Papers for Tuesday, Aug 27 2024

The temperature and polarization anisotropies of the cosmic microwave background (CMB) have been used to set constraints on decaying dark matter models down to keV masses. In this work, we extend these limits down to the sub-keV mass range. Using principal component analysis, we estimate the lower bound on the decay lifetime for a basis of different dark matter masses and Standard Model final states, from which the bound on an arbitrary model can be calculated. We validate our principal component analysis using Markov chain Monte Carlo methods and Planck 2018 data. We perform a separate analysis for models decaying into photons below the hydrogen ionization threshold. We demonstrate that for these models, the effect of energy deposition can be captured approximately by a single parameter, but the redshift dependence of the effect is very different from higher-energy injections; in particular, the perturbations to CMB anisotropies are more sensitive to energy deposited around the time of recombination.

It has been proposed that the accelerated expansion of the universe can be explained by the merging of our universe with baby universes, resulting in dark energy with a phantom-like equation of state. However, the evidence in favor of it did not include the full set of cosmological observables. Here we examine the implications of this model for both early and late universe cosmology using data from Planck collaboration, DESI 2024 and other experiments. We find that the pure baby universe model gives a poor fit to current data. Extending it to include a contribution from the cosmological constant, we find two allowed regions of parameter space: one close to $\Lambda$CDM, and another with $\Lambda <0$ plus the exotic dark energy component. The two regions can be significantly favored over $\Lambda$CDM, depending on the choice of supernova datasets, and they can ameliorate the Hubble tension to the level of $2\sigma$, depending on the supernova dataset. The model with $\Lambda<0$ features an equation of state $w(a)$ with a pole singularity at early times.

The origin of Hot Jupiters (HJs) is disputed between a variety of in situ and ex situ formation scenarios. One of the early proposed ex situ scenarios was the Eccentric Kozai-Lidov (EKL) mechanism combined with tidal circularization, which can produce HJs with the aid of a stellar or planetary companion. However, observations have revealed a lack of stellar companions to HJs, which challenges the importance of the binary star-driven EKL plus tides scenario. In this work, we explore so far unaccounted-for stellar evolution effects on HJ formation, in particular the effect of white dwarf (WD) formation. Gaia observations have revealed that WDs often undergo a kick during formation, which can alter a binary's orbital configuration or even unbind it. Based on this WD kick, in this letter we propose and explore two novel HJ formation pathways: 1) HJs that are presently orbiting single stars, but were initially formed in a binary that was later unbound by a WD kick; 2) Binaries that survive the WD kick can trigger enhanced EKL oscillations and lead to 2nd generation HJ formation. We demonstrate that the majority of seemingly single HJs could have formed in binary star systems. As such, HJ formation in binaries via the EKL mechanism could be one of the dominant HJ formation pathways, and our results highlight that unaccounted-for stellar evolution effects, like WD formation, can obscure the actual origin of observed exoplanet populations.

Although exoplanetary science was not initially projected to be a substantial part of the Spitzer mission, its exoplanet observations set the stage for current and future surveys with JWST and Ariel. We present a comprehensive reduction and analysis of Spitzer's 4.5 micron phase curves of 29 hot Jupiters on low-eccentricity orbits. The analysis, performed with the Spitzer Phase Curve Analysis (SPCA) pipeline, confirms that BLISS mapping is the best detrending scheme for most, but not all, observations. Visual inspection remains necessary to ensure consistency across detrending methods due to the diversity of phase curve data and systematics. Regardless of the model selection scheme - whether using the lowest-BIC or a uniform detrending approach - we observe the same trends, or lack thereof. We explore phase curve trends as a function of irradiation temperature, orbital period, planetary radius, mass, and stellar effective temperature. We discuss the trends that are robustly detected and provide potential explanations for those that are not observed. While it is almost tautological that planets receiving greater instellation are hotter, we are still far from confirming dynamical theories of heat transport in hot Jupiter atmospheres due to the sample's diversity. Even among planets with similar temperatures, other factors like rotation and metallicity vary significantly. Larger, curated sample sizes and higher-fidelity phase curve measurements from JWST and Ariel are needed to firmly establish the parameters governing day-night heat transport on synchronously rotating planets.

RV Tauri variable stars are pulsating evolved stars that are identified by a characterizing feature in their light curves: alternating deep and shallow minima. Many RV Tauri variable stars were originally classified decades ago using visual observations and photographic plates. However, recent studies suggest that there are imposters, or misclassified variables, amongst the sample of RV Tau stars. In this study, we examine 84 known and suspected RV Tau stars that appear in the variable star catalog of the All Sky Automated Survey for Supernova (ASAS-SN). For each ASAS-SN light curve, we performed a period analysis and compared our results with those of the General Catalog of Variables Stars (GCVS) and the automatic classification algorithm used by ASAS-SN. We found that the pattern of alternating minima present in RV Tau variable stars often confuses automatic classification. Our results include updated periods and classifications for our sample. Our study provides an important step towards obtaining a robust sample of RV Tau variable stars to better understand the pulsation mechanism and the evolutionary pathway of these variable stars.

The X-ray transient Swift J1727.8-1613 ended its 10-month discovery outburst on June of 2024, when it reached an optical brightness comparable to pre-discovery magnitudes. With the aim of performing a dynamical study, we launched an optical spectroscopy campaign with the GTC telescope. We detect the companion star and construct its radial velocity curve, yielding a binary orbital period of Porb = 10.8041 +- 0.0010 h and a radial velocity semi-amplitude of K2 = 390 +- 4 km/s. This results in a mass function of f(M1)=2.77 +- 0.09 Msun. Combined with constraints on the binary inclination, it sets a lower limit to the compact object mass of M1 > 3.12 +- 0.10 Msun, dynamically confirming the black hole nature of the accretor. Comparison of the average spectrum in the rest frame of the companion with synthetic stellar templates supports a K4V donor partially veiled (74%) by the accretion disc. A refined distance measurement of 3.7+- 0.3 kpc, together with the astrometric proper motion and the systemic velocity derived from the radial velocity curve (181 +-4 km/s), supports a natal kick velocity of 220 +30 -40 km/s, at the upper end of the observed distribution.

Mapping stars and gas in nearby galaxies is fundamental for understanding their growth and the impact of their environment. This issue is addressed by comparing the stellar "edges" of galaxies $D_{\rm stellar}$, defined as the outermost diameter where in situ star formation significantly drops, with the gaseous distribution parameterized by the neutral atomic hydrogen diameter measured at 1 $M_{sun}$/pc$^2$, $D_{HI}$. By sampling a broad HI mass range $10^5 M_{sun} < M_{HI} < 10^{11} M_{sun}$, we find several dwarf galaxies with $M_{HI} < 10^9 M_{sun}$ from the field and Fornax Cluster which are distinguished by $D_{\rm stellar} >> D_{HI}$. For the cluster dwarfs, the average HI surface density near $D_{\rm stellar}$ is $\sim$0.3 $M_{sun}$/pc$^2$, reflecting the impact of quenching and outside-in gas removal from ram pressure and tidal interactions. In comparison, $D_{\rm stellar}/D_{HI}$ ranges between 0.5-2 in dwarf field galaxies, consistent with the expectations from stellar feedback. Only more massive disk galaxies in the field can thus be characterized by the common assumption that $D_{\rm stellar} \lesssim D_{HI}$. We discover a break in the $D_{\rm stellar}-M_{\rm stellar}$ relation at $m_{break} \sim 4\times10^8 M_{sun}$ that potentially differentiates the low mass regime where the influence of stellar feedback and environmental processes more prominently regulates the sizes of nearby galaxies. Our results highlight the importance of combining deep optical and HI imaging for understanding galaxy evolution.

JWST has enabled the detection of the UV continuum of galaxies at z>10, evidencing a population of extremely blue, potentially dust-free galaxies. Interpreting the UV spectra of galaxies as they redden is complicated by the well-known degeneracy between stellar ages, dust, and nebular continuum. The main goal of this paper is to develop a theoretical model for the relationship between galaxy UV slopes, bursty star formation histories, dust evolution, and the contribution from nebular regions. We accomplish this via cosmological zoom-in simulations, and in specific, build a layered model where we simulate the UV slopes of galaxies with increasingly complex physics. Our main results follow. (i) Unattenuated stellar populations with no nebular emission exhibit a diverse range of intrinsic UV slopes, with values ranging from beta ~ -3 --> -2.2 due to long delays between bursts. This is manifested by an inverse correlation between the intrinsic UV slope and sSFR for early galaxies such that higher sSFR corresponds to bluer UV slopes. (ii) When including dust, our model galaxies demonstrate a rapid rise in dust obscuration between z ~ 8-10. This increase in dust mass is due to high grain-grain shattering rates, and enhanced growth per unit dust mass in very small grains, resulting in UV-detected galaxies at z ~ 12 descending into ALMA-detectable galaxies by z ~ 6. The rapid rise in dust content at z ~ 8-10 leads to a systematic reddening of the UV slopes during this redshift range. (iii) The inclusion of nebular continuum reddens the UV slope by a median factor Delta beta ~ 0.2-0.4. However, when including nebular continuum, our highest redshift galaxies (z~12) are insufficiently blue compared to observations; this may imply an evolving escape fraction from HII regions with redshift.

Enhanced emission in the months to years preceding explosion has been detected for several core-collapse supernovae (SNe). Though the physical mechanisms driving the emission remain hotly debated, the light curves of detected events show long-lived ($\geq$50 days), plateau-like behavior, suggesting hydrogen recombination may significantly contribute to the total energy budget. The Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) will provide a decade-long photometric baseline to search for this emission, both in binned pre-explosion observations after an SN is detected and in single-visit observations prior to the SN explosion. In anticipation of these searches, we simulate a range of eruptive precursor models to core-collapse SNe and forecast the discovery rates of these phenomena in LSST data. We find a detection rate of ~40-130 yr$^{-1}$ for SN IIP/IIL precursors and ~110 yr$^{-1}$ for SN IIn precursors in single-epoch photometry. Considering the first three years of observations with the effects of rolling and observing triplets included, this number grows to a total of 150-400 in binned photometry, with the highest number recovered when binning in 100-day bins for 2020tlf-like precursors and in 20-day bins for other recombination-driven models from the literature. We quantify the impact of using templates contaminated by residual light (from either long-lived or separate precursor emission) on these detection rates, and explore strategies for estimating baseline flux to mitigate these issues. Spectroscopic follow-up of the eruptions preceding core-collapse SNe and detected with LSST will offer important clues to the underlying drivers of terminal-stage mass loss in massive stars.

We present the results of Sun-as-a-star observations by the NEID Solar Telescope at WIYN Observatory, spanning January 1, 2021 through June 30, 2024. We identify 117,060 observations which are unlikely to be significantly affected by weather, hardware or major calibration issues. We describe several high-level data products being made available to the community to aid in the interpretation and inter comparisons of NEID solar observations. Solar observations demonstrate excellent performance of NEID, including radial velocity (RV) accuracy and long-term stability of better than $\simeq 0.37$ m s$^{-1}$ over $\simeq 3.5$ years, even though NEID was not originally designed or optimized for daytime observations of the Sun. Currently, intrinsic stellar variability is the primary barrier to detecting Earth-analog planets for most nearby, Sun-like stars. We present a comparison of the effectiveness of several methods proposed to mitigate the effects of solar variability on the Sun's estimated RV. We find that the Scalpels algorithm performs particularly well and substantially reduces the RMS RV of solar spectra from over 2 m s$^{-1}$ to 0.277 m s$^{-1}$. Even when training on a subset of days with NEID solar observations and testing on a held-out sample, the RMS of cleaned RV is 0.34-0.42 m s$^{-1}$. This is significantly better than previous attempts at removing solar variability and suggests that the current generation of EPRV instruments are technically capable of detecting Earth-mass planets orbiting a solar twin if provided with sufficient observing time allocations ($\sim 10^3$ nights of observations).

We present thermal emission measurements of GJ 1132b spanning 5--12 um obtained with the Mid-Infrared Instrument Low-Resolution Spectrometer (MIRI/LRS) on the James Webb Space Telescope (JWST). GJ 1132b is an M-dwarf rocky planet with Teq=584 K and an orbital period of 1.6 days. We measure a white-light secondary eclipse depth of 140+/-17 ppm, which corresponds to a dayside brightness temperature of Tp,dayside= 709+/-31 K using improved star and planet parameters. This measured temperature is only 1 sigma below the maximum possible dayside temperature of a bare rock (i.e., assuming a zero albedo planet with no heat redistribution, Tmax = 746+14/-11 K). The emission spectrum is consistent with a featureless blackbody, which agrees with a wide range of possible surface compositions. By comparing forward models to the dayside emission spectrum, we rule out Earth-thickness (P ~ 1 bar) atmospheres with at least 1% H2O, atmospheres of any modeled thickness (10^-4 -- 10^2 bar) that contain at least 1% CO2, and thick, Venus-like atmospheres (P>~100 bar) with at least 1 ppm CO2 or H2O. We therefore conclude that GJ 1132b likely does not have a significant atmosphere. This finding supports the concept of a universal 'Cosmic Shoreline' given the high level of bolometric and XUV irradiation received by the planet.

The M31 nucleus contains a supermassive black hole embedded in a massive stellar disk of apsidally-aligned eccentric orbits. It has recently been shown that this disk is slowly precessing at a rate consistent with zero. Here we demonstrate using N-body methods that apsidally-aligned eccentric disks can form with a significant (~0.5) fraction of orbits counter-rotating as the result of a gravitational wave recoil kick of merging supermassive black holes. Higher amplitude kicks map to a larger retrograde fraction in the surrounding stellar population which in turns results in slow precession. We furthermore show that disks with significant counter-rotation are more stable (that is, apsidal-alignment is most pronounced and long lasting), more eccentric, and have the highest rates of stars entering the black hole's tidal disruption radius

The values of the Hubble constant obtained from the data of the early and of the late Universe differ by 10%. Data obtained by DESI Collaboration can help solve this problem, establishing with high precision the value of the product of the Hubble constant h and comoving of sound horizon at the end of drag epoch to be equal to 101.8 Mpc. This value agrees very well with the earlier refinement of the standard recombination theory and the value of the Hubble constant H0=73.5 km/s/Mpc obtained from measurements in the late Universe.

The Cosmic Distance Duality Relation (CDDR) defines a fundamental link in modern cosmology between geometric distance measurements and redshift. Any observed deviations from the CDDR could suggest the need for new physics beyond the $\Lambda$CDM model. In this study, we model geometric distances (both luminosity and angular) using a cosmographic expansion with the Pade method, which stabilizes the behavior of the cosmographic series at high redshifts. This approach allows us to determine geometric distances independently of any specific model. By incorporating updated data from supernovae (SN), baryon acoustic oscillations (BAO), and cosmic chronometers (CC), we obtain observational constraints for parametric models that quantify the CDDR, alongside all baseline cosmographic parameters. We find that potential CDDR violations introduce new statistical correlations in the cosmographic parameters ($H_0$, $q_0$, and $j_0$). However, within this framework, we do not observe significant deviations from the CDDR, and our results remain consistent with the predictions of the $\Lambda$CDM model. This work offers a novel and straightforward method for testing the CDDR by fixing background evolution through cosmographic techniques, thereby paving the way for new geometric observational tests of possible deviations from standard cosmology.

Galaxy chemical enrichment mechanisms have primarily been constrained by [$\alpha$/Fe] and [Fe/H] measurements of individual stars and integrated light from stellar populations. However such measurements are limited at higher redshifts (z>1). Recently, we proposed an analogous diagram of the oxygen-to-argon abundance ratio, log(O/Ar), vs Ar abundance, 12+log(Ar/H), as a new diagnostic window for emission nebulae. In this Letter, using robust line flux measurements including temperature sensitive auroral lines, we present direct determination of O and Ar abundances in nine SFGs from JWST/NIRSPEC spectra at z$\sim$1.3-7.7, and two more with Keck/MOSFIRE spectra at z$\sim$2.2. Utilising their positions on the log(O/Ar) vs 12+log(Ar/H) plane, we present the first inference of galaxy chemical enrichment mechanisms from an ensemble of galaxies. The SFGs at z$\sim$1.3-3.4 are consistent with the solar neighbourhood galactic chemical enrichment models of the Milky Way Galaxy that are driven by core-collapse and Type Ia supernovae. Such enrichment mechanisms thus occur at least out to z$\sim$3.4. However, the highest-redshift SFGs (z$\sim$3.6-7.7) have very low log(O/Ar) values, revealing a different enrichment process at z>3.6. Such low log(O/Ar) values may be caused by a rapid but intermittent star-formation and/or additional sources. The new diagnostic window for SFGs enables us to reveal the unique fingerprints of galaxy chemical enrichment out to cosmic dawn.

Nowhere in the solar system are impact morphologies observed in greater variety than on the icy Galilean satellites. This is likely a consequence of the structural and thermal state of the crust at the time of impact, and perhaps impact velocity. Palimpsest-type impact features show smooth enclosed central plains surrounded by undulating plains, within which are distributed concentric arcuate ridges, and no recognizable rim. Buto Facula on Ganymede is the best resolved of any palimpsest, having mostly been imaged at 190 m/pixel and optimum lighting, allowing insight into the circumstances that form this type of impact feature. Part of an impact crater on the eastern edge of Buto Facula has been buried by undulating plains material, suggesting that at the time of their emplacement the undulating plains behaved as a low-viscosity flow advancing across the landscape around the impact zone, encroaching on landforms that it encountered. We evaluated hypotheses for the formation of undulating plains using impact and ejecta modeling. We do not attribute the source of Buto's undulating plains to "dry" impact ejecta due to the existence of impact features larger than Buto that are not surrounded by undulating plains deposits. We performed iSALE impact simulations incorporating a subsurface liquid layer (or low strength layer) at various depths. An impact into a surface with a pre-existing, 5 km-deep, 5 km-thick fluid layer results in excavation of fluid material from that layer, producing a nearly flat final surface profile that is consistent with Buto's flat profile and the distribution of its undulating plains material. A liquid layer at a depth of 20-40 km results in impact feature profiles that resemble classic impact craters. We offer tests of this shallow subsurface liquid layer hypothesis for palimpsest formation that could potentially be performed by JUICE and Europa Clipper.

Type-I solar noise storms are perhaps the most commonly observed active radio emissions from the Sun at meter-wavelengths. Noise storms have a long-lived and wideband continuum background with superposed islands of much brighter narrowband and short-lived emissions, known as type-I bursts. There is a serious paucity of studies focusing on the morphology of these two types of emissions, primarily because of the belief that coronal scattering will always wash out any features at small angular scales. However, it is important to { investigate} their spatial structures in detail to make a spatio-temporal connection with observations at extreme-ultraviolet/ X-ray bands to understand the detailed nature of these emissions. In this work, we use high angular resolution observations from the upgraded Giant Metrewave Radio Telescope to demonstrate that it is possible to detect structures with angular scales as small as $\sim 9\arcsec$, about three times smaller than the smallest structure reported to date from noise storms. Our observations also suggest while the individual type-I bursts are narrowband in nature, the bursts are probably caused by traveling disturbance(s) inducing magnetic reconnections at different coronal heights, and thus leading to correlated change in the morphology of the type-I bursts observed at a wide range of frequencies.

Using gamma-ray bursts as standard candles for cosmological parameter constraints rely on their empirical luminosity relations and low-redshift calibration. In this paper, we examine the Amati relation and its potential corrections based on the A118 sample of higher-quality gamma-ray bursts, using both Hubble data set and Pantheon+ samples as calibration samples in the redshift range of z < 1.965. In calibrating gamma-ray bursts using these two datasets, we employ Gaussian processes to obtain corresponding Hubble diagrams to avoid the dependence on cosmological models in the calibration process. We first divided the low-redshift sample of GRBs into two bins and examined the Amati relation and its potential modifications. We found that under both calibrations, the Amati relation did not show evidence of redshift evolution (68% confidence level). For the other two Amati relations that include redshift evolution terms, the central values of the redshift evolution coefficients deviated from 0, but due to the limitations of the sample size and the increase in the number of parameters, most of the redshift evolution coefficients were not able to be excluded from 0 at the 1 sigma level. Therefore, to assess their situation across the entire redshift range, we employed MCMC to globally fit three types of Amati relations. By computing AIC and BIC, we found that for the GRB A118 sample, the standard Amati relation remains the most fitting empirical luminosity formula, and no potential redshift evolution trend was observed for two different low-redshift calibrating sources.

We introduce a simple parametric model of the radio-infrared correlation (i.e., the ratio between the IR luminosity and the 1.4 GHz radio luminosity, $q_{\mbox{\tiny IR}}$) by considering the energy loss rate of high-energy cosmic ray (CR) electron governed by the radiative cooling (synchrotron, bremsstrahlung, inverse Compton scattering), ionization, and adiabatic expansion. Each process of CR electron energy loss is explicitly computed and compared to each other. We rewrite the energy loss rate of each process to be dependent on the gas surface density and redshift using the relevant scaling relations. By combining each energy loss rate, the fraction of the synchrotron energy loss rate is computed as a function of gas surface density and redshift, and used to extrapolate the well-established `local' radio-infrared correlation to the high-redshift universe. The locally established $q_{\mbox{\tiny IR}}$ is reformulated to be dependent upon the redshift and the gas surface density and applied for understanding the observed distribution of the radio-infrared correlation of high-redshift galaxies in \cite{delvecchio_etal_2021}. Our model predicts that $q_{\mbox{\tiny IR}}$ value is anti-correlated with gas surface density and the redshift dependency of $q_{\mbox{\tiny IR}}$ value changes by gas surface density of galaxies, which captures the observed trend of $q_{\mbox{\tiny IR}}$ values for stellar mass selected star forming galaxies with a minimal impact of radio-infrared selection bias.

In this work, we study the X-ray spectral and temporal properties of an ultraluminous X-ray source (ULX) in NGC 628 by using multi-epoch archival X-ray data. The physical parameters were estimated in each epoch in order to constrain the nature of the compact object in the system. Also, the optical counterpart candidates of the ULX were examined using the archival {\it Hubble} Space Telescope (HST)/Wide Field Camera 3 (WFC3) data. {\it XMM-Newton}, {\it Chandra}, and {\it Swift} data were used to create the long-term light curve (which covers a period of 22 years) and perform the spectral analysis. Lomb-Scargle periodograms of the source were constructed to examine the short-term variability in each epoch. In order to search for an optical counterpart in the HST/WFC3 images, a relative astrometric correction was initially applied to the {\it Chandra} and HST/WFC3 images. The X-ray flux of the source changes by a factor of $\sim 200$ throughout the observations. The previously detected quasi-periodic signal (in the range of 0.1$-$0.4 mHz) was confirmed by using the Lomb-Scargle method. After astrometric correction, two optical counterpart candidates were detected for the source. The obtained spectral energy distributions in the optical band for both candidates indicate that the optical emission is dominated by the irradiation of the accretion disc. Considering the best-fit model parameters of the multi-colour disc black-body model, we derived the mass of the black hole in the system as being in the range of (3$-$16) $M_{\odot}$. Nonetheless, the long-term variability and the spectral transitions in the hardness-luminosity diagram make it difficult to rule out the neutron star scenario.

In this review we show that the space experiment with gamma-ray detector with sensitivity 2 orders of magnitude higher than existing ones will make it possible to discover up to a thousand neutron star mergers, even at those moments when gravitational wave (GW) antennas are not working. At the same time, synchronous detection of neutron stars mergers by gamma-ray and GW detectors will make it possible not only to study in detail the physical processes occurring at the time of the catastrophe, but also to determine the full gamma ray beam pattern, including the average jet divergence angle and the real energy of the explosion. A gamma detector that has the required sensitivity at a relatively low flight weight is proposed. The latter, in turn, will make it possible to clarify our ideas about the genesis of double relativistic stars in the Universe.

HD 209458 b is one of the most studied exoplanets to date. Despite this, atmospheric characterisation studies yielded inconsistent species detections and abundances. Values reported for the C/O ratio range from 0.1 to 1.0. Of particular interest is the simultaneous detection of H2O and HCN reported by some studies using high-resolution ground-based observations, which would require the atmospheric C/O ratio to be fine-tuned to a narrow interval around 1. HCN has however not been detected from recent space-based observations. We aim to provide an independent study of HD 209458 b's atmosphere with high-resolution observations, in order to infer the presence of several species, including H2O and HCN. We observed four primary transits of HD 209458 b at a high resolution (R=92000) with CRIRES+ in the near infrared (band H, 1.4--1.8 um). After reducing the data with pycrires, we prepared the data using the SysRem algorithm and performed a cross-correlation (CCF) analysis of the transmission spectra. We also compared the results with those obtained from simulated datasets constructed by combining the Exo-REM self-consistent model with the petitRADTRANS package. Combining the four transits, we detect H2O with a signal-to-noise CCF metric of 8.7. This corresponds to a signal emitted at $K_p=151.3^{+31.1}_{-23.4}$ km/s and blueshifted by $-6^{+1}_{-2}$ km/s, consistent with what is expected for HD 209458 b. We do not detect any other species among C2H2, CH4, CO, CO2, H2S, HCN, and NH3. Comparing this with our simulated datasets, this result is consistent with a C/O ratio of 0.1 and an opaque cloud top pressure of 50 Pa, at a 3 times solar metallicity. This would also be consistent with recent JWST observations. However, none of the simulated results obtained with a bulk C/O ratio of 0.8, a value suggested by previous studies using GIANO-B and CRIRES, are consistent with our observations.

Compact star-forming galaxies were dominant galaxy types in the early Universe. Blueberry galaxies (BBs) represent their local analogues being very compact and having intensive star formation. Motivated by high X-ray emission recently found in other analogical dwarf galaxies, called Green Peas, we probe into the X-ray properties of BBs to determine if their X-ray emission is consistent with the empirical laws for star-forming galaxies. We performed the first X-ray observations of a small sample of BBs with the XMM-Newton satellite. Spectral analysis for detected sources and upper limits measured via Bayesian-based analysis for very low-count measurements were used to determine the X-ray properties of our galaxy sample. Clear detection was obtained only for 2 sources, with one source exhibiting an enhanced X-ray luminosity to the scaling relations. For the remaining 5 sources, only an upper limit was constrained, suggesting BBs to be rather underluminous as a whole. Our analysis shows that the large scatter cannot be easily explained by the stochasticity effects. While the bright source is above and inconsistent at almost the 99% confidence level, the upper limits of the two sources are below the expected distribution. These results indicate that the empirical relations between the star formation rate, metallicity, and X-ray luminosity might not hold for BBs with uniquely high specific star formation rates. One possible explanation could be that the BBs may not be old enough to have a significant X-ray binary population. The high luminosity of the only bright source can be then caused by an additional X-ray source, such as a hidden active galactic nucleus or more extreme ultraluminous X-ray sources.

The Rapid and accurate identification of Gamma-Ray Bursts (GRBs) is crucial for unraveling their origins. However, current burst search algorithms frequently miss low-threshold signals or lack universality for observations. In this study, we propose a novel approach utilizing transfer learning experiment based on convolutional neural network (CNN) to establish a universal GRB identification method, which validated successfully using GECAM-B data. By employing data augmentation techniques, we enhance the diversity and quantity of the GRB sample. We develop a 1D CNN model with a multi-scale feature cross fusion module (MSCFM) to extract features from samples and perform classification. The comparative results demonstrated significant performance improvements following pre-training and transferring on a large-scale dataset. Our optimal model achieved an impressive accuracy of 96.41% on the source dataset of GECAM-B, and identified three previously undiscovered GRBs by contrast with manual analysis of GECAM-B observations. These innovative transfer learning and data augmentation methods presented in this work hold promise for applications in multi-satellite exploration scenarios characterized by limited data sets and a scarcity of labeled samples in high-energy astronomy.

To test the hypothesis that megaparsec-scale giant radio galaxies (GRGs) have multiple epochs of recurrent activity which lead to their giant sizes and to understand the nature of double-double radio galaxies (DDRGs), we have built the largest sample of giant DDRGs from the LOFAR Two Metre Sky Survey (LoTSS). The sample consists of 111 sources, of which 76 are newly identified. Their redshift and total projected size ranges are between 0.06 to 1.6 and 0.7 Mpc to 3.2 Mpc, respectively. We conducted a detailed analysis to characterise their properties, including arm-length ratios, flux density ratios of pairs of lobes, and misalignment angles. These measurements allow us to study the symmetry parameters, which are influenced by the immediate and large-scale environments of DDRGs. Our study conclusively shows that the cocoons in which the inner lobes of DDRGs grow are often contaminated with surrounding material from the external medium, leading to observed asymmetries in the inner lobes. Our analysis also reveals highly misaligned DDRGs, which could be due to environmental factors and/or changes in the supermassive black hole jet ejection axes. By studying the misalignment angles, we assess the stability of the jets in these systems in relation to their environment. Additionally, we have also characterised the large-scale environments of DDRGs, identifying them within dense galaxy clusters and studying how cluster weather significantly affects their morphologies. Interestingly, we have identified two gigahertz peaked-spectrum (GPS) candidates in the unresolved cores of the DDRGs, as well as one triple-double candidate, which, if confirmed, would be only the fifth known case. Additionally, we have discovered a DDRG in a distant galaxy cluster at $z\sim$\,1.4. Abridged.

The cosmic microwave background (CMB) lensing power spectrum is a powerful probe of the late-time universe, encoding valuable information about cosmological parameters such as the sum of neutrino masses and dark energy equation of state. However, the presence of anisotropic cosmic birefringence can bias the reconstructed CMB lensing power spectrum using CMB polarization maps, particularly at small scales, and affect the constraints on these parameters. Upcoming experiments, which will be dominated by the polarization lensing signal, are especially susceptible to this bias. We identify the dominant contribution to this bias as an $N_L^{(1)}$-like noise, caused by anisotropic rotation instead of lensing. We show that, for an CMB-S4-like experiment, a scale-invariant anisotropic rotation field with a standard deviation of 0.05 degrees can suppress the small-scale lensing power spectrum ($L\gtrsim 2000$) at a comparable level to the effect of massive neutrino with $\sum_i m_{\nu_{i}}=50~\rm{meV}$, making rotation field an important source of degeneracy in neutrino mass measurement for future CMB experiments. We provide an analytic expression and a simulation-based estimator for this $N_L^{(1)}$-like noise, which allows for efficient forecasting and mitigation of the bias in future experiments. Furthermore, we investigate the impact of a non-scale-invariant rotation power spectrum on the reconstructed lensing power spectrum and find that an excess of power in the small-scale rotation power spectrum leads to a larger bias. Our work provides an effective numeric framework to accurately model and account for the bias caused by anisotropic rotation in future CMB lensing measurements.

Binary star systems are expected to follow Newtonian dynamics similarly to planetary systems. However, reports have been made of wide binary systems with separations around 0.01 pc and larger, showing potential deviations from standard Newtonian motion. This phenomenon, suggestive of the flattening of galactic rotation curves, calls for closer inspection. This study presents an analysis of wide binary stars using data from Gaia Data Release 3 (DR3), a space-based astrometry mission funded by the European Space Agency. The study compares different choices of selection criteria to examine the nature of these apparent anomalous kinematics within the solar neighbourhood. The Gaia data set furnishes detailed parameters, including radial velocity, mass, age, and binary probability in addition to standard parameters. A custom Python tool named BYNARY facilitated both data processing and analysis. This report reveals that the signs of any anomalous signals systematically diminish as the initial selection criteria are relaxed for degrouping while subsequent filtering remains stringent, leading to the complete elimination of any apparent non-Newtonian motion for binary separations within 0.5 pc. The study shows that an y observ ed anomalous behaviour in solar neighbourhood wide binaries within 130 pc must be produced either by faint companion stars orbiting primary or secondary stars, or by flyby stars. The findings emphasize the importance of the choice of selection criteria in disentangling genuine binary dynamics from external influences. The conclusions align with the predictions of Newtonian mechanics and general relativity, though they do not exclude other phenomena at larger scales.

Accurate classification of celestial objects is essential for advancing our understanding of the universe. MargNet is a recently developed deep learning-based classifier applied to SDSS DR16 dataset to segregate stars, quasars, and compact galaxies using photometric data. MargNet utilizes a stacked architecture, combining a Convolutional Neural Network (CNN) for image modelling and an Artificial Neural Network (ANN) for modelling photometric parameters. In this study, we propose enhancing MargNet's performance by incorporating attention mechanisms and Vision Transformer (ViT)-based models for processing image data. The attention mechanism allows the model to focus on relevant features and capture intricate patterns within images, effectively distinguishing between different classes of celestial objects. Additionally, we leverage ViTs, a transformer-based deep learning architecture renowned for exceptional performance in image classification tasks. We enhance the model's understanding of complex astronomical images by utilizing ViT's ability to capture global dependencies and contextual information. Our approach uses a curated dataset comprising 240,000 compact and 150,000 faint objects. The models learn classification directly from the data, minimizing human intervention. Furthermore, we explore ViT as a hybrid architecture that uses photometric features and images together as input to predict astronomical objects. Our results demonstrate that the proposed attention mechanism augmented CNN in MargNet marginally outperforms the traditional MargNet and the proposed ViT-based MargNet models. Additionally, the ViT-based hybrid model emerges as the most lightweight and easy-to-train model with classification accuracy similar to that of the best-performing attention-enhanced MargNet.

We explain the large variety of star formation laws in terms of one single, simple law that can be inferred from the definition of the star formation rate and basic algebra. The resulting equation, $\SFR = \eff\ \Mcollapsing/\tauff$, although it has been presented elsewhere, is interpreted in terms of clouds undergoing collapse { rather than being turbulence-supported, an idea that different groups have pursued this century}. Under such assumption, one can explain the constancy of $\eff$, the different intra-cloud correlations observed in Milky Way's molecular clouds, as well as the resolved and unresolved extragalactic relationships between SFR and a measurement of the mass in CO, HCN, and CO+HI. We also explain why the slope of the correlation changes when the orbital time $\tauorb$ is considered instead of the free-fall time, and why estimations of the free-fall time from extragalactic observations skew the correlation, providing a false sublinear correlation. We furthermore show that the apparent nearly linear correlation between the star formation rate and the dynamical equilibrium pressure in the midplane of the galaxies, $\PDE$, is just a consequence of $\PDE$ values being dominated by the variation of the column density of molecular gas. All in all, we argue that the star formation law is driven by the collapse of cold, dense gas, which happens to be primarily molecular in the present Universe, and that the role of stellar feedback is just to shut down the star formation process, not to shape the star formation law.

As an inflated Hot Jupiter orbiting an early-type primary star in the evolved binary HD 202772 system, HD 202772 A b's presence invites a study of how such a planet forms and evolves. As a prelude to potential atmospheric characterization with the latest generation of observatories, we present a reduction and analysis of eclipse light curve observations of HD 202772 A b acquired with the Spitzer Space Telescope using the 3.6 and 4.5 $\mu$m channels. We find eclipse depths of $680\pm68$ and $1081^{+54}_{-53}$ ppm, respectively, corresponding to day-side effective temperatures of $2130^{+102}_{-91}$ and $2611^{+46}_{-49}$ K. The corresponding Bond albedos are consistent with the distribution of albedos for Hot Jupiters observed with both Spitzer and TESS. The heat redistribution efficiencies consistent with the Bond albedo range predicted by 1-D atmospheric models in radiative-convective equilibrium are $0.71\pm0.10$ and $0.03^{+0.03}_{-0.02}$, respectively, indicating a weak day-night contrast for the former and a strong contrast for the latter. Given this, and the unique environment in which this planet resides, we recommend follow-up observations with JWST to more precisely constrain its atmospheric composition and structure, as well as its host stellar environment, to elucidate if and how the atmospheres of these close-in giants evolve with host stars in binaries past the main sequence.

A crucial aspect of galaxy evolution is the pace at which galaxies build up their mass. We can investigate this hierarchical assembly by uncovering and timing accretion events that our galaxy has experienced. In the Milky Way, accreted debris has been previously identified in the local halo, thanks to the advent of Gaia data. We aim to couple this dataset with advancements in colour-magnitude diagram fitting techniques to characterise the building blocks of the Galaxy. Here we focus on the retrograde halo, specifically Thamnos and Sequoia. We do this as part of the ChronoGal project by fitting absolute colour-magnitude diagrams (using CMDft.Gaia) of samples of stars associated with these substructures, extracted from a local 5D Gaia DR3 dataset. Comparing their derived age and metallicity distributions with those of the expected contamination, from the dominant Gaia Enceladus and low energy in-situ populations, we can unveil the stellar population signatures of the progenitors of Sequoia and Thamnos. We show that both Thamnos and Sequoia have a metal-poor population ([Fe/H]<-1.5 dex) that is distinct from the expected contamination. The age distributions allow us to see that Sequoia formed half of its stars by a lookback time of 12 Gyr, while Thamnos is on average older, having formed half its stars at 12.3 Gyr. Gaia Enceladus and the low energy in-situ populations formed half of their stars by 12.1 Gyr and 12.9 Gyr respectively. This suggests that Thamnos was accreted earlier than Gaia Enceladus and Sequoia is the most recent accretion event. We have presented, for the first time, age distributions of the retrograde halo substructures: Sequoia and Thamnos. These are derived purely photometrically using CMD fitting techniques, which also provide metallicity distributions that successfully reproduce the spectroscopic distributions, highlighting the capability of CMDft.Gaia.

We present a detailed study of the MASTER OT J030227.28+191754.5 outburst in 2021-2022, reaching an amplitude of 10.2 mag and a duration of 60 d. The detections of (1) the double-peaked optical emission lines, and (2) the early and ordinary superhumps, established that MASTER OT J030227.28+191754.5 is an extremely energetic WZ Sge-type dwarf nova (DN). Based on the superhump observations, we obtained its orbital period and mass ratio as 0.05986(1) d and 0.063(1), respectively. These are within a typical range of low-mass-ratio DNe. According to the binary parameters derived based on the thermal-tidal instability model, our analyses showed that (1) the standard disk model requires an accretion rate $\simeq$ 10$^{20}$ g s$^{-1}$ to explain its peak optical luminosity and (2) large mass was stored in the disk at the outburst onset. These cannot be explained solely by the impact of its massive ($\gtrsim$ 1.15 M$_\odot$) primary white dwarf implied by Kimura et al. (2023). Instead, we propose that the probable origin of this enormously energetic DN outburst is the even lower quiescence viscosity than other WZ Sge-type DNe. This discussion is qualitatively valid for most possible binary parameter spaces unless the inclination is low ($\lesssim 40^\circ$) enough for the disk to be bright explaining the outburst amplitude. Such low inclinations, however, would not allow detectable amplitude of early superhumps in the current thermal-tidal instability model. The optical spectra at outburst maximum showed the strong emission lines of Balmer, He I, and He II series whose core is narrower than $\sim 800$ km s$^{-1}$. Considering its binary parameters, a Keplerian disk cannot explain this narrow component, but the presumable origin is disk winds.

Radiative transfer (RT) modelling is a necessary tool in the interpretation of observations of the thermal emission of interstellar dust. It is also often part of multi-physics modelling. In this context, the efficiency of radiative transfer calculations is important, even for one-dimensional models. We investigate the use of the so-called immediate re-emission (IRE) method for fast calculation of one-dimensional spherical cloud models. We wish to determine whether weighting methods similar to those used in traditional Monte Carlo simulations can speed up the estimation of dust temperature. We present the program DIES, a parallel implementation of the IRE method, which makes it possible to do the calculations also on graphics processing units (GPUs). We tested the program with externally and internally heated cloud models, and examined the potential improvements from the use of different weighted sampling schemes. The execution times of the program compare favourably with previous programs, especially when run on GPUs. On the other hand, weighting schemes produce only limited improvements. In the case of an internal radiation source, the basic IRE method samples the re-emission well, while traditional Monte Carlo requires the use of spatial importance sampling. Some noise reduction could be achieved for externally heated models by weighting the initial photon directions. Only in optically very thin models does weighting - such as the proposed method of forced first interaction - result in noise reduction by a factor of several. The IRE method performs well for both internally and externally heated models, typically without the need for any additional weighting schemes. With run times of the order of one second for our test models, the DIES program is suitable even for larger parameter studies.

We report the evidence of a hidden black hole (BH) in a low-mass galaxy, MaNGA 9885-9102, and provide a new method to identify active BH in low mass galaxies. This galaxy is originally selected from the MaNGA survey with distinctive bipolar H$\alpha$ blobs at the minor axis. The bipolar feature can be associated with AGN activity, while the two blobs are classified as the H II regions on the BPT diagram, making the origins confusing. The Swift UV continuum shows that the two blobs do not have UV counterparts, suggesting that the source of ionization is out of the blobs. Consistent with this, the detailed photoionization models prefer to AGN rather than star-forming origin with a significance of 5.8$\sigma$. The estimated BH mass is $M_{\rm BH}\sim$7.2$\times 10^5 M_\odot$ from the $M_{\rm BH}-\sigma_*$ relationship. This work introduces a novel method for detecting the light echo of BHs, potentially extending to intermediate mass, in low metallicity environments where the traditional BPT diagram fails.

The Dark Energy Spectroscopic Instrument (DESI) Peculiar Velocity Survey aims to measure the peculiar velocities of early and late type galaxies within the DESI footprint using both the Fundamental Plane and Tully-Fisher relations. Direct measurements of peculiar velocities can significantly improve constraints on the growth rate of structure, reducing uncertainty by a factor of approximately 2.5 at redshift 0.1 compared to the DESI Bright Galaxy Survey's redshift space distortion measurements alone. We assess the quality of stellar velocity dispersion measurements from DESI spectroscopic data. These measurements, along with photometric data from the Legacy Survey, establish the Fundamental Plane relation and determine distances and peculiar velocities of early-type galaxies. During Survey Validation, we obtain spectra for 6698 unique early-type galaxies, up to a photometric redshift of 0.15. 64\% of observed galaxies (4267) have relative velocity dispersion errors below 10\%. This percentage increases to 75\% if we restrict our sample to galaxies with spectroscopic redshifts below 0.1. We use the measured central velocity dispersion, along with photometry from the DESI Legacy Imaging Surveys, to fit the Fundamental Plane parameters using a 3D Gaussian maximum likelihood algorithm that accounts for measurement uncertainties and selection cuts. In addition, we conduct zero-point calibration using the absolute distance measurements to the Coma cluster, leading to a value of the Hubble constant, $H_0 = 76.05 \pm 0.35$(statistical) $\pm 0.49$(systematic FP) $\pm 4.86$(statistical due to calibration) $\mathrm{km \ s^{-1} Mpc^{-1}}$. This $H_0$ value is within $2\sigma$ of Planck Cosmic Microwave Background results and within $1\sigma$, of other low redshift distance indicator-based measurements.

Much interest surrounds the nature of the compact remnant formed in core collapse supernovae (SNe). One means to constrain its nature is to search for signatures of power injection from the remnant in the SN observables years after explosion. In this work, we conduct a large grid of 1D nonlocal thermodynamic equilibrium radiative transfer calculations of He-star explosions under the influence of magnetar-power injection from post-explosion age of about one to ten years. Our results for SN observables vary with He-star mass, SN age, injected power, or ejecta clumping. At high mass, the ejecta coolants are primarily O and Ne, with [OI]6330A, [OII]7325A, and [OIII]5000A dominating in the optical, and with strong [NeII]12.81micron in the infrared -- this line may carry more than half the total SN luminosity. For lower He-star masses, a greater diversity of coolants appear, in particular Fe, S, Ar, or Ni from the Si- and Fe-rich regions. All models tend to rise in ionization in time, with twice-ionized species (i.e., OIII, NeIII, SIII, or FeIII) dominating at ~10yr, although this ionization is significantly reduced if clumping is introduced. Our treatment of magnetar power in the form of high-energy electrons or X-ray irradiation yields similar results -- no X-rays emerge from our ejecta even at ten years because of high-optical depth in the keV range. An uncertainty of our work concerns the power deposition profile, which is not known from first principles, although this profile could be constrained from observations. Our magnetar-powered model he8p00 with moderate clumping yields a good match to the optical and near-infrared observations of Type Ib SN2012au at both ~300d (power of 1-2x10^41erg/s) and 2269d (power of 10^40erg/s). Unless overly ionized, we find that all massive magnetar-powered ejecta should be infrared luminous at 5-10yr through strong [NeII]12.81micron line emission.

A giant radio halo (RH) is a diffuse synchrotron emission observed on the scale of megaparsecs (Mpc), typically found in the central region of merging galaxy clusters. Its large size and steep spectrum suggest that it originates from the re-energization of an aged population of cosmic ray electrons (CREs), while the secondary leptons produced in the $pp$ hadronic collision of cosmic ray protons (CRPs) may contribute to the emission. In this study, we investigate the reacceleration model including both primary and secondary CREs, assuming that the primary CRs originate from internal galaxies. In our new method, we follow the cosmological evolution of each cluster and calculate the energy spectra and one-dimensional spatial distributions of CRs. The primary CRE model with $\sim 3$ Gyr duration of reacceleration successfully reproduces the statistical properties of the RHs observed in the recent LOFAR survey, as well as the spectrum and profile of the Coma cluster. The gamma-ray and neutrino emissions produced by reaccelerated CRPs are consistent with the upper limits. However, if the CRP injection rate is high and the secondary CREs become significant, the model with the required $\sim 3$ Gyr reacceleration overproduces the number of RHs. The limit on the CRP injection rate, $L_{\rm p} \lesssim 10^{41}$ erg/s, is significantly lower than that expected from the early starburst activity or jets from active galactic nuclei. This discrepancy requires a revision of either the model of CR supply from galaxies or the turbulent reacceleration model.

We construct two new summary statistics, the scale-dependent peak height function (scale-PKHF) and the scale-dependent valley depth function (scale-VLYDF), and forecast their constraining power on PNG amplitudes $\{f_\mathrm{NL}^\mathrm{local}, f_\mathrm{NL}^\mathrm{equil},f_\mathrm{NL}^\mathrm{ortho}\}$ and standard cosmological parameters based on ten thousands of density fields drawn from \textsc{Quijote} and \textsc{Quijote-PNG} simulations at $z=0$. With the Fisher analysis, we find that the scale-PKHF and scale-VLYDF are capable of capturing a wealth of primordial information about the Universe. Specifically, the constraint on the scalar spectral index $n_s$ obtained from the scale-VLYDF (scale-PKHF) is 12.4 (8.6) times tighter than that from the power spectrum, and 3.9 (2.7) times tighter than that from the bispectrum. The combination of the two statistics yields constraints on $\{f_\mathrm{NL}^\mathrm{local}, f_\mathrm{NL}^\mathrm{equil}\}$ similar to those from the bispectrum and power spectrum combination, but provides a 1.4-fold improvement in the constraint on $f_\mathrm{NL}^\mathrm{ortho}$. After including the power spectrum, its constraining power well exceeds that of the bispectrum and power spectrum combination by factors of 1.1--2.9 for all parameters.

A technique has recently been developed for tracking short-term spectral variations in Galactic cosmic rays (GCRs) using data from a single neutron monitor (NM), by collecting histograms of the time delay between successive neutron counts and extracting the leader fraction $L$ as a proxy of the spectral index. Here we analyze $L$ from four Antarctic NMs during 2015 March to 2023 September. We have calibrated $L$ from the South Pole NM with respect to a daily spectral index determined from published data of GCR proton fluxes during 2015--2019 from the Alpha Magnetic Spectrometer (AMS-02) aboard the International Space Station. Our results demonstrate a robust correlation between the leader fraction and the spectral index fit over the rigidity range 2.97--16.6 GV for AMS-02 data, with uncertainty 0.018 in the daily spectral index as inferred from $L$. In addition to the 11-year solar activity cycle, a wavelet analysis confirms a 27-day periodicity in the GCR flux and spectral index corresponding to solar rotation, especially near sunspot minimum, while the flux occasionally exhibited a strong harmonic at 13.5 days, and that the magnetic field component along a nominal Parker spiral (i.e., the magnetic sector structure) is a strong determinant of such spectral and flux variations, with the solar wind speed exerting an additional, nearly rigidity-independent influence on flux variations. Our investigation affirms the capability of ground-based NM stations to accurately and continuously monitor cosmic ray spectral variations in the long-term future.

Polycyclic aromatic hydrocarbon (PAH) dust emission has been proposed as an effective extinction-independent star formation rate (SFR) indicator in the mid-infrared (MIR), but this may depend on conditions in the interstellar medium. The coverage of the AKARI/Infrared Camera (IRC) allows us to study the effects of metallicity, starburst intensity, and active galactic nuclei on PAH emission in galaxies with $f_{\nu}(L18W)\lesssim 19$ AB mag. Observations include follow-up, rest-frame optical spectra of 443 galaxies within the AKARI North Ecliptic Pole survey that have IRC detections from 7-24 $\mu$m. We use optical emission line diagnostics to infer SFR based on H$\alpha$ and [O II]$\lambda\lambda 3726,3729$ emission line luminosities. The PAH 6.2 $\mu$m and PAH 7.7 $\mu$m luminosities ($L(PAH\ 6.2\ \mu m)$ and $L(PAH\ 7.7\ \mu m)$, respectively) derived using multi-wavelength model fits are consistent with those derived from slitless spectroscopy within 0.2 dex. $L(PAH\ 6.2\ \mu m)$ and $L(PAH\ 7.7\ \mu m)$ correlate linearly with the 24 $\mu$m-dust corrected H$\alpha$ luminosity only for normal, star-forming ``main-sequence" galaxies. Assuming multi-linear correlations, we quantify the additional dependencies on metallicity and starburst intensity, which we use to correct our PAH SFR calibrations at $0<z<1.2$ for the first time. We derive the cosmic star formation rate density (SFRD) per comoving volume from $0.15 \lesssim z \lesssim 1$. The PAH SFRD is consistent with that of the far-infrared and reaches an order of magnitude higher than that of uncorrected UV observations at $z\sim1$. Starburst galaxies contribute $\gtrsim 0.7$ of the total SFRD at $z\sim1$ compared to main-sequence galaxies.

The presence of Dark matter (DM) within a neutron star (NS) can substantially influence the macroscopic properties. It is commonly assumed that the pressure inside an NS is isotropic, but in reality, pressure is locally anisotropic. This study explores the properties of anisotropic NS with a subfraction of DM (isotropic) trapped inside. Implementing a two-fluid formalism with three Equations of State (EOS): AP3 (a realistic nucleon-nucleon interaction model), BSk22 (modeling atomic nuclei and neutron-matter), and MPA1 (considering relativistic effects in nuclear interactions). The properties of NS, such as mass ($M$), radius ($R$), and dimensionless tidal deformability ($\Lambda$), for various DM-anisotropic configurations, have been rigorously tested against observational constraints. These constraints include data from the binary NS merger GW170817, NICER x-ray measurements, and pulsar mass-radius observations. We observe that with increasing DM subfraction, higher anisotropies could also satisfy the observational constraints. Furthermore, increasing the coupling ($g$) between DM and its mediator leads to the formation of a core-halo structure, with a DM halo surrounding the baryonic matter (BM). Specifically, for coupling values of $g = 10^{-4}$, $10^{-3.7}$, and $10^{-3.5}$, we observe that the maximum radius ($R_{max}$) decreases with increasing anisotropy, which contrasts with the behavior at $g = 10^{-5}$ and in scenarios with no DM. Our analysis indicates that binary pulsar systems could potentially constrain the extent of admixed anisotropic NS or, more optimistically, provide evidence for the existence of DM-admixed anisotropic NS.

According to the hierarchical clustering scenario, galaxies like the Milky Way form hierarchically, and many supporting evidences have been found in the Galactic halo. However, most stars in the Milky Way are disk stars. Disk stars have almost lost their spatial distribution and kinematic features at birth, retaining solely chemical signatures. Identifying such substructures using abundances of iron or light elements is difficult due to their high degeneracy. Heavy elements, especially neutron capture elements have limited sources so have lower degeneracy, but spectral line fitting of these elements is tough, requiring mid to high resolution spectra, which are currently limited in sample size. This work utilizes the overall effect of changing in abundances to introduce the information of heavy elements, especially neutron capture elements to weaken the degeneracy of [Fe/H]. The analysis suggests the presence of at least 12 clusters, proves that the Galactic Disk form hierarchically. We excluded giants to prevent the possibly impact of mixing effects and that of too large differences in surface gravity, so the conclusions were all drawn from the solar neighborhood.

Gravitational waves induced by primordial perturbations serve as crucial probes for studying the early universe, providing a significant window into potential new physics during cosmic evolution. Due to the potentially large amplitudes of primordial perturbations on small scales, the contributions of high-order cosmological perturbations are highly significant. We propose a vertex approach applicable to the study of induced gravitational waves for arbitrary higher orders. Using the vertex approach and tree diagrams, we can directly derive the explicit expressions of higher-order induced gravitational waves without involving the complex and lengthy calculations of higher-order cosmological perturbations. Correlations between different tree diagrams correspond to the loop diagrams of two-point correlation functions of induced gravitational waves. Our investigation reveals that one-particle reducible diagrams impact tensor-scalar induced gravitational waves while leaving scalar induced gravitational waves unaffected.

The structure of the jet in Cen A is likely better revealed in X-rays than in the radio band, which is usually used to investigate jet proper motions. In this paper, we analyze Chandra ACIS observations of Cen A from 2000 to 2022 and develop an algorithm for systematically fitting the proper motions of its X-ray jet knots. Most of the knots had an apparent proper motion below the detection limit. However, one knot at a transverse distance of $520~\mathrm{pc}$ had an apparent superluminal proper motion of $2.7\pm0.4~\mathrm{c}$. This constrains the inclination of the jet to be $i<41\pm6^{\circ}$, and the velocity of this knot to be $\beta>0.94\pm0.02$. This agrees well with the inclination measured in the inner jet by the EHT, but contradicts previous estimates based on jet and counterjet brightness. It also disagrees with the proper motion of the corresponding radio knot, of $0.8\pm0.1~\mathrm{c}$, which further indicates that the X-ray and radio bands trace distinct structures in the jet. There are four prominent X-ray jet knots closer to the nucleus, but only one of these is inconsistent with being stationary. A few jet knots also have a significant proper motion component in the non-radial direction. This component is typically larger closer to the center of the jet. We also detect brightness and morphology variations at a transverse distance of $100~\mathrm{pc}$ from the nucleus.

The shape of the cold interstellar molecular gas is determined by several processes, including self-gravity, tidal force, turbulence, magnetic field, and galactic shear. Based on the 3D dust extinction map derived by Vergely et al., we identify a sample of 550 molecular clouds (MCs) within $3\hspace{0.2em}\rm kpc$ of the solar vicinity in the Galactic disk. Our sample contains clouds whose size ranges from pc to kiloparsec, which enables us to study the effect of Galactic-scale processes, such as shear, on cloud evolution. We find that our sample clouds follow a power-law mass-size relation of $M\propto 32.00\hspace{0.2em}{R_{\rm{max}}}^{1.77}$, $M\propto 20.59\hspace{0.2em}{R_S}^{2.04}$ and $M\propto 14.41\hspace{0.2em}{R_V}^{2.29}$, where $R_{\rm{max}}$ is the major axis-based cloud radius, $R_S$ is the area-based radius, and $R_V$ is the volume-based radius, respectively. These clouds have a mean constant surface density of $\sim 7 \hspace{0.2em} \rm {M_{\odot}pc^{-2}}$, and follow a volume density-size relation of $\rho \propto 2.60\hspace{0.2em}{R_{\rm{max}}}^{-0.55}$. As cloud size increases, their shapes gradually transition from ellipsoidal to disk-like to bar-like structures. Large clouds tend to have a pitch angle of $28^{\circ} - 45^{\circ}$, where the angle is measured concerning the Galactic tangential direction. These giant clouds also tend to stay parallel to the Galactic disk plane and are confined within the Galactic molecular gas disk. Our results show that large molecular clouds in the Milky Way can be shaped by Galactic shear and confined in the vertical direction by gravity.

Magnetar outbursts are powered by an intense magnetic field. The phenomenon has recently drawn significant attention because of a connection to some fast radio bursts that has been reported. Understanding magnetar outbursts may provide the key to mysterious transient events. The elastic deformation of the solid crust due to magnetic field evolution accumulates over a secular timescale. Eventually, the crust fractures or responds plastically beyond a particular threshold. Determination of the critical limit is required to obtain the shear strain tensor in response to magnetic stress. In some studies, the tensor was substituted with an approximate expression determined algebraically from the magnetic stress. This study evaluated the validity of the approximation by comparing it with the strain tensor obtained through appropriate calculations. The differential equations for the elastic deformation driven by the magnetic field were solved. The results indicated that the approximation did not represent the correct strain tensor value, both in magnitude and spatial profile. Previous evolutionary calculations based on spurious criteria are likely to overestimate the magnitude of the strain tensor, and crustal failure occurs on a shorter timescale. Therefore, revisiting evolutionary calculations using the correct approach is necessary. This study is essential for developing the dynamics of crustal fractures and the magnetic-field evolution in a magnetar.

We report on the discovery of three large diffuse dwarf (LDD) galaxies located in isolated triple systems. They have effective diameters of (3.6-10.0) kpc and effective surface brightness of (26.2--27.3)$^m/$sq.arcsec. We note that the LDD galaxies tend to occur in small groups with a very low dispersion of radial velocities. The total (orbital) mass of the triplets approximately equals to their integral stellar mass within velocity measurement errors. The presence of LDD galaxies in cold multiple systems seems mysterious.

SRGe J194401.8+284452 is the brightest point-like X-ray object within the position uncertainty ellipse of an unidentified $\gamma$-ray source 4FGL J1943.9+2841. We performed multi-wavelength spectral and photometric studies to determine its nature and possible association with the $\gamma$-ray source. We firmly established its optical counterpart with the Gaia based distance of about 415 pc. Our data show that the object is a cataclysmic variable with an orbital period of about 1.5 hours, a late type donor star and an accretion disk around the white dwarf. SRGe J194401.8+284452 exhibits fast spontaneous transitions between the high and low luminosity states simultaneously in the optical and X-rays, remaining relatively stable between the transitions on scales of several months/years. This can be caused by an order of magnitude changes in the accretion rate. The brightness of the source is about 17 mag and 20 mag in the 2000 - 8000~A range and $5\times 10^{-12}$ and $5\times 10^{-13}$ erg/cm$^2$/s in the 0.3 -- 10 keV range in the high and low states, respectively. We constrained the mass of the white dwarf (0.3 -- 0.9 $M_\odot$) and its temperature in the low state (14750 $\pm$ 1250 K), the mass of the donor star ($\leq$ 0.08 $\pm$ 0.01 $M_ \odot$). In the low state, we detected regular optical pulsations with an amplitude of 0.2 mag and a period of 8 min. They are likely associated with the spin of the white dwarf, rather than with its non-radial pulsations. In the high state, the object demonstrates only stochastic optical brightness variations on time scales of 1 -- 15 minutes with amplitudes of 0.2 -- 0.6 mag. We conclude, that SRGe J194401.8+284452 based on its properties can be classified as an intermediate polar, and its association with the $\gamma$-ray source is very unlikely.

Superclusters are the largest-scale environments where a number of galaxy clusters interact with each other through minor/major mergers and grow via accretion along cosmic filaments. We focus on the A3528 complex in the core of the Shapley Supercluster. This system includes three clusters, A3528 (composed itself by two sub-clusters, namely A3528N and A3528S), A3532 and A3530, and presents a mildly active dynamical state. We study how minor mergers affect the evolution of radio galaxies and whether they are able to re-accelerate relativistic electrons in the ICM. We used observations from the uGMRT (Band 3, 4 and 5) and MeerKAT (L-band) telescopes to obtain images and spectral index maps over a wide frequency band and spatial resolutions. We compare these data with those from the SRG/eROSITA X-ray telescope. We detect faint diffuse radio emission associated with the radio galaxies. The BCGs in A3528S and A3532 show filaments of diffuse radio emission which extend for $\sim200-400$ kpc out of the radio galaxy. The spectral index of these filaments is extremely steep and almost constant ($\alpha\sim -2, -2.5$). Contrary to the radio tails in A3528N, the spectral properties of these radio filaments are not consistent with standard models of plasma ageing. We also detect roundish diffuse radio emission around the BCG in A3528S which could be classified as a radio mini-halo. The radio tail in this cluster appears longer that in earlier detections, being $\sim300$ kpc long at all frequencies. We linked the presence of extended radio emission in the form of filaments and threads in the A3528 complex with the effect of minor mergers. This is reinforced by the increasing X-ray fluctuations in correspondence with the radio extended emission in A3528S. Despite the less energy involved, our findings support the hypothesis that these events can re-energise plasma originating from radio galaxies.

Surface porosity has been found to be an important property for small bodies. Some asteroids and comets can exhibit an extremely high surface porosity in the first millimeter layer. This layer may be produced by various processes and maintained by the lack of an atmosphere. However, the influence of porosity on the spectro-photometric properties of small body surfaces is not yet fully understood. In this study, we looked into the effect of porosity on the spectro-photometric properties of Phobos regolith spectroscopic simulants; created by mixing the simulants with ultra-pure water, producing ice-dust particles, and then sublimating the water. The reflectance spectroscopic properties in the visible and near-infrared (0.5-4.2 $\mu$m) show no strong variations between the porous and compact samples. However, one simulant exhibits a bluing of the slope after increasing porosity, providing possible insights into the differences between the blue and red units observed on Phobos. In the mid-infrared range, a contrast increase of the 10-$\mu$m emissivity plateau due to silicates is observed. Photometry reveals a modification in the phase reddening behavior between the compact powder and the sublimation residue for both simulants. However, the observed behavior is different between the simulants, suggesting that the phase reddening may be dependent on the composition of the simulants. The phase curve also appears to be modified by the addition of porosity, with a higher contribution of forward scattering observed for the sublimation residue. The derivation of the Hapke parameters indicates an increase in roughness for the porous sample, but no significant modification of the opposition effect. This study aims to provide new insights into the understanding of porosity by using two Phobos simulants in the context of the upcoming JAXA/Martian Moons eXploration mission.

Classical novae are known to demonstrate a supersoft X-ray source (SSS) state following outbursts, which is associated with residual thermonuclear burning on the white dwarf (WD) surface. During its all-sky survey (eRASS1), the eROSITA telescope onboard the Spectrum-Roentgen-Gamma observatory discovered a bright new SSS, whose position is consistent with the known classical nova AT 2018bej in the Large Magellanic Cloud. There were two eROSITA spectra obtained during eRASS1 and eRASS2 monitoring epochs and one XMM-Newton grating spectrum close to the eRASS1 epoch. We aim to describe the eROSITA and XMM-Newton spectra of AT 2018bej with our local thermodynamic equilibrium (LTE) atmosphere models. We focused on the evolution of the hot WD properties between the eRASS1 and eRASS2 epochs, especially on the change of the carbon abundance. A grid of LTE model atmosphere spectra were calculated for different values of the effective temperature (from $T_{\rm eff}= 525$ to $700\,\rm kK$), surface gravity (six values) and chemical composition with five different values of carbon and nitrogen abundances. Both eRASS1 and XMM $0.3-0.6$ keV spectral analyses yield a temperature of the WD of $T_{\rm eff}{\sim}\,600\, \rm kK$ and a WD radius of $8000-8700\,\rm km$. Simultaneous fitting of the eROSITA spectra for two epochs (eRASS1 and eRASS2) with a common WD mass parameter demonstrates a decrease in $T_{\rm eff}$ accompanied by an increase in the WD radius and a decrease in the carbon abundance. However, these changes are marginal and coincide within errors. The derived WD mass is estimated to be $1.05-1.15\, M_\odot$. We traced a minor evolution of the source on a half-year timescale accompanied by a decrease in carbon abundance and concluded that LTE model atmospheres can be used to analyse the available X-ray spectra of classical novae during their SSS stage.

Solar eruptive events, including solar flares and coronal mass ejections (CMEs), are typically characterised by energetically significant X-ray emissions from flare-accelerated electrons and hot thermal plasmas. Occulted events, where the main flare is blocked by the solar limb, provide an opportunity to observe and analyse the X-ray emissions specifically associated with CMEs. This study investigates the X-ray and extreme ultraviolet (EUV) emissions associated with a large filament eruption and CME that occurred on February 15, 2022. This event was highly occulted from the three vantage points of Solar Orbiter, STEREO-A, and Earth. We utilised X-ray observations from the Spectrometer/Telescope for Imaging X-rays (STIX) and EUV observations from the Full Sun Imager (FSI) of the Extreme Ultraviolet Imager (EUI) on-board Solar Orbiter, supplemented by multi-viewpoint observations from STEREO-A/EUVI. This enabled a comprehensive analysis of the X-ray emissions in relation to the filament structure observed in EUV. We used STIX's imaging and spectroscopy capabilities to characterise the X-ray source associated with the eruption. Our analysis reveals that the X-ray emissions associated with the occulted eruption originated from an altitude exceeding 0.3Rsun above the main flare site. The X-ray time-profile showed a sharp increase and exponential decay, and consisted of both a hot thermal component at 17 MK and non-thermal emissions (>11.4 keV) characterised by an electron spectral index of 3.9. Imaging analysis showed an extended X-ray source that coincided with the EUV emission as observed from EUI, and was imaged until the source grew to a size larger than the imaging limit of STIX (180 arcsec). The findings demonstrate that STIX combined with EUI provides a unique and powerful tool for examining the energetic properties of the CME component of solar energetic eruptions.

We present a systematic study of the thermodynamic properties of the intracluster medium (ICM), including the plasma conditions in the ICM, in the merging galaxy cluster A3667 using Suzaku. We analyze the X-ray spectra of the ICM between the northwestern and the southeastern radio relics across a prominent cold front in A3667 to search for the ICM in a non-equilibrium ionization (NEI) state. We find that the ICM inside the cold front exhibits an NEI state with an ionization parameter of $(2.5\pm0.7)\times10^{11}\,\mathrm{s\,cm^{-3}}$ at the $90\%$ confidence level, which is lower than that expected from a collisional ionization equilibrium (CIE) state of the ICM (i.e., $n_\mathrm{e}t > 10^{12}\,\mathrm{s\,cm^{-3}}$). The timescale calculated from the ionization parameter of the ICM inside the cold front is $ 4.6\pm1.2\,\mathrm{Myr}$, which is much shorter than the thermal equilibration timescale. A weak transonic sloshing motion might explain the possibility that the ICM inside the cold front of A3667 is in the NEI state. In addition, the NEI state of the ICM in A3667 seems to be associated with the elongated radio halo bridging the region between the northwestern radio relic and the prominent cold front in A3667 across the cluster center.

The tidal Love numbers of self-gravitating compact objects describe their response to external tidal perturbations, such as those from a companion in a binary system, offering valuable insights into their internal structure. For static tidal fields, asymptotically flat black holes in vacuum exhibit vanishing Love numbers in general relativity, even though this property is sensitive to the presence of an external environment. In this work we study the tidal deformability of black holes surrounded by thin accretion disks, showing that the Love numbers could be large enough to mask any effect of modified gravity and to intrinsically limit tidal tests of black-hole mimickers. Furthermore, we investigate the measurability of the tidal parameters with next-generation gravitational wave experiments, like LISA and Einstein Telescope. Our findings suggest that these parameters could be measured with high precision, providing a powerful tool to probe the environment around coalescing binary systems.

We study the statistics and properties of microlensing events that can be detected by the Transiting Exoplanet Survey Satellite(TESS) based on Monte Carlo simulations. We simulate potential microlensing events from a sample of the TESS Candidate Target List(CTL) stars by assuming different observational time spans(or different numbers of sectors for each star) and a wide range of lens masses, i.e., $M_{\rm l}\in [0.1M_{\oplus},~2 M_{\odot}]$. On average, the microlensing optical depth and the event rate for CTL stars are $\simeq 0.2\times 10^{-9}$, and $\Gamma_{\rm{TESS}}\simeq0.6\times10^{-9}$ per star per day, respectively. The microlensing optical depth decreases by increasing the CTL priority, whereas the efficiency for detecting their microlensing signals enhances with the priority. Additionally, we simulate the microlensing events from the TESS Full-Frame Images(FFIs) stars extracted from the \texttt{TESS}-\texttt{SPOC} pipeline. The optical depth and event rate for these stars are on average $\simeq 1$-$3\times 10^{-9}$, and $\Gamma_{\rm{TESS}}\simeq 1$-$4\times 10^{-9}$ per star per day, and their highest values occur for sector $12$. The total number of microlensing events for the CTL stars is $N_{\rm e, \rm{tot}}\sim0.03$, whereas for the FFIs' stars number of events per star during $27.4$-day observing windows is $\hat{N}_{\rm e, \rm{tot}}\simeq1.4 \times 10^{-6}$. Based on four criteria we extract the detectable microlensing events and evaluate the detection efficiencies. The highest efficiency for detecting microlensing events from the TESS data occurs for the lens mass $\log_{10}[M_{\rm l}(M_{\odot})] \in [-4.5$,~$-2.5]$, i.e., super-Earth to Jupiter-mass Free Floating Planets(FFPs). The detectable microlensing events from the TESS stars are significantly affected by both finite-source and parallax effects.

The Pan-STARRS 3$\pi$ survey has detected hundreds of thousands of variable stars thanks to its coverage and 4-year time span, even though the sampling of the light curves is relatively sparse. These light curves contain only 10-15 detections in each of the five filters (g,r,i,z,y). During the K2 mission, the Kepler space telescope observed with a high sampling frequency, although only for about 80 days in each of its campaigns. Crossmatching and investigating the RR Lyrae stars observed by both K2 and Pan-STARRS can serve as a valuable tool to validate the classification and period determination of the survey. We used the Sesar catalogue of RR Lyrae stars detected by Pan-STARRS. After determining the overlap, we also considered the Gaia DR3 RR Lyrae catalogue data for these stars wherever it was available. The frequencies of the light variations were calculated by applying the Lomb-Scargle periodogram method on the K2 light curves that were prepared with autoEAP photometry. The calculated frequencies of the stars then were compared with those given in the Sesar catalogue and the Gaia DR3 RR Lyrae catalogue. We found that for the majority of the stars, the classification (95.6%) and the frequency determination (90.1%) of the PanSTARRS RR Lyrae stars were consistent within 0.03 d-1 with those that we derived from the K2 autoEAP light curves. For a significant subset of the sample, 7.4%, however, an offset of 1 or 2 d-1 was found in the frequencies. These are the result of the sampling of the detections, because Pan-STARRS observations are affected by diurnal cycles, whereas Kepler carried out measurements quasicontinuously. We found that RRc subtypes are significantly more affected (25.3%) than RRab subtypes (3.7%), which is most likely caused by RRc stars having less sharp light curve features. Validation via space-based data will be important for future ground-based surveys, as well.

The Atacama Large Millimeter/submillimeter Array (ALMA) offers new diagnostic capabilities for studying the Sun, providing complementary insights through high spatial and temporal resolution at millimeter wavelengths. ALMA acts as a linear thermometer for atmospheric gas, aiding in understanding the solar atmosphere's structure, dynamics, and energy balance. Given the Sun's complex emission patterns and rapid evolution, high-cadence imaging is essential for solar observations. Snapshot imaging is required, though it limits available visibility data, making full exploitation of ALMA's capabilities non-trivial. Challenges in processing solar ALMA data highlight the need for revising and enhancing the solar observing mode. The ALMA development study High-Cadence Imaging of the Sun demonstrated the potential benefits of high cadence observations through a forward modelling approach. The resulting report provides initial recommendations for improved post-processing solar ALMA data and explores increasing the observing cadence to sub-second intervals to improve image reliability.

We present Super-RDI, a unique framework for the application of reference star differential imaging (RDI) to Keck/NIRC2 high-contrast imaging observations with the vortex coronagraph. Super-RDI combines frame selection and signal-to-noise ratio (S/N) optimization techniques with a large multi-year reference point spread function (PSF) library to achieve optimal PSF subtraction at small angular separations. We compile a $\sim$7000 frame reference PSF library based on a set of 288 new Keck/NIRC2 $L'$ sequences of 237 unique targets acquired between 2015 and 2019 as part of two planet-search programs, one focusing on nearby young M dwarfs and the other targeting members of the Taurus star-forming region. For our dataset, synthetic companion injection-recovery tests reveal that frame selection with the mean-squared error (MSE) metric combined with KLIP-based PSF subtraction using 1000-3000 frames and $<$500 principal components yields the highest average S/N for injected synthetic companions. We uniformly reduce targets in the young M-star survey with both Super-RDI and angular differential imaging (ADI). For the typical parallactic angle rotation of our dataset ($\sim$10$^\circ$), Super-RDI performs better than a widely used implementation of ADI at separations $\lesssim$0.4" ($\approx$5 $\lambda$/$D$) gaining an average of 0.25 mag in contrast at 0.25" and 0.4 mag in contrast at 0.15". This represents a performance improvement in separation space over RDI with single-night reference star observations ($\sim$100 frame PSF libraries) applied to a similar Keck/NIRC2 dataset in previous work. We recover two known brown dwarf companions and provide detection limits for 155 targets in the young M-star survey. Our results demonstrate that increasing the PSF library size with careful selection of reference frames can improve the performance of RDI with the Keck/NIRC2 vortex coronagraph in $L'$.

{Sulfur Monoxide is known to be a good shock tracer in molecular clouds and protostar environments, but its abundance is difficult to reproduce even using state-of-the-art astrochemical models. } {We investigate the properties of the observed SO emission in the protoplanetary disk of AB Auriga, a Herbig Ae star of 2.4 $\Msun$ located at 156 pc. The AB Aur system is unique because it exhibits a dust trap and at least one young putative planet orbiting at about 30 au from the central star.} {We reduced ALMA archival data (projects 2019.1.00579.S, 2021.1.00690.S and 2021.1.01216.S) and analyzed the three detected SO lines (SO $6_5-5_4$, $6_7-5_6$ and $5_6-4_5$). We also present C$^{17}$O and C$^{18}$O 2-1 data to complement the interpretation of the SO data.} {For the three SO lines, the maximum SO emission in the ring is not located in the dust trap. Moreover, the inner radius of the SO ring is significantly larger than the CO emission inner radius, $\sim 160$ au versus $\sim 90$ au. The SO emission traces gas located in part beyond the dust ring. This emission likely originates from shocks at the interface of the outer spirals, observed in CO and scattered light emission, and the molecular and dust ring. SO is also detected inside the cavity, at a radius $\sim 20-30$ au and with a rotation velocity compatible with the proto-planet P1. We speculate that this SO emission originates from accretion shocks onto the circumplanetary disk of the putative proto-planet P1.} {These observations confirm that SO is a good tracer of shocks in protoplanetary disks and could be a powerful, new tool to detect embedded (proto-)planets.} Conclusions. These observations confirm that SO is a good tracer of shocks in protoplanetary disks and could be a powerful, new tool to detect embedded (proto-)planets.

The mass of galaxy clusters (GCs) can be determined by calculating the hydrostatic equilibrium equation. In this work, we revisit the derivation of hydrostatic mass of GCs within modified emergent Newtonian gravity (MENG) with generalized uncertainty principle (GUP) correction, Eddington-inspired Born-Infeld (EiBI) theory, and beyond Horndeski gravity (BHG). We apply the formulations on the masses of 10 GCs. We compare our results with the Newtonian mass of GCs and their baryonic masses $M_{bar}$. We find that all formulations could match with the Newtonian mass and the impact of the modified theories of gravity used in this work can be neglected. The noteworthy impact for MENG starts if we set $\beta_0=-1.656\times10^{110}$, $\kappa=5\times10^{40}$ for EiBI theory, and $\Upsilon=-0,1655\times10^{69}$ for BHG. In contrast, for the comparison with baryonic mass, only the MENG formulation produces the best linear fit with the slope $\frac{M_{MENG}}{M_{bar}}=0.755\pm0.061$. It means that in this work, MENG is the best theory of gravity that alleviates the mass discrepancy of hydrostatic mass and baryonic mass of GCs.

Direct imaging of exoplanets relies on complex wavefront sensing and control architectures. In addition to fast adaptive optics systems, most of the future high-contrast imaging instruments will soon be equipped with focal plane wavefront sensing algorithms. These techniques use the science detector to estimate the static and quasi-static aberrations induced by optical manufacturing defects and system thermal variations. Pair-wise probing (PWP) has been the most widely used, especially for space-based application and will be tested at contrast levels of ~1e-9 on-sky along with the future coronagraph instrument onboarding the Roman Space Telescope. This algorithm leans on phase diversities applied on the deformable mirror that are recorded in pairs. A minimum of two pairs of probes are required per bandwidth. An additional unprobed image is also recorded to verify the convergence rate of the correction. Before PWP, Borde & Traub proposed a similar algorithm that takes advantage of the unprobed image in the estimation process to get rid of the pair diversity requirement. In this work, we theoretically show that this latter technique should be more efficient than PWP when the convergence time is not limited by photon noise. We then present its performance and practical limitations on coronagraphic testbeds at JPL and exhibit a first on-sky control of non-common path aberrations with such method on VLT/SPHERE.

Detecting a compact subsolar object would have profound implications in physics, the reach of which depends on the nature of the object. Here we explore such consequences for a putative subsolar-mass gravitational wave event detected by the LIGO-Virgo-KAGRA Collaboration. We forecast that the nature of a subsolar binary (made of light neutron stars, primordial black holes, or more exotic compact objects) can be inferred with a great statistical confidence level already during the ongoing fourth observing run, based on the large tidal deformability effects on the signal. The detection of a primordial black hole would have implications for cosmology and dark matter scenarios, while the measurement of the tidal deformability of a subsolar neutron star could rule out or confirm the existence of strange stars made of quarks.

Supermassive Black Holes (BHs) are known to efficiently grow through gas accretion, but even sustained and intense mass build-up through this mechanism struggles to explain the assembly of the most massive BHs observed in the local Universe. Using the Chandra Deep-Wide Field Survey (CDFWS) in the Boötes field, we measure BH-galaxy assembly in massive galaxies ($M_\star\gtrsim10^{10}\, \rm M_\odot$) through the AGN fraction and specific Black Hole accretion rate (sBHAR) distribution as a function of redshift and stellar mass. We determine stellar masses and star formation rates for a parent sample of optically selected galaxies as well as those with X-ray detections indicating the presence of an AGN through Spectral Energy Distribution (SED) fitting. We derive a redshift-dependent mass completeness limit and extract X-ray information for every galaxy as to identify AGN in various types of galaxies at low and high redshift providing a comprehensive picture of the AGN population in massive galaxies. While X-ray AGN samples are dominated by moderately massive host galaxies of $M_{\star} \geqslant 10^{10}\, \rm M_{\odot}$, we do not find a strong stellar mass dependence in AGN fraction (to limits in sBHAR) or sBHAR distribution above this mass threshold, indicating a bias towards massive galaxies in the observed samples. A mild increase in AGN fraction is seen towards high stellar mass for low sBHAR AGN, as well as a suppression of high sBHAR events in the most massive galaxies. We derive BH-galaxy growth tracks over time, which reveal that while most BH mass has been accumulated since $z=4$ for lower mass BHs, the assembly of the most massive BHs is more complex, with little to no relative mass gain since $z=4$, implying that rapid and intense growth episodes prior to $z=4$ were necessary to form these massive BHs.

We present an analysis of the first two XRISM/Resolve spectra of the well-known Seyfert-1.5 active galactic nucleus in NGC 4151, obtained in December 2023. Our work focuses on the nature of the narrow Fe K$_{\alpha}$ emission line at 6.4 keV, the strongest and most common X-ray line observed in AGN. The total line is found to consist of three components. Even the narrowest component of the line is resolved with evident Fe K$_{\alpha,1}$ (6.404 keV) and K$_{\alpha,2}$ (6.391 keV) contributions in a 2:1 flux ratio, fully consistent with neutral gas with negligible bulk velocity. Subject to the limitations of our models, the narrowest and intermediate-width components are consistent with emission from optically thin gas, suggesting that they arise in a disk atmosphere and/or wind. Modeling the three line components in terms of Keplerian broadening, they are readily associated with (1) the inner wall of the ``torus,'' (2) the innermost optical ``broad line region'' (or, ``X-ray BLR''), and (3) a region with a radius of $r\simeq 100~GM/c^{2}$ that may signal a warp in the accretion disk. Viable alternative explanations of the broadest component include a fast wind component and/or scattering; however, we find evidence of variability in the narrow Fe K$_{\alpha}$ line complex on time scales consistent with small radii. The best-fit models are statistically superior to simple Voigt functions, but when fit with Voigt profiles the time-averaged lines are consistent with a projected velocity broadening of FWHM$=1600^{+400}_{-200}~{\rm km}~{\rm s}^{-1}$. Overall, the resolution and sensitivity of XRISM show that the narrow Fe K line in AGN is an effective probe of all key parts of the accretion flow, as it is currently understood. We discuss the implications of these findings for our understanding of AGN accretion, future studies with XRISM, and X-ray-based black hole mass measurements.

We present an initial analysis of the XRISM first-light observation of the supernova remnant (SNR) N132D in the Large Magellanic Cloud. The Resolve microcalorimeter has obtained the first high-resolution spectrum in the 1.6-10 keV band, which contains K-shell emission lines of Si, S, Ar, Ca, and Fe. We find that the Si and S lines are relatively narrow, with a broadening represented by a Gaussian-like velocity dispersion of $\sigma_v \sim 450$ km s$^{-1}$. The Fe He$\alpha$ lines are, on the other hand, substantially broadened with $\sigma_v \sim 1670$ km s$^{-1}$. This broadening can be explained by a combination of the thermal Doppler effect due to the high ion temperature and the kinematic Doppler effect due to the SNR expansion. Assuming that the Fe He$\alpha$ emission originates predominantly from the supernova ejecta, we estimate the reverse shock velocity at the time when the bulk of the Fe ejecta were shock heated to be $-1000 \lesssim V_{\rm rs}~[{\rm km s}^{-1}] \lesssim 3300$ (in the observer frame). We also find that Fe Ly$\alpha$ emission is redshifted with a bulk velocity of $\sim 890$ km s$^{-1}$, substantially larger than the radial velocity of the local interstellar medium surrounding N132D. These results demonstrate that high-resolution X-ray spectroscopy is capable of providing constraints on the evolutionary stage, geometry, and velocity distribution of SNRs.

The presence of the $\approx 10^6$ K gas in the circumgalactic medium of the Milky Way has been well established. However, the location and the origin of the newly discovered hot gas at `super-virial' temperatures of $\approx 10^7$ K have been puzzling. This hot gas has been detected in both absorption and emission; here we focus on the emitting gas only. We show that both the `virial' and the `super-virial' temperature gas as observed in \emph{emission} occupy disk-like extraplanar regions, in addition to the diffuse virial temperature gas filling the halo of the Milky Way. We perform idealized hydrodynamical simulations to show that the $\approx 10^7$ K emitting gas is likely to be produced by stellar feedback in and around the Galactic disk. We further show that the emitting gas at both super-virial and virial temperatures in the extraplanar regions is metal enriched and is not in hydrostatic equilibrium with the halo but is continuously evolving.

The recent detection of large column density absorption lines from highly ionized gas in a few directions through the circumgalactic medium (CGM) of the Milky Way (MW) has been puzzling. The inferred temperature from these absorption lines far exceeds the virial temperature of the MW, and the column densities are also too large to be easily explained. In this paper, we propose a novel idea to explain these observations and claim that they may not have originated from the CGM, but from a totally different type of source, namely, stellar ejecta from supernovae (SNe) above the Galactic disk that happen to lie in the line of sight to the background quasars. About $\sim 20\%$ of massive OB stars (progenitors of core-collapse supernovae) are known to be runaway stars that have high ejection velocities near the Galactic plane and can end up exploding as SNe above the Galactic disk. We show that the associated reverse shock in the supernova remnant in the early non-radiative phase can heat the ejecta to temperatures of $\gtrsim 10^7\,{\rm K}$ and can naturally explain the observed high column density of ions in the observed `super-virial' phase along with $\alpha$-enriched super-solar abundance that is typical of core-collapse supernovae. However, SNe from runaway stars has a covering fraction of $\lesssim 0.7 \%$ and thus can only explain the observations along limited sightlines.

Current observational evidence reveals that fast radio bursts (FRBs) exhibit bandwidths ranging from a few dozen MHz to several GHz. Traditional FRB searches primarily employ matched filter methods on time series collapsed across the entire observational bandwidth. However, with modern ultra-wideband receivers featuring GHz-scale observational bandwidths, this approach may overlook a significant number of events. We investigate the efficacy of sub-banded searches for FRBs, a technique seeking bursts within limited portions of the bandwidth. These searches aim to enhance the significance of FRB detections by mitigating the impact of noise outside the targeted frequency range, thereby improving signal-to-noise ratios. We conducted a series of Monte Carlo simulations, for the $400$-MHz bandwidth Parkes 21-cm multi-beam (PMB) receiver system and the Parkes Ultra-Wideband Low (UWL) receiver, simulating bursts down to frequency widths of about $100$\,MHz. Additionally, we performed a complete reprocessing of the high-latitude segment of the High Time Resolution Universe South survey (HTRU-S) of the Parkes-Murriyang telescope using sub-banded search techniques. Simulations reveal that a sub-banded search can enhance the burst search efficiency by $67_{-42}^{+133}$ % for the PMB system and $1433_{-126}^{+143}$ % for the UWL receiver. Furthermore, the reprocessing of HTRU led to the confident detection of eighteen new bursts, nearly tripling the count of FRBs found in this survey. These results underscore the importance of employing sub-banded search methodologies to effectively address the often modest spectral occupancy of these signals.

We have investigated the possibility of molecular cloud formation via the Collision-induced Magnetic Reconnection (CMR) mechanism of the cold neutral medium (CNM). Two atomic gas clouds with conditions typical of the CNM were set to collide at the interface of reverse magnetic fields. The cloud-cloud collision triggered magnetic reconnection and produced a giant 20pc filamentary structure which was not seen in the control models without CMR. The cloud, with rich fiber-like sub-structures, developed a fully molecular spine at 5Myr. Radiative transfer modeling of dust emission at far infrared wavelengths showed that the middle part of the filament contained dense cores over a span of 5pc. Some of the cores were actively forming stars and typically exhibited both connecting fibers in dust emission and high-velocity gas in CO line emission, indicative of active accretion through streamers. Supersonic turbulence was present in and around the CMR-filament due to inflowing gas moving at supersonic velocities in the collision mid-plane. The shocked gas was condensed and transported to the main filament piece by piece by reconnected fields, making the filament and star formation a bottom-up process. Instead of forming a gravitationally bounded cloud which then fragments hierarchically (top-down) and forms stars, the CMR process creates dense gas pieces and magnetically transports them to the central axis to constitute the filament. Since no turbulence is manually driven, our results suggest that CMR is capable of self-generating turbulence. Finally, the resulting helical field should show field-reversal on both sides of the filament from most viewing angles.

Non-LTE radiative transfer is a key tool for modern astrophysics: it is the means by which many key synthetic observables are produced, thus connecting simulations and observations. Radiative transfer models also inform our understanding of the primary formation layers and parameters of different spectral lines, and serve as the basis of inversion tools used to infer the structure of the solar atmosphere from observations. The default approach for computing the radiation field in multidimensional solar radiative transfer models has long remained the same: a short characteristics, discrete ordinates method, formal solver. In situations with complex atmospheric structure and multiple transitions between optically-thick and -thin regimes these solvers require prohibitively high angular resolution to correctly resolve the radiation field. Here, we present the theory of radiance cascades, a technique designed to exploit structure inherent to the radiation field, allowing for efficient reuse of calculated samples, thus providing a very high-resolution result at a fraction of the computational cost of existing methods. We additionally describe our implementation of this method in the DexRT code, and present initial results of the synthesis of a snapshot of a magnetohydrodynamic model of a solar prominence formed via levitation-condensation. The approach presented here provides a credible route for routinely performing multidimensional radiative transfer calculations free from so-called ray effects, and scaling high-quality non-LTE models to next-generation high-performance computing systems with GPU accelerators.

We present the discovery of 118 new ultracool dwarf candidates, discovered using a new machine learning tool, named \texttt{SMDET}, applied to time series images from the Wide-field Infrared Survey Explorer. We gathered photometric and astrometric data to estimate each candidate's spectral type, distance, and tangential velocity. This sample has a photometrically estimated spectral class distribution of 28 M dwarfs, 64 L dwarfs, and 18 T dwarfs. We also identify a T subdwarf candidate, two extreme T subdwarf candidates, and two candidate young ultracool dwarfs. Five objects did not have enough photometric data for any estimations to be made. To validate our estimated spectral types, spectra were collected for 2 objects, yielding confirmed spectral types of T5 (estimated T5) and T3 (estimated T4). Demonstrating the effectiveness of machine learning tools as a new large-scale discovery technique.

Solar physicists routinely utilize observations of Ar-like Fe IX and Cl-like Fe X emission to study a variety of solar structures. However, unidentified lines exist in the Fe IX and Fe X spectra, greatly impeding the spectroscopic diagnostic potential of these ions. Here, we present measurements using the Lawrence Livermore National Laboratory EBIT-I electron beam ion trap in the wavelength range 238-258 A. These studies enable us to unambiguously identify the charge state associated with each of the observed lines. This wavelength range is of particular interest because it contains the Fe IX density diagnostic line ratio 241.74 A/244.91 A, which is predicted to be one of the best density diagnostics of the solar corona, as well as the Fe X 257.26 A magnetic-field-induced transition. We compare our measurements to the Fe IX and Fe X lines tabulated in CHIANTI v10.0.1, which is used for modeling the solar spectrum. In addition, we have measured previously unidentified Fe X lines that will need to be added to CHIANTI and other spectroscopic databases.

Pulsating stars in eclipsing binary systems offer a unique possibility to empirically identify pulsation modes using the geometric effect of the eclipses on the pulsation signals. Here we explore the $\delta$ Scuti type pulsations in the eclipsing binary system KIC~3858884 with the aim of identifying the dominant modes using various photometric methods. We used the \textit{Kepler} short cadence photometry data. Refined binary model and pulsation parameters were determined using an iterative separation of the eclipsing binary and pulsation signals. We used various methods to identify the host stars of the dominant pulsations. Échelle diagram diagnostics were employed to locate the frequencies with most influence from the eclipses. Direct Fitting methods assuming spherical harmonic surface patterns were explored to determine an orientation for the symmetry axis and to infer surface mode numbers $\ell$ and $m$. General surface patterns were reconstructed using dynamic eclipse mapping, and provided ancillary mode number estimates. We established the secondary star as the main source of the pulsations. Seven peaks, including the two strongest modes, were found to show modulations during the secondary eclipses. For the first time ever, we could detect two hidden modes with amplitude intensification during the eclipses. Only one frequency appears to originate from the primary. We successfully reconstructed surface patterns and determined mode numbers for most of the selected frequencies with both of our methods. One radial and three sectoral modes, among them hidden modes. The two hidden modes were identified as (3,$\pm$1) and (2,$\pm$1). One additional radial mode turned out to be a combination frequency. Partial disagreement between EM and DF results may indicate that the strongest modes deviate from strict spherical harmonics.

Cosmic demographics -- the statistical study of populations of astrophysical objects -- has long relied on *multivariate statistics*, providing methods for analyzing data comprising fixed-length vectors of properties of objects, as might be compiled in a tabular astronomical catalog (say, with sky coordinates, and brightness measurements in a fixed number of spectral passbands). But beginning with the emergence of automated digital sky surveys, ca. ~2000, astronomers began producing large collections of data with more complex structure: light curves (brightness time series) and spectra (brightness vs. wavelength). These comprise what statisticians call *functional data* -- measurements of populations of functions. Upcoming automated sky surveys will soon provide astronomers with a flood of functional data. New methods are needed to accurately and optimally analyze large ensembles of light curves and spectra, accumulating information both along and across measured functions. Functional data analysis (FDA) provides tools for statistical modeling of functional data. Astronomical data presents several challenges for FDA methodology, e.g., sparse, irregular, and asynchronous sampling, and heteroscedastic measurement error. Bayesian FDA uses hierarchical Bayesian models for function populations, and is well suited to addressing these challenges. We provide an overview of astronomical functional data, and of some key Bayesian FDA modeling approaches, including functional mixed effects models, and stochastic process models. We briefly describe a Bayesian FDA framework combining FDA and machine learning methods to build low-dimensional parametric models for galaxy spectra.