Papers for Friday, Dec 27 2024

Recent observations by Mistele et al. show that the circular velocity curves of isolated galaxies remain flat out to the largest radii probed so far, i.e. $\approx 1$ Mpc. The velocity decline beyond the expected virial radius is not observed. These results imply that the galaxy halo is in thermal equilibrium even at large radii where particles did not have time to relax. The galaxies must have already formed in the isothermal state. How is this possible? In the present note we try to understand the formation of galaxies with warm dark matter in the expanding universe.

Previous studies of the response of the Moon to gravitational waves have been carried out using analytical or semi-analytical models assuming ideal lunar structures. Such models are advantageous for their high-speed calculation but fail to account for the extremely heterogeneous subsurface and/or interior structures of the Moon. Numerical calculations are needed, but it is challenging to model the topography and lateral heterogeneity of the Moon. In addition, the computational cost is great especially when performing the GW simulation for a long time. As a first step towards overcoming the above difficulties, we employ a two-dimensional finite element method to numerically simulate the lunar response to gravitational waves. We verify our method by comparing our numerical results with those semi-analytical solutions. Based on such comparison, we also analyze the limitation of the two-dimensional simulation. Our work breaks a new way towards the precise simulation of realistic lunar response to gravitational waves in the future and lays down a solid foundation for three-dimensional numerical simulations.

Recent high-angular-resolution observations indicate the need for core growth to form high-mass stars. To understand the gas dynamics at the core scale in the very early evolutionary stages before being severely affected by feedback, we have conducted Atacama Large Millimeter/submillimeter Array (ALMA) observations toward a 70 $\mu$m dark massive clump, G337.541-00.082 as part of the Global and Local infall in the ASHES sample (GLASHES) program. Using dense gas tracers such as N$_2$H$^+$ ($J = 1-0$) and HNC ($J = 3-2$), we find signs of infall from the position-velocity diagram and more directly from the blue asymmetry profile in addition to the clump-scale velocity gradient. We estimate infall velocities from intermediate and low-mass cores to be 0.28-1.45 km s$^{-1}$, and infall rates to be on the order of 10$^{-4}$ to 10$^{-3}$ $M_\odot$ yr$^{-1}$, both are higher than those measured in low-mass star-forming regions by more than a factor of five and an order of magnitude, respectively. We find a strong correlation of the infall velocity with the nonthermal velocity dispersion, suggesting that infall may contribute significantly to the observed line width. Consistent with clump-fed scenarios, we show that the mass infall rate is larger for larger core masses and shorter distances to the clump center. Such high infall rates in cores embedded in IRDCs can be considered as strong signs of core growth, allowing high-mass star formation from intermediate-mass cores that would not initially form high-mass stars at their current mass.

We discuss a new source of gravitational waves (GWs) from first-order phase transitions. The collisions of bubbles of the new phase can efficiently produce particles that couple to the background field undergoing the transition, thereby transferring a significant fraction of the released vacuum energy into a distribution of inhomogeneous and dynamic particle populations that persist long after the bubbles have disappeared. We study the GWs produced by such particle distributions, showing that GWs arise from the quadrupolar anisotropy in the radiation emitted from the bubble collisions, and present a semi-analytical calculation of the two-point correlation function for the associated energy distributions. We find that this new contribution can qualitatively modify the overall GW signal from such phase transitions, creating a distinct shift in the spectral slope at low frequencies that could be observed by future GW experiments. It is therefore important to take this new contribution into account for any transition where the background field has significant self-coupling or couplings to other fields that could lead to efficient particle production at bubble collision.

We characterize the overdensity (spike) of fermionic dark matter (DM) particles around a supermassive black hole (SMBH) within a general relativistic analysis. The initial DM halo distribution is obtained by solving the equilibrium equations of a self-gravitating system of massive fermions at a finite temperature, according to the Ruffini-Argüelles-Rueda (RAR) model. The final fermionic DM spike is calculated around a Schwarzschild SMBH. We explore two possible interpretations for the origin and evolution of the SMBH seed. One corresponds to the traditional scenario proposed by Gondolo & Silk (1999), where a small BH mass of baryonic origin sits at the halo's center and grows adiabatically. The other one, from DM origin, where the dense and degenerate fermion core predicted by the RAR model grows adiabatically by capturing baryons until its gravitational collapse, providing a heavy SMBH seed, whose specific value depends on the fermion mass. We study different initial fermionic DM profiles that the theory allows. We show that overall dilute (i.e., Boltzmannian) fermionic DM develops the well-known spike with density profile $\rho \sim r^{-3/2}$. Instead, for semi-degenerate fermions with a dense and compact core surrounded by a diluted halo, we find novel spike profiles that depend on the particle mass and nature. In the more general case, fermionic spikes do not develop a simple power-law profile. Furthermore, the SMBH does not always imply an enhancement of the surrounding DM density; it might also deplete it. Thus, the self-consistent inclusion of the DM candidate nature and mass in determining the structure and distribution of DM in galaxies, including the DM spikes around SMBHs, is essential for the specification of DM astrophysical probes such as BH mergers, gravitational waves, or stellar orbits.

Here we briefly describe each of the modules that constitute Skynet's new "Astrophotography of the Multi-Wavelength Universe!", or MWU!, curriculum.

Sub-Neptunes occupy an intriguing region of planetary mass-radius space, where theoretical models of interior structure predict that they could be water-rich, where water is in steam and supercritical state. Such planets are expected to evolve according to the same principles as canonical H$_2$-He rich planets, but models that assume a water-dominated atmosphere consistent with the interior have not been developed yet. Here, we present a state of the art structure and evolution model for water-rich sub-Neptunes. Our set-up combines an existing atmosphere model that controls the heat loss from the planet, and an interior model that acts as the reservoir of energy. We compute evolutionary tracks of planetary radius over time. We find that planets with pure water envelopes have smaller radii than predicted by previous models, and the change in radius is much slower (within $\sim$10\%). We also find that water in the deep interior is colder than previously suggested, and can transition from plasma state to superionic ice, which can have additional implications for their evolution. We provide a grid of evolutionary tracks that can be used to infer the bulk water content of sub-Neptunes. We compare the bulk water content inferred by this model and other models available in the literature, and find statistically significant differences between models when the uncertainty on measured mass and radius are both smaller than 10\%. This study shows the importance of pursuing efforts in the modeling of volatile-rich planets, and how to connect them to observations.

We present Atacama Compact Array (ACA) Band-3 observations of the protocluster SPT2349$-$56, an extreme system hosting ${\gtrsim}\,12$ submillimeter galaxies (SMGs) at $z\,{=}\,4.3$, to study its integrated molecular gas content via CO(4-3) and long-wavelength dust continuum. The $\sim$30-hour integration represents one of the longest exposures yet taken on a single pointing with the ACA 7-m. The low-resolution ACA data ($21.0''\,{\times}\,12.2''$) reveal a 75% excess CO(4-3) flux compared to the sum of individual sources detected in higher-resolution Atacama Large Millimeter Array (ALMA) data ($1.0''\,{\times}\,0.8''$). Our work also reveals a similar result by tapering the ALMA data to $10''$. In contrast, the 3.2mm dust continuum shows little discrepancy between ACA and ALMA. A single-dish [CII] spectrum obtained by APEX/FLASH supports the ACA CO(4-3) result, revealing a large excess in [CII] emission relative to ALMA. The missing flux is unlikely due to undetected faint sources but instead suggests that high-resolution ALMA observations might miss extended and low-surface-brightness gas. Such emission could originate from the circum-galactic medium (CGM) or the pre-heated proto-intracluster medium (proto-ICM). If this molecular gas reservoir replenishes the star formation fuel, the overall depletion timescale will exceed 400Myr, reducing the requirement for the simultaneous SMG activity in SPT2349$-$56. Our results highlight the role of an extended gas reservoir in sustaining a high star formation rate (SFR) in SPT2349$-$56, and potentially establishing the ICM during the transition phase to a mature cluster.

The recent and brightest GRB 221009A observed by LHAASO marked the first detection of the onset of TeV afterglow, with a total of 7 GRBs exhibiting very high energy (VHE) afterglow radiation. However, consensus on VHE radiation of GRBs is still lacking. Multi-wavelength studies are currently a primary research method for investigating high-energy $\gamma$-ray astronomy. The limited sample of VHE GRBs, combined with their transient nature, hinders the progress of physical studies of GRBs. This paper aims to obtain useful information for GRB research through the properties of blazars, which share significant similarities with GRBs. By fitting high-quality and simultaneous multiwavelength spectral energy distributions with a one-zone leptonic model, the study explores the similarity of radiation properties of blazars and GRBs. A tight correlation between synchrotron and synchrotron self-Compton (SSC) emission luminosities suggests that blazars and GRBs share similar radiation mechanisms, to be specific, synchrotron radiation produces the observed X-ray photons, which also serve as targets for electrons in the SSC process. We hope that ground-based experiments can observe more GRBs in sub-TeV to confirm these findings.

The Herzberg Extensible Adaptive optics Real-Time Toolkit (HEART) is a complete framework written in C and Python for building next-generation adaptive optics (AO) system real-time controllers, with the performance needed for extremely large telescopes. With numerous HEART-based RTCs now in their design or build phases, each with different AO algorithms, target hardware, and observatory requirements, continuous automated builds and tests are a cornerstone of our development effort. In this paper we describe the many levels of testing that we perform, from low-level unit tests of individual functions, to more complex component and system-level tests that verify both numerical correctness and execution performance. Incorporating extensive testing into HEART since its inception has allowed us to continuously (and confidently) refactor and extend it to both meet the changing needs of local on-sky experiments, as well as those of the several major facility instruments that we are developing.

Mechanical properties of minerals in planetary materials are not only interesting from a fundamental point of view but also critical to the development of future space missions. Here we present nanoindentation experiments to evaluate the hardness and reduced elastic modulus of olivine, (Mg,Fe)2SiO4, in meteorite NWA 12008, a lunar basalt. Our experiments suggest that the olivine grains in this lunaite are softer and more elastic than their terrestrial counterparts. This may be attributed to macroscopic effects, like increased porosity, or even to modifications at the chemical bond scale. We have performed high-pressure X-ray diffraction (HP-XRD) measurements to probe the elastic compressibility properties on this meteorite and, for comparison purposes, on three ordinary chondrites. The HP-XRD results suggest that the axial compressibility of the orthorhombic $b$ lattice parameter of olivine is higher in NWA 12008 and also in the highly-shocked Chelyabinsk meteorite, relative to terrestrial olivine. The origin of the observed differences may be the consequence of a combination of factors reflecting their complex history. The combined study by nanoindentations and HP-XRD of the mechanical and elastic properties of meteorites and returned samples opens up a new avenue to characterize these materials that will be crucial for future extraterrestrial resource utilization purposes.

Persistence effects in HgCdTe infrared detectors cause significant artifacts that can impact the quality of science observations for up to many hours after exposure to bright/saturating sources. This problem will have a substantially greater impact on viable observing modes for infrared cameras on future ELTs due to the leap in sensitivities that is expected. In this paper we present new results from an updated test system that was previously used to prototype ``on-detector guide windows'' to provide fast T/T feedback to AO systems, interleaved with simultaneous (slow) full-frame readouts for science. We now explore the possibility of continuously resetting these small regions of the detector that are illuminated with a compact source as a strategy for mitigating persistence, using two different detectors. While our results generally show promise for this observing strategy, we found for one of our detectors that the combination of fast localized resets with intense illumination can introduce a potentially problematic persistent change in local reset levels.

Globular clusters (GCs) are fossil tracer of formation and evolution of galaxies and studying their properties provides crucial insights into past formation and interaction events. We study the properties of GC candidates in 1.25 $\times$ 1.03 sq. degrees area centred on NGC 5018 group of galaxies using deep, wide field and multi-band observations obtained with VST. We derived photometric catalogues of compact and extended sources and identified GC candidates using a set of photometric and morphometric selection parameters. A GC candidates catalogue is inspected using a statistical background decontamination technique. The 2D distribution map of GC candidates reveals an overdensity of sources on brightest member of group, NGC 5018. No significant GC overdensities are observed in other bright galaxies of group. We report discovery of a candidate local nucleated LSB dwarf galaxy possibly in tidal interaction with NGC 5018. The 2D map also reveals an intra-group GC population aligning with bright galaxies and along intra-group light component. Radial density profile of GC candidates in NGC 5018 follows galaxy surface brightness profile. Colour profile of GC candidates centred on this galaxy shows no evidence of well-known colour bimodality, which is instead observed in intra-group population. From GC luminosity function (GCLF) analysis, we find a low specific frequency $S_{\!\rm N}=0.59 \pm 0.27$ for NGC 5018, consistent with previous results. This relatively low $S_{\!\rm N}$ and lack of colour bimodality might be due to a combination of observational data limitations and the post-merger status of NGC 5018, which might host a population of relatively young GCs. For intra-group GC population, we obtain a lower limit of $S_{\!\rm {N,gr}}\sim0.6$. Using GCLF as a distance indicator, we estimate that NGC 5018 is located $38.0 \pm 7.9$ Mpc away, consistent with values in the literature..

Most Kreutz family sungrazing comets are discovered only days before perihelion, severely limiting observational opportunities to study their physical nature and decay. Kreutz sungrazer C/2024 S1 (ATLAS) was discovered a month before reaching its perihelion distance of 0.008 au, allowing physical observations from both ground- and space-based telescopes. We present observations from 0.9 au to 0.4 au using the Nordic Optical Telescope showing that 1) nucleus disintegration was on-going already at 0.7 au pre-perihelion, 2) the activity varied unpredictably with distance and 3) the nucleus radius was $<$600 m (red geometric albedo 0.04 assumed). We also use coronagraphic observations from the STEREO-A spacecraft to study C/2024 S1 at heliocentric distances $\lesssim$0.1 au. We find that the coma scattering cross-section peaked near 0.075 au and faded progressively, by a factor $\sim$20, towards the last observation at 0.02 au. We interpret the near-perihelion fading as a result of the sublimation of refractory coma grains, beginning at blackbody temperatures $\sim$1000 K, consistent with olivine composition. The comet was not detected after perihelion. We consider processes operating to destroy the nucleus when near perihelion, concluding that rotational instability and sublimation losses work together towards this end, even before entry of the comet into the Roche lobe of the Sun.

We investigate the bar formation process using $N$-body simulations across the Toomre's parameter $Q_{min}$ and central mass concentration (CMC), focusing principally on the formation timescale. Of importance is that, as suggested by cosmological simulations, disk galaxies have limited time of $\sim 8$ Gyr in the Universe timeline to evolve secularly, starting when they became physically and kinematically steady to prompt the bar instability. By incorporating this time limit, bar-unstable disks are further sub-divided into those that establish a bar before and after that time, namely the normal and the slowly bar-forming disks. Simulations demonstrate that evolutions of bar strengths and configurations of the slowly bar-forming and the bar-stable cases are nearly indistinguishable prior to $8$ Gyr, albeit dynamically distinct, while differences can be noticed afterwards. Differentiating them before $8$ Gyr is possible by identifying the proto-bar, a signature of bar development visible in kinematical maps such as the Fourier spectrogram and the angular velocity field, which emerges in the former group $1-2$ Gyr before the fully developed bar, whereas it is absent in the latter group until $8$ Gyr and such bar-stable disk remains unbarred until at least $10$ Gyr. In addition, we find complicated interplays between $Q_{min}$ and CMC in regulating the bar formation. Firstly, disk stabilization requires both high $Q_{min}$ and CMC. Either high $Q_{min}$ or high CMC only results in slow bar formation. Secondly, some hot disks can form a bar more rapidly than the colder ones in a specific range of $Q_{min}$ and CMC.

The geometry and kinematics of the broad-line region (BLR) in AGNs are still unclear, which is crucial for studying the physics and evolution of supermassive black holes (SMBHs) and AGNs. The broad-line profile provides valuable information on BLR geometry and kinematics. In this work, we explore the evolution of line profiles in variable AGNs based on the BLR model of Czerny \& Hryniewicz, where the BLR is driven by the radiation pressure acting on dust at the surface layers of the accretion disk. The line profiles in the low-Eddington-ratio regime show a double-peak profile, which will become a single peak at high Eddington ratios. The high metallicity of $Z\gtrsim 5Z_{\odot}$ is required to reproduce the observational anti-correlation between the peak separation of broad lines and the Eddington ratio for a sample of AGNs. For the broad lines in variable AGNs, it will take several years to several decades to change their line profile if the disk luminosity suffers strong variation in a much shorter timescale. More monitoring of the broad line and continuum in strongly variable AGNs can shed special light on BLR physics.

Robotic telescope networks play an important role in capturing early and bright optical afterglows, providing critical insights into the energetics and emission mechanisms of GRBs. In this study, we analyze GRB 230204B, an exceptionally energetic and multi-pulsed long GRB, detected by the Fermi GBM and MAXI detectors, with an isotropic equivalent gamma-ray energy exceeding 10$^{54}$ erg. Time-resolved spectral analysis reveals a transition in the prompt emission from hard (sub-photospheric dominated) spectra during early pulses to softer (synchrotron radiation dominated) spectra in later pulses, indicative of a hybrid jet composition. We report the discovery and characterization of the optical afterglow using the MASTER and BOOTES robotic telescope networks, alongside long-term radio observations extending to 335 days post-burst with the ATCA. At ~1.3 ks post-burst, the optical luminosity was exceptionally high, surpassing even other bright GRBs, such as GRB 221009A (the ``BOAT"). Multi-wavelength modeling, incorporating data from MASTER, BOOTES, DOT, Swift/XRT, and radio observations, was conducted using an external ISM forward-shock top-hat jet model with afterglowpy. The results reveal a narrow and highly collimated jet with a circumburst density of n$_{0}$ ~ 28.12 cm$^{-3}$, kinetic energy E$_{K}$ ~ 4.18 x 10$^{55}$ erg, and a relatively low value of $\epsilon_{B}$ = 2.14 x 10$^{-6}$, indicating shock-compression of the magnetic field in the surrounding interstellar medium. We constrained a low radiative efficiency of ~ 4.3 %. This study highlights the indispensable contribution of robotic networks to early afterglow observations and advances our understanding of GRB 230204B unique characteristics and underlying jet physics.

A potential difference of 1.3 Giga-Volts (GV) was inferred across a thundercloud using data from the GRAPES-3 muon telescope (G3MT). This was the first-ever estimation of gigavolt potential in thunderstorms, confirming prediction of C.T.R. Wilson almost a century ago. To infer the thundercloud potential required acceleration of muons in atmospheric electric field to be incorporated in the Monte Carlo simulation software CORSIKA. The G3MT records over 4 billion muons daily that are grouped into 169 directions covering 2.3 sr sky. This enabled changes as small as 0.1% in the muon flux on minute timescale, caused by thunderstorms to be accurately measured. But that requires high statistics simulation of muon fluxes in thunderstorm electric fields. The CORSIKA offers a choice of several generators for low- (FLUKA, GHEISHA, and UrQMD) and high-energy (SIBYLL, EPOS-LHC, and QGSJETII) hadronic interactions. Since it is unclear which combination of the low- and high-energy generators provides the correct description of hadronic interactions, all nine combinations of generators were explored, and they yielded thundercloud potentials ranging from 1.3 GV to 1.6 GV for the event recorded on 1 December 2014. The result of SIBYLL-FLUKA combination yielded the lowest electric potential of 1.3 GV was reported. Furthermore, another seven major thunderstorm events recorded between April 2011 and December 2020 were analyzed to measure the dependence of their thundercloud potential on the hadronic interaction generators. It is observed that the low-energy generators produce larger variation ($\sim$14%) in thundercloud potential than the high-energy generators ($\sim$8%). This probably reflects the fact that the GeV muons are predominantly produced in low-energy ($<$80 GeV) interactions, which effectively magnifies the differences in the meson production cross-sections among the low-energy generators.

Magnetic helicity is a key geometrical parameter to describe the structure and evolution of solar coronal magnetic fields. The accumulation of magnetic helicity is correlated with the non-potential magnetic field energy, which is released in the solar eruptions. Moreover, the relative magnetic helicity fluxes can be estimated only relying on the line-of-sight magnetic field (e.g., Demoulin and Berger Sol. Phys. 215, 203, 2003). The payload Full-disk MagnetoGraph (FMG) on the Advanced Space-based Solar Observatory (ASO-S) currently has been supplying the continuous evolution of line-of-sight magnetograms for the solar active regions, which can be used to estimate the magnetic helicity flux. In this study, we useeight hours line-of-sight magnetograms of NOAA 13273, when at which the Sun-Earth direction speed of the satellite is zero to avoid the oscillation of the magnetic field caused by the Doppler effect on polarization measurements. We obtain the helicity flux by applying Fast Fourier Transforms (FFT) and local correlation tracking (LCT) methods to obtain the horizontal vector potential field and the motions of the ine-sf-sight polarities. We also compare the helicity flux derived using data from the Heliosesmic and Magnetic Imager (HMI) on board the Solar Dynamics Observatory (SDO) and the same method. It is found that the flux has the same sign and the correlation between measurements is 0.98. The difference of the absolute magnetic helicity normalized to themagnetic flux is less than 4%. This comparison demonstrates the reliability of ASO-S/FMG data and that it can be reliably used in future studies

An instability criterion in the MHD with the open boundary of magnetic field is proposed in this paper. We use a series of linear force-free extrapolation field, in which the normal part of magnetic field is fixed, to obtain the linear fitting coefficient called relative alpha by using the co-joined value of magnetic free energy and magnetic flux at the open boundary ($E_f \Phi ^2$) and the square of relative magnetic helicity ($H_R^2$). We calculate this coefficient of the magnetic field above active regions NOAA~8210 and NOAA~11429 obtained by the photospheric-data-driven magnetohydrodynamics (MHD) model. It is found that the fitting coefficient is a good proxy of the criterion to indicate the occurrence of instability after which the magnetic reconnection happens and caused the fast release of magnetic energy. We also applied this method to the continuous evolution of three-dimension magnetic field of NOAA~11158 based on the measurement of photospheric vector magnetic field of SDO/HMI by the Non-linear Force-Free (NLFF) extrapolation method. The calculated coefficient when the major flare happened based on the extrapolation data is very close to the expected ones, which perfectly reflects the occurrence of instability and the difference is even less than 7\%. This relative alpha is very helpful to evaluate how far it is from the instability in the MHD and quantitatively estimate the occurrence of solar eruption in the space weather forecast.

Type Ia supernovae are triggered by accretion onto a white dwarf from a companion which is most likely Roche lobe-filling at the time of the explosion. The collision between the ejecta and a surviving companion carves out a conical wake, which could manifest as an asymmetry when the ejecta reaches the remnant phase. We simulate the companion interaction using the Athena++ hydrodynamics solver to determine the ejecta structure for a double-degenerate type Ia supernova. Ejecta in the wake is of lower density and higher velocity than the unperturbed ejecta. We then evolve the ejecta for several thousand years using the expanding-grid code Sprout. The forward shock within the wake is initially indented, but becomes spherical after roughly a thousand years due to transverse motion of shocked ejecta that fills the wake. The reverse shock travels quickly within the wake, leading to an off-center convergence of the reverse shock and leaving the remnant with an asymmetrical core. This also draws material from the interstellar medium deep into the remnant, eventually reaching the center. Large Rayleigh-Taylor plumes are found around the edge of the wake, creating a toroidal structure composed primarily of ejecta. Estimates of the thermal X-ray emission show that such remnants exhibit observable asymmetries for thousands of years.

We present a study on rotation measure (RM) of the quasar 3C 273. This analysis aims to discern the magnetic field structure and its temporal evolution. The quasar 3C 273 is one of the most studied active galactic nuclei due to its high brightness, strong polarization, and proximity, which enables resolving the transverse structure of its jet in detail. We used polarized data from 2014, collected at six frequencies (5, 8, 15, 22, 43, 86 GHz) with the Very Long Baseline Array, to produce total and linear polarization intensity images, as well as RM maps. Our analysis reveals a well-defined transverse RM gradient across the jet, indicating a helical, ordered magnetic field that threads the jet and likely contributes to its collimation. Furthermore, we identified temporal variations in the RM magnitude when compared with prior observations. These temporal variations show that the environment around the jet is dynamic, with changes in the density and magnetic field strength of the sheath that are possibly caused by interactions with the surrounding medium.

The High Energy cosmic-Radiation Detection (HERD) facility is a dedicated high energy astronomy and particle physics experiment planned to be installed on the Chinese space station, aiming to detect high-energy cosmic rays (\si{\giga\electronvolt} $\sim$ \si{\peta\electronvolt}) and high-energy gamma rays (> \SI{500}{\mega\electronvolt}). The Plastic Scintillator Detector (PSD) is one of the sub-detectors of HERD, with its main function of providing real-time anti-conincidence signals for gamma-ray detection and the secondary function of measuring the charge of cosmic-rays. In 2023, a prototype of PSD was developed and tested at CERN PS\&SPS. In this paper, we investigate the position response of the PSD using two reconstruction algorithms: the classic dual-readout ratio and the deep learning method (KAN \& MLP neural network). With the latter, we achieved a position resolution of 2 mm ($1\sigma$), which is significantly better than the classic method.

Using the cross-correlation data from the first three observing runs of the LIGO-Virgo-KAGRA Collaboration, we search for a gravitational-wave background (GWB) from primordial black holes, arising from the superposition of compact binary coalescence events. We consider both early and late binary formation mechanisms, and perform Bayesian parameter inference, investigating different prior distributions of the model parameters. From the non-detection of the GWB, we provide constraints on the fraction of primordial black holes contributing to the present dark matter energy density.

We present the results of determining the parameters characterizing the shape and orientation of residual velocity ellipsoids from the Gaia DR3 red giants and subgiants. We show the distribution of velocity dispersions in the Galactic plane obtained from three components of the spatial velocity, as well as the coordinate distribution of the intersection points of the velocity ellipsoid axes with the celestial sphere, in particular the deviations of the longitudes and latitudes of the vertices of stellar regions located within spheres with a radius of 1 kpc centered in the Galactic mid-plane. The area of the Galactic disk under study is in the range of Galactocentric coordinates 0 < R < 15 kpc and $120^\circ < \theta < 240^\circ$. We show that the vertex deviations in some regions of the Galactic mid-plane can reach $30^\circ$ in longitude, and $15^\circ$ in latitude. This indicates the presence of kinematic distortions of the stellar velocity field, especially noticeable in the angular range of $150^\circ < \theta < 210^\circ$ at a distance of approximately 13 kpc. We propose the angles of deviation of longitudes and latitudes of the ellipsoid axes of residual stellar velocities to be considered as kinematic signatures of various Galactic deformations determined from real fields of spatial velocities. We present the distribution of parameters characterizing the shapes of velocity ellipsoids, as well as their distribution of the semi-axes length ratios. We note a local feature in this distribution and in the distribution of the elongation measurements of the ellipsoids. We perform a comparison of the results obtained from the tensor of deformation velocities and from the observed spatial velocities.

SPARC4 is a new astronomical instrument developed entirely by Brazilian institutions, currently installed on the 1.6-m Perkin-Elmer telescope of the Pico dos Dias Observatory. It allows the user to perform photometric or polarimetric observations simultaneously in the four SDSS bands (g, r, i, and z). In this paper, we describe the control system developed for SPARC4. This system is composed of S4ACS, S4ICS, and S4GUI softwares and associated hardware. S4ACS is responsible for controlling the four EMCCD scientific cameras (one for each instrument band). S4ICS controls the sensors and motors responsible for the moving parts of SPARC4. Finally, S4GUI is the interface used to perform observations, which includes the choice of instrument configuration and image acquisition parameters. S4GUI communicates with the instrument subsystems and with some observatory facilities, needed during the observations. Bench tests were performed for the determination of the overheads added by SPARC4 control system in the acquisition of photometric and polarimetric series of images. In the photometric mode, SPARC4 allows the acquisition of a series of 1400 full-frame images, with a deadtime of 4.5 ms between images. Besides, several image series can be concatenated with a deadtime of 450 ms plus the readout time of the last image. For the polarimetric mode, measurements can be obtained with a deadtime of 1.41 s plus the image readout time between subsequent waveplate positions. For both photometric and polarimetric modes, the user can choose among operating modes with image readout times between 5.9 ms and 1.24 s, which ultimately defines the instrument temporal performance.

In this study, we present the results of the relationship between spectral type (ST) and the projected stellar rotational velocity ($vsini$), utilising a sample of approximately 50,000 single stars across a range of evolutionary stages. The STs of the stars included in this study span a broad range, from O0 to M9. We examine the stars in our data set, which has been divided into two groups according to ST and luminosity class (LC). The groups were conducted an investigation into the relationship between the mean $vsini$ ($\langle vsini \rangle$) and STs, as well as the dependence of $\langle vsini \rangle$ on STs and LCs. The rationale for investigating the two subgroups separately is to take into account for the evolutionary status of the stars and ascertain the impact on stellar rotation. The results demonstrate a notable decline in $\langle vsini \rangle$ as the spectral type progresses from early to late types. In particular, we found a significant decrease in $\langle vsini \rangle$ values, amounting to approximately 100 km/s, between hot stars (STs O0 to F2) and cool stars (STs F2 to M9). Moreover, a reduction in $\langle vsini \rangle$ is discernible as stars evolve, with this trend being most pronounced in evolutionary stages beyond the subgiant phase.

Model-independent bounds on the Hubble constant $H_0$ are important to shed light on cosmological tensions. We work out a model-independent analysis based on the sum rule, which is applied to late- and early-time data catalogs to determine $H_0$. Through the model-independent Bézier interpolation of the observational Hubble data (OHD) and assuming a flat universe, we reconstruct the dimensionless distances of the sum rule and apply them to strong lensing data to derive constraints on $H_0$. Next, we extend this method to the high-redshift domain including, in other two separated analyses, gamma-ray burst (GRB) data sets from the well-established Amati and Combo correlations. In all three analyses, our findings agree at $1\sigma$ level with the $H_0$ determined from type Ia supernovae (SNe Ia), and only at $2\sigma$ level with the measurement derived from the cosmic microwave background (CMB) radiation. Our method evidences that the bounds on $H_0$ are significantly affected by strong lensing data, which favor the local measurement from SNe Ia. Including GRBs causes only a negligible decrease in the value of $H_0$. This may indicate that GRBs can be used to trace the expansion history and, in conjunction with CMB measurements, may heal the Hubble tension and accommodate to the flat $\Lambda$CDM paradigm purported by CMB data.

Massive stars play an important role in the Universe. Unlike low-mass stars, the formation of these objects located at great distances is still unclear. It is expected to be governed by some combination of self-gravity, turbulence, and magnetic fields. In this work, we aim to study the chemical and physical conditions of dense clumps at different evolutionary stages. We performed observations towards 5 regions of massive star and stellar cluster formation (L1287, S187, S231, DR 21(OH), NGC 7538) with the IRAM-30m telescope. We covered the 2 and 3$-$4 mm wavelength bands and analysed the lines of HCN, HNC, HCO$^+$, HC$_3$N, HNCO, OCS, CS, SiO, SO$_2$, and SO. Using astrodendro algorithm on the 850 $\mu$m dust emission data from the SCUBA Legacy catalogue, we determined the masses, H$_2$ column densities, and sizes of the clumps. Furthermore, the kinetic temperatures, molecular abundances, and dynamical state were obtained. The Red Midcourse Space Experiment Source survey (RMS) was used to determine the clump types. A total of 20 clumps were identified. Three clumps were found to be associated with the Hii regions, 10 with young stellar objects (YSOs), and 7 with submillimetre emission. The clumps have typical sizes of about 0.2 pc and masses from 1 to 10$^{2}\,M_\odot$, kinetic temperatures ranging from 20 to 40 K and line widths of $\rm H^{13}CO^{+} (1-0)$ approximately 2 $\rm km\,s^{-1}$. We found no significant correlation in the line width$-$size and the line width$-$mass relationships. However, a strong correlation is observed in mass$-$size relationships. The virial analysis indicated that three clumps are gravitationally bound. Furthermore, we suggested that magnetic fields of about 1 mG provide additional support for clump stability. The molecular abundances relative to H$_2$ are approximately $10^{-10}-10^{-8}$.

I review the properties and discuss some of the puzzling aspects of the unique binary system composed of the luminous hot subdwarf HD 49798 and a white dwarf with mass of 1.2 solar masses and spin period of 13.2 s. This is one of the few massive white dwarfs with a dynamically measured mass and the one with the shortest spin period. It emits pulsed X-rays with a very soft spectrum, powered by accretion from the tenuous stellar wind of its companion of sdO spectral type. The current level of mass accretion cannot provide enough angular momentum to explain the small, but precisely measured, spin-up rate of 72 nanoseconds per year, which is instead best interpreted as the result of the radial contraction of this young white dwarf. The higher mass transfer rate expected during the future evolutionary stages of HD 49798 will drive the white dwarf above the Chandrasekhar limit, but the final fate, a type Ia SN explosion or the collapse to a millisecond pulsar, is uncertain.

The flat spectrum radio quasar PKS 1510-089 is one of the most active blazars across the entire electromagnetic spectrum, displaying periods of flaring activity. This study explores its spectral variability over a decade. By employing the non-thermal dominance parameter, we analyze the H$\beta$ and $\lambda5100\text{ Å}$ continuum light curves, as well as the full width at half maximum of the H$\beta$ emission line, to identify whether the primary source of the continuum emission is the accretion disk or the jet during activity periods. Our results shows an anti-correlation between the full width at half maximum and the luminosity of the H$\beta$ emission line across all datasets. This indicates, that variations in H$\beta$ luminosity consistently reflects the canonical broad-line region, irrespective of whether the primary ionizing source is the accretion disk or the jet. The anti-correlation persisted when comparing the full width at half maximum of H$\beta$ against the luminosity at $\lambda5100\text{ Å}$ in the disk dominance regime. These findings, along with the observation that flaring events in the $\lambda5100\text{ Å}$ continuum, attributed to the jet, coincide with flares in the H$\beta$ emission line, suggest that the base of the jet is located within the broad-line region. Based on the 219 spectra within the disk dominance regime, we estimated a mean black hole mass of $M_{BH}=2.85\pm0.37\times10^{8}\: M_{\odot}$.

Red supergiant stars (RSGs) represent the final evolutionary phase of the majority of massive stars and hold a unique role in testing the physics of stellar models. Eighty eight RSGs in the Small Magellanic Cloud (SMC) were recently found to have an ultra-violet excess that was attributed to a B-type companion. We present follow-up Hubble Space Telescope (HST) Space Telescope Imaging Spectrograph (STIS) ultra-violet (1700 -- 3000\,Å) spectroscopy for 16 of these stars to investigate the nature of the UV excess and confirm the presence of a hot companion. In all cases we are able to confirm that the companion is a main-sequence B-type star based on the near-UV continuum. We determine effective temperatures, radii and luminosities from fitting the UV continuum with TLUSTY models and find stellar parameters in the expected range of SMC B-type stars. We display these results on a Hertzsprung--Russell diagram and assess the previously determined stellar parameters using UV photometry alone. From this comparison we conclude that UV photometric surveys are vital to identify such companions and UV spectroscopy is similarly vital to characterise the hot companions. From a comparison with IUE spectra of 32 Cyg, a well known RSG binary system in the Galaxy, four targets display evidence of being embedded in the wind of the RSG, like 32 Cyg, although none to the more extreme extent of VV Cep. The ages of six targets, determined via the stellar parameters of the hot companions, are found to be in tension with the ages determined for the RSG. A solution to this problem could be binary mass-transfer or red straggler stars.

The study of high-energy gamma-ray emission from gamma-ray bursts (GRBs) involves complex synchrotron radiation and synchrotron self-Compton scattering (SSC) mechanisms with multiple parameters exhibiting a wide distribution. Recent advancements in GRB research, particularly the observation of very high energy (VHE, $\rm >100~GeV$) radiation, have ushered in a new era of multiwavelength exploration, offering fresh perspectives and limitations for understanding GRB radiation mechanisms. This study aimed to leverage VHE observations to refine constraints on synchrotron + SSC radiation from electrons accelerated by forward shocks. By analyzing two external environments - the uniform interstellar medium and stratified stellar wind medium, we conducted spectral and variability fitting for five specific bursts (GRB~180720B, GRB~190114C, GRB~190829A, GRB~201216C, and GRB~221009A) to identify the optimal parameters characterizing these events. A comparative analysis of model parameter distributions with and without VHE radiation observations reveals that the magnetic energy equipartition factor $\epsilon_B$ is more concentrated with VHE emissions. This suggests that VHE emissions may offer greater constraints on this microphysical parameter. Additionally, we found that the energy budget between VHE and keV-MeV $\gamma$-ray emissions under the SSC radiation exhibits an almost linear relationship, which may serve as a tool to differentiate radiation mechanisms. We anticipate future statistical analyses of additional VHE bursts to validate our findings.

Using a three-flavor Nambu--Jona-Lasinio model to describe the charge-parity violating effects through axion field, we investigate the axion effects on quark matter and quark-matter cores in massive hybrid stars. The properties of quark matter vary with the scaled axion field $a/f_a$ in a periodic manner, with a period of $2\pi$. Within the range from 0 to $\pi $, axion field decrease the baryon chemical potential of the first-order phase transition, leading to an increase in normalized pressure and stiffening of the quark matter equation of state. The effect of axions on hybrid star matter that includes the hadron-quark phase transition is contrary to expectations. The axion field shifts the onset of the hadron-quark mixed phase to lower densities but slightly softens the equation of state of the mixed phase matter, which also results in a slight decrease in the maximum mass and corresponding radius of the hybrid stars. However, we also find that the lowering of the onset of the mixed phase significantly increases the radius and mass of the quark-matter core in the hybrid star. Therefore, our results indicate with axion effects, a sizable quark-matter core can appear in $2M_{\odot}$ massive neutron stars.

The construction of high-resolution shock-capturing schemes is vital in producing highly accurate gravitational waveforms from neutron star binaries. The entropy based flux limiting (EFL) scheme is able to perform fast converging binary neutron star merger simulations reaching up to fourth-order convergence in the gravitational waveform phase. In these results the EFL method was used only in the dynamical evolution of initial data constructed with the Lorene library. Here, we extend the use of the EFL method to the construction of eccentricity reduced initial data for neutron star binaries and present several new BNS simulations resulting from such initial data and show for the first time up to optimal fifth-order convergence in the gravitational waveform phase.

A scalar field coupled conformally and disformally to matter affects both the linear memory effect for binary systems on hyperbolic orbits, as well as the kick velocity for binaries on bound or unbound orbits. We study these corrections in detail, their order of magnitude, and discuss their detectability. In particular, we find that the disformal interaction does not contribute to the memory effect and the emitted power spectrum at zero frequency. The conformal interaction corrects the GR linear memory and the quadrupole emitted power at zero frequency resulting in a breaking of the GR memory-power spectrum relationship. On the other hand, disformal interactions give rise to a change of momentum for the centre of mass. Hence, measuring both the linear memory effect and the kicks for hyperbolic orbits would give access to the conformal and disformal couplings of nearly massless scalars to matter.

As next-generation gravitational-wave (GW) observatories approach unprecedented sensitivities, the need for robust methods to analyze increasingly complex, overlapping signals becomes ever more pressing. Existing matched-filtering approaches and deep-learning techniques can typically handle only one or two concurrent signals, offering limited adaptability to more varied and intricate superimposed waveforms. To overcome these constraints, we present the UnMixFormer, an attention-based architecture that not only identifies the unknown number of concurrent compact binary coalescence GW events but also disentangles their individual waveforms through a multi-decoder architecture, even when confronted with five overlapping signals. Our UnMixFormer is capable of capturing both short- and long-range dependencies by modeling them in a dual-path manner, while also enhancing periodic feature representation by incorporating Fourier Analysis Networks. Our approach adeptly processes binary black hole, binary neutron star, and neutron star-black hole systems over extended time series data (16,384 samples). When evaluating on synthetic data with signal-to-noise ratios (SNR) ranging from 10 to 50, our method achieves 99.89% counting accuracy, a mean overlap of 0.9831 between separated waveforms and templates, and robust generalization ability to waveforms with spin precession, orbital eccentricity, and higher modes, marking a substantial advance in the precision and versatility of GW data analysis.

We give the connection formulae for ordinary differential equations with 5 and 6 (and in principle can be generalized to more) regular singularities from the data of instanton partition functions of quiver gauge theories. We check the consistency of these connection formulae by numerically computing the quasinormal modes (QNMs) of Reissner-Nordström de Sitter (RN-dS) blackhole. Analytic expressions are obtained for all the families of QNMs, including the photon-sphere modes, dS modes, and near-extremal modes. We also argue that a similar method can be applied to the dS-Kerr-Newman blackhole.

The composition and equation of state (EoS) of dense matter relevant to compact stars are quite inconclusive. However, certain observational constraints on the structural properties of compact stars help us constrain the EoS to a fair extent. Moreover, gravitational asteroseismology gives us a notion of the interior of the compact star viz., its composition and EoS. The next generation gravitational wave (GW) detectors are likely to detect several oscillation mode frequencies of the GWs, specially the fundamental ($f$) mode frequency. In this work we compute the $f$ and the first pressure ($p_1$) mode frequencies ($f_f$ and $f_{p1}$, respectively) with different compositions viz., hadronic matter, quark matter, and hybrid star (HS) matter. For HSs, we also study the gravity ($g$) mode frequency ($f_g$) arising due to the discontinuity in density. For each phase we also study the correlation between the oscillation frequencies of the 1.4 $M_{\odot}$ and 2.01 $M_{\odot}$ compact stars with the other different properties. We find that the different possible composition of compact stars substantially affects the different oscillation mode frequencies. However, the mass-scaled angular $f$ mode frequency ($\omega_f M$) varies universally with compactness ($C$) for all hadronic, quark and hybrid stars. The $f$ mode frequency ($f_{f_{1.4}}$) of the canonical 1.4 $M_{\odot}$ compact star, obtained with the different composition, is quite correlated with the canonical radius ($R_{1.4}$) and tidal deformability ($\Lambda_{1.4}$) while $f_{p_{1.4}}$ is well correlated with the slope parameter ($L$) of the symmetry energy. We also show that $f_{g_{1.4}}$ of the HSs varies almost linearly with $\Lambda_{1.4}$. Should $g$ modes be detected, they could not only support the existence of HSs, but the magnitude of $f_g$ could be useful to understand the strength of quark repulsion in HSs.