Papers for Tuesday, Apr 12 2022

The MIllimeter Sardinia radio Telescope Receiver based on Array of Lumped elements kids, MISTRAL, is a millimetric ($\simeq 90GHz$) multipixel camera being built for the Sardinia Radio Telescope. It is going to be a facility instrument and will sample the sky with 12 arcsec angular resolution, 4 arcmin field of view, through 408 Kinetic Inductance Detectors (KIDs). The construction and the beginning of commissioning is planned to be in 2022. MISTRAL will allow the scientific community to propose a wide variety of scientific cases including protoplanetary discs study, star forming regions, galaxies radial profiles, and high angular resolution measurements of the Sunyaev Zel'dovich (SZ) effect with the investigation of the morphology of galaxy cluster and the search for the Cosmic Web.

We study dark matter-helium scattering in the early Universe and its impact on constraints from cosmic microwave background (CMB) anisotropy measurements. We describe possible theoretical frameworks for dark matter-nucleon interactions via a scalar, pseudoscalar, or vector mediator; such interactions give rise to hydrogen and helium scattering, with cross sections that have a power-law dependence on relative velocity. Within these frameworks, we consider three scenarios: dark matter coupling to only neutrons, to only protons, and to neutrons and protons with equal strength. For these various cases, we use \textit{Planck} 2018 temperature, polarization, and lensing anisotropy data to place constraints on dark matter scattering with hydrogen and/or helium for dark matter masses between 10 keV and 1 TeV. For any model that permits both helium and hydrogen scattering with a non-negative power-law velocity dependence, we find that helium scattering dominates the constraint for dark matter masses well above the proton mass. Furthermore, we place the first CMB constraints on dark matter that scatters dominantly/exclusively with helium in the early Universe.

We discuss a new set of $\sim$ 500 numerical n-body calculations designed to constrain the masses and bulk densities of Styx, Nix, Kerberos, and Hydra. Comparisons of different techniques for deriving the semimajor axis and eccentricity of the four satellites favor methods relying on the theory of Lee & Peale (2006), where satellite orbits are derived in the context of the restricted three body problem (Pluto, Charon, and one massless satellite). In each simulation, we adopt the nominal satellite masses derived in Kenyon & Bromley (2019a), multiply the mass of at least one satellite by a numerical factor $f \ge 1$, and establish whether the system ejects at least one satellite on a time scale $\le$ 4.5 Gyr. When the total system mass is large ($f \gg 1$), ejections of Kerberos are more common. Systems with lower satellite masses ($ f \approx$ 1) usually eject Styx. In these calculations, Styx often `signals' an ejection by moving to higher orbital inclination long before ejection; Kerberos rarely signals in a useful way. The n-body results suggest that Styx and Kerberos are more likely to have bulk densities comparable with water ice, $\rho_{SK} \lesssim$ 2 g cm$^{-3}$, than with rock. A strong upper limit on the total system mass, $M_{SNKH} \lesssim 9.5 \times 10^{19}$ g, also places robust constraints on the average bulk density of the four satellites, $\rho_{SNKH} \lesssim$ 1.4 g cm$^{-3}$. These limits support models where the satellites grow out of icy material ejected during a major impact on Pluto or Charon.

We present a model-independent search for the gravitational wave background from cosmic domain walls (DWs) in the NANOGrav 12.5 years dataset and International PTA Data Release 2. DWs that annihilate at temperatures $\sim 20-50~\text{MeV}$ with tensions $\sim (40-100~\text{TeV})^3$ provide as good a fit to both datasets as the astrophysical background from supermassive black hole mergers. DWs may decay into the Standard Model (SM) or a dark sector. In the latter case we predict an abundance $\Delta N_{\text{eff}}$ of dark radiation well within the reach of upcoming CMB surveys. Complementary signatures at colliders and laboratories can arise if couplings to the SM are present. As an example, we discuss heavy axion scenarios, where DW annihilation may interestingly be induced by QCD confinement.

Radio interferometers, which are designed to observe astrophysical objects in the universe, can also be used to study the Earth's ionosphere. Radio interferometers like the Giant Metrewave Radio Telescope (GMRT) detect variations in ionospheric total electron content (TEC) on a much wider spatial scale at a relatively higher sensitivity than traditional ionospheric probes like the Global Navigation Satellite System (GNSS). The hybrid configuration of the GMRT (compact core and extended arms) and its geographical location make this interferometer an excellent candidate to explore the sensitive regions between the northern crest of the Equatorial Ionization Anomaly (EIA) and the magnetic equator. For this work, a bright radio source, 3C68.2, is observed from post-midnight to post-sunrise ($\sim$\,9 hours) to study the ionospheric activities at solar-minima. This study presents data reduction and processing techniques to measure differential TEC ($\delta\rm{TEC}$) between the set of antennas with an accuracy of $1\times10^{-3}$ TECU. Furthermore, using these $\delta\rm{TEC}$ measurements, we have demonstrated techniques to compute the TEC gradient over the full array and micro-scale variation in two-dimensional TEC gradient surface. These variations are well equipped to probe ionospheric plasma, especially during the night-time. Our study, for the first time, reports the capability of the GMRT to detect ionospheric activities. Our result validates, compared to previous studies with VLA, LOFAR and MWA, the ionosphere over the GMRT is more active, which is expected due to its location near the magnetic equator.

Solar System observations that serve as analogs for exoplanet remote sensing data can provide important opportunities to validate ideas and models related to exoplanet environments. Critically, and unlike true exoplanet observations, Solar System analog data benefit from available high-quality ground- or orbiter-derived "truth" constraints that enable strong validations of exoplanet data interpretation tools. In this work, we first present a versatile atmospheric retrieval suite, capable of application to reflected light, thermal emission, and transmission observations. The tool -- dubbed rfast -- is designed, in part, to enable exoplanet mission concept feasibility studies. Following model validation, the retrieval tool is applied to a range of Solar System analog observations for exoplanet environments. Retrieval studies using Earth reflected light observations from NASA's EPOXI mission provide a key proof-of-concept for under-development exo-Earth direct imaging concept missions. Inverse modeling applied to an infrared spectrum of Earth from the Mars Global Surveyor Thermal Emission Spectrometer achieves good constraints on atmospheric gases, including many biosignature gases. Finally, retrieval analysis applied to a transit spectrum of Titan derived from the Cassini Visual and Infrared Mapping Spectrometer provides a proof-of-concept for interpreting more feature-rich transiting exoplanet observations from NASA's James Webb Space Telescope (JWST). In the future, Solar System analog observations for exoplanets could be used to verify exoplanet models and parameterizations, and future exoplanet analog observations of any Solar System worlds from planetary science missions should be encouraged.

Dark Matter (DM) annihilation in our Galaxy may produce a linearly polarized synchrotron signal. We use, for the first time, synchrotron polarization to constrain the DM annihilation cross section by comparing theoretical predictions with the latest polarization maps obtained by the Planck satellite collaboration. We find that synchrotron polarization is typically more constraining than synchrotron intensity by about one order of magnitude, independently of uncertainties in the modeling of electron and positron propagation, or of the Galactic magnetic field. Our bounds compete with Cosmic Microwave Background limits in the case of leptophilic DM.

Large and less-biased samples of star-forming galaxies are essential to investigate galaxy evolution. H{\alpha} emission line is one of the most reliable tracers of star-forming galaxies because its strength is directly related to recent star formation. However, it is observationally expensive to construct large samples of H{\alpha} emitters by spectroscopic or narrow-band imaging survey at high-redshifts. In this work, we demonstrate a method to extract H{\alpha} fluxes of galaxies at z = 2.1-2.5 from Ks broad-band photometry of ZFOURGE catalog. Combined with 25-39 other filters, we estimate the emission line fluxes by SED fitting with stellar population models that incorporate emission-line strengths. 2005 galaxies are selected as H{\alpha} emitters by our method and their fluxes show good agreement with previous measurements in the literature. On the other hand, there are more H{\alpha} luminous galaxies than previously reported. The discrepancy can be explained by extended H{\alpha} profiles of massive galaxies and a luminosity dependence of dust attenuation, which are not taken into account in the previous work. We also find that there are a large number of low-mass galaxies with much higher specific star formation rate (sSFR) than expected from the extrapolated star formation main sequence. Such low-mass galaxies exhibit larger ratios between H{\alpha} and UV fluxes compared to more massive high sSFR galaxies. This result implies that a "starburst" mode may differ among galaxies: low-mass galaxies appear to assemble their stellar mass via short-duration bursts while more massive galaxies tend to experience longer-duration (> 10 Myr) bursts.

The pathway to forming the iron-rich planet Mercury remains mysterious. Mercury's core makes up 70% of the planetary mass, which implies a significant enrichment of iron relative to silicates, while its mantle is strongly depleted in oxidized iron. The high core mass fraction is traditionally ascribed to evaporative loss of silicates, e.g. following a giant impact, but the high abundance of moderately volatile elements in the mantle of Mercury is inconsistent with reaching temperatures much above 1,000 K during its formation. Here we explore the nucleation of solid particles from a gas of solar composition that cools down in the hot inner regions of the protoplanetary disc. The high surface tension of iron causes iron particles to nucleate homogeneously (i.e., not on a more refractory substrate) under very high supersaturation. The low nucleation rates lead to depositional growth of large iron pebbles on a sparse population of nucleated iron nano-particles. Silicates in the form of iron-free MgSiO$_3$ nucleate at similar temperatures but obtain smaller sizes due to the much higher number of nucleated particles. This results in a chemical separation of large iron particles from silicate particles with ten times lower Stokes numbers. We propose that such conditions lead to the formation of iron-rich planetesimals by the streaming instability. In this view, Mercury formed by accretion of iron-rich planetesimals with a sub-solar abundance of highly reduced silicate material. Our results imply that the iron-rich planets known to orbit the Sun and other stars are not required to have experienced mantle-stripping impacts. Instead their formation could be a direct consequence of temperature fluctuations in protoplanetary discs and chemical separation of distinct crystal species through the ensuing nucleation process.

The ideal exoplanets to search for life are those within a star's habitable zone. However, even within the habitable zone planets can still develop uninhabitable climate states. Sustaining a temperate climate over geologic ($\sim$Gyr) timescales requires a planet contain sufficient internal energy to power a planetary-scale carbon cycle. A major component of a rocky planet's energy budget is the heat produced by the decay of radioactive elements, especially $^{40}$K, $^{232}$Th, $^{235}$U and $^{238}$U. As the planet ages and these elements decay, this radiogenic energy source dwindles. Here we estimate the probability distribution of the amount of these heat producing elements (HPEs) that enter into rocky exoplanets through Galactic history, by combining the system-to-system variation seen in stellar abundance data with the results from Galactic chemical evolution models. Using these distributions, we perform Monte-Carlo thermal evolution models that maximize the mantle cooling rate. This allows us to create a pessimistic estimate of lifetime a rocky, stagnant-lid exoplanet can support a global carbon cycle and temperate climate as a function of its mass and when it in Galactic history. We apply this framework to a sample of 17 likely rocky exoplanets with measured ages, 7 of which we predict are likely to be actively degassing today despite our pessimistic assumptions. For the remaining planets, including those orbiting TRAPPIST-1, we cannot confidently assume they currently contain sufficient internal heat to support mantle degassing at a rate sufficient to sustain a global carbon cycle or temperate climate without additional tidal heating or undergoing plate tectonics.

We present a new code named pyHIIextractor, which detects and extracts the main features (positions and radii) of clumpy ionized regions, i.e. candidate HII regions, using H{\alpha} emission line images. Our code is optimized to be used on the dataproducts provided by the Pipe3D pipeline (or dataproducts with such a format), applied to high spatial resolution Integral Field Spectroscopy data (like that provided by the AMUSING++ compilation, using MUSE). The code provides the properties of both the underlying stellar population and the emission lines for each detected H ii candidate. Furthermore, the code delivers a novel estimation of the diffuse ionized gas (DIG) component, independent of its physical properties, which enables a decontamination of the properties of the HII regions from the DIG. Using simulated data, mimicking the expected observations of spiral galaxies, we characterise pyHIIextractor and its ability to extract the main properties of the H ii regions (and the DIG), including the line fluxes, ratios and equivalent widths. Finally, we compare our code with other such tools adopted in the literature, which have been developed or used for similar purposes: pyhIIexplorer, SourceExtractor, HIIphot, and astrodendro. We conclude that pyHIIextractor exceeds the performance of previous tools in aspects such as the number of recovered regions and the distribution of sizes and fluxes (an improvement that is especially noticeable for the faintest and smallest regions). pyHIIextractor is therefore an optimals tool to detect candidate HII regions, offering an accurate estimation of their properties and a good decontamination of the DIG component.

Methane has been proposed as an exoplanet biosignature. Imminent observations with the James Webb Space Telescope may enable methane detections on potentially habitable exoplanets, so it is essential to assess in what planetary contexts methane is a compelling biosignature. Methane's short photochemical lifetime in terrestrial planet atmospheres implies that abundant methane requires large replenishment fluxes. While methane can be produced by a variety of abiotic mechanisms such as outgassing, serpentinizing reactions, and impacts, we argue that, in contrast to an Earth-like biosphere, known abiotic processes cannot easily generate atmospheres rich in CH$_4$ and CO$_2$ with limited CO due to the strong redox disequilibrium between CH$_4$ and CO$_2$. Methane is thus more likely to be biogenic for planets with 1) a terrestrial bulk density, high mean-molecular-weight and anoxic atmosphere, and an old host star; 2) an abundance of CH$_4$ that implies surface fluxes exceeding what could be supplied by abiotic processes; and 3) atmospheric CO$_2$ with comparatively little CO.

Are WO-type Wolf Rayet (WR) stars in the final stage of massive star evolution before core-collapse? Although WC- and WO-type WRs have very similar spectra, WOs show a much stronger O VI $\lambda \lambda$3811,34 emission-line feature. This has usually been interpreted to mean that WOs are more oxygen rich than WCs, and thus further evolved. However, previous studies have failed to model this line, leaving the relative abundances uncertain, and the relationship between the two types unresolved. To answer this fundamental question, we modeled six WCs and two WOs in the LMC using UV, optical, and NIR spectra with the radiative transfer code CMFGEN in order to determine their physical properties. We find that WOs are not richer in oxygen; rather, the O VI feature is insensitive to the abundance. However, the WOs have a significantly higher carbon and lower helium content than the WCs, and hence are further evolved. A comparison of our results with single-star Geneva and binary BPASS evolutionary models show that while many properties match, there is more carbon and less oxygen in the WOs than either set of evolutionary model predicts. This discrepancy may be due to the large uncertainty in the $^{12}$C+$^4$He$\rightarrow^{16}$O nuclear reaction rate; we show that if the Kunz et al. rate is decreased by a factor of 25-50%, then there would be a good match with the observations. It would also help explain the LIGO/VIRGO detection of black holes whose masses are in the theoretical upper mass gap.

We complete the analysis of all 2018 prime-field microlensing planets identified by the KMTNet AnomalyFinder. Among the 10 previously unpublished events with clear planetary solutions, 8 are either unambiguously planetary or are very likely to be planetary in nature: OGLE-2018-BLG-1126, KMT-2018-BLG-2004, OGLE-2018-BLG-1647, OGLE-2018-BLG-1367, OGLE-2018-BLG-1544, OGLE-2018-BLG-0932, OGLE-2018-BLG-1212, and KMT-2018-BLG-2718. Combined with the 4 previously published new AnomalyFinder events and 12 previously published (or in preparation) planets that were discovered by eye, thismakes a total of 24 2018 prime-field planets discovered or recovered by AnomalyFinder. Together with a paper in preparation on 2018 sub-prime planets, this work lays the basis for the first statistical analysis of the planet mass-ratio function based on planets identified in KMTNet data. By systematically applying the heuristic analysis of Hwang et al. (2022) to each event, we identify the small modification in their formalism that is needed to unify the so-called close/wide and inner/outer degeneracies, as conjectured by

Despite that many methods have been developed to find young massive planets in protoplanetary discs, it is still challenging to directly detect low-mass planets that are embedded in discs. On the other hand, the core-accretion theory suggests that there could be a large population of embedded low-mass young planets at the Kelvin-Helmholtz (KH) contraction phase. We adopt both 1-D models and 3-D simulations to calculate the envelope structure around low-mass cores (several to tens of $M_{\oplus}$) with different luminosities, and derive their thermal fluxes at radio wavelengths. We find that, when the disc is optically thin at these wavelengths, the observation can see through the disc and probe the denser envelope within the planet's Hill sphere. When the optically thin disc is observed with the resolution reaching one disc scale height, the radio thermal flux from the planetary envelope around a 10 M$_{\oplus}$ core is more than 10\% higher than the flux from the background disc. The emitting region can be extended and elongated. Finally, our model suggests that the au-scale excess detected by ALMA at 52 au in TW Hydrae disc is consistent with the envelope of an embedded 10-20 $M_{\oplus}$ planet, which can explain the detected flux, the spectral index dip, and the tentative spirals. The observation is also consistent with the planet undergoing pebble accretion. Future ALMA and ngVLA observations may directly reveal more such low-mass planets, enabling us to study core growth and even reconstruct the planet formation history using the "protoplanet" population.

We use full-disk, SOHO/EIT 195 $\AA$ calibrated images to measure latitudinal and day to day variations of area and average photon fluxes of the near equatorial coronal holes. In addition, energy emitted by the coronal holes with their temperature and strength of magnetic field structures are estimated. By analyzing data of 2001-2008, we find that variations of average area (A), photon flux (F), radiative energy (E) and temperature (T) of coronal holes are independent of latitude. Whereas inferred strength of magnetic field structure of the coronal holes is dependent on the latitudes and varies from low near the equator to high near both the poles. Typical average values of estimated physical parameters are: $A \sim 3.8(\pm0.5)\times10^{20}~cm^{2}, F \sim 2.3(\pm0.2)\times10^{13}~photons\;cm^{-2}\;sec^{-1}, E \sim 2.32(\pm0.5)\times 10^{3}~ergscm^{-2}sec^{-1} \ and \ T \sim 0.94(\pm0.1)\times10^{6} ~$ K. Average strength of magnetic field structure of coronal hole at the corona is estimated to be $\sim$ $0.08 \pm 0.02$ Gauss. If coronal holes are anchored in the convection zone, one would expect they should rotate differentially. Hence, thermal wind balance and isorotation of coronal holes with the solar plasma implies the temperature difference between the equator and both the poles. Contrary to this fact, variation of thermal structure of near equatorial coronal holes is independent of latitude leading to a conclusion that coronal holes must rotate rigidly that are likely to be anchored initially below the tachocline confirming our previous study (ApJ, 763, 137, 2013).

Gamma-ray bursts (GRB) are short and intense bursts of $\sim$100 keV$-$1MeV photons, usually followed by long-lasting decaying afterglow emission in a wide range of electromagnetic wavelengths from radio to X-ray and, sometimes, even to GeV gamma-rays. These emissions are believed to originate from a relativistic jet, which is driven due to the collapse of special massive stars and the mergers of compact binaries (i.e., double neutron stars or a neutron star and a black hole). This chapter first briefly introduces the basic observational facts of the GRB phenomena, including the prompt emission, afterglow emission, and host galaxies. Secondly, a general theoretical understanding of the GRB phenomena is described based on a relativistic jet's overall dynamical evolution, including the acceleration, propagation, internal dissipation, and deceleration phases. Here a long-lasting central engine of the GRBs can substantially influence the dynamical evolution of the jet. In addition, a supernova/kilonova emission can appear in the optical afterglow of some nearby GRBs, which can provide an important probe to the nature of the GRB progenitors. Finally, as luminous cosmological phenomena, it is expected to use GRBs to probe the early universe and to constrain the cosmological parameters.

We present \textsc{MG-MAMPOSSt}, a license-free code to constrain modified gravity models by reconstructing the mass profile of galaxy clusters with the kinematics of the cluster's member galaxies. We describe the main features of the code and we show the capability of the method when the kinematic information is combined with lensing data. We discuss recent results and forecasts on two classes of models currently implemented in the code, characterized by different screening mechanisms, namely, chameleon and Vainshtein screening. We further explore the impact of possible systematics in view of application to the data from upcoming surveys. This proceedings summarizes the results presented at the ALTECOSMOFUN workshop in September 2021.

We demonstrate the effectiveness of a Bayesian evidence-based analysis for diagnosing and disentangling the sky-averaged 21-cm signal from instrumental systematic effects. As a case study, we consider a simulated REACH pipeline with an injected systematic. We demonstrate that very poor performance or erroneous signal recovery is achieved if the systematic remains unmodelled. These effects include sky-averaged 21-cm posterior estimates resembling a very deep or wide signal. However, when including parameterised models of the systematic, the signal recovery is dramatically improved in performance. Most importantly, a Bayesian evidence-based model comparison is capable of determining whether or not such a systematic model is needed as the true underlying generative model of an experimental dataset is in principle unknown. We, therefore, advocate a pipeline capable of testing a variety of potential systematic errors with the Bayesian evidence acting as the mechanism for detecting their presence.

The Laser Interferometer Space Antenna (LISA), which is currently under construction, is designed to measure gravitational wave signals in the milli-Hertz frequency band. It is expected that tens of millions of Galactic binaries will be the dominant sources of observed gravitational waves. The Galactic binaries producing signals at mHz frequency range emit quasi monochromatic gravitational waves, which will be constantly measured by LISA. To resolve as many Galactic binaries as possible is a central challenge of the upcoming LISA data set analysis. Although it is estimated that tens of thousands of these overlapping gravitational wave signals are resolvable, and the rest blurs into a galactic foreground noise; extracting tens of thousands of signals using Bayesian approaches is still computationally expensive. We developed a new end-to-end pipeline using Gaussian Process Regression to model the log-likelihood function in order to rapidly compute Bayesian posterior distributions. Using the pipeline we are able to solve the Lisa Data Challange (LDC) 1-3 consisting of noisy data as well as additional challenges with overlapping signals and particularly faint signals.

The limiting mass of cold white dwarfs was first calculated by E. Stoner in an approximate model of a uniform star and was soon reduced by ~20% in papers by S. Chandrasekhar and L. D. Landau based on an exact solution of the equations for the stellar equilibrium. Here we examine uniform models of white dwarfs taking general relativity effects and the influence of finite temperature into account. Solutions are obtained in the form of finite analytic formulas and, for masses differing by no more than ~20% from the exact solutions, found by numerical integration of the differential equations for the stellar equilibrium.

We demonstrate that the Bayesian evidence can be used to find a good approximation of the true likelihood function of a dataset, a goal of the likelihood-free inference (LFI) paradigm. As a concrete example, we use forward modelled sky-averaged 21-cm signal antenna temperature datasets where we artificially inject noise structures of various physically motivated forms. We find that the Gaussian likelihood performs poorly when the noise distribution deviates from the Gaussian case e.g. heteroscedastic radiometric or heavy-tailed noise. For these non-Gaussian noise structures, we show that the generalised normal likelihood is on a similar Bayesian evidence scale with comparable sky-averaged 21-cm signal recovery as the true likelihood function of our injected noise. We therefore propose the generalised normal likelihood function as a good approximation of the true likelihood function if the noise structure is a priori unknown.

We present a directed search for continuous gravitational wave (CW) signals emitted by spinning neutron stars located in the inner parsecs of the Galactic Center (GC). Compelling evidence for the presence of a numerous population of neutron stars has been reported in the literature, turning this region into a very interesting place to look for CWs. In this search, data from the full O3 LIGO--Virgo run in the detector frequency band $[10,2000]\rm~Hz$ have been used. No significant detection was found and 95$\%$ confidence level upper limits on the signal strain amplitude were computed, over the full search band, with the deepest limit of about $7.6\times 10^{-26}$ at $\simeq 142\rm~Hz$. These results are significantly more constraining than those reported in previous searches. We use these limits to put constraints on the fiducial neutron star ellipticity and r-mode amplitude. These limits can be also translated into constraints in the black hole mass -- boson mass plane for a hypothetical population of boson clouds around spinning black holes located in the GC.

We explore the properties of central galaxies living in voids using the EAGLE cosmological hydrodynamic simulations in conjunction with the \cite{paillas2017} void catalogue. Based on the distance to the closest void, we define four galaxy samples: inner void, outer void, wall, and skeleton. We find that inner void galaxies with host halo masses $<10^{12}M_\odot$ have lower stellar mass and stellar mass fractions than those in denser environments, and the fraction of galaxies with star formation (SF) activity and atomic hydrogen (HI) gas decreases with increasing distance to their closest void, in agreement with observations. To mitigate the influence of stellar mass, we compare inner void galaxies to subsamples of fixed stellar (halo) mass. Inner void galaxies with $M_{*}= 10^{[9.0-9.5]}M_\odot$ have similar SF activity and HI gas fractions, but the lowest quenched galaxy fraction. Inner void galaxies with $M_{*}= 10^{[9.5-10.5]}M_\odot$ have the lowest HI gas fraction, the highest quenched fraction, and lowest gas metallicities compared to galaxies in denser environments. On the other hand, inner void galaxies with $M_{*}>10^{10.5}M_\odot$ have comparable SF activity and HI gas fractions to their analogues in denser environments. They retain the highest metallicity gas that might be linked to physical processes only fostered in underdense regions. Furthermore, inner void galaxies have the lowest fraction of positive gas-phase metallicity gradients, which are typically associated with external processes or feedback events, suggesting they have more quiet merger histories than galaxies in denser environments. Our findings provide a unique insight into how galaxies are affected by their large-scale environment.

Multiple observations made by several different telescopes have shown asymmetry between the number of spiral galaxies rotating in opposite directions in different parts of the sky. One of the immediate questions regarding the possible asymmetry of the spin directions is whether the distribution forms a cosmological-scale axis. This paper analyzes and compares 10 different datasets published in the past decade, collected by SDSS, Pan-STARRS, and Hubble Space Telescope. The datasets contain spiral galaxies separated by their spin direction, and the distribution can show dipole axes. The analysis shows that the directions of the most probable dipole axes are consistent in datasets that have similar average redshift, but different between datasets that have different average redshift. The analysis also shows that the location of the most probable axis correlates with the average redshift of the galaxies in the datasets. That is, the location of the most probable axis shifts when the redshift gets higher, and the correlation is statistically significant. This provides a certain indication of a drift in a possible axis formed by the distribution of galaxy spin directions, or a cosmological scale structure that peaks at a certain distance from Earth.

Coronal magnetic fields are well known to be one of the crucial parameters defining coronal physics and space weather. However, measuring the global coronal magnetic fields remains challenging. The polarization properties of coronal radio emissions are sensitive to coronal magnetic fields. While they can prove to be useful probes of coronal and heliospheric magnetic fields, their usage has been limited by technical and algorithmic challenges. We present a robust algorithm for precise polarization calibration and imaging of low-radio frequency solar observations and demonstrate it on data from the Murchison Widefield Array, a Square Kilometer Array (SKA) precursor. This algorithm is based on the {\it Measurement Equation} framework, which forms the basis of all modern radio interferometric calibration and imaging. It delivers high dynamic range and fidelity full Stokes solar radio images with instrumental polarization leakages $<1\%$, on par with general astronomical radio imaging, and represents the state-of-the-art. Opening up this rewarding, yet unexplored, phase space will enable multiple novel science investigations and offer considerable discovery potential. Examples include detection of low-level of circularly polarization from thermal coronal emission to estimate large-scale quiescent coronal fields; polarization of faint gyrosynchrotron emissions from coronal mass ejections for robust estimation of plasma parameters; and detection of the first-ever linear polarization at these frequencies. This method has been developed with the SKA in mind and will enable a new era of high fidelity spectro-polarimetric snapshot solar imaging at low-radio frequencies.

Gaia19bxc is a transient detected on 2019 May 9 by the Gaia Photometric Science Alerts Team. I analyzed the past public Zwicky Transient Facility (ZTF) data and found that Gaia19bxc has a period of 0.04473647(3) d and two different maxima in one cycle. This object also showed high and low states in the ZTF data. Based on the high amplitude (2.0 mag) of the short-term variations, short period, almost zero color indices between the different ZTF bands and the absence of a longer period, I classified it to be a likely polar. There has been no established polar below the period minimum of cataclysmic variables (CVs) and Gaia19bxc could be the first such an object. CVs below the period minimum usually have a secondary star with a stripped evolved core and Gaia19bxc is expected to have a similar secondary. If this is indeed the case, Gaia19bxc could become a highly magnetized exotic ultracompact binary during its secular evolution.

We present ALMA archival data for 219-235 GHz continuum and line observations toward the hot molecular core (HMC) W49N MCN-a (UCHII region J1) at a resolution of ~0."3. The dust continuum emission, showing an elongated structure of 1."40x0."95 (PA=43.5deg) perpendicular to the outflow seen in SiO and SO, represents a rotating flattened envelope, or torus, with a radius of 7,800au inclined at 47.5deg or larger. The emissions from CH3CN and 11 molecular lines exhibit a consistent velocity gradient as a result of rotation. The magnitude of each velocity gradient is different, reflecting that each line samples a specific radial region. This allows us to derive a rotation curve as Vrot prop R^0.44+-0.11 for 2,400au < R < 14,000au, giving the dynamical mass as Mdyn = 57.0+24.5-17.1 (R [au]/3, 000)^1.88 Msun. The envelope mass independently estimated from the dustemission is 910Msun (for Tdust =180K) for R<7,800au and 32Msun (for Tdust=300K) for R<1,700 au. The dynamical mass formula agrees well with these mass estimates within an uncertainty of a factor of three in the latter. The envelope is self-gravitating and is unstable to form spiral arms and fragments, allowing rapid accretion to the inner radii with a rate of order 10^-2 Msun yr^-1, although inward motion was not detected. The envelope may become a non self-gravitating Keplerian disk at R<(300-1,000) au. The formula is also consistent with the total mass ~10^4 Msun of the entire HMC 0.15 pc (31,000 au) in radius. Multiple transitions of CH3CN, HNCO and CH3OH provide the rotation temperatures, suggesting that the central source of MCN-a has an intrinsic bolometric luminosity of ~10^6 Lsun. These results have revealed the structure and kinematics of MCN-a at its intermediate radii. With no broad-line H30alpha emission detected, MCN-a may be in the earliest phase of massive star formation.

A recent result from Fermilab suggests that the measured W-boson mass deviates from the prediction of the Standard Model with a significance of $>7\sigma$, and there may exist new physics beyond the SM. It is proposed that the inert two Higgs doublet model (i2HDM) can well explain the new W-boson mass. The preferred dark matter mass is between 54 and 74 GeV. More interestingly, it is found that part of the parameter space of this model can explain both the Galactic center GeV gamma-ray excess detected by Fermi-LAT and the GeV antiproton excess detected by AMS-02 through a $SS\rightarrow WW^*$ annihilation. In this paper, we aim to test this model using the Fermi-LAT observation of Milky Way dwarf spheroidal (dSph) galaxy. We mainly focus on the four nearest confirmed dSphs, which are among the dSphs with the largest J-factors. We find that our constraints are above the favored parameters and can not exclude such a model, suggesting i2HDM is a promising model that can interpret the W-boson mass anomaly, GeV excess, and antiproton excess.

The multi-messenger gravitational-wave (GW) observation for binary neutron star merger events could provide a rather useful tool to explore the evolution of the universe. In particular, for the third-generation GW detectors, i.e., the Einstein Telescope (ET) and the Cosmic Explorer (CE), proposed to be built in Europe and the U.S., respectively, lots of GW standard sirens with known redshifts could be obtained, which would exert great impacts on the cosmological parameter estimation. The total neutrino mass could be measured by cosmological observations, but such a measurement is model-dependent and currently only gives an upper limit. In this work, we wish to investigate whether the GW standard sirens observed by ET and CE could help improve the constraint on the neutrino mass, in particular in the interacting dark energy (IDE) models. We find that the GW standard siren observations from ET and CE can only slightly improve the constraint on the neutrino mass in the IDE models, compared to the current limit. The improvements in the IDE models are weaker than in the standard cosmological model. Although the limit on neutrino mass can only be slightly updated, the constraints on other cosmological parameters can be significantly improved by using the GW observations.

The evolution of magnetism in late-type dwarfs remains murky, as we can only weakly predict levels of activity for M dwarfs of a given mass and age. We report results from our spectroscopic survey of M dwarfs in the Southern Continuous Viewing Zone (CVZ) of the Transiting Exoplanet Survey Satellite (TESS). As the TESS CVZs overlap with those of the James Webb Space Telescope, our targets constitute a legacy sample for studies of nearby M dwarfs. For 122 stars, we obtained at least one $R\approx 2000$ optical spectrum with which we measure chromospheric $\mathrm{H}\alpha$ emission, a proxy for magnetic field strength. The fraction of active stars is consistent with what is expected for field M dwarfs; as in previous studies, we find that late-type M dwarfs remain active for longer than their early type counterparts. While the TESS light curves for $\approx$20% of our targets show modulations consistent with rotation, TESS systematics are not well enough understood for confident measurements of rotation periods ($P_{\mathrm{rot}}$) longer than half the length of an observing sector. We report periods for 12 stars for which we measure $P_{\mathrm{rot}} {\lower0.8ex\hbox{$\buildrel <\over\sim$}}$ 15 d or find confirmation for the TESS-derived $P_{\mathrm{rot}}$ in the literature. Our sample of 21 $P_{\mathrm{rot}}$, which includes periods from the literature, is consistent with our targets being spun-down field stars. Finally, we examine the $\mathrm{H}\alpha$-to-bolometric luminosity distribution for our sample. Two stars are rotating fast enough to be magnetically saturated, but are not, hinting at the possibility that fast rotators may appear inactive in $\mathrm{H}\alpha$.

Blazars are a class of AGN, one of their jets is pointed towards the earth. Here, we report about the multi-wavelength study for blazar S5 1803+78 between MJD 58727 to MJD 59419. We analysed $\gamma$-ray data collected by Fermi-LAT, X-ray data collected by Swift-XRT \& NuSTAR, optical photons detected by Swift-UVOT \& TUBITAK observatory in Turkey. Three flaring states are identified by analysing the $\gamma$-ray light curve. A day scale variability is observed throughout the flares with the similar rise and decay times suggesting a compact emission region located close to the central engine. Cross-correlation studies are carried out between $\gamma$-ray, radio, and X-ray bands, and no significant correlation is detected. The $\gamma$-ray and optical emission are significantly correlated with zero time lag suggesting a co-spatial origin of them. A significant positive correlation between the R-I index and the V magnitude is observed. The broadband spectral energy distributions (SEDs) modeling was performed for all the flaring episodes as well as for one quiescent state for comparison. SEDs are best fitted with the synchrotron-self Compton (SSC) model under a one-zone leptonic scenario. The SED modeling shows that to explain the high flaring state strong Doppler boosting is required.

We have conducted an extensive search for nitrogen-, oxygen- and sulfur-bearing heterocycles toward Taurus Molecular Cloud 1 (TMC-1) using the deep, broadband centimeter-wavelength spectral line survey of the region from the GOTHAM large project on the Green Bank Telescope. Despite their ubiquity in terrestrial chemistry, and the confirmed presence of a number of cyclic and polycyclic hydrocarbon species in the source, we find no evidence for the presence of any heterocyclic species. Here, we report the derived upper limits on the column densities of these molecules obtained by Markov Chain Monte Carlo (MCMC) analysis and compare this approach to traditional single-line upper limit measurements. We further hypothesize why these molecules are absent in our data, how they might form in interstellar space, and the nature of observations that would be needed to secure their detection.

UNIONS is an ongoing collaboration that will provide the largest deep photometric survey of the Northern sky in four optical bands to date. As part of this collaboration, CFIS is taking $r$-band data with an average seeing of 0.65 arcsec, which is complete to magnitude 24.5 and thus ideal for weak-lensing studies. We perform the first weak-lensing analysis of CFIS $r$-band data over an area spanning 1700 deg$^2$ of the sky. We create a catalogue with measured shapes for 40 million galaxies, corresponding to an effective density of 6.8 galaxies per square arcminute, and demonstrate a low level of systematic biases. This work serves as the basis for further cosmological studies using the full UNIONS survey of 4800 deg$^2$ when completed. Here we present ShapePipe, a newly developed weak-lensing pipeline. This pipeline makes use of state-of-the-art methods such as Ngmix for accurate galaxy shape measurement. Shear calibration is performed with metacalibration. We carry out extensive validation tests on the Point Spread Function (PSF), and on the galaxy shapes. In addition, we create realistic image simulations to validate the estimated shear. We quantify the PSF model accuracy and show that the level of systematics is low as measured by the PSF residuals. Their effect on the shear two-point correlation function is sub-dominant compared to the cosmological contribution on angular scales <100 arcmin. The additive shear bias is below 5x$10^{-4}$, and the residual multiplicative shear bias is at most $10^{-3}$ as measured on image simulations. Using COSEBIs we show that there are no significant B-modes present in second-order shear statistics. We present convergence maps and see clear correlations of the E-mode with known cluster positions. We measure the stacked tangential shear profile around Planck clusters at a significance higher than $4\sigma$.

In this work, we point out that the Q/U Stokes parameters and E/B mode polarizations are the four components of a unique quaternion, which describes at the same time the directions and the parity states of spherical linear polarizations. We then point out that, with this polarization quaternion, the mathematical form of all Q/U and E/B transforms are greatly simplified, to an extent that requires only one quaternion multiplication for each transform. A preliminary application of the polarization quaternion is shown as an example to detect peculiar pixel domain patterns within the E- and B-families, which are the former and latter halves of the polarization quaternion.

The stellar surface density within the inner 1 kpc ($\Sigma_{1}$) has become a popular tool for understanding the growth of galaxies and its connection with the quenching of star formation. The emerging picture suggests that building a central dense core is a necessary condition for quenching. However, it is not clear whether changes in $\Sigma_{1}$ trace changes in stellar kinematics and the growth of dispersion-dominated bulges. In this paper, we combine imaging from the Sloan Digital Sky Survey with stellar kinematics from the Sydney-AAO Multi-object Integral-field unit (SAMI) and Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) surveys to quantify the correlation between $\Sigma_{1}$ and the proxy for stellar spin parameter within one effective radius ($\lambda_{re}$) for 1599 nearby galaxies. We show that, on the star-forming main sequence and at fixed stellar mass, changes in $\Sigma_{1}$ are mirrored by changes in $\lambda_{re}$. While forming stars, main sequence galaxies remain rotationally-dominated systems, with their $\Sigma_{1}$ increasing but their stellar spin staying either constant or slightly increasing. The picture changes below the main sequence, where $\Sigma_{1}$ and $\lambda_{re}$ are no longer correlated. Passive systems show a narrower range of $\Sigma_{1}$, but a wider range of $\lambda_{re}$ compared to star-forming galaxies. Our results indicate that, from a structural point of view, passive galaxies are a more heterogeneous population than star-forming systems, and may have followed a variety of evolutionary paths. This also suggests that, if dispersion-dominated bulges still grow significantly at $z\sim$0, this generally takes place during, or after, the quenching phase.

Solar cycle is modeled as a forced and damped harmonic oscillator and the amplitudes, frequencies, phases and decay factors of such a harmonic oscillator are estimated by non-linear fitting the equation of sinusoidal and transient parts to the sunspot and irradiance (proxy for the sunspot) data for the years 1700-2008. We find that:(i) amplitude and frequency (or period of $\sim$11 yr) of the sinusoidal part remain constant for all the solar cycles; (ii) the amplitude of the transient part is phase locked with the phase of the sinusoidal part; (iii) for all the cycles, the period and decay factor (that is much less than 1) of the transient part remain approximately constant. The constancy of the amplitudes and the frequencies of the sinusoidal part and a very small decay factor from the transient part suggests that the solar activity cycle mainly consists of a persistent oscillatory part that might be compatible with long-period ($\sim$22 yr) Alfven oscillations. For all the cycles, with the estimated physical parameters (amplitudes, phases and periods) and, by an autoregressive model, we forecast (especially for coming solar cycle 25) and backcast (to check whether Maunder minimum type solar activity exists or not) the solar cycles. We find that amplitude of coming solar cycle 25 is almost same as the amplitude of the previous solar cycle 24. We also find that sun might not have experienced a deep Maunder minimum (MM) type of activity during 1645-1700 AD corroborating some of the paleoclimatic inferences and, MM type of activity will not be imminent in near future, until at least 200 years.

The detection of multi-TeV gamma-rays from the afterglow phase of GRB 190829A by High Energy Stereoscopic System (H.E.S.S.) telescope is an addition to the already existing list of two more GRBs observed in the very high energy (VHE) gamma-rays in recent years. Jets of blazars and GRBs have many similarities and the photohadronic model is very successful in explaining the VHE gamma-ray spectra from the high energy blazars. Recently, the photohadronic model has been successfully applied to study the sub-TeV gamma-rays from the afterglow phases of GRB 180720B and GRB 190114C. We employed this model again to explain the VHE spectra observed for the two consecutive nights from GRB 190829A. We show that the spectra of GRB 190829A can be due to the interactions of high energy protons with the synchrotron self-Compton photons in the forward shock region of the GRB jet, similar to the low emission state of the VHE flaring events of high energy blazars. We speculate that, if in future, it is possible to observe the VHE gamma-ray spectra from nearby GRBs in their afterglow phases, then some of them could only be explained by employing two different spectral indices. If confirmed, such VHE spectra could be interpreted as a result of the interactions of the high energy protons with the photons, both from the synchrotron background and the synchrotron self-Compton background in the forward shock region.

We discuss the feasibility of direct multipixel imaging of exoplanets with the solar gravitational lens (SGL) in the context of a realistic deep space mission. For this, we consider an optical telescope, placed in the image plane that forms in the strong interference region of the SGL. We consider an Earth-like exoplanet located in our immediate stellar neighborhood and model its characteristics using our own Earth. We estimate photon fluxes from such a compact, extended, resolved exoplanet. This light appears in the form of an Einstein ring around the Sun, seen through the solar corona. The solar corona background contributes a significant amount of stochastic noise and represents the main noise source for observations utilizing the SGL. We estimate the magnitude of this noise. We compute the resulting signal-no-noise ratios and related integration times that are needed to perform imaging measurements under realistic conditions. We conclude that an imaging mission is challenging but feasible, using technologies that are either already available or in active development. Under realistic conditions, megapixel imaging of Earth-like exoplanets in our galactic neighborhood requires only weeks or months of integration time, not years as previously thought.

We study the Jeans-type gravitational instability for a self-gravitating medium composed of two species, baryonic (bright) and dark matter particles, using a hybrid quantum-classical fluid approach. Baryonic matter is treated classically, which is appropriate for most astrophysical environments, e.g., Bok globules, while dark matter is treated through a quantum hydrodynamic approach allowing for possible nonlinearities. These nonlinearities may arise in bosonic dark matter due to attractive or repulsive short-range self-interaction (attractive interaction being more relevant for axions) or from the Pauli exclusion principle for fermionic dark matter, e.g., massive neutrinos. This allows us to explore, in a very broad context, the impact of a dark matter background on the Jeans process for different scenarios discussed in the literature. In the simplest case, it is shown that the effect of a dark matter background on the Jeans mass depends only on the ratio of densities and velocity dispersions of baryonic and dark matter particles. Taking advantage of that, we confront the established stability criterion with Bok globule stability observations and show that the model adequately accounts for the data with dark matter parameters close to those predicted independently from numerical simulations.

Core collapse of massive stars leads to the different fate for various physical factors, which gives different spectra of the emitted neutrinos. We focus on the supernova relic neutrinos (SRNs) as a probe to investigate the stellar collapse fate. We present the SRN fluxes and event rate spectra at a detector for three resultant states after stellar core collapse, the typical mass neutron star, the higher mass neutron star, or the failed supernova forming a black hole, based on different nuclear equations of state. Then possible SRN fluxes are formed as mixtures of the three components. We also show the expected sensitivities at the next-generation water-based Cherenkov detectors, SK-Gd and Hyper-Kamiokande, as constraining the mixture fractions. This study provides a practical example of extracting astrophysical constraints through the SRN measurement.

Molecular hydrogen normally has only weak, quadrupole transitions between its rovibrational states, but in a static electric field it acquires a dipole moment and a set of allowed transitions. Here we use published ab initio calculations of the static electrical response tensors of the H2 molecule to construct the perturbed rovibrational eigensystem and its ground state absorptions. We restrict attention to two simple field configurations that are relevant to condensed hydrogen molecules in the interstellar medium: a uniform electric field, and the field of a point-like charge. The energy eigenstates are mixtures of vibrational and angular momentum eigenstates so there are many transitions that satisfy the dipole selection rules. We find that mixing is strongest amongst the states with high vibrational excitation, leading to hundreds of absorption lines across the optical and near infrared. These spectra are very different to that of the field-free molecule, so if they appeared in astronomical data they would be difficult to assign. Furthermore in a condensed environment the excited states likely have short lifetimes to internal conversion, giving the absorption lines a diffuse appearance. We therefore suggest electrified H2 as a possible carrier of the Diffuse Interstellar Bands (DIBs). We further argue that in principle it may be possible to account for all of the DIBs with this one carrier. However, despite electrification the transitions are not very strong and a large column of condensed H2 would be required, making it difficult to reconcile this possibility with our current understanding of the ISM.

We present the new ProFuse R package, a simultaneous spectral (ultraviolet to far infrared) and spatial structural decomposition tool that produces physical models of galaxies and their components. This combines the functionality of the recently released ProFound (for automatic source extraction), ProFit (for extended source profiling) and ProSpect (for stellar population modelling) software packages. The key novelty of ProFuse is that it generates images using a self-consistent model for the star formation and metallicity history of the bulge and disk separately, and uses target images across a range of wavelengths to define the model likelihood and optimise our physical galaxy reconstruction. The first part of the paper explores the ProFuse approach in detail, and compares results to published structural and stellar population properties. The latter part of the paper applies ProFuse to 6,664 z < 0.06 GAMA galaxies. Using re-processed ugriZYJHKs imaging we extract structural and stellar population properties for bulges and disks in parallel. As well as producing true stellar mass based mass-size relationships, we further extend this correlation to explore the third dimensions of age and gas phase metallicity. The disks in particular demonstrate strong co-dependency between mass-size-age in a well defined plane, where at a given disk stellar mass younger disks tend to be larger. These findings are in broad agreement with work at higher redshift suggesting disks that formed earlier are physically smaller.

We investigate mantle stripping giant impacts (GI) between super-Earths with masses between 1 M$_{\oplus}$ and 20 M$_{\oplus}$. We infer new scaling laws for the mass of the largest fragment and its iron mass fraction, as well as updated fitting coefficients for the critical specific impact energy for catastrophic disruption, $Q_{RD}^{*}$. With these scaling laws, we derive equations that relate the impact conditions, i.e., target mass, impact velocity and impactor-to-target mass ratio, to the mass and iron mass fraction of the largest fragment. This allows one to predict collision outcomes without performing a large suite of simulations. Using these equations we present the maximum and minimum planetary iron mass fraction as a result of collisional stripping of its mantle for a given range of impact conditions. We also infer the radius for a given mass and composition using interior structure models and compare our results to observations of metal-rich exoplanets. We find good agreement between the data and the simulated planets suggesting that GI could have played a key role in their formation. Furthermore, using our scaling laws we can further constrain the impact conditions that favour their masses and compositions. Finally, we present a flexible and easy-to-use tool that allows one to predict mass and composition of a planet after a GI for an arbitrary range of impact conditions which in turn allows to assess the role of GI in observed planetary systems.

Aims. We aim to investigate how polarimetric observations can improve our understanding of the nature and diversity of M/X-type asteroids. Methods. Polarimetric observations of the selected M/X-type asteroids were carried out at the Tohoku 0.6-m telescope at Haleakala Observatory, Hawaii (simultaneously in BVR filters), the 2-m telescope of the Bulgarian National Astronomical Observatory in Rozhen (in R filter), and the 2.15-m telescope of the Complejo Astron\'omico El Leoncito (CASLEO), Argentina (in V filter). We analysed the polarimetric characteristics of M/X-type asteroids along with the available data obtained by other techniques. Results. New polarimetric observations of 22 M/X-type asteroids combined with published observations provide a data set of 41 asteroids for which the depth of a negative polarisation branch and/or inversion angle were determined. We found that the depth of the negative polarisation branch tends to increase with decreasing steepness of the near-infrared spectra. Asteroids with a deeper negative polarisation branch tend to have a higher radar circular polarisation ratio. We show that, based on the relationship of the depth of the negative polarisation branch and inversion angle, two main sub-types can be distinguished among M-type asteroids. We suggest that these groups may be related to different surface compositions similar to (1) irons and stony-irons and (2) enstatite and iron-rich carbonaceous chondrites.

The survey of northern hemisphere were made at a frequency of 111 MHz. The total accumulation time at each point of the sky was at least one hour. 75 sources of pulse emission were detected. More then 80% of these sources are known pulsars observed in the side lobes of the antenna. From one to several hundreds pulses were detected in twelve known pulsars. For four pulsars (J0157+6212, J1910+5655, J2337+6151, J2354+6155) the narrowness of the strongest pulses and the ratio of peak flux densities in the strongest pulse and in the average profile indicate that they can be pulsars with giant pulses. One new rotating radio transient (RRAT) J0812+8626 with a dispersion measure DM=40.25 pc/cm3 was detected.

We present here the kinematics of the EUV wave associated with a GOES M1.0-class solar flare, which originates in NOAA AR 12740. The event is thoroughly observed with Atmospheric Imaging Assembly (AIA) onboard Solar Dynamics Observatory (SDO) with high spatio-temporal resolutions. This event displays many features of EUV waves, which are very decisive for the understanding of the nature of EUV waves. These features include: a fast-mode wave, a pseudo wave, a slow-mode wave and stationary fronts, probably due to mode conversion. One fast-mode wave also propagates towards the coronal hole situated close to the north pole and the wave speed does not change when it encounters the coronal hole. We intend to provide self-consistent interpretations for all these different features.

The direct detection and imaging of exoplanets requires the use of high-contrast adaptive optics(AO). In these systems quasi-static aberrations need to be highly corrected and calibrated. In order to achieve this, the pupil-modulated point-diffraction interferometer (m-PDI), was presented in an earlier paper. This present paper focuses on m-PDI concept validation through three experiments. First, the instrument's accuracy and dynamic range are characterised by measuring the spatial transfer function at all spatial frequencies and at different amplitudes. Then, using visible monochromatic light, an adaptive optics control loop is closed on the system's systematic bias to test for precision and completeness. In a central section of the pupil with 72% of the total radius the residual error is 7.7nm-rms. Finally, the control loop is run using polychromatic light with a spectral FWHM of 77nm around the R-band. The control loop shows no drop in performance with respect to the monochromatic case, reaching a final Strehl ratio larger than 0.7.

The detection of lines in emission in planetary atmospheres provides direct evidence of temperature inversion. We confirm the trend of ultra-hot Jupiters orbiting A-type stars showing temperature inversions on their daysides, by detecting metals emission lines in the dayside of KELT-20b. We first detect the planetary emission by using the G2 stellar mask of the HARPS-N pipeline, which is mainly composed of neutral iron lines, as a template. Using neutral iron templates, we perform a retrieval of the atmospheric temperature-pressure profile of the planet, confirming a thermal inversion. Then we create models of planetary emission of different species using the retrieved inverted temperature-pressure profile. By using the cross-correlation technique, we detect FeI, FeII and CrI at signal-to-noise ratio levels of 7.1, 3.9 and 3.6, respectively. The latter is detected for the first time in emission in the atmosphere of an exoplanet. Contrary to FeI, FeII and CrI are detected only after the occultation and not before, hinting for different atmospheric properties in view on the pre- and post- occultation orbital phases. A further retrieval of the temperature-pressure profile performed independently on the pre- and post- occultation phases, while not highly significant, points to a steeper thermal inversion in the post-occultation.

We present the identifications of a flux-limited sample of highly variable X-ray sources on long time-scales from the second catalogue of the XMM$-$Newton SLew survey (XMMSL2). The carefully constructed sample, comprising 265 sources (2.5 per cent) selected from the XMMSL2 clean catalogue, displayed X-ray variability of a factor of more than 10 in 0.2$-$2 keV compared to the ROSAT All Sky Survey. Of the sample sources, 94.3 per cent are identified. The identification procedure follows a series of cross-matches with astronomical data bases and multiwavelength catalogues to refine the source position and identify counterparts to the X-ray sources. Assignment of source type utilizes a combination of indicators including counterparts offset, parallax measurement, spectral colours, X-ray luminosity, and light-curve behaviour. We identified 40 per cent of the variables with stars, 10 per cent with accreting binaries, and at least 30.4 per cent with active galactic nuclei. The rest of the variables are identified as galaxies. It is found that the mean effective temperatures of the highly variable stars are lower than those of less variable stars. Our sample of highly variable AGN tend to have lower black hole masses, redshifts, and marginally lower soft X-ray luminosities compared to the less variable ones, while no difference was found in the Eddington ratio distributions. Five flaring events are tidal disruption events published previously. This study has significantly increased the number of variable sources in XMMSL2 with identifications and provides greater insight on the nature of many of the sources, enabling further studies of highly variable X-ray sources.

Duo to the large magnetic Reynolds number, the magnetic helicity originating from the solar interior can be carried away through the photosphere into the corona. However, the relationship between the accumulated magnetic helicity flux through the photosphere and the magnetic helicity in the corona is still unclear. By selecting 36 newly emerging active regions in the 23rd solar cycle, we apply optical flow methods to derive the accumulated magnetic helicity through the photosphere ($H_m^p$) by using the sequential longitudinal magnetograms, use nonlinear force-free field extrapolation to obtain the 3D coronal magnetic field, and adopt finite volume methods to calculate the instantaneous relative magnetic helicity in the corona ($H_m^c$) by using vector magnetograms. It is found that the local correlation tracking (LCT)-based $H_m^p$ is larger than $H_m^c$ in $1"$, and that the Differential Affine Velocity Estimator-based $H_m^p$ is more consistent with $H_m^c$ than the LCT-based $H_m^p$. $H_m^p$ is more consistent with $H_m^c$ in evaluation from $2"$ than from $1"$. Moreover, $H_m^c - H_m^p$ systematically shows consistency with the Hemispheric Helicity Rule (over 55\%), no matter which resolution and method are used. These estimations suggest that the consistency of $H_m^c$ and $H_m^p$ is partly dependent on the resolution of the magnetograms and the calculation methods.

We develop and apply a bespoke fitting routine to a large volume of solar wind electron distribution data measured by Parker Solar Probe (PSP) over its first five orbits, covering radial distances from 0.13 to 0.5 au. We characterise the radial evolution of the electron core, halo and strahl populations in the slow solar wind during these orbits. The fractional densities of these three electron populations provide evidence for the growth of the combined suprathermal halo and strahl populations from 0.13 to 0.17 au. Moreover, the growth in the halo population is not matched by a decrease of the strahl population at these distances, as has been reported for previous observations at distances greater than 0.3 au. We also find that the halo is negligible at small heliocentric distances. The fractional strahl density remains relatively constant ~1 % below 0.2 au, suggesting that the rise in the relative halo density is not solely due to the transfer of strahl electrons into the halo.

Galactic winds are a crucial player in galaxy formation and evolution, but observations of them have proven extraordinarily difficult to interpret, leaving large uncertainties even in basic quantities such as mass outflow rates. Part of this uncertainty arises from the relatively simplistic models to which complex wind observations are often fit, which inevitably discard much of the available information. Here we present an analysis of the wind of the nearby dwarf starburst galaxy M82 using a semi-analytic model that is able to take advantage of the full three-dimensional information present in position-position-velocity data cubes measured in the Hi 21 cm line, the CO 2-1 line, and the Ha line. Our best-fitting model produces position-dependent spectra in good agreement with the observations, and shows that the total wind mass flux in the atomic and molecular phases is approx 10 M_sun yr-1 (corresponding to a mass loading factor of about 2 - 3), with less than a factor of two uncertainty; the mass flux in the warm ionised phase is more poorly constrained, and may be comparable to or smaller than this. At least over the few kpc off the plane for which we trace the outflow, it appears to be a wind escaping the galaxy, rather than a fountain that falls back. Our fits require that clouds of cool gas entrained into the wind expand only modestly, suggesting they are magnetically confined. Finally, we demonstrate that attempts to model the wind using simplifying assumptions such as instantaneous acceleration and a constant terminal wind speed can yield significantly erroneous results.

This paper compares the population of BBH mergers detected by LIGO/Virgo with selected long GRB world models convolved with a delay function (LGRBs are used as a tracer of stellar mass BH formation). The comparison involves the redshift distribution and the fraction of LGRBs required to produce the local rate of BBH mergers. We find that BBH mergers and LGRBs cannot have the same formation history, unless BBHs mergers have a long coalescence time of several Gyr. This would imply that BHs born during the peak of long GRB formation at redshift z = 2 - 3 merge within the horizon of current GW interferometers. We also show that LGRBs are more numerous than BBH mergers, so that most of them do not end their lives in BBH mergers. We interpret these results as an indication that BBH mergers and LGRBs constitute two distinct populations of stellar mass BHs, with LGRBs being more frequent than BBH mergers. We speculate that the descendants of LGRBs may resemble galactic High-Mass X-Ray Binaries more than BBH mergers. We finally discuss the possible existence of a sub-population of fast-spinning LGRBs descendants among BBH mergers, showing that this population, if it exists, is expected to become dominant beyond redshift z = 1, leading to a change in the observed properties of BBH mergers.

In this work we report the results of a nineteen-month Fast Radio Burst observational campaign carried out with the North-South arm of the Medicina Northern Cross radio telescope at 408~MHz in which we monitored four repeating sources: FRB20180916B, FRB20181030A, FRB20200120E and FRB20201124A. We present the current state of the instrument and the detection and characterisation of three bursts from FRB20180916B. Given our observing time, our detections are consistent with the event number we expect from the known burst rate ($2.7 \pm 1.9$ above our 10$\sigma$, 38~Jy~ms detection threshold) in the 5.2 day active window of the source, further confirming the source periodicity. We detect no bursts from the other sources. We turn this result into a 95\% confidence level lower limit on the slope of the differential fluence distribution $\alpha$ to be $\alpha > 2.1$ and $\alpha > 2.2$ for FRB20181030A and FRB20200120E respectively. Given the known rate for FRB20201124A, we expect $1.0 \pm 1.1$ bursts from our campaign, consistent with our non-detection.

Redshift evolution and peculiar velocities break the isotropy of cosmological surveys with respect to the directions parallel and transverse to the line of sight, limiting the accuracy of the Fourier representation to small areas and redshift ranges. In contrast to the Fourier space power spectrum, the full information about the two-point function of tracers of large-scale structure is encapsulated in the redshift-dependent angular power spectrum $C_\ell^{ij} (z_i,z_j)$ for the tracer species $i$ and $j$ at the redshift slices $z_i$ and $z_j$, expressed in harmonic space. In this paper we derive semi-analytical expressions for the multi-tracer Fisher matrix of angular power spectra, in real and in redshift space, which are exact in the linear regime of structure formation. Our expressions can be used to forecast the constraining power of galaxy surveys with many tracers and a large number of redshift slices, for which the derivation of the Fisher matrix from numerically evaluated covariance matrices may not be feasible or practical.

We study elastic properties of solid Yukawa systems. Elastic moduli and effective shear modulus of body-centered cubic (bcc) and face-centered cubic (fcc) lattices are obtained from electrostatic energies of deformed crystals. For the bcc lattice our results are well consistent with previous calculations and improve them, while results for the fcc lattice are new. We have also obtained an analytical expression of the elastic moduli for the weak polarization and constructed a convenient approximation for the higher polarization.

PKS 0735+178 is a bright radio and $\gamma$-ray blazar that is possibly associated with multiple neutrino events observed by the IceCube, Baikal, Baksan, and KM3NeT neutrino telescopes. The source was found to undergo a major flaring activity in $\gamma$-ray, X-ray, ultraviolet (UV) and optical bands. We present a long-term detailed study of this peculiar blazar to investigate the temporal and spectral changes in the multi-wavelength emission when the neutrino events were observed. The analysis of Swift-XRT snapshots reveal a flux variability of more than a factor 2 in about $5\times10^3$ seconds during the observation on December 17, 2021. In the $\gamma$-ray band, the source was in its historical highest flux level at the time of the arrival of the neutrinos. The observational comparison between PKS 0735+178 and other neutrino source candidates, such as TXS 0506+056, PKS 1424+240, and GB6 J1542+6129, shows that all these sources share similar spectral energy distributions, very high radio and $\gamma$-ray powers, and parsec scale jet properties. Moreover, PKS 0735+178, like all the others, is a masquerading BL Lac. We perform comprehensive modelling of the multiwavelength emission from PKS 0735+178 within one-zone lepto-hadronic models considering both internal and external photon fields and estimate the expected accompanying neutrino flux. The most optimistic scenario invokes a jet with luminosity close to the Eddington value and the interactions of $\sim$ PeV protons with an external UV photon field. This scenario predicts $\sim 0.067$ muon and antimuon neutrinos over the observed 3-week flare. Our results are consistent with the detection of one very-high-energy neutrino like IceCube-211208A.

We are in an era of large catalogs and, thus, statistical analysis tools for large data sets, such as machine learning, play a fundamental role. One example of such a survey is the Sloan Moving Object Catalog (MOC), which lists the astrometric and photometric information of all moving objects captured by the Sloan field of view. One great advantage of this telescope is represented by its set of five filters, allowing for taxonomic analysis of asteroids by studying their colors. However, until now, the color variation produced by the change of phase angle of the object has not been taken into account. In this paper, we address this issue by using absolute magnitudes for classification. We aim to produce a new taxonomic classification of asteroids based on their magnitudes that is unaffected by variations caused by the change in phase angle. We selected 9481 asteroids with absolute magnitudes of Hg, Hi and Hz, computed from the Sloan Moving Objects Catalog using the HG12 system. We calculated the absolute colors with them. To perform the taxonomic classification, we applied a unsupervised machine learning algorithm known as fuzzy C-means. This is a useful soft clustering tool for working with {data sets where the different groups are not completely separated and there are regions of overlap between them. We have chosen to work with the four main taxonomic complexes, C, S, X, and V, as they comprise most of the known spectral characteristics. We classified a total of 6329 asteroids with more than 60% probability of belonging to the assigned taxonomic class, with 162 of these objects having been characterized by an ambiguous classification in the past. By analyzing the sample obtained in the plane Semimajor axis versus inclination, we identified 15 new V-type asteroid candidates outside the Vesta family region.

The linear MHD Kelvin-Helmholtz instability (KHI) in an anisotropic plasma concerning the direction of an external magnetic field is examined in detail. For this purpose, the MHD equations are used to describe the motion of plasma as a fluid, which is derived from 16 moments of Boltzmann-Vlasov kinetic equations for collisionless plasma. In addition, the heat flux along the magnetic field is taken into account. The growing rates of KHI are calculated as functions of the anisotropic plasma properties for a shear flow along the magnetic field at supersonic velocities. On the other hand, the quasi-transverse propagation of surface waves between flows with varying velocities is thoroughly examined for both zero-width and finite-width transition layers. In contrast to the tangential discontinuity, it is proved that the limiting breadth of the transition layer constrains the KHI excitation as the wavenumber grows. The instability under investigation could be one of the main ways of dissipation of large-scale low-frequency Alf\'en wave turbulence existing in the solar wind plasma.

A recent analysis of data from the ESA Gaia mission demonstrated that the kinematics of stars in the Milky Way can be modelled without invoking the presence of dark matter whatsoever. Indeed, the higher-than-Keplerian velocities observed in outer stars can be ascribed to the properties of the general relativistic (GR) metric assumed to describe the Galaxy. Here we generalize the concept, and derive the most general exact GR model for a stationary axi-symmetric dynamically cold dust structure. We explicitly show how deviations from the commonly adopted Newtonian dynamics are indeed manifest even at low velocities and low densities. We provide for the first time a detailed description of the frequency shift experienced by photons travelling from any emission site within an external disc galaxy to the observer, relating the outcome of the shift measurement to the gravitational properties of the galaxy. Finally, we devise a novel, groundbreaking observational test potentially able to fully characterize the GR metric under the minimal set of assumptions mentioned above. The proposed experiment exploits the effects non-diagonal GR terms have on the frequency shift of photons, ultimately providing a test to evaluate whether dark matter is actually required by disc galaxy kinematics.

With the prevalence of wide-field, time-domain photometric sky surveys, the number of eclipsing white dwarf systems being discovered is increasing dramatically. An efficient method to follow these up will be key to determining any population trends and finding any particularly interesting examples. We demonstrate that multi-band eclipse photometry of binaries containing a white dwarf and an M~dwarf can be used to determine the masses and temperatures of the white dwarfs to better than 5 per cent. For the M~dwarfs we measure their parameters to a precision of better than 6 per cent with the uncertainty dominated by the intrinsic scatter of the M~dwarf mass-radius relationship. This precision is better than what can typically be achieved with low-resolution spectroscopy. The nature of this method means that it will be applicable to LSST data in the future, enabling direct characterisation without follow-up spectroscopy. Additionally, we characterise three new post-common-envelope binaries from their eclipse photometry, finding two systems containing hot helium-core white dwarfs with low-mass companions (one near the brown dwarf transition regime) and a possible detached cataclysmic variable at the lower edge of the period gap.

We extend the analysis of Bok et al. (2020) in which the HI content of isolated galaxies from the AMIGA sample and selected paired galaxies from ALFALFA were examined as a potential driver of galaxy location on the WISE mid-infrared SFR-Mstar sequence. By further characterizing the isolated and pair galaxy samples, i.e. in terms of optical galaxy morphology, a more detailed and quantitative description of local galaxy environment by way of the local number density (eta) and tidal strength (Q) parameters, star formation efficiency (SFE(HI)), and HI integrated profile asymmetries, we present plausible pathways for the broadening of the pair sample HI deficiency distribution towards both high and low deficiencies compared to the narrower isolated galaxy sample distribution (i.e. sigma_pairs = 0.34 versus sigma_AMIGA = 0.28). We associate the gas-rich tail of the pair deficiency distribution with the highest Q values, large profile asymmetries, and low SFEs. From this we infer that merger activity is enhancing gas supplies, as well as disrupting the efficiency of star formation, via strong gravitational torques. The gas-poor wing of the deficiency distribution appears to be populated with galaxies in denser environments (with larger eta values on average), more akin to groups. Despite our gas-rich selection criterion, there is a small population of early-type galaxies in the pair sample, which primarily fall in the positive deficiency wing of the distribution. These results suggest that a combination of a denser galaxy environment, early-type morphology, and higher stellar mass, is contributing to the broadening of the deficiency distribution towards larger deficiencies.

Global hydrodynamic simulations of internal solar dynamics have focused on replicating the conditions for solar-like differential rotation and meridional circulation using the results of helioseismic inversions as a constraint. Inferences of meridional circulation, however, have provided controversial results showing the possibility of one, two, or multiple cells along the radius. To resolve this controversy and develop a more robust understanding of global flow regimes in the solar interior, we apply a "forward-modeling" approach to the analysis of helioseismic signatures of meridional circulation profiles obtained from numerical simulations. We employ the global acoustic modeling code GALE to simulate the propagation of acoustic waves through regimes of mean mass flows generated by global hydrodynamic and magnetohydrodynamic models: EULAG, the Pencil Code, and the Rayleigh code. These models are used to create synthetic dopplergram data products, used as inputs for local time-distance helioseismology techniques. Helioseismic travel-time signals from solutions obtained through global numerical simulations are compared directly with inferences from solar observations, in order to set additional constraints on global model parameters in a direct way. We show that even though these models are able to replicate solar-like differential rotation, the resulting rotationally-constrained convection develops a multi-cell global meridional circulation profile that is measurably inconsistent with local time-distance inferences of solar observations. However, we find that the development of rotationally-unconstrained convection close to the model surface is able to maintain solar-like differential rotation, while having a significant impact on the helioseismic travel-time signal, replicating solar observations within one standard deviation of the error due to noise.

Context. Active galactic nuclei (AGN) of galaxies play an important role in the life and evolution of galaxies due to the impact they exert on certain properties and the evolutionary path of galaxies. It is well known that infrared (IR) emission is useful for selecting galaxies with AGNs, although it has been observed that there is contamination by star-forming galaxies. Aims. In this work we investigate galaxy properties hosting AGNs identified at mid and near-IR wavelengths. The sample of AGNs selected at IR wavelengths was confirmed using optical spectroscopy and X-ray photometry. We study the near-UV, optical, near and mid-IR (MIR) properties, as well as [O III] {\lambda}5007 luminosity, black hole mass and morphology properties of optical and IR colour selected AGNs. Methods. We selected AGN candidates using two mid-IR colour selection techniques, a power-law emission method and a combination of mid and near-IR selection techniques. We confirm the AGN selection with two line diagnostic diagrams that use the ratio [O III]/H\b{eta} and the emission line width {\sigma} [O III] (kinematics-excitation diagram, KEx) and the host galaxy stellar mass (mass-excitation diagram, MEx), as well as X-ray photometry. Results. According to the diagnostic diagrams, the methods with the greatest success in selecting AGNs are those that use a combination of a mid and near-IR selection technique and a power-law emission. The method that use a combination of mid and near-IR observation selects a large number of AGNs, and is reasonably efficient in both the success rate (61%) and total number of AGN recovered. We also find that the KEx method presents contamination of SF galaxies within the AGN selection box. According to morphological studies based on the S\'ersic index, AGN samples have higher percentages of galaxy morphologies with bulge+disk components compared to galaxies without AGNs.

After large galaxies merge, their central supermassive black holes are expected to form binary systems whose orbital motion generates a gravitational wave background (GWB) at nanohertz frequencies. Searches for this background utilize pulsar timing arrays, which perform long-term monitoring of millisecond pulsars (MSPs) at radio wavelengths. We use 12.5 years of Fermi Large Area Telescope data to form a gamma-ray pulsar timing array. Results from 35 bright gamma-ray pulsars place a 95\% credible limit on the GWB characteristic strain of $1.0\times10^{-14}$ at 1 yr$^{-1}$, which scales as the observing time span $t_{\mathrm{obs}}^{-13/6}$. This direct measurement provides an independent probe of the GWB while offering a check on radio noise models.

Given that the reionization history of cosmic hydrogen is yet to be stringently constrained, it is worth checking the prospects of doing so using physically motivated models and available observational data. For this purpose, we use an extended version of the explicitly photon-conserving semi-numerical model of reionization, $\texttt{SCRIPT}$, which also includes thermal evolution of the intergalactic medium (IGM). The model incorporates the effects of inhomogeneous recombination and radiative feedback self-consistently and is characterized by five free parameters (two for the redshift-dependent ionization efficiency, two for the ionizing escape fraction, and another for reionization temperature increment). We constrain these free parameters by simultaneously matching with various observational probes, e.g., estimates of the ionized hydrogen fraction, the CMB scattering optical depth and the galaxy UV luminosity function. In addition, we include the low-density IGM temperature measurements obtained from Lyman-$\alpha$ absorption spectra at $z \sim 5.5$, a probe not commonly used for Bayesian analysis of reionization parameters. We find that the interplay of the various data sets, particularly inclusion of the temperature data, leads to tightening of the parameter constraints. Our default models prefer a late end of reionization (at $z \lesssim 6$), in agreement with other recent studies. We can also derive constraints on the duration of reionization, $\Delta z=1.81^{+0.51}_{-0.67}$ and the midpoint of reionization, $z_{\mathrm{mid}}=7.0^{+0.30}_{-0.40}$. The constraints can be further tightened by including other available and upcoming data sets.

GRS 1915+105 can show type-C quasi-periodic oscillations (QPOs) in the power density spectrum. A high-frequency QPO (HFQPO) at 67 Hz has been observed in this source, albeit less often than the type-C QPOs. Besides these features, GRS 1915+105 sometimes shows a broad bump in the power spectrum at around 30-150 Hz. We study the power spectra of GRS 1915+105 with the Rossi X-ray Timing Explorer when the source was in the $\chi$ class. We find that the rms amplitude of the bump depends strongly upon both the frequency of the type-C QPO and the hardness ratio, and is correlated with the corona temperature and anti-correlated with the radio flux at 15 GHz. The characteristic frequency of the bump is better correlated with a combination of the frequency of the type-C QPO and the hardness ratio than with the frequency of the type-C QPO alone. The rms amplitude of the bump generally increases with energy from ~1-2% at ~3 keV to ~10-15% at ~30 keV. We suggest that the bump and the high-frequency QPO may be the same variability component but the properties of the corona affect the coherence of this variability, leading either to a HFQPO when the spectrum is in the relatively soft $\gamma$ class, or to a bump when the spectrum is in the hard $\chi$ class. Finally, we discuss the anti-correlation between the rms amplitude of the bump and the radio flux in the context of the relation between the corona and the jet.

We study the radio power of the core and its relation to the optical properties of the host galaxy in samples of high excitation (HERG) and low excitation (LERG) Fanaroff-Riley type II (FRII) radio galaxies. The radio galaxy sample is divided into two groups of core/non-core FRII, based on the existence of strong, weak or lack of single radio core component. We show that FRII LERGs with radio emission of the core have significantly higher [O III] line luminosities compared to the non-core LERG FRIIs. There is no significant difference between the hosts of the core and non-core FRIIs of LERG type in galaxy sizes, concentration indices, star formation rates, 4000-\AA\ break strengths, colours, black hole masses and black hole to stellar masses. We show that the results are not biased by the stellar masses, redshifts and angular sizes of the radio galaxies. We argue that the detection of higher [O III] luminosities in the core FRIIs may indicate the presence of higher amounts of gas, very close to the AGN nucleus in the core FRIIs compared to the non-core FRIIs or may result from the interaction of the radio jets with this gas. The core and non-core FRIIs of the HERG type show no significant differences perhaps due to our small sample size. The effect of relativistic beaming on the radio luminosities and the contribution of restating AGN activity have also been considered.

Utilizing the data and tools provided through the NASA Ames PAH IR Spectroscopic Database (PAHdb), we study the PAH component of over 900 Spitzer-IRS galaxy spectra. Employing a database-fitting approach, the average PAH size, the PAH size distribution, and PAH ionization fraction are deduced. In turn, we examine their connection with the properties of the host galaxy. We found that PAH population within galaxies consists of middle-sized PAHs with an average number of carbon atoms of $\bar{N_{C}}$ = 55, and a charge state distribution of $\sim$40% ionized - 60% neutral. We describe a correlation between the 6.2/11.2 $\mu$m PAH ratio with the ionization parameter ($\gamma\equiv(G_{0}/n_{\rm e})(T_{\rm gas} / 1\ \mathrm{K})^{0.5}$), a moderate correlation between the 8.6/11.2 $\mu$m PAH ratio and specific star-formation rate, and a weak anti-correlation between $\gamma$ and M$_{*}$. From the PAHdb decomposition we provide estimates for the 3.3 $\mu$m PAH band, not covered by Spitzer observations, and establish a correlation between the 3.3/11.2 $\mu$m PAH ratio with N$_{\mathrm{C}}$. We further deliver a library of mid-IR PAH template spectra parameterized on PAH size and ionization fraction, which can be used in galaxy spectral energy distribution fitting codes for the modeling of the mid-IR PAH emission component in galaxies.

We identify a new cosmological signal, the Doppler-boosted Cosmic Infrared Background (DB-CIB), arising from the peculiar motion of the galaxies whose thermal dust emission source the cosmic infrared background (CIB). This new observable is an independent probe of the cosmic velocity field, highly analogous to the well-known kinematic Sunyaev-Zel'dovich (kSZ) effect. Interestingly, DB-CIB does not suffer from the 'kSZ optical depth degeneracy', making it immune from the complex astrophysics of galaxy formation. We forecast that the DB-CIB effect is detectable in the cross-correlation of CCAT-Prime and DESI-like experiments. We show that it also acts as a new CMB foreground which can bias future kSZ cross-correlations, if not properly accounted for.

Using a suite of $N$-body simulations we study the angular clustering of galaxies, halos, and dark matter in $\mathrm{\Lambda \text{CDM}}$ and Modified Gravity (MG) scenarios. We consider two general categories of such MG models, one is the $f(R)$ gravity, and the other is the normal branch of the Dvali-Gabadadze-Porrati brane world (nDGP). To measure angular clustering we construct a set of observer-frame lightcones and resulting mock sky catalogs. We focus on the area-averaged angular correlation functions, $W_J$, and the associated reduced cumulants, $S_J\equiv W_J/W_2^{(J-1)}$, and robustly measure them up to the 9th order using counts-in-cells (CIC). We find that $0.15 < z < 0.3$ is the optimal redshift range to maximize the MG signal in our lightcones. Analyzing various scales for the two types of statistics, we identify up to 20\% relative departures in MG measurements from general relativity (GR), with varying signal significance. For the case of halos and galaxies, we find that $3$rd order statistics offer the most sensitive probe of the different structure formation scenarios, with both $W_3$ and the reduced skewness, $S_3$, reaching from $2\sigma$ to $4\sigma$ significance at angular scales $\theta \sim 0.13 ^\circ$. The MG clustering of the smooth dark matter field is characterized by even stronger deviations ($\stackrel{>}{{}_\sim} 5\sigma$) from GR, albeit at a bit smaller scales of $\theta\sim0.08^\circ$, where baryonic physics is already important. Finally, we stress out that our mock halo and galaxy catalogs are characterized by rather low surface number densities when compared to existing and forthcoming state-of-the-art photometric surveys. This opens up exciting potential for testing GR and MG using angular clustering in future applications, with even higher precision and significance than reported here.

Sterile neutrinos can be produced in the early universe via interactions with their active counterparts. For small active-sterile mixing angles, thermal equilibrium with the standard model plasma is not reached and sterile neutrinos are only produced via flavor oscillations. We study in detail this regime, taking into account matter potentials and decoherence effects caused by elastic scatterings with the plasma. We find that resonant oscillations occurring at temperatures $T\lesssim 10\,\mathrm{GeV}$ lead to a significant enhancement of the sterile neutrino production rate. Taking this into account, we improve constraints on the active-sterile mixing from Big Bang nucleosynthesis and the cosmic microwave background, excluding mixing angles down to $\theta_s\sim 10^{-10}-10^{-16}$ for sterile neutrino masses in the $10\,\mathrm{MeV}$ to $10\,\mathrm{GeV}$ range. We observe that if sterile neutrinos predominantly decay into metastable hidden sector particles, this process provides a novel dark matter production mechanism, consistent with the sterile neutrino origin of light neutrino masses via the seesaw mechanism.

Neutrino interactions with protons and neutrons probe their deep structure and may reveal new physics. The higher the neutrino energy, the sharper the probe. So far, the neutrino-nucleon ($\nu N$) cross section is known across neutrino energies from a few hundred MeV to a few PeV. Soon, ultra-high-energy (UHE) cosmic neutrinos, with energies above 100 PeV, could take us farther. So far, they have evaded discovery, but upcoming UHE neutrino telescopes endeavor to find them. We present the first detailed measurement forecasts of the UHE $\nu N$ cross section, geared to IceCube-Gen2, one of the leading detectors under planning. We use state-of-the-art ingredients in every stage of our forecasts: in the UHE neutrino flux predictions, the neutrino propagation inside Earth, the emission of neutrino-induced radio signals in the detector, their propagation and detection, and the treatment of backgrounds. After 10 years, if at least a few tens of UHE neutrino-induced events are detected, IceCube-Gen2 could measure the $\nu N$ cross section at center-of-mass energies of $\sqrt{s} \approx 10-100$ TeV for the first time, with a precision comparable to that of its theory prediction.

Measures of discrepancy between probability distributions (statistical distance) are widely used in the fields of artificial intelligence and machine learning. We describe how certain measures of statistical distance can be implemented as numerical diagnostics for simulations involving charged-particle beams. Related measures of statistical dependence are also described. The resulting diagnostics provide sensitive measures of dynamical processes important for beams in nonlinear or high-intensity systems, which are otherwise difficult to characterize. The focus is on kernel-based methods such as Maximum Mean Discrepancy, which have a well-developed mathematical foundation and reasonable computational complexity. Several benchmark problems and examples involving intense beams are discussed. While the focus is on charged-particle beams, these methods may also be applied to other many-body systems such as plasmas or gravitational systems.

Two experiments from the Fermilab, E989 and CDF II, have reported two anomalies for muon anomalous magnetic moment ($g$-2) and $W$-boson mass that may indicate the new physics at the low energy scale. Here we examine the possibility of a common origin of these two anomalies in the Next-to-Minimal Supersymmetric Standard Model. Considering various experimental and astrophysical constraints such as the Higgs mass, collider data, B-physics, dark matter relic density and direct detection experiments, we find that a neutralino in the mass range of $\sim 160-270$ GeV is a viable solution. Moreover, the favored parameter region can be effectively probed by the ongoing direct detection experiments like LZ, PandaX-4T and XENON-nT. The velocity averaged annihilation cross section of the dark matter particles, however, is suppressed.

This paper highlights a methodological approach designed to enhance the search for extraterrestrial intelligence (SETI) by hypothesizing that a transmission technosignature would likely have two features: 1) be wideband in the microwave or higher frequency range that originates from a hub within a supposed ET interplanetary navigation/communication (nav/comm) network, and 2) contain x-ray pulsar-based navigation (XNAV) metadata. Potential contributions to the field include improved accuracy in finding transmission technosignatures and other technosignatures in the electromagnetic spectrum, a common standard in reaching a Schelling Point (a mutual realization of how we and ETs can find each other), and operationalizing models such as the Drake Equation.

We perform a three dimensional lattice simulation of the electroweak symmetry breaking process through a two-step phase transition, where one of the two steps is a first order phase transition. Our results show that: 1) when the electroweak symmetry breaking is driven by the beyond Standard Model sector around $\sim \mathcal{O}(10^{2-3})$ GeV, the gravitational wave spectra produced from the phase transitions are of broken power-law double-peak shapes; 2) when the electroweak symmetry breaking is induced by a first-order phase transition of a high-scale global U(1) theory, cosmic strings can form and then disappear through particle radiation, and the yielded gravitational wave spectra are of plateau shapes. The two scenarios can be distinguished through probing gravitational wave spectra. Our study suggests that the stochastic gravitational waves provide an alternative way to probe the beyond Standard Model sector relevant to the electroweak symmetry breaking pattern in the early Universe.

Finsler geometry is an important extension of Riemann geometry, in which to each point of the spacetime manifold an arbitrary internal variable is associated. Interesting Finsler geometries, with many physical applications, are the Randers and Kropina type geometries, respectively. A subclass of Finsler geometries is represented by the osculating Finsler spaces, in which the internal variable is a function of the base manifold coordinates only. In an osculating Finsler geometry one introduces the Barthel connection, which has the remarkable property that it is the Levi-Civita connection of a Riemannian metric. In the present work we consider the gravitational and cosmological implications of a Barthel-Kropina type geometry. We assume that in this geometry the Ricci type curvatures are related to the matter energy-momentum tensor by the standard Einstein equations. The generalized Friedmann equations in the Barthel-Kropina geometry are obtained by considering that the background Riemannian metric is of Friedmann-Lemaitre-Robertson-Walker type. The matter energy balance equation is also derived. The cosmological properties of the model are investigated in detail, and it is shown that the model admits a de Sitter type solution, and that an effective dark energy component can also be generated. Several cosmological solutions are also obtained by numerically integrating the generalized Friedmann equations. A comparison of two specific classes of models with the observational data and with the standard $\Lambda$CDM model is also performed, and it turns out that the Barthel-Kropina type models give a satisfactory description of the observations.

We present a phenomenological and numerical study of strong Alfv\'enic turbulence in a magnetically dominated collisionless relativistic plasma with a strong background magnetic field. In contrast with the non-relativistic case, the energy in such turbulence is contained in magnetic and electric fluctuations. We argue that such turbulence is analogous to turbulence in a strongly magnetized non-relativistic plasma in the regime of broken quasi-neutrality. Our 2D particle-in-cell numerical simulations of turbulence in a relativistic pair plasma find that the spectrum of the total energy has the scaling $k^{-3/2}$, while the difference between the magnetic and electric energies, the so-called residual energy, has the scaling $k^{-2.4}$. The electric and magnetic fluctuations at scale $\ell$ exhibit dynamic alignment with the alignment-angle scaling close to $\cos\phi_\ell\propto \ell^{1/4}$. At scales smaller than the (relativistic) plasma inertial scale, the energy spectrum of relativistic inertial Alfv\'en turbulence steepens to $k^{-3.5}$.

The isomerization of hydrogen cyanide to hydrogen isocyanide on icy grain surfaces is investigated by an accurate composite method (jun-Cheap) rooted in the coupled cluster ansatz and by density functional approaches. After benchmarking density functional predictions of both geometries and reaction energies against jun-Cheap results for the relatively small model system HCN -- (H2O)2 the best performing DFT methods are selected. A large cluster containing 20 water molecules is then employed within a QM/QM$'$ approach to include a realistic environment mimicking the surface of icy grains. Our results indicate that four water molecules are directly involved in a proton relay mechanism, which strongly reduces the activation energy with respect to the direct hydrogen transfer occurring in the isolated molecule. Further extension of the size of the cluster up to 192 water molecules in the framework of a three-layer QM/QM'/MM model has a negligible effect on the energy barrier ruling the isomerization. Computation of reaction rates by transition state theory indicates that on icy surfaces the isomerization of HNC to HCN could occur quite easily even at low temperatures thanks to the reduced activation energy that can be effectively overcome by tunneling.

Small, expendable drop probes are an attractive method for making measurements in the lower atmosphere of Venus, \,augmenting the capabilities of orbiters or aerial platforms that must remain in the benign temperature region above 50 km altitude. However, probe miniaturization is impeded by the need to provide thermal and pressure protection for conventional payloads. This paper determines the minimum mass limits for an insulated pressure vessel probe that operates all the way to the Venusian surface. Scaling laws for the probe performance and mass of major system components are explicitly derived using a simple model that captures the relevant physics. Streamlining the probe is found to be a highly effective strategy for lowering the system mass, but it also reduces the time available for data collection and transmission. Tradeoffs, guidelines and design charts are presented for an array of miniaturized probes. Total system masses on the order of 5 kilograms are plausible with streamlined probes if the desired science measurements can be performed faster than a standard Venus descent timeline.

Recently, the CDF collaboration has reported the precise measurement of the W boson mass, $M_W = 80433.5\pm 9.4 \,$MeV, based on $8.8$ fb$^{-1}$ of $\sqrt{s}=1.96$ TeV $p\bar{p}$ collision data from the CDF II detector at the Fermilab Tevatron. This is about $7\sigma$ away from the Standard Model prediction, $M_{W}^{\rm SM}=80357 \pm 6 \,$MeV. Such a large discrepancy may be partially due to exotic particles that radiatively alter the relation between the W and Z boson masses. In this Letter, we study singlet extensions of the Standard Model focusing on the shift of the W boson mass. In the minimal extension with a real singlet field, using the bounds from the electroweak oblique parameters, B meson decays, LEP, and LHC, we find that the W boson mass shift is at most a few MeV, and therefore it does not alleviate the tension between the CDF II result and the SM prediction. We then examine how much various bounds are relaxed when the singlet is allowed to decay invisibly and find that the increase of the W boson mass does not exceed $5$ MeV due to the bound from the Higgs signal strength. We also discuss phenomenological and cosmological implications of the singlet extensions such as the muon $g-2$ anomaly, axion/hidden photon dark matter, and self-interacting dark radiation as a possible alleviation of the Hubble tension.

In this paper we have considered an interacting model of dark energy and have looked into the evolution of the dark sectors. By solving the perturbation equations numerically, we have studied the imprints on the growth of matter as well as dark energy fluctuations. It has been found that for higher rate of interaction strength for the coupling term, visible imprints on the dark energy density fluctuations are observed at the early epochs of evolution.

We calculate one-loop correction to the two-point functions of curvature perturbation in single-field inflation generated by cubic self-interaction. Incorporating the observed red-tilted spectrum of curvature perturbation, the relevant one-loop correction takes a finite value and inversely proportional to the spectral tilt. Requiring one-loop correction to be much smaller than the tree-level contribution leads to an upper bound on primordial non-Gaussianity. While observationally allowed region of non-Gaussian parameter space is found to be entirely included by the region where one-loop correction is smaller than the tree-level contribution, if we tighten the theoretical requirement to require that the former should be smaller than, say, 10% or 1% of the latter, significant parts of the observationally allowed region are ruled out. If future observations conclude non-Gaussianity falls in such a region, then it would be important to incorporate higher-order corrections to the spectrum in order to achieve precise cosmology. In some extreme cases, it might indicate breakdown of the cosmological perturbation theory in the context of single-field inflation.