Papers for Wednesday, Apr 19 2023

Deep VLT/MUSE optical integral field spectroscopy has recently revealed an abundant population of ultra-faint galaxies ($M_{UV} = -$15; 0.01 $L_{\star}$) at $z=$2.9$-$6.7 due to their strong Lyman-$\alpha$ emission. The implied Lyman-$\alpha$ equivalent widths are in excess of 100-200 Angstrom, challenging existing models of normal star formation and implying extremely young ages, small stellar masses, and a very low amount of metal enrichment. We use JWST/NIRSpec's microshutter array to follow-up 45 of these galaxies (11h in G235M/F170LP and 7h in G395M/F290LP), as well as 45 lower-equivalent width Lyman-$\alpha$ emitters. Our spectroscopy covers the range 1.7$-$5.1 micron in order to target strong optical emission lines: H$\alpha$, [OIII], H$\beta$, and [NII]. Individual measurements as well as stacks reveal line ratios consistent with a metal poor nature (2$-$30% $Z_{\odot}$) and intense ionizing radiation fields. The galaxies with the highest equivalent widths of Lyman-$\alpha$, in excess of 120 Angstrom, have lower gas-phase metallicities than those with lower equivalent widths. This implies a selection based on Lyman-$\alpha$ equivalent width is an efficient technique for identifying younger, less chemically enriched systems.

We discuss our recent constraints on the coupling of Very-Light Axion-Like Particles (of masses $<$$ 10^{-12} \ \mathrm{eV}$) to electromagnetism from $Chandra$ observations of the cluster-hosted Active Galactic Nuclei (AGN) H1821+643 and NGC1275. In both cases, the inferred high-quality AGN spectra excluded all photon-ALP couplings $g_\mathrm{a\gamma} > (6.3 - 8.0) \times 10^{-13} \ {\mathrm{GeV}}^{-1}$ at the $99.7\%$ level, respectively, based on the non-detection of spectral distortions attributed to photon-ALP inter-conversion along the cluster line-of-sight. Finally, we present the prospects of tightening current bounds on such ALPs by up to a factor of 10 with next-generation X-ray observatories such as $Athena$, $AXIS$ and $LEM$ given their improved spectral and spatial resolution and collecting area compared to current missions.

The transient X-ray source MAXI J1810-222 was discovered in 2018 and has been active ever since. A long combined radio and X-ray monitoring campaign was performed with ATCA and Swift respectively. It has been proposed that MAXI J1810-222 is a relatively distant black hole X-ray binary, albeit showing a very peculiar outburst behaviour. Here, we report on the spectral study of this source making use of a large sample of NICER observations performed between 2019 February and 2020 September. We detected a strong spectral absorption feature at $\sim$1 keV, which we have characterised with a physical photoionisation model. Via a deep scan of the parameters space, we obtained evidence for a spectral-state dependent outflow, with mildly relativistic speeds. In particular, the soft and intermediate states point to a hot plasma outflowing at 0.05-0.15 $c$. This speeds rule-out thermal winds and, hence, they suggest that such outflows could be radiation pressure or (most likely) magnetically-driven winds. Our results are crucial to test current theoretical models of wind formation in X-ray binaries.

We analyse the gas-phase metallicity properties of a sample of 66 low stellar mass (log M*/M_sun < 8.5) galaxies at 3<z<10, observed with JWST/NIRSpec as part of the JADES programme in its deep GOODS-S tier. By combining this sample with more massive galaxies at similar redshifts from other programmes, we study the scaling relations between stellar mass (M*), oxygen abundance (O/H), and star-formation rate (SFR) across three orders of magnitude in mass out to the epoch of early galaxy assembly. We find evidence for a shallower slope at the low-mass-end of the mass-metallicity relation (MZR), with 12 + log(O/H) = (7.88 +- 0.03) + (0.17 +- 0.04) log(M*/10^8 M_sun), in good agreement with the MZR probed by local analogues of high-redshift systems like `green pea' and `blueberry' galaxies. The inferred slope is well matched by models including `momentum-driven' SNe winds, suggesting that feedback mechanisms in dwarf galaxies (and at high-z) might be different from those in place at higher masses. The evolution in the normalisation is instead observed to be relatively mild compared to previous determinations of the MZR at z~3 (~0.1-0.2 dex on average). We also find a progressive deviation from the local fundamental metallicity relation (FMR) as a function of redshift, especially at z>6, with galaxies significantly less enriched (by ~0.4 dex on average) than predicted given their M* and SFR. These observations are consistent with an enhanced stochasticity in the accretion rate from the cosmic web, and/or with an increased efficiency in metal removals by outflows, prompting us to reconsider the nature of the relationship between M*, O/H, and SFR in the early Universe.

A population of massive galaxies invisible or very faint in deep optical/near-infrared surveys, but brighter at longer wavelengths has been uncovered in the past years. However, the nature of these optically dark/faint galaxies (OFGs, among other nomenclatures) is highly uncertain. In this work, we investigate the drivers of dust attenuation in the JWST era. Particularly, we study the role of stellar mass, size, and orientation in obscuring star-forming galaxies (SFGs) at $3 < z < 7.5$, focusing on understanding why galaxies like OFGs are so faint at optical/near-infrared wavelengths. We find that stellar mass is the primary proxy of dust attenuation among those studied. Effective radius and axis ratio do not show a clear link with dust attenuation, with the effect of orientation close to random. However, there is a subset of highly dust attenuated ($A_V > 1$, typically) SFGs, of which OFGs are a specific case. For this subset, we find that the key distinctive feature is their compact size (for massive systems with $\log (M_{*}/M_{\odot}) > 10$), exhibiting 30% smaller effective radius than the average SFGs at the same stellar mass and redshift. On the contrary, they do not exhibit a preference for low axis ratios (i.e., edge-on disks). The results in this work show stellar mass as a primary proxy of dust attenuation and compact stellar light profiles behind thick dust columns obscuring typical massive SFGs.

Extremely large telescopes (ELTs) present an unparalleled opportunity to study the magnetism, atmospheric dynamics, and chemistry of very low mass stars (VLMs), brown dwarfs, and exoplanets. Instruments such as the Giant Magellan Telescope - Consortium Large Earth Finder (GMT/GCLEF), the Thirty Meter Telescope's Multi-Objective Diffraction-limited High-Resolution Infrared Spectrograph (TMT/MODHIS), and the European Southern Observatory's Mid-Infrared ELT Imager and Spectrograph (ELT/METIS) provide the spectral resolution and signal-to-noise (S/N) necessary to Doppler image ultracool targets' surfaces based on temporal spectral variations due to surface inhomogeneities. Using our publicly-available code, $Imber$, developed and validated in Plummer & Wang (2022), we evaluate these instruments' abilities to discern magnetic star spots and cloud systems on a VLM star (TRAPPIST-1); two L/T transition ultracool dwarfs (VHS J1256$-$1257 b and SIMP J0136+0933); and three exoplanets (Beta Pic b and HR 8799 d and e). We find that TMT/MODHIS and ELT/METIS are suitable for Doppler imaging the ultracool dwarfs and Beta Pic b over a single rotation. Uncertainties for longitude and radius are typically $\lesssim 10^{\circ}$, and latitude uncertainties range from $\sim 10^{\circ} \ \rm{to} \ 30^{\circ}$. TRAPPIST-1's edge-on inclination and low $\upsilon \sin i$ provide a challenge for all three instruments while GMT/GCLEF and the HR 8799 planets may require observations over multiple rotations. We compare the spectroscopic technique, photometry-only inference, and the combination of the two. We find combining spectroscopic and photometric observations can lead to improved Bayesian inference of surface inhomogeneities and offers insight into whether ultracool atmospheres are dominated by spotted or banded features.

We investigate dust attenuation and its dependence on viewing angle for 308 star-forming galaxies at $1.3\leq z\leq2.6$ from the MOSFIRE Deep Evolution Field (MOSDEF) survey. We divide galaxies with a detected H$\alpha$ emission line and coverage of H$\beta$ into eight groups by stellar mass, star formation rate (SFR), and inclination (i.e., axis ratio), then stack their spectra. From each stack, we measure Balmer decrement and gas-phase metallicity, then we compute median \AV and UV continuum spectral slope ($\beta$). First, we find that none of the dust properties (Balmer decrement, \AV, $\beta$) vary with axis ratio. Second, both stellar and nebular attenuation increase with increasing galaxy mass, showing little residual dependence on SFR or metallicity. Third, nebular emission is more attenuated than stellar emission, and this difference grows even larger at higher galaxy masses and SFRs. Based on these results, we propose a three-component dust model where attenuation predominantly occurs in star-forming regions and large, dusty star-forming clumps, with minimal attenuation in the diffuse ISM. In this model, nebular attenuation primarily originates in clumps, while stellar attenuation is dominated by star-forming regions. Clumps become larger and more common with increasing galaxy mass, creating the above mass trends. Finally, we argue that a fixed metal yield naturally leads to mass regulating dust attenuation. Infall of low-metallicity gas increases SFR and lowers metallicity, but leaves dust column density mostly unchanged. We quantify this idea using the Kennicutt-Schmidt and fundamental metallicity relations, showing that galaxy mass is indeed the primary driver of dust attenuation.

We present a dynamical study of the eclipsing intermediate polar XY Ari based on time-resolved near-infrared spectroscopy obtained with the EMIR spectrograph on the 10.4-m Gran Telescopio Canarias. Using main sequence template spectra taken with the same instrument setup as the target spectra, we measure a radial velocity amplitude of the late K-type donor star $K_2=256 \pm 2$ km s$^{-1}$. We also obtain the rotational broadening of its photospheric lines $v_\mathrm{rot} \sin i = 141 \pm 3$ km s$^{-1}$. From these and the eclipse geometry, we derive a donor-to-white dwarf mass ratio $q = M_2/M_1 = 0.62 \pm 0.02$, an orbital inclination $i = 80.8^{\circ} \pm 0.5^{\circ}$ and dynamical masses $M_{1} = 1.21 \pm 0.04 \, \mathrm{M}_{\odot}$ and $M_2 = 0.75 \pm 0.04 \, \mathrm{M}_{\odot}$ ($1 \sigma$). This result places the white dwarf in XY Ari as one of the three most massive known in a cataclysmic variable. Comparison with a number of white dwarf mass estimates from X-ray spectral modelling indicates the necessity of a major revision of the cooling models currently assumed for magnetic accretion in intermediate polars, as most of the X-ray white dwarf masses lie significantly below the dynamical mass value.

Reliable indirect diagnostics of LyC photon escape from galaxies are required to understand which sources were the dominant contributors to reionization. While multiple escape fraction ($f_{\rm esc}$) indicators have been proposed to trace favourable conditions for LyC leakage from the interstellar medium of low-redshift ``analog'' galaxies, it remains unclear whether these are applicable at high redshifts where LyC emission cannot be directly observed. Using a library of 14,120 mock spectra of star-forming galaxies with redshifts $4.64 \leq z \leq 10$ from the SPHINX$^{20}$ cosmological radiation hydrodynamics simulation, we develop a framework for the physics that leads to high $f_{\rm esc}$. We investigate LyC leakage from our galaxies based on the criteria that successful LyC escape diagnostics must \textit{i)} track a high specific star formation rate, \textit{ii)} be sensitive to stellar population age in the range $3.5-10$~Myr representing the times when supernova first explode to when LyC production significantly drops, and \textit{iii)} include a proxy for neutral gas content and gas density in the interstellar medium. ${\rm O}_{32}$, $\Sigma_{\rm SFR}$, M$_{\rm UV}$, and H$\beta$ equivalent width select for one or fewer of our criteria, rendering them either necessary but insufficient or generally poor diagnostics. In contrast, UV slope ($\beta$), and ${\rm E(B-V)}$ match two or more of our criteria, rendering them good $f_{\rm esc}$ diagnostics (albeit with significant scatter). Using our library, we build a quantitative model for predicting $f_{\rm esc}$ based on $\beta$, ${\rm E(B-V)}$, H$\beta$, M$_{\rm UV}$, ${\rm R_{23}}$, and ${\rm O_{32}}$. When applied to bright $z > 6$ Ly$\alpha$ emitters observed with JWST, we find that the majority of them have $f_{\rm esc} \lesssim 10\%$.

We present deep optical imaging and photometry of four objects classified as "Almost-Dark" galaxies in the ALFALFA survey because of their gas-rich nature and extremely faint or missing optical emission in existing catalogs. They have HI masses of $10^7$-$10^9$ $M_{\odot}$ and distances of $\sim$9-100 Mpc. Observations with the WIYN 3.5m telescope and One Degree Imager reveal faint stellar components with central surface brightnesses of $\sim$24-25 $\mathrm{mag}\,\mathrm{arcsec}^{-2}$ in the g-band. We also present the results of HI synthesis observations with the Westerbork Synthesis Radio Telescope. These Almost-Dark galaxies have been identified as possible tidal dwarf galaxies (TDGs) based on their proximity to one or more massive galaxies. We demonstrate that AGC 229398 and AGC 333576 likely have the low dark matter content and large effective radii representative of TDGs. They are located much farther from their progenitors than previously studied TDGs, suggesting they are older and more evolved. AGC 219369 is likely dark matter dominated, while AGC 123216 has a dark matter content that is unusually high for a TDG, but low for a normal dwarf galaxy. We consider possible mechanisms for the formation of the TDG candidates such as a traditional major merger scenario and gas ejection from a high velocity fly-by. Blind HI surveys like ALFALFA enable the detection of gas-rich, optically faint TDGs that can be overlooked in other surveys, thereby providing a more complete census of the low-mass galaxy population and an opportunity to study TDGs at a more advanced stage of their life cycle.

Particle lateral distributions were investigated in cosmic ray air showers with energy $E_0 \simeq 10^{18}$ eV registered at the Yakutsk array with surface and underground scintillation detectors with $\simeq 1 \times \sec\theta$~GeV threshold during the period of continuous observations from 1986 to 2016. The analysis covers events with arrival direction zenith angles $\theta \le 60^{\circ}$ within five intervals with step $\Delta\cos\theta = 0.1$. Experimental values were compared to simulation results obtained with the use of CORSIKA code within the framework of QGSJet01 hadron interaction model. The whole dataset points at probable cosmic ray composition which is close to protons.

We report ALMA high-resolution observations of water emission lines $p-{\rm{H_2O}} (2_{02}-1_{11}$), $o-{\rm{H_2O}} (3_{21}-3_{12})$, $p-{\rm{H_2O}} (4_{22}-4_{13})$, in the strongly lensed galaxy HATLASJ113526.2-01460 at redshift z $\sim$ 3.1. From the lensing-reconstructed maps of water emission and line profiles, we infer the general physical properties of the ISM in the molecular clouds where the lines arise. We find that the water vapor lines $o-{\rm{H_2O}} (3_{21}-3_{12})$, $p-{\rm{H_2O}} (4_{22}-4_{13})$ are mainly excited by FIR pumping from dust radiation in a warm and dense environment, with dust temperatures ranging from 70 K to $\sim 100$ K, as suggested by the line ratios. The $p-{\rm{H_2O}} (2_{02}-1_{11})$ line instead, is excited by a complex interplay between FIR pumping and collisional excitation in the dense core of the star-forming region. This scenario is also supported by the detection of the medium-level excitation of CO resulting in the line emission CO (J=8-7). Thanks to the unprecedented high resolution offered by the combination of ALMA capabilities and gravitational lensing, we discern the different phases of the ISM and locate the hot molecular clouds into a physical scale of $\sim$ 500 pc. We discuss the possibility of J1135 hosting an AGN in its accretion phase. Finally, we determine the relation between the water emission lines and the total IR luminosity of J1135, as well as the SFR as a function of water emission intensities, comparing the outcomes to local and high-$z$ galactic samples from the literature.

We report results of an observing campaign to study the dwarf M flare star EV Lacertae. Between October 2021 and January 2022 we obtained concurrent B band photometry and low resolution spectroscopy of EV Lac on 39 occasions during 10 of which we observed flares with amplitude greater than 0.1 magnitude. Spectra were calibrated in absolute flux using concurrent photometry and flare-only spectra obtained by subtracting mean quiescent spectra. We measured B band flare energies between Log E = 30.8 and 32.6 erg. In the brightest flares we measured temporal development of flare flux in H I and He I emission lines and in the adjacent continuum and found that flux in the continuum subsided more rapidly than in the emission lines. Although our time resolution was limited, in our brightest flare we saw flux in the continuum clearly peaking before flux in the emission lines. We observed a progressive decrease in flare energy from H\b{eta} to H{\delta}. On average we found 37% of B band flare energy appeared in the H\b{eta} to H{\epsilon} emission lines with the remainder contributing to a rise in continuum flux. We measured black-body temperatures for the brightest flares between 10,500 +- 700 K and 19,500 +- 500 K and found a linear relationship between flare temperature and continuum flux at 4170 {\AA}. Balmer lines in flare-only spectra were well fitted by Gaussian profiles with some evidence of additional short-lived blue-shifted emission at the flare peak.

Does radiometry (e.g., signal-to-noise ratio) limit the performance of near-IR subcloud imaging of our sister planet's surface at night? It does not. We compute subcloud radiometry using above-cloud observations, an assumed ground temperature, sub-cloud absorption and emission modeling, and Rayleigh scattering simulations. We thus confirm both archival and recent studies that deployment of a modest subcloud camera does enable high-resolution surface imaging.

We present a new approach to measuring the thickness of a partially face-on stellar disk, using dust geometry. In a moderately-inclined disk galaxy, the fraction of reddened stars is expected to be 50% everywhere, assuming that dust lies in a thin midplane. In a thickened disk, however, a wide range of radii project onto the line of sight. Assuming stellar density declines with radius, this geometrical projection leads to differences in the numbers of stars on the near and far sides of the thin dust layer. The fraction of reddened stars will thus differ from the 50% prediction, with a deviation that becomes larger for puffier disks. We map the fraction of reddened red giant branch (RGB) stars across M31, which shows prominent dust lanes on only one side of the major axis. The fraction of reddened stars varies systematically from 20% to 80%, which requires that these stars have an exponential scale height h_z that is 0.14+/-0.015 times the exponential scale length (h_r~5.5kpc). M31's RGB stars must therefore have h_z=770+/-80pc, which is far thicker than the Milky Way's thin disk, but comparable to its thick disk. The lack of a significant thin disk in M31 is unexpected, but consistent with its interaction history and high disk velocity dispersion. We suggest that asymmetric reddening be used as a generic criteria for identifying ``thick disk'' dominated systems, and discuss prospects for future 3-dimensional tomographic mapping of the gas and stars in M31.

Context. Open clusters (OCs) are fundamental benchmarks where theories of star formation and stellar evolution can be tested and validated. Coma Ber and Latyshev 2 (Group X) are the second and third OCs closest to the Sun, making them excellent targets to search for low-mass stars and ultra-cool dwarfs. In addition, this pair will experience a flyby in 10-16 Myr which makes it a benchmark to test OCs pair interactions. Aims. We aim at analysing the membership, luminosity, mass, phase-space (i.e., positions and velocities), and energy distributions for Coma Ber and Latyshev 2 and test the hypothesis of the mixing of their populations at the encounter time. Methods. We develop a new phase-space membership methodology and apply it to Gaia data. With the recovered members we infer the phase-space, luminosity and mass distributions using publicly available Bayesian inference codes. Then, with a publicly available orbit integration code and members' positions and velocities, we integrate their orbits 20 Myr into the future. Results. In Coma Ber, we identify 302 candidate members distributed in the core and tidal tails. The tails are dynamically cold and asymmetrically populated. The stellar system called Group X is made of two structures: the disrupted OC Latyshev 2 (186 candidate members) and a loose stellar association called Mecayotl 1 (146 candidate members), both of them will fly by Coma Ber in $11.3\pm0.5$ Myr and $14.0\pm0.6$ Myr, respectively, and each other in $8.1\pm1.3$ Myr. Conclusions. We study the dynamical properties of the core and tails of Coma Ber and also confirm the existence of the OC Latyshev 2 and its neighbour stellar association Mecayotl 1. Although these three systems will experience encounters we find no evidence supporting the mixing of their populations.

Gamma-ray bursts (GRBs) have been traditionally classified based on their duration. The increasing number of extended emission (EE) GRBs, lasting typically more than 2 seconds but with properties similar to those of a short GRBs, challenges the traditional classification criteria. In this work, we use the t-Distributed Stochastic Neighbor Embedding (t-SNE), a machine learning technique, to classify GRBs. We present the results for GRBs observed until July 2022 by the {\itshape Swift}/BAT instrument in all its energy bands. We show the effects of varying the learning rate and perplexity parameters as well as the benefit of pre-processing the data by a non-parametric noise reduction technique. % named {\sc FABADA}. Consistently with previous works, we show that the t-SNE method separates GRBs in two subgroups. We also show that EE GRBs reported by various authors under different criteria tend to cluster in a few regions of our t-SNE maps, and identify seven new EE GRB candidates by using the gamma-ray data provided by the automatic pipeline of {\itshape Swift}/BAT and the proximity with previously identified EE GRBs.

We present the first 3D radiation-hydrodynamic simulations of the formation of planetary nebulae (PNe) emerging from 3D spiral patterns. We use the GUACHO code to create 3D spiral structures as a consequence of the distortions on the geometry of the intrinsically isotropic wind of an asymptotic giant branch (AGB) star produced by a companion star in a circular orbit. We found that the orbital period of the binary producing the 3D spiral pattern has consequences on the formation and shaping of the PN itself. Stellar systems with longer period create less entwined 3D spirals, producing PNe with rounder inner cavities, and prevent the expansion of jet towards the polar directions. The spiral fitting procedure used in the literature to predict the binary's orbital period may be misleading in the case of proto-PNe and PNe as spiral patterns are diluted by their own thermal expansion down to the average AGB density profile within a few hundred years and are further disrupted by the action of jets. By adopting a phase of jet ejections between the AGB and post-AGB stages, we are able to recover the morphologies of proto-PNe and PNe that exhibit ring-like structures in their halos.

In this paper, we compare the chemistry and the emission spectra of nitrogen-dominated cool, warm, and hot ultra-short-period (USP) super-Earth atmospheres in and out of chemical equilibrium at various surface pressure scenarios ranging from 0.1 to 10 bar. We link the one-dimensional VULCAN chemical kinetic code, in which thermochemical kinetic and vertical transport and photochemistry are taken into account, to the one-dimensional radiative transfer model, PETITRADTRANS, to predict the emission spectra of these planets. The radiative-convective temperature-pressure profiles were computed with the HELIOS code. Then, using PANDEXO noise simulator, we explore the observability of the differences produced by disequilibrium processes with the JWST. Our grids show how different surface pressures can significantly affect the temperature profiles, the atmospheric abundances, and consequently the emission spectra of these planets. We find that the divergences due to disequilibrium processes would be possible to observe in cooler planets by targeting HCN, C2H4, and CO, and in warmer planets by targeting CH4 with HCN, using the NIRSpec and MIRI LRS JWST instruments. These species are also found to be sensitive indicators of the existence of surfaces on nitrogen-dominated USP super-Earths, providing information regarding the thickness of these atmospheres.

Magnetic reconnection is a key mechanism involved in solar eruptions and is also a prime possibility to heat the low corona to millions of degrees. Here, we present ultra-high-resolution extreme ultraviolet observations of persistent null-point reconnection in the corona at a scale of about 390 km over one hour observations of the Extreme-Ultraviolet Imager on board Solar Orbiter spacecraft. The observations show formation of a null-point configuration above a minor positive polarity embedded within a region of dominant negative polarity near a sunspot. The gentle phase of the persistent null-point reconnection is evidenced by sustained point-like high-temperature plasma (about 10 MK) near the null-point and constant outflow blobs not only along the outer spine but also along the fan surface. The blobs appear at a higher frequency than previously observed with an average velocity of about 80 km/s and life-times of about 40 s. The null-point reconnection also occurs explosively but only for 4 minutes, its coupling with a mini-filament eruption generates a spiral jet. These results suggest that magnetic reconnection, at previously unresolved scales, proceeds continually in a gentle and/or explosive way to persistently transfer mass and energy to the overlying corona.

An interesting result recently reported for Mrk 421 is the detection of a significant excess at hard X-ray energies, which could provide useful information for investigating particle acceleration and emission mechanisms in the relativistic jet. Considering a two-injection scenario, we develop a self-consistent one-zone leptonic model to understand the origin of the hard X-ray excess in Mrk 421 during the period of extremely low X-ray and very high energy (VHE) flux in 2013 January. In the model, two populations of mono-energetic ultrarelativistic electrons are injected into the emission region that is a magnetized plasmoid propagating along the blazar jet. We numerically calculate the emitting electron energy distribution by solving a kinetic equation that incorporates both shock acceleration and stochastic acceleration processes. Moreover, we infer analytic expressions relating electrons acceleration, cooling, escape and injection to the observed spectra and variability. For the injection luminosity in particular, we derive a new approximate analytical expression for the case of continual injection with a mono-energetic distribution. Based on a comparison between the theoretical predictions and the observed SED, we conclude that the hard X-ray excess observed in Mrk 421 may be due to the synchrotron radiation emitted by an additional electrons population, which is co-spatial with an electron population producing simultaneous Optical/UV, soft X-ray, and $\gamma$-ray emissions. The stochastic acceleration may play a major role in producing the observed X-ray spectrum.

The Sunyaev Zel'dovich (SZ) effect is sensitive to the pressure of ionized gas inside galaxy clusters. The gas pressure responds to changes in the gravitational potential of the cluster, which is dominated by the host dark matter halo. Changes in halo concentration therefore impact the SZ signal, with implications for cosmological and other analyses of SZ-selected clusters. We investigate the concentration-SZ relation in theory and simulations. We find that the impact of concentration on the inner SZ profile ($R \lesssim 0.75 R_{200c}$) can be captured with standard polytropic gas models. However, we find that such models do a poor job of reproducing the outer SZ profiles ($R \gtrsim 0.75 R_{200c}$) and the relation between the integrated SZ signal, $Y$, and concentration. This disagreement results from a sharp truncation of the gas pressure profile near the splashback radius, likely caused by virial shocks. We develop a simple description of the truncation that leads to a good match with simulated SZ profiles out to several $R_{200c}$ for clusters of varying mass and concentration, and that also accurately predicts the concentration-$Y$ relationship. Finally, we determine how inference of the linear bias parameter and splashback radius for SZ-selected clusters can be biased by ignoring the concentration dependence of the SZ signal.

The present study aims at characterizing the habitability conditions of exoplanets with an Earth-like magnetosphere inside the habitable zone of M stars and F stars like tau Boo, caused by the direct deposition of the stellar wind on the exoplanet surface if the magnetosphere shielding is inefficient. In addition, the radio emission generated by exoplanets with a Earth-like magnetosphere is calculated for different space weather conditions. The study is based on a set of MHD simulations performed by the code PLUTO reproducing the space weather conditions expected for exoplanets orbiting the habitable zone of M stars and F stars type tau Boo. Exoplanets hosted by M stars at 0.2 au are protected from the stellar wind during regular and CME-like space weather conditions if the star rotation period is slower than 3 days, that is to say, faster rotators generate stellar winds and interplanetary magnetic fields large enough to endanger the exoplanet habitability. Exoplanets hosted by a F stars type tau Boo at >= 2.5 au are protected during regular space weather conditions, but a stronger magnetic field compared to the Earth is mandatory if the exoplanet is close to the inner edge of the star habitable zone (2.5 au) to shield the exoplanet surface during CME-like space weather conditions. The range of radio emission values calculated in the simulations are consistent with the scaling proposed by [Zarka 2018] during regular and common CME-like space weather conditions. If the radio telescopes measure a relative low radio emission signal with small variability from an exoplanet, that may indicate favorable exoplanet habitability conditions with respect to the space weather states considered and the intrinsic magnetic field of the exoplanet.

IGR J17407-2808 is an enigmatic and poorly studied X-ray binary that was recently observed quasi-simultaneously with NuSTAR and XMM-Newton. In this paper we report the results of this observational campaign. During the first 60 ks of observation, the source was caught in a relatively low emission state, characterised by a modest variability and an average flux of ~8.3E-13 erg/cm^2/s (4-60 keV). Afterwards, IGR J17407-2808 entered a significantly more active emission state that persisted for the remaining ~40 ks of the NuSTAR observation. During this state, IGR J17407-2808 displayed several fast X-ray flares, featuring durations of ~1-100 s and profiles with either single or multiple peaks. The source flux in the flaring state reached values as high as 2E-9 erg/cm^2/s (4-60 keV), leading to a measured dynamic range during the NuSTAR and XMM-Newton campaign of >~ 10^3. We also analysed available archival photometric near-infrared data of IGR J17407-2808 to improve the constraints available so far on the the nature of the donor star hosted in this system. Our analysis shows that the donor star can be either a rare K or M-type sub-subgiant or an K type main sequence star, or sub-giant star. Our findings support the classification of IGR J17407-2808 as a low-mass X-ray binary. We discuss the source X-ray behaviour as recorded by NuSTAR and XMM-Newton in view of this revised classification.

The accretion of dark matter (DM) into astrophysical black holes slowly increases their mass. The rate of this mass accretion depends on the DM model and the model parameters. If this mass accretion effect can be measured accurately enough, it is possible to rule out some DM models, and, with the sufficient technology and the help of other DM constraints, possibly confirm one model. We propose a DM probe based on accreting pulsar-black hole binaries, which provide a high-precision measurement on binary orbital phase shifts induced by DM accretion into black holes, and can help rule out DM models and study the nature of DM.

We present a possible evolutionary pathway to form planetary nebulae (PNe) with close neutron star (NS)-white dwarf (WD) binary central stars. By employing a comprehensive binary population synthesis technique we find that the evolution involves two common envelope evolution (CEE) phases and a core collapse supernova explosion between them that forms the NS. Later the lower mass star engulfs the NS as it becomes a red giant, a process that leads to the second CEE phase and to the ejection of the envelope. This leaves a hot horizontal branch star that evolves to become a helium WD and an expanding nebula. Both the WD and the NS power the nebula. The NS in addition might power a pulsar wind nebula inside the expanding PN. From our simulations we find that the Galactic formation rate of NS-WD PNe is $1.8 \times 10^{-5} {\rm yr}^{-1}$ while the Galactic formation rate of all PNe is $0.42 {\rm yr}^{-1}$. There is a possibility that one of the observed Galactic PNe might be a NS-WD PN, and a few NS-WD PNe might exist in the Galaxy. The central binary systems might be sources for future gravitational wave detectors like LISA, and possibly of electromagnetic telescopes.

Tidal disruption events (TDEs) occur when a star passes close to a massive black hole, so that the tidal forces of the black hole exceed the binding energy of a star and cause it to be ripped apart. Part of the matter will fall onto the black hole, causing a strong increase in the luminosity. Such events are often seen in the optical or the X-ray (or both) or even at other wavelengths such as in the radio, where the diversity of observed emission is still poorly understood. The XMM-Newton catalogue of approximately a million X-ray detections covering 1283$^2$ degrees of sky contains a number of these events. Here I will show the diverse nature of a number of TDEs discovered in the catalogue and discuss their relationship with quasi periodic eruptions.

Type Ia supernovae (SNe Ia) are among the most energetic events in the Universe. They are excellent cosmological distance indicators due to the remarkable homogeneity of their light curves. However, the nature of the progenitors of SNe Ia is still not well understood. In the single-degenerate model, a carbon-oxygen white dwarf (CO WD) could grow its mass by accreting material from an asymptotic giant branch (AGB) star, leading to the formation of SNe Ia when the mass of the WD approaches to the Chandrasekhar-mass limit, known as the AGB donor channel. In this channel, previous studies mainly concentrate on the wind-accretion pathway for the mass-increase of the WDs. In the present work, we employed an integrated mass-transfer prescription for the semidetached WD+AGB systems, and evolved a number of WD+AGB systems for the formation of SNe Ia through the Roche-lobe overflow process or the wind-accretion process. We provided the initial and final parameter spaces of WD+AGB systems for producing SNe Ia. We also obtained the density distribution of circumstellar matter at the moment when the WD mass reaches the Chandrasekhar-mass limit. Moreover, we found that the massive WD+AGB sample AT 2019qyl can be covered by the final parameter space for producing SNe Ia, indicating that AT 2019qyl is a strong progenitor candidate of SNe Ia with AGB donors.

Using the data set of The Three Hundred project, i.e. 324 hydrodynamical resimulations of cluster-sized haloes and the regions of radius 15 $h^{-1}$Mpc around them, we study galaxy pairs in high-density environments. By projecting the galaxies' 3D coordinates onto a 2D plane, we apply observational techniques to find galaxy pairs. Based on a previous theoretical study on galaxy groups in the same simulations, we are able to classify the observed pairs into "true" or "false", depending on whether they are gravitationally bound or not. We find that the fraction of true pairs (purity) crucially depends on the specific thresholds used to find the pairs, ranging from around 30 to more than 80 per cent in the most restrictive case. Nevertheless, in these very restrictive cases, we see that the completeness of the sample is low, failing to find a significant number of true pairs. Therefore, we train a machine learning algorithm to help us to identify these true pairs based on the properties of the galaxies that constitute them. With the aid of the machine learning model trained with a set of properties of all the objects, we show that purity and completeness can be boosted significantly using the default observational thresholds. Furthermore, this machine learning model also reveals the properties that are most important when distinguishing true pairs, mainly the size and mass of the galaxies, their spin parameter, gas content and shape of their stellar components.

We present a detailed multi-wavelength study of star formation in the dwarf galaxy NGC 4395 which hosts an active galactic nucleus (AGN). From our observations with the Ultra-Violet Imaging Telescope, we have compiled a catalogue of 284 star forming (SF) regions, out of which we could detect 120 SF regions in H$\alpha$ observations. Across the entire galaxy, we found the extinction corrected star formation rate (SFR) in the far ultra-violet (FUV) to range from 2.0 $\times$ 10$^{-5}$ M$_\odot$yr$^{-1}$ to 1.5 $\times$ 10$^{-2}$ M$_\odot$yr$^{-1}$ with a median of 3.0 $\times$ 10$^{-4}$ M$_\odot$yr$^{-1}$ and age to lie in the range of $\sim$ 1 to 98 Myr with a median of 14 Myr. In H$\alpha$ we found the SFR to range from 7.2 $\times$ 10$^{-6}$ M$_\odot$yr$^{-1}$ to 2.7 $\times$ 10$^{-2}$ M$_\odot$yr$^{-1}$ with a median of 1.7 $\times$ 10$^{-4}$ M$_\odot$yr$^{-1}$ and age to lie between 3 to 6 Myr with a median of 5 Myr. The stellar ages derived from H$\alpha$ show a gradual decline with galactocentric distance. We found three SF regions close to the center of NGC~4395 with high SFR both from H$\alpha$ and UV which could be attributed to feedback effects from the AGN. We also found six other SF regions in one of the spiral arms having higher SFR. These are very close to supernovae remnants which could have enhanced the SFR locally. We obtained a specific SFR (SFR per unit mass) for the whole galaxy 4.64 $\times$ 10$^{-10}$ yr$^{-1}$.

We use time-domain simulations of Jupiter observations to test and develop a beam reconstruction pipeline for the Simons Observatory Small Aperture Telescopes. The method relies on a map maker that estimates and subtracts correlated atmospheric noise and a beam fitting code designed to compensate for the bias caused by the map maker. We test our reconstruction performance for four different frequency bands against various algorithmic parameters, atmospheric conditions and input beams. We additionally show the reconstruction quality as function of the number of available observations and investigate how different calibration strategies affect the beam uncertainty. For all of the cases considered, we find good agreement between the fitted results and the input beam model within a ~1.5% error for a multipole range l = 30 - 700.

Time variability of astronomical sources provides crude information on their typical size and on the implied physical mechanisms. PKS1830-211 is a remarkable radio-bright lensed quasar with a foreground molecular absorber at z=0.89. Small-scale morphological changes in the core-jet structure of the quasar -- which is magnified by the lensing -- result in a varying illumination of the absorber screen, which in turn causes variations in the absorption profile. We aim to study the time variations of the system [...] in order to obtain constraints on both the quasar activity and small-scale structures in the ISM of the absorber. We used ALMA to monitor the submm continuum emission, together with the absorption spectra of the H2O and CH molecules, with 17 visits spread over six months in 2016. [...] From the continuum data, we followed the evolution of the flux density, flux-density ratio, spectral index, and differential polarization between the two lensed images of the quasar; all quantities show significant variations related to the intrinsic activity of the quasar. We propose a simple parametric model of a core plus a ballistic plasmon to account for the continuum evolution, from which we constrain a time delay of 25+/-3~days between lensed images. The spectral lines reveal significant variations in the foreground absorption. A PCA highlights apparent wavy time variations, possibly linked to the helical jet precession period of the quasar. From the deep averaged spectra towards the SW image, we detect the absorption of 13CH and estimate an abundance ratio of 12CH/13CH~150. We also measure the oxygen isotopic ratios, 16O/18O=65.3+/-0.7 and 18O/17O=11.5+/-0.5. Finally, we find a remarkable continuous shallow trough in the water absorption spanning a velocity interval of nearly 500 km/s. This broad absorption could be the signature of an extra-planar molecular component. [Abridged]

We quantify the evolution of matter and galaxy clustering in cosmological hydrodynamical simulations via correlation and bias functions of matter and galaxies. We use simulations TNG100 and TNG300 with epochs from $z=5$ to $z=0$. We calculate spatial correlation functions of galaxies, $\xi(r)$, for simulated galaxies and dark matter (DM) particles to characterise the evolving cosmic web. We find that bias parameters decrease during the evolution, confirming earlier results. At low and medium luminosities, bias parameters of galaxies, $b_0$, are equal, suggesting that dwarf galaxies reside in the same filamentary web as brighter galaxies. Bias parameters of the lowest luminosity galaxies estimated from CFs are lower relative to CFs of particle density-limited clustered samples of DM. We find that bias parameters $b_0$, estimated from CFs of clustered DM, agree with the expected values from the fraction of particles in the clustered population, $b=1/F_c$. The cosmic web contains filamentary structures of various densities, and fractions of matter in the clustered and the unclustered populations are both less than unity. Thus the CF amplitude of the clustered matter is always higher than for all matter, i.e. bias parameter must be $b>1$. Differences between CFs of galaxies and clustered DM suggest that these functions describe different properties of the cosmic web.

We introduce a quasi-periodic restricted Hamiltonian to describe the secular motion of a small-mass planet in a multi-planetary system. In particular, we refer to the motion of $\upsilon$-And $b$ which is the innermost planet among those discovered in the extrasolar system orbiting around the $\upsilon$-Andromedae A star. We preassign the orbits of the Super-Jupiter exoplanets $\upsilon$-And $c$ and $\upsilon$-And $d$ in a stable configuration. The Fourier decompositions of their secular motions are reconstructed by using the Frequency Analysis and are injected in the equations describing the orbital dynamics of $\upsilon$-And $b$ under the gravitational effects exerted by those two external exoplanets (expected to be major ones in such an extrasolar system). We end up with a $2+3/2$ degrees of freedom Hamiltonian model; its validity is confirmed by the comparison with several numerical integrations of the complete $4$-body problem. Furthermore, the model is enriched by taking into account also the relativistic effects on the secular motion of the innermost exoplanet. We focus on the problem of the stability of $\upsilon$-And $b$ as a function of the parameters that mostly impact on its orbit, i.e. the initial values of its inclination and the longitude of its node. We study the evolution of its eccentricity, crucial to exclude orbital configurations with high probability of (quasi)collision with the central star in the long-time evolution of the system. Moreover, we also introduce a normal form approach, that further reduces our Hamiltonian model to a system with $2$ degrees of freedom, which is integrable because it admits a constant of motion related to the total angular momentum. This allows us to quickly preselect the domains of stability for $\upsilon$-And $b$, with respect to the set of the initial orbital configurations that are compatible with the observations.

We measured the average Compton profile of 461 clusters detected jointly by the South Pole Telescope (SPT) and Planck. The number of clusters included in this analysis is about one order of magnitude larger than in previous analyses. We propose an innovative method developed in Fourier space to combine optimally the Planck and SPT-SZ data, allowing us to perform a clean deconvolution of the point spread and transfer functions while simultaneously rescaling by the characteristic radial scale $R_{\rm 500}$ with respect to the critical density. The method additionally corrects for the selection bias of SPT clusters in the SPT-SZ data. We undertake a generalised Navarro-Frenk-White (NFW) fit to the profile with only one parameter fixed, allowing us to constrain the other four parameters with excellent precision. The best-fitting profile is in good agreement with the Universal Pressure Profile based on REXCESS in the inner region and with the Planck Intermediate Paper V profile based on Planck and the XMM archive in the outer region. We investigate trends with redshift and mass, finding no indication of redshift evolution but detecting a significant difference in the pressure profile of the low vs. high mass subsamples, in the sense that the low mass subsample has a profile that is more centrally-peaked than that of the high mass subsample. [abridged]

Searching for afterglows not associated with any gamma-ray bursts (GRBs) is a longstanding goal of transient surveys. These surveys provide the very chance of discovering the so-called orphan afterglows. Recently, a promising orphan afterglow candidate, AT2021any, was found by the Zwicky Transient Facility. Here we perform multi-wavelength fitting of AT2021any with three different outflow models, namely the top-hat jet model, the isotropic fireball model, and the structured Gaussian jet model. Although the three models can all fit the observed light curve well, it is found that the structured Gaussian jet model presents the best result, and thus is preferred by observations. In the framework of the Gaussian jet model, the best-fit Lorentz factor is about 68, which indicates that AT2021any should be a failed GRB. The half-opening angle of the jet and the viewing angle are found to be 0.104 and 0.02, respectively, which means that the jet is essentially observed on-axis. The trigger time of the GRB is inferred to be about 1000 s before the first detection of the orphan afterglow. An upper limit of 21.4% is derived for the radiative efficiency, which is typical in GRBs.

PSR J1933$-$6211 is a 3.5-ms pulsar in a 12.8-d orbit with a white dwarf (WD). Its high proper motion and low dispersion measure result in such significant interstellar scintillation that high signal-to-noise detections require long observing durations or fortuitous timing. We turn to the sensitive MeerKAT telescope and, combined with historic Parkes data, leverage PSR J1933$-$6211's kinematic and relativistic effects to constrain its 3D orbital geometry and the component masses. We obtain precise proper motion and parallax estimates, and measure their effects as secular changes in the Keplerian orbital parameters: a variation in orbital period of $7(1) \times 10^{-13}$ s s$^{-1}$ and a change in projected semi-major axis of $1.60(5) \times 10^{-14}$ s s$^{-1}$. A self-consistent analysis of all kinematic and relativistic effects yields a distance of $1.6^{+0.2}_{-0.3}$ kpc, an orbital inclination, $i = 55(1)$ deg and a longitude of the ascending node, $\Omega = 255^{+8}_{-14}$ deg. The probability densities for $\Omega$ and $i$ and their symmetric counterparts, ($180-i$, $360-\Omega$), are seen to depend on the fiducial orbit used to measure the time of periastron passage. We investigate this unexpected dependence and rule out software-related causes using simulations. Nevertheless, we constrain the pulsar and WD masses to $1.4^{+0.3}_{-0.2}$ M$_\odot$ and $0.43(5)$ M$_\odot$ respectively. These strongly disfavour a helium-dominated WD. The orbital similarities between PSRs J1933$-$6211 and J1614$-$2230 suggest they underwent Case A Roche lobe overflow, an extended evolution while the companion star is still on the Main Sequence. However, with a mass of $\sim 1.4$ M$_\odot$, PSR J1933$-$6211 has not accreted significant matter. This highlights the low accretion efficiency of the spin-up process and suggests that observed neutron star masses are mostly a result of supernova physics.

Gamma-ray bursts categorically produce broadband afterglow emission, but in some cases, emission in the optical band is dimmer than expected based on the contemporaneously observed X-ray flux. This phenomenon, aptly dubbed "optical darkness", has been studied extensively in long GRBs (associated with the explosive deaths of massive stars), with possible explanations ranging from host environment extinction to high redshift to possibly unique emission mechanisms. However, investigations into optical darkness in short GRBs (associated with the mergers of compact object binaries) have thus far been limited. This work implements a procedure for determining the darkness of GRBs based on spectral indices calculated using temporally-matched Swift-XRT data and optical follow-up observations; presents a complete and up-to-date catalog of known short GRBs that exhibit optical darkness; and outlines some of the possible explanations for optically dark short GRBs. In the process of this analysis, we developed versatile and scalable data processing code that facilitates reproducibility and reuse of our pipeline. These analysis tools and resulting complete sample of dark short GRBs enable a systematic statistical study of the phenomenon and its origins, and reveal that optical darkness is indeed quite rare in short GRBs, and highly dependent on observing response time and observational effects.

We report the detection of CO(6-5) and CO(7-6) and their underlying continua from the host galaxy of quasar J100758.264+211529.207 (P\=oniu\=a'ena) at z=7.5419, obtained with the NOrthern Extended Millimeter Array (NOEMA). P\=oniu\=a'ena belongs to the HYPerluminous quasars at the Epoch of ReionizatION (HYPERION) sample of 17 $z>6$ quasars selected to be powered by supermassive black holes (SMBH) which experienced the fastest mass growth in the first Gyr of the Universe. The one reported here is the highest-redshift measurement of the cold and dense molecular gas to date. The host galaxy is unresolved and the line luminosity implies a molecular reservoir of $\rm M(H_2)=(2.2\pm0.2)\times 10^{10}$ $\rm M_\odot$, assuming a CO spectral line energy distribution typical of high-redshift quasars and a conversion factor $\alpha=0.8$ $\rm M_{\odot} (K\,km \, s^{-1} \,pc^{2})^{-1} $. We model the cold dust spectral energy distribution (SED) to derive a dust mass of M$_{\rm dust} =(2.1\pm 0.7)\times 10^8$ $\rm M_\odot$, and thus a gas to dust ratio $\sim100$. Both the gas and dust mass are not dissimilar from the reservoir found for luminous quasars at $z\sim6$. We use the CO detection to derive an estimate of the cosmic mass density of $\rm H_2$, $\Omega_{H_2} \simeq 1.31 \times 10^{-5}$. This value is in line with the general trend suggested by literature estimates at $ z < 7 $ and agrees fairly well with the latest theoretical expectations of non-equilibrium molecular-chemistry cosmological simulations of cold gas at early times.

We review the recent theoretical and observational developments concerning the interaction of the Large Magellanic Cloud (LMC) with the Milky Way and its neighbourhood. An emerging picture is that the LMC is a fairly massive companion (10-20% of the Milky Way mass) and just passed the pericentre of its orbit, likely for the first time. The gravitational perturbation caused by the LMC is manifested at different levels. The most immediate effect is the deflection of orbits of stars, stellar streams or satellite galaxies passing in the vicinity of the LMC. Less well known but equally important is the displacement (reflex motion) of central regions of the Milky Way about the centre of mass of both galaxies. Since the Milky Way is not a rigid body, this displacement varies with the distance from the LMC, and as a result, the Galaxy is deformed and its outer regions (beyond a few tens kpc) acquire a net velocity with respect to its centre. These phenomena need to be taken into account at the level of precision warranted by current and future observational data, and improvements on the modelling side are also necessary for an adequate interpretation of these data.

The high-precision astrometric mission GAIA recently reported the remarkable discovery of a Sun-like star closely orbiting a dark object, with a semi-major axis and period of $1.4\, \rm{AU}$ and $187.8$ days respectively. While the plausible expectation for the central dark object is a black hole, the evolutionary mechanism leading to the formation of such a two-body system is highly challenging. Here, we challenge the scenario of a central black hole and show that the observed orbital dynamics can be explained under fairly general assumptions if the central dark object is a stable clump of bosonic particles of spin-0, or spin-1, known as a boson star. We further explain how future astrometric measurements of similar systems will provide an exciting opportunity to probe the fundamental nature of compact objects and test compact alternatives to black holes.

[Abridged] Galaxy clusters are the most massive gravitationally-bound systems in the universe and are widely considered to be an effective cosmological probe. We propose the first Machine Learning method using galaxy cluster properties to derive unbiased constraints on a set of cosmological parameters, including Omega_m, sigma_8, Omega_b, and h_0. We train the machine learning model with mock catalogs including "measured" quantities from Magneticum multi-cosmology hydrodynamical simulations, like gas mass, gas bolometric luminosity, gas temperature, stellar mass, cluster radius, total mass, velocity dispersion, and redshift, and correctly predict all parameters with uncertainties of the order of ~14% for Omega_m, ~8% for sigma_8, ~6% for Omega_b, and ~3% for h_0. This first test is exceptionally promising, as it shows that machine learning can efficiently map the correlations in the multi-dimensional space of the observed quantities to the cosmological parameter space and narrow down the probability that a given sample belongs to a given cosmological parameter combination. In the future, these ML tools can be applied to cluster samples with multi-wavelength observations from surveys like CSST in the optical band, Euclid and Roman in the near-infrared band, and eROSITA in the X-ray band to constrain both the cosmology and the effect of the baryonic feedback.

In-situ plasma measurements as well as remote mapping of energetic neutral atoms around Jupiter provide indirect evidence that an enhancement of neutral gas is present near the orbit of the moon Europa. Simulations suggest that such a neutral gas torus can be sustained by escape from Europa's atmosphere and consists primarily of molecular hydrogen, but the neutral gas torus has not yet been measured directly through emissions or in-situ. Here we present observations by the Cosmic Origins Spectrograph of the Hubble Space Telescope (HST/COS) from 2020 and 2021, which scanned the equatorial plane between 8 and 10 planetary radii west of Jupiter. No neutral gas emissions are detected. We derive upper limits on the emissions and compare these to modelled emissions from electron impact and resonant scattering using a Europa torus Monte Carlo model for the neutral gases. The comparison supports the previous findings that the torus is dilute and primarily consists of molecular hydrogen. A detection of sulfur ion emissions radially inward of the Europa orbit is consistent with emissions from the extended Io torus and with sulfur ion fractional abundances as previously detected.

We propose a cross-internal linear combination (cross-ILC) approach to measure the small-scale cosmic microwave background (CMB) anisotropies robustly against the contamination from astrophysical signals. In particular, we focus on the mitigation of systematics from cosmic infrared background (CIB) and thermal Sunyaev-Zeldovich (tSZ) signals in kinematic SZ (kSZ) power spectrum and CMB lensing. We show the cross-spectrum measurement between two CMB maps created by nulling the contributions from CIB (CIB-free map) and tSZ (tSZ-free map) to be robust for kSZ as the approach significantly suppresses the total contribution of CIB and tSZ signals. Similarly, for CMB lensing, we use the approach introduced by Madhavacheril & Hill (2018) but with a slight modification by using the tSZ-free and CIB-free maps in the two legs of the quadratic estimator. By cross-correlating the CMB lensing map created using this technique with galaxy surveys, we show that the biases from both CIB/tSZ are negligible. We also compute the impact of unmodeled CIB/tSZ residuals on kSZ and cosmological parameters finding that the kSZ measured using the standard ILC to be significantly biased. The kSZ estimate from the cross-ILC remains less affected by CIB/tSZ making it crucial for CMB surveys such as the South Pole Telescope (SPT), Simons Observatory (SO) and CMB-S4. With the cross-ILC method, we find the total kSZ power spectrum can be measured at very high significance: $35\sigma$ by SPT, $22\sigma$ by SO, and $80\sigma$ by CMB-S4. We forecast constraints on the epoch of reionization using the kSZ power spectrum and find that the duration of reionization, currently unconstrained by {\it Planck}, can be constrained to $\sigma(z_{\rm dur})$= 1.5 (or) 0.5 depending on the choice of $\tau_{\rm re}$ prior. The data products and codes can be downloaded from this https://github.com/sriniraghunathan/cross_ilc_methods_paper.

We conduct the first study of how the relative quenching probability of galaxies depends on environment over the redshift range $0.5 < z < 3$, using data from the UKIDSS Ultra-Deep Survey. By constructing the stellar mass functions for quiescent and post-starburst (PSB) galaxies in high, medium and low density environments to $z = 3$, we find an excess of quenched galaxies in dense environments out to at least $z \sim 2$. Using the growth rate in the number of quenched galaxies, combined with the star-forming galaxy mass function, we calculate the probability that a given star-forming galaxy is quenched per unit time. We find a significantly higher quenching rate in dense environments (at a given stellar mass) at all redshifts. Massive galaxies (M$_* > 10^{10.7}$ M$_{\odot}$) are on average 1.7 $\pm$ 0.2 times more likely to quench per Gyr in the densest third of environments compared to the sparsest third. Finally, we compare the quiescent galaxy growth rate to the rate at which galaxies pass through a PSB phase. Assuming a visibility timescale of 500 Myr, we find that the PSB route can explain $\sim$ 50\% of the growth in the quiescent population at high stellar mass (M$_* > 10^{10.7}$ M$_{\odot}$) in the redshift range $0.5 < z < 3$, and potentially all of the growth at lower stellar masses.

The universe's current acceleration is a pretty recent phenomenon in cosmological time scales. This means that the modes that have left our horizon since the beginning of the contemporary acceleration phase, have not really reached the exact IR limit. Noting this observation, we reconsider the possibility of having a ghost condensate as dark energy with a sixth-order dispersion relation. Looking at the three-point function of such a theory, we obtain the constraints on the coefficient of the sixth-order dispersion relation to avoid strong coupling. Such a ghost condensate if coupled to the standard model fields, induces a constant Lorentz-violating spin-dependent force, which can gravitate or anti-gravitate.

To numerically evolve the full Einstein equations (or modifications thereof), simulations of cosmological spacetimes must rely on a particular formulation of the field equations combined with a specific gauge/frame choice. Yet truly physical results cannot depend on the given formulation or gauge/frame choice. In this paper, we present a resolution of the gauge problem and, as an example, numerically implement it to evaluate our previous work on contracting spacetimes.

Recently, the IceCube collaboration observed a neutrino excess in the direction of NGC 1068 with high statistical significance. This constitutes the second detection of an astrophysical neutrino point source after the discovery of a variable emission originating from the blazar TXS~0506+056. Neutrinos emitted by these sources traverse huge, well-determined distances on their way to Earth. This makes them a promising tool to test new physics in the neutrino sector. We consider secret interactions with the cosmic neutrino background and discuss their impact on the flux of neutrino point sources. The observation of emission from NGC 1068 and TXS 0506+056 can then be used to put limits on the strength of the interaction. We find that our ignorance of the absolute neutrino masses has a strong impact and, therefore, we present limits in two benchmark scenarios with the sum of the neutrino masses around their lower and upper limits.

We present global predictions of the ground state mass of atomic nuclei based on a novel Machine Learning (ML) algorithm. We combine precision nuclear experimental measurements together with theoretical predictions of unmeasured nuclei. This hybrid data set is used to train a probabilistic neural network. In addition to training on this data, a physics-based loss function is employed to help refine the solutions. The resultant Bayesian averaged predictions have excellent performance compared to the testing set and come with well-quantified uncertainties which are critical for contemporary scientific applications. We assess extrapolations of the model's predictions and estimate the growth of uncertainties in the region far from measurements.

Generalized disformal transformations enable us to construct the generalized disformal Horndeski theories, which form the most general class of ghost-free scalar-tensor theories to this date. We extend the effective field theory (EFT) of cosmological perturbations to incorporate these generalized disformal Horndeski theories. The main difference from the conventional EFT is that our extended EFT involves operators with higher spatial derivatives of the lapse function. Our EFT also accommodates the generalized disformal transformation of U-DHOST theories.

Finsler geometry is a natural and fundamental generalization of Riemann geometry, and is a tool to research Lorentz invariance violation. We find the connection between the most general modified dispersion relation and a pseudo-Finsler structure, and then we calculate the arrival time delay of astroparticles with different modified dispersion relations in the framework of Finsler geometry. The result suggests that the time delay is irrelevant with the exact form of the modified dispersion relation. If the modified term becomes 0 when $E=p$, there is no arrival time difference, otherwise the time delays only depend on the Lorentz violation scale and the order at which the Lorentz invariance breaks.

We study the properties of hot and dense neutron star matter under the presence of strong magnetic fields using two Skyrme interactions, namely the LNS and the BSk21 ones. Asking for $\beta$--stability and charge neutrality, we construct the equation of state of the system and analyze its composition for a range of densities, temperatures and magnetic field intensities of interest for the study of supernova and proto-neutron star matter, with a particular interest on the degree of spin-polarization of the different components. The results show that system configurations with larger fractions of spin up protons and spin down neutrons and electrons are energetically favored over those with larger fractions of spin down protons and spin up neutrons and electrons. The effective mass of neutrons and protons is found to be in general larger for the more abundant of their spin projection component, respectively, spin down neutrons and spin up protons. The effect of the magnetic field on the Helmhotz total free energy density, pressure and isothermal compressibility of the system is almost negligible for all the values of the magnetic field considered.

Cosmological implication of a generalized model of two scalar and two vector fields, in which both scalar fields are non-minimally coupled to each vector field, is studied in this paper. In particular, we will seek an anisotropic power-law inflationary solution to this model. Furthermore, the stability of the obtained solution will be examined by using the dynamical system approach. As a result, we will show that this solution turns out to be stable and attractive during the inflationary phase as expected due to the existence of the unusual couplings between two scalar and two vector fields. Remarkably, we will point out that the existence of phantom field will lead to an instability of the corresponding anisotropic power-law inflation.

We study the possibility of constraining a scenario with high scale first order phase transition (FOPT) responsible for the cogenesis of baryon and dark matter using gravitational wave (GW) (non)-observations. While the FOPT at high scale is responsible for generating baryon asymmetry through leptogenesis and dark matter via the \textit{mass-gain} mechanism, the resulting GW spectrum falls within the ongoing LIGO-VIRGO experimental sensitivity. The dark matter is preferred to be in the non-thermal ballpark with sub-GeV masses and the criteria of successful dark matter relic rules out a large portion of the parameter space consistent with high scale FOPT and successful leptogenesis. Some part of the parameter space allowed from dark matter and leptogenesis criteria also gives rise to a large signal-to-noise ratio at ongoing experiments and hence can be disfavoured in a conservative way from the non-observation of such stochastic GW background. Future data from ongoing and planned experiments will offer a complementary and indirect probe of the remaining parameter space which is typically outside the reach of any direct experimental probe.

Any violation of the fundamental principles of general relativity (GR), including the violations of the equivalence principle and parity/Lorentz symmetries, could induce possible derivations in the gravitational wave (GW) propagations so they can be tested/constrained directly by the GW data. In this letter, we present a universal parametrization for characterizing possible derivations from GW propagations in GR. This parametrization provides a general framework for exploring possible modified GW propagations arising from a large number of modified theories of gravity. With this parameterization, we construct the modified GW waveforms generated by the coalescence of compact binaries with the effects of the gravitational parity/Lorentz violations, then analyze the open data of compact binary merging events detected by LIGO/Virgo/KAGRA collaboration. We do not find any signatures of gravitational parity/Lorentz violations, thereby allowing us to place several of the most stringent constraints on parity/Lorentz violations in gravity and a first constraint on the Lorentz-violating damping effect in GW. This also represents the most comprehensive tests on the modified GW propagations.

Particles measured in large gradual solar energetic particle (SEP) events are believed to be predominantly accelerated at shocks driven by coronal mass ejections (CMEs). Ion charge state and composition analyses suggest that the origin of the seed particle population for the mechanisms of particle acceleration at CME-driven shocks is not the bulk solar wind thermal material, but rather a suprathermal population present in the solar wind. This suprathermal population could result from remnant material accelerated in prior solar flares and/or preceding CME-driven shocks. In this work, we examine the distribution of this suprathermal particle population in the inner heliosphere by combining a magnetohydrodynamic (MHD) simulation of the solar wind and a Monte-Carlo simulation of particle acceleration and transport. Assuming that the seed particles are uniformly distributed near the Sun by solar flares of various magnitudes, we study the longitudinal distribution of the seed population at multiple heliocentric distances. We consider a non-uniform background solar wind, consisting of fast and slow streams that lead to compression and rarefaction regions within the solar wind. Our simulations show that the seed population at a particular location (e.g., 1 au) is strongly modulated by the underlying solar wind configuration. Corotating interaction regions (CIRs) and merged interactions regions (MIRs) can strongly alter the energy spectra of the seed particle populations. In addition, cross-field diffusion plays an important role in mitigating strong variations of the seed population in both space and energy.