Papers for Monday, Dec 06 2021

We present a multiwavelength study of the H II region Sh 2-301 (S301) using deep optical data, near-infrared data, radio continuum data and other archival data at longer wavelengths. A cluster of young stellar objects (YSOs) is identified in the north-east (NE) direction of S301. The H{\alpha} and radio continuum images trace the distribution of the ionized gas surrounding a massive star ALS 207, and the S301 H II region is bounded by an arc-like structure of gas and dust emission in the south-eastern direction. The north-western part of S301 seems to be devoid of gas and dust emission, while the presence of molecular material between the NE cluster and the central massive star ALS 207 is found. The distribution of warm dust emission, ionized gas, and neutral hydrogen together suggests a blistered morphology of the S301 H II region powered by ALS 207, which appears to be located near the edge of the cloud. The location of the NE cluster embedded in the cold molecular cloud is found opposite to the blistered morphology. There is a noticeable age difference investigated between the massive star and the NE cluster. This age difference, pressure calculation, photodissociation regions (PDRs), and the distribution of YSOs favour the positive feedback of the massive star ALS 207 in S301. On a wider scale of S301, the H II region and the young stellar cluster are depicted toward the central region of a hub-filamentary system, which is evident in the infrared images.

The unusually low velocity dispersion and large size of the Crater II dwarf spheroidal pose a challenge to our understanding of dwarf galaxies in the Lambda Cold Dark Matter (LCDM) cosmogony. The low velocity dispersion suggests either a dark halo mass much lower than the minimum expected from hydrogen cooling limit arguments, or one that is in the late stages of extreme tidal stripping. The tidal interpretation has been favoured in recent work, and is supported by the small pericentric distances allowed by orbits consistent with available estimates of distance, proper motions, and radial velocity. We use N-body simulations to examine this interpretation in some detail, assuming a Navarro-Frenk-White (NFW) profile for Crater II's progenitor halo. Our main finding is that, although the low velocity dispersion can indeed result from the effect of tides, the large size of Crater II is inconsistent with this hypothesis. This is because galaxies stripped to match the observed velocity dispersion are also reduced to sizes much smaller than the observed half-light radius of Crater II. Unless the size of Crater II has been substantially overestimated, reconciling this system with LCDM requires that either (i) it is not bound and near equilibrium (unlikely, given its crossing time is shorter than the time elapsed since pericentre), or that (ii) its progenitor halo deviates from the assumed NFW profile. The latter alternative may signal that baryons can affect the inner halo cusp even in extremely faint dwarfs or, more intriguingly, may signal effects associated with the intimate nature of the dark matter, such as finite self-interactions, or other such deviations from the canonical LCDM paradigm.

Large-scale astronomical surveys have the potential to capture data on large numbers of strongly gravitationally lensed supernovae (LSNe). To facilitate timely analysis and spectroscopic follow-up before the supernova fades, an LSN needs to be identified soon after it begins. To quickly identify LSNe in optical survey datasets, we designed ZipperNet, a multi-branch deep neural network that combines convolutional layers (traditionally used for images) with long short-term memory (LSTM) layers (traditionally used for time series). We tested ZipperNet on the task of classifying objects from four categories -- no lens, galaxy-galaxy lens, lensed type Ia supernova, lensed core-collapse supernova -- within high-fidelity simulations of three cosmic survey data sets -- the Dark Energy Survey (DES), Rubin Observatory's Legacy Survey of Space and Time (LSST), and a Dark Energy Spectroscopic Instrument (DESI) imaging survey. Among our results, we find that for the LSST-like dataset, ZipperNet classifies LSNe with a receiver operating characteristic area under the curve of 0.97, predicts the spectroscopic type of the lensed supernovae with 79\% accuracy, and demonstrates similarly high performance for LSNe 1-2 epochs after first detection. We anticipate that a model like ZipperNet, which simultaneously incorporates spatial and temporal information, can play a significant role in the rapid identification of lensed transient systems in cosmic survey experiments.

We present "Extending the Satellites Around Galactic Analogs Survey" (xSAGA), a method for identifying low-$z$ galaxies on the basis of optical imaging, and results on the spatial distributions of xSAGA satellites around host galaxies. Using spectroscopic redshift catalogs from the SAGA Survey as a training data set, we have optimized a convolutional neural network (CNN) to identify $z < 0.03$ galaxies from more distant objects using image cutouts from the DESI Legacy Imaging Surveys. From the sample of $> 100,000$ CNN-selected low-$z$ galaxies, we identify $>20,000$ probable satellites located between 36-300 projected kpc from NASA-Sloan Atlas central galaxies in the stellar mass range $9.5 < \log(M_\star/M_\odot) < 11$. We characterize the incompleteness and contamination for CNN-selected samples, and apply corrections in order to estimate the true number of satellites as a function of projected radial distance from their hosts. Satellite richness depends strongly on host stellar mass, such that more massive host galaxies have more satellites, and on host morphology, such that elliptical hosts have more satellites than disky hosts with comparable stellar masses. We also find a strong inverse correlation between satellite richness and the magnitude gap between a host and its brightest satellite. The normalized satellite radial distribution between 36-300 kpc does not depend strongly on host stellar mass, morphology, or magnitude gap. The satellite abundances and radial distributions we measure are in reasonable agreement with predictions from hydrodynamic simulations. Our results deliver unprecedented statistical power for studying satellite galaxy populations, and highlight the promise of using machine learning for extending galaxy samples of wide-area surveys.

Stellar-mass BHs (sBHs) are predicted to be embedded in active galactic nuclei (AGN) disks due to gravitational drag and in-situ star formation. However, we find that due to a high gas density in an AGN disk environment, compact objects may rapidly grow to intermediate-mass BHs and deplete matter from the AGN disk unless accretion is suppressed by some feedback process(es). These consequences are inconsistent with AGN observations and the dynamics of the Galactic center. Here we consider mechanical feedback mechanisms for the reduction of gas accretion. Rapidly accreting sBHs launch winds and/or jets via the Blandford-Znajek mechanism, which produce high-pressure shocks and cocoons. Such a shock and cocoon can spread laterally in the plane of the disk, eject the outer regions of a circum-sBH disk (CsBD) and puncture a hole in the AGN disk with horizontal size comparable to the disk scale-height. Since the depletion timescale of the bound CsBD is much shorter than the resupply timescale of gas to the sBH, the time-averaged accretion rate onto sBHs is reduced by this process by a factor of $\sim 10$--$100$. This feedback mechanism can therefore help alleviate the sBH over-growth and AGN-disk depletion problems. On the other hand, we find that cocoons of jets can unbind a large fraction of the gas accreting in the disks of less massive SMBHs, which may help explain the dearth of high-Eddington ratio AGNs with SMBH mass $\lesssim10^5{\rm M_\odot}$.

Recent works have shown that weak lensing magnification must be included in upcoming large-scale structure analyses, such as for the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST), to avoid biasing the cosmological results. In this work we investigate whether including magnification has a positive impact on the precision of the cosmological constraints, as well as being necessary to avoid bias. We forecast this using an LSST mock catalog, a halo model to calculate the galaxy power spectra and a multiplicative factor to account for systematic effects in the magnification measurement. We find that including magnification has little effect on the cosmological parameter constraints for an LSST galaxy clustering analysis. In particular, we find that for the LSST gold sample ($i < 25.3$) including weak lensing magnification only improves the galaxy clustering constraint on $\Omega_{\rm{m}}$ by a factor of 1.03, and when using a deeper LSST mock sample ($i<26.5$) by a factor of 1.3. Since magnification predominantly contributes to the clustering measurement and provides similar information to that of cosmic shear, this mild improvement would be further reduced for a combined galaxy clustering and shear analysis. We also confirm that not modelling weak lensing magnification will catastrophically bias the cosmological results from LSST. Magnification must therefore be included in LSST large-scale structure analyses even though it does not significantly enhance the cosmological constraints.

Murchikova et al 2019 discovered a disk of cool ionized gas within 20,000 Schwarzschild radii of the Milky Way's Galactic Center black hole Sagittarius A*. They further demonstrated that the ionizing photon flux in the region is enough to keep the disk ionized, but there is not ample excess of this radiation. This raised the possibility that some neutral gas could also be in the region shielded within the cool ionized clumps. Here we present ALMA observations of a broad 1.3 millimeter hydrogen recombination line H30alpha: n = 31 -> 30, conducted during the flyby of the S0-2 star by Sgr A*. We report that the velocity-integrated H30alpha line flux two month prior to the S0-2 pericenter passage is about 20% larger than it was one month prior to the passage. The S0-2 is a strong source of ionizing radiation moving at several thousand kilometers per second during the approach. Such a source is capable of ionising parcels of neural gas along its trajectory, resulting in variation of the recombination line spectra from epoch to epoch. We conclude that there are at least (6.6 +- 3.3) x 10^{-6} Msun of neutral gas within 20,000 Schwarzschild radii of Sgr A*.

We apply a conventional accretion disk model to the FU Ori-type objects HBC 722 and Gaia 17bpi. Our base model is a steady-state, thin Keplerian disk featuring a modified Shakura-Sunyaev temperature profile, with each annulus radiating as an area-weighted spectrum given by a NextGen atmosphere at the appropriate temperature. We explore departures from the standard model by altering the temperature distribution in the innermost region of the disk to account for "boundary region"-like effects. We consider the overall spectral energy distribution (SED) as well as medium- and high-resolution spectra in evaluating best-fit models to the data. Parameter degeneracies are studied via a Markov Chain Monte Carlo (MCMC) parameter estimation technique. Allowing all parameters to vary, we find accretion rates for HBC 722 of $\dot{M} = 10^{-4.90} M_\odot \textrm{ yr}^{-1}\; {}^{+0.99}_{-0.40} \textrm{ dex}$ and for Gaia 17bpi of $\dot{M} = 10^{-6.70} M_\odot \textrm{ yr}^{-1}\; {}^{+0.46}_{-0.36} \textrm{ dex}$; the corresponding maximum disk temperatures are $7100_{-500}^{+300}$ K and $7900_{-400}^{+900}$ K, respectively. While the accretion rate of HBC 722 is on the same order as other FU Ori-type objects, Gaia 17bpi has a lower rate than previously reported as typical, commensurate with its lower luminosity. Alternate models that fix some disk or stellar parameters are also presented, with tighter confidence intervals on the remaining fitted parameters. In order to improve upon the somewhat large credible intervals for the $\dot{M}$ values, and make progress on boundary layer characterization, flux-calibrated ultraviolet spectroscopy is needed.

Context. The S\'ersic profile is a widely-used model to describe the surface brightness distribution of galaxies. Spiral galaxies, however, are qualitatively different from a S\'ersic model. Aims. The goal of this study is to assess how accurately can the total flux and half-light radius of a galaxy with spiral arms be recovered, when fitted with a S\'ersic profile. Methods. I selected a sample of bulge-dominated galaxies with spiral arms. Using photometric data from the Hyper Suprime-Cam survey, I estimated the contribution of the spiral arms to their total flux. Then I generated simulated images of galaxies with similar characteristics, fitted them with a S\'ersic model and quantified the error on the determination of the total flux and half-light radius. Results. Spiral arms can introduce biases on the photometry of galaxies, in a way that depends on the underlying smooth surface brightness profile, the location of the arms, and the depth of the photometric data. A set of spiral arms accounting for 10% of the flux of a bulge-dominated galaxy typically causes the total flux and the half-light radius to be overestimated by 15% and 30%, respectively. This bias, however, is much smaller if the galaxy is disk-dominated. Conclusions. Galaxies with a prominent bulge and a non-zero contribution from spiral arms are the most susceptible to biases in the total flux and half-light radius, when fitted with a S\'ersic profile. If photometric measurements with high accuracy are required, then measurements over finite apertures are to be preferred over global estimates of the flux.

We connect galaxy properties with their orbital classification by analysing a sample of galaxies with stellar mass $M_{\star} \geq 10^{8.5}h^{-1}M_\odot$ residing in and around massive and isolated galaxy clusters with mass $M_{200} > 10^{15}h^{-1}M_\odot$ at redshift $z=0$. The galaxy population is generated by applying the semi-analytic model of galaxy formation SAG on the cosmological simulation MultiDark Planck 2. We classify galaxies considering their real orbits (3D) and their projected phase-space position using the ROGER code (2D). We define five categories: cluster galaxies, galaxies that have recently fallen into a cluster, backsplash galaxies, infalling galaxies, and interloper galaxies. For each class, we analyse the $g-r$ colour, the specific star formation rate (sSFR), and the stellar age, as a function of the stellar mass. For the 3D classes, we find that cluster galaxies have the lowest sSFR, and are the reddest and the oldest, as expected from environmental effects. Backsplash galaxies have properties intermediate between the cluster and recent infaller galaxies. For each 2D class, we find an important contamination by other classes. We find it necessary to separate the galaxy populations in red and blue to perform a more realistic analysis of the 2D data. For the red population, the 2D results are in good agreement with the 3D predictions. Nevertheless, when the blue population is considered, the 2D analysis only provides reliable results for recent infallers, infalling galaxies and interloper galaxies.

We reexamine the relationship between energy spectral indices and element abundance enhancements in solar energetic particle (SEP) events at energies of a few MeV/amu. We find correlated behavior only in the largest gradual SEP4 events when all ions are accelerated from the ambient coronal plasma by shock waves driven by fast, wide coronal mass ejections (CMEs). This correlated abundance behavior can track complex time variations in the spectral indices during an event. In other (SEP3) events, CME-driven shock waves, days apart, sample seed particles from a single pool of suprathermal impulsive ions contributed earlier. Of the smaller, Fe-rich, impulsive SEP events, previously related to magnetic reconnection in solar jets, over half are subsequently reaccelerated by CME-driven shock waves (SEP2) causing typical ion intensities to have a 64% correlation with shock speed. In these SEP2 events, onset of shock acceleration is signaled by a new component in the abundances, large proton excesses. The remaining SEP1 events lack evidence of shock acceleration. However, for all these events (SEP1 - SEP3) with abundances determined by magnetic reconnection, spectra and abundances are decoupled.

Observations of YZ LMi show enhanced emission along the stream trajectory beyond impact at disk rim during outbursts as well as when the quiescent disk is large. We investigated whether these features can be explained in terms of either gas stream overflow or penetration within the frameworks of the disk-instability (DIM) and the mass-transfer instability (MTIM) models of outbursting disks. Gas stream overflow is not possible because the vertical scaleheight of the stream is significantly lower than that of the outer disk and because there is no combination of parameters which enables stream overflow on a larger disk while preventing it on a smaller disk. Stream penetration requires the gas stream to be denser than the outer disk regions. This requirement cannot be met by a low-viscosity DIM disk because its density is significantly larger than that of the gas stream over the whole range of mass transfer rates where the thermal-viscous instability occurs. On the other hand, the high-viscosity MTIM disk has much lower densities which decrease with increasing radius, easily allowing for gas stream penetration during outbursts (when mass transfer rate and stream density increase) as well as in large quiescent disks. The observed features are not consistent with DIM, but can be plausibly explained by MTIM. These results suggest that the outbursts of YZ LMi are the response of a high-viscosity disk to bursts of enhanced mass transfer rate. In this case, the outburst decline timescale of (2-3) d implies a viscosity parameter in the range alpha=3-4.

We present the photometric analysis of Gaia19bld, a high-magnification ($A\approx60$) microlensing event located in the southern Galactic plane, which exhibited finite source and microlensing parallax effects. Due to a prompt detection by the Gaia satellite and the very high brightness of $I = 9.05~$mag at the peak, it was possible to collect a complete and unique set of multi-channel follow-up observations, which allowed us to determine all parameters vital for the characterisation of the lens and the source in the microlensing event. Gaia19bld was discovered by the Gaia satellite and was subsequently intensively followed up with a network of ground-based observatories and the Spitzer Space Telescope. We collected multiple high-resolution spectra with Very Large Telescope (VLT)/X-Shooter to characterise the source star. The event was also observed with VLT Interferometer (VLTI)/PIONIER during the peak. Here we focus on the photometric observations and model the light curve composed of data from Gaia, Spitzer, and multiple optical, ground-based observatories. We find the best-fitting solution with parallax and finite source effects. We derived the limit on the luminosity of the lens based on the blended light model and spectroscopic distance. We compute the mass of the lens to be $1.13 \pm 0.03~M_{\odot}$ and derive its distance to be $5.52^{+0.35}_{-0.64}~\mathrm{kpc}$. The lens is likely a main sequence star, however its true nature has yet to be verified by future high-resolution observations. Our results are consistent with interferometric measurements of the angular Einstein radius, emphasising that interferometry can be a new channel for determining the masses of objects that would otherwise remain undetectable, including stellar-mass black holes.

In near-infrared bands, co-adding and tiling of astronomical imaging datasets require a sufficiently high calibration quality (flat fielding, background subtraction). Here we present a complete workflow for obtaining imaging mosaics with the MMT and Magellan Infrared Spectrograph (MMIRS) operated at the 6.5-m MMT in Arizona and open-source add-on tools developed for the MMIRS pipeline for preparation and data reduction of mosaic observations. We describe pre-observing actions, such as design of dithering patterns and mosaic layouts and post-processing steps to perform absolute astrometric and photometric calibration, and also generate HiPS maps to display the final data product in Aladin / Aladin Lite.

This paper presents time series observations and analysis of broadband night sky airglow intensity 4 September 2018 through 30 April 2020. Data were obtained at 5 sites spanning more than 8500 km during the historically deep minimum of Solar Cycle 24 into the beginning of Solar Cycle 25. New time series observations indicate previously unrecognized significant sources of broadband night sky brightness variations, not involving corresponding changes in the Sun's 10.7cm solar flux, occur during deep solar minimum. Even during a deep solar minimum the natural night sky is rarely, if ever, constant in brightness. Changes with time scales of minutes, hours, days, and months are observed. Semiannual night sky brightness variations are coincident with changes in the orientation of Earth's magnetic field relative to the interplanetary magnetic field. Solar wind plasma streams from solar coronal holes arriving at Earth's bow shock nose are coincident with major night sky brightness increase events. Sites more than 8500 km along the Earth's surface experience nights in common with either very bright or very faint night sky airglow emissions. The reason for this observational fact remains an open question. It is plausible, terrestrial night airglow and geomagnetic indices have similar responses to the solar energy input into Earth's magnetosphere. Our empirical results contribute to a quantitative basis for understanding and predicting broadband night sky brightness variations. They are applicable in astronomical, planetary science, space weather, light pollution, biological, and recreational studies.

Most black holes possess accretion disks. Models of such disks inform observations and constrain the properties of the black holes and their surrounding medium. Here, we study isothermal shocks in a thin black hole accretion flow. Modelling infinitesimal molecular viscosity allows the use of multiple-scales matched asymptotic methods. We thus derive the first explicit calculations of isothermal shock stability. We find that the inner shock is always unstable, and the outer shock is always stable. The growth/decay rates of perturbations depend only on an effective potential and the incoming--outgoing flow difference at the shock location. We give a prescription of accretion regimes in terms of angular momentum and black hole radius. Accounting for angular momentum dissipation implies unstable outer shocks in much of parameter space, even for realistic viscous Reynolds numbers of the order $\approx 10^{20}$.

The MaNGA Stellar Library (MaStar) is a large collection of high-quality empirical stellar spectra designed to cover all spectral types and ideal for use in the stellar population analysis of galaxies observed in the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey. The library contains 59,266 spectra of 24,130 unique stars with spectral resolution $R\sim1800$ and covering a wavelength range of $3,622-10,354$ \r{A}. In this work, we derive five physical parameters for each spectrum in the library: effective temperature ($T_{eff}$), surface gravity ($\log g$), metalicity ($[Fe/H]$), micro-turbulent velocity ($\log(v_{micro})$), and alpha-element abundance ($[\alpha/Fe]$). These parameters are derived with a flexible data-driven algorithm that uses a neural network model. We train a neural network using the subset of 1,675 MaStar targets that have also been observed in the Apache Point Observatory Galactic Evolution Experiment (APOGEE), adopting the independently-derived APOGEE Stellar Parameter and Chemical Abundance Pipeline (ASPCAP) parameters for this reference set. For the regions of parameter space not well represented by the APOGEE training set ($7,000 \leq T \leq 30,000$ K), we supplement with theoretical model spectra. We present our derived parameters along with an analysis of the uncertainties and comparisons to other analyses from the literature.

Rapidly rotating neutron stars are promising sources for existing and upcoming gravitational-wave interferometers. While relatively dim, these systems are expected to emit continuously, allowing for signal to be accumulated through persistent monitoring over year-long timescales. If, at some point during the observational window, the source comes to lie behind a dense collection of stars, transient gravitational lensing may occur. Such events, though rare, would modulate the waveform, induce phase drifts, and ultimately affect parameter inferences concerning the nuclear equation of state and/or magnetic field structure of the neutron star. Importantly, the source wavelength will typically exceed the Schwarzschild radius of the individual perturbers in this scenario, implying that (micro-)lensing occurs in the diffractive regime where geometric optics does not apply. In this paper, we make use of numerical tools that borrow from Picard-Lefschetz theory to efficiently evaluate the relevant Fresnel-Kirchhoff integrals for $n \gtrsim 10^{2}$ microlenses. Modulated strain profiles are constructed both in general and for particular neutron star trajectories relative to some simulated macrolenses.

We analyzed the evolution of 3GHz luminosity function of radio selected galaxies through cosmic time, up to about 1 Gyr after Big Bang, with making use of the VLA-COSMOS3GHz large project. We estimated the luminosity function for SFGs and AGNs in the VLA-COSMOS field respectively with two independent methods, the non-parametric C-method for accuracy to avoid the effect of cosmic variance and parametric MCMC method assuming pure luminosity evolution. We obtained pure luminosity evolution parameter as L*=(1+z)(3.36-0.31)z for SFG and L* = (1+z)(2.91-0.99)z for AGN. The resultant model agreed well with observed number counts down to sub-muJy level. Additionally, we compared the model with cosmic star formation rate density history with multiple observations in various wavelength and confirmed that our result is consistent with literature up to z=3. We further analyzed the statistical effect of gravitational lensing due to the dark matter halos and found that the effect is maximized at S_3GHz = 100muJy with about 1% contribution for SFG and less than 0.5% for AGN for the flux range considered in this study. In order to evaluate the effect of lensing magnification on parameter estimation in the SKA observations, we performed Fisher analysis and find that while the bias is limited, there is a distinguishable effect on SFG LF evolution parameter with SKA I MID wide survey, and all-sky and wide survey for AGN LF evolution parameters.

Due to their scarcity, microlensing events in the Galactic disk are of great interest and high-cadence photometric observations, supplemented by spectroscopic follow-up, are necessary for constraining the physical parameters of the lensing system. In particular, a precise estimate of the source characteristics is required to accurately measure the lens distance and mass. We conducted a spectroscopic follow-up of microlensing event Gaia19bld to derive the properties of the microlensing source and, ultimately, to estimate the mass and distance of the lens. We obtained low- and high-resolution spectroscopy from multiple sites around the world during the course of the event. The spectral lines and template matching analysis has led to two independent, consistent characterizations of the source. We found that the source is a red giant located at about 8.5 kpc from the Earth. Combining our results with the photometric analysis has led to a lens mass of Ml=1.1 M at a distance of D l = 5.5 kpc. We did not find any significant blend light in the spectra (with an upper detection limit of V < 17 mag), which is in agreement with photometric observations. Therefore, we cannot exclude the possibility that the lens is a main-sequence star. Indeed, we predict in this scenario a lens brightness of V about 20 mag, a value that would make it much fainter than the detection limit.

Asteroid pairs are genetically related asteroids that recently separated ($<$few million years), but still reside on similar heliocentric orbits. A few hundred of these systems have been identified, primarily in the asteroid main-belt. Here we studied a newly discovered pair of near-Earth objects (NEOs): 2019 PR2 and 2019 QR6. Based on broadband photometry, we found these asteroids to be spectrally similar to D-types, a type rare amongst NEOs. We recovered astrometric observations for both asteroids from the Catalina Sky Survey from 2005, which significantly improved their fitted orbits. With these refinements we ran backwards orbital integrations to study formation and evolutionary history. We found that neither a pure gravitational model nor a model with the Yarkovsky effect could explain their current orbits. We thus implemented two models of comet-like non-gravitational forces based on water or CO sublimation. The first model assumed quasi-continuous, comet-like activity after separation, which suggested a formation time of the asteroid pair $300^{+120}_{-70}$ years ago. The second model assumed short-term activity for up to one heliocentric orbit ($\sim$13.9 years) after separation, which suggested that the pair formed 272$\pm$7 years ago. Image stacks showed no activity for 2019~PR2 during its last perihelion passage. These results strongly argue for a common origin that makes these objects the youngest asteroid pair known to date. Questions remain regarding whether these objects derived from a parent comet or asteroid, and how activity may have evolved since their separation.

We present the environmental properties around high-$z$ radio galaxies (HzRGs) at $z\sim4$, which have been poorly investigated because of their rarity. We use the largest samples of HzRGs and $g$-dropout galaxy overdense regions at $z\sim4$, which were constructed from Hyper Suprime-Cam Subaru Strategic Program, to characterize the HzRG environments statistically. We measure the $g$-dropout galaxy overdensities around 21 HzRGs whose rest-frame 1.4 GHz radio luminosities ($L_{1.4\mathrm{GHz}}$) are $10^{26-27}$ W Hz$^{-1}$. We find that the overdensities around the faint HzRGs with $L_{1.4\mathrm{GHz}}\sim10^{26.0-26.5}$ W Hz$^{-1}$ tend to be higher than that of the $g$-dropout galaxies. On the other hand, no significant difference of density environments is found between the luminous HzRGs with $ L_{1.4\mathrm{GHz}}\sim10^{26.5-27.0} $ W Hz$^{-1}$ and the $g$-dropout galaxies. The HzRGs are found to occupy more massive halos than $g$-dropout galaxies through a cross-correlation between the HzRGs and $g$-dropout galaxies. This trend is more pronounced in the faint HzRGs. These results are consistent with a scenario where HzRGs get older and more massive as the radio-luminosity decreases. The HzRGs are expected to trace the progenitors of local cluster halos from their calculated halo mass. In addition, we find that surrounding galaxies tend to distribute along the radio-jet major axis of the HzRGs at angular distances less than $\lesssim500$ physical kpc. Our findings imply the onset of the filamentary structures around the HzRGs at $z\sim4$.

We report the performance of a radiation hardened SRI 4K$\times$2K CMOS image sensor at temperatures down to 140 K. A biasing scheme to suppress readout glow during long exposures is highly effective so that the 0.077 $\mathrm{m\mathrm{e^-}/s}$ dark current floor is reached at 160 K, rising to 1 $\mathrm{m}\mathrm{e^-}\mathrm{s}$ at 184 K. We examine the trade-off between readout speed and read noise, finding that 1.43 $\mathrm{e^-}$ median read noise is achieved using line-wise digital correlated double sampling at 700 $\mathrm{kpix/s/ch}$ which corresponds to 1.5 $\mathrm{s}$ readout time. The 15 $\mathrm{k}\mathrm{e^-}$ well capacity in high gain mode is extended to 120 $\mathrm{k}\mathrm{e^-}$ in dual gain mode. Continued collection of photo-generated charge during readout makes possible a further dynamic range extension beyond $10^6\,\mathrm{e^-}$ effective well capacity with only $1\%$ loss of exposure efficiency by combining short and long exposures. A quadratic fit to correct for non-linearity reduces gain correction residuals from $1.5\%$ to $0.2\%$ in low gain mode and to $0.4\%$ in high gain mode. Crosstalk to adjacent pixels is only $0.4\%$ vertically, $0.6\%$ horizontally and $0.1\%$ diagonally. These characteristics plus the relatively large ($10\mathrm{\mu}\mathrm{m}$) pixel size, quasi 4-side buttability, electronic shutter and subarray readout make this sensor an excellent choice for wide field imaging in space, even at FUV wavelengths where sky background is very low.

In this work we combine information from solar-like oscillations, high-resolution spectroscopy and Gaia astrometry to derive stellar ages, chemical abundances and kinematics for a group of seven metal-poor Red Giants and characterise them in a multidimensional chrono-chemo-dynamical space. Chemical abundance ratios were derived through classical spectroscopic analysis employing 1D LTE atmospheres on Keck/HIRES spectra. Stellar ages, masses and radii were calculated with grid-based modelling, taking advantage of availability of asteroseismic information from Kepler. The dynamical properties were determined with Galpy using Gaia EDR3 astrometric solutions. Our results suggest that underestimated parallax errors make the effect of Gaia parallaxes more important than different choices of model grid or -- in the case of stars ascending the RGB -- mass-loss prescription. Two of the stars in this study are identified as potentially evolved halo blue stragglers. Four objects are likely members of the accreted Milky Way halo, and their possible relationship with known accretion events is discussed.

Solar wind turbulence is anisotropic with respect to the mean magnetic field. Anisotropy leads to ambiguity when interpreting in-situ turbulence observations in the solar wind because an apparent change in the measurements could be due either to the change of intrinsic turbulence properties or to a simple change of the spacecraft sampling direction. We demonstrate the ambiguity using the spectral index and magnetic compressibility in the inertial range observed by the Parker Solar Probe during its first seven orbits ranging from 0.1 to 0.6 AU. To unravel the effects of the sampling direction, we assess whether the wavevector anisotropy is consistent with a two-dimensional (2D) plus slab turbulence transport model and determine the fraction of power in the 2D versus slab component. Our results confirm that the 2D plus slab model is consistent with the data and the power ratio between 2D and slab components depends on radial distance, with the relative power in 2D fluctuations becoming smaller closer to the Sun.

The Gaia Alert System issued an alert on 2020 August 28, on Gaia 20eae when its light curve showed a $\sim$4.25 magnitude outburst. We present multi-wavelength photometric and spectroscopic follow-up observations of this source since 2020 August and identify it as the newest member of the FUor/EXor family of sources. We find that the present brightening of Gaia 20eae is not due to the dust clearing event but due to an intrinsic change in the spectral energy distribution. The light curve of Gaia 20eae shows a transition stage during which most of its brightness ($\sim$3.4 mag) has occurred at a short timescale of 34 days with a rise-rate of 3 mag/month. Gaia 20eae has now started to decay at a rate of 0.3 mag/month. We have detected a strong P Cygni profile in H$\alpha$ which indicates the presence of winds originating from regions close to the accretion. We find signatures of very strong and turbulent outflow and accretion in Gaia 20eae during this outburst phase. We have also detected a red-shifted absorption component in all the Ca II IR triplet lines consistent with signature of hot in-falling gas in the magnetospheric accretion funnel. This enables us to constrain the viewing angle with respect to the accretion funnel. Our investigation of Gaia 20eae points towards magnetospheric accretion being the phenomenon for the current outburst.

Blazars are a class of active galactic nuclei which host relativistic jets oriented close to the observer's line of sight. Blazars have very complex variability properties. Flares, namely flux variations around the mean value with a well-defined shape and duration, are one of the identifying properties of the blazar phenomenon. Blazars are known to exhibit multi-wavelength flares, but also "orphan" flares, namely flux changes that appear only in a specific energy range. Various models, sometimes at odds with each other, have been proposed to explain specific flares even for a single source, and cannot be synthesized into a coherent picture. In this paper, we propose a unified model for explaining orphan and multi-wavelength flares from blazars in a common framework. We assume that the blazar emission during a flare consists of two components: (i) a quasi-stable component that arises from the superposition of numerous but comparatively weak dissipation zones along the jet, forming the background (low-state) emission of the blazar, and (ii) a transient component, which is responsible for the sudden enhancement of the blazar flux, forming at a random distance along the jet by a strong energy dissipation event. Whether a multi-wavelength or orphan flare is emitted depends on the distance from the base of the jet where the dissipation occurs. Generally speaking, if the dissipation occurs at a small/large distance from the supermassive black hole, the inverse Compton/synchrotron radiation dominates and an orphan gamma-ray/optical flare tends to appear. On the other hand, we may expect a multi-wavelength flare if the dissipation occurs at an intermediate distance. We show that the model can successfully describe the spectral energy distribution of different flares from the flat spectrum radio quasar 3C 279 and the BL Lac object PKS 2155-304.

The origin of ultra-high energy neutrinos still lacks observational evidence, besides, the physical mechanism is also unclear. There is an association of a PeV neutrino event (IceCube-191001A) with an optical tidal disruption event (TDE, AT2019dsg) which was detected 6 months ahead from IceCube-191001A. The numerical simulations and observations suggested that a TDE can produce ultrafast outflows, which will interact with clouds near the supermassive black hole. In this paper, we study the interactions between TDE outflows and clouds and the possible production. In the shock waves generated by the outflow-cloud interactions, protons can be accelerated to $\sim$ 60 PeV with the outflow velocity 0.07 c and kinetic luminosity $10^{45}\rm erg/s$. PeV neutrinos can be produced through hadronic reactions. The calculation illustrates that the expected PeV neutrino event from AT2019dsg is 0.014 for a power-law proton energy distribution of $\Gamma=1.5$ and 0.0016 for $\Gamma=1.9$. The GeV --TeV $\gamma$-rays through hadronic processes are lower than the present observed limits. Outflows escaped from the TDE center colliding with clouds, which also can naturally explain the half-year delay between the neutrino event and TDE.

We follow up on the surprising recent announcement by Vernstrom et al. (2021) of the detection of the synchrotron cosmic web. We attempt to reproduce their detection with new observations with the Phase II, extended configuration of the Murchison Widefield Array at 118.5 MHz. We reproduce their detection methodology by stacking pairs of nearby luminous red galaxies (LRGs) -- used as tracers for clusters and galaxy groups -- contained in our low frequency radio observations. We show our observations are significantly more sensitive than those used in Vernstrom et al., and that our angular sensitivity is sufficient. And yet, we make no statistically significant detection of excess radio emission along the bridge spanning the LRG pairs. This non-detection is true both for the original LRG pair catalogue as used in Vernstrom et al., as well as for other larger catalogues with modified selection criteria. Finally, we return to the original data sets used in Vernstrom et al., and find that whilst we clearly reproduce the excess X-ray emission from ROSAT, we find no evidence of intercluster filamentary emission in the original 118.5 MHz MWA survey data. In the interests of understanding this result, as part of this paper we release images of the 14 fields used in this study, the final stacked images, as well as key components of our stacking and modelling code.

We present the results of the Fermi-Large Area Telescope (LAT) light curve (LC) modelling of selected blazars: six flat spectrum radio quasars (FSRQs) and five BL Lacertae (BL Lacs). All objects have densely sampled and long-term LCs, over 10 years. For each blazar we generated three LCs with 7, 10, and 14 days binning, using the latest LAT 8-year source catalog and binned analysis provided within the fermipy package. The LCs were modelled with several tools: the Fourier transform, the Lomb-Scargle periodogram (LSP), the autoregressive moving average (ARMA), the fractional autoregressive integrated moving average, the continuous-time autoregressive moving average (CARMA) processes, the Hurst exponents ($H$), the $\mathcal{A}-\mathcal{T}$ plane, and the wavelet scalogram. The power law indices $\beta$ calculated from the Fourier and LSP modelling are consistent with each other. Many objects yield $\beta\simeq 1$, with PKS 2155$-$304 even flatter, but some are significantly steeper, e.g. Mrk 501 and B2 1520+31. The power law PSD is indicative of a self-affine stochastic process characterised by $H$, underlying the observed variability. Several algorithms for the $H$ estimations are employed. For some objects we found $H>0.5$, indicating long-term memory. We confirm a quasi-periodic oscillation (QPO) in the PKS 2155$-$304 data with the period of $612\pm 42$ days at a $3\sigma$ significance level, but do not detect any QPOs in other objects. The ARMA results give in general higher orders for 7 days binned LCs and lower orders for 10 and 14 days binned LCs, implying temporal variations in the LCs are consistently captured by the fitted models. CARMA modelling leads to featureless PSDs. The recently introduced $\mathcal{A}-\mathcal{T}$ plane allows us to successfully classify the PSDs based on the LCs alone and clearly separates the FSRQ and BL Lac types of blazars.

The origin of magnetic fields in clusters of galaxies is still a matter of debate. Observations for intracluster magnetic fields over a wide range of redshifts are crucial to constrain possible scenarios for the origin and evolution of the fields. (Differences of Faraday rotation measures (RMs) of an embedded double radio sources, i.e. a pair of lobes of mostly Fanaroff--Riley type II radio galaxies, are free from the Faraday rotation contributions from the interstellar medium inside the Milky Way and the intergalactic medium between radio galaxies and us, and hence provide a novel way to estimate average magnetic field within galaxy clusters.) We have obtained a sample of 627 pairs whose RMs and redshifts are available in the most updated RM catalogues and redshift databases. The RM differences of the pairs are derived. The statistically large RM differences for pairs of redshifts $z>0.9$ indicate that intracluster magnetic fields is as strong as about 4~$\mu$G. Such strong magnetic fields in the intracluster medium at the half age of the Universe, comparable to intracluster field strength in nearby galaxy clusters, pose a challenge on the theories for origin of cosmic magnetic fields.

This work analyzes the Hall magnetohydrodynamics (HMHD) and magnetohydrodynamics (MHD) numerical simulations of a flaring solar active region as a testbed while idealizing the coronal Alfv\'en speed to be of two orders of magnitude lesser. HMHD supports faster magnetic reconnection and shows richer complexity in magnetic field line evolution compared to the MHD. The magnetic reconnections triggering the flare are explored by numerical simulations augmented with relevant multi-wavelength observations. The initial coronal magnetic field is constructed by non-force-free extrapolation of photospheric vector magnetic field. Magnetic structure involved in the flare is identified to be a flux rope, with its overlying magnetic field lines constituting the quasi-separatrix layers (QSLs) along with a three-dimensional null point and a null line. Compared to the MHD simulation, the HMHD simulation shows a higher and faster ascend of the rope together with the overlying field lines, which further reconnect at the QSL located higher up in the corona. The foot points of the field lines match better with the observations for the HMHD case with the central part of the flare ribbon located at the chromosphere. Additionally, field lines are found to rotate in a circular pattern in the HMHD, whereas no such rotation is seen in the MHD results. Interestingly, plasma is also observed to be rotating in a co-spatial chromospheric region, which makes the HMHD simulation more credible. Based on the aforementioned agreements, HMHD simulation is found to agree better with observations and, thus, opens up a novel avenue to explore.

Formation process(es) of galactic bulges are not yet clarified although several mechanisms have been proposed. In a previous study, we suggested one possibility that galactic bulges have been formed from the cold gas inflowing through surrounding hot halo gas in massive dark matter halos at high redshifts. It was shown that this scenario leads to the bulge-to-total stellar mass ratio increasing with the galaxy mass, in agreement with the well-known observed trend. We here indicate that it also reproduces recent observational results that the mean stellar age of the bulge increases with the galaxy mass while the age gradient across the bulge decreases. We infer that this formation path applies mainly to high-mass galaxies and the bulges in lower-mass galaxies have different origins such as secular formation from the disc material.

We derive general expressions for the multi-tracer Fisher matrix, both assuming that the cross-spectra are constrained by the auto-spectra, and also allowing for independent degrees of freedom in the cross-spectra. We show that, just like the ratios of power spectra, the independent degrees of freedom of the cross-spectra are also not constrained by cosmic variance. Moreover, whereas the uncertainties in the ratios of power spectra decrease with the number density of the tracers as $\sim 1/\sqrt{\bar{n}}$, the uncertainties in the independent degrees of freedom of the cross-spectra decrease even faster, as $\sim 1/\bar{n}$. We also derive simple expressions for the optimal number of tracers in a survey.

The space missions designed to visit small bodies of the Solar System boosted the study of the dynamics around non-spherical bodies. In this vein, we study the dynamics around a class of objects classified by us as Non-Spherical Symmetric Bodies, including contact binaries, triaxial ellipsoids, spherical bodies with a mass anomaly, among others. In the current work, we address the results for a body with a mass anomaly. We apply the pendulum model to obtain the width of the spin-orbit resonances raised by non-asymmetric gravitational terms of the central object. The Poincare surface of section technique is adopted to confront our analytical results and to study the system's dynamics by varying the parameters of the central object. We verify the existence of two distinct regions around an object with a mass anomaly: a chaotic inner region that extends beyond the corotation radius and a stable outer region. In the latter, we identify structures remarkably similar to those of the classical restrict and planar 3-body problem in the Poincare surface of sections, including asymmetric periodic orbits associated with 1:1+p resonances. We apply our results to a Chariklo with a mass anomaly, obtaining that Chariklo rings are probably related to first kind periodic orbits and not with 1:3 spin-orbit resonance, as proposed in the literature. We believe that our work presents the first tools for studying mass anomaly systems.

The nature of dark matter (DM) is still an open question in Physics. Gamma-ray and neutrino telescopes have been searching for DM signatures for several years and no detection has been obtained so far. In their quest, these telescopes have gathered a wealth of observations that, if properly combined and analyzed, can improve on the constraints to the nature of DM set by individual instruments. In this contribution, we present two open-source analysis tools aimed at performing the before mentioned combined analysis: gLike, a general-purpose ROOT-based code framework for the numerical maximization of joint likelihood functions, and LklCom, a Python-based tool combining likelihoods from different instruments to produce combined exclusion limits on the DM annihilation cross-section.

The Major Atmospheric Gamma Imaging Cherenkov (MAGIC) telescope system consists of two imaging atmospheric Cherenkov telescopes (IACTs) and is located on the Canary island of La Palma. IACTs are excellent tools to inspect the very-high-energy (few tens of GeV and above) gamma-ray sky by capturing images of the air showers, originated by the absorption of gamma rays and cosmic rays by the atmosphere, through the detection of Cherenkov photons emitted in the shower. One of the main factors determining the sensitivity of IACTs to gamma-ray sources, in general, is how well reconstructed the properties (type, energy, and incoming direction) of the primary particle triggering the air shower are. We present how deep convolutional neural networks (CNNs) are being explored as a promising method for IACT full-event reconstruction. The performance of the method is evaluated on observational data using the standard MAGIC Analysis and Reconstruction Software, MARS, and CTLearn, a package for IACT event reconstruction through deep learning.

We use a sample of 706 galaxies, selected as [OII]$\lambda$3727 ([OII]) emitters in the Survey for High-$z$ Absorption Red and Dead Sources (SHARDS) on the CANDELS/GOODS-N field, to study the differential attenuation of the nebular emission with respect to the stellar continuum. The sample includes only galaxies with a counterpart in the infrared and $\mathrm{log}_{10}(M_{*}/\mathrm{M}_{\odot})$ $>$ 9, over the redshift interval 0.3 $\lesssim$ $z$ $\lesssim$ 1.5. Our methodology consists in the comparison of the star formation rates inferred from [OII] and H$\alpha$ emission lines with a robust quantification of the total star-forming activity (${SFR}_{\mathrm{TOT}}$) that is independently estimated based on both infrared and ultraviolet (UV) luminosities. We obtain $f$$=$$E(B-V)_{\mathrm{stellar}}$/$E(B-V)_{\mathrm{nebular}}$ $=$ 0.69$^{0.71}_{0.69}$ and 0.55$^{0.56}_{0.53}$ for [OII] and H$\alpha$, respectively. Our resulting $f$-factors display a significant positive correlation with the UV attenuation and shallower or not-significant trends with the stellar mass, the $SFR_{\mathrm{TOT}}$, the distance to the main sequence, and the redshift. Finally, our results favour an average nebular attenuation curve similar in shape to the typical dust curve of local starbursts.

Laboratory results of the optical properties of vapor-deposited water ice, specifically the refractive index and extinction coefficient, are available mainly for a selective set of wavelengths and a limited number of deposition temperatures. Experimental limitations are the main reason for the lack of broadband data, which is unfortunate as these quantities are needed to interpret and predict astronomical and planetary observations. The goal of this work is to address these lacking data, using an experimental broadband method that is capable of rapidly providing reliable water ice data across the entire UV-visible range. This approach combines the simultaneous use of a monochromatic HeNe laser and a broadband Xe-arc lamp to record interference fringes of water ice during deposition at astronomically relevant ice temperatures. The ice thickness is typically more than 20 $\mu$m. Analyzing the period and intensity patterns combining both the monochromatic and broadband interference patterns allows the determination of the wavelength-dependent refractive index and extinction coefficient. We present accurate refractive index and extinction coefficient graphs for wavelengths between 250 and 750 nm and ices deposited between 30 and 160 K. From our data, we find a possible structural change in the ice in the 110-130 K region that has not been reported before. We also discuss that the data presented in this paper can be used to interpret astronomical observations of icy surfaces.

Despite the huge improvements guaranteed by future GRAVITY observations of S0-2 star, these will not be able to unveil the fundamental nature, whether black hole or wormhole, of the central supermassive object. Nevertheless, observing stars orbiting closer to the central gravitational source could allow to distinguish between the black hole and wormhole nature of this object at more than 5$\sigma$. Firstly, we have used publicly available astrometric and spectroscopic measurements of S0-2 star to constrain the metric around the supermassive object without finding any evidence either favouring or ruling out the wormhole nature. Secondly, we have designed a mock catalogue of future observations of S0-2 star mirroring the accuracy and precision of GRAVITY. Afterwards, we firstly tested our methodology showing that our procedure recovers the input model, and subsequently we demonstrated that the constraining power of such a dataset is not enough to distinguish between black hole and wormhole. Additionally, we evaluate whether future observation of the star S62, whose pericentre is located at only 16 AU away from the central object, would allow to shed light on its nature. Nevertheless, we conclude that the star S62 will not serve to achieve our aim. Finally, we built some toy models representing stars orbiting much closer the central object than S0-2 and S62 stars. We used these toy models to investigate which are the ideal orbital features and observational strategies to achieve our aim of unveiling the fundamental nature of the central supermassive object, demonstrating that a star with a period of the order of $\sim 5$ years and a pericentre distance of $\sim 5$ AU could identify the nature of the central object at almost 5$\sigma$ accuracy.

Using the galaxy clusters from \theth, we define a new parameter: $\lambda_{DS}$ to describe the dynamical state of clusters, which assumes a double-Gaussian distribution in logarithm scale for our mass-complete cluster sample at $z=0$ from the dark-matter-only (DMO) run. Therefore, the threshold for distinguishing relaxed and unrelaxed clusters is naturally determined by the crossing point of the double-Gaussian fitting which has a value of $\lambda_{DS} = 3.424$. By applying $\lambda_{DS}$ with the same parameters from the DMO run to the hydro-dynamically simulated clusters (\gadgetx\ run and \simba\ run), we investigate the effect of baryons on the cluster dynamical state. We find a weak baryon-model dependence for the $\lambda_{DS}$ parameter. Finally, we study the evolution of $\lambda_{DS}$ along with clusters mass accretion history. We notice an upper limit of halo mass change $\frac{\Delta M_{200}}{M_{200}} \sim 0.12$ that don't alter the cluster dynamical state, i.e. from relaxed to unrelaxed. We define relaxation period (from the most relaxed state to disturb and relaxed again) which reflects how long the dynamical state of a cluster restores its relaxation state, and propose a correlation between this relaxation period and the strength of halo mass change $\frac{\Delta M_{200}}{M_{200}}$. With the proposed fitting to such correlation, we verify the relaxation period can be estimated from $\frac{\Delta M_{200}}{M_{200}}$ (including multi mass change peaks) with considerably small error.

The relative rarity of giant planets around low mass stars compared with solar-type stars is a key prediction from core accretion planet formation theory. In this paper we report on the discovery of four gas giant planets that transit low mass late K and early M dwarfs. The planets HATS-74Ab (TOI 737b), HATS-75b (TOI 552b), HATS-76b (TOI 555b), and HATS-77b (TOI 730b), were all discovered from the HATSouth photometric survey and followed-up using TESS and other photometric facilities. We use the new ESPRESSO facility at the VLT to confirm and systems and measure their masses. We find that that planets have masses of 1.46 +- 0.14 MJ , 0.491 +- 0.039 MJ , 2.629 +- 0.089 MJ and 1.374 +0.100-0.074 MJ , respectively, and radii of 1.032 +- 0.021 RJ , 0.884 +- 0.013 RJ , 1.079 +- 0.031 RJ , and 1.165 +- 0.021 RJ, respectively. The planets all orbit close to their host stars with orbital periods ranging from 1.7319 d to 3.0876 d. With further work we aim to test core accretion theory by using these and further discoveries to quantify the occurrence rate of giant planets around low mass host stars.

The collisional pumping of CH3OH and OH masers in non-dissociative C-type shock waves is studied. The chemical processes responsible for the evolution of molecule abundances in the shock wave are considered in detail. The large velocity gradient approximation is used to model radiative transfer in molecular lines. We present calculations of the optical depth in maser transitions of CH3OH and OH for a grid of C-type shock models that vary in cosmic ray ionization rate, gas density and shock speed. We show that pre-shock gas densities $n_{H,tot} = 2 \times 10^4-2 \times 10^5$ cm$^{-3}$ are optimal for pumping of methanol maser transitions. A complete collisional dissociation of methanol at the shock front takes place for shock speeds $u_{s} \gtrsim 25$ km s$^{-1}$. At high pre-shock gas density $n_{H,tot} = 2 \times 10^{6}$ cm$^{-3}$, the collisional dissociation of methanol takes place at shock speeds just above the threshold speed $u_{s} \approx 15-17.5$ km s$^{-1}$ corresponding to sputtering of icy mantles of dust grains. We show that the methanol maser transition E $4_{-1} \to 3_0$ at 36.2 GHz has the optical depth higher than that of the transition A$^+$ $7_0 \to 6_1$ at 44.1 GHz at high cosmic ray ionization rate $\zeta_{H_2} \gtrsim 10^{-15}$ s$^{-1}$ and pre-shock gas density $n_{H,tot} = 2 \times 10^4$ cm$^{-3}$. These results can be applied to the interpretation of observational data on methanol masers near supernova remnants and in molecular clouds of the Central Molecular Zone. At the same time, a necessary condition for the operation of 1720 MHz OH masers is a high ionization rate of molecular gas, $\zeta_{H_2} \gtrsim 10^{-15}$ s$^{-1}$. We find that physical conditions conducive to the operation of both hydroxyl and methanol masers are cosmic ray ionization rate $\zeta_{H_2} \approx 10^{-15} - 3 \times 10^{-15}$ s$^{-1}$, and a narrow range of shock speeds.

A powerful test of fundamental physics consists on probing the variability of fundamental constants in Nature. Although they have been measured on Earth laboratories and in our Solar neighbourhood with extremely high precision, it is crucial to carry out these tests at the distant Universe, as any significant variation of these quantities would immediately hint at new physics. We perform a cosmological measurement of the speed of light using the latest Type Ia Supernova and cosmic chronometer observations at the redshift range $0<z<2$. Our method relies on the numerical reconstruction of the data in order to circumvent {\it a priori} assumptions of the underlying cosmology. We confirm the constancy of the speed of light at such redshift range, reporting a $\sim 5$\% precision measurement of $c = (3.22 \pm 0.16) \; \mathrm{km \; s}^{-1}$ in $z \simeq 1.60$ at a $1\sigma$ confidence level

The LIGO-Virgo-KAGRA Collaboration recently detected gravitational waves (GWs) from the merger of black-hole-neutron-star (BHNS) binary systems GW200105 and GW200115. No coincident electromagnetic (EM) counterparts were detected. While the mass ratio and BH spin in both systems were not sufficient to tidally disrupt the NS outside of the BH event horizon, other, magnetospheric mechanisms for EM emission exist in this regime and depend sensitively on the NS magnetic field strength. Combining GW measurements with EM flux upper limits, we place upper limits on the NS surface magnetic field strength above which magnetospheric emission models would have generated an observable EM counterpart. We consider fireball models powered by the black-hole battery mechanism, where energy is output in gamma-rays over $\lesssim1$~second. Consistency with no detection by Fermi-GBM or INTEGRAL SPI-ACS constrains the NS surface magnetic field to $\lesssim10^{15}$~G. Hence, joint GW detection and EM upper limits rule out the theoretical possibility that the NSs in GW200105 and GW200115, and the putative NS in GW190814, retain $\gtrsim10^{15}$~G dipolar magnetic fields until merger. They also rule out formation scenarios where strongly magnetized magnetars quickly merge with BHs. We alternatively rule out operation of the BH-battery powered fireball mechanism in these systems. This is the first multi-messenger constraint on NS magnetic fields in BHNS systems and a novel approach to probe fields at this point in NS evolution. This demonstrates the constraining power that multi-messenger analyses of BHNS mergers have on BHNS formation scenarios, the magnetic-field evolution in NSs, and the physics of BHNS magnetospheric interactions.

Non-thermal radio emission is detected in the region of the gamma-ray source FHES J1723.5-0501. The emission has an approximately circular shape 0.8 degrees in diameter. The observations confirm its nature as a new supernova remnant, G17.8+16.7. We derive constraints on the source parameters using the radio data and gamma-ray observations of the region. The distance to the object is possibly in the range 1.4-3.5 kpc. An SNR age of the order of 10 kyr is compatible with the radio and GeV features, but an older or younger SNR cannot be ruled out. A simple one-zone leptonic model naturally explains the multi-wavelength non-thermal fluxes of the source at its location outside the Galactic plane.

OAO 1657-415 is an atypical supergiant X-ray binary among wind-fed and disk-fed systems, showing alternate spin-up/spin-down intervals lasting on the order of tens of days. We study different torque states of OAO 1657-415 based on the spin history monitored by {\it Fermi}/GBM, together with fluxes from {\it Swift}/BAT and {\it MAXI}/GSC. Its spin frequency derivatives are well correlated with {\it Swift}/BAT fluxes during rapid spin-up episodes, anti-correlated with {\it Swift}/BAT fluxes during rapid spin-down episodes, and not correlated in between. The orbital profile of spin-down episodes is reduced by a factor of 2 around orbital phases of 0.2 and 0.8 compared to that of spin-up episodes. The orbital hardness ratio profile of spin-down episodes is also lower than that of spin-up episodes around phases close to the mid-eclipse, implying that there is more material between the neutron star and the observer for spin-down episodes than for spin-up episodes around these phases. These results indicate that the torque state of the neutron star is connected with the material flow on orbital scale and support the retrograde/prograde disk accretion scenario for spin-down/spin-up torque reversal.

GW190521 is the most massive binary black hole (BBH) merger observed to date. Due to its peculiar properties, the origin of this system is still a matter of debate: several hints may favor a dense stellar environment as a birthplace. Here, we investigate the possible formation of GW190521-like systems via three-body encounters in young massive star clusters (YSCs) by means of direct N-body simulations.

This paper documents the seventeenth data release (DR17) from the Sloan Digital Sky Surveys; the fifth and final release from the fourth phase (SDSS-IV). DR17 contains the complete release of the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey, which reached its goal of surveying over 10,000 nearby galaxies. The complete release of the MaNGA Stellar Library (MaStar) accompanies this data, providing observations of almost 30,000 stars through the MaNGA instrument during bright time. DR17 also contains the complete release of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) survey which publicly releases infra-red spectra of over 650,000 stars. The main sample from the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), as well as the sub-survey Time Domain Spectroscopic Survey (TDSS) data were fully released in DR16. New single-fiber optical spectroscopy released in DR17 is from the SPectroscipic IDentification of ERosita Survey (SPIDERS) sub-survey and the eBOSS-RM program. Along with the primary data sets, DR17 includes 25 new or updated Value Added Catalogs (VACs). This paper concludes the release of SDSS-IV survey data. SDSS continues into its fifth phase with observations already underway for the Milky Way Mapper (MWM), Local Volume Mapper (LVM) and Black Hole Mapper (BHM) surveys.

With the launch of the JWST, we will obtain more precise data for exoplanets than ever before. However, this data can only inform and revolutionize our understanding of exoplanets when placed in the larger context of planet-star formation. Therefore, gaining a deeper understanding of their host stars is equally important and synergistic with the upcoming JWST data. We present detailed chemical abundance profiles of 17 FGK stars that will be observed in exoplanet-focused Cycle 1 JWST observer programs. The elements analyzed (C, N, O, Na, Mg, Si, S, K, and Fe) were specifically chosen as being informative to the composition and formation of planets. Using archival high-resolution spectra from a variety of sources, we perform an LTE equivalent width analysis to derive these abundances. We look to literature sources to correct the abundances for non-LTE effects, especially for O, S, and K, where the the corrections are large (often $> 0.2~\textrm{dex}$). With these abundances and the ratios thereof, we will begin to paint clearer pictures of the planetary systems analyzed by this work. With our analysis, we can gain insight into the composition and extent of migration of Hot Jupiters, as well as the possibility of carbon-rich terrestrial worlds.

We report on the detection of a large, extended HI cloud complex in the GAMA G23 field, located at a redshift of $z\,\sim\,0.03$, observed as part of the MeerHOGS campaign (a pilot survey to explore the mosaicing capabilities of MeerKAT). The cloud complex, with a total mass of $10^{10.0}\,M_\odot$, lies in proximity to a large galaxy group with $M_\mathrm{dyn}\sim10^{13.5}\,M_\odot$. We identify seven HI peak concentrations, interconnected as a tenuous 'chain' structure, extending $\sim 400\,\mathrm{kpc}$ from east-to-west, with the largest (central) concentration containing $10{^{9.7}}\,M_\odot$ in HI gas distributed across $50\,\mathrm{kpc}$. The main source is not detected in ultra-violet, optical or infrared imaging. The implied gas mass-to-light ($M_\mathrm{HI}$/$L_\mathrm{r}$) is extreme ($>$1000) even in comparison to other 'dark clouds'. The complex has very little kinematic structure ($110\,\mathrm{km}\,\mathrm{s}^{-1}$), making it difficult to identify cloud rotation. Assuming pressure support, the total mass of the central concentration is $>10^{10.2}\,M_\odot$, while a lower limit to the dynamical mass in the case of full rotational support is $10^{10.4}\,M_\odot$. If the central concentration is a stable structure, it has to contain some amount of unseen matter, but potentially less than is observed for a typical galaxy. It is, however, not clear whether the structure has any gravitationally stable concentrations. We report a faint UV--optical--infrared source in proximity to one of the smaller concentrations in the gas complex, leading to a possible stellar association. The system nature and origins is enigmatic, potentially being the result of an interaction with or within the galaxy group it appears to be associated with.

The mass function of globular cluster (GC) populations is a fundamental observable that encodes the physical conditions under which these massive stellar clusters formed and evolved. The high-mass end of star cluster mass functions are commonly described using a Schechter function, with an exponential truncation mass $M_{c,*}$. For the GC mass functions in the Virgo galaxy cluster, this truncation mass increases with galaxy mass ($M_{*}$). In this paper we fit Schechter mass functions to the GCs in the most massive galaxy group ($M_{\mathrm{200}} = 5.14 \times 10^{13} M_{\odot}$) in the E-MOSAICS simulations. The fiducial cluster formation model in E-MOSAICS reproduces the observed trend of $M_{c,*}$ with $M_{*}$ for the Virgo cluster. We therefore examine the origin of the relation by fitting $M_{c,*}$ as a function of galaxy mass, with and without accounting for mass loss by two-body relaxation, tidal shocks and/or dynamical friction. In the absence of these mass-loss mechanisms, the $M_{c,*}$-$M_{*}$ relation is flat above $M_* > 10^{10} M_{\odot}$. It is therefore the disruption of high-mass GCs in galaxies with $M_{*}\sim 10^{10} M_{\odot}$ that lowers the $M_{c,*}$ in these galaxies. High-mass GCs are able to survive in more massive galaxies, since there are more mergers to facilitate their redistribution to less-dense environments. The $M_{c,*}-M_*$ relation is therefore a consequence of both the formation conditions of massive star clusters and their environmentally-dependent disruption mechanisms.

Temperate terrestrial exoplanets are likely common objects, but their discovery and characterization is very challenging. Concepts for optimized space missions to overcome these challenges are being studied. The LIFE initiative focuses on the development of a space-based mid-infrared (MIR) nulling interferometer probing the thermal emission of a large sample of exoplanets. We derive first estimates for the signal-to-noise (S/N), spectral resolution (R), and wavelength requirements for LIFE. Using an Earth-twin exoplanet as reference case, we quantify how well planetary/atmospheric properties can be constrained from MIR spectra of different quality. We simulate LIFE observations of an Earth-twin orbiting a G2V star at 10 pc from the Sun with different levels of exozodiacal dust emissions. We combine a cloud-free 1D radiative transfer model and the nested sampling algorithm to retrieve planetary/atmospheric properties from input spectra of different wavelength coverage, R, and S/N. We find that H2O, CO2, and O3 are detectable if S/N$\geq$10 (uncertainty $\leq\pm1.0$ dex). We find upper limits for N2O (abundance $\leq10^{-3}$). CO, N2, and O2 are unconstrained in all cases. The limit for a CH4 detection is R $= 50$, S/N $=10$. We further correctly determine the exoplanet radius (uncertainty $\leq\pm10\%$), surface temperature (uncertainty $\leq\pm20$K), and surface pressure (uncertainty $\leq\pm0.5$ dex). With the current LIFE design, the observation time required to reach the specified S/N amounts to $\sim7$ weeks with 4x2m apertures. We conclude that a minimum wavelength coverage of $4-18.5\mu$m, a R of 50 and an S/N of 10 is required. With the current assumptions, the atmospheric characterization of several Earth-like exoplanets at a distance of 10 pc and within a reasonable amount of observing time will require apertures $\geq2$ meters.

We introduce a novel method for identifying the mass composition of ultra-high-energy cosmic rays using deep learning. The key idea of the method is to use a chain of two neural networks. The first network predicts the type of a primary particle for individual events, while the second infers the mass composition of an ensemble of events. We apply this method to the Monte-Carlo data for the Telescope Array Surface Detectors readings, on which it yields an unprecedented low error of 7% for 4-component approximation. The statistical error is shown to be inferior to the systematic one related to the choice of the hadronic interaction model used for simulations.

We explore a new mechanism for reproducing the Dark Matter (DM) abundance: scatterings of one DM particle on light Standard Model particles. Strong bounds on its decays can be satisfied if DM undergoes freeze-in and has a mass around or below the pion mass. This happens, for example, in theories with a right-handed neutrino interacting with charged fermions through a leptoquark exchange. These leptoquarks can be linked to the ones motivated by the B-physics anomalies if assumptions about the flavour structure are made. DM signals are unusual, with interesting possibilities for direct and indirect detection. Achieving thermal freeze-out instead requires models with more than one DM flavour, and couplings parametrically smaller than what needed by the usual pair annihilations.

The existing search tools for exploring the NASA Astrophysics Data System (ADS) can be quite rich and empowering (e.g., similar and trending operators), but researchers are not yet allowed to fully leverage semantic search. For example, a query for "results from the Planck mission" should be able to distinguish between all the various meanings of Planck (person, mission, constant, institutions and more) without further clarification from the user. At ADS, we are applying modern machine learning and natural language processing techniques to our dataset of recent astronomy publications to train astroBERT, a deeply contextual language model based on research at Google. Using astroBERT, we aim to enrich the ADS dataset and improve its discoverability, and in particular we are developing our own named entity recognition tool. We present here our preliminary results and lessons learned.

With first-principles kinetic simulations, we show that a large-scale Alfv\'en wave (AW) propagating in an inhomogeneous background decays into kinetic Alfv\'en waves (KAWs), triggering ion and electron energization. We demonstrate that the two species can access unequal amounts of the initial AW energy, experiencing differential heating. During the decay process, the electric field carried by KAWs produces non-Maxwellian features in the particle VDFs, in accordance with space observations. The process we present solely requires the interaction of a large-scale AW with a magnetic shear and may be relevant for several astrophysical and laboratory plasmas.

Our understanding of quantum field theory rests largely on explicit and controlled calculations in perturbation theory. Because of this, much recent effort has been devoted to improve our grasp of perturbative techniques on cosmological spacetimes. While scattering amplitudes in flat space at tree level are obtained from simple algebraic operations, things are harder for cosmological observables. Indeed, computing cosmological correlation functions or the associated wavefunction coefficients requires evaluating a growing number of nested time integrals already at tree level, which is computationally challenging. Here, we present a new "differential" representation of the cosmological wavefunction in de Sitter spacetime that obviates this problem for a large class of phenomenologically relevant theories. Given any tree-level Feynman-Witten diagram, we give simple algebraic rules to write down a seed function and a differential operator that transforms it into the desired wavefunction coefficient for any scale-invariant, parity-invariant theory of massless scalars and gravitons with general boost-breaking interactions. In particular, this applies to large classes of phenomenologically relevant theories such as those described by the effective field theory of inflation or solid inflation. Trading nested bulk time integrals for derivatives on boundary kinematical data provides a great computational advantage, especially for processes involving many vertices.

The temperature dependence of stellar electron-capture (EC) rates is investigated, with a focus on nuclei around $N=50$, just above $Z=28$, which play an important role during the collapse phase of core-collapse supernovae (CCSN). Two new microscopic calculations of stellar EC rates are obtained from a relativistic and a non-relativistic finite-temperature quasiparticle random-phase approximation approaches, for a conventional grid of temperatures and densities. In both approaches, EC rates due to Gamow-Teller transitions are included. In the relativistic calculation contributions from first-forbidden transitions are also included, and add strongly to the EC rates. The new EC rates are compared with large-scale shell model calculations for the specific case of $^{86}$Kr, providing insight into the finite-temperature effects on the EC rates. At relevant thermodynamic conditions for core-collapse, the discrepancies between the different calculations of this work are within about one order of magnitude. Numerical simulations of CCSN are performed with the spherically-symmetric GR1D simulation code to quantify the impact of such differences on the dynamics of the collapse. These simulations also include EC rates based on two parametrized approximations. A comparison of the neutrino luminosities and enclosed mass at core bounce shows that differences between simulations with different sets of EC rates are relatively small ($\approx 5\%$), suggesting that the EC rates used as inputs for these simulations have become well constrained.

Proportional electroluminescence (EL) is the physical effect used in two-phase dark matter detectors, to optically record in the gas phase the ionization signal produced by particle scattering in the liquid phase. In our previous works the presence of a new EL mechanism in noble gases, namely that of neutral bremsstrahlung (NBrS), was demonstrated both theoretically and experimentally, in addition to the ordinary EL mechanism due to excimer emission. In this work we show that the similar theoretical approach can apply to noble liquids, namely to liquid helium, neon, argon, krypton and xenon. In particular, the photon yields and spectra for NBrS EL in noble liquids have for the first time been calculated, using the electron energy and transport parameters obtained in the framework of Cohen-Lekner and Atrazhev theory. The relevance of the results obtained to the development of noble liquid detectors for dark matter searches and neutrino experiments is discussed.

We investigate a model of gravitational collapse of matter inhomogeneities where the latter are modelled as Bianchi type IX (BIX) spacetimes. We found that this model contains, as limiting cases, both the standard spherical collapse model and the Zeldovich solution for a 1-dimensional perturbation. We study how these models are affected by small anisotropic perturbations within the BIX potential. For the spherical collapse case, we found that the model is equivalent to a closed FLRW Universe filled with matter and two perfect fluids representing the anisotropic contributions. From the linear evolution up to the turnaround, the anisotropies effectively shift the value of the FLRW spatial curvature, because the fluids have effective Equation of State (EoS) parameters $w \approx -1/3$. Then we estimate the impact of such anisotropies on the number density of haloes using the Press-Schechter formalism. If a fluid description of the anisotropies is still valid after virialization, the averaged over time EoS parameters are $w\approx 1/3$. Using this and demanding hydrostatic equilibrium, we find a relation between the mass $M$, the average radius $R$ and the pressure $p$ of the virialized final structure. When we consider perturbations of the Zeldovich solution, our qualitative analysis suggests that the so called \textit{pancakes} exhibit oscillatory behavior, as would be expected in the case of a vacuum BIX spacetime.