Papers for Tuesday, Aug 29 2023

Due to the massive increase in astronomical images (such as James Webb and Solar Dynamic Observatory), automatic image description is essential for solar and astronomical. Zernike moments (ZMs) are unique due to the orthogonality and completeness of Zernike polynomials (ZPs); hence valuable to convert a two-dimensional image to one-dimensional series of complex numbers. The magnitude of ZMs is rotation invariant, and by applying image normalization, scale and translation invariants can be made, which are helpful properties for describing solar and astronomical images. In this package, we describe the characteristics of ZMs via several examples of solar (large and small scale) features and astronomical images. ZMs can describe the structure and morphology of objects in an image to apply machine learning to identify and track the features in several disciplines.

Galactic aberration (GA) is a small effect in proper motions of celestial objects with an amplitude of about 5 $\mu$as/yr already noticeable in highly accurate astrometric observations such as VLBI and Gaia. However accurate accounting for this effect faces difficulty caused by the uncertainty in the GA amplitude (GA constant). Its estimates derived from VLBI and Gaia data processing differ significantly, so it would be very desirable to involve another independent method to solve the problem of inconsistency between these two methods. Such a method, that we consider in this paper, is using determination of the Galactic rotation parameters by methods of stellar astronomy. The result obtained in this study showed that the GA constant estimate obtained from stellar astronomy is closer to the estimate obtained from Gaia.

The primary objective of this thesis is to identify the key dynamical mechanisms responsible for the superthin stellar discs. We use HI 21cm radio-synthesis observations and stellar photometry to construct detailed dynamical models of a sample of superthin galaxies to determine the primary mechanism responsible for the existence of superthin stellar discs in these galaxies. Our study is based on a sample of superthin galaxies with $\rm 10< a/b < 16$ for which H1 21cm radio-synthesis data were already available in the literature. In addition, we had the two thinnest galaxies in our sample with $\rm a/b \sim 21$, for which we carried out Giant Meterwave Radio Telescope (GMRT) 21cm radio-synthesis observations. To identify the physical mechanism primarily responsible for the superthin vertical structure, we carry out a Principal Component Analysis of the following dynamical parameters: 1) Dark matter dominance at inner galactocentric radii given by $V_{\rm{rot}}/{(R_{c}/R_{d})}$, 2) the ratio of the vertical-to-radial stellar velocity dispersion $(\sigma_{z,s}/\sigma_{R,s})$, 3) Disc dynamical stability against local axisymmetric perturbations $Q_{RW}$, 4) Specific angular momentum of the disc $(j_{*})$, along with $a/b$ for all superthins and the extremely thin galaxies. We note that the first two principal components explain $\sim$ 80$\%$ of the variation in the data, and the major contributions are from $a/b$, $Q_{RW}$ and $V_{\rm{rot}}/{(R_{c}/R_{d})}$. This possibly indicates that high values of the disc dynamical stability and dark matter dominance at inner galactocentric radii are fundamentally responsible for the superthin stellar discs.

The abundance of protoplanetary bodies ejected from their parent star system is presently poorly-constrained. With only two existing optical observations of interstellar objects in the $10^{8} - 10^{10}$ kg mass range and a small number of robust microlensing observations of free-floating planets (FFPs) in the $10^{24} - 10^{25}$ kg mass range, there is a large range of masses for which there are no existing measurements of the unbound population. The three primary microlensing surveys currently searching for FFPs operate at a cadence greater than 15 minutes, which limits their ability to observe events associated with bodies with a mass much below an Earth mass. We demonstrate that existing high-cadence observations of M31 with the Subaru Hyper Suprime-Cam place the best direct constraints at present on the abundance of unbound objects at sub-terrestrial masses, with peak sensitivity at $10^{-4}~M_\oplus$ for Milky Way lenses and $10^{-1}~M_\oplus$ for lenses in M31. For a fiducial $\frac{dn}{dM}\propto M^{-2}$ mass distribution, we find that the abundance of unbound objects is constrained to $n_\text{unbound} < 4.4 \times 10^{7} ~\rm{pc}^{-3}$ for masses within 1 dex of $10^{-4}~M_\oplus$. Additionally, we compute limits on an artificial ``monochromatic'' distribution of unbound objects and compare to existing literature, demonstrating that the assumed spatial distribution of lenses has very significant consequences for the sensitivity of microlensing surveys. Our limits place robust constraints on the unbound population in the sub-terrestrial mass range and motivate new observational strategies for microlensing surveys.

We present radio observations of 23 optically-discovered tidal disruption events (TDEs) on timescales of about 500-3200 days post-discovery. We detect 9 new TDEs that did not have detectable radio emission at earlier times, indicating a late-time brightening after several hundred (and up to 2300 days); an additional 6 TDEs exhibit radio emission whose origin is ambiguous or may be attributed to the host galaxy or an AGN. We also report new rising components in two TDEs previously detected in the radio (iPTF16fnl and AT2019dsg) at ~1000 days. While the radio emission in some of the detected TDEs peaked on a timescale of ~2-4 years, more than half of the sample still shows rising emission. The range of luminosities for the sample is 10^37-10^39 erg/s, about two orders of magnitude below the radio luminosity of the relativistic TDE Sw1644+57. Our data set indicates that about 40% of all optical TDEs are detected in the radio hundreds to thousands of days after discovery, and that this is probably more common than early radio emission peaking at ~100 days. Using an equipartition analysis, we find evidence for a delayed launch of the radio-emitting outflows, with delay timescales of ~500-2000 days, inferred velocities of ~0.02-0.15c, and kinetic energies of ~10^47-10^49 erg. We rule out off-axis relativistic jets as a viable explanation for this population, and conclude delayed outflows are a more likely explanation, such as from delayed disk formation. Finally, we find comparable densities in the circumnuclear environments of these TDEs as for those with early radio emission, and find the TDEs still rising in luminosity are consistent with free expansion. We conclude that late radio emission marks a fairly ubiquitous but heretofore overlooked phase of TDE evolution.

During the journey from the cloud to the disc, the chemical composition of the protostellar envelope material can be either preserved or processed to varying degrees depending on the surrounding physical environment. This works aims to constrain the interplay of solid (ice) and gaseous methanol (CH$_3$OH) in the outer regions of protostellar envelopes located in the Coronet cluster in Corona Australis (CrA), and assess the importance of irradiation by the Herbig Ae/Be star R CrA. CH$_3$OH is a prime test-case as it predominantly forms as a consequence of the solid-gas interplay (hydrogenation of condensed CO molecules onto the grain surfaces) and it plays an important role in future complex molecular processing. We present 1.3 mm Submillimeter Array (SMA) and Atacama Pathfinder Experiment (APEX) observations towards the envelopes of four low-mass protostars in the Coronet. Eighteen molecular transitions of seven species are identified. We calculate CH$_3$OH gas-to-ice ratios in this strongly irradiated cluster and compare them with ratios determined towards protostars located in less irradiated regions such as the Serpens SVS 4 cluster in Serpens Main and the Barnard 35A cloud in the $\lambda$ Orionis region. The CH$_3$OH gas-to-ice ratios in the Coronet vary by one order of magnitude (from 1.2$\times$10$^{-4}$ to 3.1$\times$10$^{-3}$) which is similar to less irradiated regions as found in previous studies. We find that the CH$_3$OH gas-to-ice ratios estimated in these three regions are remarkably similar despite the different UV radiation field intensities and formation histories. This result suggests that the overall CH$_3$OH chemistry in the outer regions of low-mass envelopes is relatively independent of variations in the physical conditions and hence that it is set during the prestellar stage.

In this study, we modify the semi-analytic model Galacticus in order to accurately reproduce the observed properties of dwarf galaxies in the Milky Way. We find that reproducing observational determinations of the halo occupation fraction and mass-metallicity relation for dwarf galaxies requires us to include H$_2$ cooling, an updated UV background radiation model, and to introduce a model for the metal content of the intergalactic medium. By fine-tuning various model parameters and incorporating empirical constraints, we have tailored the model to match the statistical properties of Milky Way dwarf galaxies, such as their luminosity function and size$-$mass relation. We have validated our modified semi-analytic framework by undertaking a comparative analysis of the resulting galaxy-halo connection. We predict a total of $300 ^{+75} _{-99}$ satellites with an absolute $V$-band magnitude (M$_{V}$) less than $0$ within $300$ kpc from our Milky Way-analogs. The fraction of subhalos that host a galaxy at least this bright drops to $50\%$ by a halo peak mass of $\sim 8.9 \times 10^{7}$ M$_{\odot}$, consistent with the occupation fraction inferred from the latest observations of Milky Way satellite population.

Detecting the end of the cosmic-ray (CR) electron spectrum would provide important new insights. While we know that Milky Way sources can accelerate electrons up to at least $\sim$1~PeV, the observed CR electron spectrum at Earth extends only up to 5~TeV (possibly 20~TeV), a large discrepancy. The question of the end of the CR electron spectrum has received relatively little attention, despite its importance. We take a comprehensive approach, showing that there are multiple steps at which the observed CR electron spectrum could be cut off. At the highest energies, the accelerators may not have sufficient luminosity, or the sources may not allow sufficient escape, or propagation to Earth may not be sufficiently effective, or present detectors may not have sufficient sensitivity. For each step, we calculate a rough range of possibilities. Although all of the inputs are uncertain, a clear vista of exciting opportunities emerges. We outline strategies for progress based on CR electron observations and auxiliary multi-messenger observations. In addition to advancing our understanding of CRs in the Milky Way, progress will also sharpen sensitivity to dark matter annihilation or decay.

In the present study, we identified an extremely low-mass white dwarf as a companion to a low luminous blue straggler star within the Galactic globular cluster NGC 362. To conduct the analysis, we utilized data obtained from various sources, including AstroSat Ultra Violet Imaging Telescope, UVOT, and the 2.2-m ESO telescope. By examining the spectral energy distribution of the blue straggler star candidate, we successfully identified an extremely low mass white dwarf as its binary companion. We determined the effective temperature, radius, and luminosity for both, the extremely low-mass white dwarf and the blue straggler star candidate. Specifically, the extremely low-mass white dwarf exhibits an effective temperature of 15000 K, a radius of 0.17 Rsun, a luminosity of 1.40 Lsun, and a mass range of 0.19-0.20 Msun. Furthermore, the position of the straggler star within the cluster suggests their formation via the Case A/B mass-transfer mechanism in a low-density environment.

Lava planets are rocky exoplanets that orbit so close to their host star that their day-side is hot enough to melt silicate rock. Their short orbital periods ensure that lava planets are tidally locked into synchronous rotation, with permanent day and night hemispheres. Such asymmetric magma oceans have no analogs in the Solar System and will exhibit novel fluid dynamics. Here we report numerical simulations of lava planet interiors showing that solid-liquid fractionation in the planetary interior has a major impact on the compositional structure and evolution of the planet. We explored two styles of dynamics that depend primarily on the interior thermal state : 1) a hot fully molten interior, and 2) a mostly solid interior with a shallow day-side magma ocean. In the hot interior scenario, the atmosphere reflects the planet's bulk silicate composition and the night-side crust is gravitationally unstable and constantly replenished. In the cool interior scenario, the distilled atmosphere will lack Na, K and FeO, and the night-side mantle is entirely solid, with a cold surface. These two end-member cases can be distinguished with observations from the James Webb Space Telescope, offering an avenue to probe the diversity of terrestrial exoplanet evolutions.

We demonstrate how a changing and negative equation-of-state (EoS) for dark matter can alleviate cosmic tensions and explain the integrated Sachs-Wolfe (ISW) void anomaly. We discuss the effect of the model on the cosmic expansion history, growth of structure and the ISW. We show that a negative EoS at late times is able to produce a larger Hubble constant and smaller $\sigma_{8}$, which can explain both cosmological tensions. Furthermore, the model uniquely predicts larger ISW at low redshift, a prediction which is in agreement with observations of larger ISW from voids. The preference for a negative EoS for dark matter at late times is indicative of a unified dark sector and degenerate with models of dark matter and dark energy interaction. Future measurements of the ISW from cosmic voids can provide a unique test for this solution to tensions in cosmology, should they continue to persist.

We present a thermal observation of Callisto's leading hemisphere obtained using the Atacama Large Millimeter/submillimeter Array (ALMA) at 0.87 mm (343 GHz). The angular resolution achieved for this observation was $\sim$$0.16^{\prime\prime}$, which for Callisto at the time of this observation ($D\sim 1.05^{\prime\prime}$) was equivalent to $\sim$6 elements across the surface. Our disk-integrated brightness temperature of 116 $\pm$ 5 K (8.03 $\pm$ 0.40 Jy) is consistent with prior disk-integrated observations. Global surface properties were derived from the observation using a thermophysical model (de Kleer et al. 2021) constrained by spacecraft data. We find that models parameterized by two thermal inertia components more accurately fit the data than single thermal inertia models. Our best-fit global parameters adopt a lower thermal inertia of 15-50 $\text{J}\:\text{m}^{-2}\:\text{K}^{-1}\:\text{s}^{-1/2}$ and a higher thermal inertia component of 1200-2000 $\text{J}\:\text{m}^{-2}\:\text{K}^{-1}\:\text{s}^{-1/2}$, with retrieved millimeter emissivities of 0.89-0.91. We identify several thermally anomalous regions, including spots $\sim$3 K colder than model predictions co-located with the Valhalla impact basin and a complex of craters in the southern hemisphere; this indicates the presence of materials possessing either a higher thermal inertia or a lower emissivity. A warm region confined to the mid-latitudes in these leading hemisphere data may be indicative of regolith property changes due to exogenic sculpting.

The population of strongly irradiated Jupiter-sized planets has no equivalent in the Solar System. It is characterised by strongly bloated atmospheres and atmospheric large-scale heights. Recent space-based observations of SO2 photochemistry demonstrated the knowledge that can be gained from detailed atmospheric studies of these unusual planets about Earth's uniqueness. Aims. Here we explore the atmosphere of WASP-172b a similar planet in temperature and bloating to the recently studied HD~149026~b. In this work, we characterise the atmospheric composition and subsequently the atmospheric dynamics of this prime target. Methods. We observed a particular transit of WASP-172b in front of its host star with ESO's ESPRESSO spectrograph and analysed the spectra obtained before during and after transit. Results. We detect the absorption of starlight by WASP-172b's atmosphere by sodium (5.6sigma), hydrogen (19.5sigma) and obtained a tentative detection of iron (4.1sigma). We detect strong - yet varying - blue shifts, relative to the planetary rest frame, of all of these absorption features. This allows for a preliminary study of the atmospheric dynamics of WASP-172b. Conclusions. With only one transit, we were able to detect a wide variety of species, clearly tracking different atmospheric layers with possible jets. WASP-172b is a prime follow-up target for a more in-depth characterisation both for ground and space-based observatories. If the detection of Fe is confirmed, this may suggest that radius inflation is an important determinant for the detectability of Fe in hot Jupiters, as several non-detections of Fe have been published for planets that are hotter but less inflated than WASP-172b.

The Magellanic Cloud system represents a unique laboratory for study of both interacting dwarf galaxies and the ongoing process of the formation of the Milky Way and its halo. We focus on one aspect of this complex, 3 body interaction - the dynamical perturbation of the Small Magellanic Cloud (SMC) by the Large Magellanic Cloud (LMC), and specifically potential tidal effects on the SMC's eastern side. Using Gaia astrometry and the precise radial velocities and multi-element chemical abundances from APOGEE-2 DR17, we explore the well-known distance bimodality on the eastern side of the SMC. Through estimated stellar distances, proper motions, and radial velocities, we characterize the kinematics of the two populations in the bimodality and compare their properties with those of SMC populations elsewhere. Moreover, while all regions explored by APOGEE seem to show a single chemical enrichment history, the metallicity distribution function (MDF), of the "far" stars on the eastern periphery of the SMC is found to resemble that for the more metal-poor fields of the western periphery, whereas the MDF for the "near" stars on the eastern periphery resembles that for stars in the SMC center. The closer eastern periphery stars also show radial velocities (corrected for SMC rotation and bulk motion) that are, on average, approaching us relative to all other SMC populations sampled. We interpret these trends as evidence that the near stars on the eastern side of the SMC represent material pulled out of the central SMC as part of its tidal interaction with the LMC.

Hubble Space Telescope images obtained in 2018 are combined with archival HST data taken in 1995 to detect changes and measure proper motions in the HH 80/81 shock complex which is powered by the fastest known jet driven by a forming star, the massive object IRAS 18162-2048. Some persistent features close to the radio jet axis have proper motions grater than 1,000 km/s away from IRAS 18162-2048. About 3 to 5 parsecs downstream from the IRAS source and beyond HH 80/81, H-alpha emission traces the rim of a parsec-scale bubble blown by the jet. Lower speed motions are seen in [Sii] away from the jet axis; these features have a large component of motion at right-angles to the jet. We identify new HH objects and H2 shocks in the counterflow opposite HH 80/81. The northeastern counterflow to HH 80/81 exhibits an extended but faint complex of 2.12 um H2 shocks. The inner portion of the outflow is traced by dim 1.64 um [Feii] emission. The full extent of this outflow is at least 1,500" (about 10 pc in projection at a distance of 1.4 kpc). We speculate about the conditions responsible for the production of the ultra-fast jet and the absence of prominent large-scale molecular outflow lobes.

Solar Wind Charge eXchange X-ray (SWCX) emission in the heliosphere and Earth's exosphere is a hard to avoid signal in soft X-ray observations of astrophysical targets. On the other hand, the X-ray imaging possibilities offered by the SWCX process has led to an increasing number of future dedicated space missions for investigating the solar wind-terrestrial interactions and magnetospheric interfaces. In both cases, accurate modelling of the SWCX emission is key to correctly interpret its signal, and remove it from observations, when needed. In this paper, we compile solar wind abundance measurements from ACE for different solar wind types, and atomic data from literature, including charge exchange cross-sections and emission probabilities, used for calculating the compound cross-section $\alpha$ for the SWCX X-ray emission. We calculate $\alpha$ values for charge-exchange with H and He, relevant to soft X-ray energy bands (0.1 - 2.0 keV) for various solar wind types and solar cycle conditions.

Forward-modeling observables from galaxy simulations enables direct comparisons between theory and observations. To generate synthetic spectral energy distributions (SEDs) that include dust absorption, re-emission, and scattering, Monte Carlo radiative transfer is often used in post-processing on a galaxy-by-galaxy basis. However, this is computationally expensive, especially if one wants to make predictions for suites of many cosmological simulations. To alleviate this computational burden, we have developed a radiative transfer emulator using an artificial neural network (ANN), ANNgelina, that can reliably predict SEDs of simulated galaxies using a small number of integrated properties of the simulated galaxies: star formation rate, stellar and dust masses, and mass-weighted metallicities of all star particles and of only star particles with age <10 Myr. Here, we present the methodology and quantify the accuracy of the predictions. We train the ANN on SEDs computed for galaxies from the IllustrisTNG project's TNG50 cosmological magnetohydrodynamical simulation. ANNgelina is able to predict the SEDs of TNG50 galaxies in the ultraviolet (UV) to millimetre regime with a typical median absolute error of ~7 per cent. The prediction error is the greatest in the UV, possibly due to the viewing-angle dependence being greatest in this wavelength regime. Our results demonstrate that our ANN-based emulator is a promising computationally inexpensive alternative for forward-modeling galaxy SEDs from cosmological simulations.

Very-high-energy and ultra-high-energy neutrinos are messengers of energetic sources in the universe. Sub-orbital and satellite-based neutrino telescopes employ detectors of the atmospheric Cherenkov emission from extensive air showers (EASs) generated by charged particles. These Cherenkov detectors can be pointed below or above the Earth's limb. Cherenkov emissions produced from directions below the limb are from upward-going EASs produced in the atmosphere sourced by Earth-skimming neutrinos. When the Cherenkov telescope is pointed slightly above the Earth's limb, signals from EASs are initiated by cosmic ray interactions in the atmosphere. For sub-orbital detectors, muons produced from cosmic rays in the atmosphere can directly hit the Cherenkov telescope. Using a semi-analytic technique with cascade equations for atmospheric particle fluxes, we quantify the atmospheric muon flux that reaches sub-orbital telescopes like Extreme Universe Space Observatory Super Pressure Balloon 2 (EUSO-SPB2). We assess this potential background to the EAS signals. The calculation technique may also provide an understanding of the evolution of the muon content in individual EAS.

We present Fermi Gamma-ray Burst Monitor (Fermi-GBM) and Swift Burst Alert Telescope (Swift-BAT) searches for gamma-ray/X-ray counterparts to gravitational wave (GW) candidate events identified during the third observing run of the Advanced LIGO and Advanced Virgo detectors. Using Fermi-GBM on-board triggers and sub-threshold gamma-ray burst (GRB) candidates found in the Fermi-GBM ground analyses, the Targeted Search and the Untargeted Search, we investigate whether there are any coincident GRBs associated with the GWs. We also search the Swift-BAT rate data around the GW times to determine whether a GRB counterpart is present. No counterparts are found. Using both the Fermi-GBM Targeted Search and the Swift-BAT search, we calculate flux upper limits and present joint upper limits on the gamma-ray luminosity of each GW. Given these limits, we constrain theoretical models for the emission of gamma-rays from binary black hole mergers.

Waste gas products from technological civilizations may accumulate in an exoplanet atmosphere to detectable levels. We propose nitrogen trifluoride (NF3) and sulfur hexafluoride (SF6) as ideal technosignature gases. Earth life avoids producing or using any N-F or S-F bond-containing molecules and makes no fully fluorinated molecules with any element. NF3 and SF6 may be universal technosignatures owing to their special industrial properties, which unlike biosignature gases, are not species-dependent. Other key relevant qualities of NF3 and SF6 are: their extremely low water solubility, unique spectral features, and long atmospheric lifetimes. NF3 has no non-human sources and was absent from Earth's pre-industrial atmosphere. SF6 is released in only tiny amounts from fluorine-containing minerals, and is likely produced in only trivial amounts by volcanic eruptions. We propose a strategy to rule out SF6's abiotic source by simultaneous observations of SiF4, which is released by volcanoes in an order of magnitude higher abundance than SF6. Other fully fluorinated human-made molecules are of interest, but their chemical and spectral properties are unavailable. We summarize why life on Earth-and perhaps life elsewhere-avoids using F. We caution, however, that we cannot definitively disentangle an alien biochemistry byproduct from a technosignature gas.

The famous recurrent nova (RN) T Coronae Borealis (T CrB) has had observed eruptions peaking at a visual magnitude of 2.0 in the years 1866 and 1946, while a third eruption is now expected for the year 2024.4+-0.3. Each RN has very similar light curves of eruptions that come with a fairly even-spacing in time, for which T CrB has a recurrence timescale near 80 years. So it is reasonable to look backwards in time for prior eruptions, around 1786, and so on back. I have investigated two long-lost suggestions that T CrB was seen in eruption in the years 1217 and 1787. (1) In a catalog published in 1789, the Reverend Francis Wollaston reports an astrometric position for a star that is exactly on top of T CrB. From his letters, these observations were made on at least four occassions with both a large and small telescope, within a few days before 1787 December 28. Wollaston's limiting magnitude for his astrometry is near 7.8 mag, so T CrB would have to have been in eruption. With other transients strongly rejected, the only way that Wollaston could get the coordinates was to have measured the coordinates of T CrB itself during an eruption. (2) The 1217 event has an eyewitness report written by Abbott Burchard of Upsberg as a fast-rising stellar point-source ("stella") in Corona Borealis that "shone with great light", lasted for "many days", and was ascribed as being a "wonderful sign". This event cannot be a report of a comet, because Burchard used the term for a star ("stella") and not for a comet, and because Burchard had the omen being very positive, with such being impossible for comets that are universally the worst of omens. The reported event is just as expected for a prior eruption of T CrB, and all other possibilities are strongly rejected, so the case for the 1217 eruption of T CrB is strong.

We report on the observation of a rare neutrino event detected by Baikal-GVD in April 2021. The event GVD210418CA is the highest-energy cascade observed by Baikal-GVD so far from the direction below the horizon. The estimated cascade energy is $224\pm75$~TeV. The evaluated signalness parameter of GVD210418CA is 97.1\% using an assumption of the E$^{-2.46}$ spectrum of astrophysical neutrinos. The arrival direction of GVD210418CA is near the position of the well-known radio blazar TXS~0506+056, with the angular distance being within a 90\% directional uncertainty region of the Baikal-GVD measurement. The event was followed by a radio flare observed by the RATAN-600 radio telescope, further strengthening the case for the neutrino-blazar association.

Massive black holes (MBHs) are crucial in shaping their host galaxies. How the MBH co-evolves with its host galaxy is a pressing problem in astrophysics and cosmology. The valuable information carried by the binary MBH is encoded in the gravitational waves (GWs), which will be detectable by the space-borne GW detector LISA. In the GW data analysis, usually, only the dominant $(2,2)$ mode of the GW signal is considered in the parameter estimation for LISA. However, including the higher harmonics in parameter estimation can break the degeneracy between the parameters, especially for the inclination angle and luminosity distance. This may enable the identification of GW signals without electromagnetic counterparts, known as ''dark sirens''. Thus, incorporating higher harmonics will be beneficial to resolve the Hubble tension and constrain the cosmological model. In this paper, we investigate the role of higher harmonics in the parameter estimation for GWs emitted by binary MBHs. We demonstrate that including $(3,3)$ mode can lead to a $10^3$-times improvement in angular resolution and a $10^4$-times improvement in luminosity distance. Meanwhile, our results indicate that considering higher harmonics increases the probability of identifying over 70% host galaxies from $10^{-2}\,\rm{Gpc}^3$ cosmological volume threshold (corresponding $10^5$ host galaxies), while the probability less than 8% for only the $(2,2)$ mode. Thus, our results underscore the importance of including higher modes in the GW signal from binary MBHs, for LISA at least $(3,3)$ mode.

Core-collapse explosions of massive stars leave behind neutron stars, with a known diversity that includes the "Central Compact Objects" (CCOs). Typified by the neutron star discovered near the centre of the Cas A supernova remnant (SNR), CCOs have been observed to shine only in X-rays. To address their supernova progenitors, we perform a systematic study of SNRs that contain a CCO and display X-ray emission from their shock-heated ejecta. We make use of X-ray data primarily using the Chandra X-ray observatory, complemented with XMM-Newton. This study uses a systematic approach to the analysis of each SNR aimed at addressing the supernova progenitor as well as the explosion properties (energy and ambient density). After fitting for the ejecta abundances estimated from a spatially resolved spectroscopic study, we compare the data to six nucleosynthesis models making predictions on supernova ejecta yields in core-collapse explosions. We find that the explosion models commonly used by the astrophysics community do not match the ejecta yields for any of the SNRs, suggesting additional physics, e.g. multi-dimensional explosion models or updated progenitor structures, are required. Overall we find low-mass ($\leq$25 solar masses) progenitors among the massive stars population and low-energy explosions ($<$10$^{51}$ ergs). We discuss degeneracies in our model fitting, particularly how altering the explosion energy affects the estimate of the progenitor mass. Our systematic study highlights the need for improving on the theoretical models for nucleosynthesis predictions as well as for sensitive, high-resolution spectroscopy observations to be acquired with next-generation X-ray missions.

We present the first detailed chemical-abundance analysis of stars from the dwarf-galaxy stellar stream Wukong/LMS-1 covering a wide metallicity range ($-3.5 < \rm[Fe/H] \lesssim -1.3$). We find abundance patterns that are effectively indistinguishable from the bulk of Indus and Jhelum, a pair of smaller stellar streams proposed to be dynamically associated with Wukong/LMS-1. We discovered a carbon-enhanced metal-poor star ($\rm[C/Fe] > +0.7$ and $\rm[Fe/H] \sim -2.9$) in Wukong/LMS-1 with strong enhancements in Sr, Y, and Zr, which is peculiar given its solar-level [Ba/Fe]. Wukong/LMS-1 stars have high abundances of $\alpha$ elements up to $\rm[Fe/H] \gtrsim -2$, which is expected for relatively massive dwarfs. Towards the high-metallicity end, Wukong/LMS-1 becomes $\alpha$-poor, revealing that it probably experienced fairly standard chemical evolution. We identified a pair of N- and Na-rich stars in Wukong/LMS-1, reminiscent of multiple populations in globular clusters. This indicates that this dwarf galaxy contained at least one globular cluster that was completely disrupted in addition to two intact ones previously known to be associated with Wukong/LMS-1, which is possibly connected to similar evidence found in Indus. From these $\geq$3 globular clusters, we estimate the total mass of Wukong/LMS-1 to be ${\approx}10^{10} M_\odot$, representing ${\sim}1$% of the present-day Milky Way. Finally, the [Eu/Mg] ratio in Wukong/LMS-1 continuously increases with metallicity, making this the first example of a dwarf galaxy where the production of $r$-process elements is clearly dominated by delayed sources, presumably neutron-star mergers.

Mini-EUSO is the first mission of the JEM-EUSO program on board the International Space Station. It was launched in August 2019 and it is operating since October 2019 being located in the Russian section (Zvezda module) of the station and viewing our planet from a nadir-facing UV-transparent window. The instrument is based on the concept of the original JEM-EUSO mission and consists of an optical system employing two Fresnel lenses of 25 cm each and a focal surface composed of 36 Multi-Anode Photomultiplier tubes, 64 channels each, for a total of 2304 channels with single photon counting sensitivity and an overall field of view of 44$\times$44$^\circ$. Mini-EUSO can map the night-time Earth in the near UV range (predominantly between 290 nm and 430 nm), with a spatial resolution of about 6.3 km and different temporal resolutions of 2.5 $\mu$s, 320 $\mu$s and 41 ms. Mini-EUSO observations are extremely important to better assess the potential of a space-based detector in studying Ultra-High Energy Cosmic Rays (UHECRs) such as K-EUSO and POEMMA. In this contribution we focus the attention on the results of the UV measurements and we place them in the context of UHECR observations from space, namely the estimation of exposure.

Very young massive planets are sufficiently luminous by their internal heat of formation to permit detailed studies, including spectroscopy of their atmospheres with large telescopes at sufficient resolution ($\lambda / \Delta \lambda \gtrsim 1000$) to identify major constituents to inform models of planet formation and early evolution. We obtained 1-2.4$\mu$m ($YJHK$) spectra of the planetary-mass "b" companion of 2MASS~J04372171+2651014, a 1-3 Myr-old M dwarf member of the Taurus star-forming region, and one of the youngest such objects discovered to date. These indicate the presence of CO and possibly H$_2$O and CH$_4$ in the atmosphere, all suggesting a $T_{\rm eff}$ of around 1200K, characteristic of a L-T transition spectral type and consistent with previous estimates based on its luminosity and age. The absence or attenuation of spectral features at shorter wavelengths suggests the presence of micron-size dust, consistent with the object's red color. The spectrum of 2M0437b resembles those of the HR 8799 planets, especially the innermost "b" planet, with the exception of a pronounced flux deficit in the $H$-band of uncertain origin.

The orbits of trans-Neptunian objects (TNOs) can indicate the existence of an undiscovered planet in the outer solar system. Here, we used N-body computer simulations to investigate the effects of a hypothetical Kuiper Belt planet (KBP) on the orbital structure of TNOs in the distant Kuiper Belt beyond ~50 au. We used observations to constrain model results, including the well-characterized Outer Solar System Origins Survey (OSSOS). We determined that an Earth-like planet (m ~ 1.5-3 Earth masses) located on a distant (semimajor axis a ~ 250-500 au, perihelion q ~ 200 au) and inclined (i ~ 30 deg) orbit can explain three fundamental properties of the distant Kuiper Belt: a prominent population of TNOs with orbits beyond Neptune's gravitational influence (i.e., detached objects with q > 40 au), a significant population of high-i objects (i > 45 deg), and the existence of some extreme objects with peculiar orbits (e.g., Sedna). Furthermore, the proposed KBP is compatible with the existence of identified Gyr-stable TNOs in the 2:1, 5:2, 3:1, 4:1, 5:1, and 6:1 Neptunian mean motion resonances. These stable populations are often neglected in other studies. We predict the existence of an Earth-like planet and several TNOs on peculiar orbits in the outer solar system, which can serve as observationally testable signatures of the putative planet's perturbations.

A weak absorption feature at 2.07 \textmu m on Europa's trailing hemisphere has been suggested to arise from radiolytic processing of an endogenic salt, possibly sourced from the interior ocean. However, if the genesis of this feature requires endogenic material to be present, one might expect to find a correlation between its spatial distribution and the recently disrupted chaos terrains. Using archived near-infrared observations from Very Large Telescope/SINFONI with a $\sim$1 nm spectral resolution and a linear spatial resolution $\sim$130 km, we examine the spatial distribution of this feature in an effort to explore this endogenic formation hypothesis. We find that while the presence of the 2.07 \textmu m feature is strongly associated with the irradiation pattern on Europa's trailing hemisphere, there is no apparent association between the presence or depth of the absorption feature and Europa's large-scale chaos terrain. This spatial distribution suggests that the formation pathway of the 2.07 \textmu m feature on Europa is independent of any endogenous salts within the recent geology. Instead, we propose that the source of this feature may simply be a product of the radiolytic sulfur cycle or arise from some unidentified parallel irradiation process. Notably, the 2.07 \textmu m absorption band is absent from the Pwyll crater ejecta blanket, suggesting that radiolytic processing has not had enough time to form the species responsible and placing a lower limit on the irradiation timescale. We are unable to find a plausible spectral match to the 2.07 \textmu m feature within the available laboratory data.

We investigate which scalar quantity or quantities can best predict the loss of equilibrium and subsequent eruption of magnetic flux ropes in the solar corona. Our models are initialized with a potential magnetic arcade, which is then evolved by means of two effects on the lower boundary: firstly a gradual shearing of the arcade, modelling differential rotation on the solar surface, and secondly supergranular diffusion. These result in flux cancellation at the polarity inversion line and the formation of a twisted flux rope. We use three model setups: full magnetohydrodynamics (MHD) in cartesian coordinates, and the magnetofrictional model in both cartesian and polar coordinates. The flux ropes are translationally-invariant, allowing for very fast computational times and thus a comprehensive parameter study, comprising hundreds of simulations and thousands of eruptions. Similar flux rope behavior is observed using either magnetofriction or MHD, and there are several scalar criteria that could be used as proxies for eruptivity. The most consistent predictor of eruptions in either model is the squared current in the axial direction of the rope, normalised by the relative helicity, although a variation on the previously proposed 'eruptivity index' is also found to perform well in both the magnetofrictional and MHD simulations.

Near-Earth asteroids Ryugu and Bennu, were visited, characterised, and sampled by the Hayabusa2 and OSIRIS-REx missions: remote sensing data and sample return analysis showed that both asteroids have primitive, hydrated and organic-rich compositions. The dark families of the inner main belt (IMB) that belong to the spectroscopic C-complex have been claimed to be the sources of both Ryugu and Bennu. Hence, there has been large effort to characterise them. Here we used the Gaia Data Release 3 (DR3) asteroid reflectance spectra to investigate the 11 known IMB C-complex families (Chaldaea, Chimaera, Clarissa, Erigone, Eulalia, Klio, Polana, Primordial, Sulamitis, Svea, Tamara). For each family, we extracted the family members that have known geometric visible albedo values and Gaia DR3 data and we created an average reflectance spectrum per family between 370 and 950 nm. The average DR3 reflectance spectra of each family were compared with the previous literature data and to Bennu's and Ryugu's spectra. We found that DR3 reflectance spectra of the IMB C-complex families are in general consistent with previous findings with the only exception of the Svea family. We also showed that the Polana and the Eulalia families can be distinguished in the wavelength region 370 - 500 nm. Among all the IMB C-complex families, we determined that the average reflectance spectra of the Eulalia and Polana families are the most similar to those of Bennu and Ryugu, respectively. In particular, Eulalia family's average spectrum is a good match to Bennu's in the wavelength range 450 - 800 nm, while beyond 800 nm the spectrum of Bennu is bluer than that of Eulalia. Moreover, the spectrum of the Polana family has the smallest discrepancy against the spectrum of Ryugu, although this match is formally unsatisfactory (reduced chi^2 ~ 1.9).

The Baikal-GVD alert system was launched at the beginning of 2021. There are alerts for muon neutrinos (long upward-going track-like events) and all-flavour neutrinos (high-energy cascades). The system is able to get a preliminary response to external alerts with a temporal delay of about 3-10 minutes. The Baikal-GVD data processing and the results of the follow-up procedure are described. We report on the analysis of the coincidence in time and direction between the Baikal-GVD cascade GVD20211208CA with an estimated energy of 43 TeV and the announced alert IceCube211208A possibly associated with a flaring state of the blazar PKS 0735+178.

We present counts-level fits to the keV-GeV data of the early afterglow of the brightest gamma-ray burst detected to date, GRB 221009A. We discuss the complexity of the data reduction due to the unprecedented brightness and the location in the Galactic plane. We find the energy spectrum to be well described as a smoothly broken power law with a break around 10 keV and no indications for additional features towards GeV energies. An interpretation as synchrotron emission from forward-shock accelerated and subsequently cooled electrons yields three possible types of solutions: (1) a slow cooling solution with low magnetic fields (few percent of a Gauss) but poorly constrained minimum injected electron energy (<100 GeV), (2) a fast cooling solution with stronger magnetic fields (few percent to few Gauss) and minimum injected electron energy 10-100 GeV or (3) the transition between both regimes with low magnetic fields and minimum injected electron energy around 100 GeV. Limited statistics at GeV energies prevent conclusions extrapolating towards a cut-off or the onset of a new component at higher energies.

We extend the analytic theory presented by Sternberg et al. 2014 and Bialy \& Sternberg (2016) for the production of atomic hydrogen (HI) via FUV photodissociation at the boundaries of dense interstellar molecular (H$_2$) clouds, to also include the effects of penetrating (low-energy) cosmic-rays for the growth of the total HI column densities. We compute the steady-state abundances of the HI and H$_2$ in one-dimensional gas slabs in which the FUV photodissociation rates are reduced by depth-dependent H$_2$ self-shielding and dust absorption, and in which the cosmic-ray ionization rates are either constant or reduced by transport effects. The solutions for the HI and H$_2$ density profiles and the integrated HI columns, depend primarily on the ratios $I_{\rm UV}/Rn$ and $\zeta/Rn$, where $I_{\rm UV}$ is the intensity of the photodissociating FUV field, $\zeta$ is the H$_2$ cosmic-ray ionization rate, $n$ is the hydrogen gas density, and $R$ is the dust-surface H$_2$ formation rate coefficient. We present computations for a wide range of FUV field strengths, cosmic-ray ionization rates, and dust-to-gas ratios. We develop analytic expressions for the growth of the HI column densities. For Galactic giant molecular clouds (GMCs) with multiphased (warm/cold) HI envelopes, the interior cosmic-ray zones will dominate the production of the HI only if $\zeta \gtrsim 4.5\times 10^{-16} \times (M_{\rm GMC}/10^6 \ M_{\odot})^{-1/2}$~s$^{-1}$, where $M_{\rm GMC}$ is the GMC mass, and including attenuation of the cosmic-ray fluxes. For most Galactic GMCs and conditions, FUV photodissociation dominates over cosmic-ray ionization for the production of the HI column densities. Furthermore, the cosmic-rays do not affect the HI-to-H$_2$ transition points.

The planetary obliquity plays a significant role in determining physical properties of planetary surfaces and climate. As direct detection is constrained due to the present observation accuracy, kinetic theories are helpful to predict the evolution of the planetary obliquity. Here the coupling effect between the eccentric Kozai-Lidov (EKL) effect and the equilibrium tide is extensively investigated, the planetary obliquity performs to follow two kinds of secular evolution paths, based on the conservation of total angular momentum. The equilibrium timescale of the planetary obliquity $t_{\mathrm{eq}}$ varies along with $r_{t}$, which is defined as the initial timescale ratio of the tidal dissipation and secular perturbation. We numerically derive the linear relationship between $t_{\mathrm{eq}}$ and $r_{t}$ with the maximum likelihood method. The spin-axis orientation of S-type terrestrials orbiting M-dwarfs reverses over $90^\circ$ when $r_{t} > 100$, then enter the quasi-equilibrium state between $40^\circ$ and $60^\circ$, while the maximum obliquity can reach $130^\circ$ when $r_{t} > 10^4 $. Numerical simulations show that the maximum obliquity increases with the semi-major axis ratio $a_1$/$a_2$, but is not so sensitive to the eccentricity $e_2$. The likelihood of obliquity flip for S-type terrestrials in general systems with $a_2 < 45$ AU is closely related to $m_1$. The observed potential oblique S-type planets HD 42936 b, GJ 86 Ab and $\tau$ Boot Ab are explored to have a great possibility to be head-down over the secular evolution of spin.

We study the long-term behavior of the bright gamma-ray blazar PKS 0402-362. Over a span of approximately 12.5 years, from August 2008 to January 2021, we gathered Fermi-LAT temporal data and identified three distinct periods of intense $\gamma$-ray activity. Notably, the second period exhibited the highest brightness ever observed in this particular source. We observed most of the $\gamma$-ray flare peaks to be asymmetric in profile suggesting a slow cooling time of particles or the varying Doppler factor as the main cause of these flares. The $\gamma$-ray spectrum is fitted with power-law and log-parabola models, and in both cases, the spectral index is very steep. The lack of time lags between optical-IR and $\gamma$-ray emissions indicates the presence of a single-zone emission model. Using this information, we modeled the broadband SEDs with a simple one-zone leptonic model using the publicly available code `GAMERA'. The particle distribution index is found to be the same as expected in diffusive shock acceleration suggesting it as the main mechanism of particle acceleration to very high energy up to 4 - 6 GeV. During the different flux phases, we observed that the thermal disk dominates the optical emission, indicating that this source presents a valuable opportunity to investigate the connection between the disk and the jet.

Fitting model spectral energy distributions (SED) to galaxy photometric data is a widely used method to recover galaxy parameters from galaxy surveys. However, the parameter space used to describe galaxies is wide and interdependent, and distinctions between real and spurious correlations that are found between these parameters can be difficult to discern. In this work, we use the SED fitting code BAGPIPES to investigate degeneracies between galaxy parameters and the effect of the choice of different sets of photometric bands. In particular, we focus on optical to infrared wavelength coverage, and on two parameters describing the galaxies' dust attenuation law: $A_V$ and $\delta$, which characterize dust column density and the slope of a flexible dust attenuation law, respectively. We demonstrate that 1) a degeneracy between the residual (the difference between truth and recovered value) $A_V$ and star formation rate exists, but this is lifted when WISE bands are included; 2) BAGPIPES is able to accurately recover the input $A_V$ and $\delta$ distributions and relations (differences in slope of less than 1.7$\sigma$ for a flat relation, less than 1.2$\sigma$ for an observationally-motivated relation from Salim et al. 2018) and is not introducing spurious correlations between these parameters. Our findings suggest that the information needed to constrain $A_V$ and $\delta$ well enough individually exists in the data, especially when IR is added. This indicates that recent works finding a correlation between $A_V$ and $\delta$ are not being misled by fitting degeneracies from their SED fitting code.

We present the results of determining the parameters of the spiral arms of the Galaxy using the stars Gaia DR3, whose absolute magnitude is $M_G$ < 4, and which allow tracing spiral arms at large distances from the Sun. As tracers of spiral arms, we use the centroids of stellar spherical regions with a radius of 0.5 kpc, in which the deformation velocities along the coordinate axis R are insignificant. These kinematic tracers cover the Galactic plane within the Galactocentric coordinate ranges 140{\deg} < ${\theta}$ < 220{\deg} and 4 kpc < R < 14 kpc. The numerical values of the pitch angles of the spirals and their Galactocentric distances to the point of intersection of the spiral with the direction of the Galactic center - the Sun are in good agreement with the results of other authors. By extrapolating beyond the data we have, we present a schematic four-arm global pattern, consisting of the Scutum-Centaurus, Sagitarius-Carina, Perseus, Norma-Outer arms, as well as the local arm Orion. The uncertainties of the determined spiral parameters confirm that the structures identified are not false, but are reliable from the statistical point of view.

Interstellar lines observed toward stellar targets change slowly over long timescales, mainly due to the proper motion of the background target relative to the intervening clouds. On longer timescales, the cloud's slowly changing physical and chemical conditions can also cause variation. We searched for systematic variations in the absorption profiles of the diffuse interstellar bands (DIBs) and interstellar atomic and molecular lines by comparing the high-quality data set from the ESO diffuse interstellar bands extensive exploration survey (EDIBLES) to older archival observations, bridging typical timescales of 10 years with a maximum timescale of 22 years. We found good archival observations for 64 EDIBLES targets. Our analysis focused on 31 DIBs, 7 atomic, and 5 molecular lines. We considered various systematic effects and applied a robust Bayesian test to establish which absorption features could display significant variations. While systematic effects greatly complicate our search, we find evidence for variations in the profiles of the $\lambda\lambda$4727 and 5780 DIBs in a few sightlines. Toward HD~167264, we find a new \ion{Ca}{i} cloud component that appears and becomes stronger after 2008. The same sightline furthermore displays marginal but systematic changes in the column densities of the atomic lines originating from the leading cloud component in the sightline. Similar variations are seen toward HD~147933. Our high-quality spectroscopic observations and archival data show that it is possible to probe interstellar time variations on time scales of typically a decade. Even though systematic uncertainties, as well as the generally somewhat lower quality of older data, complicate matters, we can conclude that time variations can be made visible, both in atomic lines and DIB profiles for a few targets, but that generally, these features are stable along many lines of sight.

The James Webb Space Telescope (JWST) ushers in a new era of astronomical observation and discovery, offering unprecedented precision in a variety of measurements such as photometry, astrometry, morphology, and shear measurement. Accurate point spread function (PSF) models are crucial for many of these measurements. In this paper, we introduce a hybrid PSF construction method called HybPSF for JWST NIRCam imaging data. HybPSF combines the WebbPSF software, which simulates the PSF for JWST, with observed data to produce more accurate and reliable PSF models. We apply this method to the SMACS J0723 imaging data and construct supplementary structures from residuals obtained by subtracting the WebbPSF PSF model from the data. Our results show that HybPSF significantly reduces discrepancies between the PSF model and the data compared to WebbPSF. Specifically, the PSF shape parameter ellipticity and size comparisons indicate that HybPSF improves precision by a factor of approximately 10 for \$R^2\$ and \$50\%\$ for \$e\$. This improvement has important implications for astronomical measurements using JWST NIRCam imaging data.

The O I 1355{\AA} spectral line is one of the only optically thin lines that are both routinely observed and thought to be formed in the chromosphere. We present analysis of a variety of observations of this line with the Interface Region Imaging Spectrograph (IRIS), and compare it with other IRIS diagnostics as well as diagnostics of the photospheric magnetic field. We utilize special deep exposure modes on IRIS and provide an overview of the statistical properties of this spectral line for several different regions on the Sun. We analyze the spatio-temporal variations of the line intensity, and find that it is often significantly enhanced when and where magnetic flux of opposite polarities cancel. Significant emission occurs in association with chromospheric spicules. Because of the optically thin nature of the O I line, the non-thermal broadening can provide insight into unresolved small-scale motions. We find that the non-thermal broadening is modest, with typical values of 5-10 km/s, and shows some center-to-limb variation, with a modest increase towards the limb. The dependence with height of the intensity and line broadening off-limb is compatible with the line broadening being dominated by the superposition of Alfv\'en waves on different structures. The non-thermal broadening shows a modest but significant enhancement above locations that are in between photospheric magnetic flux concentrations in plage, i.e., where the magnetic field is likely to be more inclined with respect to the line-of-sight. Our measurements provide strict constraints on future theoretical models of the chromosphere.

(Abridged) Very high-energy (VHE, $E>100$ GeV) observations of the blazar Mrk 501 with MAGIC in 2014 have provided evidence for an unusual narrow spectral feature at about 3 TeV during an extreme X-ray flaring activity. The one-zone synchrotron-self Compton scenario, widely used in blazar broadband spectral modeling, fails to explain the narrow TeV component. Motivated by this rare observation, we propose an alternative model for the production of narrow features in the VHE spectra of flaring blazars. These spectral features may result from the decay of neutral pions ($\pi^0$ bumps) that are in turn produced via interactions of protons (of tens of TeV energy) with energetic photons, whose density increases during hard X-ray flares. We explore the conditions needed for the emergence of narrow $\pi^0$ bumps in VHE blazar spectra during X-ray flares reaching synchrotron energies $\sim100$ keV using time-dependent radiative transfer calculations. We focus on high-synchrotron peaked (HSP) blazars, which consist the majority of VHE detected extragalactic sources. We find that synchrotron-dominated flares with peak energies $\gtrsim100$ keV can be ideal periods for the search of $\pi^0$ bumps in the VHE spectra of HSP blazars. Application of the model to the SED of Mrk 501 on MJD 56857.98 shows that the VHE spectrum of the flare is described well by the sum of an SSC component and a distinct $\pi^0$ bump centered at 3 TeV. Spectral fitting of simulated SSC+$\pi^0$ spectra for the Cherenkov Telescope Array (CTA) show that a $\pi^0$ bump could be detected at a 5$\sigma$ significance level with a 30-min exposure.

Identifying the astrophysical sites of the $r$-process, one of the primary mechanisms by which heavy elements are formed, is a key goal of modern astrophysics. The discovery of the brightest gamma-ray burst of all time, GRB 221009A, at a relatively nearby redshift, presented the first opportunity to spectroscopically test the idea that $r$-process elements are produced following the collapse of rapidly rotating massive stars. Here we present spectroscopic and photometric $\textit{James Webb Space Telescope}$ (JWST) observations of GRB 221009A obtained $+168$ and $+170$ rest-frame days after the initial gamma-ray trigger, and demonstrate they are well-described by a supernova (SN) and power-law afterglow, with no evidence for an additional component from $r$-process emission, and that the SN component strongly resembles the near-infrared spectra of previous SNe, including SN 1998bw. We further find that the SN associated with GRB 221009A is slightly fainter than the expected brightness of SN 1998bw at this phase, concluding that the SN is therefore not an unusual GRB-SN. We infer a nickel mass of $\approx0.09$ M$_{\odot}$, consistent with the lack of an obvious SN detection in the early-time data. We find that the host galaxy of GRB 221009A has a very low metallicity of $\approx0.12$ Z$_{\odot}$ and our resolved host spectrum shows that GRB 221009A occurred in a unique environment in its host characterized by strong H$_2$ emission lines consistent with recent star formation, which may hint at environmental factors being responsible for its extreme energetics.

From 2011 to 2021, LAMOST has released a total of 76,167 quasar data. We try to search for gravitationally lensed QSOs by limiting coordinate differences and redshift differences of these QSOs. The name, brightness, spectrum, photometry and other information of each QSO will be visually checked carefully. Special attention should be paid to check whether there are groups of galaxies, gravitationally lensed arcs, Einstein crosses, or Einstein rings near the QSOs. Through careful selection, we select LAMOST J160603.01+290050.8 (A) and LAMOST J160602.81+290048.7 (B) as a candidate and perform an initial analysis. The component A and B are 3.36 arc seconds apart and they display blue during photometric observations. The redshift values of component A and B are 0.2\% different, their Gaia$\_$g values are 1.3\% different, and their ugriz values are 1.0\% or less different. For the spectra covering from 3,690 {\AA} to 9,100 {\AA}, the emission lines of C\,II, Mg, H\,$\gamma$, O\,III, and H\,$\beta$ are present for both component A and B and the ratio of flux(B) to flux(A) from LAMOST is basically a constant, around 2.2. We accidentally find a galaxy group near the component A and B. If the center of dark matter in the galaxy group is at the center between component A and B, the component A and B are probably gravitationally lensed QSOs. We estimate that the Einstein mass is 1.46 $\times$ $10^{11}$ $M_{\odot}$ and the total mass of the lens is 1.34 $\times$ $10^{13}$ $M_{\odot}$. The deflection angle is 1.97 arc seconds at positions A and B and the velocity dispersion is 261\,$km\,s^{-1}$. Theoretically, this candidate could be a pair of fold images of a strong lensing system by a galaxy group, and we will investigate the possibility when the redshifts of nearby galaxies are available.

We present rest-frame UV \textit{Hubble Space Telescope} imaging of the largest and most complete sample of 23 long duration gamma-ray burst (GRB) host galaxies between redshifts 4 and 6. Of these 23, we present new WFC3/F110W imaging for 19 of the hosts, which we combine with archival WFC3/F110W and WFC3/F140W imaging for the remaining four. We use the photometry of the host galaxies from this sample to characterize both the rest-frame UV luminosity function (LF) and the size-luminosity relation of the sample. We find that when assuming the standard Schechter-function parameterization for the UV LF, the GRB host sample is best fit with $\alpha = -1.30^{+0.30}_{-0.25}$ and $M_* = -20.33^{+0.44}_{-0.54}$ mag, which is consistent with results based on $z\sim5$ Lyman-break galaxies. We find that $\sim68\%$ of our size-luminosity measurements fall within or below the same relation for Lyman-break galaxies at $z\sim4$. This study observationally confirms expectations that at $z\sim5$ Lyman-break and GRB host galaxies should trace the same population and demonstrates the utility of GRBs as probes of hidden star-formation in the high-redshift universe. Under the assumption that GRBs unbiasedly trace star formation at this redshift, our non-detection fraction of 7/23 is consistent at the $95\%$-confidence level with $13 - 53\%$ of star formation at redshift $z\sim5$ occurring in galaxies fainter than our detection limit of $M_{1600 A} \sim -18.3$ mag.

In the smooth mass distribution model, the critical curve represents a line with magnification divergence on the image plane in a strong gravitational lensing system. Considering the microlensing effects caused by discrete masses, the magnification map in the source plane exhibits a complex structure, which offers a promising way for detecting dark matter. However, simulating microlensing near the critical curve poses challenges due to magnification divergence and the substantial computational demands involved. To achieve the required simulation accuracy, direct inverse ray-shooting would require significant computational resources. Therefore we applied a GPU-based code optimized with interpolation method to enable efficient computation on a large scale. Using the GPU of NVIDIA Tesla V100S PCIe 32GB, it takes approximately 7000 seconds to calculate the effects of around 13,000 microlenses for a simulation involving 1013 emitted rays. Then we generated 80 magnification maps, and select 800 light curves for a statistical analysis of microcaustic density and peak magnification.

Precursors provide important clues to the nature of gamma-ray burst (GRB) central engines and can be used to contain GRB physical processes. In this letter, we study the self-organized criticality in precursors of long GRBs in the third Swift/BAT Catalog. We investigate the differential and cumulative size distributions of 100 precursors, including peak flux, duration, rise time, decay time, and quiescent time with the Markov Chain Monte Carlo technique. It is found that all of the distributions can be well described by power-law models and understood within the physical framework of a self-organized criticality system. In addition, we inspect the cumulative distribution functions of the size differences with a q-Gaussian function. The scale-invariance structures of precursors further strengthen our findings. Particularly, similar analyses are made in 127 main bursts. The results show that both precursors and main bursts can be attributed to an self-organized criticality system with the spatial dimension S = 3 and driven by the similar magnetically dominated process.

The division of Gamma-ray bursts (GRBs) into different classes, other than the "short" and "long", has been an active field of research. We investigate whether GRBs can be classified based on a broader set of parameters, including prompt and plateau emission ones. Observational evidence suggests the existence of more GRB sub-classes, but results so far are either conflicting or not statistically significant. The novelty here is producing a machine-learning-based classification of GRBs using their observed X-rays and optical properties. We used two data samples: the first, composed of 203 GRBs, is from the Neil Gehrels Swift Observatory (Swift/XRT), and the latter, composed of 134 GRBs, is from the ground-based Telescopes and Swift/UVOT. Both samples possess the plateau emission (a flat part of the light curve happening after the prompt emission, the main GRB event). We have applied the Gaussian Mixture Model (GMM) to explore multiple parameter spaces and sub-class combinations to reveal if there is a match between the current observational sub-classes and the statistical classification. With these samples and the algorithm, we spot a few micro-trends in certain cases, but we cannot conclude that any clear trend exists in classifying GRBs. These microtrends could point towards a deeper understanding of the physical meaning of these classes (e.g., a different environment of the same progenitor or different progenitors). However, a larger sample and different algorithms could achieve such goals. Thus, this methodology can lead to deeper insights in the future.

FIRST is a post Extreme Adaptive-Optics (ExAO) spectro-interferometer operating in the Visible (600-800 nm, R~400). Its exquisite angular resolution (a sensitivity analysis of on-sky data shows that bright companions can be detected down to 0.25lambda/D) combined with its sensitivity to pupil phase discontinuities (from a few nm up to dozens of microns) makes FIRST an ideal self-calibrated solution for enabling exoplanet detection and characterization in the future. We present the latest on-sky results along with recent upgrades, including the integration and on-sky test of a new spectrograph (R~3,600) optimized for the detection of H-alpha emission from young exoplanets accreting matter.

Based on the APOGEE survey we conducted a search for carbon-deficient red giants (CDGs). We found 103 new CDGs, increasing the number in the literature by more than a factor of 3. CDGs are very rare, representing $0.03$~per cent of giants. They appear as an extended tail off the normal carbon distribution. We show that they are found in all components of the Galaxy, contrary to previous findings. The location of CDGs in the Hertzsprung-Russell diagram (HRD) shows that they are primarily intermediate-mass stars ($2-4~\rm{M}_{\odot}$). Their extended distribution may indicate that CDGs can also sometimes have $M < 2.0~\rm{M}_{\odot}$. We attempted to identify the evolutionary phases of the CDGs using stellar model tracks. We found that the bulk of the CDGs are likely in the subgiant branch or red clump phase, whereas other CDGs may be in the red giant branch or early asymptotic giant branch phases. Degeneracy in the HRD makes exact identification difficult. We examined their C, N, and O compositions and confirmed previous studies showing that the envelope material has undergone extensive hydrogen burning through the CN(O) cycle. The new-CDGs have [C+N+O/Fe] that generally sum to zero, indicating that they started with scaled-solar composition. However, the previously known-CDGs generally have [C+N+O/Fe$] > 0.0$, indicating that some He-burning products were added to their envelopes. As to the site(s) in which this originally occurred, we do not find a convincing solution.

The full phase space information on the kinematics of a huge number of stars provided by the Gaia third Data Release raises the demand for a better understanding of the 3D stellar dynamics. In this paper, we investigate the possible regimes of motion of stars in the axisymmetric approximation of a Galactic potential model. The model consists of three components: the axisymmetric disk, the central spheroidal bulge and the spherical halo of dark matter. The axisymmetric disk is divided into stellar and gaseous disk subcomponents, each one modeled by three Miyamoto-Nagai profiles. The physical and structural parameters of the Galaxy components are adjusted by observational kinematic constraints. The phase space of the two-degrees-of-freedom model is studied by means of the Poincar\'e and dynamical mapping, the dynamical spectrum method and the direct numerical integrations of the Hamiltonian equations of motion. For the chosen physical parameters, the nearly-circular and low-altitude stellar behaviour is composed of two weakly coupled simple oscillations, radial and vertical motions. The amplitudes of the vertical oscillations of these orbits are gradually increasing with the growing Galactocentric distances, in concordance with the exponential mass decay assumed. However, for increasing planar eccentricities and the altitudes over the equatorial disk, new regimes of stellar motion emerge as a result of the beating between the radial and vertical oscillation frequencies, which we refer to as e-z resonances. The corresponding resonant motion produces characteristic sudden increase or decrease of the amplitude of the vertical oscillation, bifurcations in the dynamical spectra and the chains of islands of stable motion in the phase space. The results obtained can be useful in the understanding and interpretation of the features observed in the stellar 3D distribution around the Sun.

Recently, some fast radio bursts (FRBs) have been reported to exhibit complex and diverse variations in Faraday rotation measurements (RM) and polarization, suggesting that dynamically evolving magnetization environments may surround them. In this paper, we investigate the Faraday conversion (FC) effect in a binary system involving an FRB source and analyze the polarization evolution of FRBs. For an strongly magnetized high-mass companion binary (HMCB), when an FRB with $\sim100\%$ linear polarization passes through the radial magnetic field of the companion star, the circular polarization (CP) component will be induced and oscillate symmetrically around the point with the CP degree equal to zero, the rate and amplitude of the oscillation decrease as the frequency increases. The very strong plasma column density in the HMCBs can cause CP to oscillate with frequency at a very drastic rate, which may lead to depolarization. Near the superior conjunction of the binary orbit, the DM varies significantly due to the dense plasma near the companion, and the significant FC also occurs in this region. As the pulsar moves away from the superior conjunction, the CP gradually tends towards zero and then returns to its value before incidence. We also investigate the effect of the rotation of the companion star. We find that a sufficiently significant RM reversal can be produced at large magnetic inclinations and the RM variation is very diverse. Finally, we apply this model to explain some polarization observations of PSR B1744-24A and FRB 20201124A.

We investigate the orbital evolution of planetesimals in the inner disk in the presence of nebula gas and a (proto-) cold Jupiter. By varying the mass, eccentricity, and semi-major axis of the planet, we study the dependence of the relative velocities of the planetesimals on these parameters. For classic small planetesimals ($10^{16}-10^{20} $g) whose mutual gravitational interaction is negligible, gas drag introduces a size-dependent alignment of orbits and keeps the relative velocity low for similar-size bodies, while preventing orbital alignment for different-size planetesimals. Regardless of the location and the mass ratio of the planetesimals, increasing the mass and eccentricity or decreasing the orbital distance of the planet always leads to higher relative velocities of planetesimals. However, for massive planetesimals, the interplay of viscous stirring, gas damping, and secular perturbation results in lower velocity dispersion of equal-size planetesimals when the planet is more massive or when it is located on a closer or more eccentric orbit. The random velocities of such planetesimals remain almost unperturbed when the planet is located beyond Jupiter's current orbit, or when it is less massive or less eccentric than Jupiter. Unlike small planetesimals, such large planetesimals can grow in a runaway fashion as in the unperturbed case. Our results imply that the presence of a cold Jupiter does not impede the formation of inner rocky planets through planetesimal accretion, provided that the planetesimals are initially large.

Fast Radio Bursts (FRBs) are extragalactic transients of (sub-)millisecond duration that show wide-ranging spectral, temporal, and polarimetric properties. The polarimetric analysis of FRBs can be used to probe intervening media, study the emission mechanism, and test possible progenitor models. In particular, low frequency depolarisation of FRBs can identify dense, turbulent, magnetised, ionised plasma thought to be near the FRB progenitor. An ensemble of repeating FRBs has shown low-frequency depolarisation. The depolarisation is quantified by the parameter $\sigma_{\rm RM}$, which correlates with proxies for both the turbulence and mean magnetic field strength of the putative plasma. However, while many non-repeating FRBs show comparable scattering (and hence inferred turbulence) to repeating FRBs, it is unclear whether their surrounding environments are comparable to those of repeating FRBs. To test this, we analyse the spectro-polarimetric properties of five one-off FRBs and one repeating FRB, detected and localised by the Australian Square Kilometer Array Pathfinder. We search for evidence of depolarisation due to $\sigma_{\rm RM}$ and consider models where the depolarisation is intrinsic to the source. We find no evidence (for or against) the sample showing spectral depolarisation. Under the assumption that FRBs have multipath propagation-induced depolarisation, the correlation between our constraint on $\sigma_{\rm RM}$ and RM is consistent with repeating FRBs only if the values of $\sigma_{\rm RM}$ are much smaller than our upper limits. The observations provide further evidence for differences in the environments and sources of one-off and repeating FRBs.

Shear instabilities can be the source of significant amounts of turbulent mixing in stellar radiative zones. Past attempts at modeling their effects (either theoretically or using numerical simulations) have focused on idealized geometries where the shear is either purely vertical or purely horizontal. In stars, however, the shear can have arbitrary directions with respect to gravity. In this work, we use direct numerical simulations to investigate the nonlinear saturation of shear instabilities in a stably stratified fluid, where the shear is sinusoidal in the horizontal direction, and either constant or sinusoidal in the vertical direction. We find that, in the parameter regime studied here (non-diffusive, fully turbulent flow), the mean vertical shear does not play any role in controlling the dynamics of the resulting turbulence unless its Richardson number is smaller than one (approximately). As most stellar radiative regions have a Richardson number much greater than one, our result implies that the vertical shear can essentially be ignored in the computation of the vertical mixing coefficient associated with shear instabilities for the purpose of stellar evolution calculations, even when it is much larger than the horizontal shear (as in the solar tachocline, for instance).

Many experimental quantities show a power-law distribution $p(x)\propto x^{-\alpha}$. In astrophysics, examples are: size distribution of dust grains or luminosity function of galaxies. Such distributions are characterized by the exponent $\alpha$ and by the extremes $x_\text{min}$ $x_\text{max}$ where the distribution extends. There are no mathematical tools that derive the three unknowns at the same time. In general, one estimates a set of $\alpha$ corresponding to different guesses of $x_\text{min}$ $x_\text{max}$. Then, the best set of values describing the observed data is selected a posteriori. In this paper, we present a tool that finds contextually the three parameters based on simple assumptions on how the observed values $x_i$ populate the unknown range between $x_\text{min}$ and $x_\text{max}$ for a given $\alpha$. Our tool, freely downloadable, finds the best values through a non-linear least-squares fit. We compare our technique with the maximum likelihood estimators for power-law distributions, both truncated and not. Through simulated data, we show for each method the reliability of the computed parameters as a function of the number $N$ of data in the sample. We then apply our method to observed data to derive: i) the slope of the core mass function in the Perseus star-forming region, finding two power-law distributions: $\alpha=2.576$ between $1.06\,M_{\sun}$ and $3.35\,M_{\sun}$, $\alpha=3.39$ between $3.48\,M_{\sun}$ and $33.4\,M_{\sun}$; ii) the slope of the $\gamma$-ray spectrum of the blazar J0011.4+0057, extracted from the Fermi-LAT archive. For the latter case, we derive $\alpha=2.89$ between 1,484~MeV and 28.7~GeV; then we derive the time-resolved slopes using subsets 200 photons each.

White dwarfs (WDs) in active galactic nucleus (AGNs) discs might migrate to the inner radii of the discs and form restricted three-body systems with two WDs moving around the central supermassive black hole (SMBH) in close orbits. These systems could be dynamical unstable, which can lead to very close encounters or direct collisions. In this work, we use N-body simulations to study the evolution of such systems with the different initial orbital separation $p$, relative orbital inclination $\Delta{i}$ and SMBH mass $M$. It is found that the close encounters of WDs mainly occur at $1.1R_{\rm H} \lesssim p \lesssim 2\sqrt{3}R_{\rm H}$, where $R_{\rm H}$ is the mutual Hill radius. For $p<1.1R_{\rm H}$, the majority of WDs move in horseshoe or tadpole orbits, and only few of them with small initial orbital phase difference undergo close encounters. For $p=3.0R_{\rm H}$, WD-WD collisions occur in most of the samples within a time of $10^5P_1$, and considerable collisions occur within a time of $t<62P_1$ for small orbital radii, where $P_1$ is the orbital period. The peak of the closest separation distribution increase and the WD-WD collision fraction decreases with an increase of the relative inclination. The closest separation distribution is similar in cases with the different SMBH mass, but the WD-WD collision fraction decreases as the mass of SMBHs increases. According to our estimation, the event rate of the cosmic WD-WD collision in AGN discs is about $300{\rm Gpc^{-3}yr^{-1}}$, roughly $1\%$ of the one of the observed type Ia supernova. The corresponding electromagnetic emission signals can be observed by large surveys of AGNs.

Glycolaldehyde (CH$_{2}$OHCHO) is the simplest monosaccharide sugar in the interstellar medium, and it is directly involved in the origin of life via the 'RNA world' hypothesis. We present the first detection of glycolaldehyde (CH$_{2}$OHCHO) towards the hot molecular core G358.93-0.03 MM1 using the Atacama Large Millimeter/Submillimeter Array (ALMA). The calculated column density of CH$_{2}$OHCHO towards G358.93-0.03 MM1 is (1.52$\pm$0.9)$\times$10$^{16}$ cm$^{-2}$ with an excitation temperature of 300$\pm$68.5 K. The derived fractional abundance of CH$_{2}$OHCHO with respect to H$_{2}$ is (4.90$\pm$2.92)$\times$10$^{-9}$, which is consistent with that estimated by existing two-phase warm-up chemical models. We discuss the possible formation pathways of CH$_{2}$OHCHO within the context of hot molecular cores and hot corinos and find that CH$_{2}$OHCHO is likely formed via the reactions of radical HCO and radical CH$_{2}$OH on the grain surface of G358.93-0.03 MM1.

In the inner region of the disc of an active galactic nucleus (AGN), the collision of two white dwarfs (WDs) through Jacobi capture might be inevitable, leading to a Type Ia supernova (SN Ia) explosion. This transient event, influenced by the disc gas and the gravity of the supermassive black hole (SMBH), exhibits distinct characteristics compared to normal SNe Ia. The energy of the explosion is mainly stored in the ejecta in the form of kinetic energy. Typically, the ejecta is not effectively decelerated by the AGN disc and rushes rapidly out of the AGN disc. However, under the influence of the SMBH, most of the ejecta falls back toward the AGN disc. As the fallback ejecta becomes more dispersed, it interacts with the disc gas, converting its kinetic energy into thermal energy. This results in a high-energy transient characterized by a rapid initial rise followed by a decay with $L\propto t^{-2.8}$. The time-scale of the transient ranges from hours to weeks, depending on the mass of the SMBH. This process generates high-energy radiation spanning from hard X-rays to the soft $\gamma$ range. Additionally, the subsequent damage to the disc may result in changing-look AGNs. Moreover, the falling back of SNe Ia ejecta onto the AGN disc significantly increases the metallicity of the AGN and can even generate heavy elements within the AGN discs.

Cosmic string loops are non-linear density fluctuations which form in the early universe and could play an important role in explaining many phenomena which are in tension with the standard $\Lambda$CDM model. Hence, the details of the accretion process onto cosmic string loops should be studied in detail. Most previous works view loops as point masses and ignore the impact of a finite loop size. In this work, we utilize the Zel'dovich approximation to calculate the non-linear mass sourced by a static extended loop with a time-averaged density profile derived from the trajectory of the loop oscillation, and compare the result with what is obtained for a point-mass source. We find that the finite size of a loop mainly affects the evolution of turnaround shells during the early stages of accretion, converging to the point mass result after a critical redshift, $z^{(II)/(III)}_{c}$. For $z>z^{(II)/(III)}_{c}$, the total accreted mass surrounding a loop is suppressed relative to the point mass case and has a growth rate proportional to $(1+z)^{-3/2}$. As an immediate extension, we also qualitatively analyse the accretion onto moving point masses and onto moving extended loops. In addition to the reduction in the nonlinear mass, the loop finite size also changes the shape of the turnaround surface at early stages of accretion.

We present a study aiming at detecting VMS in local star-forming region from the imprint they leave on the integrated UV and optical light. We analyzed a sample of 27 star-forming regions and galaxies in the local Universe. We selected sources with a metallicity close to that of the LMC. We defined empirical criteria to distinguish sources dominated by VMS and Wolf-Rayet stars (WR), using template spectra of VMS- and WR-dominated regions. We subsequently built population synthesis models with an updated treatment of VMS. We show that the UV range alone is not sufficient to distinguish between VMS- and WR-dominated sources. The region of the WR bumps in the optical breaks the degeneracy. In particular, the morphology of the blue bump at 4640-4686 A is a key diagnostic. Beyond the prototypical R136 region we identify two galaxies showing clear signatures of VMS. In two other galaxies or regions the presence of VMS can be suspected, as already discussed in the literature. The stellar population is clearly dominated by WR stars in seven other sources. The most recent BPASS population synthesis models can neither account for the strong HeII 1640 emission, nor for the shape of the blue bump in VMS- and WR-dominated sources. Our models that include VMS more realistically reproduce the UV-optical spectra of VMS-dominated sources. We conclude that VMS are present in some local star-forming regions, but that separating them from WR-dominated populations requires optical spectroscopy with a high signal-to-noise ratio. A high equivalent width of HeII 1640 is not a sufficient condition for identifying VMS. Populations synthesis models need to take VMS into account by incorporating not only evolutionary tracks, but also dedicated spectral libraries. Finally, we stress that the treatment of WR stars needs to be improved as well.

Atmospheric characterisation of exoplanets from the ground is an actively growing field of research. In this context we have created the ATMOSPHERIX consortium: a research project aimed at characterizing exoplanets atmospheres using ground-based high resolution spectroscopy. This paper presents the publicly-available data analysis pipeline and demonstrates the robustness of the recovered planetary parameters from synthetic data. Simulating planetary transits using synthetic transmission spectra of a hot Jupiter that were injected into real SPIRou observations of the non-transiting system Gl 15 A, we show that our pipeline is successful at recovering the planetary signal and input atmospheric parameters. We also introduce a deep learning algorithm to optimise data reduction which proves to be a reliable, alternative tool to the commonly used principal component analysis. We estimate the level of uncertainties and possible biases when retrieving parameters such as temperature and composition and hence the level of confidence in the case of retrieval from real data. Finally, we apply our pipeline onto two real transits of HD~189733 b observed with SPIRou and obtain similar results than in the literature. In summary, we have developed a publicly available and robust pipeline for the forthcoming studies of the targets to be observed in the framework of the ATMOSPHERIX consortium, which can easily be adapted to other high resolution instruments than SPIRou (e.g. VLT-CRIRES, MAROON-X, ELT-ANDES)

In a companion paper, we introduced a publicly-available pipeline to characterise exoplanet atmospheres through high-resolution spectroscopy. In this paper, we use this pipeline to study the biases and degeneracies that arise in atmospheric characterisation of exoplanets in near-infrared ground-based transmission spectroscopy. We inject synthetic planetary transits into sequences of SPIRou spectra of the well known M dwarf star Gl 15 A, and study the effects of different assumptions on the retrieval. We focus on (i) mass and radius uncertainties, (ii) non isothermal vertical profiles and (iii) identification and retrieval of multiple species. We show that the uncertainties on mass and radius should be accounted for in retrievals and that depth-dependent temperature information can be derived from high-resolution transmission spectroscopy data. Finally, we discuss the impact of selecting wavelength orders in the retrieval and the issues that arise when trying to identify a single species in a multi-species atmospheric model. This analysis allows us to understand better the results obtained through transmission spectroscopy and their limitations in preparation to the analysis of actual SPIRou data.

Recent Event Horizon Telescope observations of M87* and Sgr A* strongly suggests the presence of supermassive black hole at their respective cores. In this work, we use the semi-analytic Radiatively Inefficient Accretion Flows (RIAF) model to investigate the resulting images of Joshi-Malafarina-Narayan (JMN-1) naked singularity and the Schwarzschild BH. We aim at choosing the JMN-1 naked singularity model and compare the synchrotron images with the Schwarzschild solution to search any distinct features which can distinguish the two objects and find alternative to the black hole solution. We perform general relativistic ray-tracing and radiative transfer simulations using Brahma code to generate synchrotron emission images utilising thermal distribution function for emissivity and absorptivity. We investigate effects in the images by varying inclination angle, disk width and frequency. The shadow images simulated by the JMN-1 model closely resemble those generated by the Schwarzschild black hole. When we compare these images, we find that the disparities between them are minimal. We conduct simulations using various plasma parameters, but the resulting images remain largely consistent for both scenarios. This similarity is evident in the horizontal cross-sectional brightness profiles of the two instances. Notably, the JMN-1 model exhibits slightly higher intensity in comparison to the Schwarzschild black hole. We conclude that JMN-1 presents itself as a viable substitute for the black hole scenario. This conclusion is not solely grounded in the fact that they are indistinguishable from their respective shadow observations, but also in the consideration that JMN-1 emerges as an end state of a continual gravitational collapse. This paradigm not only allows for constraints on spacetime but also provides a good probe for the nature of the central compact object.

Photometric observations from the last decade have revealed additional low-amplitude periodicities in many classical pulsators that are likely due to pulsations in non-radial modes. One group of multi-mode RR Lyrae stars, the so-called 0.61 stars, is particularly interesting. In these stars, the radial first overtone is accompanied by additional signals with period ratios around 0.61. The most promising explanation for these signals is pulsation in non-radial modes of degrees 8 and 9. If the theory behind the additional signals in the 0.61 stars is substantiated, it would allow us to use non-radial modes to study classical pulsators. We aim to perform asteroseismic modeling of selected 0.61 stars with independently determined physical parameters to test whether this assumption behind the modeling leads to correct results. Namely, we test whether the additional signals are indeed due to non-radial modes of the proposed moderate degrees. We selected a number of and RR Lyrae stars that are also 0.61 stars and have good observational constraints on their other physical parameters. We assume that the nature of those modes is correctly explained with non-radial modes of degrees 8 or 9. Using this assumption and observational constraints on physical parameters, we performed asteroseismic modeling to test whether the observed periods and period ratios can be reproduced. For the majority of selected targets, we obtained a good match between observed and calculated periods and period ratios. For a few targets however, the results obtained are ambiguous and not straightforward to interpret.

Bolometric light curves play an important role in understanding the underlying physics of various astrophysical phenomena, as they allow for a comprehensive modeling of the event and enable comparison between different objects. However, constructing these curves often requires the approximation and extrapolation from multicolor photometric observations. In this study, we introduce vector Gaussian processes as a new method for reduction of supernova light curves. This method enables us to approximate vector functions, even with inhomogeneous time-series data, while considering the correlation between light curves in different passbands. We applied this methodology to a sample of 29 superluminous supernovae (SLSNe) assembled using the Open Supernova Catalog. Their multicolor light curves were approximated using vector Gaussian processes. Subsequently, under the black-body assumption for the SLSN spectra at each moment of time, we reconstructed the bolometric light curves. The vector Gaussian processes developed in this work are accessible via the Python library gp-multistate-kernel on GitHub. Our approach provides an efficient tool for analyzing light curve data, opening new possibilities for astrophysical research.

Context. V1309 Sco is an example of a red nova, a product of the merger between non-compact stars. V1309 Sco is particularly important within the class of red novae due to the abundance of the progenitor binary before the merger. Aims. We aim to investigate the spatio-kinematic and chemical properties of the circumstellar environment, including deriving the physical conditions and establishing the origins of the different circumstellar components. Methods. We use radioactive transfer modelling of molecular emission in sub-mm spectra to examine the properties of the molecular gas, and use forbidden line diagnostics from optical spectra to constrain electron density and temperature using forbidden line diagnostics. We compare line intensities from shock models to observations to look for and constrain shocks. Results. We derive a new kinematical distance of 5.6 kpc to the source. The detection of ro-vibrational H2 and sub-mm HCO+ emission in 2016 and 2019, respectively, indicate active shock interactions within the circumstellar environment. The velocity profiles of both H2 and HCO+, as well as the moment-1 maps of sub-mm CO and 29-SiO, indicate a bipolar structure that may be asymmetric. The sub-mm and optical molecular emission exhibits temperatures of 35-113 and 200 K, respectively, whilst the atomic gas is much hotter, with temperatures of 5-15 kK, which may be due to shock heating. Conclusions. The detection of a bipolar structure in V1309 Sco indicates further similarities with the structure of another Galactic red nova, V4332 Sgr. It provides evidence that bipolar structures may be common in red novae. All collected data are consistent with V1309 Sco bring a kinematically and chemically complex system.

Accretion plays a central role in the physics that governs the evolution and dispersal of protoplanetary disks. The primary goal of this paper is to analyze the stability over time of the mass accretion rate onto TW Hya, the nearest accreting solar-mass young star. We measure veiling across the optical spectrum in 1169 archival high-resolution spectra of TW Hya, obtained from 1998--2022. The veiling is then converted to accretion rate using 26 flux-calibrated spectra that cover the Balmer jump. The accretion rate measured from the excess continuum has an average of $2.51\times10^{-9}$~M$_\odot$~yr$^{-1}$ and a Gaussian distribution with a FWHM of 0.22 dex. This accretion rate may be underestimated by a factor of up to 1.5 because of uncertainty in the bolometric correction and another factor of 1.7 because of excluding the fraction of accretion energy that escapes in lines, especially Ly$\alpha$. The accretion luminosities are well correlated with He line luminosities but poorly correlated with H$\alpha$ and H$\beta$ luminosity. The accretion rate is always flickering over hours but on longer timescales has been stable over 25 years. This level of variability is consistent with previous measurements for most, but not all, accreting young stars.

In radio astronomy, visibility data, which are measurements of wave signals from radio telescopes, are transformed into images for observation of distant celestial objects. However, these resultant images usually contain both real sources and artifacts, due to signal sparsity and other factors. One way to obtain cleaner images is to reconstruct samples into dense forms before imaging. Unfortunately, existing visibility reconstruction methods may miss some components of the frequency data, so blurred object edges and persistent artifacts remain in the images. Furthermore, the computation overhead is high on irregular visibility samples due to the data skew. To address these problems, we propose PolarRec, a reconstruction method for interferometric visibility data, which consists of a transformer-conditioned neural fields pipeline with a polar coordinate representation. This representation matches the way in which telescopes observe a celestial area as the Earth rotates. We further propose Radial Frequency Loss function, using radial coordinates in the polar coordinate system to correlate with the frequency information, to help reconstruct complete visibility. We also group visibility sample points by angular coordinates in the polar coordinate system, and use groups as the granularity for subsequent encoding with a Transformer encoder. Consequently, our method can capture the inherent characteristics of visibility data effectively and efficiently. Our experiments demonstrate that PolarRec markedly improves imaging results by faithfully reconstructing all frequency components in the visibility domain while significantly reducing the computation cost.

Fuzzy Dark Matter (FDM) has recently gained attention as a motivated candidate for the dark matter (DM) content of the Universe, as opposed to the commonly assumed cold DM (CDM), since the soliton profile intrinsic to FDM models was found to be particularly well suited to reproduce observed galaxy mass profiles. While FDM as a single DM component has been strongly constrained by multiple probes, there remained a mass window between $10^{-25}\,\mathrm{eV}$ and $10^{-24}\,\mathrm{eV}$ in which it can comprise a large portion ($\gtrsim\mathcal{O}(10\%)$) of the total DM. In this work, we consider gravitational lensing measurements in the strong lensing regime, which are one of the only means to directly constrain the distribution and profile of DM in astronomical bodies. Using a simple model that combines a soliton FDM component with a Navarro-Frenk-White (NFW) profile, we explore under what conditions DM halos with this hybrid profile are able to reproduce the observed Einstein radii of several known lenses. We find that FDM with a particle mass of $\lesssim 10^{-24}\,{\rm eV}$ cannot explain the observations if it makes up more than $\sim10\%$ of the total DM, effectively closing the lingering FDM mass window.

Crater mapping using neural networks and other automated methods has increased recently with automated Crater Detection Algorithms (CDAs) applied to planetary bodies throughout the solar system. A recent publication by Benedix et al. (2020) showed high performance at small scales compared to similar automated CDAs but with a net positive diameter bias in many crater candidates. I compare the publicly available catalogs from Benedix et al. (2020) and Lee & Hogan (2021) and show that the reported performance is sensitive to the metrics used to test the catalogs. I show how the more permissive comparison methods indicate a higher CDA performance by allowing worse candidate craters to match ground-truth craters. I show that the Benedix et al. (2020) catalog has a substantial performance loss with increasing latitude and identify an image projection issue that might cause this loss. Finally, I suggest future applications of neural networks in generating large scientific datasets be validated using secondary networks with independent data sources or training methods.

Magnetic flux ropes are a key component of coronal mass ejections, forming the core of these eruptive phenomena. However, determining whether a flux rope is present prior to eruption onset and, if so, the rope's handedness and the number of turns that any helical field lines make is difficult without magnetic field modelling or in-situ detection of the flux rope. We present two distinct observations of plasma flows along a filament channel on 4 and 5 September 2022 made using the \textit{Solar Orbiter} spacecraft. Each plasma flow exhibited helical motions in a right-handed sense as the plasma moved from the source active region across the solar disk to the quiet Sun, suggesting that the magnetic configuration of the filament channel contains a flux rope with positive chirality and at least one turn. The length and velocity of the plasma flow increased from the first to the second observation, suggesting evolution of the flux rope, with the flux rope subsequently erupting within $\sim$5~hours of the second plasma flow. The erupting flux rope then passed over the \textit{Parker Solar Probe} spacecraft during its Encounter 13, enabling \textit{in-situ} diagnostics of the structure. Although complex and consistent with the flux rope erupting from underneath the heliospheric current sheet, the \textit{in-situ} measurements support the inference of a right-handed flux rope from remote-sensing observations. These observations provide a unique insight into the eruption and evolution of a magnetic flux rope near the Sun.

Besides the neutrino source detected by IceCube, NGC 1068, the association of the IceCube-170922A neutrino with the blazar in a flaring state among several wavelengths (from radio up to high-energy (HE) gamma-rays), the site and mechanisms of production of HE neutrino remains in discussion. Extragalactic sources such as Quasars, Blazars, Radio galaxies, and Gamma-ray bursts have been proposed as progenitors of HE neutrinos. In this work, we study the Blazars reported by Fermi-LAT in the 4LAC catalog, which are embedded inside the 90\% error of the best-fit position from the neutrinos reported by IceCube. We propose a one-zone lepto-hadronic scenario to describe the broadband Spectral Energy Distribution and then estimate the number of neutrinos to compare with those in the direction of each source. A brief discussion is provided of the results.

We investigate whether the oblate, spheroidal morphology of common dwarf spheroidal galaxies (dSph) may result from the slow relaxation of stellar orbits within a halo of Wave Dark Matter ($\psi$DM) when starting from an initial disk of stars. Stellar orbits randomly walk over a Hubble time, perturbed by the pervasive "granular" interference pattern of $\psi$DM, that fully modulates the dark matter density on the de Broglie scale. Our simulations quantify the level of stellar disk thickening over the Hubble time, showing that distribution of stars is predicted to become an oblate spheroid of increasing radius, that plausibly accounts for the morphology of dSph galaxies. We predict a low level of residual rotation remains after a Hubble time at the 1-3 km/s level, depending on orientation, that compares with recent claims of rotation for some well studied local dSph galaxies. This steady internal dynamical evolution may be witnessed directly with JWST for well resolved dwarf galaxies, appearing more oblate with look back time and tending to small disks of young stars at high redshift.

The data recently released by the North American Nanohertz Observatory for Gravitational Waves provides compelling evidence supporting the existence of a stochastic signal that aligns with a gravitational-wave background. We show this signal can be the scalar-induced gravitational waves from the Higgs inflation model with the parametric amplification mechanism. Such a gravitational-wave background naturally predicts the substantial existence of planet-mass primordial black holes, which can be planet 9 in our solar system and the lensing objects for the ultrashort-timescale microlensing events observed by the Optical Gravitational Lensing Experiment. The future observations of stochastic gravitational wave background by pulsar timing arrays and planet-mass primordial black holes provide such a possibility to give further confirmations on Higgs inflation, which unifies two fundamental aspects of theoretical physics, particle physics and cosmology.

We present the first study on the gravitational impact of supernova feedback in an isolated soliton and a spherically symmetric dwarf SFDM halo of virial mass $1\times 10^{10}\mathrm{M_\odot}$. We use a boson mass $m=10^{-22}\mathrm{eV/c^2}$ and a soliton core $r_c \approx 0.7$kpc, comparable to typical half-light radii of Local Group dwarf galaxies. We simulate the rapid gas removal from the center of the soliton by a concentric external time-dependent Hernquist potential. We explore two scenarios of feedback blowouts: i) a massive single burst, and ii) multiple consecutive blowouts injecting the same total energy to the system, including various magnitudes for the blowouts in both scenarios. In all cases, we find one single blowout has a stronger effect on reducing the soliton central density. Feedback leads to central soliton densities that oscillate quasi-periodically for an isolated soliton and stochastically for a SFDM halo. The range in the density amplitude depends on the strength of the blowout, however we observe typical variations of a factor of $\geqslant$2. One important consequence of the stochastic fluctuating densities is that, if we had no prior knowledge of the system evolution, we can only know the configuration profile at a specific time within some accuracy. By fitting soliton profiles at different times to our simulated structures, we found the (1-$\sigma$) scatter of their time-dependent density profiles. For configurations within the 1$\sigma$ range, we find the inferred boson mass is typically less than 20\% different from the real value used in our simulations. Finally, we compare the observed dynamical masses of field dwarf galaxies in our Local Group with the implied range of viable solitons from our simulations and find good agreement.

We analyze new JWST NIRCam and NIRSpec data on the redshift 9.11 galaxy MACS1149-JD1. Our NIRCam imaging data reveal that JD1 comprises three spatially distinct components. Our spectroscopic data indicate that JD1 appears dust-free but is already enriched, $12 + \log {\rm (O/H) } = 7.875^{+0.042}_{-0.045}$. We also find that the Carbon and Neon abundances in JD1 are below the solar abundance ratio. Particularly the Carbon under-abundance is suggestive of recent star formation where Type~II supernovae have already enriched the ISM in Oxygen but intermediate mass stars have not yet enriched the ISM in Carbon. A recent burst of star formation is also revealed by the star formation history derived from NIRCam photometry. Our data do not reveal the presence of a significant amount of old populations, resulting in a factor of $\sim7\times$ smaller stellar mass than previous estimates. Thus, our data support the view that JD1 is a young object.

As we enter the era of TESS and JWST, instrumentation that can carry out radial velocity measurements of exoplanet systems is in high demand. We will address this demand by upgrading the UC Lick Observatory's 2.4-meter Automated Planet Finder (APF) telescope with an adaptive optics (AO) system. The AO upgrade will be directly integrated into the APF telescope by replacing the telescope's static secondary mirror with a 61-actuator adaptive secondary mirror (ASM) to minimize the disturbance to the spectrograph optics. This upgrade is enabled by The Netherlands Organization for Applied Scientific Research's (TNO) large-format deformable mirror technology, which will be constructed using a new style of high-efficiency hybrid-variable reluctance actuator. We outline the technical design and manufacturing plan for the proposed APF AO upgrade and simulate the improvement to the science yield using HCIpy. Our simulations predict the AO upgrade will reduce the PSF instabilities due to atmospheric turbulence, concentrating the light on the spectrograph slit by a multiplicative factor of more than two (doubling the telescope's observing efficiency) for targets as dim as I = 14. When completed, the APF adaptive secondary mirror will be among the first pairings of an ASM with a radial velocity spectrograph and become a pathfinder for similar AO systems in telescopes of all sizes.

At $z \lesssim 1$, shock heating caused by large-scale velocity flows and possibly violent feedback from galaxy formation, converts a significant fraction of the cool gas ($T\sim 10^4$ K) in the intergalactic medium (IGM) into warm-hot phase (WHIM) with $T >10^5$K, resulting in a significant deviation from the previously tight power-law IGM temperature-density relationship, $T=T_0 (\rho / {\bar{\rho}})^{\gamma -1}$. This study explores the impact of the WHIM on measurements of the low-$z$ IGM thermal state, $[T_0,\gamma]$, based on the $b$-$N_{H I}$ distribution of the Lyman-$\alpha$ forest. Exploiting a machine learning-enabled simulation-based inference method trained on Nyx hydrodynamical simulations, we demonstrate that [$T_0$, $\gamma$] can still be reliably measured from the $b$-$N_{H I}$ distribution at $z=0.1$, notwithstanding the substantial WHIM in the IGM. To investigate the effects of different feedback, we apply this inference methodology to mock spectra derived from the IllustrisTNG and Illustris simulations at $z=0.1$. The results suggest that the underlying $[T_0,\gamma]$ of both simulations can be recovered with biases as low as $|\Delta \log(T_0/\text{K})| \lesssim 0.05$ dex, $|\Delta \gamma | \lesssim 0.1$, smaller than the precision of a typical measurement. Given the large differences in the volume-weighted WHIM fractions between the three simulations (Illustris 38\%, IllustrisTNG 10\%, Nyx 4\%) we conclude that the $b$-$N_{H I}$ distribution is not sensitive to the WHIM under realistic conditions. Finally, we investigate the physical properties of the detectable Lyman-$\alpha$ absorbers, and discover that although their $T$ and $\Delta$ distributions remain mostly unaffected by feedback, they are correlated with the photoionization rate used in the simulation.

We probe dark matter-electron scattering using high-energy neutrino observations from the Sun. Dark matter (DM) interacting with electrons can get captured inside the Sun. These captured DM may annihilate to produce different Standard Model (SM) particles. Neutrinos produced from these SM states can be observed in IceCube and DeepCore. Although there is no excess of neutrinos from the Solar direction, we find that the current data-sets of IceCube and DeepCore set the strongest constraint on DM-electron scattering cross section in the DM mass range $10\,$GeV to $10^5\,$GeV. Our work implies that future observations of the Sun by neutrino telescopes have the potential to discover DM-electron interaction.

Using the Earth as a neutrino converter, tau neutrino fluxes from astrophysical point sources can be detected by tau-lepton-induced extensive air showers (EASs). Both muon neutrino and tau neutrino induced upward-going EAS signals can be detected by terrestrial, sub-orbital and satellite-based instruments. The sensitivity of these neutrino telescopes can be evaluated with the nuSpaceSim package, which includes the nuPyProp simulation package. The nuPyProp package propagates neutrinos ($\nu_\mu$, $\nu_\tau$) through the Earth to produce the corresponding charged leptons (muons and tau-leptons). We use nuPyProp to quantify the uncertainties from Earth density models, tau depolarization effects and photo-nuclear electromagnetic energy loss models in the charged lepton exit probabilities and their spectra. The largest uncertainties come from electromagnetic energy loss modeling, with as much as a 20-50% difference between the models. We compare nuPyProp results with other simulation package results.

We have attempted to mitigate the challenge of connecting the neutron star (NS) properties with the nuclear matter parameters that describe equations of state (EoSs). The outcomes of correlations of NS properties with individual nuclear matter parameters are at variance. A Principal Component Analysis is employed as a tool to uncover the connection between multiple nuclear matter parameters and the tidal deformability as well as the radius of neutron stars within the mass range of $1.2-1.8M_\odot$. The essential EoSs for neutron star matter at low densities have been derived using both uncorrelated uniform distributions and minimally constrained joint posterior distributions of nuclear matter parameters. For higher densities ($\rho > 0.32$fm$^{-3}$), the EoSs have been established through a suitable parameterization of the sound speed, which consistently maintains causality and gradually approaches the conformal limit. Our analysis reveals that in order to account for over 90\% of the variability in NS properties, it is crucial to consider two or more principal components, emphasizing the significance of employing multivariate analysis. To explain the variability in tidal deformability needs a greater number of principal components compared to those for the radius at a given NS mass. The contributions from iso-vector nuclear matter parameters to the tidal deformability and radius of NS decrease by $\sim$ 25\% with the increase in mass of NS 1.2$M_\odot$ to 1.8$M_\odot$ and accordingly those of iso-scalar nuclear matter parameters increase.

We investigate the properties of the hadron-quark mixed phase, often termed the \textit{pasta} phase, expected to exist in the cores of massive neutron stars. To construct the equations of state (EoS), we combine an analytical representation based on the APR EoS for hadronic matter with the MIT bag model featuring vector interactions for quark matter. For modeling the mixed phase, we utilize the compressible liquid drop model that consistently accounts for finite-size and Coulomb effects. Unlike most previous analyses that treated surface tension as a constant free parameter and neglected curvature tension, we employ microphysical calculations using the multiple reflection expansion formalism to determine these parameters, while also ensuring their self-consistency with the EoS. We construct an extensive set of mixed hybrid EoSs by varying model parameters, solve the stellar structure equations to obtain neutron star mass-radius relationships, and select the models that satisfy current astrophysical constraints. Our findings closely align with calculations using a constant surface tension in terms of EoS stiffness and resulting stellar structure. However, they reveal significant differences in the types of geometric structures and their prevalence ranges within the mixed phase. Specifically, curvature effects enhance the emergence of tubes and bubbles at high densities despite the large value of surface tension, while suppressing the existence of drops and rods at low densities.

This paper studies the quantum corrections to the Higgs inflation model in the context of the Einstein-Cartan (E-C) gravity in the large-$ N $ limit with $N$ being the number of real scalar components in Higgs. Recently, it is realized that the Higgs inflation in the E-C formalism smoothly connects those in the metric and the Palatini formalisms in the presence of a non-minimal coupling between the Higgs fields and the Nieh-Yan term. This motivates us to investigate the quantum corrections to the E-C Higgs inflation and to clarify how the Ricci curvature squared $ R^2 $ induced by the quantum corrections succeeds in Ultraviolet (UV)-extending the Higgs inflation in metric formalism while it fails in the Palatini case. We show that a generalized $ R^2 $-term required for the renormalization in the E-C formalism induces a new scalar degree of freedom (DoF), the scalaron, which gradually decouples with the system due to its increasing mass as approaching the Palatini limit. The presence of the scalaron extends the UV cutoff at vacuum of the original model except for the parameter space close to the Palatini limit. This UV-extension is expected to solve the strong coupling problem that may exist during (p)reheating in the absence of the scalaron.

We use Bayesian analysis in order to constrain the equation of state for nuclear matter from astrophysical data related to the recent measurements from the NICER mission, LIGO/Virgo collaboration, and probability distributions of mass and radius from other 12 sources, including thermonuclear busters, and quiescent low-mass X-ray binaries. For this purpose, we base our study on a relativistic hadronic mean field model including an $\omega-\rho$ interaction. Our results indicate optimal ranges for some bulk parameters at the saturation density, namely, effective mass, incompressibility, and symmetry energy slope ($L_0$). For instance, we find $L_0 = 50.79^{+15.16}_{-9.24}$ MeV (Case 1) and $L_0 = 75.06^{+8.43}_{-4.43}$ MeV (Case 2) in a $68\%$ confidence interval for the 2 cases analyzed (different input ranges for $L_0$ related to the PREX-II data). The respective parametrizations are in agreement with important nuclear matter constraints, as well as observational neutron star data, such as the dimensionless tidal deformability of the GW170817 event. From the mass-radius curves obtained from these best parametrizations, we also find the ranges of $11.97 \mbox{ km}\leqslant R_{1.4}\leqslant 12.73 \mbox{ km}$ (Case 1) and $12.34 \mbox{ km}\leqslant R_{1.4}\leqslant 13.06 \mbox{ km}$ (Case 2) for the radius of the $1.4M_\odot$ neutron star.

We generate three families of extended covariant density functionals of nuclear matter which have varying slope of symmetry energy and skewness at nuclear saturation density, but otherwise share the same basic parameters (symmetry energy, compressibility, saturation parameters, etc.) with the standard DDME2, DD2, and MPE functionals. Tables of the parameters of these new density functionals are given, that can be straightforwardly used in DDME2, DD2, and MPE parametrization-based codes. Furthermore, we provide tables of a large number of equation of state (81 for each family) which can be used in astrophysical simulations to assess the impact of variations of not-well-known slope of symmetry energy and skewness of nuclear system on the astrophysics of compact objects. We also provide tables of computed integral parameters (mass, radius, and tidal deformability) which can be used, e.g., for modeling gravitational wave forms. Finally, for the extended DDME2-based parametrization, we implement a first-order phase transition to quark matter to obtain a family of equation of state which accommodates a phase transition to quark matter. Analogous Tables of the EoS and integral parameters are provided for this case as well.

The origin of tiny neutrino mass is an unsolved puzzle leading to a variety of phenomenological aspects beyond the Standard Model (BSM). Among several interesting attempts, $U(1)$ gauge extension of Standard Model (SM) is a simple and interesting set-up where the so-called seesaw mechanism is incarnated by the addition of three generations of right-handed neutrinos followed by the breaking of $U(1)$ and electroweak symmetries. Such scenarios are anomaly free in nature appearing with a neutral BSM gauge boson ($Z^\prime$). In addition to that, there comes another open question regarding the existence of a non-luminous, hitherto unidentified object called Dark Matter (DM) originating from the measurement of its relic density. To explore properties of $Z^\prime$, we focus on chiral and flavored scenarios where $Z^\prime-$neutrinos interaction could be probed in the context of cosmic explosions like gamma-ray burst (GRB221009A, so far the highest energy), blazar (TXS 0506+056) and Active galaxy (NGC1068) respectively. The neutrino antineutrino annihilation produces electron-positron pair which could energize GRB through energy deposition. Taking the highest energy GRB under consideration and estimating the energy deposition rates we constrain $Z^\prime$ mass $(M_{Z^\prime})$ and the additional $U(1)$ coupling $(g_X)$ for chiral and flavored scenarios in the Schwarzchild, Hartle-Thorne and modified gravity frameworks. On the other hand, adding viable and alternative DM candidates in these models we study neutrino-DM scattering mediated by $Z^\prime$ in the $t-$ channel and estimate constraints on $g_X-M_{Z^\prime}$ plane using observed data of high energy neutrinos from cosmic blazar and active galaxy at the IceCube experiment. We compare our results with bounds obtained from different scattering, beam-dump and $g-2$ experiments.

We consider the relic density and positivity bounds for freeze-in scalar dark matter with general Higgs-portal interactions up to dimension-8 operators. When dimension-4 and dimension-6 Higgs-portal interactions are proportional to mass squares for Higgs or scalar dark matter in certain microscopic models such as massive graviton, radion or general metric couplings with conformal and disconformal modes, we can take the dimension-8 derivative Higgs-portal interactions to be dominant for determining the relic density via the 2-to-2 thermal scattering of the Higgs fields after reheating. We show that there is a wide parameter space for explaining the correct relic density from the freeze-in mechanism and the positivity bounds can curb out the dimension-8 derivative Higgs-portal interactions nontrivially in the presence of the similar dimension-8 self-interactions for Higgs and dark matter.

Fifty years after the first mention of light pollution in Science, the journal recently elevated this topic to the cover of its 16 June 2023 issue, highlighting the large impact on human and ecological health, circadian rhythms, migratory patterns, and more. We offer here the term noctalgia to express "sky grief" for the accelerating loss of the home environment of our shared skies - representing loss of science, heritage, millennia-old sky traditions, place-based language, and more - and summarize next steps to address the protection of our nighttime and daytime skies.