Papers for Monday, Aug 14 2023

In this work, we perform the re-calibration of PS1 photometry by correcting for position-dependent systematic errors using the spectroscopy-based Stellar Color Regression method (SCR), the photometry-based SCR method (SCR$'$), and the Gaia XP synthetic photometry method (XPSP). We confirm the significant large-scale and small-scale spatial variation of magnitude offsets for all the $grizy$ filters. We show that the PS1 photometric calibration precisions in the $grizy$ filters are around 5--7\,mmag when averaged over 14$'$ regions. We note a much larger calibration error up to 0.04 mag in the Galactic plane, which is probably caused by the systematic errors of the PS1 magnitudes in crowded fields. The results of the three methods are consistent with each other within 1--2\,mmag or better for all the filters. We provide two-dimensional maps and a python package ({\url{https://doi.org/10.12149/101283}}) to correct for position-dependent magnitude offsets of PS1, which can be used for high-precision investigations and as a reference to calibrate other surveys.

Simulations of the effects of stellar fly-bys on planetary systems in star-forming regions show a strong dependence on subtle variations in the initial spatial and kinematic substructure of the regions. For similar stellar densities, the more substructured star-forming regions disrupt up to a factor of two more planetary systems. We extend this work to look at the effects of substructure on stellar binary populations. We present $N$-body simulations of substructured, and non-substructured (smooth) star-forming regions in which we place different populations of stellar binaries. We find that for binary populations that are dominated by close ($<$100au) systems, a higher proportion are destroyed in substructured regions. However, for wider systems ($>$100au), a higher proportion are destroyed in smooth regions. The difference is likely due to the hard-soft, or fast-slow boundary for binary destruction. Hard (fast/close) binaries are more likely to be destroyed in environments with a small velocity dispersion (kinematically substructured regions), whereas soft (slow/wide) binaries are more likely to be destroyed in environments with higher velocity dispersions (non-kinematically substructured regions). Due to the vast range of stellar binary semimajor axes in star-forming regions ($10^{-2} - 10^4$au) these differences are small and hence unlikely to be observable. However, planetary systems have a much smaller initial semimajor axis range (likely $\sim$1 -- 100au for gas giants) and here the difference in the fraction of companions due to substructure could be observed if the star-forming regions that disrupt planetary systems formed with similar stellar densities.

An equivalent width (EW) based classification may cause the erroneous judgement to the flat spectrum radio quasars (FSRQs) and BL Lacerate objects (BL Lac) due to the diluting the line features by dramatic variations in the jet continuum flux. To help address the issue, the present paper explore the possible intrinsic classification on the bias of a random forest supervised machine learning algorithm. In order to do so, we compile a sample of 1680 Fermi blazars that have both gamma-rays and radio-frequencies data available from the 4LAC-DR2 catalog, which includes 1352 training and validation samples and 328 forecast samples. By studying the results for all of the different combinations of 23 characteristic parameters, we found that there are 178 optimal parameters combinations (OPCs) with the highest accuracy ($\simeq$ 98.89\%). Using the combined classification results from the nine combinations of these OPCs to the 328 forecast samples, we predict that there are 113 true BL Lacs (TBLs) and 157 false BL Lacs (FBLs) that are possible intrinsically FSRQs misclassified as BL Lacs. The FBLs show a clear separation from TBLs and FSRQs in the $\gamma$-ray photon spectral index, $\Gamma_{\rm ph}$, and X-band radio flux, ${\rm log}{F_{R}}$, plot. Phenomenally, existence a BL Lac to FSRQ (B-to-F) transition zone is suggested, where the FBLs are in the stage of transition from BL Lacs to FSRQs. Comparing the LSP Changing-Look Blazars (CLBs) reported in the literatures, the majority of LSP CLBs are located at the B-to-F zone. We argue that the FBLs located at B-to-F transition zone are the most likely Candidates of CLBs.

Large-scale outflows driven by supermassive black holes are thought to play a fundamental role in suppressing star formation in massive galaxies. However, direct observational evidence for this hypothesis is still lacking, particularly in the young universe where star formation quenching is remarkably rapid, thus requiring effective removal of gas as opposed to slow gas heating. While outflows of ionized gas are commonly detected in massive distant galaxies, the amount of ejected mass is too small to be able to suppress star formation. Gas ejection is expected to be more efficient in the neutral and molecular phases, but at high redshift these have only been observed in starbursts and quasars. Using deep spectroscopy from JWST, here we show the presence of an outflow of neutral and ionized gas in a massive galaxy observed during the rapid quenching of its star formation, at a redshift of z=2.445. The outflowing mass is mostly in the neutral phase, and the mass outflow rate is larger than the residual star formation rate, indicating that the gas ejection is likely to have a strong impact on the evolution of the galaxy. We do not detect X-ray or radio activity; however the presence of a supermassive black hole is suggested by the properties of the ionized gas emission lines. We thus conclude that supermassive black holes are able to rapidly suppress star formation in massive galaxies by efficiently ejecting neutral gas.

Deriving properties of cometary nuclei from coma data is of significant importance for our understanding of cometary activity and has implications beyond. Ground-based data represent the bulk of measurements available for comets. Yet, to date these observations only access a comet's gas and dust coma at rather large distances from the surface and do not directly observe its surface. In contrast, spacecraft missions are one of the only tools that gain direct access to surface measurements. However, such missions are limited to roughly one per decade. We can overcome these challenges by recognising that the coma contains information about the nucleus's properties. In particular, the near-surface gas environment is most representative of the nucleus. It can inform us about the composition, regionality of activity, and sources of coma features and how they link to the topography, morphology, or other surface properties. The inner coma data is a particularly good proxy because it has not yet, or only marginally, been contaminated by coma chemistry or secondary gas sources. Additionally, the simultaneous observation of the innermost coma with the surface provides the potential to make a direct link between coma measurements and the nucleus. If we hope to link outer coma measurements obtained by Earth-based telescopes to the surface, we must first understand how the inner coma measurements are linked to the surface. Numerical models that describe the flow from the surface into the immediate surroundings are needed to make this connection. This chapter focuses on the advances made to understand the flow of the neutral gas coma from the surface to distances up to a few tens of nuclei radii. We describe both simple/heuristic models and state-of-the-art physically consistent models. Model limitations and what they each are best suited for are discussed.

Spectra of the highest redshift galaxies taken with JWST are now allowing us to see into the heart of the reionization epoch. Many of these observed galaxies exhibit strong damping wing absorption redward of their Lyman-$\alpha$ emission. These observations have been used to measure the redshift evolution of the neutral fraction of the intergalactic medium and sizes of ionized bubbles. However, these estimates have been made using a simple analytic model for the intergalactic damping wing. We explore the recent observations with models of inhomogeneous reionization from the Sherwood-Relics simulation suite. We carry out a comparison between the damping wings calculated from the simulations and from the analytic model. We find that although the agreement is good on the red side of the Lyman-$\alpha$ emission, there is a discrepancy on the blue side due to residual neutral hydrogen present in the simulations, which saturates the intergalactic absorption. For this reason, we find that it is difficult to reproduce the claimed observations of large bubble sizes at z ~ 7, which are driven by a detection of transmitted flux blueward of the Lyman-$\alpha$ emission. We suggest instead that the observations can be explained by a model with smaller ionized bubbles and larger intrinsic Lyman-$\alpha$ emission from the host galaxy.

The Milky Way is thought to host a huge population of interstellar objects (ISOs), numbering approximately $10^{15}\mathrm{pc}^{-3}$ around the Sun, which are formed and shaped by a diverse set of processes ranging from planet formation to galactic dynamics. We define a novel framework: firstly to predict the properties of this Galactic ISO population by combining models of processes across planetary and galactic scales, and secondly to make inferences about the processes modelled, by comparing the predicted population to what is observed. We predict the spatial and compositional distribution of the Galaxy's population of ISOs by modelling the Galactic stellar population with data from the APOGEE survey and combining this with a protoplanetary disk chemistry model. Selecting ISO water mass fraction as an example observable quantity, we evaluate its distribution both at the position of the Sun and averaged over the Galactic disk; our prediction for the Solar neighbourhood is compatible with the inferred water mass fraction of 2I/Borisov. We show that the well-studied Galactic stellar metallicity gradient has a corresponding ISO compositional gradient. We also demonstrate the inference part of the framework by using the current observed ISO composition distribution to constrain the parent star metallicity dependence of the ISO production rate. This constraint, and other inferences made with this framework, will improve dramatically as the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) progresses and more ISOs are observed. Finally, we explore generalisations of this framework to other Galactic populations, such as that of exoplanets.

We present results of a wideband high-resolution polarization study of Hydra A, one of the most luminous FR I radio galaxies known and amongst the most well-studied. The radio emission from this source displays extremely large Faraday rotation measures (RM), ranging from -12300 rad m$^{-2}$ to 5000 rad m$^{-2}$, the majority of which are believed to originate from magnetized thermal gas external to the radio tails. The radio emission from both tails strongly depolarizes with decreasing frequency. The depolarization, as a function of wavelength, is commonly non-monotonic, often showing oscillatory behavior, with strongly non-linear rotation of the polarization position angle with $\lambda^2$. A simple model, based on the RM screen derived from the high frequency, high resolution data, predicts the lower frequency depolarization remarkably well. The success of this model indicates the majority of the depolarization can be attributed to fluctuations in the magnetic field on scales $< 1500$ pc, suggesting the presence of turbulent magnetic field/electron density structures on sub-kpc scales within a Faraday rotating (FR) medium.

In this paper, we have revisited a class of coupled dark energy models where dark energy interacts with dark matter via phenomenological interactions. We included correction terms on the perturbation equations taking into account the perturbation of the Hubble rate, absent in previous works. We also consider more recent data sets such as cosmic microwave background anisotropies from \textit{Planck} 2018, type I-a supernovae measurements from Pantheon+ and data from baryon acoustic oscillations, and redshift space distortions. We analyzed the influence of the SH0ES Cepheid host distances on the results and for one model the discrepancy of $H_0$ is reduced to $1.3\sigma$ when compared to $\Lambda$CDM and $4.6\sigma$ when compared to the SH0ES team.

In this paper we present details of the construction of a wideband, cryogenic receiver and its successful commissioning on the Arecibo 12m telescope. The cryogenic receiver works in the 2.5-14 GHz frequency range. The telescope is operated by the Arecibo Observatory, and is located within the premises of the Observatory. We upgraded the current narrow band, room temperature receivers of the telescope with the new wideband receiver. The current receiver is built around a Quadruple-Ridged Flared Horn (QRHF) developed by Akgiray et al. (2013). To mitigate strong radio frequency interference (RFI) below 2.7 GHz, we installed a highpass filter before the first stage low noise amplifier (LNA). The QRHF, highpass filter, noise coupler and LNA are located inside a cryostat and are cooled to 15 K. The measured receiver temperature is 25 K (median value) over 2.5 GHz to 14 GHz. The system temperature measured at zenith is about 40 K near 3.1 and 8.6 GHz and the zenith antenna gains are 0.025 and 0.018 K/Jy at the two frequencies respectively. In the next stage of the development, we plan to upgrade the highpass filter in order to achieve better RFI rejection near 2.5 GHz, improve the aperture efficiency at 8.6 GHz and upgrade the IF system to increase the upper frequency of operation from 12 GHz to 14 GHz.

Broad lines in active galaxies are primarily emitted by photoionization processes that are driven by the incident continuum arising from a complex geometrical structure circumscribing the supermassive black hole. A model of the broad-band spectral energy distribution (SED) effective in ionizing the gas-rich broad line emitting region (BLR) is needed to understand the various radiative processes that eventually lead to the emission of emission lines from diverse physical conditions. Photoionization codes are a useful tool to investigate two aspects - the importance of the shape of the SED, and the physical conditions in the BLR. In this work, we focus on the anisotropy of continuum radiation from the very centre a direct consequence of the development of a funnel-like structure at regions very close to the black hole. Accounting for the diversity of Type-1 active galactic nuclei (AGNs) in the context of the main sequence of quasars permits us to locate the super Eddington sources along the sequence and constrain the physical conditions of their line-emitting BLR.

NGC 5548 has been hailed as an archetypical type-1 active galactic nuclei (AGN) and serves as a valuable laboratory to study the long-term variation of its broad-line region (BLR). In this work, we re-affirm our finding on the connection between the continuum variability in the optical regime and the corresponding H$\beta$ response to it, in order to realize the increase, albeit with a gradual saturation, in the H$\beta$ emitting luminosity with increasing AGN continuum. This effect is also known as the Pronik-Chuvaev effect after the authors who first demonstrated this effect using long-term monitoring of another well-studied Type-1 AGN - Mrk 6. We employ a homogeneous, multi-component spectral fitting procedure over a broad range of spectral epochs that is then used to create the continuum and H$\beta$ light curves. We focus on the epoch range 48636-49686 MJD, different from our previous analysis. We again notice a clear signature of shallowing in the trend between the H$\beta$ and the continuum luminosities. We attempt to recover this H$\beta$ emission trend as a response to a significant continuum flux increase using CLOUDY photoionization simulations and employ a suitable broad-band spectral energy distribution for this source. We explore the wide range in the physical parameters space for modelling the H$\beta$ emission from the BLR appropriate for this source. We employ a constant density, single cloud model approach in this study and successfully recover the observed shallowing of the H$\beta$ emission with respect to the rising AGN continuum and provide constraints on the local BLR densities and the location of the H$\beta$ emitting BLR which agrees with the H$\beta$ time-lags reported from the long-term reverberation mapping monitoring. On the contrary, we do not find a significant breathing effect in the location of the H$\beta$ line-emitting BLR for this epoch in NGC 5548.

Emission-line studies in the active galactic nuclei (AGNs), particularly those utilizing high spatial resolution, provide the most accurate method to determine critical quantities of the central engine and of the gas a few tens of parsecs away. Using seeing-limited data with spectroscopy, we have explored the extended narrow-line region for a sample of active galactic nuclei (AGNs) with strong, forbidden emission lines that have high-ionization potentials (IP $\gtrsim$ 100 eV). We have studied the optical and near-infrared spectra for these AGNs, extracted and compared their spectral energy distributions, and put constraints on the physical conditions of the region producing the coronal lines. We have realized a novel black hole mass scaling relation with one such prominent coronal line - [Si VI] 1.963 microns, over the 10$^6$ - 10$^8$ M$_{\odot}$ interval, that suggests photoionization by the continuum produced by the accretion disk is the primary physical process at play here. We perform a detailed parameter space study to optimize the emission from these coronal lines in terms of fundamental black hole parameters and test predictions that can be used to measure the kinematics of the extended X-ray emission gas. With the successful launch and first light of the JWST, we are well-poised to refine our findings using the superb angular resolution of the telescope that will allow us to map the inner few parsecs to the central supermassive black holes. This opens up the study of the higher ionization lines that will be spatially resolved by JWST, expanding our sample to tens of hundreds of AGNs, and putting firmer constraints on the physical conditions in the coronal line region.

Solar radio bursts are strongly affected by radio-wave scattering on density inhomogeneities, changing their observed time characteristics, sizes, and positions. The same turbulence causes angular broadening and scintillation of galactic and extra-galactic compact radio sources observed through the solar atmosphere. Using large-scale simulations of radio-wave transport, the characteristics of anisotropic density turbulence from $0.1 \, R_\odot$ to $1$ au are explored. For the first time, a profile of heliospheric density fluctuations is deduced that accounts for the properties of extra-solar radio sources, solar radio bursts, and in-situ density fluctuation measurements in the solar wind at $1$ au. The radial profile of the spectrum-weighted mean wavenumber of density fluctuations (a quantity proportional to the scattering rate of radio-waves) is found to have a broad maximum at around $(4-7) \, R_\odot$, where the slow solar wind becomes supersonic. The level of density fluctuations at the inner scale (which is consistent with the proton resonance scale) decreases with heliocentric distance as $\langle\delta{n_i}^2 \rangle (r) \simeq 2 \times 10^7 \, (r/R_\odot-1)^{-3.7}$ cm$^{-6}$. Due to scattering, the apparent positions of solar burst sources observed at frequencies between $0.1$ and $300$ MHz are computed to be essentially cospatial and to have comparable sizes, for both fundamental and harmonic emission. Anisotropic scattering is found to account for the shortest solar radio burst decay times observed, and the required wavenumber anisotropy is $q_\parallel/q_\perp =0.25-0.4$, depending on whether fundamental or harmonic emission is involved. The deduced radio-wave scattering rate paves the way to quantify intrinsic solar radio burst characteristics.

Exoplanet imaging has thus far enabled studies of wide-orbit ($>$10 AU) giant planet ($>$2 Jupiter masses) formation and giant planet atmospheres, with future 30 meter-class Extremely Large Telescopes (ELTs) needed to image and characterize terrestrial exoplanets. However, current state-of-the-art exoplanet imaging technologies placed on ELTs would still miss the contrast required for imaging Earth-mass habitable-zone exoplanets around low-mass stars by ~100x due to speckle noise--scattered starlight in the science image due to a combination of aberrations from the atmosphere after an adaptive optics (AO) correction and internal to the telescope and instrument. We have been developing a focal plane wavefront sensing technology called the Fast Atmospheric Self-coherent camera Technique (FAST) to address both of these issues; in this work we present the first results of simultaneous first and second stage AO wavefront sensing and control with a Shack Hartmann wavefront sensor (SHWFS) and FAST, respectively, using two common path deformable mirrors. We demonstrate this "multi-WFS single conjugate AO" real-time control at up to 200 Hz loop speeds on the Santa Cruz Extreme AO Laboratory (SEAL) testbed, showing a promising potential for both FAST and similar high-speed diffraction-limited second-stage wavefront sensing technologies to be deployed on current and future observatories, helping to remove speckle noise as the main limitation to ELT habitable exoplanet imaging.

The MIDAS (Micro-Imaging Dust Analysis System) atomic force microscope on board the Rosetta comet orbiter investigated and measured the 3D topography of a few hundred nm to tens of $\mu$m sized dust particles of 67P/Churyumov-Gerasimenko with resolutions down to a few nanometers, giving insights into the physical processes of our early Solar System. We analyze the shapes of the cometary dust particles collected by MIDAS on the basis of a recently updated particle catalog with the aim to determine which structural properties remained pristine. We develop a set of shape descriptors and metrics such as aspect ratio, elongation, circularity, convexity, and particle surface/volume distribution, which can be used to describe the distribution of particle shapes. Furthermore, we compare the structure of the MIDAS dust particles and the clusters in which the particles were deposited to those found in previous laboratory experiments and by Rosetta/COSIMA. Finally, we combine our findings to calculate a pristineness score for MIDAS particles and determine the most pristine particles and their properties. We find that the morphological properties of all cometary dust particles at the micrometer scale are surprisingly homogeneous despite originating from diverse cometary environments (e.g., different collection targets that are associated with cometary activities/source regions and collection velocities/periods). We next find that the types of clusters found by MIDAS show good agreement with those defined by previous laboratory experiments, however, there are some differences to those found by Rosetta/COSIMA. Based on our result, we rate 19 out of 1082 MIDAS particles at least moderately pristine, i.e., they are not substantially flattened by impact, not fragmented, and/or not part of a fragmentation cluster.

In addition to the intrinsic clustering of galaxies themselves, the spatial distribution of galaxies observed in surveys is modulated by the presence of weak lensing due to matter in the foreground. This effect, known as magnification bias, is a significant contaminant to analyses of galaxy-lensing cross-correlations and must be carefully modeled. We present a method to estimate the magnification bias in spectroscopically confirmed galaxy samples based on finite differences of galaxy catalogs while marginalizing over errors due to finite step size. We use our estimator to measure the magnification biases of the CMASS and LOWZ samples in the SDSS BOSS galaxy survey, analytically taking into account the dependence on galaxy shape for fiber and PSF magnitudes, finding $\alpha_{\rm CMASS} = 2.71 \pm 0.02$ and $\alpha_{\rm LOWZ} = 2.45 \pm 0.02$ and quantify modeling uncertainties in these measurements. Finally, we quantify the redshift evolution of the magnification bias within the CMASS and LOWZ samples, finding a difference of up to a factor of three between the lower and upper redshift bounds for the former. We discuss how to account for this evolution in modeling and its interaction with commonly applied redshift-dependent weights. Our method should be readily-applicable to upcoming surveys and we make our code publicly available as part of this work.

The first JWST observations of TRAPPIST-1 c showed a secondary eclipse depth of 421+/-94 ppm at 15 um, which is consistent with a bare rock surface or a thin, O2-dominated, low CO2 atmosphere (Zieba et al. 2023). Here, we further explore potential atmospheres for TRAPPIST-1 c by comparing the observed secondary eclipse depth to synthetic spectra of a broader range of plausible environments. To self-consistently incorporate the impact of photochemistry and atmospheric composition on atmospheric thermal structure and predicted eclipse depth, we use a two-column climate model coupled to a photochemical model, and simulate O2-dominated, Venus-like, and steam atmospheres. We find that a broader suite of plausible atmospheric compositions are also consistent with the data. For lower pressure atmospheres (0.1 bar), our O2-CO2 atmospheres produce eclipse depths within 1$\sigma$ of the data, consistent with the modeling results of Zieba et al. (2023). However, for higher-pressure atmospheres, our models produce different temperature-pressure profiles and are less pessimistic, with 1-10 bar O2, 100 ppm CO2 models within 2.0-2.2$\sigma$ of the measured secondary eclipse depth, and up to 0.5% CO2 within 2.9$\sigma$. Venus-like atmospheres are still unlikely. For thin O2 atmospheres of 0.1 bar with a low abundance of CO2 ($\sim$100 ppm), up to 10% water vapor can be present and still provide an eclipse depth within 1$\sigma$ of the data. We compared the TRAPPIST-1 c data to modeled steam atmospheres of $\leq$ 3 bar, which are 1.7-1.8$\sigma$ from the data and not conclusively ruled out. More data will be required to discriminate between possible atmospheres, or to more definitively support the bare rock hypothesis.

We use the dispersion measure (DM) and redshift measurements of 24 localized fast radio bursts (FRBs) to compare cosmological models and investigate the Hubble tension. Setting a flat prior on the DM contribution from the Milky Way's halo, $\mathrm{DM_{halo}^{MW}}\in[5,\;80]\;\mathrm{pc\;cm^{-3}}$, the best fit for flat $\Lambda$CDM is obtained with a Hubble constant $H_0=95.8^{+7.8}_{-9.2}\;\mathrm{km\;s^{-1}\;Mpc^{-1}}$ and a median matter density $\Omega_{\mathrm{m}}\approx0.66$. The best fit for the $R_{\mathrm{h}}=ct$ universe is realized with $H_0=94.2^{+5.6}_{-6.2}\;\mathrm{km\;s^{-1}\;Mpc^{-1}}$. We emphasize that the $H_0$ measurement depends sensitively on the $\mathrm{DM_{halo}^{MW}}$ prior. Since flat $\Lambda$CDM has one more free parameter, $R_{\mathrm{h}}=ct$ is favored by the Bayesian Information Criterion (BIC) with a likelihood of $\sim73\%$ versus $\sim27\%$. Through simulations, we find that if the real cosmology is $\Lambda$CDM, a sample of $\sim1,150$ FRBs in the redshift range $0<z<3$ would be sufficient to rule out $R_{\mathrm{h}}=ct$ at a $3\sigma$ confidence level, while $\sim550$ FRBs would be necessary to rule out $\Lambda$CDM if the real cosmology is instead $R_{\mathrm{h}}=ct$. The required sample sizes are different, reflecting the fact that the BIC imposes a severe penalty on the model with more free parameters. We further adopt a straightforward method of deriving an upper limit to $H_{0}$, without needing to consider the poorly known probability distribution of the DM contributed by the host galaxy. The theoretical DM contribution from the intergalactic medium ($\mathrm{DM_{IGM}}$) at any $z$ is proportional to $H_0$. Thus, requiring the extragalactic $\mathrm{DM_{ext}}$ to be larger than $\mathrm{DM_{IGM}}$ delimits $H_0$ to the upside. Assuming flat $\Lambda$CDM, we have $H_0<89.0\;\mathrm{km\;s^{-1}\;Mpc^{-1}}$ at a 95\% confidence level.

Complex Fe-K emission/absorption line features are commonly observed in the 6--11 keV band from Active Galactic Nuclei (AGN). These features are formed in various physical components surrounding the black holes. The Narrow-Line Seyfert 1 (NLS1) galaxy Mrk 766, in particular, exhibits characteristic blue-shifted Fe-K absorption lines caused by the ultra-fast outflow (UFO), and a broad Fe-K emission line, as well as variable absorbers partially covering the X-ray emitting region. We re-analyze the Mrk 766 archival data of XMM-Newton, NuSTAR, and Swift to investigate the origin of the Fe-K line feature and the 0.3--79 keV energy spectral variation. We have found that the spectral variation in $\lesssim$10 keV is primarily explained by the variable partial covering of the central X-ray source by multi-layer absorbing clouds. The Fe-K line feature consists of the blue-shifted absorption lines due to the UFO, a narrow emission line from the distant material, a broad emission line from the inner-disk reflection, and a slightly broadened weak emission line at around 6.4--6.7 keV whose equivalent width is $\sim$0.05 keV. The last one is presumably due to the resonance scattering in the UFO out of the line-of-sight, as predicted by a Monte Carlo simulation based on the hydrodynamical UFO modeling. We suggest that the seemingly complex Fe-K line features and the X-ray energy spectra of Mrk 766 are explained by a moderately extended central X-ray source around a Schwarzschild black hole, an optically thick accretion disk with a truncated inner-radius, the UFO, multi-layer partial covering clouds, and a torus.

Context: Polycyclic Aromatic Hydrocarbons, largely known as PAHs, are widespread in the universe and have been identified in a vast array of astronomical observations from the interstellar medium to protoplanetary discs. They are likely to be associated with the chemical history of the universe and the emergence of life on Earth. However, their abundance on exoplanets remains unknown. Aims: We aim to investigate the feasibility of PAH formation in the thermalized atmospheres of irradiated and non-irradiated hot Jupiters around Sun-like stars. Methods: To this aim, we introduced PAHs in the 1-D self-consistent forward modeling code petitCODE. We simulated a large number of planet atmospheres with different parameters (e.g. carbon to oxygen ratio, metallicity, and effective planetary temperature) to study PAH formation. By coupling the thermochemical equilibrium solution from petitCODE with the 1-D radiative transfer code, petitRADTRANS, we calculated the synthetic transmission and emission spectra for irradiated and non-irradiated planets, respectively, and explored the role of PAHs on planet spectra. Results: Our models show strong correlations between PAH abundance and the aforementioned parameters. In thermochemical equilibrium scenarios, an optimal temperature, elevated carbon to oxygen ratio, and increased metallicity values are conducive to the formation of PAHs, with the carbon to oxygen ratio having the largest effect.

The presence of methanol among the common ice components in interstellar clouds and protostellar envelopes has been confirmed by the James Webb Space Telescope. Methanol is often detected in the gas phase toward lines of sight shielded from UV radiation. We measured the volumetric density of methanol ice, grown under simulated interstellar conditions, and the infrared spectroscopy at different deposition temperatures and during the warm-up. The IR band strengths are provided and the experimental spectra are compared to those computed with a model. The transition from amorphous to crystalline methanol ice was also explored. Finally, we propose new observations of methanol ice at high resolution to probe the methanol ice structure.

We report detections of scintillation arcs for pulsars in globular clusters M5, M13 and M15 for the first time using the Five-hundred-meter Aperture Spherical radio Telescope (FAST). From observations of these arcs at multiple epochs, we infer that screen-like scattering medium exists at distances $4.1_{-0.3}^{+0.2}$ kpc, $6.7_{-0.2}^{+0.2}$ kpc and $1.3_{-1.0}^{+0.7}$ kpc from Earth in the directions of M5, M13 and M15, respectively. This means M5's and M13's scattering screens are located at $3.0_{-0.2}^{+0.1}$ kpc and $4.4_{-0.1}^{+0.1}$ kpc above the galactic plane, whereas, M15's is at $0.6_{-0.5}^{+0.3}$ kpc below the plane. We estimate the scintillation timescale and decorrelation bandwidth for each pulsar at each epoch using the one-dimensional auto-correlation in frequency and time of the dynamic spectra. We found that the boundary of the Local Bubble may have caused the scattering of M15, and detected the most distant off-plane scattering screens to date through pulsar scintillation, which provides evidence for understanding the medium circulation in the Milky Way.

This paper is written as a follow-up observations to reinterpret the radial velocity (RV) of HD 36384, where the existence of planetary systems is known to be ambiguous. In giants, it is, in general, difficult to distinguish the signals of planetary companions from those of stellar activities. Thus, known exoplanetary giant hosts are relatively rare. We, for many years, have obtained RV data in evolved stars using the high-resolution, fiber-fed Bohyunsan Observatory Echelle Spectrograph (BOES) at the Bohyunsan Optical Astronomy Observatory (BOAO). Here, we report the results of RV variations in the M giant HD 36384. We have found two significant periods of 586d and 490d. Considering the orbital stability, it is impossible to have two planets at so close orbits. To determine the nature of the RV variability variations, we analyze the HIPPARCOS photometric data, some indicators of stellar activities, and line profiles. A significant period of 580d was revealed in the HIPPARCOS photometry. H{\alpha} EW variations also show a meaningful period of 582d. Thus, the period of 586d may be closely related to the rotational modulations and/or stellar pulsations. On the other hand, the other significant period of 490d is interpreted as the result of the orbiting companion. Our orbital fit suggests that the companion was a planetary mass of 6.6 MJ and is located at 1.3 AU from the host.

The minimum initial mass required for a star to explode as an Fe core collapse supernova, typically denoted $M_\text{mas}$, is an important quantity in stellar evolution because it defines the border between intermediate mass and massive stellar evolutionary paths. The precise value of $M_\text{mas}$ carries implications for models of galactic chemical evolution and the calculation of star formation rates. Despite the fact that stars with super solar metallicities are commonplace within spiral and some giant elliptical galaxies, there are currently no studies of this mass threshold in super metal-rich models with $Z>0.05$. Here, we study the minimum mass necessary for a star to undergo an Fe core collapse supernova when its initial metal content falls in the range $2.5\times 10^{-3} \leq Z \leq 0.10$. Although an increase in initial $Z$ corresponds to an increase in the Fe ignition threshold for $Z \approx 1\times 10^{-3}$ to $Z\approx0.04$, we find that there is a steady reversal in trend that occurs for $Z > 0.05$. Our super metal-rich models thus undergo Fe core collapse at lower initial masses than those required at solar metallicity. Our results indicate that metallicity--dependent curves extending to $Z=0.10$ for the minimum Fe ignition mass should be utilised in galactic chemical evolution simulations to accurately model supernovae rates as a function of metallicity, particularly for simulations of metal-rich spiral and elliptical galaxies.

The interior of Enceladus, a medium sized icy moon of Saturn hosts hydrothermal activity and exhibits tidal heating and related geyser-like activity. There are major disagreements in the existing literature on the porosity of the interior, due to the different theoretical assumptions on which porosity related calculations were based. We present an application of experimental equations - derived for Earth - for icy planetary objects and Enceladus in particular. We chose a set of boundary values for our initial parameters from measured porosity values of chondrite samples as references, and calculated the porosity related values of Enceladus using various approaches. We present a comprehensive investigation of the effects of using these different porosity calculation methods on icy moons. With our most realistic approach we also calculated the same values for Earth and Mars for comparison. Our result for Enceladus is a minimum porosity of about 5\% at the centre of the body. For the total pore volume we estimated $1.51*10^7 km^3$ for Enceladus, $2.11*10^8 km^3$ for Earth and $1.62*10^8 km^3$ for Mars. Using the same method, we estimated the total pore surface area. From this we derived that the pore surface under a given $1 km^2$ area of the surface on Enceladus is about $1.37*10^9 km^2$, while for Earth this value is only $5.07*10^7 km^2$.

We present optical photometric and spectroscopic observations of the 02es-like type Ia supernova (SN) 2022ywc. The transient occurred in the outskirts of an elliptical host galaxy and showed a striking double-peaked light curve with an early excess feature detected in the ATLAS orange and cyan bands. The early excess is remarkably luminous with an absolute magnitude $\sim -19$, comparable in luminosity to the subsequent radioactively-driven second peak. The spectra resemble the hybrid 02es-like SN 2016jhr, that is considered to be a helium shell detonation candidate. We investigate different physical mechanisms that could power such a prominent early excess and rule out massive helium shell detonation, surface $^{56}$Ni distribution and ejecta-companion interaction. We conclude that SN ejecta interacting with circumstellar material (CSM) is the most viable scenario. Semi-analytical modelling with MOSFiT indicates that SN ejecta interacting with $\sim 0.05\,$M$_{\odot}$ of CSM at a distance of $\sim 10^{14}$ cm can explain the extraordinary light curve. A double-degenerate scenario may explain the origin of the CSM, either by tidally-stripped material from the secondary white dwarf, or disk-originated matter launched along polar axes following the disruption and accretion of the secondary white dwarf. A non-spherical CSM configuration could suggest that a small fraction of 02es-like events viewed along a favourable line of sight may be expected to display a very conspicuous early excess like SN 2022ywc.

Neutron stars accreting from massive binary companions come in a wide range of types. Systems with an OB supergiant donor are often divided between persistently and transiently accreting systems, respectively called Supergiant X-ray Binaries (SgXBs) and Supergiant Fast X-ray Transients (SFXTs). The origin of this dichotomy in accretion behaviour is typically attributed to systematic differences in the massive stellar wind, the binary orbit, or magnetic field configuration, but direct observational evidence for these hypotheses remains sparse. Here, we present the results of a pilot exploration of a novel approach to this long-standing question, turning to the mm band to probe the outer regions of the stellar wind beyond the binary orbit. Specifically, we present 100-GHz NOEMA observations of a SgXB, X1908+075, and a SFXT, IGR J18410-0535. We detect the SFXT as a point source at $63.4 \pm 9.6$ $\mu$Jy, while the SgXB is not detected. The spectrum of IGR J18410-0535 is constrained to be flat or inverted by comparing with quasi-simultaneous $5.5$+$9$ GHz radio observations, ruling out non-thermal flaring and consistent with thermal wind emission. Additional X-ray measurements further constrain the wind mass loss rate and velocity of the SgXB. We compare our targets with each other and earlier wind estimates, and reflect on future opportunities using this novel observational approach to characterize stellar winds in X-ray binaries.

Context. The long-term evolution of an atmosphere and the remote detectability of its chemical constituents are susceptible to how the atmospheric gas responds to stellar irradiation. The response remains poorly characterized for water and its dissociation products, however, this knowledge is relevant to our understanding of hypothetical water-rich exoplanets. Aims: Our work investigates the effect of photoelectrons, namely, the non-thermal electrons produced by photoionizing stellar radiation on the heating and ionization of extended atmospheres dominated by the dissociation products of water. Methods: We used a Monte Carlo model and up-to-date collision cross sections to simulate the slowing down of photoelectrons in O-H mixtures for a range of fractional ionizations and photoelectron energies. Results: We find that that the fraction of energy of a photoelectron that goes into heating is similar in a pure H gas and in O-H mixtures, except for very low fractional ionizations, whereby the O atom remains an efficient sink of energy. The O-H mixtures will go on to produce more electrons because the O atom is particularly susceptible to ionization. We quantified all that information and present it in a way that can be easily incorporated into photochemical-hydrodynamical models. Conclusions: Neglecting the role of photoelectrons in models of water-rich atmospheres will result in overestimations of the atmospheric heating and, foreseeably, the mass-loss rates as well. It will also underestimate the rate at which the atmospheric gas becomes ionized, which may have implications for the detection of extended atmospheres with Lyman-{\alpha} transmission spectroscopy. Our simulations for the small exoplanets {\pi} Men c and TRAPPIST-1 b reveal that they respond very differently to irradiation from their host stars, with water remaining in molecular form at lower pressures in the latter case.

The modelling of astrophysical systems such as binary neutron star mergers or the formation of magnetars from the collapse of massive stars involves the numerical evolution of magnetised fluids at extremely large Reynolds numbers. This is a major challenge for (unresolved) direct numerical simulations which may struggle to resolve highly dynamical features as, e.g. turbulence, magnetic field amplification, or the transport of angular momentum. Sub-grid models offer a means to overcome those difficulties. In a recent paper we presented MInIT, an MHD-instability-induced-turbulence mean-field, sub-grid model based on the modelling of the turbulent (Maxwell, Reynolds, and Faraday) stress tensors. While in our previous work MInIT was assessed within the framework of the magnetorotational instability, in this paper we further evaluate the model in the context of the Kelvin-Helmholtz instability (KHI). The main difference with other sub-grid models (as e.g. the alpha-viscosity model or the gradient model) is that in MInIT we track independently the turbulent energy density at sub-grid scales, which is used, via a simple closure relation, to compute the different turbulent stresses relevant for the dynamics. The free coefficients of the model are calibrated using well resolved box simulations of magnetic turbulence generated by the KHI. We test the model against these simulations and show that it yields order-of-magnitude accurate predictions for the evolution of the turbulent Reynolds and Maxwell stresses.

The detection of gamma-ray emission from accreting pulsars in X-ray binaries (XRBs) has long been sought after. For some high-mass X-ray binaries (HMXBs), marginal detections have recently been reported. Regardless of whether these will be confirmed or not, future telescopes operating in the gamma-ray band could offer the sensitivity needed to achieve solid detections and possibly spectra. In view of future observational advances, we explored the expected emission above 10 GeV from XRBs, based on the Cheng & Ruderman model, where gamma-ray photons are produced by the decay of pion-0 originated by protons accelerated in the magnetosphere of an accreting pulsar fed by an accretion disc. We improved this model by considering, through Monte Carlo simulations, the development of cascades inside of and outside the accretion disc, taking into account pair and photon production processes that involve interaction with nuclei, X-ray photons from the accretion disc, and the magnetic field. We produced grids of solutions for different input parameter values of the X-ray luminosity (L_x), magnetic field strength (B), and for different properties of the region where acceleration occurs. We found that the gamma-ray luminosity spans more than five orders of magnitude, with a maximum of ~1E35 erg/s. The gamma-ray spectra show a large variety of shapes: some have most of the emission below ~100 GeV, others are harder (emission up to 10-100 TeV). We compared our results with Fermi/LAT and VERITAS detections and upper-limits of two HMXBs: A0535+26 and GROJ1008-57. More consequential comparisons will be possible when more sensitive instruments will be operational in the coming years.

The enigmatic nature of 55 Cancri e has defied theoretical explanation. Any explanation needs to account for the observed variability of its secondary eclipse depth, which is at times consistent with zero in the visible/optical range of wavelengths -- a phenomenon that does not occur with its also variable infrared eclipses. Yet, despite this variability its transit depth remains somewhat constant in time and is inconsistent with opaque material filling its Hill sphere. The current study explores the possibility of a thin, transient, secondary atmosphere on 55 Cancri e that is sourced by geochemical outgassing. Its transient nature derives from the inability of outgassing to be balanced by atmospheric escape. As the outgassed atmosphere escapes and is replenished, it rapidly adjusts to radiative equilibrium and the temperature fluctuations cause the infrared eclipse depths to vary. Atmospheres of pure carbon dioxide or carbon monoxide produce sufficient Rayleigh scattering to explain the observed optical/visible eclipse depths, which vanish in the absence of an atmosphere and the presence of a dark rocky surface. Atmospheres of pure methane are ruled out because they produce insufficient Rayleigh scattering. Upcoming observations by the James Webb Space Telescope will potentially allow the atmospheric temperature and surface pressure, as well as the surface temperature, to be measured.

This paper presents new optical and near-UV spectra of 11 extremely powerful jetted quasars, with radio to optical flux density ratio $>$ 10$^3$, that concomitantly cover the low-ionization emission of \mgii\ and \hb\ as well as the \feii\ blends in the redshift range $0.35 \lesssim z \lesssim 1$. We aim to quantify broad emission line differences between radio-loud (RL) and radio-quiet (RQ) quasars by using the 4D eigenvector 1 parameter space and its Main Sequence (MS) and to check the effect of powerful radio ejection on the low ionization broad emission lines. The \hb\ and \mgii\ emission lines were measured by using non-linear multicomponent fittings as well as by analysing their full profile. We found that broad emission lines show large redward asymmetry both in \hb\ and \mgii. The location of our RL sources in a UV plane looks similar to the optical one, with weak \feiiuv\ emission and broad \mgii. We supplement the 11 sources with large samples from previous work to gain some general inferences. We found that, compared to RQ, our extreme RL quasars show larger median \hb\ full width at half maximum (FWHM), weaker \feii\ emission, larger \mbh, lower \lledd, and a restricted bf space occupation in the optical and UV MS planes. The differences are more elusive when the comparison is carried out by restricting the RQ population to the region of the MS occupied by RL sources, albeit an unbiased comparison matching \mbh\ and \lledd\ suggests that the most powerful RL quasars show the highest redward asymmetries in \hb.

We derive a minimal basis of kernels furnishing the perturbative expansion of the density contrast and velocity divergence in powers of the initial density field that is applicable to cosmological models with arbitrary expansion history, thereby relaxing the commonly adopted Einstein-de-Sitter (EdS) approximation. For this class of cosmological models, the non-linear kernels are at every order given by a sum of terms, each of which factorizes into a time-dependent growth factor and a wavenumber-dependent basis function. We show how to reduce the set of basis functions to a minimal amount, and give explicit expressions up to order $n=5$. We find that for this minimal basis choice, each basis function individually displays the expected scaling behaviour due to momentum conservation, being non-trivial at $n\geq 4$. This is a highly desirable property for numerical evaluation of loop corrections. In addition, it allows us to match the density field to an effective field theory (EFT) description for cosmologies with an arbitrary expansion history, which we explicitly derive at order four. We evaluate the differences to the EdS approximation for $\Lambda$CDM and $w_0w_a$CDM, paying special attention to the irreducible cosmology dependence that cannot be absorbed into EFT terms for the one-loop bispectrum. Finally, we provide algebraic recursion relations for a special generalization of the EdS approximation that retains its simplicity and is relevant for mixed hot and cold dark matter models.

Since the release of the Gravitational Wave Transient Catalogue GWTC-2.1 by the LIGO-Virgo collaboration, sub-threshold gravitational wave (GW) candidates are publicly available. They are expected to be released in real-time as well, in the upcoming O4 run. Using these GW candidates for multi-messenger studies complement the ongoing efforts to identify neutrino counterparts to GW events. This in turn, allows us to schedule electromagnetic follow-up searches more efficiently. However, the definition and criteria for sub-threshold candidates are pretty flexible. Finding a multi-messenger counterpart via archival studies for these candidates will help to set up strong bounds on the GW parameters which are useful for defining a GW signal as sub-threshold, thereby increasing their significance for scheduling follow-up searches. Here, we present the current status of this ongoing work with the IceCube Neutrino Observatory. We perform a selection of the sub-threshold GW candidates from GWTC-2.1 and conduct an archival search for sub-TeV neutrino counterparts detected by the dense infill array of the IceCube Neutrino Observatory, known as "DeepCore". For this, an Unbinned Maximum Likelihood (UML) method is used. We report the 90% C.L. sensitivities of this sub-TeV neutrino dataset for each selected sub-threshold GW candidate, considering the spatial and temporal correlation between the GW and neutrino events within a 1000 s time window.

The origin of extremely fast variability is one of the long-standing questions in the gamma-ray astronomy of blazars. While many models explain the slower, lower energy variability, they cannot easily account for such fast flares reaching hour-to-minute time scales. Magnetic reconnection, a process where magnetic energy is converted to the acceleration of relativistic particles in the reconnection layer, is a candidate solution to this problem. In this work, we employ state-of-the-art particle-in-cell simulations in a statistical comparison with observations of a flaring episode of a well-known blazar, Mrk 421, at very high energy (VHE, E > 100 GeV). We tested the predictions of our model by generating simulated VHE light curves that we compared quantitatively with methods that we have developed for a precise evaluation of theoretical and observed data. With our analysis, we can constrain the parameter space of the model, such as the magnetic field strength of the unreconnected plasma, viewing angle and the reconnection layer orientation in the blazar jet. Our analysis favours parameter spaces with magnetic field strength 0.1 G, rather large viewing angles (6-8 degrees), and misaligned layer angles, offering a strong candidate explanation for the Doppler crisis often observed in the jets of high synchrotron peaking blazars.

Short-duration gamma-ray bursts (sGRBs) are commonly attributed to the mergers of double neutron stars (NSs) or the mergers of a neutron star with a black hole (BH). While the former scenario was confirmed by the event GW170817, the latter remains elusive. Here, we consider the latter scenario in which, a NS is tidally disrupted by a fast spinning low-mass BH and the accretion onto the BH launches a relativistic jet and hence produces a sGRB. The merging binary's orbit is likely misaligned with the BH's spin. Hence, the Lense-Thirring precession around the BH may cause a hyper-accreting thick disk to precess in a solid-body manner. We propose that a jet, initially aligned with the BH spin, is deflected and collimated by the wind from the disk, therefore being forced to precess along with the disk. This would result in a quasi-periodic oscillation or modulation in the gamma-ray light curve of the sGRB, with a quasi-period of $\sim 0.01-0.1$ s. The appearance of the modulation may be delayed respective to the triggering of the light curve. This feature, unique to the BH-NS merger, may have already revealed itself in a few observed sGRBs (such as GRB 130310A), and it carries the spin-obit orientation information of the merging system. Identification of this feature would be a new approach to reveal spin-orbit-misaligned merging BH-NS systems, which are likely missed by the current gravitational-wave searching strategy principally targeting aligned systems.

X-ray binaries (XRBs) are thought to regulate cosmic thermal and ionisation histories during the Epoch of Reionisation and Cosmic Dawn ($z\sim 5-30$). Theoretical predictions of the X-ray emission from XRBs are important for modeling such early cosmic evolution. Nevertheless, the contribution from Be-XRBs, powered by accretion of compact objects from decretion disks around rapidly rotating O/B stars, has not been investigated systematically. Be-XRBs are the largest class of high-mass XRBs (HMXBs) identified in local observations and are expected to play even more important roles in metal-poor environments at high redshifts. In light of this, we build a physically motivated model for Be-XRBs based on recent hydrodynamic simulations and observations of decretion disks. Our model is able to reproduce the observed population of Be-XRBs in the Small Magellanic Cloud with appropriate initial conditions and binary stellar evolution parameters. We derive the X-ray output from Be-XRBs as a function of metallicity in the (absolute) metallicity range $Z\in [10^{-4},0.03]$. We find that Be-XRBs can contribute a significant fraction ($\sim 60\%$) of the total X-ray budget from HMXBs observed in nearby galaxies for $Z\sim 0.0003-0.02$. A similar fraction of observed ultra-luminous ($\gtrsim 10^{39}\ \rm erg\ s^{-1}$) X-ray sources can also be explained by Be-XRBs. Moreover, the predicted metallicty dependence in our fiducial model is consistent with observations, showing a factor of $\sim 8$ increase in X-ray luminosity per unit star formation rate from $Z=0.02$ to $Z=0.0003$.

Context. Understanding complex phenomena and unsolved problems in modern astronomy requires wider-bandwidth observations. The current technique for designing and fabricating an astronomical instrument potentially provides such observations with higher efficiency and precision than in the past. Higher-order modes in an instrument associated with wider bandwidths have been reported, which may degrade observation precision. Aims. To reduce the unfavorable degradation, we need to quantify the higher-order propagation modes, though their power is too difficult to measure directly. Instead of the direct mode measurement, we aim at developing a method based on measurable radiation patterns from an instrument of interest. Method. Assuming a linear system, whose radiated field is determined as a superposition of the mode coefficients in an instrument, we obtain a coefficient matrix connecting the inside modes and the outside radiated field and calculate the pseudo-inverse matrix. To understand the estimation accuracy of the proposed method, we demonstrate two cases with numerical simulations, axially-corrugated horn case and offset Cassegrain antenna case, and investigate the effect of random errors on the accuracy. Results. Both cases showed the estimated mode coefficients with a precision of 10e-6 with respect to the maximum mode amplitude and 10e-3 degrees in phase, respectively. The calculation errors were observed when the random errors were smaller than 0.01 percent of the maximum radiated field amplitude. The demonstrated method works independently of the details of a system. Conclusions. The method can quantify the propagation modes inside an instrument and will be applicable to most of linear components and antennas. This method can be employed for a general purpose, such as diagnosis of feed alignment and higher-performance feed design.

Aims. The evolutionary status of the blue supergiant Sher 25 and its membership to the massive cluster NGC 3603 are investigated. Methods. A hybrid non-LTE (local thermodynamic equilibrium) spectrum synthesis approach is employed to analyse a high-resolution optical spectrum of Sher 25 and five similar early B-type comparison stars in order to derive atmospheric parameters and elemental abundances. Fundamental stellar parameters are determined by considering stellar evolution tracks, Gaia Data Release 3 (DR3) data and complementary distance information. Interstellar reddening and the reddening law along the sight line towards Sher 25 are constrained employing UV photometry for the first time in addition to optical and infrared data. The distance to NGC 3603 is reevaluated based on Gaia DR3 data of the innermost cluster O-stars. Results. The spectroscopic distance derived from the quantitative analysis implies that Sher 25 lies in the foreground of NGC 3603, which is found to have a distance of $d_\mathrm{NGC 3603}$ = 6250$\pm$150 pc. A cluster membership is also excluded as the hourglass nebula is unaffected by the vigorous stellar winds of the cluster stars and from the different excitation signatures of the hourglass nebula and the nebula around NGC 3603. Sher 25 turns out to have a luminosity of log L/L$_\odot$ = 5.48$\pm$0.14, equivalent to that of a $\sim$27 $M_\odot$ supergiant in a single-star scenario, which is about half of the mass assumed so far, bringing it much closer in its characteristics to Sk-69{\deg}202, the progenitor of SN 1987A. Sher 25 is significantly older than NGC 3603. Further arguments for a binary (merger) evolutionary scenario of Sher 25 are discussed.

The spin parameter of the black hole in the accreting X-ray binary LMC X-1 has been measured in a number of studies to be $a_*\gtrsim 0.9$. These measurements were claimed to take into account both statistical and systematic (model-dependent) uncertainties. We perform new measurements using a recent simultaneous observation of LMC X-1 by NICER and NuSTAR, providing a data set of very high quality. We use the disk continuum method together with improved models for coronal Comptonization. With the standard relativistic blackbody disk model and optically thin Comptonization, we obtain values of $a_*$ similar to those obtained before. We then consider modifications to the standard disk model. Using a phenomenological color correction set to 2, we find lower values of $a_*\approx 0.64$--0.84. We then consider disks thicker than the standard one, i.e., with some dissipation in surface layers, as expected if partially supported by magnetic pressure. To account for that, we assume the disk is covered by a warm and optically thick layer, Comptonizing the emission of the underlying disk. Our model with the lowest $\chi^2$ yields then a low range of the spin, $a_*\approx 0.40^{+0.41}_{-0.32}$. That last model is also in agreement with the inverse disk temperature-luminosity relation found in this source. We conclude that determinations of the spin using the continuum method is highly sensitive to the assumptions about the disk structure.

Recent associations of high-energy neutrinos with active galactic nuclei (AGN) have revived the interest in leptohadronic models of radiation from astrophysical sources. The rapid increase in the amount of acquired multi-messenger data will require soon fast numerical models that may be applied to large source samples. We develop a time-dependent leptohadronic code, LeHaMoC, that offers several notable benefits compared to other existing codes, such as versatility and speed. LeHaMoC solves the Fokker-Planck equations of photons and relativistic particles (i.e. electrons, positrons, protons, and neutrinos) produced in a homogeneous magnetized source that may also be expanding. The code utilizes a fully implicit difference scheme that allows fast computation of steady-state and dynamically evolving physical problems. We first present test cases where we compare the numerical results obtained with LeHaMoC against exact analytical solutions and numerical results computed with ATHE$\nu$A, a well-tested code of similar philosophy but different numerical implementation. We find a good agreement (within 10-30%) with the numerical results obtained with ATHE$\nu$A without evidence of systematic differences. We then demonstrate the capabilities of the code through illustrative examples. First, we fit the spectral energy distribution from a jetted AGN in the context of a synchrotron-self Compton model and a proton-synchrotron model using Bayesian inference. Second, we compute the high-energy neutrino signal and the electromagnetic cascade induced by hadronic interactions in the corona of NGC 1068. LeHaMoC is easily customized to model a variety of high-energy astrophysical sources and has the potential to become a widely utilized tool in multi-messenger astrophysics.

Messier 104, NGC 4594, also known as the Sombrero Galaxy, has been extensively studied, especially its structure and stellar halo. Its abundance of globular clusters has given rise to many theories and much speculation (Ford H. C. et al 1996). But other objects in the vicinity of such a spectacular galaxy are sometimes ignored. While studying HST images available on the HST Legacy website of the halo of M104 (HST proposal 9714, PI: Keith Noll), the author observed at 12:40:07.829 -11:36:47.38 (in j2000) an object about 4 arcseconds in diameter. A study with VO tools suggests that the object is a SBc galaxy with AGN (Seyfert).

Open clusters are excellent tools to probe the history of the Galactic disk and properties of star formation. In this work, we present a study of an old age open cluster Berkley 39 using the observations from UVOT instrument of the Neil Gehrels Swift observatory. Making use of a machine learning algorithm, ML-MOC, we have identified a total of 861 stars as cluster members out of which 17 are blue straggler stars. In this work, we present a characterisation of 2 blue straggler stars. To estimate the fundamental parameters of blue straggler stars and their companions (if any), we constructed spectral energy distributions using UV data from swift/UVOT and GALEX, optical data from Gaia DR3, and infrared (IR) data from 2MASS, Spitzer/IRAC, and WISE. We find excess flux in UV in one blue straggler star, implying the possibility of a hot companion.

We perform a multi-tracer full-shape analysis in Fourier space based on the effective field theory of large-scale structure (EFTofLSS) using the complete Sloan Digital Sky Survey IV (SDSS-IV) extended Baryon Oscillation Spectroscopic Survey (eBOSS) DR16 luminous red galaxy (LRG) and emission line galaxy (ELG) samples. We study in detail the impact of the volume projection effect and different prior choices when doing the full-shape analysis based on the EFTofLSS model. We show that adopting a combination of Jeffreys prior and Gaussian prior can mitigate the volume effect and avoid exploring unphysical regions in the parameter space at the same time, which is crucial when jointly analysing the eBOSS LRG and ELG samples. We validate our pipeline using 1000 eBOSS EZmocks. By performing a multi-tracer analysis on mocks with comparable footprints, we find that cosmological constraints can be improved by $\sim10-40\%$ depending on whether we assume zero stochastic terms in the cross power spectrum, which breaks the degeneracy and boosts the constraints on the standard deviation of matter density fluctuation $\sigma_8$. Combining with the Big Bang Nucleosynthesis (BBN) prior and fixing the spectral tilt $n_s$ to Planck value, our multi-tracer full-shape analysis measures $H_0=70.0\pm2.3~{\rm km}~{\rm s}^{-1} {\rm Mpc}^{-1}$, $\Omega_m=0.317^{+0.017}_{-0.021}$, $\sigma_8=0.787_{-0.062}^{+0.055}$ and $S_8=0.809_{-0.078}^{+0.064}$, consistent with the Planck 2018 results. In particular, the constraint on $\sigma_8$ is improved beyond that obtained from the single tracer analysis by $18\%$, or by $27\%$ when assuming zero stochastic terms in the cross power spectrum.

How multiple close-in super-Earths form around stars with masses lower than that of the Sun is still an open issue. Several recent modeling studies have focused on planet formation around M-dwarf stars, but so far no studies have focused specifically on K dwarfs, which are of particular interest in the search for extraterrestrial life. We aim to reproduce the currently known population of close-in super-Earths observed around K-dwarf stars and their system characteristics. We performed 48 high-resolution N-body simulations of planet formation via planetesimal accretion using the existing GENGA software running on GPUs. In the simulations we varied the initial disk mass and the solid and gas surface density profiles. Each simulation began with 12000 bodies with radii of between 200 and 2000 km around two different stars, with masses of 0.6 and 0.8 $M_{\odot}$. Most simulations ran for 20 Myr, with several simulations extended to 40 or 100 Myr. The mass distributions for the planets with masses between 2 and 12 $M_\oplus$ show a strong preference for planets with masses $M_p<6$ $M_\oplus$ and a lesser preference for planets with larger masses, whereas the mass distribution for the observed sample increases almost linearly. However, we managed to reproduce the main characteristics and architectures of the known planetary systems and produce mostly long-term angular-momentum-deficit-stable, nonresonant systems, but we require an initial disk mass of 15 $M_\oplus$ or higher and a gas surface density value at 1 AU of 1500 g cm$^{-2}$ or higher. Our simulations also produce many low-mass planets with $M<2$ $M_\oplus$, which are not yet found in the observed population, probably due to the observational biases. The final systems contain only a small number of planets, which could possibly accrete substantial amounts of gas, and these formed after the gas had mostly dissipated.

Currently, an upgrade consisting of seven densely instrumented strings in the center of the volume of the IceCube detector with new digital optical modules (DOMs) is being built. On each string, DOMs will be regularly spaced with a vertical separation of 3 m between depths of 2160 m and 2430 m below the surface of the ice, which is a denser configuration than the existing DOMs of IceCube detector. For a precise calibration of the IceCube Upgrade it is important to understand the properties of the ice, both inside and surrounding the deployment holes. The camera system together with the LED illumination system was developed and produced at Sungkyunkwan university and are installed in almost every DOM to measure these properties. For these calibration measurements, a new simulation framework, which produces expected images from various geometric and optical variables has been developed. Images produced from the simulation will be used to develop an analysis framework for the IceCube Upgrade camera calibration system and for the design of the IceCube Gen2 camera system.

We present new optical transmission spectra for two hot Jupiters: WASP-25b (M = 0.56~M$_J$; R = 1.23 R$_J$; P =~3.76 days) and WASP-124b (M = 0.58~M$_J$; R = 1.34 R$_J$; P = 3.37 days), with wavelength coverages of 4200 - 9100\AA\ and 4570 - 9940\AA, respectively. These spectra are from the ESO Faint Object Spectrograph and Camera (v.2) mounted on the New Technology Telescope (NTT) and Inamori-Magellan Areal Camera & Spectrograph on Magellan Baade. No strong spectral features were found in either spectra, with the data probing 4 and 6 scale heights, respectively. \texttt{Exoretrievals} and \texttt{PLATON} retrievals favor stellar activity for WASP-25b, while the data for WASP-124b did not favor one model over another. For both planets the retrievals found a wide range in the depths where the atmosphere could be optically thick ($\sim0.4\mu$ - 0.2 bars for WASP-25b and 1.6 $\mu$ -- 32 bars for WASP-124b) and recovered a temperature that is consistent with the planets' equilibrium temperatures, but with wide uncertainties (up to $\pm$430$^\circ$K). For WASP-25b, the models also favor stellar spots that are $\sim$500-3000$^\circ$K cooler than the surrounding photosphere. The fairly weak constraints on parameters are owing to the relatively low precision of the data, with an average precision of 840 and 1240 ppm per bin for WASP-25b and WASP-124b, respectively. However, some contribution might still be due to an inherent absence of absorption or scattering in the planets' upper atmospheres, possibly because of aerosols. We attempt to fit the strength of the sodium signals to the aerosol-metallicity trend proposed by McGruder et al. 2023, and find WASP-25b and WASP-124b are consistent with the prediction, though their uncertainties are too large to confidently confirm the trend.

Magnetic nulls are locations where the magnetic field vanishes. Nulls are the location of magnetic reconnection, and they determine to a large degree the magnetic connectivity in a system. We describe a novel approach to understanding movement, appearance, and disappearance of nulls in magnetic fields. This approach is based on the novel concept of isotropes, or lines where the field direction is constant. These lines are streamlines of a vector field whose flux is sourced by the topological indices of nulls, and can be conceptualized as corresponding "lines of force" between nulls. We show how this topological approach can be used to generate analytical expressions for the location of nulls in the presence of external fields for dipoles and for a field defined by the Hopf fibration.

We answer frequently asked questions (FAQs) about the Hellings and Downs correlation curve -- the "smoking-gun" signature that pulsar timing arrays (PTAs) have detected gravitational waves (GWs). Many of these questions arise because of intuition based on how ground-based interferometers like LIGO respond to GWs. These have arms that are short (km scale) compared to the wavelengths of the GWs that they detect (hundreds to thousands of km). In contrast, PTAs respond to GWs whose wavelengths (tens of light-years) are much shorter than their arms (a typical PTA pulsar is thousands of light-years from Earth). To elucidate this, we calculate the exact response of a "one-arm, one-way" detector to a passing GW, and compare it in the "short-arm" (LIGO-like) and "long-arm" (PTA) limits. This provides qualitative and quantitative answers to many questions about the Hellings and Downs correlation. The resulting "FAQ sheet" should help in understanding the "evidence for GWs" recently announced by several PTA collaborations.

Periodic patterns of logarithmic oscillations can arise in primordial curvature perturbation spectra and in the associated stochastic gravitational-wave background via different mechanisms. We show that, in the presence of log oscillations, the spectral shape of the stochastic background has a unique parametrization independent of its physical origin. We also show that this log-periodic modulation can be generated in any scenario beyond Einstein gravity endowed with an approximate discrete scale invariance, a symmetry typical of deterministic fractal spacetimes that could emerge in quantum gravity under certain conditions. We discuss how a log-oscillatory spectral shape arises from concrete inflationary models beyond Einstein gravity and the prospects for detection in Einstein Telescope and other next-generation gravitational-wave observatories. We find that these instruments will be sensitive to log-periodic features if the detection is made with a high signal-to-noise ratio (SNR) and that the error scales as $1/{\rm SNR}$.

In this work, we explore the intriguing possibility of connecting self-interacting dark matter (SIDM) with the recently observed exceptionally bright and long-duration Gamma Ray Burst (GRB221009A). The proposed minimal scenario involves a light scalar mediator, simultaneously enabling dark matter (DM) self-interaction and explaining the observed very high energy (VHE) photons from GRB221009A reported by LHAASO's data. The scalar's mixing with the standard model (SM) Higgs boson allows for its production at the GRB site, which will then propagate escaping attenuation by the extra-galactic background light (EBL). These scalars, if highly boosted, have the potential to explain LHAASO's data. Moreover, the same mixing also facilitates DM-nucleon or DM-electron scatterings at terrestrial detectors, linking SIDM phenomenology to the GRB221009A events. This manuscript presents the parameter space meeting all constraints and offers an exciting opportunity to explore SIDM in future direct search experiments using insights from the GRB observation.