Papers for Tuesday, Jul 06 2021

Young transiting exoplanets (< 100 Myr) provide crucial insight into atmospheric evolution via photoevaporation. However, transmission spectroscopy measurements to determine atmospheric composition and mass loss are challenging due to the activity and prominent stellar disk inhomogeneities present on young stars. We observed a full transit of V1298 Tau c, a 23 Myr, 5.59$R_\oplus$ planet orbiting a young K0-K1.5 solar analogue with GRACES on Gemini-North. We were able to measure the Doppler tomographic signal of V1298 Tau c using the Ca II infrared triplet (IRT) and find a projected obliquity of $\lambda = 5^\circ \pm 15^\circ$. The tomographic signal is only seen in the chromospherically driven core of the Ca II IRT, which may be the result of star-planet interactions. Additionally, we find that excess absorption of the H-alpha line decreases smoothly during the transit. While this could be a tentative detection of hot gas escaping the planet, we find this variation is consistent with similar timescale observations of other young stars that lack transiting planets over similar timescales. We show this variation can also be explained by the presence of starspots with surrounding facular regions. More observations both in- and out-of the transits of V1298 Tau c are required to determine the nature of the Ca II IRT and H-alpha line variations.

To indirectly study the internal structure of giant clumps in main sequence galaxies at $z \sim 1-3$, we target very turbulent and gas-rich local analogues from the DYNAMO sample with the Hubble Space Telescope, over a wavelength range of $\sim 200-480$ nm. We present a catalog of 58 clumps identified in six DYNAMO galaxies, including the WFC3/UVIS F225W, F336W, and F467M photometry where the ($225-336$) and ($336-467$) colours are sensitive to extinction and stellar population age respectively. We measure the internal colour gradients of clumps themselves to study their age and extinction properties. We find a marked colour trend within individual clumps, where the resolved colour distributions show that clumps generally have bluer ($336-467$) colours (denoting very young ages) in their centers than at their edges, with little variation in the ($225-336$) colour associated with extinction. Furthermore, we find that clumps whose colours suggest they are older, are preferentially located closer toward the centers of their galaxies, and we find no young clumps at small galactocentric distances. Both results are consistent with simulations of high-redshift star forming systems that show clumps form via violent disk instability, and through dynamic processes migrate to the centers of their galaxies to contribute to bulge growth on timescales of a few 100 Myr, while continually forming stars in their centers. When we compare the DYNAMO clumps to those in these simulations, we find the best agreement with the long-lived clumps.

Infrared excesses around white dwarf stars indicate the presence of various astrophysical objects of interest, including companions and debris disks. In this second paper of a series, we present follow-up observations of infrared excess candidates from Gaia and unWISE discussed in the first paper, Paper I. We report space-based infrared photometry at 3.6 and 4.5 micron for 174 white dwarfs from the Spitzer Space Telescope and ground-based near-infrared J, H, and K photometry of 235 white dwarfs from Gemini Observatory with significant overlap between Spitzer and Gemini observations. This data is used to confirm or rule-out the observed unWISE infrared excess. From the unWISE-selected candidate sample, the most promising infrared excess sample comes from both colour and flux excess, which has a Spitzer confirmation rate of 95%. We also discuss a method to distinguish infrared excess caused by stellar or sub-stellar companions from potential dust disks. In total, we confirm the infrared excess around 61 white dwarfs, 10 of which are likely to be stellar companions. The remaining 51 bright white dwarf with infrared excess beyond two microns has the potential to double the known sample of white dwarfs with dusty exoplanetary debris disks. Follow-up high-resolution spectroscopic studies of a fraction of confirmed excess white dwarfs in this sample have discovered emission from gaseous dust disks. Additional investigations will be able to expand the parameter space from which dust disks around white dwarfs are found.

We report multi-epoch radial velocities, rotational velocities, and atmospheric parameters for 37 T-type brown dwarfs observed with Keck/NIRSPEC. Using a Markov Chain Monte Carlo forward-modeling method, we achieve median precisions of 0.5 km s$^{-1}$ and 0.9 km s$^{-1}$ for radial and rotational velocities, respectively. All of the T dwarfs in our sample are thin disk brown dwarfs. We confirm previously reported moving group associations for four T dwarfs. However, the lack of spectral indicators of youth in two of these sources suggests that these are chance alignments. We confirm two previously un-resolved binary candidates, the T0+T4.5 2MASS J11061197+2754225 and the L7+T3.5 2MASS J21265916+7617440, with orbital periods of 4 yr and 12 yr, respectively. We find a kinematic age of 3.5$\pm$0.3 Gyr for local T dwarfs, consistent with nearby late-M dwarfs (4.1$\pm$0.3 Gyr). Removal of thick disk L dwarfs in the local ultracool dwarf sample gives a similar age for L dwarfs (4.2$\pm$0.3 Gyr), largely resolving the local L dwarf age anomaly. The kinematic ages of local late-M, L, and T dwarfs can be accurately reproduced with population simulations incorporating standard assumptions of the mass function, star formation rate, and brown dwarf evolutionary models. A kinematic dispersion break is found at the L4$-$L6 subtypes, likely reflecting the terminus of the stellar Main Sequence. We provide a compilation of precise radial velocities for 172 late-M, L, and T dwarfs within $\sim$20 pc of the Sun.

Understanding the impact of halo properties beyond halo mass on the clustering of galaxies (namely galaxy assembly bias) remains a challenge for contemporary models of galaxy clustering. We explore the use of machine learning to predict the halo occupations and recover galaxy clustering and assembly bias in a semi-analytic galaxy formation model. For stellar-mass selected samples, we train a Random Forest algorithm on the number of central and satellite galaxies in each dark matter halo. With the predicted occupations, we create mock galaxy catalogues and measure the clustering and assembly bias. Using a range of halo and environment properties, we find that the machine learning predictions of the occupancy variations with secondary properties, galaxy clustering and assembly bias are all in excellent agreement with those of our target galaxy formation model. Internal halo properties are most important for the central galaxies prediction, while environment plays a critical role for the satellites. Our machine learning models are all provided in a usable format. We demonstrate that machine learning is a powerful tool for modelling the galaxy-halo connection, and can be used to create realistic mock galaxy catalogues which accurately recover the expected occupancy variations, galaxy clustering and galaxy assembly bias, imperative for cosmological analyses of upcoming surveys.

We apply novel, recently developed plasma ray-tracing techniques to model the propagation of radio photons produced by axion dark matter in neutron star magnetospheres and combine this with both archival and new data for the galactic centre magnetar PSR J1745-2900. The emission direction to the observer and the magnetic orientation are not constrained for this object leading to parametric uncertainty. Our analysis reveals that ray-tracing greatly reduces the signal sensitivity to this uncertainty, contrary to previous calculations where there was no emission at all in some directions. Based on a Goldreich-Julian model for the magnetosphere and a Navarro-Frank-White model for axion density in the galactic centre, we obtain the most robust limits on the axion-photon coupling, to date. These are comparable to those from the CAST solar axion experiment in the mass range $\sim 4.2-60\,\mu{\rm eV}$. If the dark matter density is larger, as might predicted by a "spike" model, the limits could be much stronger. The dark matter density in the region of the galactic centre is now the biggest uncertainty in these calculations.

Dark matter halos that reach the HI-cooling mass without prior star formation or external metal pollution represent potential sites for the formation of small Population III galaxies at high redshifts. Such objects are expected to attain total stellar masses of at most $10^6$ solar masses and will therefore typically be extremely faint. Gravitational lensing may in rare cases boost their fluxes to detectable levels, but to find even a small number of such objects requires very large sky areas to be surveyed. Because of this, a small, wide-field telescope can in principle offer better detection prospects than a large telescope with a smaller field of view. Here, we derive the Pop III galaxy properties - in terms of comoving number density, stellar initial mass function and total stellar mass - required to allow gravitational lensing to lift such objects at redshift z = 5-16 above the detection thresholds of blind surveys carried out with the James Webb space telescope (JWST), the Roman space telescope (RST) or Euclid. We find that the prospects for photometric detections of Pop III galaxies are promising, and that they are better for RST than for JWST and Euclid. However, the Pop III galaxies favoured by current simulations have number densities too low to allow spectroscopic detections based on the strength of the HeII1640 emission line in any of the considered surveys unless very high star formation efficiencies (10 per cent) are envoked.

We present the serendipitous discovery of a galaxy group in the XMM-LSS field with MIGHTEE Early Science observations. Twenty galaxies are detected in HI in this $z\sim0.044$ group, with a $3\sigma$ column density sensitivity of $N_\mathrm{HI} = 1.6\times10^{20}\,\mathrm{cm}^{-2}$. This group has not been previously identified, despite residing in a well-studied extragalactic legacy field. We present spatially-resolved HI total intensity and velocity maps for each of the objects, which reveal environmental influence through disturbed morphologies. The group has a dynamical mass of $\log_{10}(M_\mathrm{dyn}/\mathrm{M}_\odot) = 12.32$, and is unusually gas-rich, with an HI-to-stellar mass ratio of $\log_{10}(f_\mathrm{HI}^\mathrm{*}) = -0.2$, which is 0.7 dex greater than expected. The group's high HI content, spatial, velocity, and identified galaxy type distributions strongly suggest that it is in the early stages of its assembly. The discovery of this galaxy group is an example of the importance of mapping spatially-resolved HI in a wide range of environments, including galaxy groups. This scientific goal has been dramatically enhanced by the high sensitivity, large field-of-view, and wide instantaneous bandwidth of the MeerKAT telescope.

The sunspot umbra harbors the coolest plasma on the solar surface due to the presence of strong magnetic fields. The routinely used atomic lines to observe the photosphere have weak signals in the umbra and are often swamped by molecular lines. This makes it harder to infer the properties of the umbra, especially in the darkest regions. The lines of the Ti I multiplet at 2.2 $\mu$m are formed mainly at temperatures $\le\!4500$ K and are not known to be affected by molecular blends in sunspots. Since the first systematic observations in the 1990's, these lines have been seldom observed due to the instrumental challenges involved at these longer wavelengths. We revisit these lines and investigate their formation in different solar features. We synthesize the Ti I multiplet using a snapshot from 3D MHD simulation of a sunspot and explore the properties of two of its lines in comparison with two commonly used iron lines at 630.25 nm and $1.5648\,\mu$m. We find that the Ti I lines have stronger signals than the Fe I lines in both intensity and polarization in the sunspot umbra and in penumbral spines. They have little to no signal in the penumbral filaments and the quiet Sun, at $\mu=1$. Their strong and well-split profiles in the dark umbra are less affected by stray light. Consequently, inside the sunspot it is easier to invert these lines and to infer the atmospheric properties, compared to the iron lines. The Cryo-NIRSP instrument at the DKIST will provide the first ever high resolution observations in this wavelength range. In this preparatory study, we demonstrate the unique temperature and magnetic sensitivities of the Ti multiplet, by probing the Sun's coolest regions which are not favourable for the formation of other commonly used spectral lines. We thus expect such observations to advance our understanding of sunspot properties.

Hyperactive comets have high water production rates, with inferred sublimation areas of order the surface area of the nucleus. Comets 46P/Wirtanen and 103P/Hartley 2 are two examples of this cometary class. Based on observations of comet Hartley 2 by the Deep Impact spacecraft, hyperactivity appears to be caused by the ejection of water-ice grains and/or water-ice rich chunks of nucleus into the coma. These materials increase the sublimating surface area, and yield high water production rates. The historic close approach of comet Wirtanen to Earth in 2018 afforded an opportunity to test Hartley 2 style hyperactivity in a second Jupiter-family comet. We present high spatial resolution, near-infrared spectroscopy of the inner coma of Wirtanen. No evidence for the 1.5- or 2.0-$\mu$m water-ice absorption bands is found in six 0.8-2.5 $\mu$m spectra taken around perihelion and closest approach to Earth. In addition, the strong 3.0-$\mu$m water-ice absorption band is absent in a 2.0-5.3 $\mu$m spectrum taken near perihelion. Using spectroscopic and sublimation lifetime models we set constraints on the physical properties of the ice grains in the coma, assuming they are responsible for the comet's hyperactivity. We rule out pure water-ice grains of any size, given their long lifetime. Instead, the hyperactivity of the nucleus and lack of water-ice absorption features in our spectra can be explained either by icy grains on the order of 1 $\mu$m in size with a small amount of low albedo dust (greater than 0.5% by volume), or large chunks containing significant amounts of water ice.

Relativistic magnetic reconnection is a powerful agent through which magnetic energy can be tapped in astrophysics, energizing particles that then produce observed radiation. In some systems, the highest energy photons come from particles Comptonizing an ambient radiation bath supplied by an external source. If the emitting particle energies are high enough, this inverse Compton (IC) scattering enters the Klein-Nishina regime, which differs from the low-energy Thomson IC limit in two significant ways. First, radiative losses become inherently discrete, with particles delivering an order-unity fraction of their energies to single photons. Second, Comptonized photons may pair-produce with the ambient radiation, opening up another channel for radiative feedback on magnetic reconnection. We analytically study externally illuminated highly magnetized reconnecting systems for which both of these effects are important. We identify a universal (initial magnetization-independent) quasi-steady state in which gamma-rays emitted from the reconnection layer are absorbed in the upstream region, and the resulting hot pairs dominate the energy density of the inflow plasma. However, a true pair cascade is unlikely, and the number density of created pairs remains subdominant to that of the original plasma for a wide parameter range. Future particle-in-cell simulation studies may test various aspects. Pair-regulated Klein-Nishina reconnection may explain steep spectra (quiescent and flaring) from flat-spectrum radio quasars and black hole accretion disc coronae.

This is the second part of a thorough investigation of the redshift-space effects that affect void properties and the impact they have on cosmological tests. Here, we focus on the void-galaxy cross-correlation function, specifically, on the projected versions that we developed in a previous work. The pillar of the analysis is the one-to-one relationship between real and redshift-space voids above the shot-noise level identified with a spherical void finder. Under this mapping, void properties are affected by three effects: (i) a systematic expansion as a consequence of the distortions induced by galaxy dynamics, (ii) the Alcock-Paczynski volume effect, which manifests as an overall expansion or contraction depending on the fiducial cosmology, and (iii) a systematic off-centring along the line of sight as a consequence of the distortions induced by void dynamics. We found that correlations are also affected by an additional source of distortions: the ellipticity of voids. This is the first time that distortions due to the off-centring and ellipticity effects are detected and quantified. With a simplified test, we verified that the Gaussian streaming model is still robust provided all these effects are taken into account, laying the foundations for improvements in current models in order to obtain unbiased cosmological constraints from spectroscopic surveys. Besides this practical importance, this analysis also encodes key information about the structure and dynamics of the Universe at the largest scales. Furthermore, some of the effects constitute cosmological probes by themselves, as is the case of the void ellipticity.

Solar coronal rain is classified generally into two categories: flare-driven and quiescent coronal rain. The latter is observed to form along both closed and open magnetic field structures. Recently, we proposed that some of the quiescent coronal rain events, detected in the transition region and chromospheric diagnostics, along loop-like paths could be explained by the formation mechanism for quiescent coronal rain facilitated by interchange magnetic reconnection between open and closed field lines. In this study, we revisited 38 coronal rain reports from the literature. From these earlier works, we picked 15 quiescent coronal rain events out of the solar limb, mostly suggested to occur in active region closed loops due to thermal nonequilibrium, to scrutinize their formation mechanism. Employing the extreme ultraviolet images and line-of-sight magnetograms, the evolution of the quiescent coronal rain events and their magnetic fields and context coronal structures is examined. We find that 6, comprising 40%, of the 15 quiescent coronal rain events could be totally or partially interpreted by the formation mechanism for quiescent coronal rain along open structures facilitated by interchange reconnection. The results suggest that the quiescent coronal rain facilitated by interchange reconnection between open and closed field lines deserves more attention.

The Tracking Ultraviolet Set-up (TUS) is the world's first orbital imaging detector of Ultra High Energy Cosmic Rays (UHECR) and it operated in 2016-2017 as part of the scientific equipment of the Lomonosov satellite. The TUS was developed and manufactured as a prototype of the larger project K-EUSO with the main purpose of testing the efficiency of the method for measuring the ultraviolet signal of extensive air shower (EAS) in the Earth's night atmosphere. Despite the low spatial resolution ($\sim5\times5$ km$^2$ at sea level), several events were recorded which are very similar to EAS as for the signal profile and kinematics. Reconstruction of the parameters of such events is complicated by a short track length, an asymmetry of the image, and an uncertainty in the sensitivity distribution of the TUS channels. An advanced method was developed for the determination of event kinematic parameters including its arrival direction. In the present article, this method is applied for the analysis of 6 EAS-like events recorded by the TUS detector. All events have an out of space arrival direction with zenith angles less than 40{\deg}. Remarkably they were found to be over the land rather close to United States airports, which indicates a possible anthropogenic nature of the phenomenon. Detailed analysis revealed a correlation of the reconstructed tracks with direction to airport runways and Very High Frequency (VHF) omnidirectional range stations. The method developed here for reliable reconstruction of kinematic parameters of the track-like events, registered in low spatial resolution, will be useful in future space missions, such as K-EUSO.

The High Altitude Water Cherenkov (HAWC) observatory and the High Energy Stereoscopic System (H.E.S.S.) are two leading instruments in the ground-based very-high-energy gamma-ray domain. HAWC employs the water Cherenkov detection (WCD) technique, while H.E.S.S. is an array of Imaging Atmospheric Cherenkov Telescopes (IACTs). The two facilities therefore differ in multiple aspects, including their observation strategy, the size of their field of view and their angular resolution, leading to different analysis approaches. Until now, it has been unclear if the results of observations by both types of instruments are consistent: several of the recently discovered HAWC sources have been followed up by IACTs, resulting in a confirmed detection only in a minority of cases. With this paper, we go further and try to resolve the tensions between previous results by performing a new analysis of the H.E.S.S. Galactic plane survey data, applying an analysis technique comparable between H.E.S.S. and HAWC. Events above 1 TeV are selected for both datasets, the point spread function of H.E.S.S. is broadened to approach that of HAWC, and a similar background estimation method is used. This is the first detailed comparison of the Galactic plane observed by both instruments. H.E.S.S. can confirm the gamma-ray emission of four HAWC sources among seven previously undetected by IACTs, while the three others have measured fluxes below the sensitivity of the H.E.S.S. dataset. Remaining differences in the overall gamma-ray flux can be explained by the systematic uncertainties. Therefore, we confirm a consistent view of the gamma-ray sky between WCD and IACT techniques.

This paper examines long-term temporal and spatial fluctuations in the solar rotation (more than four solar cycles) by investigating radio emission escapes from various layers of the solar atmosphere during the years 1967-2010. The flux modulation approach can also be used to investigate variations in solar rotation, which is a contentious topic in solar physics. The current study makes use of a time series of radio flux data at different frequencies (245-15400 MHz) obtained at Sagamore Hill Solar Radio Observatory in Massachusetts, USA, and other observatories from 1967 to 2010. The periodicity present in the temporal variation of time series is estimated through Lomb Scargle Periodogram (LSP). The rotation period estimated for five radio emissions (606, 1415, \& 2695 MHz; from corona, and 4995 \& 8800 MHz; from transition region) through statistical approach shows continuous temporal and spatial variation throughout the years. The smoothed rotation period shows the presence of $\sim$ 22-yrs periodic and $\sim$ 11-yrs components in it. The 22-year component could be linked to the reversal of the solar magnetic field (Hale's) cycle, while the 11-yrs component is most likely related to the sunspot (Schwabe's) cycle. Besides these two components, random components are also prominently present in the analyzed data. The cross-correlation between the sunspot number and the rotation period obtained shows a strong correlation with 11-yrs Schwabe's and 22-yr Hale cycle. The corona rotates faster or slower than transition region in different epoch. The swap of faster rotation speed between corona and transition region also follows the 22-yrs cycle.

It is not yet fully understood how planet formation affects the properties of host stars, in or out of a cluster; however, abundance trends can help us understand these processes. We present a detailed chemical abundance analysis of six stars in Praesepe, a planet-hosting open cluster. Pr0201 is known to host a close-in (period of 4.4 days) giant planet (mass of 0.54$\rm M_{\rm{J}}$), while the other five cluster members in our sample (Pr0133, Pr0081, Pr0208, Pr0051, and Pr0076) have no detected planets according to RV measurements. Using high-resolution, high signal-to-noise echelle spectra obtained with Keck/HIRES and a novel approach to equivalent width measurements (XSpect-EW), we derived abundances of up to 20 elements spanning a range of condensation temperatures (Tc). We find a mean cluster metallicity of [Fe/H] = +0.21$\pm$0.02 dex, in agreement with most previous determinations. We find most of our elements show a [X/Fe] scatter of $\sim$0.02-0.03 dex and conclude that our stellar sample is chemically homogeneous. The Tc slope for the cluster mean abundances is consistent with zero and none of the stars in our sample exhibit individually a statistically significant Tc slope. Using a planet engulfment model, we find that the planet-host, Pr0201, shows no evidence of significant enrichment in its refractory elements when compared to the cluster mean that would be consistent with a planetary accretion scenario.

Four stars pulsating simultaneously with a dominant period $P_D$$\in$(0.28,0.39) d and an {\it additional} period $P_A$$\in$(0.20,0.27) d have been identified from among the more than 3000 RR Lyrae stars observed by the Kepler space telescope during NASA's K2 Mission. All four stars are located in the direction of the Galactic Bulge and have period ratios, $P_A$/$P_D$, significantly smaller than those of most double-mode RR Lyrae (RRd) stars: $P_A$/$P_D$$\in$(0.694,0.710) vs. $P_1$/$P_0$$\in$(0.726,0.748). Three of the stars are faint ($<$$V$$>$=18--20 mag) and distant and are among the `peculiar' RRd (pRRd) stars discovered by Prudil et al. (2017); the fourth star, EPIC 216764000 (=V1125 Sgr), is a newly discovered pRRd star several magnitudes brighter than the other three stars. In this paper the high-precision long-cadence K2 photometry is analyzed in detail and used to study the cycle-to-cycle light variations. The pulsational characteristics of pRRd stars are compared with those of `classical' and `anomalous' RRd (cRRd, aRRd) stars. The conclusion by Prudil et al. that pRRd stars form a separate group of double-mode pulsators and are not simply very-short-period cRRd stars is confirmed. V1127 Aql and AH Cam are identified as other probable members of the class of pRRd stars.

Ultra-precise astrometry from the Gaia mission is expected to lead to astrometric detections of more than 20,000 exoplanets in our Galaxy. One of the factors that could hamper such detections is the astrometric jitter caused by the magnetic activity of the planet host stars. In our previous study, we modeled astrometric jitter for the Sun observed equator-on. In this work, we generalize our model and calculate the photocenter jitter as it would be measured by the Gaia and Small-JASMINE missions for stars with solar rotation rate and effective temperature, but with various values of the inclination angle of the stellar rotation axis. In addition, we consider the effect of metallicity and of nesting of active regions (i.e. the tendency of active regions to emerge in the vicinity of each other). We find that, while the jitter of stars observed equator-on does not have any long-term trends and can be easily filtered out, the photocenters of stars observed out of their equatorial planes experience systematic shifts over the course of the activity cycle. Such trends allow the jitter to be detected with continuous measurements, in which case it can interfere with planet detectability. An increase in the metallicity is found to increase the jitter caused by stellar activity. Active-region nesting can further enhance the peak-to-peak amplitude of the photocenter jitter to a level that could be detected by Gaia.

With the advent of next generation high resolution telescopes, our understanding of how the magnetic field is organized in the internetwork photosphere is likely to advance.We aim to evaluate the extent to which we can retrieve information about the magnetic vector in the internetwork (IN) photosphere using inversions. We use snapshots produced from high resolution 3D magnetohydrodynamic (MHD) simulations and employ the Stokes Inversions based on Response functions (SIR) code to produce synthetic observables in the near infrared spectral window observed by the GREGOR Infrared Spectrograph (GRIS), which contains the highly magnetically sensitive photospheric Fe I line pair at 15648.52 A and 15652.87 A. We perform nearly 14 million inversions to test how well the true MHD atmospheric parameters can be constrained. Finally, we degrade the synthetic Stokes vectors spectrally and spatially to GREGOR resolutions and examine how this influences observations, considering the impact of stray light, spatial resolution and signal-to-noise (S-to-N). We find the depth-averaged parameters can be recovered by the inversions of the undegraded profiles, and by adding gradients to magnetic field strength, inclination and line of sight velocity we show an improvement in the chi squared value is achieved. We evaluate the extent to which we can constrain these parameters at various optical depths, with the kinematic and thermodynamic parameters sensitive deeper in the atmosphere than the magnetic parameters. We find the S-to-N and spatial resolution play a significant role in determining how the atmosphere appears and the magnetic and kinematic parameters are invariant upon inclusion of unpolarized stray light. We studied a linear polarization feature which resembles those recently observed by GRIS, appearing as loop-like structures with similar magnetic flux density.

The synchrotron mechanism has the radiation limit of about 160 MeV, and it is not possible to explain the very high energy (VHE) photons that are emitted by high-energy objects. Inverse Compton scattering as a traditional process is applied for the explanation of the VHE emission. In this paper, jitter radiation, the relativistic electron radiation in the random and small-scale magnetic field, is proposed to be a possible mechanism to produce VHE photons. The jitter radiation frequency is associated with the perturbation field. The spectral index of the jitter radiation is dominated by the kinetic turbulence. We utilize the jitter radiation to explain the gamma-ray burst (GRB 190114C and GRB 180720B) VHE emissions that were recently detected by the Imaging Atmospheric Cherenkov Telescopes. We suggest that this mechanism can be applied to other kinds of VHE sources.

In Sun and sun-like stars, it is believed that the cycles of the large-scale magnetic field are produced due to the existence of differential rotation and helicity in the plasma flows in their convection zones (CZs). Hence, it is expected that for each star, there is a critical dynamo number for the operation of a large-scale dynamo. As a star slows down, it is expected that the large-scale dynamo ceases to operate above a critical rotation period. In our study, we explore the possibility of the operation of the dynamo in the subcritical region using the Babcock--Leighton type kinematic dynamo model. In some parameter regimes, we find that the dynamo shows hysteresis behavior, i.e., two dynamo solutions are possible depending on the initial parameters -- decaying solution if started with weak field and strong oscillatory solution (subcritical dynamo) when started with a strong field. However, under large fluctuations in the dynamo parameter, the subcritical dynamo mode is unstable in some parameter regimes. Therefore, our study supports the possible existence of subcritical dynamo in some stars which was previously shown in a mean-field dynamo model with distributed $\alpha$ and MHD turbulent dynamo simulations.

Coronal upflows at the edges of active regions (AR), which are a possible source of slow solar wind, have been found to connect with dynamics in the transition region. To infer at what scale transition region dynamics connect to AR upflows, we investigate the statistical properties of the small-scale dynamics in the transition region underneath the upflows at the edge of AR NOAA 11934. With observations from the Interface Region Imaging Spectragraph (IRIS), we found that the Si IV 1403\,\AA\ Doppler map consists of numerous blue-shifted and red-shifted patches mostly with sizes less than 1\,$Mm^2$. The blue-shifted structures in the transition region tend to be brighter than the red-shifted ones, but their nonthermal velocities have no significant difference. With the SWAMIS feature tracking procedure, in IRIS slit-jaw 1400\,\AA\ images we found that dynamic bright dots with an average size of about 0.3\,$Mm^2$ and lifetimes mostly less than 200\,s spread all over the region. Most of the bright dots appear to be localised, without clear signature of propagation of plasma to a long distance on the projection plane. Surge-like motions with speeds about 15 km/s could be seen in some events at the boundaries of the upflow region, where the magnetic field appear to be inclined. We conclude that the transition region dynamics connecting to coronal upflows should occur in very fine scale, suggesting that the corresponding coronal upflows should also be highly-structured. It is also plausible that the transition region dynamics might just act as stimulation at the coronal base that then drives the upflows in the corona.

Earth's magnetosphere is vital for today's technologically dependent society. To date, numerous design studies have been conducted and over a dozen science missions have own to study the magnetosphere. However, a majority of these solutions relied on large monolithic satellites, which limited the spatial resolution of these investigations, as did the technological limitations of the past. To counter these limitations, we propose the use of a satellite swarm carrying numerous and distributed payloads for magnetospheric measurements. Our mission is named APIS (Applications and Potentials of Intelligent Swarms), which aims to characterize fundamental plasma processes in the Earth's magnetosphere and measure the effect of the solar wind on our magnetosphere. We propose a swarm of 40 CubeSats in two highly-elliptical orbits around the Earth, which perform radio tomography in the magnetotail at 8-12 Earth Radii (RE) downstream, and the subsolar magnetosphere at 8-12RE upstream. In addition, in-situ measurements of the magnetic and electric fields, plasma density will be performed by on-board instruments. In this article, we present an outline of previous missions and designs for magnetospheric studies, along with the science drivers and motivation for the APIS mission. Furthermore, preliminary design results are included to show the feasibility of such a mission. The science requirements drive the APIS mission design, the mission operation and the system requirements. In addition to the various science payloads, critical subsystems of the satellites are investigated e.g., navigation, communication, processing and power systems. We summarize our findings, along with the potential next steps to strengthen our design study.

The landscape of high- and ultra-high-energy astrophysics has changed in the last decade, in large part owing to the inflow of high-quality data collected by present cosmic-ray, gamma-ray, and neutrino observatories. At the dawn of the multimessenger era, the interpretation of these observations within a consistent framework is important to elucidate the open questions in this field. CRPropa 3.2 is a Monte Carlo code for simulating the propagation of high-energy particles in the Universe. This new version represents a step further towards a more complete simulation framework for multimessenger studies. Some of the new developments include: cosmic-ray acceleration, support for particle interactions within astrophysical sources, full Monte Carlo treatment of electromagnetic cascades, improved ensemble-averaged Galactic propagation, and a number of technical enhancements. Here we present some of these novel features and some applications to gamma- and cosmic-ray propagation.

Observations of the redshifted 21-cm line of neutral hydrogen (HI) are a new and powerful window of observation that offers us the possibility to map the spatial distribution of cosmic HI and learn about cosmology. BINGO (Baryon Acoustic Oscillations [BAO] from Integrated Neutral Gas Observations) is a new unique radio telescope designed to be one of the first to probe BAO at radio frequencies. BINGO has two science goals: cosmology and astrophysics. Cosmology is the main science goal and the driver for BINGO's design and strategy. The key of BINGO is to detect the low redshift BAO to put strong constraints in the dark sector models. Given the versatility of the BINGO telescope, a secondary goal is astrophysics, where BINGO can help discover and study Fast Radio Bursts (FRB) and other transients, Galactic and extragalactic science. In this paper, we introduce the latest progress of the BINGO project, its science goals, describing the scientific potential of the project in each science and the new developments obtained by the collaboration. We introduce the BINGO project and its science goals and give a general summary of recent developments in construction, science potential and pipeline development obtained by the BINGO collaboration in the past few years. We show that BINGO will be able to obtain competitive constraints for the dark sector, and also that will allow for the discovery of several FRBs in the southern hemisphere. The capacity of BINGO in obtaining information from 21-cm is also tested in the pipeline introduced here. There is still no measurement of the BAO in radio, and studying cosmology in this new window of observations is one of the most promising advances in the field. The BINGO project is a radio telescope that has the goal to be one of the first to perform this measurement and it is currently being built in the northeast of Brazil. (Abridged)

The measurement of the diffuse $21$-cm radiation from the hyperfine transition of neutral hydrogen (HI signal) in different redshifts is an important tool for modern cosmology. However, detecting this faint signal with non-cryogenic receivers in single-dish telescopes is a challenging task. The BINGO (Baryon Acoustic Oscillations from Integrated Neutral Gas Observations) radio telescope is an instrument designed to detect baryonic acoustic oscillations (BAO) in the cosmological HI signal, in the redshift interval $0.127 \le z \le 0.449$. This paper describes the BINGO radio telescope, including the current status of the optics, receiver, observational strategy, calibration and the site. BINGO has been carefully designed to minimize systematics, being a transit instrument with no moving dishes and 28 horns operating in the frequency range $980 \le \nu \le 1260$ MHz. Comprehensive laboratory tests were conducted for many of the BINGO subsystems and the prototypes of the receiver chain, horn, polarizer, magic tees and transitions have been successfully tested between 2018-2020. The survey was designed to cover $\sim 13\%$ of the sky, with the primary mirror pointing at declination $\delta=-15^{\circ}$. The telescope will see an instantaneous declination strip of $14.75^{\circ}$. The results of the prototype tests closely meet those obtained during the modelling process, suggesting BINGO will perform according to our expectations. After one year of observations with a 60% duty cycle, BINGO should achieve an expected sensitivity of $102 \mu K$ for 28 horns and 30 redshift bins, considering one polarization and be able to measure the HI power spectrum in a competitive time frame.

The BINGO telescope was designed to measure the fluctuations of the 21-cm radiation arising from the hyperfine transition of neutral hydrogen and aims to measure the Baryon Acoustic Oscillations (BAO) from such fluctuations, therefore serving as a pathfinder to future deeper intensity mapping surveys. The requirements for the Phase 1 of the projects consider a large reflector system (two 40 m-class dishes in a crossed-Dragone configuration), illuminating a focal plane with 28 horns to measure the sky with two circular polarisations in a drift scan mode to produce measurements of the radiation in intensity as well as the circular polarisation. In this paper we present the optical design for the instrument. We describe the intensity and polarisation properties of the beams and the optical arrangement of the horns in the focal plane to produce a homogeneous and well-sampled map after the end of Phase 1. Our analysis provides an optimal model for the location of the horns in the focal plane, producing a homogeneous and Nyquist sampled map after the nominal survey time. We arrive at an optimal configuration for the optical system, including the focal plane positioning and the beam behavior of the instrument. We present an estimate of the expected side lobes both for intensity and polarisation, as well as the effect of band averaging on the final side lobes. The cross polarisation leakage values for the final configuration allow us to conclude that the optical arrangement meets the requirements of the project. We conclude that the chosen optical design meets the requirements for the project in terms of polarisation purity, area coverage as well as homogeneity of coverage so that BINGO can perform a successful BAO experiment. We further conclude that the requirements on the placement and r.m.s. error on the mirrors are also achievable so that a successful experiment can be conducted.(Abridged)

The large-scale distribution of neutral hydrogen (HI) in the Universe is luminous through its 21-cm emission. The goal of the Baryon Acoustic Oscillations from Integrated Neutral Gas Observations - BINGO radio telescope is to detect Baryon Acoustic Oscillations (BAO) at radio frequencies through 21-cm Intensity Mapping (IM). The telescope will span the redshift range 0.127 $< z <$ 0.449 with an instantaneous field-of-view of $14.75^{\circ} \times 6.0^{\circ}$. In this work, we investigate different constructive and operational scenarios of the instrument by generating sky maps as they should be produced by the instrument. In doing this we use a set of end-to-end IM mission simulations. The maps will additionally also be used to evaluate the efficiency of a component separation method (GNILC). We have simulated the kind of data that would be produced in a single-dish IM experiment like BINGO. According to the results obtained we have optimized the focal plane design of the telescope. In addition, the application of the GNILC method on simulated data shows that it is feasible to extract the cosmological signal across a wide range of multipoles and redshifts. The results are comparable with the standard PCA method.

Observing the neutral Hydrogen (HI) distribution across the Universe via redshifted 21-cm line Intensity Mapping (IM) constitutes a powerful probe for cosmology. However, this redshifted 21cm signal is obscured by the foreground emission. This paper addresses the capabilities of the BINGO survey to separate such signals. Specifically, this paper (a) looks in detail at the different residuals left over by foreground components, (b) shows that a noise-corrected spectrum is unbiased and (c) shows that we understand the remaining systematic residuals by analyzing non-zero contributions to the three point function. We use the Generalized Needlet Internal Linear Combination (GNILC), which we apply to sky simulations of the BINGO experiment for each redshift bin of the survey. We present our recovery of the redshifted 21-cm signal from sky simulations of the BINGO experiment including foreground components. We test the recovery of the 21-cm signal through the angular power spectrum at different redshifts, as well as the recovery of its non-Gaussian distribution through a bispectrum analysis. We find that non-Gaussianities from the original foreground maps can be removed down to, at least, the noise limit of the BINGO survey with such techniques. Our bispectrum analysis yields strong tests of the level of the residual foreground contamination in the recovered 21-cm signal, thereby allowing us to both optimize and validate our component separation analysis. (Abridged)

BINGO (Baryon Acoustic Oscillations from Integrated Neutral Gas Observations.) is a radio telescope designed to survey from 980 MHz to 1260 MHz, observe the neutral Hydrogen (HI) 21-cm line and detect BAO (Baryon Acoustic Oscillation) signal with Intensity Mapping technique. Here we present our method to generate mock maps of the 21-cm Intensity Mapping signal covering the BINGO frequency range and related test results. (Abridged)

The 21-cm line of neutral hydrogen (HI) opens a new avenue in our exploration of the Universe's structure and evolution. It provides complementary data with different systematics, which aim to improve our current understanding of the $\Lambda$CDM model. Among several radio cosmological surveys designed to measure this line, BINGO is a single dish telescope mainly designed to detect Baryon Acoustic Oscillations (BAO) at low redshifts ($0.127 < z < 0.449$). Our goal is to assess the capabilities of the fiducial BINGO setup to constrain the cosmological parameters and analyse the effect of different instrument configurations. We will use the 21-cm angular power spectra to extract information about the HI signal and the Fisher matrix formalism to study BINGO projected constraining power. We use the Phase 1 fiducial configuration of the BINGO telescope to perform our cosmological forecasts. In addition, we investigate the impact of several instrumental setups and different cosmological models. Combining BINGO with Planck temperature and polarization data, we project a $1\%$ and a $3\%$ precision measurement at $68\%$ CL for the Hubble constant and the dark energy (DE) equation of state (EoS), respectively, within the wCDM model. Assuming a CPL parametrization, the EoS parameters have standard deviations given by $\sigma_{w_0} = 0.30$ and $\sigma_{w_a} = 1.2$. We find that BINGO can also help breaking degeneracies in alternative models, which improves the cosmological constraints significantly. Moreover, we can access information about the HI density and bias, obtaining $\sim 8.5\%$ and $\sim 6\%$ precision, respectively, assuming they vary with redshift at three independent bins. The fiducial BINGO configuration will be able to extract significant information from the HI distribution and provide constraints competitive with current and future cosmological surveys. (Abridged)

Previous X-ray studies of the Perseus Cluster, consisting of 85 Suzaku pointings along eight azimuthal directions, revealed a particularly steep decrease in the projected temperature profile near the virial radius (~r200) towards the northwest (NW). To further explore this shock candidate, another 4 Suzaku observations on the NW edge of the Perseus Cluster have been obtained. These deeper data were designed to provide the best possible control of systematic uncertainties in the spectral analysis. Using the combined Suzaku observations, we have carefully investigated this interesting region by analyzing the spectra of various annuli and extracting projected thermodynamic profiles. We find that the projected temperature profile shows a break near r200, indicating a shock with M = 1.9+-0.3. Corresponding discontinuities are also found in the projected emission measure and the density profiles at the same location. This evidence of a shock front so far away from the cluster center is unprecedented, and may provide a first insight into the properties of large-scale virial shocks which shape the process of galaxy cluster growth.

In this work we present a mathematical model for the propagation of the shock waves that occur in graded density profiles. These waves can occur in a wide range of astrophysical events, such as collisions in planetary and stellar atmospheres, common envelope explosions and peculiar type Ia supernovae. The behaviour of the shock wave and its evolution can be modelled using type II self similar solutions. In such solutions the evolution of the shock wave is determined by boundary conditions at the shock front and a singular point in the shocked region. We show how the evolution can be determined for different equations of state and density profiles, and compare these results to numerical simulations. These findings are also applied to a variety of astrophysical phenomena to further test their validity.

Coronal mass ejections (CMEs) represent one type of the major eruption from the Sun. Their interplanetary counterparts, the interplanetary CMEs (ICMEs), are the direct manifestations of these structures when they propagate into the heliosphere and encounter one or more observing spacecraft. The ICMEs generally exhibit a set of distinctive signatures from the in-situ spacecraft measurements. A particular subset of ICMEs, the so-called Magnetic Clouds (MCs), is more uniquely defined and has been studied for decades, based on in-situ magnetic field and plasma measurements. By utilizing the latest multiple spacecraft measurements and analysis tools, we report a detailed study of the internal magnetic field configuration of an MC event observed by both the Solar Orbiter (SO) and Wind spacecraft in the solar wind near the Sun-Earth line. Both two-dimensional (2D) and three-dimensional (3D) models are applied to reveal the flux rope configurations of the MC. Various geometrical as well as physical parameters are derived and found to be similar within error estimates for the two methods. These results quantitatively characterize the coherent MC flux rope structure crossed by the two spacecraft along different paths. The implication for the radial evolution of this MC event is also discussed.

We present novel cosmological constraints based on a joint analysis of our HII galaxies (HIIG) Hubble relation with the full Planck Cosmic Microwave Background anisotropy spectrum and the Baryon Acoustic Oscillations (BAO) probes. The HII galaxies span a large redshift range $(0.088 \le z \le 2.5)$, reaching significantly higher redshifts than available SNIa and hence they probe the cosmic expansion at earlier times. Our independent constraints compare well with those based on the Pantheon compilation of SNIa data, which we also analyse. We find our results to be in agreement with the conformal $\Lambda$CDM model within 1$\sigma$. We also use our HIIG data to examine the behaviour of the dark energy equation of state parameter under the CPL parameterisation, $w = w_0+w_a \frac{z}{1+z}$, and find consistent results with those based on SNIa, although the degeneracy in the parameter space as well as the individual parameter uncertainties, when marginalizing one over the other, are quite large.

We calculate the neutrino signal from Population III supermassive star collapse using a neutrino transfer code originally developed for core collapse supernovae and massive star collapse. Using this code, we are able to investigate the supermassive star mass range thought to undergo neutrino trapping ($\sim 10^4$ M$_\odot$), a mass range which has been neglected by previous works because of the difficulty of neutrino transfer. For models in this mass range, we observe a neutrino-sphere with a large radius and low density compared to typical massive star neutrino-spheres. We calculate the neutrino light-curve emitted from this neutrino-sphere. The resulting neutrino luminosity is significantly lower than the results of a previous analytical model. We briefly discuss the possibility of detecting a neutrino burst from a supermassive star or the neutrino background from many supermassive stars and conclude that the former is unlikely with current technology, unless the SMS collapse is located as close as 1 Mpc, while the latter is also unlikely even under very generous assumptions. However, the supermassive star neutrino background is still of interest as it may serve as a source of noise in proposed dark matter direct detection experiments.

Brightest cluster galaxies (BCGs), particularly those at the centers of cool-core clusters, can exhibit star formation over spatial extents of up to $\gtrsim$100\,kpc at inferred rates of up to $\gtrsim100\rm\,M_\sun\,yr^{-1}$. Is their star formation also extended over time, as might be expected if fuelled by cooling of the surrounding hot intracluster gas -- a residual cooling flow -- as demonstrated hitherto only for the BCG in the Perseus cluster? Here, to infer the formation history of relatively young stars in the BCG of MACS\,J0329.7$-$0211, we fit model single-stellar-populations to the spectral energy distributions (spanning near-UV to near-IR) measured along different sightlines towards its young stellar population. Employing a Markov Chain Monte Carlo method, we show that star formation in this BCG has persisted at a relatively constant rate of $\sim2{\rm\,M_\sun\,yr^{-1}}$ (factors of 10--40 below the rates previously inferred using simpler methods and/or ad hoc assumptions) over the past $\sim$400\,Myr, beyond which any star formation falls below the observational detection threshold. Such persistent star formation from a residual cooling flow can contribute up to $\sim$10\% of the original stellar mass of this BCG if its progenitor was among the most massive red nuggets known at $z\sim$2 having masses of $\sim1\times10^{11}\rm\,M_\sun$, but only a few percent of its overall growth in stellar mass to $\sim8\times10^{11}\rm\,M_\sun$ at $z=0.45$. Although constituting only a minor pathway for the stellar growth of this BCG, persistent star formation from a residual cooling flow can nevertheless contribute significantly to the enormous number of globular clusters found around BCGs in the local Universe.

Radiation pressure on dust is thought to play a crucial role in the formation process of massive stars by acting against gravitational collapse onto the central protostar. However, dust properties in dense regions irradiated by the intense radiation of massive protostars are poorly constrained. Previous studies usually assume the standard interstellar dust model to constrain the maximum mass of massive stars formed by accretion, which appears to contradict with dust evolution theory. In this paper, using the fact that stellar radiation exerts on dust simultaneous radiation pressure and radiative torques, we study the effects of grain rotational disruption by radiative torques (RATs) on radiation pressure and explore its implications for massive star formation. For this paper, we focus on the protostellar envelope and adopt a spherical geometry. We find that original large grains of micron-sizes presumably formed in very dense regions can be rapidly disrupted into small grains by RATs due to infrared radiation from the hot dust shell near the sublimation front induced by direct stellar radiation. Owing to the modification in the size distribution by rotational disruption, the radiation pressure opacity can be decreased by a factor of $\sim 3$ from the value expected from the original dust model. However, to form massive stars via spherical accretion, the dust-to-gas mass ratio needs to be reduced by a factor of $\sim 5$ as previously found.

We present a deep far and near-ultraviolet (FUV and NUV) wide-field imaging survey of galaxies in the Bootes Void using Ultra-Violet Imaging Telescope onboard {\em AstroSat}. Our data reach $5\sigma$ limiting magnitudes for point sources at 23.0 and 24.0 AB mag in FUV and NUV respectively. We report a total of six star-forming galaxies residing in the Bootes Void alongside the full catalog, and of these, three are newly detected in our FUV observation. Our void galaxy sample spans a range of UV colors $(-0.35\, \leq$ FUV$-$NUV $\leq\, 0.68)$ and absolute magnitudes $(-14.16\, \leq\, \mathrm{M_{NUV}}\, \leq\, -18.65)$. In addition, {\em Sloan Digital Sky Survey} and {\em Two-micron All Sky Survey} archival data are being used to study UV, optical, and infrared color-magnitude relations for our galaxies in the void. We investigate the nature of bi-modal color distribution, morphologies, and star formation of the void galaxies. Most of the galaxies in our sample are fainter and less massive than L$^{\ast}$ galaxies, with M$_\mathrm{r} > -20$. Our analysis reveals a dominant fraction of bluer galaxies over the red ones in the void region probed. The internal and Galactic extinction corrected FUV star formation rates (SFRs) in our void galaxy catalog varies in a large range of $0.05$ to $51.01$ M$_{\odot} yr^{-1}$, with a median $3.96$ M$_{\odot} yr^{-1}$. We find a weak effect of the environment on the SFRs of galaxies. Implications of our findings are discussed.

The upgraded Giant Metrewave Radio Telescope (GMRT) has been used to map the cluster A2065 at z = 0.0726. We report the discovery of a remnant radio galaxy at the peripheral cluster region. The spatially resolved radio emission from the remnant radio galaxy shows an elongated, bar-shaped structure, whose size is $\approx$ 52$^{\prime\prime}$ $\times$ 110$^{\prime\prime}$ ($\simeq$ 72 $\times$ 152 kpc$^2$). Our study with the multiwavelength GMRT data and \textit{Chandra} data shows that across the remnant radio galaxy there is a hint of a surface-brightness edge in the hot X-ray gas. We detect tentative flattening of the radio spectral index as the old plasma at the near end of the surface-brightness edge is reinvigorated by the passage of possible shock front and shows the expected change in radio emission characteristics. We suggest that the remnant radio galaxy has been seeded by the lobes of the active galactic nucleus (AGN), hosted by the WISEA J152228.01$+$274141.3 source, demonstrating the connection between AGNs and remnant radio sources. Although the number of known remnant radio sources is beginning to increase, we emphasize the need for better data to understand the physics and nature of poorly understood remnant radio sources.

We present results from a polarization study of the radio-intermediate quasar, III Zw 2, at a redshift of 0.0893, with the upgraded Giant Metrewave Radio Telescope (uGMRT) at 685 MHz and the Karl Jansky Very Large Array (VLA) at 5 and 34 GHz. We detect a kpc-scale outflow, exhibiting transverse magnetic (B-) fields. The curved jet terminates in a bow-shock-like radio structure with inferred B-fields aligned with the lobe edges. We suggest that the radio outflow in III Zw 2 is a combination of a collimated jet along with a wind-like component. This "wind" component could be a magnetized accretion disk wind or the outer layers of a broadened jet or a combination of both. The current data cannot differentiate between these possibilities. We also detect kpc-scale lobe emission that is misaligned with the primary lobes in the uGMRT images. The spectral indices and the electron lifetimes in the misaligned lobe are similar to the primary lobe, suggesting that the misaligned lobe is not a relic. We propose that changing spectral states of the accretion disk, and the subsequent intermittent behaviour of the outflow, along with the close interplay between the jet and "wind" could explain the radio-intermediate nature of III Zw 2. Our study shows that radio-intermediate quasars are promising sources for understanding the role of jets and winds in galaxy evolution and demonstrates the power of radio polarization studies towards achieving this.

Aimed at progress in MeV gamma-ray astronomy which has not yet been well-explored, Compton telescope missions with a variety of detector concepts have been proposed so far. One of the key techniques for these future missions is an event reconstruction algorithm that is able to determine the scattering orders of multiple Compton scattering events and to identify events in which gamma rays escape from the detectors before they deposit all of their energies. We propose a new algorithm that can identify whether the gamma rays escape from the detectors or not, in addition to the scattering order determination. This algorithm also corrects incoming gamma-ray energies for escape events. The developed algorithm is based on the maximum likelihood method, and we present a general formalism of likelihood functions describing physical interactions inside the detector. We also introduce several approximations in the calculation of the likelihood functions for efficient computation. For validation, we have applied the algorithm to simulation data of a Compton telescope using a liquid argon time projection chamber, which is a new type of Compton telescope proposed for the GRAMS mission, and have confirmed that it works successfully for up to 8-hit events. The proposed algorithm can be used for next-generation MeV gamma-ray missions featured by large-volume detectors, e.g., GRAMS and AMEGO.

Accurate photo-z measurements are important to construct a large-scale structure map of X-ray Universe in the ongoing SRG/eROSITA All-Sky Survey. We present machine learning Random Forest-based models for probabilistic photo-z predictions based on information from 4 large photometric surveys (SDSS, Pan-STARRS, DESI Legacy Imaging Survey, and WISE). Our models are trained on the large sample of $\approx$580000 quasars and galaxies selected from the SDSS DR14 spectral catalog and take into account Galactic extinction and uncertainties in photometric measurements for target objects. On the Stripe82X test sample we obtained photo-z accuracy for X-ray sources: $NMAD=0.034$ (normalized median absolute deviation) and $n_{>0.15}=0.088$ (catastrophic outliers fraction), which is almost $\sim2$ times better than best photo-z results available in the literature.

During solar minimum, the Sun is relatively inactive with few sunspots observed on the solar surface. Consequently, we observe a smaller number of highly energetic events such as solar flares or coronal mass ejections (CMEs), which are often associated with active regions on the photosphere. Nonetheless, our magnetofrictional simulations during the minimum period suggest that the solar corona is still dynamically evolving in response to the large-scale shearing velocities on the solar surface. The non-potential evolution of the corona leads to the accumulation of magnetic free energy and helicity, which is periodically shed in eruptive events. We find that these events fall into two distinct classes: One set of events are caused by eruption and ejection of low-lying coronal flux ropes, and they could explain the origin of occasional CMEs during solar minimum. The other set of events are not driven by destabilisation of low-lying structures but rather by eruption of overlying sheared arcades. These could be associated with streamer blowouts or stealth CMEs. The two classes differ significantly in the amount of magnetic flux and helicity shed through the outer coronal boundary. We additionally explore how other measurables such as current, open magnetic flux, free energy, coronal holes, and the horizontal component of the magnetic field on the outer model boundary vary during the two classes of event. This study emphasises the importance and necessity of understanding the dynamics of the coronal magnetic field during solar minimum.

Studying Class 0 objects is very important, as it allows to characterize dynamical processes at the onset of the star formation process, and to determine the physical mechanisms responsible for the outcome of the collapse. Observations of dense gas tracers allow the characterization of key kinematics of the gas directly involved in the star-formation process, such as infall, outflow or rotation. This work aims at investigating the molecular line velocity profiles of the Class 0 protostellar object B335 and attempts to put constraints on the infall motions happening in the circumstellar gas of the object.} Observations of C$^{17}$O (1-0), C$^{18}$O (1-0) and $^{12}CO$ (2-1) transitions are presented and the spectral profiles are analyzed at envelope radii between 100 and 860 au. C$^{17}$O emission presents a double peaked line profile distributed in a complex velocity field. Both peaks present an offset of 0.2 to 1 km s$^{-1}$ from the systemic velocity of the source in the probed area. The optical depth of the C$^{17}$O emission has been estimated and found to be less than 1, suggesting that the two velocity peaks trace two distinct velocity components of the gas in the inner envelope. After discarding possible motions that could produce the complex velocity pattern, such as rotation and outflow, it is concluded that infall is producing the velocity field. Because inside-out symmetric collapse cannot explain those observed profiles, it is suggested that those are produced by non-isotropic accretion from the envelope into the central source along the outflow cavity walls.

We report the discovery of a new unidentified extended $\gamma$-ray source in the Galactic plane named LHAASO J0341+5258 with a pre-trial significance of 8.2 standard deviations above 25 TeV. The best fit position is R.A.$=55.34^{\circ}\pm0.11^{\circ}$ and Dec$=52.97^{\circ}\pm0.07^{\circ}$. The angular size of LHAASO J0341+5258 is $0.29^\circ \pm 0.06^\circ_{stat} \pm0.02^\circ_{sys}$. The flux above 25 TeV is about $20\%$ of the flux of Crab Nebula. Although a power-law fit of the spectrum from 10 TeV to 200 TeV with the photon index $\alpha=2.98 \pm 0.19_{stat} \pm 0.02_{sys}$ is not excluded, the LHAASO data together with the flux upper limit at 10 GeV set by the Fermi LAT observation, indicate a noticeable steepening of an initially hard power-law spectrum %($\alpha \leq 1.75$) spectrum with a cutoff at $\approx 50$~TeV. We briefly discuss the origin of UHE gamma-rays. The lack of an energetic pulsar and a young SNR inside or in the vicinity of LHAASO J0341+5258 challenge, but do not exclude both the leptonic and hadronic scenarios of gamma-ray production.

The US National Park Service (NPS) Night Skies Program measured changes in sky brightness resulting from a countywide lighting retrofit project. The retrofit took place in Chelan County, a gateway community to North Cascades National Park and Lake Chelan National Recreation Area in Washington State. The county retrofitted all 3,693 county-owned high pressure sodium (HPS) street lamps to full cutoff LEDs. This number is about 60% of the County's total outdoor street and area lights. About 80% of the newly installed lights were 3000K in color temperature and 20% were 4000K. The 4000K LEDs were used to meet Washington State Department of Transportation guidelines. To measure sky brightness, we used the NPS night sky camera system before the retrofit started in 2018 and after its completion in 2019. These images were photometrically calibrated and mosaicked together to provide hemispherical images in V band. For comparison with our ground-based measurement, we obtained the satellite imagery taken by Visible Infrared Imaging Radiometer Suite (VIIRS) onboard the Suomi National Polar-orbiting Partnership satellite. Our measurements show that the post-retrofit skyglow became brighter and extended higher in the sky, but upward radiance, as measured by the day-night band radiometer, decreased. These divergent results are likely explained by a substantial increase in light emitted at wavelengths shorter than 500 nm, and a relative decrease in zenith light emission due to better shielded luminaires. These results also demonstrate that earlier models relating VIIRS day-night band data to skyglow will, at a minimum, require substantial revision to account for the different characteristics of solid state luminaires.

Core-collapse supernovae are fascinating astrophysical objects for multimessenger studies. Gravitational waves (GWs) are expected to play a role in the supernova explosion mechanism, but their modelling is also challenging due to the stochastic nature of the dynamics and the vast possible progenitors, and moreover, the GW detection from these objects is still elusive with the already advanced detectors. Low-energy neutrinos will be emitted enormously during the core-collapse explosion and can help for the gravitational wave counterpart search. In this work we develop a multi-messenger strategy to search for such astrophysical objects by exploiting a global network of both low-energy neutrino and gravitational wave detectors. First, we discuss how to improve the detection potential of the neutrino sub-network by exploiting the temporal behaviour of a neutrino burst from a core-collapse supernova. We show that with the proposed approach neutrino detectors can gain at least $10\%$ of detection efficiency at the distance where their efficiency drops. Then, we combine the information provided by GW and neutrino in a multimessenger strategy. In particular, we obtain an increase of the probability to detect the GW signal from a CCSN at $60$ kpc from zero when using only GW analysis to $\sim 33\%$ with our combined GW-$\nu$ approach. Keywords: multimessenger, supernova, core-collapse, low-energy neutrino, gravitational wave.

The $\gamma$-ray emission from stars is induced by the interaction of cosmic rays with stellar atmospheres and photon fields. This emission is expected to come in two components: a stellar disk emission, where $\gamma$-rays are mainly produced in atmospheric showers generated by hadronic cosmic rays, and an extended halo emission, where the high density of soft photons in the surroundings of stars create a suitable environment for $\gamma$-ray production via inverse Compton (IC) scattering by cosmic-ray electrons. Besides the Sun, no other disk or halo from single stars has ever been detected in $\gamma$-rays. However, by assuming a cosmic-ray spectrum similar to that observed on Earth, the predicted $\gamma$-ray emission of super-luminous stars, like e.g. Betelgeuse and Rigel, could be high enough to be detected by the Fermi Large Area Telescope (LAT) after its first decade of operations. In this work, we use 12 years of Fermi-LAT observations along with IC models to study 9 super-luminous nearby stars, both individually and via stacking analysis. Our results show no significant $\gamma$-ray emission, but allow us to restrict the stellar $\gamma$-ray fluxes to be on average $<3.3 \times 10^{-11}$ ph cm$^{-2}$ s$^{-1}$ at a 3$\sigma$ confidence level, which translates to an average local density of electrons in the surroundings of our targets to be less than twice of that observed for the Solar System.

In the standard picture of cosmic ray transport the propagation of charged cosmic rays through turbulent magnetic fields is described as a random walk with cosmic rays scattering on magnetic field turbulence. This is in good agreement with the highly isotropic cosmic ray arrival directions as this diffusion process effectively isotropizes the cosmic ray distribution. High-statistics observatories like IceCube and HAWC have however observed significant deviations from isotropy down to very small angular scales. This is in strong tension with this standard picture of cosmic ray propagation. While large scale multipoles arise naturally, for example due to the earth's motion relative to the isotropic cosmic ray distribution, there is no intuitive mechanism to account for the observed anisotropies at smaller angular scales. By relaxing one of the standard assumptions of quasi linear theory and treating correlations between fluxes of cosmic rays from different directions explicitly we show that higher multipoles also are to be expected from particle propagation through turbulent magnetic fields. We present a first analytical calculation of the angular power spectrum assuming a physically motivated model of the magnetic field turbulence and find good agreement with numerical simulations.

Current and future generations of intensity mapping surveys promise dramatic improvements in our understanding of galaxy evolution and large-scale structure. An intensity map provides a census of the cumulative emission from all galaxies in a given region and redshift, including faint objects that are undetectable individually. Furthermore, cross-correlations between line intensity maps and galaxy redshift surveys are sensitive to the line intensity and clustering bias without the limitation of cosmic variance. Using the Fisher information matrix, we derive simple expressions describing sensitivities to the intensity and bias obtainable for cross-correlation surveys, focusing on cosmic variance evasion. Based on these expressions, we conclude that the optimal sensitivity is obtained by matching the survey depth, defined by the ratio of the clustering power spectrum to noise in a given mode, between the two surveys. We find that mid- to far-infrared space telescopes could benefit from this technique by cross-correlating with coming galaxy redshift surveys such as those planned for the Nancy Grace Roman Space Telescope, allowing for sensitivities beyond the cosmic variance limit.

The IceCube Neutrino Observatory observes neutrinos interacting deep within the South Pole ice. It consists of 5,160 digital optical modules embedded within a cubic kilometer of ice, over depths of 1,450\,m to 2,450\,m. At the lower center of the array is the DeepCore subdetector. Its denser sensor configuration lowers the observable energy threshold to the GeV-scale, facilitating the study of atmospheric neutrino oscillations. The precise reconstruction of neutrino direction is critical in the measurements of oscillation parameters. This work presents a method to reconstruct the zenith angle of GeV-scale events in IceCube by using a convolutional neural network and compares the result to that of the current likelihood-based reconstruction algorithm.

While most galaxies appear to host a central supermassive black hole (SMBH), they are expected to also contain a substantial population of off-center "wandering" SMBHs naturally produced by the hierarchical merger-driven process of galaxy assembly. This population has been recently characterized in an analysis of the Romulus cosmological simulations, which correct for the dynamical forces on SMBHs without artificially pinning them to halo centers. Here we predict an array of electromagnetic signatures for these wanderers. The predicted wandering population of SMBHs from Romulus broadly reproduces the observed spatial offsets of a recent sample of hyperluminous X-ray sources. We predict that the sources with the most extreme offsets are likely to arise from SMBHs within satellite galaxies. These simulations also predict a significant population of secondary active galactic nuclei (AGN) with luminosities at least 10\% that of the central AGN. The majority of galaxies at $z=4$ that host a central AGN with bolometric luminosity $L_\mathrm{bol}>10^{42} \ \mathrm{erg} \; \mathrm{s}^{-1}$ are predicted to host a companion off-center AGN of comparable brightness. We demonstrate that stacked X-ray observations of similar mass galaxies may reveal a halo of collective emission attributable to these wanderers. Finally, because wanderers dominate the population of SMBHs with masses of $\lesssim 10^7\,M_{\odot}$ in Romulus, they may dominate tidal disruption event (TDE) rates at these masses if they retain a stellar component (e.g. a nuclear star cluster). This could warrant an order of magnitude correction to current theoretically estimated TDE rates at low SMBH masses.

LHAASO J0621$+$3755 is a TeV gamma-ray halo newly identified by LHAASO-KM2A. It is likely to be generated by electrons trapped in a slow-diffusion zone around PSR J0622$+$3749 through inverse Compton scattering. When the gamma-ray spectrum of LHAASO-KM2A is fitted, the GeV fluxes derived by the commonly used one-zone normal diffusion model for electron propagation are significantly higher than the upper limits (ULs) of Fermi-LAT. In this work, we respectively adopt the one-zone superdiffusion and two-zone normal diffusion models to solve this conflict. For the superdiffusion scenario, we find that a model with superdiffusion index $\alpha\lesssim1.2$ can meet the constraints of Fermi-LAT observation. For the two-zone diffusion scenario, the size of the slow-diffusion zone is required to be smaller than $\sim50$ pc, which is consistent with theoretical expectations. Future precise measurements of the Geminga halo may further distinguish between these two scenarios for the electron propagation in pulsar halos.

Low-mass compact galaxies (ultracompact dwarfs [UCDs] and compact ellipticals [cEs]) populate the stellar size-mass plane between globular clusters and early-type galaxies. Known to be formed either in-situ with an intrinsically low mass or resulting from the stripping of a more massive galaxy, the presence of a supermassive or an intermediate-mass black hole (BH) could help discriminate between these possible scenarios. With this aim, we have performed a multiwavelength search of active BH activity, i.e. active galactic nuclei (AGN), in a sample of 937 low-mass compact galaxies (580 UCDs and 357 cEs). This constitutes the largest study of AGN activity in these types of galaxies. Based on their X-ray luminosity, radio luminosity and morphology, and/or optical emission line diagnostic diagrams, we find a total of 11 cEs that host an AGN. We also study for the first time the location of both low-mass compact galaxies (UCDs and cEs) and dwarf galaxies hosting AGN on the BH-galaxy scaling relations, finding that low-mass compact galaxies tend to be overmassive in the BH mass-stellar mass plane but not as much in the BH mass-stellar velocity dispersion correlation. This, together with available BH mass measurements for some of the low-mass compact galaxies, supports a stripping origin for the majority of these objects that would contribute to the scatter seen at the low-mass end of the BH-galaxy scaling relations. However, the differences are too large to be explained solely by this scatter, and thus our results suggest that a flattening at such low-masses is also plausible, happening at a velocity dispersion of ~20-40 km/s.

Chromospheric umbral oscillations produce periodic brightenings in the core of some spectral lines, known as umbral flashes. They are also accompanied by fluctuations in velocity, temperature, and, according to several recent works, magnetic field. In this study, we aim to ascertain the accuracy of the magnetic field determined from inversions of the Ca II 8542 \AA\ line. We have developed numerical simulations of wave propagation in a sunspot umbra. Synthetic Stokes profiles emerging from the simulated atmosphere were computed and then inverted using the NICOLE code. The atmospheres inferred from the inversions have been compared with the original parameters from the simulations. Our results show that the inferred chromospheric fluctuations in velocity and temperature match the known oscillations from the numerical simulation. In contrast, the vertical magnetic field obtained from the inversions exhibits an oscillatory pattern with a $\sim$300 G peak-to-peak amplitude which is absent in the simulation. We have assessed the error in the inferred parameters by performing numerous inversions with slightly different configurations of the same Stokes profiles. We find that when the atmosphere is approximately at rest, the inversion tends to favor solutions that underestimate the vertical magnetic field strength. On the contrary, during umbral flashes, the values inferred from most of the inversions are concentrated at stronger fields than those from the simulation. Our analysis provides a quantification of the errors associated with the inversions of the Ca II 8542 \AA\ line and suggests caution with the interpretation of the inferred magnetic field fluctuations.

We use a science-grade Skipper Charge Coupled Device (Skipper-CCD) operating in a low-radiation background environment to develop a semi-empirical model that characterizes the origin of single-electron events in CCDs. We identify, separate, and quantify three independent contributions to the single-electron events, which were previously bundled together and classified as ``dark counts'': dark current, amplifier light, and spurious charge. We measure a dark current, which depends on exposure, of (5.89+-0.77)x10^-4 e-/pix/day, and an unprecedentedly low spurious charge contribution of (1.52+-0.07)x10^-4 e-/pix, which is exposure-independent. In addition, we provide a technique to study events produced by light emitted from the amplifier, which allows the detector's operation to be optimized to minimize this effect to a level below the dark-current contribution. Our accurate characterization of the single-electron events allows one to greatly extend the sensitivity of experiments searching for dark matter or coherent neutrino scattering. Moreover, an accurate understanding of the origin of single-electron events is critical to further progress in ongoing R&D efforts of Skipper and conventional CCDs.

Parametric stochastic simulators are ubiquitous in science, often featuring high-dimensional input parameters and/or an intractable likelihood. Performing Bayesian parameter inference in this context can be challenging. We present a neural simulator-based inference algorithm which simultaneously offers simulation efficiency and fast empirical posterior testability, which is unique among modern algorithms. Our approach is simulation efficient by simultaneously estimating low-dimensional marginal posteriors instead of the joint posterior and by proposing simulations targeted to an observation of interest via a prior suitably truncated by an indicator function. Furthermore, by estimating a locally amortized posterior our algorithm enables efficient empirical tests of the robustness of the inference results. Such tests are important for sanity-checking inference in real-world applications, which do not feature a known ground truth. We perform experiments on a marginalized version of the simulation-based inference benchmark and two complex and narrow posteriors, highlighting the simulator efficiency of our algorithm as well as the quality of the estimated marginal posteriors. Implementation on GitHub.

We discuss realization of cosmic inflation in the $\nu$-gauge mediated supersymmetry breaking scenario, in which a set of 24-dimensional chiral superfields responsible for the type III seesaw mechanism play the role of the messenger fields in gauge mediation. Using the data from neutrino oscillations, we show that the model satisfies constraints from the lepton flavor violation, perturbativity of the unified gauge couplings, the observed abundance of dark matter as well as the Higgs mass of 125.1 GeV. The predicted spectrum of the cosmic microwave background radiation fits well with the observation. We also comment on the falsifiability of this scenario by future experiments.

The formation and presence of clathrate hydrates could influence the composition and stability of planetary ices and comets; they are at the heart of the development of numerous complex planetary models, all of which include the necessary condition imposed by their stability curves, some of which include the cage occupancy or host-guest content and the hydration number, but fewer take into account the kinetics aspects. We measure the temperature-dependent-diffusion-controlled formation of the carbon dioxide clathrate hydrate in the 155-210~K range in order to establish the clathrate formation kinetics at low temperature. We exposed thin water ice films of a few microns in thickness deposited in a dedicated infrared transmitting closed cell to gaseous carbon dioxide maintained at a pressure of a few times the pressure at which carbon dioxide clathrate hydrate is thermodynamically stable. The time dependence of the clathrate formation was monitored with the recording of specific infrared vibrational modes of CO2 with a Fourier Transform InfraRed (FTIR) spectrometer. These experiments clearly show a two-step clathrate formation, particularly at low temperature, within a relatively simple geometric configuration. We satisfactorily applied a model combining surface clathration followed by a bulk diffusion-relaxation growth process to the experiments and derived the temperature-dependent-diffusion coefficient for the bulk spreading of clathrate. The derived apparent activation energy corresponding to this temperature-dependent-diffusion coefficient in the considered temperature range is E_a = 24.7 +/- 9.7 kJ/mol. The kinetics parameters favour a possible carbon dioxide clathrate hydrate nucleation mainly in planets or satellites.

We analyze the properties of the circular orbits for massive particles in the equatorial plane of symmetric rotating Ellis wormholes. In particular, we obtain the orbital frequencies and the radial and vertical epicyclic frequencies, and consider their lowest parametric, forced and Keplerian resonances. These show that quasi-periodic oscillations in accretion disks around symmetric rotating Ellis wormholes have many distinct properties as compared to quasi-periodic oscillations in accretion disks around rotating Teo wormholes and the Kerr black hole. Still we can distinguish some common features which appear in wormhole spacetimes as opposed to black holes. The most significant ones include the possibility of excitation of stronger resonances such as lower order parametric and forced resonances and the localization of these resonances deep in the region of strong gravitational interaction near the wormhole throat, which will lead to further amplification of the signal.

We present a non-linear analysis of perturbations around cosmological solutions in Generalised Massive gravity. This Lorentz invariant theory is an extension of de Rham, Gabadadze, Tolley massive gravity that propagates $5$ degrees of freedom while allowing stable open FLRW cosmologies. For a minimal model that supports a self-accelerating background, we study the dynamics of non-linear perturbations. We find that the equation for the scalar graviton is distinct from the analogues in Horndeski and DHOST theories. We numerically solve the equation to find a new type of nonlinear solution for the scalar mode, and confirm the presence of a Vainshtein screening mechanism. We show that the PPN parameter approaches its GR value at solar system scales and satisfies the current bounds.

We discuss the implications of Affleck-Dine (AD) baryogenesis for different classes of baryon and lepton number violating processes: specially focussing on implications for neutron-anti-neutron ($n-\bar{n}$) oscillation. The class of AD baryogenesis scenarios we work with uses the AD field also as the inflaton which is nonminimally coupled to gravity. We find that adequate baryogenesis and no washout by the baryon number ($B$) or the lepton number ($L$) violating operators implies constraints on the observability of the process or in the case of neutrino mass with compatibility with neutrino oscillation observations. In particular for $n-\bar{n}$ oscillation, we study some of the familiar operators that connect the AD field to $n-\bar{n}$ oscillation and find that a split scalar spectrum model turns out to be most advantageous for obtaining an observable $n-\bar{n}$ while remaining consistent with AD baryogenesis. It is interesting that this spectrum is similar to a non-supersymmetyric SO(10) model for observable $n-\bar{n}$ oscillation discussed before, suggesting that this AD scenario can be embedded into a grand unified SO(10) model. A feature of this baryogenesis scenario for $n-\bar{n}$ oscillation is that it necessarily predicts processes with $\Delta B=4$ or higher, all be it with highly suppressed amplitudes.

We describe and present the first observational evidence that light propagating near a rotating black hole is twisted in phase and carries orbital angular momentum. The novel use of this physical observable as an additional tool for the previously known techniques of gravitational lensing allows us to directly measure, for the first time, the spin parameter of a black hole. With the additional information encoded in the orbital angular momentum, not only can we reveal the actual rotation of the compact object, but we can also use rotating black holes as probes to test General Relativity.

In the present paper we show the existence of a fully nonlinear dynamical mechanism for the formation of scalarized black holes which is different from the spontaneous scalarization. We consider a class of scalar-Gauss-Bonnet gravity theories within which no tachyonic instability can occur. Although the Schwarzschild black holes are linearly stable against scalar perturbations, we show dynamically that for certain choices of the coupling function they are unstable against nonlinear scalar perturbations. This nonlinear instability leads to the formation of new black holes with scalar hair. The fully nonlinear and self-consistent study of the equilibrium black holes reveals that the spectrum of solutions is more complicated possessing additional branches with scalar field that turn out to be unstable, though. The formation of scalar hair of the Schwarzschild black hole will always happen with a jump because the stable "scalarized branch" is not continuously connected to Schwarzschild one.

LIGO is considered the most sensitive and complicated gravitational experiment ever built. Its main objective is to detect the gravitational wave from the strongest events in the universe by observing if the length of its 4-kilometer arms change by a distance 10,000 times smaller than the diameter of a proton. Due to its sensitivity, LIGO is prone to the disturbance of external noises which affects the data being collected to detect the gravitational wave. These noises are commonly called by the LIGO community as glitches. The objective of this study is to evaluate the effeciency of various deep trasnfer learning models namely VGG19, ResNet50V2, VGG16 and ResNet101 to detect glitch waveform in gravitational wave data. The accuracy achieved by the said models are 98.98%, 98.35%, 97.56% and 94.73% respectively. Even though the models achieved fairly high accuracy, it is observed that all of the model suffered from the lack of data for certain classes which is the main concern found in the experiment

The "dilemma" posed by the recent announcement of the PREX-II measurement of the neutron skin of $^{208}$Pb that suggests a stiff symmetry energy near nuclear matter density $n_0$ and its impact on the EoS of massive compact stars raise the issue as to whether the EoS determined at $n_0$ necessarily gives a stringent constraint at high densities relevant to massive compact stars, the currently widely accepted "lore" in astro-nuclear field. We present the argument that the "cusp" structure in the symmetry energy at $n_{1/2}\gtrsim 2 n_0$ predicted by a topology change in dense matter encoding hadron-quark continuity could provide a strong obstruction to the validity of the "lore." The topology change predicts an EoS that is soft below and stiffened above $n\gtrsim n_{1/2}$, involving no phase transitions, and yields the macrophysical properties of neutron stars more or less consistent with the overall observations and the maximum mass $ 2.0\lesssim M/ M_\odot\lesssim 2.2$ as well as the GW data. Furthermore it describes the interior core of the massive stars constituted of baryon-charge-fractionalized quasifermions, that are neither baryonic nor quarkonic, with the "pseudo-conformal" sound speed $v^2_{pcs}/c^2\approx 1/3$ with a nonzero trace of energy-momentum tensor converging from below at $n_{1/2}$.

The high cadence plasma, electric, and magnetic field measurements by the Magnetospheric MultiScale spacecraft allow us to explore the near-Earth space plasma with an unprecedented time and spatial resolution, resolving electron-scale structures that naturally emerge from plasma complex dynamics. The formation of small-scale turbulent features is often associated to structured, non-Maxwellian particle velocity distribution functions that are not at thermodynamic equilibrium. Using measurements in the terrestrial magnetosheath, this study focuses on regions presenting bumps in the power spectral density of the parallel electric field at sub-proton scales. Correspondingly, it is found that the ion velocity distribution functions exhibit beam-like features at nearly the local ion thermal speed. Ion cyclotron waves in the ion-scale range are frequently observed at the same locations. These observations, supported by numerical simulations, are consistent with the generation of ion-bulk waves that propagate at the ion thermal speed. This represents a new branch of efficient energy transfer at small scales, which may be relevant to weakly collisional astrophysical plasmas.

Primordial black holes as a viable candidate of cold dark matter are clustering in the large-scale structures of the Universe. The stochastic gravitational-wave background, arising from an incoherent superposition of unresolved gravitational waves from primordial black hole binaries, is expected to display anisotropies across the sky. In this work, we obtain the angular power spectrum for such anisotropies for the first time and the difference between it and the one from astrophysical black hole binaries. When multi-band measurements are employed, the angular power spectrum can be a smoking gun for the existence of primordial black holes, in particular those of stellar masses that otherwise are challenging to be distinguished from the astrophysical black holes.

We investigate the possibility of gravitationally generated particle production via the mechanism of nonminimal torsion--matter coupling. An intriguing feature of this theory is that the divergence of the matter energy--momentum tensor does not vanish identically. We explore the physical and cosmological implications of the nonconservation of the energy--momentum tensor by using the formalism of irreversible thermodynamics of open systems in the presence of matter creation/annihilation. The particle creation rates, pressure, and the expression of the comoving entropy are obtained in a covariant formulation and discussed in detail. Applied together with the gravitational field equations, the thermodynamics of open systems lead to a generalization of the standard $\Lambda$CDM cosmological paradigm, in which the particle creation rates and pressures are effectively considered as components of the cosmological fluid energy--momentum tensor. We consider specific models, and we show that cosmology with a torsion--matter coupling can almost perfectly reproduce the $\Lambda$CDM scenario, while it additionally gives rise to particle creation rates, creation pressures, and entropy generation through gravitational matter production in both low and high redshift limits.