Papers for Tuesday, Jun 06 2023

The properties of low-mass dark matter halos appear to be remarkably diverse relative to predictions of cold, collisionless dark matter, even in the presence of baryons. We show that self-interacting dark matter can simultaneously explain two extreme measurements of halo diversity in different directions -- namely, the rotation curves of low-concentration halos associated with gas-rich ultra-diffuse galaxies in the field and the inner density profile of the dense substructure perturbing the strong lens galaxy SDSSJ0946+1006. We present the first cosmological zoom-in simulation featuring strong dark matter self-interactions in a galaxy group environment centered on a $10^{13}~M_{\mathrm{\odot}}$ host halo. These interactions produce kiloparsec-scale cores in low-concentration isolated halos, which could host the ultra-diffuse galaxies, while most surviving subhalos of the group-mass host are deeply core-collapsed, yielding excellent candidates for the observed dense strong-lens perturber. Our scenario can be further tested with observations of galactic systems over a wide mass range.

Supermassive black holes (SMBHs) reside at the center of every massive galaxy in the local Universe with masses that closely correlate with observations of their host galaxy implying a connected evolutionary history. The population of binary SMBHs, which form following galaxy mergers, is expected to produce a gravitational wave background (GWB) detectable by pulsar timing arrays (PTAs). PTAs are starting to see hints of what may be a GWB, and the amplitude of the emerging signal is towards the higher end of model predictions. Simulated populations of binary SMBHs can be constructed from observations of galaxies and are used to make predictions about the nature of the GWB. The greatest source of uncertainty in these observation-based models comes from the inference of the SMBH mass function, which is derived from observed host galaxy properties. In this paper, I undertake a new approach for inferring the SMBH mass function starting from a velocity dispersion function rather than a galaxy stellar mass function. I argue that this method allows for a more direct inference by relying on a larger suite of individual galaxy observations as well as relying on a more "fundamental" SMBH mass relation. I find that the resulting binary SMBH population contains more massive systems at higher redshifts than previous models. Additionally, I explore the implications for the detection of individually resolvable sources in PTA data.

The interaction between the ejecta of supernovae (SNe) of Type IIn and a dense circumstellar medium (CSM) can efficiently generate thermal UV/optical radiation and lead to the emission of neutrinos in the $1$-$10^{3}$ TeV range. We investigate the connection between the neutrino signal detectable at the IceCube Neutrino Observatory and the electromagnetic signal observable by optical wide-field, high-cadence surveys to outline the best strategy for upcoming follow-up searches. We outline a semi-analytical model that connects the optical lightcurve properties to the SN parameters and find that a large peak luminosity ($ L_{\rm{peak}}\gtrsim 10^{43}$-$10^{44}$ erg) and an average rise time ($t_{\rm{rise}}\gtrsim 10$-$40$ days) are necessary for copious neutrino emission. Nevertheless, the most promising $ L_{\rm{peak}}$ and $t_{\rm{rise}}$ are not sufficient to guarantee ideal conditions for neutrino detection. Comparable optical properties can be obtained for SN configurations that are not optimal for neutrino emission. Such ambiguous correspondence between the optical lightcurve properties and the number of IceCube neutrino events implies that relying on optical observations only, a range of expected neutrino events should be considered (e.g. the expected number of neutrino events can vary up to two orders of magnitude for some among the brightest SNe IIn observed by the Zwicky Transient Facility up to now, SN 2020usa and SN 2020in). In addition, the peak in the high-energy neutrino curve should be expected a few $t_{\rm{rise}}$ after the peak in the optical lightcurve. Our findings highlight that it is crucial to infer the SN properties from multi-wavelength observations rather than focusing on the optical band only to enhance upcoming neutrino searches.

K- and L-type asteroids are considered to be the parent bodies of CV and CO chondrites. Spectral models of L-types invoke an enrichment in CAI with respect to the chondrites in the meteorite collection. Barbarian asteroids are associated to L-type asteroids yet the relationship between these populations is still not clear. We aim to investigate the link between the K- and L-type and Barbarian asteroids and the CV and CO chondrites by means of spectral matching of a large number of reflectance spectra of objects from either population. We seek to identify matches based on observed rather than modelled spectral features. We employ a matching criterion that accounts for the residuals and the correlation of the spectral features. The only free parameter in the comparison is the degree of alteration of the asteroids with respect to the meteorites expressed via an exponential model. We derive an absolute scale of similarity between the spectra using laboratory data from irradiation experiments. CVOxA chondrites are the best match to the asteroids, in particular to K-type (7 out of 11 asteroids matched) and Barbarians (11 out of 16). CO chondrites provide convincing matches for K-types (5 out of 11) and Barbarians (7 out of 16) as well. A single non-Barbarian L-type is matched to a meteorite. Only a few asteroids are matched to CVOxB and CVRed chondrites. Barbarian asteroids are represented among CO and CVOxA chondrites without requiring an enrichment of CAI in the asteroids. Four candidate Barbarian asteroids are identified, three of which are classified as K-types. These asteroids are favourable targets for polarimetric observations. The discrepancy between L-type asteroids and CV and CO chondrites is likely related to the ambiguity of the asteroid class itself. An extension of the taxonomy to include polarimetric properties is required.

We generalize the recently proposed Stepped Partially Acoustic Dark Matter (SPartAcous) model by including additional massless degrees of freedom in the dark radiation sector. We fit SPartAcous and its generalization against cosmological precision data from the cosmic microwave background, baryon acoustic oscillations, large-scale structure, supernovae type Ia, and Cepheid variables. We find that SPartAcous significantly reduces the $H_0$ tension but does not provide any meaningful improvement of the $S_8$ tension, while the generalized model succeeds in addressing both tensions, and provides a better fit than $\Lambda\mathrm{CDM}$ and other dark sector models proposed to address the same tensions. In the generalized model, $H_0$ can be raised to $71.4~\mathrm{km/s/Mpc}$ (the 95% upper limit) if the fitted data does not include the direct measurement from the SH0ES collaboration, and to $73.7~\mathrm{km/s/Mpc}$ (95% upper limit) if it does. A version of $\texttt{CLASS}$ that has been modified to analyze this model is publicly available at $\href{https://github.com/ManuelBuenAbad/class_spartacous}{\texttt{github.com/ManuelBuenAbad/class_spartacous}}$.

The activity of at least one repeating Fast Radio Burst (FRB) source is periodically modulated. If this modulation is the result of precession of the rotation axis and throat of an accretion disc around a black hole, driven by a companion that is also the source of accreted mass, then it may be possible to constrain the mass of the black hole. The dynamics is analogous to that of superorbital periods in ordinary mass-transfer binaries in which the accreting object may be a stellar-mass black hole, a neutron star or a white dwarf, but in the FRB source it may be an intermediate mass black hole. In a semi-detached (mass-transferring) binary the orbital period is related to the mean density of the mass-losing star. Assuming a value for its density and identifying the observed modulation period as a disc precession period would determine the mass ratio and the mass of the black hole. This model and magnetar-SNR models make distinguishable predictions of the evolution of the FRB rotation measure that may soon be tested in FRB 121102.

Solar type III radio bursts are generated by beams of energetic electrons that travel along open magnetic field lines through the corona and into interplanetary space. However, understanding the source of these electrons and how they escape into interplanetary space remains an outstanding topic. Here we report multi-instrument, multi-perspective observations of an interplanetary type III radio burst event shortly after the second perihelion of the Parker Solar Probe (PSP). This event was associated with a solar jet that produced an impulsive microwave burst event recorded by the Expanded Owens Valley Solar Array (EOVSA). The type III burst event also coincided with the detection of enhanced in situ energetic electrons recorded by both PSP at 0.37 AU and WIND at 1 AU, which were located very closely on the Parker spiral longitudinally. The close timing association and magnetic connectivity suggest that the in situ energetic electrons originated from the jet's magnetic reconnection region. Intriguingly, microwave imaging spectroscopy results suggest that the escaping energetic electrons were injected into a large opening angle of about 90 degrees, which is at least nine times broader than the apparent width of the jet spire. Our findings provide an interpretation for the previously reported, longitudinally broad spatial distribution of flare locations associated with prompt energetic electron events and have important implications for understanding the origin and distribution of energetic electrons in interplanetary space.

We analyze the directional dependence of the gravitational wave (GW) emission from 15 3D neutrino radiation hydrodynamic simulations of core-collapse supernovae. We develop a new analytic technique to characterize the distribution of GW emission over all angles. We use physics-informed toy models to provide closed form expressions for the distribution of GW emission for different CCSN phases. Using these toy models, we approximate the PNS dynamics during multiple CCSN stages and obtain similar GW distributions to simulation outputs. By applying this new technique throughout the supernova duration, we construct a distribution of preferred directions of GW emission. Our findings indicate CCSNe do not have a single `optimal' viewing angle along which the strongest GWs can be detected. For nonrotating cases, this dominant viewing angle drifts isotropically throughout the supernova, set by the dynamical timescale of the protoneutron star. For rotating cases, during core bounce and the following tens of ms, the strongest GW signal is observed along the equator. During the accretion phase, comparable -- if not stronger -- GW amplitudes are generated along the axis of rotation, which can be enhanced by the low T/|W| instability. We show two dominant factors influencing the directionality of GW emission are the degree of initial rotation and explosion morphology. Lastly, looking forward, we note the sensitive interplay between GW detector site and supernova orientation, along with its effect on detecting individual polarization modes.

The distribution of small-scale magnetic fields in stellar photospheres is an important ingredient in our understanding of the magnetism of low mass stars. Their spatial distribution connects the field generated in the stellar interior with the outer corona and the large scale field, and thereby affects the space weather of planets. Unfortunately, we lack techniques that can locate them on most low-mass stars. One strategy is to localize field concentrations using the flares that occur in their vicinity. We explore a new method that adapts the spot simulation software fleck to study the modulation of flaring times as a function of active latitude. We use empirical relations to construct flare light curves similar to those available from Kepler and the Transiting Exoplanet Survey Satellite (TESS), search them for flares, and use the waiting times between flares to determine the location of active latitudes. We find that the mean and standard deviation of the waiting time distribution provide a unique diagnostic of flaring latitudes as a function of the number of active regions. Latitudes are best recovered when stars have three or less active regions that flare repeatedly, and active latitude widths below 20 deg; when either increases, the information about the active latitude location is gradually lost. We demonstrate our technique on a sample of flaring G dwarfs observed with the Kepler satellite, and furthermore suggest that combining ensemble methods for spots and flares could overcome the limitations of each individual technique for the localization of surface magnetic fields.

Understanding the astrophysical nature of the first stars still remains an unsolved problem in cosmology. The redshifted global 21-cm signal and power spectrum act as a treasure trove to probe the Cosmic Dawn era -- when the intergalactic medium was mostly neutral. Many experiments, like SARAS 3, SKA, EDGES, DARE, etc., have been proposed to probe the cosmic dawn era. However, extracting the faint cosmological signal buried inside the brighter foregrounds $\mathcal{O}(10^4)$ remains challenging. Considering the excess radio background, we have constructed all possible $T_{21}$ signals in the EDGES limit. We have used a single Artificial Neural Network for $T_{21}$ parameter extraction in the presence of the foreground and noise with Root Mean Square Error (RMSE) and R-Squared (R2) score of $(0.2 - 0.08)$ and $(0.66 - 0.94)$, respectively. Here, we also explore the parameter estimation in the presence of heating of intergalactic medium due to background radio radiation mediated by Ly$\alpha$ photons from first stars, and we found that the effect indeed has a significant impact on parameters correlation and their estimation.

The orbital distribution of the S-star cluster surrounding the supermassive black hole in the center of the Milky Way is analyzed. A tight, roughly exponential dependence of the pericenter distance r$_{p}$ on orbital eccentricity e$_{\star}$ is found, $\log ($r$_p)\sim$(1-e$_{\star}$), which cannot be explained simply by a random distribution of semi-major axes and eccentricities. No stars are found in the region with high e$_{\star}$ and large log r$_{p}$ or in the region with low e$_{\star}$ and small log r$_{p}$. G-clouds follow the same correlation. The likelihood P(log r$_p$,(1-e$_{\star}$)) to determine the orbital parameters of S-stars is determined. P is very small for stars with large e$_{\star}$ and large log r$_{p}$. S-stars might exist in this region. To determine their orbital parameters, one however needs observations over a longer time period. On the other hand, if stars would exist in the region of low log r$_{p}$ and small e$_{\star}$, their orbital parameters should by now have been determined. That this region is unpopulated therefore indicates that no S-stars exist with these orbital characteristics, providing constraints for their formation. We call this region, defined by $\log$ (r$_p$/AU) $<$ 1.57+2.6(1-e$_{\star})$, the zone of avoidance. Finally, it is shown that the observed frequency of eccentricities and pericenter distances is consistent with a random sampling of log r$_{p}$ and e$_{\star}$. However, only if one takes into account that no stars exist in the zone of avoidance and that orbital parameters cannot yet be determined for stars with large r$_{p}$ and large e$_{\star}$.

In this brief review I summarize our basic knowledge about different types of isolated neutron stars. I discuss radio pulsars, central compact objects in supernova remnants, magnetars, near-by cooling neutron stars (aka the Magnificent seven), and sources of fast radio bursts. Several scenarios of magneto-rotational evolution are presented. Recent observational data, in the first place -- discovery of long period radio pulsar, require non-trivial evolution of magnetics fields or/and spin periods of neutron stars. In some detail I discuss different models of magnetic field decay and interaction of young neutron stars with fallback matter.

We report on \ixpe observations of the Be-transient X-ray pulsar LS V +44 17/RX J0440.9+4431 at two luminosity levels during the giant outburst in January--February 2023. Considering the observed spectral variability and changes in the pulse profiles, the source was likely caught in super- and sub-critical states with significantly different emission region geometry, associated with the presence of accretion columns and hot spots, respectively. We focus here on the pulse-phase resolved polarimetric analysis and find that the observed dependencies of the polarization degree and polarization angle (PA) on pulse phase are indeed drastically different for the two observations. The observed differences, if interpreted within the framework of the rotating vector model (RVM), imply dramatic variations of the spin axis inclination and the position angle and the magnetic colatitude by tens of degrees within just a few days separating the observations. We suggest that the apparent changes in the observed PA phase dependence are predominantly related to the presence of a polarized unpulsed component in addition to the polarized radiation associated with the pulsar itself. We show that the observed PA phase dependence in both observations can then be explained with a single set of RVM parameters defining the pulsar's geometry. We also suggest that the additional polarized component is likely produced by scattering of the pulsar radiation off the equatorial disk wind.

The \textit{absolute age} of a simple stellar population is of fundamental interest for a wide range of applications but is difficult to measure in practice, as it requires an understanding of the uncertainties in a variety of stellar evolution processes as well as the uncertainty in the distance, reddening and composition. As a result, most studies focus only on the \textit{relative age} by assuming that stellar evolution calculations are accurate and using age determinations techniques that are relatively independent of distance and reddening. Here, we construct $20,000$ sets of theoretical isochrones through Monte Carlo simulation using the Dartmouth Stellar Evolution Program to measure the absolute age of the globular cluster M92. For each model, we vary a range of input physics used in the stellar evolution models, including opacities, nuclear reaction rates, diffusion coefficients, atmospheric boundary conditions, helium abundance, and treatment of convection. We also explore variations in the distance and reddening as well as its overall metallicity and $\alpha$ enhancement. We generate simulated Hess diagrams around the main-sequence turn-off region from each set of isochrones and use a Voronoi binning method to fit the diagrams to HST ACS data. We find the age of M92 to be $13.80 \pm 0.75$ Gyr. The $5.4\%$ error in the absolute age is dominated by the uncertainty in the distance to M92 ($\sim 80\%$ of the error budget); of the remaining parameters, only the total metallicity, $\alpha$ element abundance, and treatment of helium diffusion contribute significantly to the total error.

In our previous work, we searched for super-flares on different types of stars while focusing on G-type dwarfs using entire Kepler data to study statistical properties of the occurrence rate of super-flares. Using these new data, as a by-product, we found fourteen cases of super-flare detection on thirteen slowly rotating Sun-like stars with rotation periods of 24.5 to 44 days. This result supports earlier conclusion by others that the Sun may possibly have a surprise super-flare. Moreover, we found twelve and seven new cases of detection of exceptionally large amplitude super-flares on six and four main-sequence stars of G- and M-type, respectively. No large-amplitude flares were detected in A, F, or K main-sequence stars. Here we present preliminary analysis of these cases. The super-flare detection, i.e. an estimation of flare energy, is based on a more accurate method compared to previous studies. We fit an exponential decay function to flare light curves and study the relation between e-folding decay time, $\tau$, vs. flare amplitude and flare energy. We find that for slowly rotating Sun-like stars, large values of $\tau$ correspond to small flare energies and small values of $\tau$ correspond to high flare energies considered. Similarly, $\tau$ is large for small flare amplitudes and $\tau$ is small for large amplitudes considered. However, there is no clear relation between these parameters for large amplitude super-flares in the main sequence G- and M-type stars, as we could not establish clear functional dependence between the parameters via standard fitting algorithms.

We present a refined analysis of 15038 Kepler main sequence light curves to determine the stellar rotation periods. The initial period estimates come from an autocorrelation function, as has been done before. We then measure the duration of every intensity dip in the light curve, expressed as fractions of the initial rotation period estimate. These dip duration distributions are subdivided into several regions whose relation to each other helps determine which harmonic of the initial rotation period is most physically plausible. We compare our final rotation periods to those from McQuillan, Mazeh, & Aigrain (2014) and find that the great majority agree, but about 10% of their periods are doubtful (usually twice as long as is most plausible). We are still refining our method, and will later extend it to more stars to substantially increase the sample of reliable stellar rotation periods.

In the solar corona, magnetic reconnection occurs due to the finite resistivity of the plasma. At the same time, resistivity leads to ohmic heating. Therefore, the reconnecting current sheet should heat the surrounding plasma. This paper presents experimental evidence of such plasma heating caused by magnetic reconnection. We observed the effect during a C1.4 solar flare on 16 February 2003 at the active region NOAA 10278, near the solar limb. Thanks to such a location, we successfully identified all the principal elements of the flare: the flare arcade, the fluxrope, and, most importantly, the presumed position of the current sheet. By analyzing the monochromatic X-ray images of the Sun obtained by the CORONAS-F/SPIRIT instrument in the Mg XII 8.42 A spectral line, we detected a high-temperature ($T \geq$ 4 MK) emission at the predicted location of the current sheet. The high-temperature emission appeared during the CME impulsive acceleration phase. We believe that this additionally confirms that the plasma heating around the current sheet and magnetic reconnection inside the current sheet are strongly connected.

The tilt of the velocity ellipsoid is a helpful tracer of the gravitational potential of the Milky Way. In this paper, we use nearly 140,000 RC stars selected from the LAMOST and {\it Gaia} to make a detailed analysis of the tilt of the velocity ellipsoid for various populations, as defined by the stellar ages and chemical information, within 4.5\,$\leq$\,$R$\,$\leq$\,15.0\,kpc and $|Z|$\,$\leq$\,3.0\,kpc. The tilt angles of the velocity ellipsoids of the RC sample stars are accurately described as $\alpha$ = $\alpha_{0}$\,$\mathrm{arctan}$\,($Z$/$R$) with $\alpha_{0}$ = (0.68 $\pm$ 0.05). This indicates the alignment of velocity ellipsoids is between cylindrical and spherical, implying that any deviation from the spherical alignment of the velocity ellipsoids may be caused by the gravitational potential of the baryonic disk. The results of various populations suggest that the $\alpha_{0}$ displays an age and population dependence, with the thin and thick disks respectively values $\alpha_{0}$ = (0.72 $\pm$ 0.08) and $\alpha_{0}$ = (0.64 $\pm$ 0.07), and the $\alpha_{0}$ displays a decreasing trend with age (and [$\alpha$/Fe]) increases, meaning that the velocity ellipsoids of the kinematically relaxed stars are mainly dominated by the gravitational potential of the baryonic disk. We determine the $\alpha_{0}$ -- $R$ for various populations, finding that the $\alpha_{0}$ displays oscillations with $R$ for all the different populations. The oscillations in $\alpha_{0}$ appear in both kinematically hot and cold populations, indicating that resonances with the Galactic bar are the most likely origin for these oscillations.

The effect of dynamical tide ``kicks" on eccentric binary orbits is considered using the orbital mapping method. It is demonstrated that when mode damping is negligible the mode amplitude will generically grow in time for all values of orbital eccentricity and semi-major axis, even for small kicks outside the regime exhibiting diffusive growth. The origin of the small-kick growth is the change in kick size from orbit to orbit, an effect quadratic in the mode amplitude. When damping of the mode is included, the growth is shut off when the damping time is shorter than the growth time. Hence, in practice, kicks of sufficient size and long mode damping times are required for interesting levels of growth to occur. Application to the circularization of hot Jupiters is discussed. Previous investigations found that diffusive growth of the planetary f-mode in the large-kick regime would lead to rapid orbital shrinkage, but upon exiting the diffusive regime at $e \sim 0.9$ the theory would predict a large population of highly eccentric orbits. Simulations presented here show that subsequent orbital evolution relying on the small-kick regime may further decrease the eccentricity to $e \sim 0.2$ on timescales much less than the Gyrs ages of these systems.

A detailed analysis is presented of the gravitational microlensing by intervening compact objects of the black hole shadows imaged by the Event Horizon Telescope (EHT). We show how the center, size, and shape of the shadow depend on the Einstein angle relative to the true/unlensed shadow size, and how the location of the lens affects the shift, size, and asymmetry of the black hole shadow due to microlensing. Assuming a supermassive black hole (SMBH) casts a circular-shaped true shadow, microlensing can create an asymmetry of up to approximately 8\%, which is twice the asymmetry caused by the SMBH's spin and its tilt relative to us. Furthermore, the size can be enhanced by $\sim$50\% of the true shadow. Currently, the terrestrial baselines of EHT lack the resolution to detect microlensing signatures in the shadows. However, future expansions of EHT including space-based baselines at the Moon and L$_2$, could potentially enable the detection of microlensing events. For Sgr~A$^*$, an event rate of 0.0014 per year makes the microlensing phenomena difficult to observe even with space-based baselines for the stellar population in the stellar bulge and stellar disk for lens mass $\sim M_\odot$. However, continuously monitoring the shadow of Sgr~A$^*$ could offer novel insights into the compact object population surrounding the galactic center.

We present the latest results of our ongoing multiplicity study of (Community) TESS Objects of Interest, using astrometric and photometric data from the ESA-Gaia mission to detect stellar companions of these stars and characterize their properties.\linebreak A total of 134 binary, 6 hierarchical triple, and two quadruple star systems are identified among 1106 targets whose multiplicity is investigated in the course of our survey, located at distances closer than about 500pc around the Sun. The detected companions and targets are at the same distance and have a common proper motion, as expected for components of gravitationally bound stellar systems, as demonstrated by their accurate Gaia DR3 astrometry. The companions have masses from about 0.11 to 2$M_\odot$ and are most abundant in the mass range between 0.2 and 0.5$M_\odot$. The companions have projected separations from the targets between about 50 and 9700au. Their frequency is the highest and constant from about 300 up to 750au, decreasing at larger projected separations. In addition to main sequence stars, four white dwarf companions are detected in this study, whose true nature is revealed by their photometric properties.

We examine the effect of thermal conduction on the low-angular momentum hot accretion flow (HAF) around non-rotating black holes accreting mass at very low rate. While doing so, we adopt the conductive heat flux in the saturated form, and solve the set of dynamical equations corresponding to a steady, axisymmetric, viscous, advective accretion flow using numerical methods. We study the dynamical and thermodynamical properties of accreting matter in terms of the input parameters, namely energy ($\varepsilon_0$), angular momentum ($\ell_0$), viscosity parameter ($\alpha$), and saturation constant ($\Phi_{\rm s}$) regulating the effect of thermal conduction. We find that $\Phi_{\rm s}$ plays a pivotal role in deciding the transonic properties of the global accretion solutions. In general, when $\Phi_{\rm s}$ is increased, the critical point ($r_{\rm c}$) is receded away from the black hole, and flow variables are altered particularly in the outer part of the disc. To quantify the physically acceptable range of $\Phi_{\rm s}$, we compare the global transonic solutions with the self-similar solutions, and observe that the maximum saturation constant ($\Phi^{\rm max}_{\rm s}$) estimated from the global solutions exceeds the saturated thermal conduction limit ($\Phi_{\rm sc}$) derived from the self-similar formalism. Moreover, we calculate the correlation between $\alpha$ and $\Phi^{\rm max}_{\rm s}$ and find ample disagreement between global solutions and self-similar solutions. Further, using the global flow variables, we compute the Bernoulli parameter ($Be$) which remains positive all throughout the disc, although flow becomes loosely unbound for higher $\Phi_{\rm s}$. Finally, we indicate the relevance of this work in the astrophysical context in explaining the possibility of massloss/outflows from the unbound disc.

Very-high-energy $\gamma$-ray emission provides constraints on the morphology and the physics mechanisms involved in the evolution of pulsar wind nebulae (PWNe). In the Galactic plane, around $312 ^{\circ}$ of Galactic longitude, a promising region two-degree wide containing five powerful pulsars may offer a new insight on the transition between TeV-emitting PWNe and pulsar halos. Their rotational energies range from $10^{35}$ to $10^{37}$ erg s$^{-1}$ for ages between 13.6 and 62.8 kyr. Extended emission is detected with H.E.S.S. (High Energy Stereoscopic System) in their vicinity, notably around the pulsar PSR J1413-6205. We processed 124 hours of H.E.S.S observations with an algorithm improving background fitting for the study of extended sources. We applied a three-dimensional likelihood analysis technique to model the different sources in the region using a configuration that optimizes the collection area at the highest energies. This contribution focuses on the detection of a new extended source around PSR J1413-6205 over 5$\sigma$ with a hard spectrum. Preliminary results on this source show a radius of $0.12 ^{\circ}$ $\pm$ $0.01 ^{\circ}_{\rm stat}$, an index of 2.06 $\pm$ 0.20$_{\rm stat}$ and a lower limit on a cut-off energy of 17 TeV, at a 90% confidence level. The detected emission is consistent with previous PWN models.

We present an overview of the James Webb Space Telescope (JWST) Advanced Deep Extragalactic Survey (JADES), an ambitious program of infrared imaging and spectroscopy in the GOODS-S and GOODS-N deep fields, designed to study galaxy evolution from high redshift to cosmic noon. JADES uses about 770 hours of Cycle 1 guaranteed time largely from the Near-Infrared Camera (NIRCam) and Near-Infrared Spectrograph (NIRSpec) instrument teams. In GOODS-S, in and around the Hubble Ultra Deep Field and Chandra Deep Field South, JADES produces a deep imaging region of ~45 arcmin$^2$ with an average of 130 hrs of exposure time spread over 9 NIRCam filters. This is extended at medium depth in GOODS-S and GOODS-N with NIRCam imaging of ~175 arcmin$^2$ with an average exposure time of 20 hrs spread over 8-10 filters. In both fields, we conduct extensive NIRSpec multi-object spectroscopy, including 2 deep pointings of 55 hrs exposure time, 14 medium pointings of ~12 hrs, and 15 shallower pointings of ~4 hrs, targeting over 5000 HST and JWST-detected faint sources with 5 low, medium, and high-resolution dispersers covering 0.6-5.3 microns. Finally, JADES extends redward via coordinated parallels with the JWST Mid-Infrared Instrument (MIRI), featuring ~9 arcmin$^2$ with 43 hours of exposure at 7.7 microns and twice that area with 2-6.5 hours of exposure at 12.8 microns For nearly 30 years, the GOODS-S and GOODS-N fields have been developed as the premier deep fields on the sky; JADES is now providing a compelling start on the JWST legacy in these fields.

JWST has revolutionized the field of extragalactic astronomy with its sensitive and high-resolution infrared view of the distant universe. Adding to the new legacy of JWST observations, we present the first NIRCam imaging data release from the JWST Advanced Deep Extragalactic Survey (JADES) providing 9 filters of infrared imaging of $\sim$25 arcmin$^2$ covering the Hubble Ultra Deep Field and portions of Great Observatories Origins Deep Survey (GOODS) South. Utilizing 87 on-sky dual-filter hours of exposure time, these images reveal the deepest ever near-infrared view of this iconic field. We supply carefully constructed 9-band mosaics of the JADES bands, as well as matching reductions of 5 additional bands from the JWST Extragalactic Medium-band Survey (JEMS). Combining with existing HST imaging, we provide 23-band space-based photometric catalogs and photometric redshifts for $\approx47,500$ sources. To promote broad engagement with the JADES survey, we have created an interactive {\tt FitsMap} website to provide an interface for professional researchers and the public to experience these JWST datasets. Combined with the first JADES NIRSpec data release, these public JADES imaging and spectroscopic datasets provide a new foundation for discoveries of the infrared universe by the worldwide scientific community.

We describe the NIRSpec component of the JWST Deep Extragalactic Survey (JADES), and provide deep spectroscopy of 253 sources targeted with the NIRSpec micro-shutter assembly in the Hubble Ultra Deep Field and surrounding GOODS-South. The multi-object spectra presented here are the deepest so far obtained with JWST, amounting to up to 28 hours in the low-dispersion ($R\sim 30-300$) prism, and up to 7 hours in each of the three medium-resolution $R\approx 1000$ gratings and one high-dispersion grating, G395H ($R\approx2700$). Our low-dispersion and medium-dispersion spectra cover the wavelength range $0.6-5.3\mu$m. We describe the selection of the spectroscopic targets, the strategy for the allocation of targets to micro-shutters, and the design of the observations. We present the public release of the reduced 2D and 1D spectra, and a description of the reduction and calibration process. We measure spectroscopic redshifts for 178 of the objects targeted extending up to $z=13.2$. We present a catalog of all emission lines detected at $S/N>5$, and our redshift determinations for the targets. Combined with the first JADES NIRCam data release, these public JADES spectroscopic and imaging datasets provide a new foundation for discoveries of the infrared universe by the worldwide scientific community.

We present a catalog of 717 candidate galaxies at $z > 8$ selected from 125 square arcminutes of NIRCam imaging as part of the JWST Advanced Deep Extragalactic Survey (JADES). We combine the full JADES imaging dataset with data from the JEMS and FRESCO JWST surveys along with extremely deep existing observations from HST/ACS for a final filter set that includes fifteen JWST/NIRCam filters and five HST/ACS filters. The high-redshift galaxy candidates were selected from their estimated photometric redshifts calculated using a template fitting approach, followed by visual inspection from seven independent reviewers. We explore these candidates in detail, highlighting interesting resolved or extended sources, sources with very red long-wavelength slopes, and our highest redshift candidates, which extend to $z_{phot} = 18$. We also investigate potential contamination by stellar objects, and do not find strong evidence from SED fitting that these faint high-redshift galaxy candidates are low-mass stars. Over 93\% of the sources are newly identified from our deep JADES imaging, including 31 new galaxy candidates at $z_{phot} > 12$. Using 42 sources in our sample with measured spectroscopic redshifts from NIRSpec and FRESCO, we find excellent agreement to our photometric redshift estimates, with no catastrophic outliers and an average difference of $\langle \Delta z = z_{phot}- z_{spec} \rangle= 0.26$. These sources comprise one of the most robust samples for probing the early buildup of galaxies within the first few hundred million years of the Universe's history.

We use SEDz*, a code designed to chart star formation histories of 6 < z < 12 galaxies, to analyze the SEDs of 982 galaxies with deep JWST-NIRCam imaging in the GOODS-S field. We show how SEDz* assembles an SED to match observations, from component stellar-population templates, graphing the contribution of each by epoch to confirm the robustness of the technique. Very good SED fits for most SFHs demonstrates the compatibility of our templates with generations of stars in these first galaxies, even though they are derived from present-epoch stars in our Galaxy -- as expected, because the light is dominated by main-sequence A-stars, free of post-main-sequence complexity and insensitive to heavy-element compositions. We confirm earlier results in Dressler+2023: (1) Four types of star formation histories (SFHs) -- burst, stochastic, 3-epoch `contiguous,' and longer `continuous' -- cover the variety of histories: (2) Starbursts - both single and multiple -- play the lead role in this critical period of cosmic history, even though less common, longer SFHs (0.5-1.0 Gyr) produce comparable stellar mass. This buildup of stellar mass amounts to 10^8 to 10^9 Msun in an integer epoch, whether in a ~100 Myr burst or as part of a longer history. We suggest that the absence of rising SFHs could be explained by an intense dust-enshrouded phase lasting tens of Myr that precedes their dust-free SFHs. With comparable contributions to the stellar mass by the four different SFH types, the stellar mass of these ~1000 galaxies amounts to ~2 x10^12 Msun by z = 6. We find no strong dependencies of the different SFH types with the large-scale environment. However, the discovery of a compact group of 16 galaxies, 9 of which had first star formation at z=11-12, suggests that long SFHs are much more common in rare, dense environments, likely destined for the most populated places of the modern universe.

We use deep NIRSpec spectroscopic data from the JADES survey to derive the star formation histories (SFHs) of a sample of 200 galaxies at 0.6$<$z$<$11 and spanning stellar masses from $\rm 10^6$ to $\rm 10^{9.5}~M_\odot$. We find that galaxies at high-redshift, galaxies above the Main Sequence (MS) and low-mass galaxies tend to host younger stellar populations than their low-redshift, massive, and below the MS counterparts. Interestingly, the correlation between age, M$_*$ and SFR existed even earlier than Cosmic Noon, out to the earliest cosmic epochs. However, these trends have a large scatter. Indeed, there are examples of young stellar populations also below the MS, indicating recent (bursty) star formation in evolved systems. We explore further the burstiness of the SFHs by using the ratio between SFR averaged over the last 10 Myr and averaged between 10 Myr and 100 Myr before the epoch of observation ($\mathrm{SFR_{cont, 10}/SFR_{cont, 90}}$). We find that high-redshift and low-mass galaxies have particularly bursty SFHs, while more massive and lower-redshift systems evolve more steadily. We also present the discovery of another (mini-)quenched galaxy at z = 4.4 (in addition to the one at z=7.3 reported by Looser et al. 2023), which might be only temporarily quiescent as a consequence of the extremely bursty evolution. Finally, we also find a steady decline of dust reddening of the stellar population approaching the earliest cosmic epochs, although some dust reddening is still observed in some of the highest redshift and most star forming systems.

The rest-frame UV recombination emission line Ly-alpha can be powered by ionizing photons from young massive stars in star forming galaxies, but its ability to be resonantly scattered by neutral gas complicates its interpretation. For reionization era galaxies, a neutral intergalactic medium (IGM) will scatter Ly-alpha from the line of sight, making Ly-alpha a useful probe of the neutral fraction at z>6. Here, we explore Ly-alpha in JWST/NIRSpec spectra from the ongoing GTO program JADES, which targets hundreds of galaxies in the well-studied GOODS-S and GOODS-N fields. These sources are UV-faint (-20.4<MUV<-16.4), and thus represent a poorly-explored class of galaxies. The low spectral resolution (R~100) spectra of a subset of 93 galaxies in GOODS-S with z_spec>5.5 (as derived with optical lines) are fit with line and continuum models, in order to search for significant line emission. Through exploration of the R100 data, we find evidence for Ly-alpha in 15 sources. Additional analysis of the R1000 data from the same set of galaxies results in five additional detections. This sample allows us to place observational constraints on the fraction of galaxies with Ly-alpha emission in the redshift range 5.5<z<7.5, with a decrease from z=6 to z=7. We also find a positive correlation between Ly-alpha equivalent width and MUV, as seen in other samples. These results are used to estimate the neutral gas fraction at z~7, agreeing with previous results (XHI~0.5-0.9).

The physical processes that establish the morphological evolution and the structural diversity of galaxies are key unknowns in extragalactic astrophysics. Here we report the finding of the morphologically-mature galaxy JADES-GS+53.18343-27.79097, which existed within the first 700 million years of the Universe's history. This star-forming galaxy with a stellar mass of $10^{8.6}$ solar masses consists of three components, a highly-compact core with a half-light radius of 144 pc, a strongly star-forming disc with a radius of 468 pc, and a star-forming clump, which all show distinctive star-formation histories. The central stellar mass density of this galaxy is within a factor of two of the most massive present-day ellipticals, while being globally 1000 times less massive. The radial profile of the specific star-formation rate is strongly rising toward the outskirts. This evidence strongly suggests the first detection of inside-out growth of a galaxy as a proto-bulge and a star-forming disc in the Epoch of Reionization.

A scaling relation based on thin-disc accretion theory has been used by some workers to determine the mass-inflow rate onto 20 high-redshift (z) and 80 Palomar-Green quasars. Based on several assumptions, it inexplicably implies that the inflow rate is an inverse function of black-hole (BH) mass Mbh. Moreover, its results remain untested. This paper offers a simple empirical relation essentially free of assumptions, found using available data for 59 highest-z quasars and the so-called Salpeter relation. We find that the accretion rate is proportional to Mbh(1+z)3, consistent with conventional astrophysics that the accretion rate is a direct function of both Mbh and the ambient gas density. We apply it to the 20 high-z and a subset of Palomar-Green quasars. Comparative analyses show that all empirically derived accretion rates and radiative efficiencies pass the tests, but their theoretical counterparts fail in most cases. A secondary relation defines the Eddington ratio as a function of z and radiative efficiency. Consistent with the empirical relations, spline regression analysis of Kozlowski's data for 132,000 quasars at z<2.4 shows that both the Eddington ratio and radiative efficiency are functions of Mbh and z. The results show that bigger BHs accrete more efficiently the smaller ones. For BHs > a billion solar masses, we get radiative efficiency of ~0.23 at z>5.7 and ~0.84 at z<0.005. Notably, the empirical relations predict a mass-inflow rate of 0.11-0,21 solar mass/year on to the BH in M87 that matches its Bondi accretion rate determined using observed density and temperature profiles.

The detailed modelling of stellar oscillations is a powerful approach to characterising stars. However, poor treatment of systematics in theoretical models leads to misinterpretations of stars. Here we propose a more principled statistical treatment for the systematics to be applied to fitting individual mode frequencies with a typical stellar model grid. We introduce a correlated noise model based on a Gaussian Process (GP) kernel to describe the systematics given that mode frequency systematics are expected to be highly correlated. We show that tuning the GP kernel can reproduce general features of frequency variations for changing model input physics and fundamental parameters. Fits with the correlated noise model better recover stellar parameters than traditional methods which either ignore the systematics or treat them as uncorrelated noise.

We have observed 123 pulsars with periods longer than 0.1 seconds in the Meterwavelength Single-pulse Polarimetric Emission Survey. In this work a detailed study of polarization behaviour of these pulsars have been carried out. We were able to fit the rotating vector model to the polarization position angle sweeps in 68 pulsars, and in 34 pulsars the emission heights could be measured. In all cases the radio emission was constrained to arise below 10\% of the light cylinder radius. In pulsars with low spindown energy loss, $\dot{E}<10^{34}$ ergs s$^{-1}$, we found the mean fractional linear polarization of the individual times samples in single pulses to be around 0.57 (57\%) which is significantly larger than the fractional linear polarization of 0.29 (29\%) obtained from the average profiles. On the other hand the mean fractional circular polarization of the individual time samples in single pulses is around 0.08 (8\%), similar to the measurements from the average profiles. To explain the observed polarization features, we invoke the partially screened vacuum gap model of pulsars, where dense spark associated plasma clouds exist with high pair plasma multiplicity, with significant decrease of density in the regions between the clouds, that are dominated by iron ions. The coherent radio emission is excited by curvature radiation from charge bunches in these dense plasma clouds and escape as linearly polarized waves near cloud boundaries. We suggest that the circular polarization arises due to propagation of waves in the low pair multiplicity, ion dominated inter-cloud regions.

We show how the star formation activity of galaxies is progressively inhibited from the outer region to the center of the massive cluster A2142. From an extended spectroscopic redshift survey of 2239 galaxies covering a circular area of radius $\sim 11$~Mpc from the cluster center, we extract a sample of 333 galaxies with known stellar mass, star formation rate, and spectral index $D_n4000$. We use the Blooming Tree algorithm to identify the substructures of the cluster and separate the galaxy sample into substructure galaxies, halo galaxies and outskirt galaxies. The substructure and halo galaxies are cluster members, whereas the outskirt galaxies are only weakly gravitationally bound to the cluster. For the cluster members, the star formation rate per stellar mass decreases with decreasing distance $R$ from the cluster center. Similarly, the spectral index $D_n4000$ increases with $R$, indicating an increasingly average age of the stellar population in galaxies closer to the cluster center. In addition, star formation in substructure galaxies is generally more active than in halo galaxies and less active than in outskirt galaxies, proving that substructures tend to slow down the transition between field galaxies and cluster galaxies. We finally show that most actively star forming galaxies are within the cluster infall region, whereas most galaxies in the central region are quiescent.

We use the eROSITA Final Equatorial-Depth Survey (eFEDS) to measure the rest-frame 0.1-2.4 keV band X-ray luminosities of $\sim$ 600,000 DESI groups using two different algorithms in the overlap region of the two observations. These groups span a large redshift range of $0.0 \le z_g \le 1.0$ and group mass range of $10^{10.76}h^{-1}M_{\odot} \le M_h \le 10^{15.0}h^{-1}M_{\odot}$. (1) Using the blind detection pipeline of eFEDS, we find that 10932 X-ray emission peaks can be cross matched with our groups, $\sim 38 \%$ of which have signal-to-noise ratio $\rm{S}/\rm{N} \geq 3$ in X-ray detection. Comparing to the numbers reported in previous studies, this matched sample size is a factor of $\sim 6$ larger. (2) By stacking X-ray maps around groups with similar masses and redshifts, we measure the average X-ray luminosity of groups as a function of halo mass in five redshift bins. We find, in a wide halo mass range, the X-ray luminosity, $L_{\rm X}$, is roughly linearly proportional to $M_{h}$, and is quite independent to the redshift of the groups. (3) We use a Poisson distribution to model the X-ray luminosities obtained using two different algorithms and obtain best-fit $L_{\rm X}=10^{28.46\pm0.03}M_{h}^{1.024\pm0.002}$ and $L_{\rm X}=10^{26.73 \pm 0.04}M_{h}^{1.140 \pm 0.003}$ scaling relations, respectively. The best-fit slopes are flatter than the results previously obtained, but closer to a self-similar prediction.

Planets in young star clusters could shed light on planet formation and evolution since star clusters can provide accurate age estimation. However, the number of transiting planets detected in clusters was only $\sim 30$, too small for statistical analysis. Thanks to the unprecedented high-precision astrometric data provided by Gaia DR2 and Gaia DR3, many new Open Clusters(OCs) and comoving groups have been identified. The UPiC project aims to find observational evidence and interpret how planet form and evolve in cluster environments. In this work, we cross-match the stellar catalogs of new OCs and comoving groups with confirmed planets and candidates. We carefully remove false positives and obtain the biggest catalog of planets in star clusters up to now, which consists of 73 confirmed planets and 84 planet candidates. After age validation, we obtain the radius--age diagram of these planets/candidates. We find an increment of the fraction of Hot Jupiters(HJs) around 100 Myr and attribute the increment to the flyby-induced high-e migration in star clusters. An additional small bump of the fraction of HJs after 1 Gyr is detected, which indicates the formation timescale of HJ around field stars is much larger than that in star clusters. Thus, stellar environments play important roles in the formation of HJs. The hot-Neptune desert occurs around 100 Myr in our sample. A combination of photoevaporation and high-e migration may sculpt the hot-Neptune desert in clusters.

Classification of galaxies is traditionally associated with their morphologies through visual inspection of images. The amount of data to come renders this task inhuman and Machine Learning (mainly Deep Learning) has been called to the rescue for more than a decade. However, the results look mitigate and there seems to be a shift away from the paradigm of the traditional morphological classification of galaxies. In this paper, I want to show that the algorithms indeed are very sensitive to the features present in images, features that do not necessarily correspond to the Hubble or de Vaucouleurs vision of a galaxy. However, this does not preclude to get the correct insights into the physics of galaxies. I have applied a state-of-the-art ''traditional'' Machine Learning clustering tool, called Fisher-EM, a latent discriminant subspace Gaussian Mixture Model algorithm, to 4458 galaxies carefully classified into 18 types by the EFIGI project. The optimum number of clusters given by the Integrated Complete Likelihood criterion is 47. The correspondence with the EFIGI classification is correct, but it appears that the Fisher-EM algorithm gives a great importance to the distribution of light which translates to characteristics such as the bulge to disk ratio, the inclination or the presence of foreground stars. The discrimination of some physical parameters (bulge-to-total luminosity ratio, $(B -- B)T$ , intrinsic diameter, presence of flocculence or dust, arm strength) is very comparable in the two classifications.

The ground-based gravitational wave (GW) observatories discover a population of merging stellar binary black holes (BBHs), which are promising targets for multiband observations by the low-, middle-, and high-frequency GW detectors. In this paper, we investigate the multiband GW detections of BBHs and demonstrate the advantages of such observations in improving the localization and parameter estimates of the sources. We generate mock samples of BBHs by considering different formation models as well as the merger rate density constrained by the current observations (GWTC-3). We specifically consider the astrodynamical middle-frequency interferometer GW observatory (AMIGO) in the middle-frequency band and estimate that it may detect $21$-$91$ BBHs with signal-to-noise ratio $\varrho\geq8$ in a $4$-yr observation period. The multiband observations by the low-frequency detectors [Laser Interferometer Space Antenna (LISA) and Taiji] and AMIGO may detect $5$-$33$ BBHs with $\varrho_{\rm LT}\geq5$ and $\varrho_{\rm AMI}\geq5$, which can evolve to the high-frequency band within $4$ yr and can be detected by the Cosmic Explorer (CE) and Einstein Telescope (ET). The joint observations of LISA-Taiji-AMIGO-ET-CE may localize the majority of the detectable BBHs in sky areas of $7\times10^{-7}$ to $2\times10^{-3}$ deg$^2$, which is improved by a factor of $\sim120$, $\sim2.4\times10^{5}$, $\sim1.8\times10^{4}$, or $\sim1.2\times10^{4}$, comparing with those by only adopting CE-ET, AMIGO, LISA-Taiji, or LISA-Taiji-AMIGO. These joint observations can also lead to an improvement of the measurement precision of the chirp mass (symmetric mass ratio) by a factor of $\sim5.5\times10^{4}$ ($33$), $\sim16$ ($8$), $\sim120$ ($90$), or $\sim5$ ($5$), comparing with those by CE-ET, AMIGO, LISA-Taiji, or LISA-Taiji-AMIGO.

We investigate the influence of an accretion disk on the angular momentum (AM) evolution of young M dwarfs, which parameters govern the AM distribution after the disk phase, and whether this leads to a mass-independent distribution of SAM. We find that above an initial rate $\Dot{M}_\mathrm{crit} \sim 10^{-8}~\Msolpyr$ accretion "erases" the initial SAM of M dwarfs during the disk lifetime, and stellar rotation converges to values of SAM that are largely independent of initial conditions. For stellar masses $> 0.3~\mathrm{M_\odot}$, we find that observed initial accretion rates $\Dot{M}_\mathrm{init}$ are comparable to or exceed $\Dot{M}_\mathrm{crit}$. Furthermore, stellar SAM after the disk phase scales with the stellar magnetic field strength as a power-law with an exponent of $-1.1$. For lower stellar masses, $\Dot{M}_\mathrm{init}$ is predicted to be smaller than $\Dot{M}_\mathrm{crit}$ and the initial conditions are imprinted in the stellar SAM after the disk phase. To explain the observed mass-independent distribution of SAM, the stellar magnetic field strength has to range between 20~G and 500~G (700~G and 1500~G) for a 0.1~$\mathrm{M_\odot}$ (0.6~$\mathrm{M_\odot}$) star. These values match observed large-scale magnetic field measurements of young M~dwarfs and the positive relation between stellar mass and magnetic field strength agrees with a theoretically-motivated scaling relation. The scaling law between stellar SAM, mass, and the magnetic field strength is consistent for young stars, where these parameters are constrained by observations. Due to the very limited number of available data, we advocate for efforts to obtain more such measurements. Our results provide new constraints on the relation between stellar mass and magnetic field strength and can be used as initial conditions for future stellar spin models, starting after the disk phase. (shortened)

We present a distant T$-$type brown dwarf candidate at $\approx2.55$ kpc discovered in the Cosmic Evolution Early Release Science (CEERS) fields by James Webb Space Telescope (JWST) NIRCam. In addition to the superb sensitivity, we utilised 7 filters from JWST in near-IR and thus is advantageous in finding faint, previously unseen brown dwarfs. From the model spectra in new JWST/NIRCam filter wavelengths, the selection criteria of F115W-F277W$<$-0.8 and F277W-F444W$>$1.1 were chosen to target the spectrum features of brown dwarfs having temperatures from 500K to 1300K. Searching through the data from Early Release Observations (ERO) and Early Release Science (ERS), we find 1 promising candidate in the CEERS field. The result of SED fitting suggested an early T spectral type with a low effective temperature of T$_\text{eff}\approx$1300K, the surface gravity of $\log{g}\approx5.25\text{cm s}^{-2}$, and an eddy diffusion parameter of logK$_{zz}\approx7\text{cm}^2 \text{s}^{-1}$, which indicates an age of $\approx$1.8Gyr and a mass of $\approx0.05$M$_{\odot}$. In contrast to typically found T$-$dwarf within several hundred parsecs, the estimated distance of the source is $\approx2.55$kpc, showing the JWST's power to extend the search to a much larger distance.

This study searches for neutrino signals from 18 dwarf spheroidal galaxies (dSphs) using 10 years of publicly available muon-track data of the IceCube neutrino observatory. We apply an unbinned likelihood analysis on each of these dSphs to derive the significance the putative neutrino emission. To further enhance our sensitivity, we also stack all dSphs together to perform a joint analysis. However, no significant neutrino emission signal was detected in either the single-source or stacking analysis. Based on these null results, we derive constraints on the annihilation cross section of dark matter particles. Compared to the existing literature, our constraints via the channel $\chi\chi\rightarrow\mu^+\mu^-$ are comparable to the ones from the VERITAS observations of dSphs.

We present an image classification algorithm using deep learning convolutional neural network architecture, which classifies the morphologies of eclipsing binary systems based on their light curves. The algorithm trains the machine with light curve images generated from the observational data of eclipsing binary stars in contact, detached and semi-detached morphologies, whose light curves are provided by Kepler, ASAS and CALEB catalogues. The structure of the architecture is explained, the parameters of the network layers and the resulting metrics are discussed. Our results show that the algorithm, which is selected among 132 neural network architectures, estimates the morphological classes of an independent validation dataset, 705 real data, with an accuracy of 92%.

The Mega electron volt (MeV) gamma-ray observation is a promising diagnostic tool for observing the universe. However, the sensitivity of MeV gamma-ray telescopes is limited due to peculiar backgrounds, restricting the application of MeV gamma rays for observation. Identification of backgrounds is crucial for designing next-generation telescopes. Therefore, herein, we assessed the background contribution in the electron-tracking Compton camera (ETCC) on board the SMILE- 2+ balloon experiment. This assessment was performed using the Monte Carlo simulation. The results revealed that the background below 400 keV existed due to the atmospheric gamma-ray background, the cosmic-ray/secondary-particle background, and the accidental background. On the other hand, the unresolved background component, which was not likely to be relevant to direct Compton-scattering events in the ETCC, was confirmed above 400 keV. Overall, this study demonstrated that the Compton-kinematics test provides a powerful tool to remove the background and principally improves the signal-to-noise ratio at 400 keV by an order of magnitude.

The High Energy Stereoscopic System (H.E.S.S.) is an array of five imaging atmospheric Cherenkov telescopes (IACTs) to study gamma-ray emission from astrophysical objects in the Southern hemisphere. It is the only hybrid array of IACTs, composed of telescopes with different collection areas and footprints, individually optimised for a specific energy range. Collectively, the array is most sensitive to gamma rays in the range of 100 GeV to 100 TeV. The array has been in operation since 2002 and has been upgraded with new telescopes and cameras multiple times. Recent hardware upgrades and changes in the operational procedures increased the amount of observing time, which is of key importance for time-domain science. H.E.S.S. operations saw record data taking in 2020 and 2021 and we describe the current operations with specific emphasis on system performance, operational processes and workflows, quality control, and (near) real-time extraction of science results. In light of this, we will briefly discuss the early detection of gamma-ray emission from the recurrent nova RS Oph and alert distribution to the astrophysics community.

We report on a comprehensive analysis of simultaneous X-ray polarimetric and spectral data of the bright atoll source GX 9+9 with the Imaging X-ray Polarimetry Explorer (IXPE) and NuSTAR. The source is significantly polarized in the 4--8 keV band, with a degree of $2.2\% \pm 0.5\%$ (uncertainty at the 68% confidence level). The NuSTAR broad-band spectrum clearly shows an iron line, and is well described by a model including thermal disk emission, a Comptonized component, and reflection. From a spectro-polarimetric fit, we obtain an upper limit to the polarization degree of the disk of 4% (at 99% confidence level), while the contribution of Comptonized and reflected radiation cannot be conclusively separated. However, the polarization is consistent with resulting from a combination of Comptonization in a boundary or spreading layer, plus reflection off the disc, which gives a significant contribution in any realistic scenario.

Nanoflare heating through small-scale magnetic reconnection events is one of the prime candidates to explain heating of the solar corona. However, direct signatures of nanoflares are difficult to determine, and unambiguous observational evidence is still lacking. Numerical models that include accelerated electrons, and can reproduce flaring conditions, are essential in understanding how low-energetic events act as a heating mechanism of the corona, and how such events are able to produce signatures in the spectral lines that can be detected through observations. We investigate the effects of accelerated electrons in synthetic spectra from a 3D radiative magnetohydrodynamics simulation to better understand small-scale heating events and their impact on the solar atmosphere. We synthesised the chromospheric Ca II and Mg II lines and the transition region Si IV resonance lines from a quiet Sun numerical simulation that includes accelerated electrons. We calculated the contribution function to the intensity to better understand how the lines are formed, and what factors are contributing to the detailed shape of the spectral profiles. The synthetic spectra are highly affected by variations in temperature and vertical velocity. Beam heating exceeds conductive heating at the heights where the spectral lines form, indicating that the electrons should contribute to the heating of the lower atmosphere and hence affect the line profiles. However, we find that it is difficult to determine specific signatures from the non-thermal electrons due to the complexity of the atmospheric response to the heating in combination with the relatively low energy output (~1e21 erg/s). Still, our results contribute to the understanding of small-scale heating events in the solar atmosphere, and give further guidance to future observations.

ED-2 is a stellar stream identified as a compact group in integrals of motion space in a local sample of halo stars from the third Gaia data release. Here we investigate its nature and possible association with known halo substructures. We explore the current properties of ED-2 members in phase-space, and also analyse the expected distribution via orbit integration. In addition, we study the metallicity of ED-2 using APOGEE DR17 and LAMOST DR8 (and re-calibrated DR3). ED-2 forms a compact group in the $x-z$ (or $R-z$) plane, showing a pancake-like structure as it crosses the Solar neighbourhood. Dynamically it is most similar the globular clusters NGC 3201 and NGC 6101, and the stellar stream Ylgr and Phlegethon. However, its orbit is sufficiently different that none of these objects is likely to be ED-2's progenitor. We also find ED-2 to be quite metal-poor, with all of its stars $\mathrm{[Fe/H]} \leq -2.42$, with a median $\mathrm{[Fe/H]} = -2.60^{+0.20}_{-0.21}$. At this low metallicity, it is unlikely that ED-2 stems from any known globular cluster, instead, ED-2 seems to be in a similar category as the recently discovered Phoenix and C-19 stellar streams. We find that ED-2 members are scattered across the whole sky, which is due to its current orbital phase. We predict that as this object moves to its next apocentre it will acquire an on-sky morphology that is akin to cold stellar streams. Finally, since ED-2 is nearing pericentre, we predict that additional members found below the plane should have large radial velocities, close to $\sim$ 500 km/s in the present-day direction of the globular cluster NGC 6101.

In this paper we study the production of Primordial Black Holes (PBHs) from inflation in order to explain the Dark Mater (DM) in the Universe. The evaluation of the fractional PBHs abundance to DM is sensitive to the value of the threshold $\mathrm{\delta_c}$ and the exact value of $\mathrm{\delta_c}$ is sensitive to the specific shape of the cosmological fluctuations. Different mechanisms producing PBHs lead to different thresholds and hence to different fractional abundances of PBHs. In this study, we examine various classes of inflationary models proposed in the existing literature to elucidate the formation of PBHs and we evaluate numerically the associated threshold values. Having evaluated the thresholds we compute the abundances of PBHs to DM using the Press Schecter approach and the Peak Theory. Given the influence of different power spectra on the thresholds, we investigate whether these inflationary models can successfully account for a significant fraction of DM. Moreover, we provide suggested values for the critical threshold. By examining the interplay between inflationary models, threshold values, and PBH abundances, our study aims to shed light on the viability of PBHs as a candidate for DM and contributes to the ongoing discussion regarding the nature of DM in the Universe

Oscillations are abundant in the solar corona. Coronal loop oscillations are typically studied using highly idealised models of magnetic flux tubes. In order to improve our understanding of coronal oscillations, it is necessary to consider the effect of realistic magnetic field topology and density structuring. We analyse the damping of coronal oscillations using a self-consistent 3D radiation-MHD simulation of the solar atmosphere spanning from the convection zone into the corona, the associated oscillation dissipation and heating, and finally the physical processes responsible for the damping and dissipation. The simulated corona formed in such a model does not depend on any prior assumptions about the shape of the coronal loops. We find that the bundle of magnetic loops shows damped transverse oscillations in response to perturbations in two separate instances with oscillation periods of 177 s and 191 s, velocity amplitudes of 10 km/s and 16 km/s and damping times of 176 s and 198 s, respectively. The coronal oscillations lead to the development of velocity shear in the simulated corona resulting in the formation of vortices seen in the velocity field caused by the Kelvin-Helmholtz instability, contributing to the damping and dissipation of the transverse oscillations. The oscillation parameters and evolution observed are in line with the values typically seen in observations of coronal loop oscillations. The dynamic evolution of the coronal loop bundle suggests the models of monolithic and static coronal loops with constant lengths might need to be re-evaluated by relaxing the assumption of highly idealised waveguides.

We present the software package binary_c-python which provides a convenient and easy-to-use interface to the binary_c framework, allowing the user to rapidly evolve individual systems and populations of stars. binary_c-python is available on Pip and on GitLab. binary_c-python contains many useful features to control and process the output of binary_c, like by providing binary_c-python with logging statements that are dynamically compiled and loaded into binary_c. Moreover, we have recently added standardised output of events like Roche-lobe overflow or double compact-object formation to binary_c, and automatic parsing and managing of that output in binary_c-python. binary_c-python uses multiprocessing to utilise all the cores on a particular machine, and can run populations with HPC cluster workload managers like HTCondor and Slurm, allowing the user to run simulations on large computing clusters. We provide documentation that is automatically generated based on docstrings and a suite of Jupyter notebooks. These notebooks consist of technical tutorials on how to use binary_c-python and use-case scenarios aimed at doing science. Much of binary_c-python is covered by unit tests to ensure reliability and correctness, and the test coverage is continually increased as the package is improved.

Aims. The main goal of this work is to date young open clusters using $\delta$ Sct stars. Seismic indices such as the large separation and the frequency at maximum power can help to constrain the models to better characterise the stars. We propose a reliable method to identify some radial modes, which gives us greater confidence in the constrained models. Methods. We extract the frequency content of a sample of $\delta$ Sct stars belonging to the same open cluster. We estimate the low-order large separation by means of different techniques and the frequency at maximum power for each member of the sample. We use a grid of models built with the typical parameters of $\delta$ Sct stars, including mass, metallicity and rotation as independent variables, and determine the oscillation modes. We select the observed frequencies whose ratios match those of the models. Once we find a range of radial modes matching the observed frequencies, mainly the fundamental mode, we add it to the other seismic parameters to derive the stellar age. Assuming star groups have similar chemistry and age, we estimate their mean age by computing a weighted probability density function fit to the age distribution of the seismically constrained models. Results. We estimate the age of Trumpler 10 to be $30_{-20}{+30}$ Myr, and that of Praesepe to be $580 \pm 230$ Myr. In this latter case, we find two apparent populations of $\delta$ Sct stars in the same cluster, one at $510 \pm 140$ Myr and another at $890 \pm 140$ Myr. This may be due to two different formation events, different rotational velocities of the members in our sample of stars (as rapid rotation may modify the observed large separation), or to membership of unresolved binary systems.

In this work, we consider polar perturbations and we calculate the frequency and damping time of the quadrupolar fundamental f -mode of compact objects, constructed using a wide range of model-independent hybrid equations of state that include quark matter. We give special attention to the impact of the hadron-quark conversion speed that, in the slow case, gives rise to a branch of slow stable hybrid stars. Moreover, we study the validity of universal relationships proposed in the literature and find out that none of them remains valid when slow stable hybrid stars are taken into account. This fact could constrain the applicability of asteroseismology methods with fundamental modes designed to estimate the properties of pulsating compact objects. We hope that this result could be tested with the start up of the third-generation gravitational wave observatories, which might shed some light on the f -mode emission from compact objects.

One of the commonly used non-parametric morphometric statistics for galaxy profiles and images is the asymmetry statistic. With an eye to current and upcoming large neutral hydrogen (HI) surveys, we develop a 3D version of the asymmetry statistic that can be applied to datacubes. This statistic is more resilient to variations due to the observed geometry than 1D asymmetry measures, and can be successfully applied to lower spatial resolutions (3-4 beams across the galaxy major axis) than the 2D statistic. We have also modified the asymmetry definition from an `absolute difference' version to a `squared difference' version that removes much of the bias due to noise contributions for low signal-to-noise observations. Using a suite of mock asymmetric cubes we show that the background-corrected, squared difference 3D asymmetry statistic can be applied to many marginally resolved galaxies in large wide-area HI surveys such as WALLABY on the Australian SKA Pathfinder (ASKAP).

Early dust grain growth in protostellar envelopes infalling on young discs has been suggested in recent studies, supporting the hypothesis that dust particles start to agglomerate already during the Class 0/I phase of young stellar objects (YSOs). If this early evolution were confirmed, it would impact the usually assumed initial conditions of planet formation, where only particles with sizes $\lesssim 0.25 \mu$m are usually considered for protostellar envelopes. We aim to determine the maximum grain size of the dust population in the envelope of the Class 0/I protostar L1527 IRS, located in the Taurus star-forming region (140 pc). We use Atacama Large millimetre/sub-millimetre Array (ALMA) and Atacama Compact Array (ACA) archival data and present new observations, in an effort to both enhance the signal-to-noise ratio of the faint extended continuum emission and properly account for the compact emission from the inner disc. Using observations performed in four wavelength bands and extending the spatial range of previous studies, we aim to place tight constraints on the spectral ($\alpha$) and dust emissivity ($\beta$) indices in the envelope of L1527 IRS. We find a rather flat $\alpha \sim$ 3.0 profile in the range 50-2000 au. Accounting for the envelope temperature profile, we derive values for the dust emissivity index, 0.9 < $\beta$ < 1.6, and reveal a tentative, positive outward gradient. This could be interpreted as a distribution of mainly ISM-like grains at 2000 au, gradually progressing to (sub-)millimetre-sized dust grains in the inner envelope, where at R=300 au, $\beta$ = 1.1 +/- 0.1. Our study supports a variation of the dust properties in the envelope of L1527 IRS. We discuss how this can be the result of in-situ grain growth, dust differential collapse from the parent core, or upward transport of disc large grains.

Rotating Radio Transients (RRATs) are neutron stars that emit sporadic radio bursts. We detected 1955 single pulses from RRAT J1913+1330 using the 19-beam receiver of the Five-hundred-meter Aperture Spherical Radio Telescope (FAST). These pulses were detected in 19 distinct clusters, with 49.4% of them occurring sequentially with a waiting time of one rotation period. The energy distribution of these individual pulses exhibited a wide range, spanning three orders of magnitude, reminiscent of repeating fast radio bursts (FRBs). Furthermore, we observed abrupt variations in pulse profile, width, peak flux, and fluence between adjacent sequential pulses. These findings suggest that this RRAT could be interpreted as pulsars with extreme pulse-to-pulse modulation. The presence of sequential pulse trains during active phases, along with significant pulse variations in profile, fluence, flux, and width, should be intrinsic to a subset of RRATs. Our results indicate that J1913+1330 represents a peculiar source that shares certain properties with populations of nulling pulsars, giant pulses, and FRBs from different perspectives. The dramatic pulse-to-pulse variation observed in J1913+1330 could be attributed to unstable pair creation above the polar cap region and the variation of the site where streaming pairs emit coherently. Exploring a larger sample of RRATs exhibiting similar properties to J1913+1330 has the potential to significantly advance our understanding of pulsars, RRATs, and FRBs.

Supernova remnants (SNRs) may regulate star formation in galaxies. For example, SNR-driven shocks may form new molecular gas or compress pre-existing clouds and trigger the formation of new stars. To test this scenario, we measure the deuteration of $N_2H^+$, $D_{frac}^{N_2H^+}$, a well-studied tracer of pre-stellar cores, across the Infrared Dark Cloud (IRDC) G034.77-00.55, known to be experiencing a shock interaction with the SNR W44. We use N$_2$H$^+$ and N$_2$D$^+$ J=1-0 single pointing observations obtained with the 30m antenna at the Instituto de Radioastronomia Millimetrica to infer $D_{frac}^{N_2H^+}$ toward five positions across the cloud, namely a massive core, different regions across the shock front, a dense clump and ambient gas. We find $D_{frac}^{N_2H^+}$ in the range 0.03-0.1, several orders of magnitude larger than the cosmic D/H ratio ($\sim$10$^{-5}$). Across the shock front, $D_{frac}^{N_2H^+}$ is enhanced by more than a factor of 2 ($D_{frac}^{N_2H^+}\sim$0.05-0.07) with respect to the ambient gas ($\leq$0.03) and similar to that measured generally in pre-stellar cores. Indeed, in the massive core and dense clump regions of this IRDC we measure $D_{frac}^{N_2H^+}$}$\sim$0.1. We find enhanced deuteration of $N_2H^+$ across the region of the shock, at a level that is enhanced with respect to regions of unperturbed gas. It is possible that this has been induced by shock compression, which would then be indirect evidence that the shock is triggering conditions for future star formation. However, since unperturbed dense regions also show elevated levels of deuteration, further, higher-resolution studies are needed to better understand the structure and kinematics of the deuterated material in the shock region, e.g., if it still in relatively diffuse form or already organised in a population of low-mass pre-stellar cores.

This thesis is dedicated to the study of stochastic processes; non-deterministic physical phenomena that can be well described by classical physics. The stochastic processes we are interested in are akin to Brownian Motion and can be described by an overdamped Langevin equation comprised of a deterministic drift term and a random noise term. In Part I we examine stochastic processes in the Mesoscale. For us this means that the Langevin equation is driven by thermal noise, with amplitude proportional to the temperature $T$ and there exists a genuine equilibrium thermal state. We apply a technique known as the Functional Renormalisation Group (FRG) which allows us to coarse-grain in temporal scales. We describe how to obtain effective equations of motion for the 1- and 2-point functions of a particle evolving in highly non-trivial potentials and verify their accuracy by comparison to direct numerical simulations. In this way we outline a novel procedure for describing the behaviour of stochastic processes without having to resort to time consuming numerical simulations. In Part II we turn to the Early Universe and in particular examine stochastic processes occurring during a period of accelerated expansion known as inflation. This inflationary period is driven by a scalar field called the inflaton which also obeys a Langevin equation in the Stochastic Inflation formalism. We use this to study the formation of Primordial Black Holes during a period of Ultra Slow-Roll. We finish this thesis by applying the techniques developed in Part I to a spectator field during inflation. FRG techniques can compute cosmologically relevant observables such as the power spectrum and spectral tilt. We also extend the FRG formalism to solve first-passage time problems and verify that it gives the correct prediction for the average time taken for a field (or particle) to overcome a barrier in the potential.

The search for biosignatures on exoplanets connects the fields of biology and biochemistry to astronomical observation, with the hope that we might detect evidence of active biological processes on worlds outside the solar system. Here we focus on a complementary aspect of exoplanet characterisation connecting astronomy to prebiotic chemistry: the search for molecules associated with the origin of life, prebiosignatures. Prebiosignature surveys in planetary atmospheres offer the potential to both constrain the ubiquity of life in the galaxy and provide important tests of current prebiotic syntheses outside of the laboratory setting. Here, we quantify the minimum abundance of identified prebiosignature molecules that would be required for detection by transmission spectroscopy using JWST. We consider prebiosignatures on five classes of terrestrial planet: an ocean planet, a volcanic planet, a post-impact planet, a super-Earth, and an early Earth analogue. Using a novel modelling and detection test pipeline, with simulated JWST noise, we find the detection thresholds of hydrogen cyanide (HCN), hydrogen sulfide (H2S), cyanoacetylene (HC3N), ammonia (NH3), methane (CH4), acetylene (C2H2), sulfur dioxide (SO2), nitric oxide (NO), formaldehyde (CH2O), and carbon monoxide (CO) in a variety of low mean molecular weight (<5) atmospheres. We test the dependence of these detection thresholds on M dwarf target star and the number of observed transits, finding that a modest number of transits (1-10) are required to detect prebiosignatures in numerous candidate planets, including TRAPPIST-1e with a high mean molecular weight atmosphere. We find that the NIRSpec G395M/H instrument is best suited for detecting most prebiosignatures.

The short-period AP Dor eclipsing binary's first in-depth and multiband photometric solutions are presented. We made use of our eight nights of ground-based at a southern hemisphere observatory, and twelve sectors of TESS observations. We extracted eight and 1322 minima from our observations and TESS, respectively. We suggested a new linear ephemeris based on the trend of orbital period variations using the Markov chain Monte Carlo (MCMC) approach. The PHysics Of Eclipsing BinariEs (PHOEBE) Python code and the MCMC approach were used for the light curve analysis. This system did not require a starspot for the light curve solutions. We calculated the absolute parameters of the system using Gaia DR3 parallax method. The orbital angular momentum (J_0) of the AP Dor indicates that this system is located in a region of contact binaries. According to our results, this system is an overcontact binary system with a mass ratio of 0.584, a fillout factor of 48%, and an inclination of 53deg. The positions of AP Dor stars on the Hertzsprung-Russell (HR) diagram are represented.

With JWST we can now characterize the atmospheres of planets on longer orbital planets, but this moves us into a regime where we cannot assume that tidal forces from the star have eroded planets' obliquities and synchronized their rotation rates. These rotation vectors may be tracers of formation and evolution histories and also enable a range of atmospheric circulation states. Here we delineate the orbital space over which tidal synchronization and alignment assumptions may no longer apply and present three-dimensional atmospheric models of a hypothetical warm Jupiter over a range of rotation rates and obliquities. We simulate the secondary eclipses of this planet for different possible viewing orientations and times during its orbital, seasonal cycle. We find that the eclipse depth can be strongly influenced by rotation rate and obliquity through the timing of the eclipse relative to the planet's seasonal cycle, and advise caution in attempting to derive properties such as albedo or day-night transport from this measurement. We predict that if warm Jupiters beyond the tidal limit have intrinsic diversity in their rotation vectors, then it will manifest itself as dispersion in their secondary eclipse depths. We explore eclipse mapping as a way to uniquely constrain the rotation vector of warm Jupiters but find that the associated signals are likely at the edge of JWST performance. Nevertheless, as JWST begins to measure the secondary eclipses of longer orbital period planets, we should expect to observe the consequences of a wider range of rotation states and circulation patterns.

New instrumental and telescopes covering the optical and ultra-violet spectral regions have revealed a range of small-scale dynamic features, many which may be related. For example, the range of spicule-like features hints towards a spectrum of features and not just two types; however, direct observational evidence in terms of tracking spicules across multiple wavelengths are needed in order to provide further insight into the dynamics of the Sun's outer atmosphere. This paper uses H $\alpha$ data obtained with the CRisp Imaging SpectroPolarimeter instrument on the Swedish 1-m Solar Telescope, and in the transition region using the Interface Region Imaging Spectrograph with the SJI 1400 {\AA} channel plus spectral data via the Si IV 1394 {\AA} line to track spicules termed Rapid Blue-shifted Excursions (RBEs). The RBEs as seen in the H $\alpha$ blue-wing images presented here can be sub-divided into two categories; a single or multi-threaded feature. Based on the H $\alpha$ spectra, the features can be divided into events showing broadening and line core absorption, events showing broadening and line core emission, events with a pure blue shifted H $\alpha$ profile without any absorption in the red wing, broadened line profile with the absorption in the blue stronger compared to the red wing. From the RBE-like events which have a Si IV 1394 {\AA} line profile, 78% of them show a Si IV line flux increase. Most of these features show a second broadened Si IV component which is slightly blue-shifted.

The TESS follow-up of a large number of known transiting exoplanets provide unique opportunity to study their physical properties more precisely. Being a space-based telescope, the TESS observations are devoid of any noise component resulting from the interference of Earth's atmosphere. TESS also provides a better probability to observe subsequent transit events owing to its longer uninterrupted time-series observations compared to the ground-based telescopes. Especially, for the exoplanets around bright host-stars, TESS time-series observations provides high SNR lightcurves, which can be used for higher precision studies for these exoplanets. In this work, I have studied the TESS transit photometric follow-up observations of 28 exoplanets around bright stars with $V_{mag}\le$10. The already high SNR lightcurves from TESS have been further processed with a critical noise treatment algorithm, using the wavelet denoising and the Gaussian-process regression techniques, to effectively reduce the noise components both correlated and uncorrelated in time, which were then used to estimate the physical properties of these exoplanets. The study has resulted in very precise values for the physical properties of the target exoplanets, with the improvements in precision being significant for most of the cases compared to the previous studies. Also, since a comparatively large number of transit lightcurves from TESS observations were used to estimate these physical properties for each of the target exoplanets, which removes any bias due to the lack of sufficient datasets, these updated physical properties can be considered extremely accurate and reliable for future studies.

As a potential consequence of Lorentz invariance violation~(LIV), threshold anomalies open a window to study LIV. Recently the Large High Altitude Air Shower~(LHAASO) observatory reported that more than 5000 photons from GRB 221009A have been observed with energies above 500~GeV and up to $18~\text{TeV}$. In the literature, it is suggested that this observation may have tension with the standard model result because extragalactic background light~(EBL) can prevent photons around 18~TeV photons from reaching the earth and that LIV induced threshold anomalies might be able to explain the observation. In this work we further study this proposal with more detailed numerical calculation for different LIV scales and redshifts of the sources. We find that GRB 221009A is a rather unique opportunity to search LIV, and a LIV scale $E_\text{LIV} \lesssim E_\text{Planck}\approx 1.22\times 10^{19}~\text{GeV}$ is feasible to the observation of GRB 221009A on 9 October, 2022.

We investigate the connection between the full- and flat-sky angular power spectra. First, we revisit this connection established on the geometric and physical grounds, namely that the angular correlations on the sphere and in the plane (flat-sky approximation) correspond to each other in the limiting case of small angles and a distant observer. To establish the formal conditions for this limit, we first resort to a simplified shape of the 3D power spectrum, which allows us to obtain analytic results for both the full- and flat-sky angular power spectra. Using a saddle-point approximation, we find that the flat-sky results are obtained in the limit when the comoving distance and wave modes $\ell$ approach infinity at the same rate. This allows us to obtain an analogous asymptotic expansion of the full-sky angular power spectrum for general 3D power spectrum shapes, including the LCDM Universe. In this way, we find a robust limit of correspondence between the full- and flat-sky results. These results also establish a mathematical relation, i.e., an asymptotic expansion of the ordinary hypergeometric function of a particular choice of arguments that physically corresponds to the flat-sky approximation of a distant observer. This asymptotic form of the ordinary hypergeometric function is obtained in two ways: relying on our saddle-point approximation and using some of the known properties of the hypergeometric function.

KIDSpec, the Kinetic Inductance Detector Spectrometer, is a proposed optical to near IR Microwave Kinetic Inductance Detector (MKID) spectrograph. MKIDs are superconducting photon counting detectors which are able to resolve the energy of incoming photons and their time of arrival. KIDSpec will use these detectors to separate incoming spectral orders from a grating, thereby not requiring a cross-disperser. In this paper we present a simulation tool for KIDSpec's potential performance upon construction to optimise a given design. This simulation tool is the KIDSpec Simulator (KSIM), a Python package designed to simulate a variety of KIDSpec and observation parameters. A range of astrophysical objects are simulated: stellar objects, an SDSS observed galaxy, a Seyfert galaxy, and a mock galaxy spectrum from the JAGUAR catalogue. Multiple medium spectral resolution designs for KIDSpec are simulated. The possible impact of MKID energy resolution variance and dead pixels were simulated, with impacts to KIDSpec performance observed using the Reduced Chi-Squared (RCS) value. Using dead pixel percentages from current instruments, the RCS result was found to only increase to 1.21 at worst for one of the designs simulated. SNR comparisons of object simulations between KSIM and X-Shooter's ETC were also simulated. KIDSpec offers a particular improvement over X-Shooter for short and faint observations. For a Seyfert galaxy ($m_{R}=21$) simulation with a 180s exposure, KIDSpec had an average SNR of 4.8, in contrast to 1.5 for X-Shooter. Using KSIM the design of KIDSpec can be optimised to improve the instrument further.

Determining the habitability and interpreting future atmospheric observations of exoplanets requires understanding the atmospheric dynamics and chemistry from a 3-D perspective. Previous studies have shown significant spatial variability in the ozone layer of synchronously rotating M-dwarf planets, assuming an Earth-like initial atmospheric composition. We use a 3-D Coupled Climate-Chemistry model to understand this distribution of ozone and identify the mechanism responsible for it. We document a previously unreported connection between the ozone production regions on the photochemically active dayside hemisphere and the nightside devoid of stellar radiation and thus photochemistry. We find that stratospheric dayside-to-nightside overturning circulation can advect ozone-rich air to the nightside. On the nightside, ozone-rich air subsides at the locations of two quasi-stationary Rossby gyres, resulting in an exchange between the stratosphere and troposphere and the accumulation of ozone at the gyre locations. We identify the hemispheric contrast in radiative heating and cooling as the main driver of this ozone circulation. Dynamically-driven chemistry also impacts other tracer species in the atmosphere (gaseous and non-gaseous phase) as long as chemical lifetimes exceed dynamical lifetimes. These findings illustrate the 3-D nature of planetary atmospheres, predicting spatial and temporal variability that will impact spectroscopic observations of exoplanet atmospheres.

Within the $f(Q)$-gravity framework we perform a phenomenological study of the cosmological observables in light of the degeneracy between neutrinos physics and the modified gravity parameter and we identify specific patterns which allow to break such degeneracy. We also provide separately constraints on the total mass of the neutrinos, $\Sigma m_{\nu}$, and on the effective number of neutrino species, $N_{\rm eff}$, using cosmic microwave background (CMB), baryon acoustic oscillation (BAO), redshift space distortion (RSD), supernovae (SNIa), galaxy clustering (GC) and weak gravitational lensing (WL) measurements. The strongest upper bound on the total mass of the neutrinos is found for the combination of CMB+BAO+RSD+SNIa and it is $\Sigma m_\nu<0.277$ eV at 95\% C.L. For the same combination of data we find $N_{\rm eff}=2.93^{+0.31}_{-0.34}$ at 95\% C.L. We also find that all combinations of data we consider, prefer a stronger gravitational interaction than $\Lambda$CDM. Finally, we consider the $\chi^2$ and deviance information criterion statistics and find the $f(Q)+\Sigma m_\nu$ model to be statistically supported by data over the standard scenario. On the contrary $f(Q)+N_{\rm eff}$ is supported by CMB+BAO+RSD+SNIa but a moderate evidence against it is found with GC and WL data.

The rapid development of gravitational wave astronomy provides the unique opportunity of exploring the dynamics of the Universe using clustering properties of coalescing binary black hole mergers. Gravitational wave data, along with information coming from future galaxy surveys, have the potential of shedding light about many open questions in Cosmology, including those regarding the nature of dark matter and dark energy. In this work we explore which combination of gravitational wave and galaxy survey datasets are able to provide the best constraints both on modified gravity theories and on the nature of the very same binary black hole events. In particular, by using the public Boltzmann code \texttt{Multi\_CLASS}, we compare cosmological constraints on popular $\Lambda$CDM extensions coming from gravitational waves alone and in conjunction with either deep and localized or wide and shallow galaxy surveys. We show that constraints on extensions of General Relativity will be at the same level of existing limits from gravitational waves alone or one order of magnitude better when galaxy surveys are included. Furthermore, cross-correlating both kind of galaxy survey with gravitational waves datasets will allow to confidently rule in or out primordial black holes as dark matter candidate in the majority of the allowed parameter space.

We present a detailed model atmosphere analysis of 14001 DA white dwarfs from the Montreal White Dwarf Database with ultraviolet photometry from the GALEX mission. We use the 100 pc sample, where the extinction is negligible, to demonstrate that there are no major systematic differences between the best-fit parameters derived from optical only data and the optical + UV photometry. GALEX FUV and NUV data improve the statistical errors in the model fits, especially for the hotter white dwarfs with spectral energy distributions that peak in the UV. Fitting the UV to optical spectral energy distributions also reveals UV-excess or UV-deficit objects. We use two different methods to identify outliers in our model fits. Known outliers include objects with unusual atmospheric compositions, strongly magnetic white dwarfs, and binary white dwarfs, including double degenerates and white dwarf + main-sequence systems. We present a list of 89 newly identified outliers based on GALEX UV data; follow-up observations of these objects will be required to constrain their nature. Several current and upcoming large scale spectroscopic surveys are targeting $>10^5$ white dwarfs. In addition, the ULTRASAT mission is planning an all-sky survey in the NUV band. A combination of the UV data from GALEX and ULTRASAT and optical data on these large samples of spectroscopically confirmed DA white dwarfs will provide an excellent opportunity to identify unusual white dwarfs in the solar neighborhood.

The first asteroid simulants workshop was held in late 2015. These materials are needed for tests of technologies and mission operational concepts, for training astronauts , for medical studies, and a variety of other purposes. The new program is based on lessons learned from the earlier lunar simulants program. It aims to deliver families of simulants for major spectral classes of asteroids both in cobble and regolith form, beginning with one type of carbonaceous chondrite and rapidly expanding to provide four to six more asteroid classes. These simulants will replicate a selected list of asteroid properties, but not all known properties, in order to serve the greatest number of users at an affordable price. They will be benchmarked by a variety of data sets including laboratory analysis of meteorites, observation of bolides, remote sensing of asteroids, data from asteroid missions, and scientific modeling. A variety of laboratory tests will verify the as-manufactured simulants are accurately and repeatedly providing the specified characteristics.

The Gaia-Sausage/Enceladus (GS/E) structure is an accretion remnant in the Milky Way's halo that constitutes a large fraction of the nearby stellar halo. We study GS/E using high-purity samples of kinematically selected stars from APOGEE DR16 and Gaia. Employing a novel modelling framework to account for kinematic selection biases using distribution functions, we fit density profiles to these GS/E samples and measure their masses. We find that GS/E is described by a shallow density profile in the inner Galaxy, with a break between 15-25 kpc beyond which the profile becomes very steep in the outer Galaxy. We also find that GS/E is triaxial, with axis ratios 1:0.55:0.45 (nearly prolate), and the major axis is oriented about 80~degrees from the Sun-Galactic center line and 16 degrees above the plane. We measure a stellar mass for GS/E of $1.45 ^{+0.92}_{-0.51}\,\mathrm{(stat.)}\,^{+0.08}_{-0.49} \mathrm{(sys.)}\ \times10^{8}$ M$_{\odot}$, which is lower than previously estimated. We also fit a density profile to the entire Milky Way stellar halo, finding a mass in the range of $6-8 \times 10^{8}$ M$_{\odot}$, and implying that GS/E could make up as little as 10-20 per cent of the mass of the Milky Way stellar halo. Our findings challenge recent works, which have found greater masses for GS/E, and often find that it dominates the mass budget of the stellar halo. Our lower stellar mass combined with standard stellar-mass-to-halo-mass relations implies that GS/E constituted a minor 1:8-mass-ratio merger at the time of its accretion.

The possible existence of primordial black holes (PBHs) is an open question in modern cosmology. Among the probes to test it, gravitational waves (GW) coming from their mergers constitute a powerful tool. In this work, we study how stellar mass PBH binaries could affect measurements of the clustering of merger events in future GW surveys. We account for PBH binaries formed both in the early and late Universe and show that the power spectrum modification they introduce can be detected at $\sim 2\sigma-3\sigma$ (depending on some assumptions) whenever PBH mergers make up at least $\sim 60\%$ of the overall number of detected events. By adding cross-correlations with galaxy surveys, this threshold is lowered to $\sim 40\%$. In the case of a poor redshift determination of GW sources, constraints are degraded by about a factor of 2. Assuming a theoretical model for the PBH merger rate, we can convert our results to constraints on the fraction of dark matter in PBHs, $f_{\rm PBH}$. Finally, we perform a Bayesian model selection forecast and confirm that the analysis we develop could be able to detect $\sim30M_\odot$ PBHs if they account for $f_{\rm PBH}\sim 10^{-4}-10^{-3}$, depending on the model uncertainty considered, being thus competitive with other probes.

Classical data analysis requires computational efforts that become intractable in the age of Big Data. An essential task in time series analysis is the extraction of physically meaningful information from a noisy time series. One algorithm devised for this very purpose is singular spectrum decomposition (SSD), an adaptive method that allows for the extraction of narrow-banded components from non-stationary and non-linear time series. The main computational bottleneck of this algorithm is the singular value decomposition (SVD). Quantum computing could facilitate a speedup in this domain through superior scaling laws. We propose quantum SSD by assigning the SVD subroutine to a quantum computer. The viability for implementation and performance of this hybrid algorithm on a near term hybrid quantum computer is investigated. In this work we show that by employing randomised SVD, we can impose a qubit limit on one of the circuits to improve scalibility. Using this, we efficiently perform quantum SSD on simulations of local field potentials recorded in brain tissue, as well as GW150914, the first detected gravitational wave event.

LHAASO is an instrument designed for detecting cosmic rays (CRs) and gamma rays at TeV to PeV energies. The decays of heavy dark matter particles in the Galactic halo may produce high-energy electrons that can be detected by LHAASO. The main background for the LHAASO's CR electron measurements is the hadron residuals due to mis-identification of the particle species. In this paper, we estimate the LHAASO's electron background using the known all-particle CR spectrum and the hadron rejection efficiency of LHAASO. With the estimated background, we predict the capability of LHAASO to constrain DM decay lifetime at 95% confidence level for various channels. We find that, if neglecting systematic uncertainties, the CR electron measurement by LHAASO can improve the current best results by up to on order of magnitude for DM masses between 100 - 1000TeV. However, indirect measurements of CR electrons by ground-based experiments suffer from uncertainties included in the calculation, the projected constraints will be largely weakened. So for using the CR electron observation of LHAASO to constrain the DM parameters, the key point is whether the systematic error can be effectively reduced.

We propose a method for including the effects of non-Gaussian velocity field distribution in the estimation of growth rate of the magnetic energy in a random flow with finite memory time. The method allows a reduction to the Gaussian case that was investigates earlier. For illustration we consider the multivariate Laplace distribution and compare it against the Gaussian one.

An emerging paradigm is being embraced in the conceptualization of future planetary exploration missions. Ambitious objectives and increasingly demanding mission constraints stress the importance associated with faster surface mobility. Driving speeds approaching or surpassing 1 m/s have been rarely used and their effect on performance is today unclear. This study presents experimental evidence and preliminary observations on the impact that increasing velocity has on the tractive performance of planetary rovers. Single-wheel driving tests were conducted using two different metallic, grousered wheels-one rigid and one flexible-over two different soils, olivine sand and CaCO3-based silty soil. Experiments were conducted at speeds between 0.01-1 m/s throughout an ample range of slip ratios (5-90%). Three performance metrics were evaluated: drawbar pull coefficient, wheel sinkage, and tractive efficiency. Results showed similar data trends among all the cases investigated. Drawbar pull and tractive efficiency considerably decreased for speeds beyond 0.2 m/s. Wheel sinkage, unlike what published evidence suggested, increased with increasing velocities. The flexible wheel performed the best at 1m/s, exhibiting 2 times higher drawbar pull and efficiency with 18% lower sinkage under low slip conditions. Although similar data trends were obtained, a different wheel-soil interactive behavior was observed when driving over the different soils. Overall, despite the performance reduction experienced at higher velocities, a speed in the range of 0.2-0.3 m/s would enable 5-10 times faster traverses, compared to current rovers driving capability, while only diminishing drawbar pull and efficiency by 7%. The measurements collected and the analysis presented here lay the groundwork for initial stages in the development of new locomotion subsystems for planetary surface exploration. At the same time...

We propose a novel experimental method for probing light dark matter candidates. We show that an electro-optical material's refractive index is modified in the presence of a coherently oscillating dark matter background. A high-precision resonant Michelson interferometer can be used to read out this signal. The proposed detection scheme allows for the exploration of an uncharted parameter space of dark matter candidates over a wide range of masses -- including masses exceeding a few tens of microelectronvolts, which is a challenging parameter space for microwave cavity haloscopes.

We investigate stationary accretion of the collisionless Vlasov gas onto the Kerr black hole, occurring in the equatorial plane. At infinity the gas obeys the Maxwell-J\"{u}ttner distribution, restricted to the equatorial plane. In the vicinity of the black hole, the motion of the gas is governed by the spacetime geometry. We compute accretion rates of the rest-mass, the energy, and the angular momentum, as well as the particle number surface density, focusing on the dependence of these quantities on the asymptotic temperature of the gas and the black hole spin. The accretion slows down the rotation of the black hole. We present preliminary results for a Vlasov gas accretion onto a Kerr black hole moving with a velocity parallel to the equatorial plane.

Recent indirect searches of dark matter using gamma-ray, radio, and cosmic-ray data have provided some stringent constraints on annihilating dark matter. In this article, we propose a new indirect method to constrain annihilating dark matter. By using the data of the G2 cloud near the Galactic supermassive black hole Sgr A*, we can get stringent constraints on the parameter space of dark matter mass and the annihilation cross section, especially for the non-leptophilic annihilation channels $b\bar{b}$ and $W^{\pm}$. For the thermal annihilation cross section, the lower bounds of dark matter mass can be constrained up to TeV order for the non-leptophilic channels with the standard spike index $\gamma_{\rm sp}=7/3$.

Momentum space distributions of photons coming out of any light emitting materials/devices provide critical information about its underlying physical origin. Conventional methods of determining such properties impose specific instrumentational difficulties for probing samples kept within a low temperature cryostat. There were past studies to measure one dimensional (1D) coherence function which could then be used for extracting momentum space information as well as reports of measurements of just two dimensional (2D) coherence function. However, all of those are associated with additional experimental complexities. So, here we propose a simpler, modified Michelson interferometer based optical setup kept at room temperature outside the cryostat to initially measure the 2D coherence function of emitted light, which can then be used to directly estimate the 2D in-plane momentum space distribution by calculating its fast Fourier transform. We will also discuss how this experimental method can overcome instrumentational difficulties encountered in past studies.

We present numerical solutions to Einstein's equations describing large spherical cosmic voids constituted by two components; dark matter and baryons, with a non-vanishing initial relative velocity, in an asymptotically homogeneous background compatible with the $\Lambda$CDM concordance model. We compute numerically the evolution of such configurations in the dark matter frame, with a hypothetical homogeneous distribution of baryons, but respecting the values dictated by the concordance model for the average baryon-to-dark matter density ratio. We reproduce the well known formation of overdensities at the edge of the void, and recover the Lemaitre-Tolman-Bondi solutions in the comoving limit of our simulations. We compute the average growth factor of matter fluctuations, and find that it departs significantly from the linear perturbative prescription even in the comoving case, where the non-linearity of inhomogeneities has an impact.

The climate crisis and the degradation of the world's ecosystems require humanity to take immediate action. The international scientific community has a responsibility to limit the negative environmental impacts of basic research. The HECAP+ communities (High Energy Physics, Cosmology, Astroparticle Physics, and Hadron and Nuclear Physics) make use of common and similar experimental infrastructure, such as accelerators and observatories, and rely similarly on the processing of big data. Our communities therefore face similar challenges to improving the sustainability of our research. This document aims to reflect on the environmental impacts of our work practices and research infrastructure, to highlight best practice, to make recommendations for positive changes, and to identify the opportunities and challenges that such changes present for wider aspects of social responsibility.