Papers for Wednesday, Jan 26 2022

Recent systematic searches for massive black holes (BHs) in local dwarf galaxies led to the discovery of a population of faint Active Galactic Nuclei (AGN). We investigate the agreement of the BH and AGN populations in the Illustris, TNG, Horizon-AGN, EAGLE, and SIMBA simulations with current observational constraints in low-mass galaxies. We find that some of these simulations produce BHs that are too massive, and that the BH occupation fraction at z=0 is not inherited from the simulation seeding modeling. In simulations, the ability of BHs and their host galaxies to power an AGN depends on the subgrid modeling of the BH and of the galaxy. The fraction of AGN in low-mass galaxies is not used to calibrate the simulations, and thus is a true prediction. AGN fractions at z=0 span two orders of magnitude at fixed galaxy stellar mass in simulations, similarly to observational constraints, but uncertainties and degeneracies affect both observations and simulations. The agreement is difficult to interpret due to differences in BH masses between simulations and observations, BH occupation fraction affected by numerical choices, and an unknown obscured fraction of AGN. Our work advocates for more thorough comparisons with observations in the low-mass regime to improve the modeling of cosmological simulations, and our understanding of BH and galaxy physics in this regime. The mass of BHs, their ability to efficiently accrete gas, and the AGN occupation fraction in low-mass galaxies have important implications for the build-up of the entire BH and galaxy populations with time.

Observations of present-day mass-loss rates for close-in transiting exoplanets provide a crucial check on models of planetary evolution. One common approach is to model the planetary absorption signal during the transit in lines like He I 10830 with an isothermal Parker wind, but this leads to a degeneracy between the assumed outflow temperature $T_0$ and the mass-loss rate $\dot{M}$ that can span orders of magnitude in $\dot{M}$. In this study, we re-examine the isothermal Parker wind model using an energy-limited framework. We show that in cases where photoionization is the only heat source, there is a physical upper limit to the efficiency parameter $\varepsilon$ corresponding to the maximal amount of heating. This allows us to rule out a subset of winds with high temperatures and large mass-loss rates as they do not generate enough heat to remain self-consistent. To demonstrate the utility of this framework, we consider spectrally unresolved metastable helium observations of HAT-P-11b, WASP-69b, and HAT-P-18b. For the former two planets, we find that only relatively weak ($\dot{M}\lesssim 10^{11.5}$ g s$^{-1}$) outflows can match the metastable helium observations while remaining energetically self-consistent, while for HAT-P-18b all of the Parker wind models matching the helium data are self-consistent. Our results are in good agreement with more detailed self-consistent simulations and constraints from high-resolution transit spectra.

The hot and dense core formed in the collapse of a massive star is a powerful source of hypothetical feebly-interacting particles such as sterile neutrinos, dark photons, axion-like particles (ALPs), and others. Radiative decays such as $a\to2\gamma$ deposit this energy in the surrounding material if the mean free path is less than the radius of the progenitor star. For the first time, we use a supernova (SN) population with particularly low explosion energies as the most sensitive calorimeters to constrain this possibility. These SNe are observationally identified as low-luminosity events with low ejecta velocities and low masses of ejected $^{56}$Ni. Their low energies limit the energy deposition from particle decays to less than about 0.1 B, where $1~{\rm B~(bethe)}=10^{51}~{\rm erg}$. For 1-500 MeV-mass ALPs, this generic argument excludes ALP-photon couplings $G_{a\gamma\gamma}$ in the $10^{-10}$-$10^{-8}~{\rm GeV}^{-1}$ range.

We present a detailed X-ray spectroscopic study of the supernova remnant DEM L71 in the Large Magellanic Cloud. Based on deep $\sim 103$ ks archival {\it Chandra} data, we perform a detailed spatially resolved spectral analysis of DEM L71. We analyze regional spectra extracted from thin-sliced regions along several different azimuthal directions of the SNR to construct radial profiles of elemental abundances for O, Ne, Mg, Si, and Fe. Our elemental abundance measurements reveal an asymmetrical spatial distribution of metal-rich ejecta gas. Especially the asymmetry on the western part of the central Fe distribution is remarkable. While the location of the contact discontinuity is generally at $\sim 5$ pc from the geometric center of the X-ray emission of DEM L71, it is uncertain in western part of the remnant. Fe is enhanced in the ejecta while O and Ne abundances are generally negligible. This finding confirms the Type Ia origin of DEM L71. We estimate an upper limit on the Sedov age of $\sim 6,660\pm 770$ yr and explosion energy of $\sim 1.74\pm 0.35\times 10^{51}$ erg for the remnant. This explosion energy estimate is consistent with a canonical explosion of a Type Ia supernova remnant.

The James Webb Space Telescope will have the power to characterize high-redshift quasars at z>6 with an unprecedented depth and spatial resolution. While the brightest quasars at such redshift (i.e., with bolometric luminosity L_bol> 10^46 erg/s) provide us with key information on the most extreme objects in the Universe, measuring the black hole (BH) mass and Eddington ratios of fainter quasars with L_bol= 10^45-10^46 erg/s opens a path to understand the build-up of more normal BHs at z>6. In this paper, we show that the Illustris, TNG100, TNG300, Horizon-AGN, EAGLE, and SIMBA large-scale cosmological simulations do not agree on whether BHs at z>4 are overmassive or undermassive at fixed galaxy stellar mass with respect to the M_BH-M_star scaling relation at z=0 (BH mass offsets). Our conclusions are unchanged when using the local scaling relation produced by each simulation or empirical relations. We find that the BH mass offsets of the simulated faint quasar population at z>4, unlike those of bright quasars, represent the BH mass offsets of the entire BH population, for all the simulations. Thus, a population of faint quasars with L_bol= 10^45-10^46 erg/s observed by JWST can provide key constraints on the assembly of BHs at high redshift. Moreover, this will help constraining the high-redshift regime of cosmological simulations, including BH seeding, early growth, and co-evolution with the host galaxies. Our results also motivate the need for simulations of larger cosmological volumes down to z=6, with the same diversity of sub-grid physics, in order to gain statistics on the most extreme objects at high redshift.

We present a dedicated search for new pulsating helium-atmosphere (DBV) white dwarfs from the Sloan Digital Sky Survey using the McDonald 2.1m Otto Struve Telescope. In total we observed 55 DB and DBA white dwarfs with spectroscopic temperatures between 19,000 and 35,000K. We find 19 new DBVs and place upper limits on variability for the remaining 36 objects. In combination with previously known DBVs, we use these objects to provide an update to the empirical extent of the DB instability strip. With our sample of new DBVs, the red edge is better constrained, as we nearly double the number of DBVs known between 20,000 and 24,000K. We do not find any new DBVs hotter than PG 0112+104, the current hottest DBV at $T_{\mathrm{eff}}\,{\approx}$ 31,000K, but do find pulsations in four DBVs with temperatures between 27,000 and 30,000K, improving empirical constraints on the poorly defined blue edge. We investigate the ensemble pulsation properties of all currently known DBVs, finding that the weighted mean period and total pulsation power exhibit trends with effective temperature that are qualitatively similar to the pulsating hydrogen-atmosphere white dwarfs.

We present twelve epochs of optical spectroscopy taken across the discovery outburst of the black hole candidate MAXI J1803-298 with the GTC and VLT telescopes. The source followed a standard outburst evolution with hard and soft states. The system displays a triangular shape in the hardness intensity diagram, consistent with that seen in high inclination black hole transients and the previously reported detection of X-ray dips. The two epochs observed during the initial hard state exhibited asymmetric emission line profiles, including a P-Cygni profile simultaneously detected in H-alpha and He I 6678, which indicates the presence of an optical wind in the system. The remaining spectra, obtained during the transition to the soft state and the subsequent decay, are instead characterized by narrower, double peaked emission lines embedded into broad absorption components. One epoch (intermediate state) also includes near-infrared coverage, revealing complex line profiles in the Paschen and Bracket series, which suggests that the outflow is still present during the outburst decay through the soft state. The growing list of low-mass X-ray binaries with optical and near-infrared outflow signatures indicates that these are common features. Furthermore, the lowest luminosity spectrum exhibits an H-alpha full-width-at-half-maximum of 1570 +- 100 km/s. This, together with previous constraints on the binary parameters, allows us to favor a compact object mass of ~ 3-10 Msun, further supporting its black hole nature.

We have measured the wavelength-dependent lags between the X-ray, UV and optical bands in the high accretion rate ($L/L_{\rm Edd}\approx40\%$) Active Galactic Nucleus Mrk 110 during two intensive monitoring campaigns in February and September 2019. We divide the observations into three intervals with different X-ray luminosities. The first interval, already published in Vincentelli et al. (2021), has the lowest X-ray luminosity and did not exhibit the U-band excess positive lag, or the X-ray excess negative lag that is seen in most AGN. However, these excess lags are seen in the two subsequent intervals of higher X-ray luminosity. Although the data are limited, the excess lags appear to scale with X-ray luminosity. Our modelling shows that lags expected from reprocessing of X-rays by the accretion disc vary hardly at all with increasing luminosity. Therefore, as the U-band excess almost certainly arises from Balmer continuum emission from the broad line region (BLR), we attribute these lag changes to changes in the contribution from the BLR. The change is easily explained by the usual increase in the inner radius of the BLR with increasing ionising luminosity.

The observed exoplanet population features a gap in the radius distribution that separates the smaller super-Earths ($\lesssim$1.7$R_\oplus$) from the larger sub-Neptunes ($\sim$1.7--4$R_\oplus$). While mass loss theories can explain many of the observed features of this radius valley, it is difficult to reconcile them with a potentially rising population of terrestrials beyond orbital periods of ~30 days. We investigate the ability of initial gas accretion to reproduce both the location of the observed radius gap and the existence of long-period terrestrials. We first update the analytic scalings of gas accretion rate accounting for the shrinking of the bound radius by hydrodynamic effects. From the cooling evolution of planetary envelope with realistic opacity and equation of state, we find that the envelope mass fraction depends only weakly with the radius shrinking factor ($M_{\rm gas}/M_{\rm core} \propto f_R^{0.31}$). Co-evolving planetary masses and disk structures, focussing on dust-free opacity, we find that gas accretion alone is able to carve out the observed radius gap, with slopes $R_{\rm gap} \propto P^{-0.11}$ and $R_{\rm gap} \propto M_\star^{0.24}$ for top-heavy; and $R_{\rm gap} \propto P^{-0.10}$ and $R_{\rm gap} \propto M_\star^{0.21}$ for bottom-heavy core mass distributions, in good agreement with observations. Our model reconciles the location of the radius gap with the existence of long period terrestrials. The peaks and valleys in the radius distribution were likely set in place primordially while post-formation processes further tune the exoplanetary population. We provide potential observational tests that may be possible with PLATO and Roman Space Telescope.

Measuring the gas mass of protoplanetary disks, the reservoir available for giant planet formation, has proven to be difficult. We currently lack a far-infrared observatory capable of observing HD, and the most common gas mass tracer, CO, suffers from a poorly constrained CO-to-H$_2$ ratio. Expanding on previous work, we investigate if N2H+, a chemical tracer of CO poor gas, can be used to observationally measure the CO-to-H$_2$ ratio and correct CO-based gas masses. Using disk structures obtained from the literature, we set up thermochemical models for three disks, TW Hya, DM Tau and GM Aur, to examine how well the CO-to-H$_2$ ratio and gas mass can be measured from N2H+ and C18O line fluxes. Furthermore, we compare these gas masses to independently gas masses measured from archival HD observations. The N2H+ (3-2)/C18O (2-1) line ratio scales with the disk CO-to-H$_2$ ratio. Using these two lines, we measure 4.6e-3 Msun < Mdisk < 1.1e-1 Msun for TW Hya, 1.5e-2 Msun < Mdisk < 9.6e-2 Msun for GM Aur and 3.1e-2 Msun < Mdisk < 9.6e-2 Msun for DM Tau. These gas masses agree with values obtained from HD within their respective uncertainties. The uncertainty on the N2H+ + C18O gas mass can be reduced by observationally constraining the cosmic ray ionization rate in disks. These results demonstrate the potential of using the combination of N2H+ and C18O to measure gas masses of protoplanetary disks.

In this study, we investigated the orientation model of Broad Absorption Line (BAL) quasars using a sample of sources that are common in Sloan Digital Sky Survey (SDSS) Data Release (DR)-16 quasar catalog and Very Large Array (VLA)-Faint Images of the Radio Sky at Twenty Centimeters (FIRST) survey. Using the radio cut-out images from the FIRST survey, we first designed a deep learning model using convolutional neural networks (CNN) to classify the quasar radio morphologies into the core-only, young jet, single lobe, or triples. These radio morphologies are further sub-classified into core-dominated and lobe-dominated sources. The CNN models can classify the sources with a high precision of >98% for all the morphological sub-classes. The average BAL fraction in the resolved core, core-dominated, and lobe-dominated quasars are consistent with the BAL fraction inferred from radio and infra-red surveys. We also present the distribution of BAL quasars as a function of quasar orientation by using the radio core-dominance as an orientation indicator. A similar analysis is performed for HiBALs, LoBALs, and FeLoBALs. All the radio morphological sub-classes and BAL sub-classes show an increase in BAL fraction at high orientation angles of the jets with respect to the line of sight. Our analysis suggests that BAL quasars are more likely to be found in viewing angles close to the equatorial plane of the quasar. However, a pure orientation model is inadequate, and a combination of orientation and evolution is probably the best way to explain the complete BAL phenomena.

We combine deep optical and radio data, from the Hyper Suprime-Cam and the Low-Frequency Array (LOFAR) respectively, to study 78 radio AGN in nearby (z<0.5) dwarf galaxies. Comparison to a control sample, matched in stellar mass and redshift, indicates that the AGN and controls reside in similar environments, show similar star-formation rates (which trace gas availability) and exhibit a comparable incidence of tidal features (which indicate recent interactions). We explore the AGN properties by combining the predicted gas conditions in dwarfs from a cosmological hydrodynamical simulation with a Monte-Carlo suite of simulated radio sources, based on a semi-analytical model for radio-galaxy evolution. In the subset of LOFAR-detectable simulated sources, which have a similar distribution of radio luminosities as our observed AGN, the median jet powers, ages and accretion rates are $\sim 10^{35}$ W, $\sim 5$ Myr and $\sim 10^{-3.4}$ M$_{\odot}$ yr$^{-1}$ respectively. The median mechanical energy output of these sources is $\sim100$ times larger than the median binding energy expected in dwarf gas reservoirs, making AGN feedback plausible. Since special circumstances (in terms of environment, gas availability and interactions) are not necessary for the presence of AGN, and the central gas masses are predicted to be an order of magnitude larger than that required to fuel the AGN, AGN triggering in dwarfs is likely to be stochastic and a common phenomenon. Together with the plausibility of energetic feedback, this suggests that AGN could be important drivers of dwarf-galaxy evolution, as is the case in massive galaxies.

Study Analysis Group 21 (SAG21) of the Exoplanet Exploration Program Analysis Group (ExoPAG) was organized to study the effect of stellar contamination on space-based transmission spectroscopy, a method for studying exoplanetary atmospheres by measuring the wavelength-dependent radius of a planet as it transits its star. Transmission spectroscopy relies on a precise understanding of the spectrum of the star being occulted. However, stars are not homogeneous, constant light sources but have temporally evolving photospheres and chromospheres with inhomogeneities like spots, faculae, and plages. This SAG has brought together an interdisciplinary team of more than 100 scientists, with observers and theorists from the heliophysics, stellar astrophysics, planetary science, and exoplanetary atmosphere research communities, to study the current needs that can be addressed in this context to make the most of transit studies from current NASA facilities like HST and JWST. The analysis produced 14 findings, which fall into three Science Themes encompassing (1) how the Sun is used as our best laboratory to calibrate our understanding of stellar heterogeneities ("The Sun as the Stellar Benchmark"), (2) how stars other than the Sun extend our knowledge of heterogeneities ("Surface Heterogeneities of Other Stars") and (3) how to incorporate information gathered for the Sun and other stars into transit studies ("Mapping Stellar Knowledge to Transit Studies").

The symbiotic X-ray binary Sct X-1 was suggested as the first known neutron star accreting from a red supergiant companion. Although known for nearly 50 years, detailed characterization of the donor remains lacking, particularly due to the extremely high reddening towards the source ($A_V\gtrsim25$ mag). Here, we present i) improved localization of the counterpart using Gaia and Chandra observations, ii) the first broadband infrared spectrum ($\approx1-5\,\mu$m; $R\approx 2000$) obtained with SpeX on the NASA Infrared Telescope Facility and iii) $J$-band light curve from the Palomar Gattini-IR survey. The infrared spectrum is characterized by i) deep water absorption features (H$_2$O index $\approx 40$%), ii) strong TiO, VO and CO features, and iii) weak/absent CN lines. We show that these features are inconsistent with known red supergiants, but suggest a M8-9 III type O-rich Mira donor star. We report the discovery of large amplitude ($\Delta J\approx3.5$ mag) periodic photometric variability suggesting a pulsation period of $621\pm36\,{\rm(systematic)}\pm8\,{\rm(statistical)}$ days, which we use to constrain the donor to be a relatively luminous Mira ($M_K=-8.6\pm0.3$ mag) at a distance of $3.6^{+0.8}_{-0.7}$ kpc. Comparing these characteristics to recent models, we find the donor to be consistent with a $\approx 3-5$ M$_\odot$ star at an age of $\approx 0.1-0.3$ Gyr. Together, we show that Sct X-1 was previously mis-classified as an evolved High Mass X-ray Binary; instead it is an intermediate mass system with the first confirmed Mira donor in an X-ray binary. We discuss the implications of Mira donors in symbiotic X-ray binaries, and highlight the potential of wide-field infrared time domain surveys and broadband infrared spectroscopy to unveil their demographics.

The magnificent bubbles at the Galactic center provide a great channel to understand the effects of feedback on galaxy evolution. The newly discovered eROSITA bubbles show enhanced X-ray emission from the shells around bubbles. Previous works assumed that the X-ray emitting gas in the shells has a single temperature component and that they trace the shock-heated lower-temperature Galactic halo gas. Here we show that the thermal structure of the eROSITA bubble shells is more complex. Using Suzaku observations we find with high confidence that the X-ray emission from the shells is best described by a two-temperature thermal model, one near Galaxy's virial temperature at $\rm kT \approx 0.2 ~keV$ and the other at super-virial temperatures ranging between $\rm kT = 0.4-1.1 ~keV$. Furthermore, we show that temperatures of the virial and super-virial components are similar in the shells and in the ambient medium, although the emission measures are significantly higher in the shells. We argue that the X-ray bright eROSITA bubble shells are the signature of compressed isothermal radiative shocks. The age of the bubbles is constrained to $70$--$\rm 130~Myr$. This expansion timescale, as well as the observed non-solar Ne/O and Mg/O ratios, favor the stellar feedback models for the formation of the Galactic bubbles, settling a long-standing debate on the origin of the Galactic bubbles.

We confirm the existence of a second Large Magellanic Cloud (LMC) star cluster, KMHK 1592, with an age that falls in the middle of the so-called LMC star cluster age gap, a long period of time (~ 4 - 11 Gyr) where no star cluster had been uncovered, except ESO121-SC03. The age (8.0+-0.5 Gyr) and the metallicity ([Fe/H]=-1.0+-0.2 dex) of KMHK1592 were derived from the fit of theoretical isochrones to the intrinsic star cluster colour-magnitude diagram sequences, which were unveiled using a robust star-by-star membership probability procedure. Because of the relative low brightness of the star cluster, deep GEMINI GMOS images were used. We discuss the pros and cons of three glimpsed scenarios that could explain the presence of both LMC age gap star clusters in the outskirts of the LMC, namely: in-situ star cluster formation, capture from the Small Magellanic Cloud, or accretion of a small dwarf galaxy.

The Large Magellanic Cloud (LMC) is the nearest laboratory for detailed studies on the formation and survival of complex organic molecules (COMs), including biologically important ones, in low-metallicity environments--typical for earlier cosmological epochs. We report the results of 1.2 mm continuum and molecular line observations of three fields in the star-forming region N105 with the Atacama Large Millimeter/submillimeter Array (ALMA). N105 lies at the western edge of the LMC bar with on-going star formation traced by H$_2$O, OH, and CH$_3$OH masers, ultracompact H II regions, and young stellar objects. Based on the spectral line modeling, we estimated rotational temperatures, column densities, and fractional molecular abundances for twelve 1.2 mm continuum sources. We identified sources with a range of chemical make-ups, including two bona fide hot cores and four hot core candidates. The CH$_3$OH emission is widespread and associated with all the continuum sources. COMs CH$_3$CN and CH$_3$OCH$_3$ are detected toward two hot cores in N105 together with smaller molecules typically found in Galactic hot cores (e.g., SO$_2$, SO, and HNCO) with the molecular abundances roughly scaling with metallicity. We report a tentative detection of the astrobiologically relevant formamide molecule (NH$_2$CHO) toward one of the hot cores; if confirmed, this would be the first detection of NH$_2$CHO in an extragalactic sub-solar metallicity environment. We suggest that metallicity inhomogeneities resulting from the tidal interactions between the LMC and the Small Magellanic Cloud (SMC) might have led to the observed large variations in COM abundances in LMC hot cores.

We study the formation of star clusters in molecular clouds by performing three-dimensional radiation hydrodynamics simulations with far ultraviolet (FUV; $6 ~{\rm eV} \leqq h \nu \leqq 13.6 ~{\rm eV}$) and extreme ultraviolet (EUV; $h\nu \geqq 13.6~{\rm eV}$) radiative feedback. We find that the FUV feedback significantly suppresses the star formation in diffuse clouds with the initial surface densities of $\Sigma_{\rm cl} \lesssim \rm 50~M_{\odot} \; pc^{-2}$. In the cases of clouds with $\Sigma_{\rm cl} \sim \rm 100-200~M_{\odot} \; pc^{-2}$, the EUV feedback plays a main role and decrease the star formation efficiencies less than $0.3$. We show that thermal pressure from PDRs or H{\sc ii} regions disrupts the clouds and makes the size of the star clusters larger. Consequently, the clouds with the mass $M_{\rm cl} \lesssim 10^{5}~\rm M_{\odot}$ and the surface density $\Sigma_{\rm cl} \lesssim 200~\rm M_{\odot}\; pc^{-2}$ remain the star clusters with the stellar densities of $\sim 100~\rm M_{\odot}\; pc^{-3}$ that nicely match the observed open clusters in the Milky Way. If the molecular clouds are massive ($M_{\rm cl} \gtrsim 10^{5}~\rm M_{\odot}$) and compact ($\Sigma \gtrsim 400~\rm M_{\odot}\; pc^{-2}$), the radiative feedback is not effective and they form massive dense cluster with the stellar densities of $\sim 10^{4}~\rm M_{\odot}\; pc^{-3}$ like observed globular clusters or young massive star clusters. Thus, we suggest that the radiative feedback and the initial conditions of molecular clouds are key factors inducing the variety of the observed star clusters.

The Large Magellanic Cloud (LMC) is the most luminous satellite galaxy of the Milky Way and owing to its companion, the Small Magellanic Cloud (SMC), represents an excellent laboratory to study the interaction of dwarf galaxies. The aim of this study is to investigate the kinematics of the outer regions of the LMC by using stellar proper motions to understand the impact of interactions, e.g. with the SMC about 250 Myr ago. {We calculate proper motions using multi-epoch $K_\mathrm{s}$-band images from the VISTA survey of the Magellanic Clouds system (VMC). Observations span a time baseline of 2$-$5 yr. We combine the VMC data with data from the \textit{Gaia} early Data Release 3 and introduce a new method to distinguish between Magellanic and Milky Way stars based on a machine learning algorithm. This new technique enables a larger and cleaner sample selection of fainter sources as it reaches below the red clump of the LMC. We investigate the impact of the SMC on the rotational field of the LMC and find hints of stripped SMC debris. The south east region of the LMC shows a slow rotational speed compared to the overall rotation. $N$-body simulations suggest that this could be caused by a fraction of stripped SMC stars, located in that particular region, that move opposite to the expected rotation.

Hot cores characterized by rich lines of complex organic molecules are considered as ideal sites for investigating the physical and chemical environments of massive star formation. We present a search for hot cores by using typical nitrogen- and oxygen-bearing complex organic molecules (C$_2$H$_5$CN, CH$_3$OCHO and CH$_3$OH), based on ALMA Three-millimeter Observations of Massive Star-forming regions (ATOMS). The angular resolutions and line sensitivities of the ALMA observations are better than 2 arcsec and 10 mJy/beam, respectively. A total of 60 hot cores are identified with 45 being newly detected, in which the complex organic molecules have high gas temperatures ($>$ 100 K) and small source sizes ($<$ 0.1 pc). So far this is the largest sample of hot cores observed with similar angular resolution and spectral coverage. The observations have also shown nitrogen and oxygen differentiation in both line emission and gas distribution in 29 hot cores. Column densities of CH$_3$OH and CH$_3$OCHO increase as rotation temperatures rise. The column density of CH$_3$OCHO correlates tightly with that of CH$_3$OH. The pathways for production of different species are discussed. Based on the spatial position difference between hot cores and UC~H{\sc ii} regions, we conclude that 24 hot cores are externally heated while the other hot cores are internally heated. The observations presented here will potentially help establish a hot core template for studying massive star formation and astrochemistry.

We present recent updates and improvements of the graphical processing unit (GPU) N-body code GENGA. Modern state-of-the-art simulations of planet formation require the use of a very high number of particles to accurately resolve planetary growth and to quantify the effect of dynamical friction. At present the practical upper limit is in the range of 30'000 - 60'000 fully interactive particles, with possibly a little more on the latest GPU devices. While the original hybrid symplectic integration method has difficulties to scale up to these numbers, we have improved the integration method by i) introducing higher level changeover functions, and ii) the code is better able to use the most recent GPU hardware efficiently for such large simulations. We added treatments of non-Newtonian forces such as general relativity, tidal interaction, rotational deformation, the Yarkovsky effect and Poynting-Robertson drag as well as a new model to treat virtual collisions of small bodies in the Solar System. We added new tools to GENGA, such as semi-active test particles that feel more massive bodies but not each other, a more accurate collision handling and a real-time openGL-visualization. We present example simulations, including a 1.5 billion year terrestrial planet formation simulation that initially started with 65536 particles, a 3.5 billion year simulation without gas giants starting with 32768 particles, the evolution of asteroid fragments in the Solar System, and the planetesimal accretion of a growing Jupiter simulation. GENGA runs on modern NVIDIA and AMD GPUs.

We use data from the first six encounters of Parker Solar Probe and employ the Partial Variance of Increments ($PVI$) method to study the statistical properties of coherent structures in the inner heliosphere with the aim of exploring physical connections between magnetic field intermittency and observable consequences such as plasma heating and turbulence dissipation. Our results support proton heating localized in the vicinity of, and strongly correlated with, magnetic structures characterized by $PVI \geq 1$. We show that on average, such events constitute $\approx 19\%$ of the dataset, though variations may occur depending on the plasma parameters. We show that the waiting time distribution ($WT$) of identified events is consistent across all six encounters following a power-law scaling at lower $WTs$. This result indicates that coherent structures are not evenly distributed in the solar wind but rather tend to be tightly correlated and form clusters. We observe that the strongest magnetic discontinuities, $PVI \geq 6$, usually associated with reconnection exhausts, are sites where magnetic energy is locally dissipated in proton heating and are associated with the most abrupt changes in proton temperature. However, due to the scarcity of such events, their relative contribution to energy dissipation is minor. Taking clustering effects into consideration, we show that smaller scale, more frequent structures with PVI between, $1\lesssim PVI \lesssim 6$, play the major role in magnetic energy dissipation. The number density of such events is strongly associated with the global solar wind temperature, with denser intervals being associated with higher $T_{p}$.

In this work we present relensing, a package written in python whose goal is to model galaxy clusters from gravitational lensing. With relensing we extent the amount of software available, which provides the scientific community with a wide range of models that help to compare and therefore validate the physical results that rely on them. We implement a free-form approach which computes the gravitational deflection potential on a adaptive irregular grid, from which one can characterize the cluster and its properties as a gravitational lens. Here, we use two alternative penalty functions to constrain strong lensing. We apply relensing to two toy models, in order to explore under which conditions one can get a better performance in the reconstruction. We find that by applying a smoothing to the deflection potential enhances the capability of this approach to recover the shape and size of the galaxy cluster's mass profile, as well as its magnification map, which translates in a better estimation of the critical and caustic curves. The power that the smoothing provides is also tested on the simulated clusters Ares and Hera, for which our results represent an improvement with respect to reconstructions that were carried out with methods of the same nature than relensing. At the same time, the smoothing also increases the stability of our implementation, and decreases the computation time. On its current state, relensing is available upon request.

The radiation of a Fast Radio Burst (FRB) reflects from the Moon and Sun. If a reflection is detected, the time interval between the direct and reflected signals constrains the source to a narrow arc on the sky. If both Lunar and Solar reflections are detected these two arcs intersect, narrowly confining the source location on the sky. Galactic FRB like FRB 200428 may be bright enough to be detected by a 25 m diameter radio telescope staring at the Moon or Sun. A previous paper calculated reflection by the Moon. Here we calculate the reflectivity of the Sun in the "flat Sun" approximation as a function of angle of incidence and frequency. At grazing incidence the reflectivity is high at frequencies $\lessapprox 200\,$MHz but low at higher frequencies; for near-normal incidence the reflectivity is high only for frequencies $\lessapprox 100\,$MHz.

Context. Classical novae are powered by thermonuclear runaways occurring on the surface of accreting white dwarfs (WDs). In the observations, the enrichments of heavy elements in nova ejecta have been detected, indicating the existence of a mixing process between the accreted matter and the matter from the outer layers of the underlying WDs prior to nova outbursts. However, the mixing fraction in classical novae is still uncertain. Aims. The purpose of this article is to investigate some elemental abundance ratios during nova outbursts, which can be used to estimate the WD-mixing fraction in classical novae. Methods. By considering different WD-mixing fractions with the stellar evolution code Modules for Experiments in Stellar Astrophysics (MESA), we carried out a series of simulations of nova outbursts, in which the initial CO WD masses range from $0.7-1.0\,M_\odot$. Results. We identified four elemental abundance ratios (i.e. $\rm (H+He)/\sum CNO$, $\rm (H+He)/Ne$, $\rm \sum CNO/Mg$ and $\rm \sum CNO/Si$) that satisfy the conditions for determining the WD-mixing fraction, in which $\rm (H+He)/\sum CNO$ is the most suitable mixing meter. We also estimated the WD-mixing fraction in some representative classical novae. Additionally, we found that a higher metallicity (i.e. higher WD-mixing fraction) prefers to be accompanied by a longer $t_{\rm 2}$ (the time of decline by two magnitudes from peak luminosity) during nova outbursts. The present results can be used to constrain the mixing process in classical novae.

The fluctuating Gunn-Peterson approximation (FGPA) is a commonly-used method to generate mock Lyman-$\alpha$ (Ly$\alpha$) forest absorption skewers at Cosmic Noon ($z\gtrsim 2$) from the matter-density field of $N$-body simulations without running expensive hydrodynamical simulations. Motivated by recent developments in 3D IGM tomography observations as well as matter density field reconstruction techniques applied to galaxy redshift samples at $z\sim 2$, we examine the possibility of observationally testing FGPA by directly examining the relationship between the Ly$\alpha$ transmission and the underlying matter density field. Specifically, we analyze the EAGLE, Illustris, IllustrisTNG and Nyx cosmological hydrodynamic simulations, that were run with different codes and sub-grid models. While FGPA is an excellent description of the IGM in lower-density regions, the slope of the transmission-density distribution at higher densities is significantly affected by feedback processes causing FGPA to break down in that regime. Even without added feedback, we find significant deviations caused by hydrodynamical effects arising from non-linear structure growth. We then proceed to make comparisons using realistic mock data assuming the sightline sampling and spectral properties of the recent CLAMATO survey, and find that it would be challenging to discern between FGPA and hydrodynamical models with current data sets. However, the improved sightline sampling from future extremely large telescopes or large volumes from multiplexed spectroscopic surveys such as Subaru PFS should allow for stringent tests of FGPA, and make it possible to detect the effect of galaxy feedback on the IGM.

The next generation of telescopes such as the SKA and the Rubin Observatory will produce enormous data sets, requiring automated anomaly detection to enable scientific discovery. Here, we present an overview and friendly user guide to the Astronomaly framework for active anomaly detection in astronomical data. Astronomaly uses active learning to combine the raw processing power of machine learning with the intuition and experience of a human user, enabling personalised recommendations of interesting anomalies. It makes use of a Python backend to perform data processing, feature extraction and machine learning to detect anomalous objects; and a JavaScript frontend to allow interaction with the data, labelling of interesting anomalous and active learning. Astronomaly is designed to be modular, extendable and run on almost any type of astronomical data. In this paper, we detail the structure of the Astronomaly code and provide guidelines for basic usage.

In the search for life in the Universe, exoplanets represent numerous natural experiments in planet formation, evolution, and the emergence of life. This raises the fascinating prospect of evaluating cosmic life on a statistical basis. One key statistic is the occurrence rate of life-bearing worlds, $f_{\rm L}$, the 'frequency of life' term in the famous Drake Equation. Measuring $f_{\rm L}$ would give profound insight into how common life is and may help to constrain origin-of-life theories. I propose $f_{\rm L}$ as the goal for the DRAKE mission (Dedicated Research for Advancing Knowledge of Exobiology): a transit spectroscopy survey of M-dwarf habitable zone terrestrial planets. I investigate how the uncertainty on the observed value of $f_{\rm L}$ scales with sample size. I determine that sampling error dominates over observational error and that the uncertainty is a function of the observed $f_{\rm L}$ value. I show that even small sample sizes can provide significant constraints on $f_{\rm L}$, boding well for the transit spectroscopy approach. I perform a feasibility study of the DRAKE mission using a nominal instrument design and mission plan. Due to low observing efficiencies, DRAKE may need to be incorporated into a wider-ranging deep-space or lunar observatory. A 50-planet survey could constrain $f_{\rm L}$ to $\leq$ 0.06 (at 95% confidence) if the sample $f_{\rm L}$ = 0, or 0.03-0.2 if the sample $f_{\rm L}$ = 0.1. This can be achieved (on average) in 10 years using a 17-m telescope with an unrestricted field-of-regard. DRAKE is a viable approach to attempting the first experimental measurement of $f_{\rm L}$.

Ultralight axions are theoretically interesting and phenomenologically rich dark sector candidates, but they are difficult to track across cosmological timescales because of their fast oscillations. We resolve this problem by developing a novel method to evolve them efficiently and accurately. We first construct an exact effective fluid which at late times matches the axion but which evolves in a simple way. We then approximate this evolution with a carefully chosen equation of state and sound speed. With our scheme we find that we can obtain subpercent accuracy for the linear theory suppression of axion density fluctuations relative to that of cold dark matter without tracking even a single complete oscillation of the axion field. We use our technique to test other approximation schemes and to provide a fitting formula for the transfer function for the matter power spectrum in linear theory in axion models. Implementing our approach in existing cosmological axion codes is straightforward and will help unleash the potential of high-precision next-generation experiments.

Detection of the gravitational-wave events revealed that there are numerous population of the black hole binaries which can merge within the age of the Universe. Although several formation channels of such binaries are known, considerable theoretical uncertainties associated in each channel defeats the robust prediction of how much each channel contributes to the total merger rate density. Given that the time evolution of the merger rate density in some channels is (exactly or nearly) independent of the BH masses, clarifying this feature from the observational data will shed some light on the nature of the black hole binaries. Based on this motivation, we formulate the methodology to perform the statistical test of whether the mass distribution of the black hole mergers evolves in time or not by means of the hypothesis testing. Our statistical test requires neither a priori specification of the mass distribution which is largely uncertain nor that of the time dependence of the merger rate. We then apply it to the mock data for some concrete shapes of the merger rate density and show that the proposed method rejects/(does not reject) the null hypothesis correctly for the large sample size. We also investigate if the catalog of the gravitational-wave events obtained during the LIGO-Virgo's third observing run has a large sample size enough to apply our hypothesis testing. We find that the number of the events is too small to draw any statistical conclusion regarding our test and the meaningful result of our hypothesis testing can be obtained only by the future detectors having much better sensitivity. These results demonstrate the effectiveness of our hypothesis testing to determine from the (future) observational data whether the merger rate density evolves in time independently of the BH masses or not.

INTEGRAL/SPI is a coded mask instrument observing since 2002 in the keV to MeV energy range, which covers the peak of the $\nu F\nu$ spectrum of most Gamma-Ray Bursts (GRBs). Since its launch in 2008, Fermi/GBM has been the primary instrument for analyzing GRBs in the energy range between $\approx$ 10 keV to $\approx$ 10 MeV. Herein, we show that SPI, covering a similar energy range, can give equivalently constraining results for some parameters if we use an advanced analysis method. Also, combining the data of both instruments reduces the allowed parameter space in spectral fits. The main advantage of SPI as compared to GBM is the energy resolution of $\approx$ 0.2\% at 1.3 MeV compared to $\approx$ 10\% for GBM. Therefore, SPI is an ideal instrument to precisely measure the curvature of the spectrum. This is important, as it has been shown in recent years that physical models rather than heuristic functions should be fit to GRB data to obtain better insights into their still unknown emission mechanism, and the curvature of the peak is unique to the different physical models. To fit physical models to SPI GRB data and get the maximal amount of information from the data, we developed a new open source analysis software {\tt PySPI}. We apply these new techniques to GRB 120711A in order to validate and showcase {\tt PySPI}'s capabilities. We show that {\tt PySPI} improves the analysis of SPI GRB data compared to the {\tt OSA} analysis. In addition, we demonstrate that the GBM and the SPI data of this GRB can be fitted well with a physical synchrotron model. This evinces that SPI can play an important role in GRB spectral model fitting.

The first interstellar object to be observed in our solar system 1I/2017 U1 'Oumuamua combines the lack of observable cometary activity with an extra-gravitational acceleration. This has given rise to several mutually exclusive explanations based on different assumptions in the material composition of 'Oumuamua. We show how a combination of observations in the infrared and optical spectra may serve to distinguish between these explanations once another object with 'Omuamua-like properties comes close enough to earth. This possibility is linked to the widely different thermal properties of the different material models that have been proposed. Developing a model for the thermal conduction and infrared signal from a fractal model we compare predictions of the infrared signal with that from standard thermal models that assume 'Oumuamua to be either a solid piece of rock/ice or a thin sheet.

Detecting molecular line emission from classical nova remnants has the potential of revealing information on the composition of the ejecta, in particular, it can deliver accurate isotopic ratios in the matter processed by a thermonuclear runaway. We conducted searches toward more than 100 classical novae for emission in lines of the CO or HCN molecules using single-dish telescopes and interferometric arrays at millimeter and submillimeter wavelengths. The survey demonstrates that classical novae, young or old, are not strong sources of molecular emission at submillimeter and millimeter wavelengths. Additionally, we mapped CO emission around Nova Persei 1901 (GK Per), earlier claimed to be circumstellar in origin. Our measurements indicate that the observed emission is from the interstellar medium. Although no molecular emission at millimeter and submillimeter wavelengths has been found in classical novae, it is still likely that some will be detected with high-sensitivity interferometers such as ALMA.

Context. To better understand the formation of high-mass stars, it is fundamental to investigate how matter accretes onto young massive stars, how it is ejected, and how all this differs from the low-mass case. The massive protocluster G31.41+0.31 is the ideal target to study all these processes because observations at millimeter and centimeter wavelengths have resolved the emission of the Main core into at least four massive dust continuum sources, named A, B, C, and D, within 1" or 0.018 pc, and have identified signatures of infall and several outflows associated with the core. Aims. We study the interplay between infall and outflow in G31.41+0.31 by investigating at a spatial resolution of a few 100 au their properties and their possible impact on the core. Methods. We carried out molecular line observations of typical high-density tracers, such as CH3CN or H2CO, and shock and outflow tracers, such as SiO, with ALMA at 1.4 mm that achieved an angular resolution of 0.09" (340 au). Results. The observations have revealed inverse P-Cygni profiles in CH3CN and H2CO toward the four sources embedded in the Main core, suggesting that all of them are undergoing collapse. The infall rates, estimated from the red-shifted absorption are on the order of 1E-2 Msun/yr. The individual infall rates imply that the accretion timescale of the Main core is an order of magnitude smaller than its rotation timescale. This confirms that rotating toroids such as the G31 Main core are non-equilibrium, transient collapsing structures that need to be constantly replenished with fresh material from a large-scale reservoir. For sources B, C, and D, the infall could be accelerating inside the sources, while for source A, the presence of a second emission component complicates the interpretation. The SiO observations have revealed the presence of at least six outflows in the G31.41+0.31 star-forming region, ...

A stellar occultation by Neptune's main satellite, Triton, was observed on 5 October 2017 from Europe, North Africa, and the USA. We derived 90 light curves from this event, 42 of which yielded a central flash detection. We aimed at constraining Triton's atmospheric structure and the seasonal variations of its atmospheric pressure since the Voyager 2 epoch (1989). We also derived the shape of the lower atmosphere from central flash analysis. We used Abel inversions and direct ray-tracing code to provide the density, pressure, and temperature profiles in the altitude range $\sim$8 km to $\sim$190 km, corresponding to pressure levels from 9 {\mu}bar down to a few nanobars. Results. (i) A pressure of 1.18$\pm$0.03 {\mu}bar is found at a reference radius of 1400 km (47 km altitude). (ii) A new analysis of the Voyager 2 radio science occultation shows that this is consistent with an extrapolation of pressure down to the surface pressure obtained in 1989. (iii) A survey of occultations obtained between 1989 and 2017 suggests that an enhancement in surface pressure as reported during the 1990s might be real, but debatable, due to very few high S/N light curves and data accessible for reanalysis. The volatile transport model analysed supports a moderate increase in surface pressure, with a maximum value around 2005-2015 no higher than 23 {\mu}bar. The pressures observed in 1995-1997 and 2017 appear mutually inconsistent with the volatile transport model presented here. (iv) The central flash structure does not show evidence of an atmospheric distortion. We find an upper limit of 0.0011 for the apparent oblateness of the atmosphere near the 8 km altitude.

One of the guaranteed features of the stochastic gravitational wave background (SGWB) is the presence of Doppler anisotropies induced by the motion of the detector with respect to the rest frame of the SGWB source. We point out that kinematic effects can be amplified if the SGWB is characterised by large tilts in its spectrum as a function of frequency, or by sizeable intrinsic anisotropies. Hence we examine the possibility to use Doppler effects as complementary probes of the SGWB frequency profile. For this purpose we work in multipole space, and we study the effect of kinematic modulation and aberration on the GW energy density parameter and on its angular power spectrum. We develop a Fisher forecast analysis and we discuss prospects for constraining parameters controlling kinematically induced anisotropies with future detector networks. As a case study, we apply our framework to a background component with constant slope in frequency, potentially detectable by a network of future ground-based interferometers. For this specific example, we show that a measurement of kinematic anisotropies with a network of Einstein Telescope and Cosmic Explorer will allow us to constrain the spectral shape with a precision of about 16$\%$. Finally, we identify cosmological and astrophysical scenarios where kinematic effects are enhanced in frequency ranges probed by current and future GW experiments.

We present the results of 15 years of VISNIR spectral monitoring studying the famously variable RW Aur A classical T Tauri star (CTTS) system. We find direct evidence for a highly excited classical T Tauri star surrounded by an IR bright, asymmetric, and time variable accretion disk. Comparison of the spectral and temporal trends found in long term SpeX monitoring determines 5 different components: (1) a stable continuum from 0.7 - 1.3 um, with color temperature 4000K, produced by the CTTS photospheric surface; (2) variable hydrogen emission lines emitted from hot excited hydrogen in the CTTS protostellar atmosphere-accretion envelope; (3) hot CO gas in the CTTS protostellar atmosphere/accretion envelope; (4) highly variable 1.8-5.0 um thermal continuum emission with color temperature ranging from 1130 to 1650K, due to a surrounding accretion disk that is spatially variable and has an inner wall at r = 0.04 AU and T=1650K, and outer edges at 1130K; and (5) transient signatures of abundant Fe II + associated SI, SiI, and SrI in the systems jet structures created by the catastrophic disruption of a planetesimal core. The Fe II signatures first appeared in 2015 as highly bifurcated jet lines, but these collapsed and disappeared into a small single peak protostellar atmosphere feature by late 2020. By contrast, nearby, coeval binary companion RW Aur B evinces only (1) a stable WTTS photospheric continuum from 0.7 - 1.3 um + (3) cold CO gas in absorption + (4) stable 1.8-5.0 um thermal continuum emission with color temperature 1650K. The temporal evolution of the RW Aur A spectral signatures we find is consistent with a dynamically excited CTTS system forming differentiated Vesta-sized planetesimals in an asymmetric accretion disk and migrating them inward to be destructively accreted.

We investigate how different mid-infrared (mid-IR) properties of galaxies trace the environment in which the galaxies are located. For this purpose, we first study the dependence of galaxy clustering on the absolute magnitude at 3.4 $\mu$m and redshift. Then, we look into the correlation between various mid-IR properties and the environment. We also explore how various infrared galaxy luminosity selections influence the galaxy clustering measurements. We use a set of W1 (3.4 $\mu$m) absolute magnitude ($M_\text{W1}$) selected samples from the Galaxy and Mass Assembly (GAMA) survey matched with mid-IR properties from the Wide-field Infrared Survey Explorer (WISE) in the redshift range $0.07 \leq z < 0.43$. We compute the galaxy two-point correlation function (2pCF) and compare the clustering lengths between subsamples binned in $M_\text{W1}$ and in redshift. We also measure the marked correlation functions (MCFs) using the luminosities in the WISE W1 to W4 (3.4 to 22 $\mu$m) bands. We also measure MCFs with different estimates of stellar mass and star formation rate used as marks and compare the results. Finally, we check how different selections applied to the sample affect the clustering measurements. We show strong clustering dependence on the W1 absolute magnitude: galaxies brighter in the W1 band are more strongly clustered than their fainter counterparts. We also observe a lack of significant redshift dependence of clustering in the redshift range $0.07 \leq z < 0.43$. We show that W1 and W2 bands closely follow the stellar mass and W3 and W4 bands closely follow the star formation rate in tracing the galaxy clustering. We demonstrate the agreement between different estimates of the stellar mass and the star formation rate in tracing the environment.

Stellar intensity interferometers correlate photons within their coherence time and could overcome the baseline limitations of existing amplitude interferometers. Intensity interferometers do not rely on phase coherence of the optical elements and thus function without high grade optics and light combining delay lines. However, the coherence time of starlight observed with realistic optical filter bandwidths (> 0.1 nm) is usually much smaller than the time resolution of the detection system (> 10 ps), resulting in a greatly reduced correlation signal. Reaching high signal to noise in a reasonably short measurement time can be achieved in different ways: either by increasing the time resolution, which increases the correlation signal height, or by increasing the photon rate, which decreases statistical uncertainties of the measurement. We present laboratory measurements employing both approaches and directly compare them in terms of signal to noise ratio. A high time-resolution interferometry setup designed for small to intermediate size optical telescopes and thus lower photon rates (diameters < some meters) is compared to a setup capable of measuring high photon rates, which is planned to be installed at Cherenkov telescopes with dish diameters of > 10 m. We use a Xenon lamp as a common light source simulating starlight. Both setups measure the expected correlation signal and work at the expected shot-noise limit of statistical uncertainties for measurement times between 10 min and 23 h. We discuss the quantitative differences in the measurement results and give an overview of suitable operation regimes for each of the interferometer concepts.

Population synthesis relies on semi-analytic formulae to determine masses of compact objects from the (helium or carbon-oxygen) cores of collapsing stars. Such formulae are combined across mass ranges that span different explosion mechanisms, potentially introducing artificial features in the compact object mass distribution. Such artifacts impair the interpretation of gravitational-wave observations. We propose a "top-down" remnant mass prescription where we remove mass from the star for each possible mass-loss mechanism, instead of relying on the fallback onto a "proto-compact-object" to get the final mass. For one of these mass-loss mechanisms, we fit the metallicity-dependent mass lost to pulsational-pair instability supernovae from numerical simulations. By imposing no mass loss in the absence of pulses, our approach recovers the existing compact object masses prescription at the low mass end and ensures continuity across the core-collapse/pulsational-pair-instability regime. Our remnant mass prescription can be extended to include other mass-loss mechanisms at the final collapse.

Neutrinos from supernova (SN) bursts can give rise to detectable number of nuclear recoil (NR) events through the coherent elastic neutrino-nucleus scattering (CE$\nu$NS) process in large scale liquid xenon detectors designed for direct dark matter search, depending on the SN progenitor mass and distance. Here we show that in addition to the direct NR events due to CE$\nu$NS process, the SN neutrinos can give rise to additional nuclear recoils due to the elastic scattering of neutrons produced through inelastic interaction of the neutrinos with the xenon nuclei. We find that the contribution of the supernova neutrino-induced neutrons ($\nu$I$n$) can significantly modify the total xenon NR spectrum at large recoil energies compared to that expected from the CE$\nu$NS process alone. Moreover, for recoil energies $\gtrsim20$ keV, dominant contribution is obtained from the ($\nu$I$n$) events. We numerically calculate the observable S1 and S2 signals due to both CE$\nu$NS and $\nu$I$n$ processes for a typical liquid xenon based detector, accounting for the multiple scattering effects of the neutrons in the case of $\nu$I$n$, and find that sufficiently large signal events, those with S1$\gtrsim$50 photo-electrons (PE) and S2$\gtrsim$2300 PE, come mainly from the $\nu$I$n$ scatterings.

The inner $\sim$200 pc region of the Galaxy contains a 4 million M$_{\odot}$ supermassive black hole (SMBH), significant quantities of molecular gas, and star formation and cosmic ray energy densities that are roughly two orders of magnitude higher than the corresponding levels in the Galactic disk. At a distance of only 8.2 kpc, the region presents astronomers with a unique opportunity to study a diverse range of energetic astrophysical phenomena, from stellar objects in extreme environments, to the SMBH and star-formation driven feedback processes that are known to influence the evolution of galaxies as a whole. We present a new survey of the Galactic center conducted with the South African MeerKAT radio telescope. Radio imaging offers a view that is unaffected by the large quantities of dust that obscure the region at other wavelengths, and a scene of striking complexity is revealed. We produce total intensity and spectral index mosaics of the region from 20 pointings (144 hours on-target in total), covering 6.5 square degrees with an angular resolution of 4$''$,at a central frequency of 1.28 GHz. Many new features are revealed for the first time due to a combination of MeerKAT's high sensitivity, exceptional $u,v$-plane coverage, and geographical vantage point. We highlight some initial survey results, including new supernova remnant candidates, many new non-thermal filament complexes, and enhanced views of the Radio Arc Bubble, Sgr A and Sgr B regions. This project is a SARAO public legacy survey, and the image products are made available with this article.

Red supergiants are the most common final evolutionary stage of stars that have initial masses between 8 and 35 times that of the Sun. During this stage, which lasts roughly 100,000 years1, red supergiants experience substantial mass loss. However, the mechanism for this mass loss is unknown. Mass loss may affect the evolutionary path, collapse and future supernova light curve of a red supergiant, and its ultimate fate as either a neutron star or a black hole. From November 2019 to March 2020, Betelgeuse - the second-closest red supergiant to Earth (roughly 220 parsecs, or 724 light years, away) - experienced a historic dimming of its visible brightness. Usually having an apparent magnitude between 0.1 and 1.0, its visual brightness decreased to 1.614 +/- 0.008 magnitudes around 7-13 February 2020 - an event referred to as Betelgeuse's Great Dimming. Here we report high-angular-resolution observations showing that the southern hemisphere of Betelgeuse was ten times darker than usual in the visible spectrum during its Great Dimming. Observations and modelling support a scenario in which a dust clump formed recently in the vicinity of the star, owing to a local temperature decrease in a cool patch that appeared on the photosphere. The directly imaged brightness variations of Betelgeuse evolved on a timescale of weeks. Our findings suggest that a component of mass loss from red supergiants is inhomogeneous, linked to a very contrasted and rapidly changing photosphere

We present high-pass filtered continuum images of the inner $3.5^\circ\times2.5^\circ$ of the Galactic center at 20 cm with $6.4''$ resolution. These mosaic images are taken with MeerKAT and reveal a large number of narrow filaments, roughly an order of magnitude increase in their numbers compared to past measurements. For the first time, we carry out population studies of the spectral index and magnetic field of the entire region. The mean spectral indices of the filaments are steeper than supernova remnants (SNRs) (-0.62) with a value of $\alpha\sim-0.83$. The variation in $\alpha$ is much larger than for the SNRs, suggesting that these characteristics have a different origin. A large-scale cosmic-ray driven wind has recently been proposed to explain the origin of filaments and the large-scale 430 pc bipolar radio and X-ray structure. This favors the possibility that the large-scale bipolar radio/X-ray structure is produced by past activity of Sgr A* rather than coordinated burst of supernovae. A trend of steeper indices is also noted with increasing distance from the Galactic plane. This could be explained either by synchrotron cooling or weak shocks accelerating cosmic-ray particles in the context of the cosmic-ray driven wind. The mean magnetic field strengths along the filaments ranges from $\sim100$ to 400 $\mu$G depending on the assumed ratio of cosmic-ray protons to electrons. Given that there is a high cosmic ray pressure in the Galactic center, the large equipartition magnetic field implies that the magnetic filed is weak in most of the interstellar volume of the Galactic center.

Magnetic reconnection is a ubiquitous astrophysical process that rapidly converts magnetic energy into some combination of plasma flow energy, thermal energy, and non-thermal energetic particles, including energetic electrons. Various reconnection acceleration mechanisms in different low-$\beta$ (plasma-to-magnetic pressure ratio) and collisionless environments have been proposed theoretically and studied numerically, including first- and second-order Fermi acceleration, betatron acceleration, parallel electric field acceleration along magnetic fields, and direct acceleration by the reconnection electric field. However, none of them have been heretofore confirmed experimentally, as the direct observation of non-thermal particle acceleration in laboratory experiments has been difficult due to short Debye lengths for \textit{in-situ} measurements and short mean free paths for \textit{ex-situ} measurements. Here we report the direct measurement of accelerated non-thermal electrons from low-$\beta$ magnetically driven reconnection in experiments using a laser-powered capacitor coil platform. We use kiloJoule lasers to drive parallel currents to reconnect MegaGauss-level magnetic fields in a quasi-axisymmetric geometry. The angular dependence of the measured electron energy spectrum and the resulting accelerated energies, supported by particle-in-cell simulations, indicate that the mechanism of direct electric field acceleration by the out-of-plane reconnection electric field is at work. Scaled energies using this mechanism show direct relevance to astrophysical observations. Our results therefore validate one of the proposed acceleration mechanisms by reconnection, and establish a new approach to study reconnection particle acceleration with laboratory experiments in relevant regimes.

14N(p,gamma)15O is one of the key reactions of nuclear astrophysics playing a role in various stellar processes and influencing energy generation of stars, stellar evolution and nucleosynthesis. For a reliable reaction rate calculation the low energy cross section of 14N(p,gamma)15O must be known with high accuracy. Owing to the unmeasurable low cross sections, theoretical calculations are unavoidable. High precision experimental cross section data are needed in a wide energy range in order to provide the necessary basis for low energy extrapolations. In the present work the total 14N(p,gamma)15O cross section was measured with a method complementary to the available data sets. The cross section was measured with activation, based on the detection of the annihilation radiation following the beta+ decay of the reaction product 15O. This method, which provides directly the astrophysically important total cross section, was never used for the 14N(p,gamma)15O cross section measurement in the studied energy range. The non-resonant cross section was measured between 550 keV and 1400 keV center-of-mass energies with total uncertainty of about 10%. The results were compared with literature data using an R-matrix analysis. It is found that the cross sections measured in this work are in acceptable agreement with the two recent measurements only if the weak transitions - not measured in those works - are included. The present data set, being largely independent from the other available data, can be used to constrain the extrapolated cross sections to astrophysical energies and helps to make the astrophysical model calculations more reliable.

The presence of light sterile neutrinos is one of the unanswered questions of particle physics. The cosmological counterpart is represented by dark radiation, i.e.\ any form of radiation present in the early Universe besides photons and standard (active) neutrinos. This short review provides a comprehensive overview of the two problems and of their connection. We review the status of neutrino oscillation anomalies, commenting on the most recent oscillation data and their mutual tensions, and we discuss the constraints from other terrestrial probes. We show the shortcomings of translating light sterile neutrinos in cosmology as additional thermalised relativistic species, produced by neutrino oscillations, and we detail alternative solutions, specifically focusing on neutrino non standard interactions, and on their link to the Hubble constant problem. The impact of a new force leading to dark radiation -- dark matter interactions is also discussed in the realm of new physics in the dark sector.

We investigate the maximum force for black holes with charge and rotation as well as for the corresponding Buchdahl stars with limiting compactness bound. In terms of gravitational potential black hole and Buchdahl star are respectively characterized by $\Phi(r)= 1/2, 4/9$. The compactness ratio in the two cases bear the relation, $(M/R)_{Buch} = (8/9) (M/R)_{BH}$ with $\alpha^2 =Q^2/M^2 \to (8/9)\alpha^2, \beta^2 =a^2/M^2 \to (8/9)^2\beta^2$. As a consequence the maximum force in the two cases are similarly related. Further the maximum force for uncharged static black hole or Buchdahl star is entirely given in terms of the fundamental constant velocity of light and the gravitational constant only in four dimension in general relativity (GR). In general this feature uniquely picks out the pure Lovelock gravity (having only one $N$th order term in action which includes GR in the linear order $N=1$) and the dimensional spectrum, $D=3N+1$, where $N$ is degree of the Lovelock polynomial action.

The pseudospectral method is a highly accurate numerical scheme suitable for turbulence simulations. We have developed an open-source pseudospectral code, \textsc{\textsf{Calliope}}, which adopts the P3DFFT library \citep{Pekurovsky2012} to perform a fast Fourier transform with the two-dimensional (pencil) decomposition of numerical grids. \textsc{\textsf{Calliope}} can solve incompressible magnetohydrodynamics (MHD), isothermal compressible MHD, and rotational reduced MHD with parallel computation using very large numbers of cores ($> 10^5$ cores for $2048^3$ grids). The code can also solve for local magnetorotational turbulence in a shearing frame using the remapping method \citep{Rogallo1981,Umurhan2004}. \textsc{\textsf{Calliope}} is currently the only pseudospectral code that can compute magnetorotational turbulence using pencil-domain decomposition. This paper presents the numerical scheme of \textsc{\textsf{Calliope}} and the results of linear and nonlinear numerical tests, including compressible local magnetorotational turbulence with the largest grid number reported to date.