Papers for Monday, Jan 06 2025

The general theory of relativity is currently established as the most precise theory of gravity supported by observations, and its application is diverse ranging from astronomy to cosmology, while its application to astrophysics has been restricted only to compact stars due to the assumption that the Newtonian approximation is sufficient for celestial bodies with medium density such as the sun. Surprisingly, the recent research of the author has implied that this long-held assumption is not valid, and that non-perturbative effects significantly change relevant results obtained by Newtonian gravity. In particular, local physical quantities inside the sun are newly predicted to exhibit power law differently from the so-called standard solar model. This surprising result is reviewed including brief discussion of physics behind the discrepancy and a new application of the new mass formula to gas planets.

Data sizes for next generation radio telescopes, such as the Square Kilometre Array (SKA), are far above that of their predecessors. The CLEAN algorithm was originally developed by Högbom [1974], long before such data sizes were thought possible and is still the most popular tool used for deconvolution in interferometric imaging. In order to facilitate these new large data sizes and reduce computation time a distributed approach to the algorithm has been investigated. The serial nature of the CLEAN algorithm, due to its matching pursuit design, makes this challenging. Splitting the image into a number of tiles which can be individually deconvolved has been investigated, but this creates discontinuities in the deconvolved image and makes it difficult to deconvolve faint sources in the presence of a point spread function associated with bright sources in other tiles. A method of feedback between each of the tiles has been developed to deal with these problems. This new approach has been tested on a simulated dataset containing multiple point sources of known intensity. When compared to a standard Högbom deconvolution the tiled feedback version produced a reconstructed image, containing sources up to 2.1 Jy, which agreed to between -0.1 Jy and +0.04 Jy of the standard method across the whole deconvolved image at a speed up to 10.66 times faster.

Debris disks, which consist of dust, planetesimals, planets, and gas, offer a unique window into the mineralogical composition of their parent bodies, especially during the critical phase of terrestrial planet formation spanning 10 to a few hundred million years. Observations from the $\textit{Spitzer}$ Space Telescope have unveiled thousands of debris disks, yet systematic studies remain scarce, let alone those with unsupervised clustering techniques. This study introduces $\texttt{CLUES}$ (CLustering UnsupErvised with Sequencer), a novel, non-parametric, fully-interpretable machine-learning spectral analysis tool designed to analyze and classify the spectral data of debris disks. $\texttt{CLUES}$ combines multiple unsupervised clustering methods with multi-scale distance measures to discern new groupings and trends, offering insights into compositional diversity and geophysical processes within these disks. Our analysis allows us to explore a vast parameter space in debris disk mineralogy and also offers broader applications in fields such as protoplanetary disks and solar system objects. This paper details the methodology, implementation, and initial results of $\texttt{CLUES}$, setting the stage for more detailed follow-up studies focusing on debris disk mineralogy and demographics.

We present a detailed study of SN 2024ahr, a hydrogen-poor superluminous supernova (SLSN-I), for which we determine a redshift of $z=0.0861$. SN 2024ahr has a peak absolute magnitude of $M_g\approx M_r\approx -21$ mag, rest-frame rise and decline times (50$\%$ of peak) of about 40 and 80 days, respectively, and typical spectroscopic evolution in the optical band. Similarly, modeling of the UV/optical light curves with a magnetar spin-down engine leads to typical parameters: an initial spin period of $\approx 3.3$ ms, a magnetic field strength of $\approx 6\times 10^{13}$ G, and an ejecta mass of $\approx 9.5$ M$_\odot$. Due to its relatively low redshift we obtained a high signal-to-noise ratio near-IR spectrum about 43 rest-frame days post-peak to search for the presence of helium. We do not detect any significant feature at the location of the He I $\,\lambda 2.058$ $\mu$m feature, and place a conservative upper limit of $\sim 0.05$ M$_\odot$ on the mass of helium in the outer ejecta. We detect broad features of Mg I $\,\lambda 1.575$ $\mu$m and a blend of Co II $\,\lambda 2.126$ $\mu$m and Mg II, $\lambda 2.136$ $\mu$m, which are typical of Type Ic SNe, but with higher velocities. Examining the sample of SLSNe-I with NIR spectroscopy, we find that, unlike SN 2024ahr, these events are generally peculiar. This highlights the need for a large sample of prototypical SLSNe-I with NIR spectroscopy to constrain the fraction of progenitors with helium (Ib-like) and without helium (Ic-like) at the time of the explosion, and hence the evolutionary path(s) leading to the rare outcome of SLSNe-I.

Context: NICER (Neutron star Interior Composition ExploreR) is the instrument of choice for the spectral analysis of type I X-ray bursts, as it provides high throughput at X-ray CCD resolution, down to 0.3 keV. Aims: This study investigates whether the energies of absorption lines detected in photospheric radius expansion (PRE) bursts correlate with the inferred blackbody radius. Previous reports suggested such a correlation, attributed to a combination of weaker gravitational redshift and higher blueshifts in bursts with larger radii. Methods: The analysis reexamines four previously studied PRE bursts and examines eight additional bursts from 4U 1820-303, evidencing PRE. Spectral evolution is tracked on the shortest possible timescales (tenth of a second) adopting two parallel continuum descriptions to characterise the photospheric expansion and line evolution. Applying the accretion-enhanced model, maximum blackbody radii of up to $\sim$ 900 km are inferred, with peak bolometric luminosities exceeding the Eddington limit of an Helium accretor. Absorption lines are assessed for significance using Monte Carlo simulations, and spectral lines are characterised using the state-of-art plasma codes available within {\sc{spex}} with a phenomenological continuum. A thorough parameter search explores Doppler shifts to avoid local minima. Results: Several significant (> 99.9%) absorption lines, including the previously reported 2.97 keV line, are detected. While no consistent correlation between line energies and blackbody radii is confirmed, bursts with larger radii exhibit up to four lines and the line strength is higher. The modelling suggests that the observed lines mostly originate from slightly redshifted (almost rest-frame) photo-/collisionally ionised gas in emission. For the burst with the largest PRE, a combination of photo-ionised plasma in both emission and absorption is preferred.

White dwarfs (WDs) are the most abundant compact objects, and recent surveys have suggested that over a third of WDs in accreting binaries host a strong (B $\gtrsim$ 1 MG) magnetic field. However, the origin and evolution of WD magnetism remain under debate. Two WD pulsars, AR Sco and J191213.72-441045.1 (J1912), have been found, which are non-accreting binaries hosting rapidly spinning (1.97-min and 5.30-min, respectively) magnetic WDs. The WD in AR Sco is slowing down on a $P/\dot{P}\approx 5.6\times 10^6$ yr timescale. It is believed they will eventually become polars, accreting systems in which a magnetic WD (B $\approx 10-240$ MG) accretes from a Roche lobe-filling donor spinning in sync with the orbit ($\gtrsim 78$ min). Here, we present multiwavelength data and analysis of Gaia22ayj, which outbursted in March 2022. We find that Gaia22ayj is a magnetic accreting WD that is rapidly spinning down ($P/\dot{P} = 6.1^{+0.3}_{-0.2}\times 10^6$ yr) like WD pulsars, but shows clear evidence of accretion, like polars. Strong linear polarization (40%) is detected in Gaia22ayj; such high levels have only been seen in the WD pulsar AR Sco and demonstrate the WD is magnetic. High speed photometry reveals a 9.36-min period accompanying a high amplitude ($\sim 2$ mag) modulation. We associate this with a WD spin or spin-orbit beat period, not an orbital period as was previously suggested. Fast (60-s) optical spectroscopy reveals a broad ``hump'', reminiscent of cyclotron emission in polars, between 4000-8000 Angstrom. We find an X-ray luminosity of $L_X = 2.7_{-0.8}^{+6.2}\times10^{32} \textrm{ erg s}^{-1}$ in the 0.3-8 keV energy range, while two VLA radio campaigns resulted in a non-detection with a $F_r < 15.8\mu\textrm{Jy}$ 3$ \sigma$ upper limit. The shared properties of both WD pulsars and polars suggest that Gaia22ayj is a missing link between the two classes of magnetic WD binaries.

Holmberg 15A (H15A), the brightest cluster galaxy of Abell 85, has an exceptionally low central surface brightness even among local massive elliptical galaxies with distinct stellar cores, making it exceedingly challenging to obtain high-quality spectroscopy to detect a supermassive black hole (SMBH) at its center. Aided by the superb sensitivity and efficiency of KCWI at the Keck II Telescope, we have obtained spatially resolved stellar kinematics over a ${\sim}100''\times 100''$ contiguous field of H15A for this purpose. The velocity field exhibits a low amplitude (${\sim}20\mathrm{~km~s}^{-1}$) rotation along a kinematic axis that is prominently misaligned from the photometric major axis, a strong indicator that H15A is triaxially shaped with unequal lengths for the three principal axes. Using 2500 observed kinematic constraints, we perform extensive calculations of stellar orbits with the triaxial Schwarzschild code, TriOS, and search over ${\sim}$40,000 galaxy models to simultaneously determine the mass and intrinsic 3D shape parameters of H15A. We determine a ratio of $p=0.89$ for the middle-to-long principal axes and $q=0.65$ for the short-to-long principal axes. Our best estimate of the SMBH mass, $M_\mathrm{BH}=(2.16^{+0.23}_{-0.18})\times 10^{10}M_{\odot}$, makes H15A -- along with NGC 4889 -- the galaxies hosting the most massive SMBHs known in the local universe. Both SMBHs lie significantly above the mean $M_\mathrm{BH}-\sigma$ scaling relation. Repeating the orbit modeling with the axisymmetrized version of TriOS produces worse fits to the KCWI kinematics and increases $M_\mathrm{BH}$ to $(2.55\pm 0.20) \times 10^{10}M_{\odot}$, which is still significantly below $M_\mathrm{BH}=(4.0\pm 0.8) \times 10^{10}M_{\odot}$ reported in a prior axisymmetric study of H15A.

Continuous gravitational waves (CWs) emission from neutron stars carries information about their internal structure and equation of state, and it can provide tests of General Relativity. We present a search for CWs from a set of 45 known pulsars in the first part of the fourth LIGO--Virgo--KAGRA observing run, known as O4a. We conducted a targeted search for each pulsar using three independent analysis methods considering the single-harmonic and the dual-harmonic emission models. We find no evidence of a CW signal in O4a data for both models and set upper limits on the signal amplitude and on the ellipticity, which quantifies the asymmetry in the neutron star mass distribution. For the single-harmonic emission model, 29 targets have the upper limit on the amplitude below the theoretical spin-down limit. The lowest upper limit on the amplitude is $6.4\!\times\!10^{-27}$ for the young energetic pulsar J0537-6910, while the lowest constraint on the ellipticity is $8.8\!\times\!10^{-9}$ for the bright nearby millisecond pulsar J0437-4715. Additionally, for a subset of 16 targets we performed a narrowband search that is more robust regarding the emission model, with no evidence of a signal. We also found no evidence of non-standard polarizations as predicted by the Brans-Dicke theory.

We present ORACLE, the first hierarchical deep-learning model for real-time, context-aware classification of transient and variable astrophysical phenomena. ORACLE is a recurrent neural network with Gated Recurrent Units (GRUs), and has been trained using a custom hierarchical cross-entropy loss function to provide high-confidence classifications along an observationally-driven taxonomy with as little as a single photometric observation. Contextual information for each object, including host galaxy photometric redshift, offset, ellipticity and brightness, is concatenated to the light curve embedding and used to make a final prediction. Training on $\sim$0.5M events from the Extended LSST Astronomical Time-Series Classification Challenge, we achieve a top-level (Transient vs Variable) macro-averaged precision of 0.96 using only 1 day of photometric observations after the first detection in addition to contextual information, for each event; this increases to $>$0.99 once 64 days of the light curve has been obtained, and 0.83 at 1024 days after first detection for 19-way classification (including supernova sub-types, active galactic nuclei, variable stars, microlensing events, and kilonovae). We also compare ORACLE with other state-of-the-art classifiers and report comparable performance for the 19-way classification task, in addition to delivering accurate top-level classifications much earlier. The code and model weights used in this work are publicly available at our associated GitHub repository (this https URL).

Common features of sub-Neptunes atmospheres observed to date include signatures of aerosols at moderate equilibrium temperatures (~500-800 K), and a prevalence of high mean molecular weight atmospheres, perhaps indicating novel classes of planets such as water worlds. Here we present a 0.83-5 micron JWST transmission spectrum of the sub-Neptune TOI-421 b. This planet is unique among previously observed counterparts in its high equilibrium temperature ($T_{eq} \approx 920$) and its Sun-like host star. We find marked differences between the atmosphere of TOI-421 b and those of sub-Neptunes previously characterized with JWST, which all orbit M stars. Specifically, water features in the NIRISS/SOSS bandpass indicate a low mean molecular weight atmosphere consistent with solar metallicity, and no appreciable aerosol coverage. Hints of SO$_2$ and CO (but not CO$_2$ or CH$_4$) also exist in our NIRSpec/G395M observations, but not at sufficient signal-to-noise to draw firm conclusions. Our results support a picture in which sub-Neptunes hotter than ~850 K do not form hydrocarbon hazes due to a lack of methane to photolyze. TOI-421 b additionally fits the paradigm of the radius valley for planets orbiting FGK stars being sculpted by mass loss processes, which would leave behind primordial atmospheres overlying rock/iron interiors. Further observations of TOI-421 b and similar hot sub-Neptunes will confirm whether haze-free atmospheres and low mean molecular weights are universal characteristics of such objects.

The isotropic stochastic gravitational wave background (SGWB) generated by a population of supermassive black hole binaries (BHBs) provides a unique window into their cosmic evolution. In addition to the isotropic power spectrum, the anisotropic component of the signal carries additional information about the supermassive black holes (BHs) and host galaxy connection. The measurement of this signal is usually carried out by angular power spectra, which is only a sufficient measure for a Gaussian and statistically isotropic distribution. In contrast, the contribution from supermassive BHBs in nano-hertz SGWB will be hosted by fewer massive galaxies, making the nano-hertz background anisotropic and non-Gaussian. As a result, the performance of angular power spectra in extracting the underlying physics is limited. In this work, we propose a novel technique called the Multi-Tracer Correlated Stacking, which enables the detection of anisotropies in the SGWB by stacking the signal from regions of the sky with tracers of BHs such as active galactic nucleus (AGNs), quasars, bright galaxies, etc., that can be mapped up to high redshift. We demonstrate this technique on a simulated supermassive BHBs distribution using an AGN catalog, which maps the underlying matter distribution approximately up to redshift $z=5$. This stacking technique uniquely distinguishes between isotropic and anisotropic distributions of SGWB source, surpassing the capabilities of angular power spectrum-based methods in detecting anisotropic signals. This highlights the effectiveness of this technique in detecting anisotropic SGWB signals and in the future, this technique can play a crucial role in its discovery.

Binospec is a wide-field optical (360 to 1000 nm) spectrograph commissioned at the MMT 6.5m telescope in 2017. In direct imaging mode Binospec addresses twin 8$^\prime$ (wide) by 15$^\prime$ (slit length) fields of view. We describe an optical fiber based integral field unit (IFU) that remaps a 12$^{\prime\prime}$ x 16$^{\prime\prime}$ contiguous region onto two pseudo slits, one in each Binospec channel. The IFU, commissioned in 2023, fits into the space of a standard slit mask frame and can be deployed as desired in a mixed program of slit masks, long slits, and IFU observations. The IFU fibers are illuminated by a hexagonal lenslet array with a 0.6$^{\prime\prime}$ pitch. A separate bundle of sky fibers consists of close-packed bare fibers arranged within an 11.8$^{\prime\prime}$ circular aperture. The 640 IFU fibers and 80 sky fibers have a core diameter of 150$\mu$m, corresponding to 0.90$^{\prime\prime}$. Three gratings are available, 270lpm with R$\sim$2000, 600lpm with R$\sim$5300, and 1000 lpm with R$\sim$6000.

We present a comprehensive photometric and spectroscopic study of the Type IIP SN 2018is. The $V$-band luminosity and the expansion velocity at 50 days post-explosion are $-$15.1$\pm$0.2 mag (corrected for A$_V$=1.34 mag) and 1400 km s$^{-1}$, classifying it as a low-luminosity SN II. The recombination phase in the $V$-band is shorter, lasting around 110 days, and exhibits a steeper decline (1.0 mag per 100 days) compared to most other low-luminosity SNe II. Additionally, the optical and near-infrared spectra display hydrogen emission lines that are strikingly narrow, even for this class. The Fe II and Sc II line velocities are at the lower end of the typical range for low-luminosity SNe II. Semi-analytical modelling of the bolometric light curve suggests an ejecta mass of $\sim$8 M$_\odot$, corresponding to a pre-supernova mass of $\sim$9.5 M$_\odot$, and an explosion energy of $\sim$0.40 $\times$ 10$^{51}$ erg. Hydrodynamical modelling further indicates that the progenitor had a zero-age main sequence mass of 9 M$_\odot$, coupled with a low explosion energy of 0.19 $\times$ 10$^{51}$ erg. The nebular spectrum reveals weak [O I] $\lambda\lambda$6300,6364 lines, consistent with a moderate-mass progenitor, while features typical of Fe core-collapse events, such as He I, [C I], and [Fe I], are indiscernible. However, the redder colours and low ratio of Ni to Fe abundance do not support an electron-capture scenario either. As a low-luminosity SN II with an atypically steep decline during the photospheric phase and remarkably narrow emission lines, SN 2018is contributes to the diversity observed within this population.

Galactic black hole (BH) X-ray binaries are known to exhibit episodic outbursts, during which accretion and spectral mode distinctively transition between low/hard and high/soft state. X-ray observations during high/soft state occasionally reveal a pronounced presence of a powerful disk wind in these systems. However, it is unexplored to date how such winds may influence disk emission in that regime. In this work, we consider an observational implication by Compton scattering of thermal disk radiation due to accretion disk winds by performing multi-dimensional Monte Carlo simulations in the context of a stratified wind of large solid angle launched over a large radial extent of the disk. Compton-scattered thermal disk spectrum is computed for a different wind property; i.e. wind density and its radial gradient. We find that the intrinsic disk radiation can be significantly down-scattered by winds of moderate-to-high density to the extent that the transmitted spectrum can substantially deviate from the conventional multi-color-disk emission by factor of up to 70\%. We thus claim that the conventional treatment of spectral hardening in the disk atmosphere may be insufficient to fully account for the observed disk continuum in the presence of strong wind scattering. It is suggested that the effect of scattering process (by $f_w$) should be incorporated to accurately evaluate an intrinsic disk spectrum besides the conventional hardening (color correction) factor (by $f_c$). We argue that BH spin measurements using thermal continuum-fitting in transient XRBs may well be mildly (if not significantly) altered by such spectral ``contamination".

The enduring technique of aperture masking interferometry, now more than 150 years old, is still widely practised today for it opens a window of high angular resolution astronomy that remains difficult to access by any competing technology. However, the requirement to apodise the pupil into a non-redundant array dramatically limits the throughput, typically to $\sim$10\% or less. This in turn has a dramatic impact on sensitivity, limiting observational reach to only bright science targets. This paper presents "Jewel Optics", a novel technology that leverages the gains in signal fidelity conferred by non-redundant Fizeau beam combination without the sensitivity penalty incurred by traditional aperture masks. Our approach fragments the pupil into several sets of sparse-array non-redundant patterns, each of which is encoded onto a unique phase wedge. After extensive searching, solutions could be found where all individual sets are fully non-redundant while fully tiling the available area of the input pupil. Each pattern is assigned a common phase wedge which diverts light from those sub-apertures onto a unique, defined region of the detector. We demonstrate a prototype designed for use in the VAMPIRES instrument at the Subaru telescope and find excellent agreement between the design and lab results. We discuss a design refinement for producing fully achromatic Jewel Optics, and finally we highlight the potential for future work with these novel optical components.

To date, only two strongly lensed type Ia supernovae (SNIa) have been discovered with an isolated galaxy acting as the lens: iPTF16geu and SN Zwicky. The observed image fluxes for both lens systems were inconsistent with predictions from a smooth macro lens model. A potential explanation for the anomalous flux ratios is microlensing: additional (de)magnification caused by stars and other compact objects in the lens galaxy. In this work, we combine observations of iPTF16geu and SN Zwicky with simulated microlensing magnification maps, leveraging their standardizable candle properties to constrain the lens galaxy mass slope, $\eta$, and the fraction of dark compact objects, $f_{\rm dc}$. The resulting mass slopes are $\eta = 1.70 \pm 0.07$ for iPTF16geu and $\eta = 1.81 \pm 0.10$ for SN Zwicky. Our results indicate no evidence for a population of dark compact objects, placing upper limits at the $95\%$ confidence level of $f_{\rm dc} < 0.25$ for iPTF16geu and $f_{\rm dc} < 0.47$ for SN Zwicky. Assuming a constant fraction of dark compact objects for both lensed SNe, we obtain $f_{\rm dc} < 0.19$. These results highlight the potential of strongly lensed SNIa to probe the innermost parts of lens galaxies and learn about compact matter.

Recent discoveries from time-domain surveys are defying our expectations for how matter accretes onto supermassive black holes (SMBHs). The increased rate of short-timescale, repetitive events around SMBHs, including the newly-discovered quasi-periodic eruptions (QPEs), are garnering further interest in stellar-mass companions around SMBHs and the progenitors to mHz frequency gravitational wave events. Here we report the discovery of a highly significant mHz Quasi-Periodic Oscillation (QPO) in an actively accreting SMBH, 1ES 1927+654, which underwent a major optical, UV, and X-ray outburst beginning in 2018. The QPO was first detected in 2022 with a roughly 18-minute period, corresponding to coherent motion on scales of less than 10 gravitational radii, much closer to the SMBH than typical QPEs. The period decreased to 7.1 minutes over two years with a decelerating period evolution ($\ddot{P} > 0$). This evolution has never been seen in SMBH QPOs or high-frequency QPOs in stellar mass black holes. Models invoking orbital decay of a stellar-mass companion struggle to explain the period evolution without stable mass transfer to offset angular momentum losses, while the lack of a direct analog to stellar mass black hole QPOs means that many instability models cannot explain all of the observed properties of the QPO in 1ES 1927+654. Future X-ray monitoring will test these models, and if it is a stellar-mass orbiter, the Laser Interferometer Space Antenna (LISA) should detect its low-frequency gravitational wave emission.

The cold and hot interstellar medium (ISM) in star forming galaxies resembles the reservoir for star formation and associated heating by stellar winds and explosions during stellar evolution, respectively. We utilize data from deep $Chandra$ observations and archival millimeter surveys to study the interconnection between these two phases and the relation to star formation activities in M51 on kiloparsec scales. A sharp radial decrease is present in the hot gas surface brightness profile within the inner 2 kpc of M51. The ratio between the total infrared luminosity ($L_{\rm IR}$) and the hot gas luminosity ($L_{\rm 0.5 - 2\,keV}^{\rm gas}$) shows a positive correlation with the galactic radius in the central region. For the entire galaxy, a twofold correlation is revealed in the $L_{\rm 0.5 - 2\,keV}^{\rm gas}$${-}$$L_{\rm IR}$ diagram, where $L_{\rm 0.5 - 2\,keV}^{\rm gas}$ sharply increases with $L_{\rm IR}$ in the center but varies more slowly in the disk. The best fit gives a steep relation of ${\rm log}(L_{\rm 0.5-2\,keV}^{\rm gas} /{\rm erg\,s^{-1}})=1.82\,{\rm log}(L_{\rm IR} /{L_{\rm \odot}})+22.26$ for the center of M51. The similar twofold correlations are also found in the $L_{\rm 0.5 - 2\,keV}^{\rm gas}$${-}$molecular line luminosity ($L^\prime_{\rm gas}$) relations for the four molecular emission lines CO(1-0), CO(2-1), HCN(1-0), and HCO$^+$(1-0). We demonstrate that the core-collapse supernovae (SNe) are the primary source of energy for heating gas in the galactic center of M51, leading to the observed steep $L_{\rm 0.5 - 2\,keV}^{\rm gas}$${-}$$L_{\rm IR}$ and $L_{\rm 0.5 - 2\,keV}^{\rm gas}$${-}$$L^\prime_{\rm gas}$ relations, as their X-ray radiation efficiencies ($\eta$ $\equiv$ $L_{\rm 0.5 - 2\,keV}^{\rm gas}$/$\dot{E}_\mathrm{SN}$) increase with the star formation rate surface densities, where $\dot{E}_\mathrm{SN}$ is the SN mechanical energy input rate.

Accurate determination of binary fractions ($f_{\rm b}$) and mass ratio ($q$) distributions is crucial for understanding the dynamical evolution of open clusters. We present an improved multiband fitting technique to enhance the analysis of binary properties. This approach enables an accurate photometric determination of $f_{\rm b}$ and $q$ distribution in a cluster. The detectable mass ratio can be down to the $q_{\rm lim}$, limited by the minimum stellar mass in theoretical models. First, we derived an empirical model for magnitudes of Gaia DR3 and 2MASS bands that match the photometry of single stars in the Pleiades. We then performed a multiband fitting for each cluster member, deriving the probability density function (PDF) of its primary mass ($\mathcal{M}_1$) and $q$ in the Bayesian framework. 1154 main-sequence (MS) single stars or unresolved MS+MS binaries are identified as members of the Pleiades. By stacking their PDFs, we conducted a detailed analysis of binary properties of the cluster. We found the $f_{\rm b}$ of this sample is $0.34 \pm 0.02$. The $q$ distribution exhibits a three-segment power-law profile: an initial increase, followed by a decrease, and then another increase. This distribution can be interpreted as a fiducial power-law profile with an exponent of -1.0 that is determined in the range of $0.3 < q < 0.8$, but with a deficiency of binaries at lower $q$ and an excess at higher $q$. The variations of $f_{\rm b}$ and $q$ with $\mathcal{M}_1$ reveal a complex binary distribution within the Pleiades, which might be attributed to a combination of primordial binary formation mechanisms, dynamical interactions, and the observational limit of photometric binaries imposed by $q_{\rm lim} (\mathcal{M}_1)$.

For the first time, we develop a simulation-based model for the Minkowski functionals (MFs) of large-scale structure, which allows us to extract the full information available from the MFs (including both the Gaussian and non-Gaussian part), and apply it to the BOSS DR12 CMASS galaxy sample. Our model is based on high-fidelity mock galaxy catalogs constructed from the \textsc{Abacus}\textsc{Summit} simulations using the halo occupation distribution (HOD) framework, which include the redshift-space distortions and Alcock-Paczynski distortions, incorporate survey realism, including survey geometry and veto masks, and account for angular plus radial selection effects. The cosmological and HOD parameter dependence of the MFs is captured with a neural network emulator trained from the galaxy mocks with various cosmological and HOD parameters. To benchmark the constraining power of the MFs, we also train an emulator for the galaxy 2-point correlation function (2PCF) using the same pipeline. Having validated our approach through successful parameter recovery tests on both internal and external mocks, including non-HOD forward models of the halo-galaxy connection, we apply our forward model to analyze the CMASS data in the redshift range $0.45<z<0.58$. We find the MFs provide stronger constraints on the cosmological parameters than the 2PCF. The combination of the two gives $\omega_{\rm cdm}=0.1172^{+0.0020}_{-0.0023}$, $\sigma_8=0.783\pm 0.026$, and $n_s=0.966^{+0.019}_{-0.015}$, which are tighter by a factor of 2.0, 1.9, and 1.6 than the 2PCF alone. The derived constraint $f\sigma_8=0.453 \pm 0.016$ is also improved by a factor of 1.9, compared to the 2PCF, and agrees well with Planck 2018 predictions and other results from a series of studies in the literature.

The Pierre Auger Observatory concluded its first phase of data taking after seventeen years of operation. The dataset collected by its surface and fluorescence detectors (FD and SD) provides us with the most precise estimates of the energy spectrum and mass composition of ultra-high energy cosmic rays yet available. We present measurements of the depth of shower maximum, the main quantity used to derive species of primary particles, determined either from the direct observation of longitudinal profiles of showers by the FD, or indirectly through the analysis of signals in the SD stations. The energy spectrum of primaries is also determined from both FD and SD measurements, where the former exhibits lower systematic uncertainty in the energy determination while the latter exploits unprecedentedly large exposure. The data for primaries with energy below 1 EeV are also available thanks to the high-elevation telescopes of FD and the denser array of SD, making measurements possible down to 6 PeV and 60 PeV, respectively.

Classical Cepheids are a cornerstone class of pulsators, fundamental to testing stellar evolution and pulsation theories. Their secular period changes, characterized through $O-C$ (Observed minus Calculated) diagrams, offer valuable insights into their evolution. While evolutionary period changes are well understood from both observational and theoretical perspectives, shorter timescale period changes (on the order of ($\sim$ 10$^{2}$-10$^{4}$ days) - known as non-evolutionary period changes are yet to be systematically explored. In this work, we present a detailed and comprehensive search for non-evolutionary period changes using $O-C$ analysis of Magellanic Cloud (MC) Cepheids, based on 20+ years of OGLE photometry data. Our sample includes both the Large Magellanic Cloud (LMC) and the Small Magellanic Cloud (SMC) Cepheids, focusing on single radial mode Cepheids (both fundamental (FU) and first overtone (FO) modes). The results are grouped into two phenomena: (a) Cepheids in binary systems (b) Non-linear period changes.

Gas and dust in the Galactic Center are subjected to energetic processing by intense UV radiation fields, widespread shocks, enhanced rates of cosmic-rays and X-rays, and strong magnetic fields. The Giant Molecular Clouds in the Galactic Center present a rich chemistry in a wide variety of chemical compounds, some of which are prebiotic. We have conducted unbiased, ultrasensitive and broadband spectral surveys toward the G+0.693-0.027 molecular cloud located in the Galactic Center, which have yielded the discovery of new complex organic molecules proposed as precursors of the "building blocks" of life. I will review our current understanding of the chemistry in Galactic Center molecular clouds, and summarize the recent detections toward G+0.693-0.027 of key precursors of prebiotic chemistry. All this suggests that the ISM is an important source of prebiotic material that could have contributed to the process of the origin of life on Earth and elsewhere in the Universe.

The carbon isotope ratio ($^{12}\mathrm{C}/^{13}\mathrm{C}$) is a potential tracer of giant planet and brown dwarf formation. The GQ Lup system, hosting the K7 T Tauri star GQ Lup A and its substellar companion GQ Lup B, offers a unique opportunity to investigate this ratio in a young system. We aim to characterise the atmosphere of GQ Lup B, determining its temperature, chemical composition, spin, surface gravity, and $^{12}\mathrm{C}/^{13}\mathrm{C}$, while also measuring the same ratio for its host star, GQ Lup A. High-resolution K-band spectra of GQ Lup were obtained using CRIRES+ at the VLT. We modelled the starlight contribution from GQ Lup A and fitted GQ Lup Bś spectrum with atmospheric models from petitRADTRANS. The $^{12}\mathrm{C}/^{13}\mathrm{C}$ ratio for GQ Lup A was derived using isotope-sensitive PHOENIX models. Atmospheric analysis of GQ Lup B revealed abundances of H$_2$O, $^{12}$CO, $^{13}$CO, HF, Na, Ca, and Ti, along with a C/O ratio of $0.50 \pm 0.01$, consistent with the solar value. The carbon isotope ratio was measured as $^{12}\mathrm{C}/^{13}\mathrm{C} = 53^{+7}_{-6}$ for GQ Lup B and $^{12}\mathrm{C}/^{13}\mathrm{C} = 51^{+10}_{-8}$ for GQ Lup A. Strong veiling of the photospheric lines of GQ Lup A was identified and accounted for. The similar $^{12}\mathrm{C}/^{13}\mathrm{C}$ values in GQ Lup A and B suggest a common origin from a shared material reservoir, supporting formation via disc fragmentation or gravitational collapse.

A black hole (BH) can tear apart a star that ventures within its tidal radius, producing a luminous flare as the stellar debris falls back, known as a tidal disruption event (TDE). While TDEs in quiescent galaxies are relatively well understood, identifying TDEs in active galactic nuclei (AGN) still remains a significant challenge. We present the discovery of AT2021aeuk, a transient exhibiting dual flares within around three years in a narrow-line Seyfert 1 galaxy. Multi-wavelength observations triggered during the second flare in 2023 revealed an extraordinary X-ray V-shaped light curve, strongly anti-correlated with the optical light curve and accompanied by a lag of $\sim$40 days. This behavior is inconsistent with both supernova and pure AGN origins. In addition, a new broad component emerges in the Balmer lines during the second flare, showing a clear reverberation signal to the continuum variation. We propose that the dual-flare may be linked to a repeating partial TDE (rpTDE), where the second flare results from a collision between the TDE stream and the inner accretion disk, triggering an optical flare while simultaneously partially destroying the X-ray corona. However, other mechanisms, such as a stellar-mass BH (sBH) merger within an accretion disk, could produce similar phenomena, which we cannot entirely rule out. The Vera C. Rubin Observatory will be a powerful tool for further investigating the nature of such events in the future.

Multi-messenger astronomy is key today to broaden our understanding of the high energy Universe. When an ultra-high energy (UHE) particle interacts in the atmosphere or underground, it initiates an extensive air shower that produces a coherent radio emission in the 50-200 MHz range. The GRAND project is an envisioned observatory of UHE particles (cosmic rays, gamma rays and especially neutrinos) that will consist of 200,000 radio antennas deployed over 20 different locations each of $\sim$ 10,000 $\rm km^2$. In its current phase, it consists of 3 prototypes running autonomously in 3 different locations: GRAND@Auger in Argentina, GRAND@Nançay in France and GRANDProto300 in China, all at commissioning stage. The first goal of these pathfinders is to demonstrate the viability of the GRAND detection concept. GRANDProto300 will also propose a rich science case, by allowing the study of the transition between galactic and extra-galactic cosmic ray sources. In the following, we present the detection concept, the preliminary designs and layout and the ongoing developments of the experiment.

Gamma-ray bursts are known to display various features on top of their canonical behavior. In this short review, we will describe and discuss two of them: the ultra-long gamma-ray bursts, which are defined by an extreme duration of their prompt phase, and the plateau phase, which is defined by a steady phase of large duration at the start of the afterglow. We will review the main properties of those two phenomena, and will discuss their possible origin, in light of the standard fireball model of gamma-ray bursts. A final section will discuss the future missions which could bring new evidences to the study of those objects.

The study of pulsar glitch phenomena serves as a valuable probe into the dynamic properties of matter under extreme high-density conditions, offering insights into the physics within neutron stars. Providing theoretical explanations for the diverse manifestations observed in different pulsars has proven to be a formidable challenge. By analyzing the distribution of glitch sizes and waiting times, along with the evolution of cumulative glitch sizes over time, we have uncovered a long-term clustering phenomenon for pulsar glitches. This perspective allows us to approach the distinct glitch representations in various pulsars from a unified standpoint, connecting the same periodicity of observational data to the randomness. Without relying on specific physical models, we utilized the coefficient of variation to numerically determine optimal clustering numbers and clustering periods for sample pulsars. Our analysis involving 27 pulsars has revealed a clear linear relationship between the glitch cluster period and characteristic age. Of interest, the cumulative distribution of functions of cluster sizes and interval times have the same patterns, which can be synchronously fitted by Gaussian processes. These results may indicate novel understandings of glitches and the resulting processes.

Recent solar observations of bipolar light bridges (BLBs) in sunspots have, in a few individual cases, revealed magnetic fields up to 8.2 kG, which is at least twice as strong as typical values measured in sunspot umbrae. However, the small number of such observations hinted that such strong fields in these bright photospheric features that separate two opposite-polarity umbrae, are a rare phenomenon. We determine the field strength in a large sample of BLBs with the aim of establishing how prevalent such strong fields are in BLBs. We apply a state-of-the-art inversion technique that accounts for the degradation of the data by the intrinsic point spread function of the telescope, to the so far largest set of spectropolarimetric observations, by Hinode/Solar Optical Telescope spectropolarimeter, of sunspots containing BLBs. We identified 98 individual BLBs within 51 distinct sunspot groups. Since 66.3% of the BLBs were observed multiple times, a total of 630 spectropolarimetric scans of these 98 BLBs were analysed. All analysed BLBs contain magnetic fields stronger than 4.5 kG at unit optical depth. The field strengths decrease faster with height than the fields in umbrae and penumbrae. BLBs display a unique continuum intensity and field strength combination, forming a population well separated from umbrae and the penumbrae. The high brightness of BLBs in spite of their very strong magnetic fields points to the presence of a so far largely unexplored regime of magnetoconvection.

We observed an optical afterglow of GRB 200131A obtaining the first photometric point 63 s after the satellite trigger. This early observation shows a steep decay, suggesting either internal engine activity or a reverse shock. By fitting this data set, we show that the early data fit well as a reverse shock component of the GRB afterglow modeled as a thin shell expanding into a constant density interstellar matter. The fitting also shows a good agreement with a catalogued Milky Way galactic extinction and leaves only little space for further extinction in the host galaxy. By judging several factors we conclude that the most likely redshift of this GRB is 0.9 +/- 0.1.

Conventional general relativity supplies the notion of a vacuum tension and thus a maximum force $F_{max}=c^4/4G\approx\ 3\times 10^{43}$ Newtons, that is realized for a black hole. In conjunction with the Wilson area rule, we are thus led to the surface confinement of the mass of a black hole analogous to the surface confinement of quarks. The central result of our paper is that PeV scale protons exist on the surface of a Kerr black hole residing at our galactic center that is in concert with the HAWC Collaboration result of a PeVatron at the galactic center.

Massive galaxies in cooling flow clusters display clear evidence of feedback from Active Galactic Nuclei (AGN). Joint X-ray and radio observations have shown that AGN radio jets push aside the surrounding hot gas and form cavities in the hot intracluster medium (ICM). These systems host complex, kiloparsec-scale, multiphase filamentary structures, from warm ionized (10,000 K) to cold molecular ($<$100 K). These striking clumpy filaments are believed to be a natural outcome of thermally unstable cooling from the hot ICM, likely triggered by feedback processes while contributing to feeding the AGN via Chaotic Cold Accretion (CCA). However, the detailed constraints on the formation mechanism of the filaments are still uncertain, and the connection between the different gas phases has to be fully unveiled. By leveraging a sample of seven X-ray bright cooling-flow clusters, we have discovered a tight positive correlation between the X-ray surface brightness and the H$\alpha$ surface brightness of the filaments over two orders of magnitude, as also found in stripped tails.

We have analyzed archival far-ultraviolet spectra of the UV-bright star III-60 in the globular cluster NGC 6723 obtained with the Far Ultraviolet Spectroscopic Explorer (FUSE) and the Cosmic Origins Spectrograph (COS). We find that the star's photospheric parameters (effective temperature $T_{\rm eff} = 44{,}800 \pm 1200$, surface gravity $\log g = 4.89 \pm 0.18$, and helium abundance $\log N({\rm He})/N({\rm H}) = -0.84 \pm 0.29$) are consistent with the values derived from its optical spectrum, suggesting that optically-derived values are generally accurate for evolved stars with $T_{\rm eff} \lesssim$ 50,000 K. Relative to the cluster's RGB stars, III-60 is enhanced in nitrogen and depleted in carbon and oxygen. The star exhibits strong P Cygni profiles in both components of the N V $\lambda 1240$ doublet, but the resonance lines of other species show no evidence of a stellar wind. The star's effective temperature and luminosity place it on the evolutionary tracks of stars evolving from the blue horizontal branch, but its high mass ($\sim 1.2 \, M_{\odot}$) indicates that it is the product of a stellar merger. Its helium, carbon, and nitrogen abundances suggest that it is following an evolutionary path similar to that of the low-carbon, intermediate helium-rich hot subdwarfs.

Direct imaging of exoplanets is crucial for advancing our understanding of planetary systems beyond our solar system, but it faces significant challenges due to the high contrast between host stars and their planets. Wavefront aberrations introduce speckles in the telescope science images, which are patterns of diffracted starlight that can mimic the appearance of planets, complicating the detection of faint exoplanet signals. Traditional post-processing methods, operating primarily in the image intensity domain, do not integrate wavefront sensing data. These data, measured mainly for adaptive optics corrections, have been overlooked as a potential resource for post-processing, partly due to the challenge of the evolving nature of wavefront aberrations. In this paper, we present a differentiable rendering approach that leverages these wavefront sensing data to improve exoplanet detection. Our differentiable renderer models wave-based light propagation through a coronagraphic telescope system, allowing gradient-based optimization to significantly improve starlight subtraction and increase sensitivity to faint exoplanets. Simulation experiments based on the James Webb Space Telescope configuration demonstrate the effectiveness of our approach, achieving substantial improvements in contrast and planet detection limits. Our results showcase how the computational advancements enabled by differentiable rendering can revitalize previously underexploited wavefront data, opening new avenues for enhancing exoplanet imaging and characterization.

The Lyman-alpha (LyA) line of atomic hydrogen encodes crucial information about both the intrinsic sources and surrounding environments of star-forming regions throughout the Universe. Due to the complexity of resonant scattering, analytic solutions remain scarce, with most studies focusing on idealized, static configurations. However, observations of LyA emitting galaxies consistently reveal signatures of outflows, imprinted through red-peak dominance in spectral line profiles. In this paper, we derive novel analytic solutions for resonant-line radiative transfer in moving media, specifically homologous-like cloud expansion and unbounded cosmological flows, which capture the main physics of velocity gradients. To validate these analytic solutions and identify regimes where diffusion-based assumptions hold, we introduce a robust Gridless Monte Carlo Radiative Transfer (GMCRT) method. By integrating optical depths exactly in the comoving frame, GMCRT updates photon frequencies continuously to account for Doppler shifts induced by velocity gradients. We demonstrate excellent consistency between GMCRT and our analytic solutions in regimes where diffusion approximations apply. At higher velocities or lower optical depths, discrepancies highlight the limitations of simplified formalisms. We also provide scaling relations for a point source in a cloud with a maximum-to-thermal velocity ratio beta = V_max / v_th, as modifying the standard dependence on line centre optical depth of (atau_0)^1/3 by additional factors, e.g. characteristic escape frequency scale as x_esc ~ beta^1/3, force multipliers as M_F ~ beta^-1/3, and trapping time as t_trap ~ beta^-2/3. Our work complements numerical simulations by improving physical intuition about nonstatic environments when interpreting LyA observations and guiding future subgrid prescriptions in galaxy formation models.

We report the first detection in interstellar medium of the 1-cyano propargyl radical, HC$_3$HCN. This species is an isomer of the 3-cyano propargyl radical (CH$_2$C$_3$N), which was recently discovered in TMC-1. The 1-cyano propargyl radical was observed in the cold dark cloud TMC-1 using data from the ongoing QUIJOTE line survey, which is being carried out with the Yebes 40m telescope. A total of seven rotational transitions with multiple hyperfine components were detected in the 31.0-50.4 GHz range. We derived a column density of (2.2$\pm$0.2)$\times$10$^{11}$ cm$^{-2}$ and a rotational temperature of 7$\pm$1\,K. The abundance ratio between HC$_3$HCN and CH$_2$C$_3$N is 1.4. The almost equal abundance of these isomers indicates that the two species may be produced in the same reaction with a similar efficiency, probably in the reaction C + CH$_2$CHCN and perhaps also in the reaction C$_2$ + CH$_3$CN and the dissociative recombination with electrons of CH$_2$C$_3$NH$^+$

We report the results of long-term time series photometry on RX J2133.7+5107 (also known as 1RXS J213344.1+510725) obtained at several observatories. Using data taken during 17 years, we determined the current value of the spin period of $570.811470$ seconds with the formal accuracy of $0.000006$ seconds and a spin-up of the white dwarf with a characteristic time of $1.483(1)\times10^5$ years. This is even faster than that reported previously and, if confirmed, makes this object have one of the fastest spin-up timescales of all known intermediate polars. We derived an improved value of the superhump period of the system to be $0^d.280130(1)$. Superhump maxima timings are moving on the phase curve from season to season, showing non-monotonic changes, without a change in superhump period.

We employ deep learning (DL) to improve photometric redshifts (photo-$z$s) in the Kilo-Degree Survey Data Release 4 Bright galaxy sample (KiDS-Bright DR4). This dataset, used as a foreground for KiDS lensing and clustering studies, is flux-limited to $r<20$ mag with mean $z=0.23$ and covers 1000 deg$^2$. Its photo-$z$s were previously derived with artificial neural networks from the ANNz2 package, trained on the Galaxy And Mass Assembly (GAMA) spectroscopy. Here we considerably improve over these previous redshift estimations by building a DL model, Hybrid-z, which combines four-band KiDS images with nine-band magnitudes from KiDS+VIKING. The Hybrid-z framework provides photo-$z$s for KiDS-Bright, with negligible mean residuals of O($10^{-4}$) and scatter at the level of $0.014(1+z)$ -- reduction by 20% over the previous nine-band derivations with ANNz2. We check our photo-$z$ model performance on test data drawn from GAMA, as well as from other KiDS-overlapping wide-angle spectroscopic surveys, namely SDSS, 2dFLenS, and 2dFGRS. We find stable behavior and consistent improvement over ANNz2 throughout. We finally apply Hybrid-z trained on GAMA to the entire KiDS-Bright DR4 sample of 1.2 million galaxies. For these final predictions, we design a method of smoothing the input redshift distribution of the training set, to avoid propagation of features present in GAMA, related to its small sky area and large-scale structure imprint in its fields. Our work paves the way towards the best-possible photo-$z$s achievable with machine learning for any galaxy type both for the final KiDS-Bright DR5 data and for future deeper imaging, such as from the Legacy Survey of Space and Time.

The dwarf galaxies comparable to the LMC and SMC, with stellar masses $7.5 <{\rm log}(M_{\ast}/M_{\odot})<9.5$, are found in a diversity of environments and have long quenching timescales. The way this phenomenon results from the dwarfs' halo properties or their locations in the large-scale structure is not well understood. We study the star-formation rates of dwarfs in the Illustris TNG50 simulation across different environments, focusing on the dwarfs having host halos with virial masses $9 < {\rm log}(M_{200}/M_{\odot}) < 11.5$, which we demonstrate are in the field as opposed to dwarf satellites in hosts with ${\rm log}(M_{200}/M_{\odot}) \geq 11.5$. Our field dwarf sample is heterogeneous, consisting of primary (central) galaxies, with smaller numbers of secondaries and dwarf galaxies that are on backsplash orbits around massive galaxies. We quantify the field dwarfs' large-scale environmental using the tidal index $\Theta_1$ and the distance to the nearest massive galaxy ${\rm d_{massive}}$, and we study how the quenched fraction and star-formation histories are correlated with these metrics. Those with ${\rm d_{massive}}>1.5$ Mpc and $\Theta_1<0$ are well isolated, with only $\sim 1\%$ being quenched. The quenched field dwarfs are in majority the backsplash dwarfs located in the neighborhood of cluster-scale halos. We discover a two-halo galactic conformity signal that arises from the tendency of the quenched dwarfs, particularly the backsplash sample, to have a quenched massive galaxy as a neighbor. We attribute the low quenched fractions of the simulated LMC/SMC analogs in the field to the locations of their low-mass hosts in the sparse large-scale environment, which predominate over the relatively small number of backsplash and pre-processed dwarfs in denser environments.

We propose models of Astrophysical Black holes (ABHs) without event horizons (EHs), as a more viable explanation for the long-term quenching phenomenon in galaxies. At the same time, the short-term quenching is explained here in terms of an efficient feedback expected in the models of stellar-mass astrophysical black holes (StMABHs). We have calculated the radiative flux from the disk in a general spherically symmetric metric background and used it to contrast the distinctive features of the BHs and ABHs scenarios. The nature of the feedbacks arising from accretion onto a BH and an ABH in the `quasar' and `radio' modes are compared and some possible observational signatures of the StMABHs are pointed out.

Light and air pollution are the two main forms of pollution, containing the highest concentration in urban areas. I set out to investigate the effects of anthropogenic air and light pollution on the night sky, how this affects the astronomical data-collecting process, and how the general public perceives both forms of pollution. This paper utilizes primary and secondary research, referring to different sources: data collected through my telescope, astrophotography, other scientists' research, and survey responses from over 40 people across 2 areas: Southern California and Guadalajara, Mexico. Through this, I have found strong correlations between increased aerosol particles from air pollutants and increased night sky brightness, increased aerosol particles amplifying cloud formation, which in turn increases light reflection, and a correlation between different areas and what they believe about pollution, with many being misinformed on air and light pollution themselves.

Dark photons are a well-motivated candidate for dark matter, but their detection becomes challenging for ultralight masses with both experimental and astrophysical probes. In this work, we propose a new approach to explore this regime through the dark inverse Compton scattering of ultralight dark photons with cosmic ray electrons and positrons. We show this process generates a potentially observable background radiation that is most prominent at frequencies below MHz. We compute this effect using the latest cosmic ray models and radio absorption maps. Comparing it to observations of the Milky Way's radio spectrum from Explorer 43, Radio Astronomy Explorer 2, and the Parker Solar Probe, we place leading constraints on the kinetic mixing of dark photon dark matter for masses $\lesssim 2 \times 10^{-17} \ \rm eV$.

Gibbons and Schiller have raised the physically interesting conjecture that forces in general relativity are bounded from above by the mathematically compact relation ${\cal F}\leq c^4/4G$. In the present compact paper we explicitly prove, using the non-linearly coupled Einstein-matter field equations, that the force function ${\cal F}\equiv 4\pi r^2 p(r)$ in {\it stable} self-gravitating horizonless matter configurations is characterized by the upper bound ${\cal F}\leq c^4/G$ [here $p(r)$ is the radial pressure inside the self-gravitating matter configuration].

Domain wall (DW) networks may have formed in the early universe following the spontaneous breaking of a discrete symmetry. Notably, several particle physics models predict the existence of current-carrying DWs, which can capture and store particles as zero modes on it. In this study, we demonstrate that gravitational waves (GWs) generated by current-carrying DWs with fermionic zeromodes exhibit a novel feature: an additional peak in the GW spectrum resembling mountains, arising from metastable topological remnants, which we term ``spherons.'' This distinct signature could be detectable in upcoming GW observatories such as LISA and ET. The results suggest that DW networks in beyond Standard Model scenarios could emit GW signals that are significantly stronger and with greater detectability than previously expected.

General relativity predicts that massless waves should scatter from the Riemann curvature of their backgrounds. These scattered waves are sometimes called $\textit{tails}$ and have never been directly observed. Here we calculate the gravitational waves scattered in the backward direction (scattering angle $\vartheta=\pi$) from the weak-field curvature of an extended massive object, finding close agreement with previous results in the long-wavelength limit. These long-wavelength results are then applied to gravitational waves in the Laser Interferometer Space Antenna (LISA) sensitivity band scattering from the Sun and planets. We estimate that scattering from the Sun and the planets from Jupiter to Neptune could contribute a $10^{-3}$ amplitude modulation to LISA observations when these objects almost intersect the line of sight between LISA and the source, leading to forward scattering with $\vartheta\simeq0$. These conditions should be realized for the Sun during the lifetime of LISA if the detection rate of long-duration sources is not much smaller than a thousand per year.

Axion-like particles (ALPs) can account for the observed dark matter (DM) of the Universe and if their masses are at the eV scale, they can decay into infrared, optical and ultraviolet photons with a decay lifetime larger than the age of the Universe. We analyze multi-wavelength data obtained from the central region of Messier 87 (M87) galaxy by several telescopes, such as, Swift, Astrosat, Kanata, Spitzer and International Ultraviolet Explorer in the infrared to ultraviolet frequencies ($\sim 2\times10^{14} \, {\rm Hz} - 3\times10^{15}$ Hz), to constrain the narrow emission lines indicative of the eV scale ALP DM decay. We derive constraints on the ALP coupling to two photons ($g_{a\gamma\gamma}$) for ALP mass range $2 \, {\rm eV} \lesssim m_a \lesssim 20 \, {\rm eV}$, assuming ALPs form the DM in the M87 halo. We find that our bounds on ALP-two-photon coupling can become stronger than the existing ones by an order of magnitude in the ALP mass range $8 \, {\rm eV} \lesssim m_a \lesssim 20 \, {\rm eV}$.

I make a brief review about the QCD phases and the equation of state inferred from the neutron star data. Along the temperature axis at low baryon density, the QCD phase transition is a smooth crossover, and it is a natural extension of our imagination to postulate a similar crossover along the density axis at low temperature. Even without phase transitions, the inferred thermodynamic properties of neutron star matter turn out to be highly nontrivial already at twice of the nuclear saturation density. I also give some discussions about the substantiation of quark matter by means of the gravitational wave signals including the multi-messenger prospect.

Scalar particles traveling faster than a subluminal gravitational wave generate gravitons via gravitational Cherenkov radiation. In this paper, we investigate graviton production by the primordial plasma within the framework of modified gravity in the early Universe, generating a relic graviton background. We find that for the minimal model, where only the speed of gravitational waves is modified and a standard model plasma minimally couples to gravity, the relic graviton background can be enhanced by several orders of magnitude, but still agrees with the Big Bang Nucleosynthesis (BBN) bound in most cases. Moreover, we also find that for Horndeski theories, such as Galileon theory, the relic background produced by the thermalized scalar field can reach significant amplitudes, exceeding the BBN bound for a region of the parameter space. By requiring the relic graviton background to remain consistent with the BBN constraint, we derive limits on the gravitational wave speed at early times in these modified gravity theories.