Papers for Wednesday, Feb 28 2024

In this paper, we derive observational constraints on an anisotropic $w$CDM model from observational data including Baryonic Acoustic Oscillations (BAOs), Cosmic Chronometer (CC), Big Bang Nucleosynthesis (BBN), Pantheon Plus (PP) compilation of Type Ia supernovae, and SH0ES Cepheid host distance anchors. We find that anisotropy is of the order $10^{-13}$, and its presence in the $w$CDM model reduces $H_0$ tension by $\sim 2\sigma$ and $\sim 1\sigma$ in the analyses with BAO+CC+BBN+PP and BAO+CC+BBN+PPSH0ES data combinations, respectively. In both analyses, the quintessence form of dark energy is favored at 95\% CL.

We report the discovery of Cerberus, an extremely red object detected with the MIRI Deep Imaging Survey (MIDIS) observations in the F1000W filter of the Hubble Ultra Deep Field. The object is detected at S/N~6, with F1000W~27 mag, and it is extremely faint in both the NIRCam data gathered by the JWST Advanced Deep Extragalactic Survey, JADES, with ~30.5 mag $5\sigma$ upper limits in individual bands, as well as in the MIDIS F560W ultra deep data ($\sim$29 mag, $5\sigma$). Analyzing the spectral energy distribution built with individual (low S/N) optical-to-mid-infrared filters and (S/N~5) stacks, we discuss the possible nature of this red NIRCam-dark source using a battery of codes. We discard the possibility of Cerberus being a Solar System body based on the $<$0.016" proper motion in the 1-year apart JADES and MIDIS observations. A sub-stellar Galactic nature is deemed unlikely, given that the Cerberus' relatively flat NIRCam-to-NIRCam and very red NIRCam-to-MIRI flux ratios are not consistent with any brown dwarf model. The extragalactic nature of Cerberus offers 3 possibilities: (1) A $z\sim0.4$ galaxy with strong emission from polycyclic aromatic hydrocarbons; the very low inferred stellar mass, $\mathrm{M}_\star=10^{5-6}$ M$_\odot$, makes this possibility highly improbable. (2) A dusty galaxy at $z\sim4$ with an inferred stellar mass $\mathrm{M}_\star\sim10^{8}$ M$_\odot$. (3) A galaxy with observational properties similar to those of the reddest little red dots discovered around $z\sim7$, but Cerberus lying at $z\sim15$, presenting a spectral energy distribution in the rest-frame optical dominated by emission from a dusty torus or a dusty starburst.

Radio jets are present in a diverse sample of AGN. However, the mechanisms of jet powering are not fully understood, and it is yet unclear to what extent they obey mass-invariant scaling relations, similar to those found for the triggering and fuelling of X-ray selected AGN. We study the incidence of eROSITA/eFEDS X-ray and LOFAR radio AGN as a function of several stellar mass normalised AGN power indicators. A new sample of radio AGN from the LOFAR-eFEDS survey is defined and we publicly release this catalogue, including host galaxy counterparts from the Legacy Survey DR9, LOFAR radio morphologies and host galaxy properties from the complete, spectroscopic (z<0.4) GAMA09 survey. The fraction of GAMA09 galaxies hosting radio, X-ray and both radio and X-ray AGN are calculated as a function of the specific black hole kinetic ($\lambda_{\rm Jet}$) and radiative ($\lambda_{\rm Edd}$) power. The incidence of eFEDS X-ray AGN as a function of $\lambda_{\rm Edd}$ shows the same mass-invariance as found in past studies. Meanwhile, radio AGN, regardless of their morphology, are more likely to be hosted in more massive galaxies, at all $\lambda_{\rm Jet}$. Across the stellar mass range, the compact radio AGN incidence follows the same power-law distribution, showing that it is not only high mass galaxies that host high power radio AGN and vice versa. On the other hand, the incidence of compact and complex radio AGN is boosted at the highest jet powers, diverging from a simple power-law. Interestingly, this increased incidence cannot be explained by more powerful radio AGN lying in more dense environments which could naturally boost their radio luminosity. Overall, we show that statistical incidence studies are a powerful method to probe disk-jet coupling for different AGN accretion modes, although future work on a more reliable determination of jet power for diverse samples of radio AGN is needed.

Fuzzy Dark Matter (FDM) comprised of ultralight ($m \sim 10^{-22}~\rm{eV}$) boson particles has received significant attention as a viable alternative to Cold Dark Matter (CDM), as it approximates CDM on large scales ($\gtrsim 1$ Mpc) while potentially resolving some of its small-scale problems via kiloparsec-scale quantum interference. However, the most basic FDM model, with one free parameter (the boson mass), is subject to a tension: small boson masses yield the desired cores of dwarf galaxies but underpredict structure in the Lyman-$\alpha$ forest, while large boson masses render FDM effectively identical to CDM. This Catch-22 problem may be alleviated by considering an axion-like particle with attractive particle self-interactions. We simulate an idealized FDM halo with self-interactions parameterized by an energy decay constant $f \sim 10^{15}~\rm{GeV}$ related to the axion symmetry-breaking conjectured to solve the strong-CP problem in particle physics. We observe solitons, a hallmark of FDM, condensing within a broader halo envelope, and find that the density profile and soliton mass depend on self-interaction strength. We propose generalized formulae to extend those from previous works to include self-interactions. We also investigate a critical mass threshold predicted for strong interactions at which the soliton collapses into a compact, unresolved state. We find that the collapse happens quickly and its effects are initially contained to the central region of the halo.

We detect a highly significant excess of X-ray (2RXS) and radio (NVSS, GMRT, VLSSr) catalog sources when stacked around MCXC galaxy clusters and groups, narrowly confined within $\lesssim100\mathrm{\,kpc}$ of the $\sim2.4 R_{500}$ virial shock radius (inferred from continuum stacking), with similar X-ray ($\sim4\sigma$ for $443$ clusters) and radio ($\sim4\sigma$ for $485$ clusters) characteristics ($>5\sigma$ joint). The excess sources show $10-100$ kpc scales, $L_X(0.1-2.4\mbox{ keV})\simeq10^{42-43}\mathrm{\,erg\,s^{-1}}$ or $\nu L_\nu(\nu=1.4\mbox{ GHz}) \simeq 10^{40-41}\mathrm{\,erg\,s^{-1}}$ luminosities, and a preferentially radial radio-polarization. The narrow localization and properties of the excess identify these sources not as AGN, often invoked speculatively for excess X-ray sources at cluster outskirts, but rather as infalling gaseous clumps interacting with the virial shock, probably galactic halos and possibly outflow remnants. The local excess of such discrete, radio-to-$\gamma$-ray sources around an object can probe its virial shock also at high redshifts and sub-cluster scales.

There are indications that stellar-origin black holes (BHs) are efficiently paired up in binary black holes (BBHs) in Active Galactic Nuclei (AGN) disc environments, which can undergo interactions with single BHs in the disc. Such binary-single interactions can potentially lead to an exceptionally high fraction of gravitational-wave mergers with measurable eccentricity in LIGO/Virgo/KAGRA. We here take the next important step in this line of studies, by performing post-Newtonian N-body simulations between migrating BBHs and single BHs set in an AGN disc-like configuration with a consistent inclusion of the central supermassive black hole (SMBH) in the equations of motion. With this setup, we study how the fraction of eccentric mergers varies in terms of the initial size of the BBH semi-major axis relative to the Hill sphere, as well as how it depends on the angle between the BBH and the incoming single BH. We find that the fraction of eccentric mergers is still relatively large, even when the interactions are notably influenced by the gravitational field of the nearby SMBH. However, the fraction as a function of the BBH semi-major axis does not follow a smooth functional shape, but instead shows strongly varying features that originate from the underlying phase-space structure. The phase-space further reveals that many of the eccentric mergers are formed through prompt scatterings. Finally, we present the first analytical solution to how the presence of an SMBH in terms of its Hill sphere affects the probability for forming eccentric BBH mergers through chaotic three-body interactions.

Strong high-ionization iron coronal lines (CLs) are a rare phenomenon observed in galaxy and quasi-stellar object spectra that are thought to be created as a result of tidal disruption event (TDE) flares. To test whether these CLs are the result of TDE activity, we search for extreme coronal line emitting galaxies (ECLEs) in the Sloan Digital Sky Survey (SDSS), measure their rate, and compare it to TDE rates from the literature. We detect sufficiently strong CLs in 14 objects, doubling the number previously found in SDSS. Using follow-up spectra from the Dark Energy Spectroscopic Instrument and Gemini Multi-Object Spectrograph, Wide-field Infrared Survey Explorer mid-infrared observations, and Liverpool Telescope optical photometry, we find that of the seven new objects, only one evolves in a manner consistent with that of the five previously discovered variable ECLEs. Using this new sample of six variable ECLEs, we calculate the galaxy-normalised rate of ECLEs in SDSS to be $R_\mathrm{G}=2.2~^{+1.3}_{-0.8}~(\mathrm{statistical})~^{+0.0}_{-1.3}~(\mathrm{systematic})\times10^{-5}~\mathrm{galaxy}^{-1}~\mathrm{year}^{-1}$. The mass-normalised rate is $R_\mathrm{M}=1.9~^{+1.1}_{-0.7}~(\mathrm{statistical})~^{+0.0}_{-1.1}~(\mathrm{systematic})\times10^{-16}~\mathrm{M_\odot^{-1}}~\mathrm{year}^{-1}$ and the volumetric rate is $R_\mathrm{V}=6.9~^{+5.6}_{-2.1}~(\mathrm{statistical})~^{+0.0}_{-3.9}~(\mathrm{systematic})\times10^{-8}~\mathrm{Mpc}^{-3}~\mathrm{year}^{-1}$. Our rates are comparable to TDE rates from the literature, supporting the suggestion that the CLs in variable ECLEs are the product of TDEs.

The elemental abundances of stars, particularly the refractory elements (e.g., Fe, Si, Mg), play an important role in connecting stars to their planets. Most Sun-like stars do not have refractory abundance measurements since obtaining a large sample of high-resolution spectra is difficult with oversubscribed observing resources. In this work we infer abundances for C, N, O, Na, Mn, Cr, Si, Fe, Ni, Mg, V, Ca, Ti, Al, and Y for solar analogs with Gaia RVS spectra (R=11,200) using the Cannon, a data-driven method. We train a linear model on a reference set of 34 stars observed by Gaia RVS with precise abundances measured from previous high resolution spectroscopic efforts (R > 30,000--110,000). We then apply this model to several thousand Gaia RVS solar analogs. This yields abundances with average upper limit precisions of 0.04--0.1 dex for 17,412 stars, 50 of which are identified planet (candidate) hosts. We subsequently test the relative refractory depletion of these stars with increasing element condensation temperature compared to the Sun. The Sun remains refractory depleted compared to other Sun-like stars regardless of our current knowledge of the planets they host. This is inconsistent with theories of various types of planets locking up or sequestering refractories. Furthermore, we find no significant abundance differences between identified close-in giant planet hosts, giant planet hosts, and terrestrial/small planet hosts and the rest of the sample within our precision limits. This work demonstrates the utility of data-driven learning for future exoplanet composition and demographics studies.

We report the discovery of 23 globular cluster (GC) candidates around the relatively isolated dwarf galaxy IC 2574 within the Messier 81 (M81) group, at a distance of 3.86 Mpc. We use observations from the HST Advanced Camera for Surveys (ACS) to analyse the imaging in the F814W and F555W broadband filters. Our GC candidates have luminosities ranging from $-5.9 \geq M_V \geq -10.4$ and half-light radii of $1.4 \leq r_h \leq 11.5$ pc. We find the total number of GCs ($N_{\mathrm{GC}})=27\pm5$ after applying completeness corrections, which implies a specific frequency of $S_N = 4.0\pm0.8$, consistent with expectations based on its luminosity. The GC system appears to have a bimodal colour distribution, with 30% of the GC candidates having redder colours. We also find 5 objects with extremely blue colours that could be young star clusters linked to an intense star formation episode that occurred in IC 2574 $\sim$1 Gyr ago. We make an independent measurement of the halo mass of IC 2574 from its kinematic data, which is rare for low mass galaxies, and find log $M_{200} = 10.93 \pm 0.08$. We place the galaxy on the well-known GC system mass-halo mass relation and find that it agrees well with the observed near-linear relation. IC 2574 has a rich GC population for a dwarf galaxy, which includes an unusually bright $\omega$ Cen-like GC, making it an exciting nearby laboratory for probing the peculiar efficiency of forming massive GCs in dwarf galaxies.

We investigate a possible close encounter between M33 and M31 in the past to understand the role of galaxy-galaxy interactions in shaping the matter distribution in galaxy outskirts. We recovered possible orbital trajectories of M33, M31 and the Milky Way in the past, which are compatible with the Early Third Data Release of the Gaia mission and with mass estimates of Local Group spirals, after tuning mass losses and the dynamical friction term with the help of N-body numerical simulations. A close encounter of M33 and M31 in the past has a low but non-negligible probability. If the two galaxies had been closer in the past, their minimum distance would be of the order of 100 kpc or larger, and this happened earlier than 3 Gyr ago. During this encounter, 35-40% of the dark matter mass of M33 might have been removed from the halo due to tidal stripping. A detailed comparison of the results of test-particle simulations with the observed disk warp or with the spatial distribution of candidate dark satellites of M33 suggests that a closer passage of M33 around M31 cannot, however, be responsible for the observed morphological features. We suggest that more recent gas accretion events, possibly from a cosmic filament, might cause the misalignment of the outer disk of M33 after the rapid inner disk formation.

We analyze 15 years of eclipse timings of the polar V808 Aur. The rapid ingress/egress of the white dwarf and bright accretion region provide timings as precise as a few tenths of a second for rapid cadence photometric data. We find that between 2015 and 2018, the eclipse timings deviated from a linear ephemeris by more than 30 s. The rapid timing change is consistent with the periastron passage of a planet in an eccentric orbit about the polar. The best fit orbital period is 11$\pm 1$ yr and we estimate a projected mass of $Msin(i)=6.8\pm 0.7$ Jupiter masses. We also show that the eclipse timings are correlated with the brightness of the polar with a slope of 1.1 s/mag. This is likely due to the change in the geometry of the accretion curtains as a function of the mass transfer rate in the polar. While an eccentric planet offers an excellent explanation to the available eclipse data for V808 Aur, proposed planetary systems in other eclipsing polars have often struggled to accurately predict future eclipse timings.

The nature of the progenitor systems and explosion mechanisms that give rise to Type Ia supernovae (SNe Ia) are still debated. The interaction signature of circumstellar material (CSM) being swept up by expanding ejecta can constrain the type of system from which it was ejected. Most previous studies have focused on finding CSM ejected shortly before the SN Ia explosion still residing close to the explosion site, resulting in short delay times until the interaction starts. We use a sample of 3627 SNe Ia from the Zwicky Transient Facility discovered between 2018 and 2020 and search for interaction signatures over 100 days after peak brightness. By binning the late-time light curve data to push the detection limit as deep as possible, we identify potential late-time rebrightening in 3 SNe Ia (SN 2018grt, SN 2019dlf, SN 2020tfc). The late-time detections occur between 550 and 1450 d after peak brightness, have mean absolute $r$-band magnitudes of -16.4 to -16.8 mag and last up to a few hundred days, significantly brighter than the late-time CSM interaction discovered in the prototype SN 2015cp. The late-time detections all occur within 0.8 kpc of the host nucleus and are not easily explained by nuclear activity, another transient at a similar sky position, or data quality issues. This suggests environment or specific progenitor characteristics playing a role in producing potential CSM signatures in these SNe Ia. By simulating the ZTF survey we estimate that <0.5 per cent of normal SNe Ia display late-time strong H $\alpha$-dominated CSM interaction. This is equivalent to an absolute rate of $8_{-4}^{+20}$ to $54_{-26}^{+91}$ Gpc$^{-3}$ yr$^{-1}$ assuming a constant SN Ia rate of $2.4\times10^{-5}$ Mpc$^{-3}$ yr$^{-1}$ for $z \leq 0.1$. Weaker interaction signatures, more similar to the strength seen in SN 2015cp, could be more common but are difficult to constrain with our survey depth.

Powerful outflows from active galactic nuclei (AGN) can significantly impact the gas reservoirs of their host galaxies. However, it is still unclear how these outflows can propagate from the very central regions of galaxies to their outskirts, and whether nuclear winds can be driven by and/or be responsible for drastic spectral transitions. In this work we test feedback propagation models on the case test of 2MASS 0918+2117 (2M0918), a z=0.149 X-ray variable AGN, which showed tentative evidence for nuclear ultra-fast outflows (UFOs) in a 2005 XMM-Newton observation. We also investigate whether UFOs can be related to the observed X-ray variability. We observed 2M0918 with XMM-Newton and NuSTAR in 2020 to confirm the presence and characterize the UFOs. We perform a kinematic analysis of the 2005 SDSS optical spectrum to reveal and measure the properties of galaxy- scale ionized outflows. Furthermore, we construct 20-year-long lightcurves of observed flux, line-of-sight column density, and intrinsic accretion rate from the spectra of the first 4 SRG/eROSITA all-sky surveys and archival observations from Chandra and XMM-Newton.We significantly detect UFOs with v$\sim$0.16c and galaxy-scale ionized outflows with velocities of $\sim$ 700 km/s. We also find that the drastic X-ray variability (factors >10) can be explained both in terms of variable obscuration and variable intrinsic luminosity.Comparing the energetics of the two outflow phases, 2M0918 is consistent with momentum-driven wind propagation. 2M0918 expands the sample of AGN with both UFOs and ionized gas winds from 5 to 6, and brings the sample of AGN hosting multiscale outflows to 19, contributing to a clearer picture of feedback physics. From the variations in accretion rate, column density, and ionization level of the obscurer, we propose a scenario that connects obscurers, an accretion enhancement, and the emergence of UFOs

Imaging spectropolarimetry is a new observing mode on the Advanced Camera for Surveys (ACS) aboard the Hubble Space Telescope (HST) that was commissioned in Cycle 30 and is available to HST observers starting in Cycle 31 (i.e., from 2023). It is a technique that is accessible from ground-based observatories, but the superb spatial resolution afforded by HST/ACS combined with the slitless nature of HST/ACS grism spectroscopy opens up the possibility of studying polarized extended emission in a way that is not currently possible even with Adaptive Optics facilities on the ground. This mode could help to study interesting targets including (but not limited to) QSOs, AGN and Radio Galaxies, ISM Dust Properties, Pre-Planetary Nebulae, Proto-Planetary and Debris Disks, Supernovae/Supernova Remnants, and Solar System objects. This research note presents the preliminary results from the calibration programs used to calibrate imaging spectropolarimetry on HST/ACS.

We report the discovery and characterisation of TIC 350842552 ("Zvrk"), an apparently isolated, rapidly-rotating ($P_\text{rot} \sim 99\ \mathrm{d}$) red giant observed by TESS in its Southern Continuous Viewing Zone. The star's fast surface rotation is independently verified by the use of p-mode asteroseismology, strong periodicity in TESS and ASAS-SN photometry, and measurements of spectroscopic rotational broadening. A two-component fit to APOGEE spectra indicates a coverage fraction of its surface features consistent with the amplitude of the photometric rotational signal. Variations in the amplitude of its photometric modulations over time suggest the evolution of its surface morphology, and therefore enhanced magnetic activity. We further develop and deploy new asteroseismic techniques to characterise radial differential rotation, and find weak evidence for rotational shear within Zvrk's convective envelope. This feature, in combination with such a high surface rotation rate, is incompatible with models of angular-momentum transport in single-star evolution. Spectroscopic abundance estimates also indicate a high lithium abundance, among other chemical anomalies. Taken together, all of these suggest a planet-ingestion scenario for the formation of this rotational configuration, various models for which we examine in detail.

Using the medium resolution spectrograph X-shooter, spectra of 235 OB and Wolf-Rayet (WR) stars in sub-solar metallicity environments have been secured. [...]This second paper focuses on the optical observations of 232 Magellanic Clouds targets. It describes the uniform reduction of the UVB (300 - 560 nm) and VIS (550 - 1020 nm) XShootU data as well as the preparation of advanced data products [...] . The data reduction of the raw data is based on the ESO CPL X-shooter pipeline. We paid particular attention to the determination of the response curves [...] We implemented slit-loss correction, absolute flux calibration, (semi-)automatic rectification to the continuum, and a correction for telluric lines. The spectra of individual epochs were corrected for the barycentric motion, re-sampled and co-added, and the spectra from the two arms were merged into a single flux calibrated spectrum covering the entire optical range with maximum signal-to-noise ratio. [...] We provide three types of data products: (i) two-dimensional spectra for each UVB and VIS exposure; (ii) one-dimensional UVB and VIS spectra before and after response-correction, as well as after applying various processing, including absolute flux calibration, telluric removal, normalisation and barycentric correction; and (iii) co-added flux-calibrated and rectified spectra over the full optical range, for which all available XShootU exposures were combined. For many of the targets, the final signal-to-noise ratio per resolution element is above 200 in both the UVB and the VIS co-added spectra. The reduced data and advanced scientific data products will be made available to the community upon publication of this paper. [...]

Time-dependent meridional circulation and differential rotation in radiative zones are central open issues in stellar evolution theory. We streamline this challenging problem using the 'downward control principle' of atmospheric science, under a geostrophic f-plane approximation. New insights emerge from this simplified approach, using pressure as the vertical coordinate. We recover the known stellar physics result that the steady-state meridional circulation decays on the length scale (N/f sqrt(Pr))H, assuming molecular viscosity is the dominant drag mechanism. Prior to steady-state, the meridional circulation and the zonal wind (differential rotation) spread together via radiative diffusion, under thermal wind balance. The corresponding (4th-order) hyperdiffusion process is reasonably well approximated by regular (2nd-order) diffusion on scales of order a pressure scale-height. We derive an inhomogeneous diffusion (equiv. advection-diffusion) equation for the zonal flow which admits closed-form time-dependent solutions in a finite depth domain, allowing rapid prototyping of the meridional circulation pattern. In the weak drag limit, we find that the time to rotational steady-state can be longer than the Eddington-Sweet time and be instead determined by the longer drag time. Unless strong enough drag operates, the internal rotation of main-sequence stars may thus never reach steady-state. Streamlined meridional circulation solutions provide leading-order internal rotation profiles for studying the role of fluid/MHD instabilities (or waves) in redistributing angular momentum in the radiative zones of stars. Despite clear geometrical limitations and simplifying assumptions, one might expect our thin-layer geostrophic approach to offer qualitatively useful results to understand deep meridional circulation in stars.

As well known, magnetic fields in space are distributed very inhomogeneously. Some-times field distributions have forms of filaments with high magnetic field values. As many ob-servations show, such a filamentation takes place in convective cells in the Sun and other astro-physical objects. This effect is associated with the frozenness of the magnetic field into a medium with high conductivity that leads to compression of magnetic field lines and forming magnetic filaments. We show analytically, based on the general analysis, that the magnetic field intensifies in the regions of downward flows in both two-dimensional and three-dimensional convective cells. These regions of the hyperbolic type for magnetic fields play a role of a specific attractor. This analysis was confirmed by numerical simulations for 2D convective cells of the roll-type. Without dissipation the magnetic field grows exponentially in time and does not depend on the aspect ratio between horizontal and vertical scale of the cell. An increase due to compression in the magnetic field in the high conductive plasma is saturated due to the natural limitation associated with dissipative effects when the maximum magnitude of the magnetic field is of the order of the root of the magnetic Reynolds number Rem. For the solar convective zone the mean kinetic energy density exceeds mean magnetic energy density at least for two orders of magnitude that allows one to use the kinematic approximation for the MHD induction equation. In this paper based on the stability analysis we explain why downward flows influence magnetic filaments from making them more flat with orientation along interfaces between convective cells.

Dwarf irregulars (dIrrs) are among the most common type of galaxy in the Universe. They typically have gas-rich, low surface-brightness, metal-poor, and relatively-thick disks. Here we summarize the current state of our knowledge of the interstellar medium (ISM), including atomic, molecular and ionized gas, along with their dust properties and metals. We also discuss star formation feedback, gas accretion, and mergers with other dwarfs that connect the ISM to the circumgalactic and intergalactic media. We highlight one of the most persistent mysteries: the nature of pervasive gas that is yet undetected as either molecular or cold hydrogen, the ``dark gas''. Here are a few highlights: 1. Significant quantities of HI are in far-outer gas disks. 2. Cold HI in dIrrs would be molecular in the Milky Way, making the chemical properties of star-forming clouds significantly different. 3. Stellar feedback has a much larger impact in dIrrs than in spiral galaxies. 4. The escape fraction of ionizing photons is significant, making dIrrs a plausible source for reionization in the early Universe. 5. Observations suggest a significantly higher abundance of hydrogen (H$_2$ or cold HI) associated with CO in star-forming regions than that traced by the CO alone.

The flat-spectrum radio quasar 3C 345 has been showing gamma-ray activity since the mid-2000s along with activity across the electromagnetic spectrum. A gamma-ray burst in 2009 was successfully linked to relativistic outflow in 43 GHz VLBI observations and has since been analyzed also using single dish measurements. A multi-wavelength follow- up VLBI observation to the 2009 flare in conjunction with 43 GHz catalogue data from the VLBA-BU-BLAZAR and BEAM-ME programs are analyzed in this study in the context of the long-term evolution of the source. We aim to probe the innermost few milli-arcseconds of the ultracompact 3C 345 jet. In the process, we analyze the long-term kinematics of the inner jet and discuss the magnetic field morphology at different scales, as well as the origin of the gamma-ray emission. New observations at 23, 43, and 86 GHz took place on ten epochs between 2017 and 2019. We calibrate the 30 data sets using the rPicard pipeline, image them in Difmap and carry out polarization calibration using the GPCAL pipeline. We complement our VLBI data with ancillay VLBI maps at multiple frequencies, as well as 43 GHz observations carried out in the framework of the BEAM-ME and VLBA-BU-BLAZAR monitoring programs. We find multiple distinct component paths in the inner jet, forming a helical geometry. The helix is anchored at a stationary feature some 0.16 mas from the 43 GHz VLBI core and has a timescale of about 8 years. The characteristic bends in the jet morphology are caused by variations in the component ejection angle. We confirm the result of previous studies that the gamma-ray emission is produced (or caused) by relativistic outflow and violent interactions within the jet.

Population III (Pop III) stars formed out of metal free gas in minihalos at $z>20$. While their ignition ended the Dark Ages and begin enrichment of the IGM, their mass distribution remains unconstrained. To date, no confirmed Pop III star has been observed and their direct detection is beyond the reach of the James Webb Space Telescope (JWST) without gravitational lensing. However, a subset of massive Pop III stars end their lives in pair instability supernova (PISN). With typical energies of $\sim10^{53}$~erg, PISN light curve peaks are bright enough to be detectable by JWST and the Roman Space Telescope. The fundamental question of this work is whether or not observed PISN can be used as a diagnostic of the Pop III IMF. In this work, we use a model of the formation of the first stars to determine the dependence of PISN rates at $z~>~5$ for a range of Pop III power law IMFs ($\alpha~=~0.2~-~2.35$) and, critically, the method by which the IMF is populated. At $z~>~15$, we predict typical rates of $10^{-2}~-~10^2$ per deg$^{2}$ per year which will produce $10^{-3}~-~0.1$/year in a single NIRCam pointing and $0.003~-~30$/year in a single Roman pointing with $0.1~-~1000$ per year detected in the HLTDS. Our work highlights that theoretical modeling of PISN rates is required if upcoming PISN studies with JWST and Roman are going to constrain the Pop III IMF.

Context.Unusually, there are still certain characteristics of the changing-look (CL) active galactic nuclei (AGNs) that remain undetected.Consequently,the trigger mechanism behind the CL phenomenon observed in partial AGNs remains unknown.Aims.We explore the light curve and spectral energy distribution (SED) of the CL blazar OQ 334 as obtained by Fermi-LAT. Methods. By examining the variability of the equivalent width (EW), we categorise the Fermi-LAT light curves of OQ 334 during the epoch of MJD 54628-58677 into seven distinct epochs, including the flat spectrum radio quasar (FSRQ) state, the transition state, and the BL Lac state. We obtained both a Fermi-LAT SED and a multi-wavelength SED for each of these distinct epochs. Results. The source exhibits a transformation from a quiescent state to a highly active state, as evidenced by the variability of the EW. The multi-wavelength SEDs display a prominent external Compton characteristic, even though the Fermi-LAT SED reveals both a FSRQ and a BL Lac state across the seven different epochs. To gain further insights, we employed a leptonic model that takes into account the soft photon fields originating from both synchrotron radiation and the external environment. By simulating the multi-wavelength SEDs for each epoch, we uncover the following results. Firstly, the energy density of the external photon fields evolves in an oscillatory manner over the seven different epochs. Also, the energy density of the external photon fields in the BL Lac state is lower than that in the FSRQ state.Conclusions. These findings suggest that the CL blazar represents a unique phase in the blazar sequence.

In this article I summarize the current state of knowledge about the composition of Titan's atmosphere, and our current understanding of the suggested chemistry that leads to that observed composition. I begin with our present knowledge of the atmospheric composition, garnered from a variety of measurements including Cassini-Huygens, the Atacama Large Millimeter/submillimeter Array (ALMA), and other ground and space-based telescopes. This review focuses on the typical vertical profiles of gases at low latitudes, rather than global and temporal variations. The main body of the paper presents a chemical description of how complex molecules are believed to arise from simpler species, considering all known 'stable' molecules - those that have been uniquely identified in the neutral atmosphere. The last section of the paper is devoted to the gaps in our present knowledge of Titan's chemical composition and how further work may fill those gaps.

The Cosmological Advanced Survey Telescope for Optical and UV Research (CASTOR) is a planned flagship space telescope, covering the blue-optical and UV part of the spectrum. Here we introduce the CASTOR image simulator, a Python GalSim package-based script capable of generating mock CASTOR images from an input catalogue. We generate example images from the CASTOR Wide, Deep, and Ultra-Deep surveys using simulated light-cones from the Santa Cruz Semi-Analytic Model. We make predictions for the performance of these surveys by comparing galaxies that are extracted from each image using Source Extractor to the input catalogue. We find that the Wide, Deep, and Ultra-Deep surveys will be complete to ~27, 29 and 30 mag, respectively, in the UV, u, and g filters, with the UV-split and u-split filters reaching a shallower depth. With a large area of ~2200 deg$^2$, the Wide survey will detect hundreds of millions of galaxies out to z~4, mostly with $M_\ast \gtrsim 10^9 M_\odot$. The Ultra-Deep survey will probe to z~5, detecting a large fraction of $M_\ast \simeq 10^8 M_\odot$ galaxies. These powerful samples will enable precision measurements of the distribution of star formation in the cosmic web, connecting the growth of stellar mass to the assembly of dark matter halos over two thirds of the history of the Universe, and other core goals of CASTOR's legacy surveys. These image simulations and the tools developed to generate them will be a vital planning tool to estimate CASTOR's performance and iterate the telescope and survey designs prior to launch.

The XMM-Newton Line Emission Analysis Program (X-LEAP) is designed to study diffuse X-ray emissions from the Milky Way (MW) hot gas, as well as emissions from the foreground solar wind charge exchange (SWCX). This paper reports an all-sky survey of spectral feature intensities corresponding to the O VII, O VIII, and iron L-shell (Fe-L) emissions. These intensities are derived from 5418 selected XMM-Newton observations with long exposure times and minimal contamination from point or extended sources. For 90% of the measured intensities, the values are within $\approx$ 2-18 photons cm$^{-2}$ s$^{-1}$ sr$^{-1}$ (line unit; L.U.), $\approx$ 0-8 L.U., and $\approx$ 0-9 L.U., respectively. We report long-term variations in O VII and O VIII intensities over 22 years, closely correlating with the solar cycle and attributed to SWCX emissions. These variations contribute $\sim30\%$ and $\sim20\%$ to the observed intensities on average and peak at $\approx$ 4 L.U. and $\approx$ 1 L.U. during solar maxima. We also find evidence of short-term and spatial variations in SWCX, indicating the need for a more refined SWCX model in future studies. In addition, we present SWCX- and absorption-corrected all-sky maps for a better view of the MW hot gas emission. These maps show a gradual decrease in oxygen intensity moving away from the Galactic center and a concentration of Fe-L intensity in the Galactic bubbles and disk.

We propose a novel deep learning framework, named SYMHnet, which employs a graph neural network and a bidirectional long short-term memory network to cooperatively learn patterns from solar wind and interplanetary magnetic field parameters for short-term forecasts of the SYM-H index based on 1-minute and 5-minute resolution data. SYMHnet takes, as input, the time series of the parameters' values provided by NASA's Space Science Data Coordinated Archive and predicts, as output, the SYM-H index value at time point t + w hours for a given time point t where w is 1 or 2. By incorporating Bayesian inference into the learning framework, SYMHnet can quantify both aleatoric (data) uncertainty and epistemic (model) uncertainty when predicting future SYM-H indices. Experimental results show that SYMHnet works well at quiet time and storm time, for both 1-minute and 5-minute resolution data. The results also show that SYMHnet generally performs better than related machine learning methods. For example, SYMHnet achieves a forecast skill score (FSS) of 0.343 compared to the FSS of 0.074 of a recent gradient boosting machine (GBM) method when predicting SYM-H indices (1 hour in advance) in a large storm (SYM-H = -393 nT) using 5-minute resolution data. When predicting the SYM-H indices (2 hours in advance) in the large storm, SYMHnet achieves an FSS of 0.553 compared to the FSS of 0.087 of the GBM method. In addition, SYMHnet can provide results for both data and model uncertainty quantification, whereas the related methods cannot.

We report the conditional occurrences between three planetary types: super-Earths (m sin i $<$ 10 M$_\oplus$, P $<$ 100 days), warm Jupiters (m sin i $>$ 95 $M_\oplus$, 10 $<$ P $<$ 100 days), and cold Jupiters (m sin i $>$ 95 M$_\oplus$, P $>$ 400 days) for sun-like stars. We find that while the occurrence of cold Jupiters in systems with super-Earths is $22.2\substack{+8.3\\-5.4}$$\%$, compared to $10$$\%$ for the absolute occurrence rate of cold Jupiters, the occurrence of super-Earths in systems with cold Jupiters is $66.0\substack{+18.0\\-16.0}$$\%$, compared to $30$$\%$ for the absolute occurrence rate of super-Earths for sun-like stars. We find that the enhancement of super-Earths in systems with cold Jupiters is evident for sun-like stars, in agreement with several previous studies. We also conduct occurrence studies between warm Jupiters and super-Earths, and between warm Jupiters and cold Jupiters, to consolidate our methods. We conduct an independent observational test to study the effects of cold Jupiters against the inner multiplicity using the well-established giant planet host star metallicity correlation for all transiting planets found to date. The conditional occurrences we find here can be used to constrain the validity of various planetary formation models. The extremely interesting correlations between the super-Earths, cold Jupiters, and warm Jupiters can also be used to understand the formation histories of these planetary types.

The afterglow of a gamma-ray burst (GRB) has been widely argued to arise from the interaction of a relativistic outflow with its ambient medium. During such an interaction, a pair of shocks are generated: a forward shock that propagates into the medium, and a reverse shock that propagates into the outflow. Extensive studies have been conducted on the emission from the forward shock viewed off-axis. Furthermore, the observation of a reverse shock in an on-axis short GRB suggests that the reverse shock can produce an electromagnetic counterpart to a gravitational wave-detected merger. In this paper, we investigate the contribution of the reverse shock to the afterglow from a top-hat jet viewed off-axis, and apply our model to some short GRBs previously modeled by an off-axis emission. We employ the Markov Chain Monte Carlo (MCMC) method to get the model parameters (i.e., the jet's half-opening angle $\theta_j$, the viewing angle $\theta_\text{obs}$, the initial Lorentz factor $\Gamma_0$, and the isotropic energy $E_\mathrm{iso}$). Our model successfully reproduces off-axis afterglow emission without a structured jet. In addition, our calculations suggest that the reverse shock may produce a prominent feature in an early afterglow, which can be potentially observed in an orphan optical afterglow.

The interstellar extinction law is important for interpreting observations and inferring the properties of interstellar dust grains. Based on the 993 prism/CLEAR spectra from the James Webb Space Telescope (JWST), we investigate $0.6-5.3\,{\rm \mu m}$ interstellar dust extinction law. We propose a pair method to obtain the reddening curves based only on JWST observed spectra. Most of the high extinction sources are toward the young star cluster Westerlund 2. The infrared (IR) $1.0-5.3\,{\rm \mu m}$ reddening curves agree with the power-law $A_\lambda \propto \lambda^{-\alpha}$ well. We determine an average value of $\alpha=1.98\pm0.15$, which is consistent with the average value of the Galaxy. We find that $\alpha$ may be variable and independent of $R_{\rm V}$. With the derived $\alpha$, we convert the reddening curves into the extinction curves and establish the non-parameterized $\alpha$-dependent extinction curves in the wavelength range of $0.6-5.3\,{\rm \mu m}$. At $\lambda<1\,{\rm \mu m}$, the derived extinction law is not well described by the parameterized power-law type curve. Our non-parameterized $\alpha$-dependent extinction curves are suitable for the extinction correction of JWST-based photometry and spectra measurements at $0.6-5.3\,{\rm \mu m}$. We also provide the extinction coefficients for the JWST NIRCam bandpasses with different $\alpha$.

The cosmological 21cm signal offers a potential probe of the early Universe and the first ionizing sources. Current experiments probe the spatially-dependent variance (Gaussianity) of the signal through the power spectrum (PS). The signal however is expected to be highly non-Gaussian due to the complex topology of reionization and X-ray heating. We investigate the non-Gaussianities of X-ray heating and reionization, by calculating the skew spectrum (SS) of the 21cm signal using MERAXES, which couples a semi-analytic galaxy population with semi-numerical reionization simulations. The SS is the cross-spectrum of the quadratic temperature brightness field with itself. We generate a set of seven simulations from $z = 30$ to $z = 5$, varying the halo mass threshold for hosting star-formation, the X-ray luminosity per star-formation rate, and the minimum X-ray energy escaping host galaxies. We find the SS is predominantly negative as a function of redshift, transitioning to positive towards the start of reionization, and peaking during the midpoint of reionization. We do not see a negative dip in the SS during reionization, likely due to the specifics of modelling ionization sources. We normalise the SS by the PS during reionization isolating the non-Gaussianities. We find a trough ($k\sim\,0.1\,\textrm{Mpc}^{-1}$) and peak ($k\sim\,0.4-1\,\textrm{Mpc}^{-1}$) in the normalised SS during the mid to late periods of reionization. These correlate to the ionization topology, and neutral islands in the IGM. We calculate the cosmic variance of the normalised SS, and find these features are detectable in the absence of foregrounds with the SKA_LOW.

We propose a new method for estimating the distances of molecular clouds traced by CO line emission. Stars from 2MASS and Gaia EDR3 are selected as on-cloud stars when they are projected on a cloud. The background on-cloud stars have redder colors on average than the foreground stars. Instead of searching for stars projected away from the cloud, we employed the TRILEGA galaxy model to mimic the stellar population without cloud extinction along the sightline toward the cloud. Our method does not require an exact boundary of a cloud. The boundaries are highly variable and depend on the sensitivity of the molecular line data. For each cloud, we compared the distributions of on-cloud stars to the TRILEGAL stellar populations in the diagram of $J-K_s$ color versus distance. The intrinsic $J-K_s$ colors of main-sequence and evolved stars from TRILEGAL were considered separately, and they were used as the baseline for subtracting the observed $J-K_s$ colors. The baseline-corrected $J-K_s$ color was deployed with the Bayesian analysis and Markov chain Monte Carlo sampling to determine the distance at which the $J-K_s$ color jump is largest. This method was successfully applied to measure the distances of 27 molecular clouds, which were selected from previously published cloud samples. By replacing TRILEGAL with the GALAXIA galaxy model, we were able to measure the distances for 21 of the 27 clouds. The distances of the 21 clouds based on the GALAXIA model agree well with those based on the TRILEGAL model. The distances of the 27 clouds estimated by this method are consistent with previous estimates. We will apply this new method to a larger region of the gaseous galactic plane, in particular, for the inner galactic region, where a region free of CO emission is hard to separate from the crowded field of clouds.

The variety of long duration gamma-ray burst (LGRB) light curves (LCs) encode a wealth of information on how LGRB engines release energy following the collapse of the progenitor star. Attempts to characterise GRB LCs focused on a number of properties, such as the minimum variability timescale, power density spectra (both ensemble average and individual), or with different definitions of variability. In parallel, a characterisation as a stochastic process was pursued by studying the distributions of waiting times, peak flux, fluence of individual peaks within GRB time profiles. Yet, the question remains as to whether the diversity of profiles can be described in terms of a common stochastic process. Here we address this issue by studying for the first time the distribution of the number of peaks in a GRB profile. We used four different GRB catalogues: CGRO/BATSE, Swift/BAT, BeppoSAX/GRBM, and Insight-HXMT. The statistically significant peaks were identified by means of well tested algorithm MEPSA and further selected by applying a set of thresholds on signal-to-noise ratio. We then extracted the corresponding distributions of number of peaks per GRB. Among the different models considered (power-law, simple or stretched exponential) only a mixture of two exponentials models all the observed distributions, suggesting the existence of two distinct behaviours: (i) an average number of 2.1+-0.1 peaks per GRB ("peak poor") and accounting for about 80% of the observed population of GRBs; (ii) an average number of 8.3+-1.0 peaks per GRB ("peak rich") and accounting for the remaining 20% of the observed population. We associate the class of peak-rich GRBs with the presence of sub-second variability, which seems to be absent among peak-poor GRBs. The two classes could result from two different regimes through which GRB engines release energy or through which energy is dissipated into gamma-rays.

Starlight polarimetry, when combined with accurate distance measurements, allows for exploration of the three-dimensional structure of local magnetic fields in great detail. We present optical polarimetric observations of stars in and close to the Southern Coalsack, taken from the Interstellar Polarization Survey (IPS). Located in five fields of view approximately $0.3^{o}$ by $0.3^{o}$ in size, these data represent the highest density of optical polarimetric observations in the Southern Coalsack to date. Using these data, combined with accurate distances and extinctions based on Gaia data, we are able to characterize the magnetic field of the Coalsack and disentangle contributions to the polarization caused by the Southern Coalsack and a background structure. For the Southern Coalsack, we find an average magnetic field orientation of $\theta\sim 75^{o}$ with respect to the Galactic north pole and an average plane-of-sky magnetic field strength of approximately $B_{POS}=10$ $\mu G$, using the Davis-Chandrasekhar-Fermi (DCF) method. These values are in agreement with some earlier estimates of the Coalsack's magnetic field. In order to study the distant structure, we introduce a simple method to separate and isolate the polarization of distant stars from foreground contribution. For the distant structure, which we estimate to be located at a distance of approximately 1.3-1.5 kpc, we find an average magnetic field orientation of $\theta\sim100^{o}$ and we estimate a field strength of $B_{POS}\sim10 \ \mu G$, although this will remain highly uncertain until the precise nature of the distant structure can be uncovered.

Results of the long-term monitoring observations of the 6.7 GHz Class II methanol masers associated with the four high-mass star-forming regions by Hitachi 32-m radio telescope are presented. We detected periodic flux variability in G06.795-0.257, G10.472+0.027, G12.209-0.102, and G13.657-0.599 with the periods of 968, 1624, 1272, and 1266 d, respectively, although the detected period is tentative due to the short monitoring term relative to the estimated period. The facts that the flux variation patterns show the symmetric sine curves and that the luminosities of the central protostar and periods of maser flux variation are consistent with the expected period-luminosity (PL) relation suggest that the mechanisms of maser flux variability of G10.472+0.027 and G12.209-0.102 can be explained by protostellar pulsation instability. From the PL relation, central stars of these two sources are expected to be very high-mass protostars with a mass of 40 M and to have the mass accretion rate of 2*10-2 M yr-1. On the other hand, G06.795-0.257 and G13.657-0.599 have the intermittent variation patterns and have luminosities that are an order of magnitude smaller than those expected from PL relation, suggesting variation mechanisms of these sources originated from binary system. Since almost all the maser features vary with the same period regardless of its geometry, periodic accretion model may be appropriate mechanisms for flux variability in G06.795-0.257 and G13.657-0.599.

The search for exomoons in time-domain photometric data has to-date generally consisted of fitting transit models that are comprised of a planet hosting a single moon. This simple model has its advantages, but it may not be particularly representative, as most of the major moons in our Solar System are found in multi-moon satellite systems. It is critical that we investigate, then, the impact of applying a single-moon model to systems containing multiple moons, as there is the possibility that utilizing an inaccurate or incomplete model could lead to erroneous conclusions about the system. To that end, in this work we produce a variety of realistic multi-moon light curves, perform standard single-moon model selection, and analyze the impacts that this model choice may have on the search for exomoons. We find that the number of moons in a system fit with a single-moon model generally has little impact on whether we find evidence for a moon in that system, and other system attributes are individually not especially predictive. However, the model parameter solutions for the moon frequently do not match any real moon in the system, instead painting a picture of a ``phantom'' moon. We find no evidence that multi-moon systems yield corresponding multi-modal posteriors. We also find a systematic tendency to overestimate planetary impact parameter and eccentricity, to derive unphysical moon densities, and to infer potentially unphysical limb darkening coefficients. These results will be important to keep in mind in future exomoon search programs.

TOI-1518b, a hot Jupiter around a late A-type star, is one of the few planetary systems that transit the edge of the stellar surface (the impact parameter $b\sim0.9 $) among hot Jupiters around hot stars (Cabot et al. 2021). The high rotation speed of the host star ($\sim85$ km s$^{-1}$) and the nearly polar orbit of the planet ($\sim 120$ deg) may cause a nodal precession. In this study, we report the nodal precession undergone by TOI-1518\,b. This system is the fourth planetary system in which nodal precession is detected. We investigate the time change in $b$ from the photometric data of TOI-1518 acquired in 2019 and 2022 with TESS and from the spectral transit data of TOI-1518b obtained in 2020 with two high-dispersion spectrographs; CARMENES and EXPRES. We find that the value of $b$ is decreasing with $db/dt=-0.0116\pm0.0036$\,year$^{-1}$, indicating that the transit trajectory is moving toward the center of the stellar surface. We also estimate the minimum value of the quadrupole mass moment of TOI-1518 $J_{2,\mathrm{min}}=4.41\times 10^{-5}$ and the logarithm of the Love number of TOI-1518 $\log{k_2}= -2.17\pm 0.33$ from the nodal precession.

Electrospheres are environments with the same origin as pulsars; a highly magnetized rotating neutron star. In pulsars, a cascade of electron-positron pair creation enriches the plasma. The plasma surrounding an electrosphere consists only of particles that have escaped from the neutron star's surface. Electrospheres with a magnetic axis aligned with the rotation axis have been well described for decades. Models of electrospheres with an oblique magnetic axis relative to the rotation axis have resisted most theoretical investigations. Some electrospheres and pulsars have been simulated using particle-in-cell codes, but the numerical constraints did not allow the use of realistic neutron star parameters. We aimed to develop a numerical simulation code optimized for understanding the physics of electrospheres and pulsars, with realistic neutron star parameters. As a first step, presented in this paper, we focused on the simulation of oblique electrospheres with realistic physical parameters. A specific code was developed for the computation of stationary solutions. The resolution of Maxwell's equations was based on spectral methods. Particle motions included their finite inertia. No hypothesis was made in relation to the force-free behavior of the electrospheric plasma. The numerical code is called Pulsar ARoMa (pulsar asymmetric rotating magnetosphere). Various numerical simulations were conducted using realistic neutron star parameters. We find that oblique electrospheres possess the same global structure as aligned force-free electrospheres, with two domes of electrons and a torus of positively charged particles. The domes are not centered on the magnetic axis; nor are they symmetric. Yet, the solutions do not exhibit a force-free behavior. The simulations performed with the Pulsar ARoMa code require modest resources and little computing time. This code will be upgraded for more ambitious investigations into pulsar physics.

The cosmic distance duality relation (CDDR), expressed as DL(z) = (1 + z)2DA(z), plays an important role in modern cosmology. In this paper, we propose a new method of testing CDDR using strongly lensed gravitational wave (SLGW) signals. Under the geometric optics approximation, we calculate the gravitational lens effects of two lens models, the point mass and singular isothermal sphere. We use functions of {\eta}1(z) = 1 + {\eta}0z and {\eta}2(z) = 1 + {\eta}0z=(1 + z) to parameterize the deviation of CDDR. By reparameterizing the SLGW waveform with CDDR and the distance-redshift relation, we include the deviation parameters {\eta}0 of CDDR as waveform parameters. We evaluate the ability of this method by calculating the parameter estimation of simulated SLGW signals from massive binary black holes. We apply the Fisher information matrix and Markov Chain Monte Carlo methods to calculate parameter estimation. We find that with only one SLGW signal, the measurement precision of {\eta}0 can reach a considerable level of 0.5-1.3% for {\eta}1(z) and 1.1-2.6% for {\eta}2(z), depending on the lens model and parameters.

The Cygnus region, which contains massive molecular and atomic clouds and young stars, is a promising Galactic neutrino source candidate. Cosmic rays transport in the region can produce neutrinos and $\gamma$-rays. Recently, the Large High Altitude Air Shower Observatory (LHAASO) detected an ultrahigh-energy $\gamma$-ray bubble (Cygnus Bubble) in this region. Using publicly available track events detected by the IceCube Neutrino Observatory in 7 years of full detector operation, we conduct searches for correlated neutrino signals from the Cygnus Bubble with neutrino emission templates based on LHAASO $\gamma$-ray observations. No significant signals were found for any employed templates. With the 7 TeV $\gamma$-ray flux template, we set a flux upper limit of 90% confidence level (C.L.) for the neutrino emission from the Cygnus Bubble to be $5.7\times10^{-13}\, \mathrm{TeV}^{-1}\mathrm{cm}^{-2}\mathrm{s}^{-1}$ at 5 TeV.

We report on the results of a search for radio pulsars in five supernova remnants (SNRs) with FAST. The observations were made using the 19-beam receiver in the Snapshot mode. The integration time for each pointing is 10 min. We discovered a new pulsar PSR J1845$-$0306 which has a spin period of 983.6 ms and a dispersion measure of 444.6$\pm$2.0 cm$^{-3}$ pc in observations of SNR G29.6+0.1. To judge the association between the pulsar and the SNR, further verification is needed. We also re-detected some known pulsars in the data from SNRs G29.6+0.1 and G29.7$-$0.3. No pulsars were detected in observations of other three SNRs.

IGR J0607.4+2205 is a transient Be X-ray binary discovered two decades ago. IGR J0607.4+2205 underwent an outburst in 2023 during which it was observed twice with \textit{NuSTAR}. The main goal of this work is to model the broadband X-ray spectrum of IGR J0607.4+2205 during the outburst and to study the variations of the spectral and timing features at different intensities. We extracted the light curve and spectrum of the source from the two \textit{NuSTAR} observations carried out during the recent outburst in the energy range of 3$-$78 keV. We used the epoch folding technique to find pulsation from the source and to study the changes in emission characteristics from the source with energy across an order of magnitude variation in source luminosity. IGR J0607.4+2205 shows pulsations with a period of $\sim$347.6 s during both the observations, with a pulse fraction of $\geq$50\%. The broadband spectrum of the source was modelled using a power-law continuum with a high-energy cutoff. During the first observation, a cyclotron absorption line at $\sim$51 keV was also present in the source with an optical depth of $\sim$1.3. However, no cyclotron line feature was detected in the second observation when the source was an order of magnitude fainter. Additionally, soft excess was detected in the second observation, which was modelled with a black body component emerging from close to the neutron star (NS). We report the first ever detection of a cyclotron line in the broadband spectrum of IGR J0607.4+2205 centred at 51$\pm$1 keV. The magnetic field strength of the NS is estimated to be $\sim$4$\times$$10^{12}$ G from the centroid energy of the absorption line. A significant change is observed in the pulse profile with luminosity during the decay of the outburst, indicating an associated change in the beaming pattern.

M-type stars are crucial for stellar activity studies since they cover two types of magnetic dynamos and particularly intriguing for habitability studies due to their abundance and long lifespans during the main-sequence stage. In this paper, we used the LAMOST DR9 catalog and the GALEX UV archive data to investigate the chromospheric and UV activities of M-type stars. All the chromospheric and UV activity indices clearly show the saturated and unsaturated regimes and the well-known activity-rotation relation, consistent with previous studies. Both the FUV and NUV activity indices exhibit a single-peaked distribution, while the {\rm H$\alpha$} and \rm {Ca \scriptsize{\uppercase\expandafter{\romannumeral2}} \normalsize H$\&$K} indices show a distinct double-peaked distribution. The gap between these peaks suggests a rapid transition from a saturated population to an unsaturated one. The smoothly varying distributions of different subtypes suggest a rotation-dependent dynamo for both early-type (partly convective) to late-type (fully convective) M stars. We identified a group of stars with high UV activity above the saturation regime (log$R^{\prime}_{\rm NUV} > -2.5$) but low chromospheric activity, and the underlying reason is unknown. By calculating the continuously habitable zone and the UV habitable zone for each star, we found about 70\% stars in the total sample and 40\% stars within 100 pc are located in the overlapping region of these two habitable zones, indicating a number of M stars are potentially habitable. Finally, we examined the possibility of UV activity studies of M stars using the China Space Station Telescope.

Context: Many barriers prevent dust to form planetesimals via coagulation in protoplanetary discs, such as bouncing, collisional fragmentation or aeolian erosion. Modelling dust and the different phenomena that can alter its evolution is therefore needed. Multiple solutions have been proposed, but still need to be confirmed. Aims: In this paper, we explore the role aeolian erosion plays in the evolution of dust. Methods: We use a monodisperse model to account for dust growth and fragmentation, implemented in a 1D model to compute the evolution of single grains and a 3D SPH code to compute the global evolution of dust and gas. We test the erosion model in our code and ensured it matches previous results. Results: With a model of disc reproducing observations, we show with both 1D and 3D studies that erosion is not significant during the evolution of dust when we take fragmentation into consideration. With a low-viscosity disc, fragmentation is less of a problem, but grain growth is also less important, preventing the formation of large objects anyway. In dust traps, close to the star, erosion is also not impactful, even when fragmentation is turned off. Conclusions: We show in this paper that aeolian erosion is negligible when radial drift, fragmentation and dust traps are taken into account and does not alter the dust evolution in the disc. However, it can have an impact on later stages, i.e. when the streaming instability forms large clumps close to the star, or when planetesimals are captured.

Galaxy clusters that host radio halos indicate the presence of population(s) of non-thermal electrons. These electrons can scatter low-energy photons of the Cosmic Microwave Background, resulting in the non-thermal Sunyaev-Zeldovich (ntSZ) effect. We measure the average ntSZ signal from 62 radio-halo hosting clusters using the $Planck$ multi-frequency all-sky maps. We find no direct evidence of the ntSZ signal in the $Planck$ data. Combining the upper limits on the non-thermal electron density with the average measured synchrotron power collected from the literature, we place lower limits on the average magnetic field strength in our sample. The lower limit on the volume-averaged magnetic field is $0.1-0.01\,\mu$G, depending on the assumed power-law distribution of electron energies. We further explore the potential improvement of these constraints from the upcoming Simons Observatory and Fred Young Submillimeter Telescope (FYST) of the CCAT-prime collaboration. We find that combining these two experiments, the constraints will improve by a factor of $3-4$, which can be sufficient to rule out some power-law models.

Despite the thousands of planets in orbit around stars known to date, the mechanisms of planetary formation, migration, and atmospheric loss remain unresolved. In this work, we confirm the planetary nature of a young Saturn-size planet transiting a solar-type star every 8.03 d, TOI-1135\,b. The age of the parent star is estimated to be in the interval of 125--1000 Myr based on various activity and age indicators, including its stellar rotation period of 5.13\,$\pm$\,0.27 d and the intensity of photospheric lithium. We obtained follow-up photometry and spectroscopy, including precise radial velocity measurements using the CARMENES spectrograph, which together with the TESS data allowed us to fully characterise the parent star and its planet. As expected for its youth, the star is rather active and shows strong photometric and spectroscopic variability correlating with its rotation period. We modelled the stellar variability using Gaussian process regression. We measured the planetary radius at 9.02\,$\pm$\,0.23 R$_\oplus$ (0.81\,$\pm$\,0.02 R$_{\mathrm{Jup}}$) and determined a 3$\sigma$ upper limit of $<$\,51.4 M$_\oplus$ ($<$\,0.16 \,M$_{\rm{Jup}}$) on the planetary mass by adopting a circular orbit. Our results indicate that TOI-1135\,b is an inflated planet less massive than Saturn or Jupiter but with a similar radius, which could be in the process of losing its atmosphere by photoevaporation. This new young planet occupies a region of the mass-radius diagram where older planets are scarse, and it could be very helpful to understanding the lower frequency of planets with sizes between Neptune and Saturn.

Observations show that faster-rotating stars tend to have stronger magnetic activity and shorter magnetic cycles. The cyclical magnetic activity of the Sun and stars is believed to be driven by the dynamo process. The success of the Babcock-Leighton (BL) dynamo in understanding the solar cycle suggests an important role that starspots could play in stellar magnetic cycles. We aim at extending the BL mechanism to solar-mass stars with various rotation rates and explore the effects of emergence properties of starspots in latitudes and tilt angles on stellar magnetic cycles. We adopt a kinematic BL-type dynamo model operating in the bulk of the convection zone. The profiles of the large-scale flow fields are from the mean-field hydrodynamical model for various rotators. The BL source term in the model is constructed based on the rotation dependence of starspots emergence. That is, faster rotators have starspots at higher latitudes with larger tilt angles.Faster rotators have poloidal flux appearing closer to about $\pm55^\circ$ latitudes, where the toroidal field generation efficiency is the strongest because of the strongest latitudinal differential rotation there. It takes a shorter time for faster rotators to transport the surface poloidal field from their emergence latitude to the $\pm 55^\circ$ latitudes of efficient $\Omega$-effect thus shortening their magnetic cycles. The faster rotators operate in a more supercritical regime due to a stronger BL $\alpha$-effect relating to the tilt angles, which leads to stronger saturated magnetic fields and a coupling of the poloidal field between two hemispheres more difficult. Thus the magnetic field parity shifts from the hemispherically asymmetric mixed mode to quadrupole, and further to dipole when a star spins down. The emergence of starspots plays an essential role in the large-scale stellar dynamo.

We present an analysis of the combined NIRSpec and MIRI spectra of dusty galaxies between 1.5 - 28 $\mu$m restframe by implementing a differential extinction model, where the strength of extinction varies across the spectrum as different layers of the obscuring dust are probed. Our model is able to recover a 2D distribution of dust temperature and extinction which allows inference of the physical nature of the dust in these environments. We show that differential extinction is necessary to reproduce the spectra of 4 highly obscured Luminous Infrared Galaxies observed with NIRSpec IFU and MIRI MRS, where simple screen or uniformly mixed dust distributions fail to fit the data. We additionally compare the extinction of HII regions in these galaxies via hydrogen recombination lines, the extinction of molecular gas via the H$_2$ lines, Polycyclic Aromatic Hydrocarbons via the 12.7/11.3 PAH ratio and the stellar continuum. We find that the molecular gas is deeply buried with the HII regions in star-forming regions, with a similar extinction to the hottest dust components. However we find the cooler dust to be less obscured, at a similar extinction to the stellar continuum and PAHs. The nuclei show a complex dust distribution with VV114 NE, NGC 3256 S, IIZw96 SW showing a deeply buried continuum source relative to the molecular gas/HII regions. Additionally, NGC 3256 S, NGC 7469 and VV114 SW show an isolated hot dust component, indicative of AGN heating, where NGC 3256 S and NGC 7469 are previously known AGN.

A detailed interpretation of the detected emission lines of environments in which propyne (or methyl acetylene, CH$_3$CCH) is observed requires a knowledge of its collisional rate coefficients with the most abundant species in the interstellar medium, He or H$_2$. We present the first three-dimensional potential energy surface (3D-PES) for the CH$_3$CCH-He molecular complex, study the dynamics of the collision, and report the first set of rate coefficients for temperatures up to 100 K for the collisional excitation of the lowest 60 ortho rotational levels and 60 para rotational levels of CH$_3$CCH by He atoms. We computed the 3D-PES with the explicitly correlated coupled-cluster with single-, double-, and perturbative triple-excitation method, in conjunction with the augmented correlation-consistent triple zeta basis set (CCSD(T)-F12a/aug-cc-pVTZ). The 3D-PES was fitted to an analytical function. Scattering computations of pure rotational (de-)excitation of CH$_3$CCH by collision with He atoms were performed and the state-to-state cross sections were computed using the close coupling method for total energies up to 100 cm$^{-1}$ and with the coupled states approximation at higher energy for both ortho and para symmetries of CH$_3$CCH. The PES obtained is caracterized by a large anisotropy and a potential well depth of 51.04 cm$^{-1}$. By thermally averaging the collisional cross sections, we determined quenching rate coefficients for kinetic temperatures up to 100 K. A strong even $\Delta j$ propensity rule at almost all collision energies exists for CH$_3$CCH-He complex. To evaluate the impact of rate coefficients in the analysis of observations, we carried out non-LTE radiative transfer computations of the excitation temperatures and we demonstrate that LTE conditions are typically not fulfilled for the propyne molecule.

Rapid and accurate evaluation of the nonlinear matter power spectrum, $P(k)$, as a function of cosmological parameters and redshift is of fundamental importance in cosmology. Analytic approximations provide an interpretable solution, yet current approximations are neither fast nor accurate relative to black-box numerical emulators. We use symbolic regression to obtain simple analytic approximations to the nonlinear scale, $k_\sigma$, the effective spectral index, $n_{\rm eff}$, and the curvature, $C$, which are required for the halofit model. We then re-optimise the coefficients of halofit to fit a wide range of cosmologies and redshifts. We then again exploit symbolic regression to explore the space of analytic expressions to fit the residuals between $P(k)$ and the optimised predictions of halofit. All methods are validated against $N$-body simulations. Our symbolic expressions for $k_\sigma$, $n_{\rm eff}$ and $C$ have root mean squared fractional errors of 0.8%, 0.2% and 0.3%, respectively, for redshifts below 3 and a wide range of cosmologies. The re-optimised halofit parameters reduce the root mean squared fractional error from 3% to below 2% for wavenumbers $k=9\times10^{-3}-9 \, h{\rm Mpc^{-1}}$. We introduce syren-halofit (symbolic-regression-enhanced halofit), an extension to halofit containing a short symbolic correction which improves this error to 1%. Our method is 2350 and 3170 times faster than current halofit and hmcode implementations, respectively, and 2680 and 64 times faster than EuclidEmulator2 (which requires running class) and the BACCO emulator. We obtain comparable accuracy to EuclidEmulator2 and the BACCO emulator when tested on $N$-body simulations. Our work greatly increases the speed and accuracy of symbolic approximations to $P(k)$, making them significantly faster than their numerical counterparts without loss of accuracy.

Massive stars can during their evolution reach the phase of critical (or very rapid, near-critical) rotation when further increase in rotation rate is no longer kinematically allowed. The mass ejection and angular momentum outward transport from such rapidly rotating star's equatorial surface may lead to formation and supports further existence of a circumstellar outflowing (stellar decretion) disk. The efficient mechanism for the outward transport of the mass and angular momentum is provided by the anomalous viscosity. The outer supersonic regions of the disks can extend up to a significantly large distance from the parent star, the exact radial extension is however basically unknown, partly due to the uncertainties in radial variations of temperature and viscosity. We study in detail the behavior of hydrodynamic quantities, i.e., the evolution of density, radial and azimuthal velocity, and angular momentum loss rate in stellar decretion disks out to extremely distant regions. We investigate the dependence of these physical characteristics on the distribution of temperature and viscosity. We also study the magnetorotational instability, which we regard to be the source of anomalous viscosity in such outflowing disks and to some extent we provide the preliminary models of the two-dimensional radially-vertically correlated distribution of the disk density and temperature. We developed our own two-dimensional hydrodynamic and magnetohydrodynamic numerical code based on an explicit Eulerian finite difference scheme on staggered grid, including full Navier-Stokes viscosity. We use semianalytic approach to investigate the radial profile of magnetorotational instability, where on the base of the numerical time-dependent hydrodynamic model we analytically study the stability of outflowing disks submerged to the magnetic field of central star.

Upcoming Cosmic Microwave Background (CMB) experiments, aimed at measuring primordial CMB B-modes, require exquisite control of Galactic foreground contamination. Minimum-variance techniques, like the Needlet Internal Linear Combination (NILC), have proven effective in reconstructing the CMB polarization signal and mitigating foregrounds across diverse sky models without suffering from mismodelling errors. Still, residual contamination may bias the recovered CMB polarization at large angular scales when confronted with the most complex foreground scenarios. By adding constraints to NILC to deproject moments of the Galactic emission, the Constrained Moment ILC (cMILC) method has proven to enhance foreground subtraction, albeit with an associated increase in overall noise variance. Faced with this trade-off between foreground bias reduction and overall variance minimization, there is still no recipe on which moments to deproject and which are better suited for blind variance minimization. To address this, we introduce the optimized cMILC (ocMILC) pipeline, which performs full optimization of the required number and set of foreground moments to deproject, pivot parameter values, and deprojection coefficients across the sky and angular scales, depending on the actual sky complexity, available frequency coverage, and experiment sensitivity. The optimal number of deprojected moments, before paying significant noise penalty, is determined through a data diagnosis inspired by the Generalized NILC (GNILC) method. Validated on B-mode simulations of the PICO space mission concept with four challenging foreground models, ocMILC exhibits lower foreground contamination compared to NILC and cMILC at all angular scales, with limited noise penalty. This multi-layer optimization enables the ocMILC pipeline to achieve unbiased posteriors of the tensor-to-scalar ratio, regardless of foreground complexity.

Modern galaxy surveys demand extensive survey volumes and resolutions surpassing current dark matter-only simulations' capabilities. To address this, many methods employ effective bias models on the dark matter field to approximate object counts on a grid. However, realistic catalogs necessitate specific coordinates and velocities for a comprehensive understanding of the Universe. In this research, we explore sub-grid modeling to create accurate catalogs, beginning with coarse grid number counts at resolutions of approximately $5.5,h^{-1}$ Mpc per side. These resolutions strike a balance between modeling nonlinear damping of baryon acoustic oscillations and facilitating large-volume simulations. Augmented Lagrangian Perturbation Theory (ALPT) is utilized to model the dark matter field and motions, replicating the clustering of a halo catalog derived from a massive simulation at $z=1.1$. Our approach involves four key stages: Tracer Assignment: Allocating dark matter particles to tracers based on grid cell counts, generating additional particles to address discrepancies. Attractor Identification: Defining attractors based on particle cosmic web environments, acting as gravitational focal points. Tracer Collapse: Guiding tracers towards attractors, simulating structure collapse. Redshift Space Distortions: Introducing redshift space distortions to simulated catalogs using ALPT and a random dispersion term. Results demonstrate accurate reproduction of monopoles and quadrupoles up to wave numbers of approximately $k=0.6,h$ Mpc$^{-1}$. This method holds significant promise for galaxy surveys like DESI, EUCLID, and LSST, enhancing our understanding of the cosmos across scales.

Compact obscured nuclei (CONs) are an extremely obscured (N$_{H2}$ >10$^{25}$ cm$^{-2}$) class of galaxy nuclei thought to exist in 20-40 per cent of nearby (ultra-)luminous infrared galaxies. While they have been proposed to represent a key phase of the active galactic nucleus (AGN) feedback cycle, the nature of these CONs - what powers them, their dynamics, and their impact on the host galaxy - remains unknown. This work analyses the large-scale optical properties of the local CON, NGC4418 (z=0.00727). We present new, targeted integral field unit observations of the galaxy with the Multi-Unit Spectroscopic Explorer (MUSE). For the first time, we map the ionised and neutral gas components of the galaxy, along with their dynamical structure, to reveal several previously unknown features of the galaxy. We confirm the presence of a previously postulated blueshifted outflow along the minor axis of NGC4418. We find this outflow to be decelerating and, for the first time, show it to extend bilaterally from the nucleus. We report the discovery of two further outflow structures: a redshifted southern outflow connected to a tail of ionised gas surrounding the galaxy and a blueshifted bubble to the north. In addition to these features, we find the [OIII] emission reveals the presence of knots across the galaxy, which are consistent with regions of the galaxy that have been photoionised by an AGN. Based on the properties of these features, we conclude that the CON in NGC4418 is most likely powered by AGN activity.

While gravitational lens inversion holds great promise to reveal the structure of the light-deflecting mass distribution, both light and dark, the existence of various kinds of degeneracies implies that care must be taken when interpreting the resulting lens models. This article illustrates how thinking in terms of the projected potential helps to gain insight into these matters. Additionally it is shown explicitly how, when starting from a discretised version of the projected potential of one particular lens model, the technique of quadratic programming can be used to create a multitude of equivalent lens models that preserve all or a subset of lens properties. This method is applied to a number of scenarios, showing the lack of grasp on the mass outside the strong lensing region, revisiting mass redistribution in between images and applying this to a recent model of the SDSS J1004+4112 cluster, as well as illustrating the generalised mass sheet degeneracy and source-position transformation. In the case of J1004 we show that this mass redistribution did not succeed at completely eliminating a dark mass clump recovered by GRALE near one of the quasar images.

The Legacy Survey of Space and Time (LSST) will revolutionize Time Domain Astronomy by detecting millions of transients. In particular, it is expected to increment the number of type Ia supernovae (SNIa) of a factor of 100 compared to existing samples up to z~1.2. Such a high number of events will dramatically reduce statistical uncertainties in the analysis of SNIa properties and rates. However, the impact of all other sources of uncertainty on the measurement must still be evaluated. The comprehension and reduction of such uncertainties will be fundamental both for cosmology and stellar evolution studies, as measuring the SNIa rate can put constraints on the evolutionary scenarios of different SNIa progenitors. We use simulated data from the DESC Data Challenge 2 (DC2) and LSST Data Preview 0 (DP0) to measure the SNIa rate on a 15 deg2 region of the Wide-Fast-Deep area. We select a sample of SN candidates detected on difference images, associate them to the host galaxy, and retrieve their photometric redshifts (z-phot). Then, we test different light curves classification methods, with and without redshift priors. We discuss how the distribution in redshift measured for the SN candidates changes according to the selected host galaxy and redshift estimate. We measure the SNIa rate analyzing the impact of uncertainties due to z-phot, host galaxy association and classification on the distribution in redshift of the starting sample. We found a 17% average lost fraction of SNIa with respect to the simulated sample. As 10% of the bias is due to the uncertainty on the z-phot alone (which also affects classification when used as a prior), it results to be the major source of uncertainty. We discuss possible reduction of the errors in the measurement of the SNIa rate, including synergies with other surveys, which may help using the rate to discriminate different progenitor models.

Lya emitters (LAEs) are particularly useful objects for the study of the Epoch of Reionization. Lya profiles can be used to estimate the amount of ionizing photons that are able to escape the galaxies, and therefore to understand what objects contributed to reionization. However, Lya is a resonant line and its complex radiative transfer effects make the interpretation of the line challenging and require the use of appropriate radiative transfer methods for anything but the simplest gas distributions. With this work we aim to study the properties of simulated LAEs, and the robustness of these inferred properties under a change in the dust model. We also explore the Lyman Continuum (LyC) escape fraction of these galaxies and compare our results with observationally calibrated methods to infer this quantity from the Lya spectrum. We use the radiative transfer code RASCAS to perform synthetic observations of 13 flux-selected galaxies from the Obelisk simulation at redshift z = 6, towards the end of the Epoch of Reionization. Each galaxy was observed in Lya, ionizing and non-ionizing continuum from 48 different viewing angles. We show that the Lya profiles emitted from a galaxy present large variations with a change in viewing angle and that the relation between peak separation and Lya escape fraction is not as strong as previously found, as we find lines of sight with both low peak separation and low escape fraction, due to their dust content. We also show that the properties of the Lya line are reasonably robust under a change in dust model. Lastly, we compare the LyC escape fractions we derive from the simulation to three observationally calibrated methods of inferring this quantity. We determine that none of these relations reproduce the scatter that we find in our sample, and that high escape fraction lines of sight have both low peak separation and low dust extinction in the UV.

We present the results of a polarimetric study from our new high-sensitivity L-band (0.8--1.7 GHz) observation of Pictor A with the MeerKAT radio telescope. We confirm the presence of the radio jet extending from the nucleus to the western hotspot of this source. Additionally, we show the radio emission expected to be coincident with previously observed X-ray emission in the radio lobes, confirming that the emission mechanism is of inverse Compton origin, as suggested by a previous study. Our spectropolarimetric analysis using the RM-Synthesis technique reveals a relatively uniform mean RM distribution across the lobes of Pictor A, with most lines-of-sight exhibiting single-peaked Faraday spectra. However, a number of the lines-of-sight exhibit single peaked spectra with a wide base or multiple peaks, suggesting the presence of multiple Faraday components or a Faraday thick structure along Pictor A's lines-of-sight. We also confirm the asymmetry in RM variability and depolarization between the two lobes of this source which were reported in a previous study.

The evolution of dwarf galaxies is dramatically affected by gaseous and dusty outflows, as they can easily deprive their interstellar medium of the material needed for the formation of new stars, simultaneously enriching their surrounding circumgalactic medium (CGM). In this Letter, we present the first evidence for extended [CII] 158 $\mu$m line and dust continuum emission in local dwarf galaxies hosting star-formation-driven outflows. By stacking the [CII], far-infrared and near-UV (NUV) emission obtained from $Herschel$ and GALEX data, we derive the average radial profiles, and compare the spatial extension of gas, dust, and stellar activity in dwarf galaxies. We find that [CII] and dust emissions are comparable to each other, and more extended than the NUV continuum. The [CII] size is in agreement with that measured for $z>4$ star-forming galaxies, suggesting that similar mechanisms could be at the origin of the observed atomic carbon reservoir around local and high-$z$ sources. The cold dust follows the [CII] emission, going beyond the stellar continuum as opposed to what is typically observed in the early Universe where measurements can be affected by the poor sensitivity and faintness of dust emission in the CGM of high-$z$ galaxies. We attribute the extended [CII] and dust continuum emission to the presence of galactic outflows. As local dwarf galaxies are considered analogs of primordial sources, we could expect that comparable feedback processes can be at the origin of the observed [CII] halos at $z>4$, dominating over other possible formation mechanisms.

The calculation of internal atmospheric (longwave) fluxes is a key component of any model of exoplanet atmospheres that requires radiative-transfer (RT) calculations. For atmospheres containing a strong scattering component such as cloud particles, most 1D multiple-scattering RT methods typically involve numerically expensive matrix inversions. This computational bottleneck is exacerbated when multitudes of RT calculations are required, such as in general circulation models (GCMs) and retrieval methods. In an effort to increase the speed of RT calculations without sacrificing too much accuracy, we investigate the applicability of approximate longwave scattering methods developed for the Earth science community to hot Jupiter atmospheres. We test the absorption approximation (AA) and variational iteration method (VIM) applied to typical cloudy hot Jupiter scenarios, using 64 stream DISORT calculations as reference solutions. We find the four-stream VIM variant is a highly promising method to explore using for hot Jupiter GCM and retrieval modelling, showing excellent speed characteristics, with typical errors $\sim$1\% for outgoing fluxes and within $\sim$50\%, but with larger errors in the deep cloud layer test case, for vertical heating rates. Other methods explored in this study were found to typically produce similar error characteristics in vertical heating rates.

We present a new class of machine-learning emulators that accurately model the cosmic shear, galaxy-galaxy lensing, and galaxy clustering real space correlation functions in the context of Rubin Observatory year one simulated data. To illustrate its capabilities in forecasting models beyond the standard $\Lambda$CDM, we forecast how well LSST Year 1 data will be able to probe the consistency between geometry $\Omega^{\rm geo}_\mathrm{m}$ and growth $\Omega^{\rm growth}_\mathrm{m}$ dark matter densities in the so-called split $\Lambda$CDM parameterization. When trained with a few million samples, our emulator shows uniform accuracy across a wide range in an 18-dimensional parameter space. We provide a detailed comparison of three neural network designs, illustrating the importance of adopting state-of-the-art Transformer blocks. Our study also details their performance when computing Bayesian evidence for cosmic shear on three fiducial cosmologies. The transformers-based emulator is always accurate within PolyChord's precision. As an application, we use our emulator to study the degeneracies between dark energy models and growth geometry split parameterizations. We find that the growth-geometry split remains to be a meaningful test of the smooth dark energy assumption.

Precise Doppler studies of extrasolar planets require fine-grained control of observational cadence, i.e. the timing of and spacing between observations. We present a novel framework for scheduling a set of Doppler campaigns with different cadence requirements at the W. M. Keck Observatory (WMKO). For a set of observing programs and allocated nights on an instrument, our software optimizes the timing and ordering of ~1000 observations within a given observing semester. We achieve a near-optimal solution in real-time using a hierarchical Integer Linear Programming (ILP) framework. Our scheduling formulation optimizes over the roughly 10^3000 possible orderings. A top level optimization finds the most regular sequence of allocated nights by which to observe each host star in the request catalog based on a frequency specified in the request. A second optimization scheme minimizes the slews and downtime of the instrument. We have assessed our algorithms performance with simulated data and with the real suite of Doppler observations of the California Planet Search in 2023.

We use a recent Pantheon+SH0ES compilation of Type Ia Supernova distance measurements at low-redshift, i.e., $0.01 \leq z \leq 0.10$, in order to investigate the directional dependency of the deceleration parameter ($q_0$) in different patches ($60^{\circ}$ size) across the sky, as a probe of the statistical isotropy of the Universe. We adopt a cosmographic approach to compute the cosmological distances, fixing $H_0$ and $M_B$ to reference values provided by the collaboration. By looking at 500 different patches randomly taken across the sky, we find a maximum $\sim 3\sigma$ CL anisotropy level for $q_0$, whose direction points orthogonally to the CMB dipole axis, i.e., $(RA^{\rm SN},DEC^{\rm SN}) = (267^{\circ},6^{\circ})$ vs $(RA^{\rm CMB},DEC^{\rm CMB}) = (167^{\circ},-7^{\circ})$. We assessed the statistical significance of those results, finding that such a signal is expected due to the limitations of the observational sample. These results support that there is no significant evidence for a departure from the cosmic isotropy assumption, one of the pillars of the standard cosmological model.

Cosmological correlators, the natural observables of the primordial universe, have been extensively studied in the past two decades using the in-in formalism pioneered by Schwinger and Keldysh for the study of dissipative open systems. Ironically, most applications in cosmology have focused on non-dissipative closed systems. We show that, for non-dissipative systems, correlators can be equivalently computed using the in-out formalism with the familiar Feynman rules. In particular, the myriad of in-in propagators is reduced to a single (Feynman) time-ordered propagator and no sum over the labelling of vertices is required. In de Sitter spacetime, this requires extending the expanding Poincar\'e patch with a contracting patch, which prepares the bra from the future. Our results are valid for fields of any mass and spin but assuming the absence of infrared divergences. We present three applications of the in-out formalism: a representation of correlators in terms of a sum over residues of Feynman propagators in the energy-momentum domain; an algebraic recursion relation that computes Minkowski correlators in terms of lower order ones; and the derivation of cutting rules from Veltman's largest time equation, which we explicitly develop and exemplify for two-vertex diagrams to all loop orders. The in-out formalism leads to a natural definition of a de Sitter scattering matrix, which we discuss in simple examples. Remarkably, we show that our scattering matrix satisfies the standard optical theorem and the positivity that follows from it in the forward limit.

In this work we present general predictions for the static observables of neutron stars (NSs) under the hypothesis of a purely nucleonic composition of the ultra-dense baryonic matter, using Bayesian inference on a very large parameter space conditioned by both astrophysical and nuclear physics constraints. The equation of states are obtained using a unified approach of the NS core and inner crust within a fully covariant treatment based on a relativistic mean-field Lagrangian density with density dependent couplings. The posterior distributions are well compatible with the ones obtained by semi-agnostic meta-modelling techniques based on non-relativistic functionals, that span a similar portion of the parameter space in terms of nuclear matter parameters, and we confirm that the hypothesis of a purely nucleonic composition is compatible with all the present observations. We additionally show that present observations do not exclude the existence of very massive neutron stars with mass compatible with the lighter partner of the gravitational event GW190814 measured by the LIGO-Virgo collaboration. Some selected representative models, that respect well all the constraints taken into account in this study, and approximately cover the residual uncertainty in our posterior distributions, will be uploaded in the CompOSE database for use by the community.

In this study, we develop a numerical method to generate images on an observer's screen, formed by radiation from hotspots on any timelike orbits outside a black hole. This method uses the calculation of fractional numbers, enabling us not only to produce the overall image but also to distinguish between primary, secondary, and higher-order images. Building upon this, we compute the images of hotspots from eight potential types of geodesic timelike orbits outside a Kerr black hole, summarizing the properties of both the overall and individual order images. Furthermore, we calculate the centroid motion and lightcurve. Notably, we observe flare phenomena across all orbit types and classify these flares into three categories based on the Doppler and gravitational redshift effects.

We present a statistical analysis of electrostatic solitary waves observed aboard Magnetospheric Multiscale spacecraft in the Earth's magnetosheath. Applying single-spacecraft interferometry to several hundred solitary waves collected in about two minute intervals, we show that almost all of them have the electrostatic potential of positive polarity and propagate quasi-parallel to the local magnetic field with plasma frame velocities of the order of 100 km/s. The solitary waves have typical parallel half-widths from 10 to 100 m that is between 1 and 10 Debye lengths and typical amplitudes of the electrostatic potential from 10 to 200 mV that is between 0.01 and 1\% of local electron temperature. The solitary waves are associated with quasi-Maxwellian ion velocity distribution functions, and their plasma frame velocities are comparable with ion thermal speed and well below electron thermal speed. We argue that the solitary waves of positive polarity are slow electron holes and estimate the time scale of their acceleration, which occurs due to interaction with ions, to be of the order of one second. The observation of slow electron holes indicates that their lifetime was shorter than the acceleration time scale. We argue that multi-spacecraft interferometry applied previously to these solitary waves is not applicable because of their too-short spatial scales. The source of the slow electron holes and the role in electron-ion energy exchange remain to be established.

We consider the effects of a bare mass term for the inflaton, when the inflationary potential takes the form $V(\phi)= \lambda \phi^k$ about its minimum with $k \ge 4$. We concentrate on $k=4$, but discuss general cases as well. Further, we assume $\lambda \phi_{\rm end}^2 \gg m_\phi^2$, where $\phi_{\rm end}$ is the inflaton field value when the inflationary expansion ends. We show that the presence of a mass term (which may be present due to radiative corrections or supersymmetry breaking) can significantly alter the reheating process, as the equation of state of the inflaton condensate changes from $w_\phi=\frac{1}{3}$ to $w_\phi=0$ when $\lambda \phi^2$ drops below $m_\phi^2$. We show that for a mass $m_\phi \gtrsim T_{\rm RH}/250$, the mass term will dominate at reheating. We compute the effects on the reheating temperature for cases where reheating is due to inflaton decay (to fermions, scalars, or vectors) or to inflaton scattering (to scalars or vectors). For scattering to scalars and in the absence of a decay, we derive a strong upper limit to the inflaton bare mass $m_\phi < 350~{\rm MeV} (T_{\rm RH}/10^{10}~{\rm GeV})^{3/5}$, as there is always a residual inflaton background which acts as cold dark matter. We also consider the effect of the bare mass term on the fragmentation of the inflaton condensate.

Characterizing the thermodynamics of turbulent plasmas is key to decoding observable signatures from astrophysical systems. In magnetohydrodynamic (MHD) turbulence, nonlinear interactions between counter-propagating Alfv\'en waves cascade energy to smaller spatial scales where dissipation heats the protons and electrons. When the thermal pressure far exceeds the magnetic pressure, linear theory predicts a spectral gap at perpendicular scales near the proton gyroradius where Alfv\'en waves become non-propagating. For simple models of an MHD turbulent cascade that assume only local nonlinear interactions, the cascade halts at this gap, preventing energy from reaching smaller scales where electron dissipation dominates, leading to an overestimate of the proton heating rate. In this work, we demonstrate that nonlocal contributions to the cascade, specifically large scale shearing and small scale diffusion, can bridge the non-propagating gap, allowing the cascade to continue to smaller scales. We provide an updated functional form for the proton-to-electron heating ratio accounting for this nonlocal energy transfer by evaluating a nonlocal weakened cascade model over a range of temperature and pressure ratios. In plasmas where the thermal pressure dominates the magnetic pressure, we observe that the proton heating is moderated compared to the significant enhancement predicted by local models.

In this study, we propose an investigation into dark photon dark matter (DPDM) within the infrared frequency band, utilizing highly sensitive infrared light detectors commonly integrated into space telescopes, such as the James Webb Space Telescope (JWST). The presence of DPDM induces electron oscillations in the reflector of these detectors. Consequently, these oscillating electrons can emit monochromatic electromagnetic waves with a frequency almost equivalent to the mass of DPDM. By employing the stationary phase approximation, we can demonstrate that when the size of the reflector significantly exceeds the wavelength of the electromagnetic wave, the contribution to the electromagnetic wave field at a given position primarily stems from the surface unit perpendicular to the relative position vector. This simplification results in the reduction of electromagnetic wave calculations to ray optics. By applying this concept to JWST, our analysis of observational data demonstrates the potential to establish constraints on the kinetic mixing between the photon and dark photon within the range [10, 500] THz. Despite JWST not being optimized for DPDM searches, our findings reveal constraints comparable to those obtained from the XENON1T experiment in the laboratory, as well as astrophysical constraints from solar emission. Additionally, we explore strategies to optimize future experiments specifically designed for DPDM searches.

We recently hypothesized that a distortion parameter exists such that its signed sum for all images of singular gravitational lensing of a source vanishes identically [1]. We found a distortion parameter (the ratio of the tangential to radial magnifications) that satisfied the hypothesis for the images of Schwarzschild lensing. We now show that another distortion parameter (the difference between tangential and radial magnifications) also supports our hypothesis when we perform computations with the primary-secondary and relativistic images. The distortion parameters, which satisfy the aesthetically appealing hypothesis, will likely aid in developing gravitational lensing theory. Finally, we discuss the conservation of distortions of images in gravitational lensing.

We reappraise the viability of asymmetric dark matter (ADM) realized as a Dirac fermion coupling dominantly to the Standard Model fermions. Treating the interactions of such a DM particle with quarks/leptons in an effective-interactions framework, we derive updated constraints using mono-jet searches from the Large Hadron Collider (LHC) and mono-photon searches at the Large Electron-Positron (LEP) collider. We carefully model the detectors used in these experiments, which is found to have significant impact. The constraint of efficient annihilation of the symmetric part of the ADM, as well as other observational constraints are synthesized to produce a global picture. Consistent with previous work, we find that ADM with mass in the range $1-100$ GeV is strongly constrained, thus ruling out its best motivated mass range. However, we find that leptophilic ADM remains allowed for $\gtrsim 10$ GeV DM, including bounds from colliders, direct detection, and stellar heating. We forecast that the Future Circular Collider for electron-positron collisions (FCC-ee) will improve sensitivity to DM-lepton interactions by almost an order of magnitude.

Recent Pulsar Timing Array (PTA) collaborations show strong evidence for a stochastic gravitational wave background (SGWB) with the characteristic Hellings-Downs inter-pulsar correlations. The signal may stem from supermassive black hole binary mergers, or early universe phenomena. The former is expected to be strongly anisotropic while primordial backgrounds are likely to be predominantly isotropic with small fluctuations. In case the observed SGWB is of cosmological origin, our relative motion with respect to the SGWB rest frame is a guaranteed source of anisotropy, leading to $\mathcal{O}(10^{-3})$ energy density fluctuations of the SGWB. These kinematic anisotropies are likely to be larger than the intrinsic anisotropies, akin to the cosmic microwave background (CMB) dipole anisotropy. We assess the sensitivity of current PTA data to the kinematic dipole anisotropy and also provide forecasts with which the magnitude and direction of the kinematic dipole may be measured in the future with an SKA-like experiment. We also discuss how the spectral shape of the SGWB and the location of pulsar observed affects the prospects of detecting the kinematic dipole with PTA. A detection of this anisotropy may even help resolve the discrepancy in the magnitude of the kinematic dipole as measured by CMB and large-scale structure observations.

We report correlations in underground seismic measurements with horizontal separations of several hundreds of meters to a few kilometers in the frequency range 0.01Hz to 40Hz. These seismic correlations could threaten science goals of planned interferometric gravitational-wave detectors such as the Einstein Telescope as well as atom interferometers such as MIGA and ELGAR. We use seismic measurements from four different sites, i.e. the former Homestake mine (USA) as well as two candidate sites for the Einstein Telescope, Sos Enattos (IT) and Euregio Maas-Rhein (NL-BE-DE) and the site housing the MIGA detector, LSBB (FR). At all sites, we observe significant coherence for at least 50% of the time in the majority of the frequency region of interest. Based on the observed correlations in the seismic fields, we predict levels of correlated Newtonian noise from body waves. We project the effect of correlated Newtonian noise from body waves on the capabilities of the triangular design of the Einstein Telescope's to observe an isotropic gravitational-wave background (GWB) and find that, even in case of the most quiet site, its sensitivity will be affected up to $\sim$20Hz. The resolvable amplitude of a GWB signal with a negatively sloped power-law behaviour would be reduced by several orders of magnitude. However, the resolvability of a power-law signal with a slope of e.g. $\alpha=0$ ($\alpha=2/3$) would be more moderately affected by a factor $\sim$ 6-9 ($\sim$3-4) in case of a low noise environment. Furthermore, we bolster confidence in our results by showing that transient noise features have a limited impact on the presented results.

If two particles moving towards a black hole collide in the vicinity of the horizon, the energy $E_{c.m.}$ in the center of mass frame can grow indefinitely if one of particles is fine-tuned. This is the Ba\~{n}ados, Silk and West (BSW) effect. One of objections against this effect consists in that for some types of a horizon fine-tuned particles cannot reach the horizon. However, this difficulty can be overcome if instead of exact fine-tuning, one of particle is nearly fine-tuned, with the value of small detuning being adjusted to the distance to the horizon. Such particles are called near-fine-tuned. We give classification of such particles and describe possible high energy scenarios of collision in which they participate. We analyze the ranges of possible motion for each type of particle and determine under which condition such particles can reach the horizon. We analyze collision energy $E_{c.m.}\,$and determine under which conditions it may grow indefinitely. We also include into consideration the forces acting on particles and find when the BSW effect with nearly-fine-tuned particles is possible with finite forces. We demonstrate that the BSW effect with particles under discussion is consistent with the principle of kinematic censorship. According to this principle, $E_{c.m.}\,$cannot be literally infinite in any event of collision (if no singularity is present), although it can be made as large as one likes.

PBH formation requires high-density regions in the (random) density field filling the primordial universe. While only the largest (and so rarest) overdensities collapse to form PBHs, the rest cause large anisotropic stresses, which are the source of GWs. We provide an overview of the theoretical aspects of the GW backgrounds associated with PBHs from large primordial fluctuations. We consider GW backgrounds associated with PBH formation, PBH reheating and unresolved PBH binaries. We present several graphical summaries and illustrations for the busy reader.

We explore the idea of restoring the full diffeomorphism (Diff) invariance in theories with only transverse diffeomorphisms (TDiff) by the introduction of additional fields. In particular, we consider in detail the case of a TDiff invariant scalar field and how Diff symmetry can be restored by introducing an additional vector field. We reobtain the corresponding dynamics and energy-momentum tensor from the covariantized action and analyze the potential and kinetic domination regimes. For the former, the theory describes a cosmological constant-type behaviour, while for the latter we show that the theory can describe an adiabatic perfect fluid whose equation of state and speed of sound is obtained in a straightforward way.