Papers for Thursday, Sep 09 2021

Planetary magnetic fields are generated by motions of electrically conducting fluids in their interiors. The dynamo problem has thus received much attention in spherical geometries, even though planetary bodies are non-spherical. To go beyond the spherical assumption, we develop an algorithm that exploits a fully spectral description of the magnetic field in triaxial ellipsoids to solve the induction equation with local boundary conditions (i.e. pseudo-vacuum or perfectly conducting boundaries). We use the method to compute the free-decay magnetic modes and to solve the kinematic dynamo problem for prescribed flows. The new method is thoroughly compared with analytical solutions and standard finite-element computations, which are also used to model an insulating exterior. We obtain dynamo magnetic fields at low magnetic Reynolds numbers in ellipsoids, which could be used as simple benchmarks for future dynamo studies in such geometries. We finally discuss how the magnetic boundary conditions can modify the dynamo onset, showing that a perfectly conducting boundary can strongly weaken dynamo action, whereas pseudo-vacuum and insulating boundaries often give similar results.

We apply the multiple point criticality principle to inflationary model building and study Coleman-Weinberg inflation when the scalar potential is quadratic in the logarithmic correction. We analyze also the impact of a non-minimal coupling to gravity under two possible gravity formulation: metric or Palatini. We compare the eventual compatibility of the results with the final data release of the Planck mission.

Accurate mass-loss rates are essential for meaningful stellar evolutionary models. For massive single stars with initial masses between 8 - 30\msun the implementation of cool supergiant mass loss in stellar models strongly affects the resulting evolution, and the most commonly used prescription for these cool-star phases is that of de Jager. Recently, we published a new \mdot\ prescription calibrated to RSGs with initial masses between 10 - 25\msun, which unlike previous prescriptions does not over estimate \mdot\ for the most massive stars. Here, we carry out a comparative study to the MESA-MIST models, in which we test the effect of altering mass-loss by recomputing the evolution of stars with masses 12-27\msun\ with the new \mdot-prescription implemented. We show that while the evolutionary tracks in the HR diagram of the stars do not change appreciably, the mass of the H-rich envelope at core-collapse is drastically increased compared to models using the de Jager prescription. This increased envelope mass would have a strong impact on the Type II-P SN lightcurve, and would not allow stars under 30\msun\ to evolve back to the blue and explode as H-poor SN. We also predict that the amount of H-envelope around single stars at explosion should be correlated with initial mass, and we discuss the prospects of using this as a method of determining progenitor masses from supernova light curves.

Pulsar timing provides a sensitive probe of small-scale structure. Gravitational perturbations arising from an inhomogeneous environment could manifest as detectable perturbations in the pulsation phase. Consequently, pulsar timing arrays have been proposed as a probe of dark matter substructure on mass scales as small as $10^{-11}$ $M_\odot$. Since the small-scale mass distribution is connected to early-Universe physics, pulsar timing can therefore constrain the thermal history prior to Big Bang Nucleosynthesis (BBN), a period that remains largely unprobed. We explore here the prospects for pulsar timing arrays to detect the dark substructure imprinted by a period of early matter domination (EMD) prior to BBN. EMD amplifies density variations, leading to a population of highly dense sub-Earth-mass dark matter microhalos. We use recently developed semianalytic models to characterize the distribution of EMD-induced microhalos, and we evaluate the extent to which the pulsar timing distortions caused by these microhalos can be detected. Broadly, we find that sub-0.1-$\mu$s timing noise residuals are necessary to probe EMD. However, with 10-ns residuals, a pulsar timing array with just 70 pulsars could detect the evidence of an EMD epoch with 20 years of observation time if the reheat temperature is of order 10 MeV. With 40 years of observation time, pulsar timing arrays could probe EMD reheat temperatures as high as 150 MeV.

Variability of a galaxy's core radio source can be a significant consequence of AGN accretion. However, this variability has not been well studied, particularly at high radio frequencies. As such, we report on a campaign monitoring the high radio frequency variability of 20 nearby, cool-core brightest cluster galaxies. Our highest cadence observations are at 15 GHz and are from the Owens Valley Radio Observatory (OVRO). They have a median time interval of 7 days and mostly span between 8 and 13 years. We apply a range of variability detection techniques to the lightcurves of the sources to analyse changes in their flux density on week to decade long timescales. Over the full period in which each source was observed, $\chi^{2}$ tests suggest that 13/20 are inconsistent with the flat lightcurve of a non-varying source. Variability amplitude tests suggest that 12/20 sources are variable on 300 day timescales, while 19/20 are variable on 3000 day timescales. At least half of the sources also show 20 per cent peak to trough variability on 3~year timescales, while at least a third vary by 60 per cent on 6~year timescales. Significant variability is therefore a common feature of these sources. We also show how the variability relates to spectral properties at frequencies of up to 353 GHz using data from the Korean VLBI network (KVN), the NIKA2 instrument of the IRAM 30m telescope, and the SCUBA-2 instrument of the James Clerk Maxwell Telescope.

We derive solutions to transit light curves of exoplanets orbiting rapidly-rotating stars. These stars exhibit significant oblateness and gravity darkening, a phenomenon where the poles of the star have a higher temperature and luminosity than the equator. Light curves for exoplanets transiting these stars can exhibit deviations from those of slowly-rotating stars, even displaying significantly asymmetric transits depending on the system's spin-orbit angle. As such, these phenomena can be used as a protractor to measure the spin-orbit alignment of the system. In this paper, we introduce a novel semi-analytic method for generating model light curves for gravity-darkened and oblate stars with transiting exoplanets. We implement the model within the code package starry and demonstrate several orders of magnitude improvement in speed and precision over existing methods. We test the model on a TESS light curve of WASP-33, whose host star displays rapid rotation ($v \sin i_* = 86.4$ km/s). We subtract the host's $\delta$-Scuti pulsations from the light curve, finding an asymmetric transit characteristic of gravity darkening. We find the projected spin orbit angle is consistent with Doppler tomography and constrain the true spin-orbit angle of the system as $\varphi=108.3^{+19.0}_{-15.4}$~$^{\circ}$. We demonstrate the method's uses in constraining spin-orbit inclinations of such systems photometrically with posterior inference. Lastly, we note the use of such a method for inferring the dynamical history of thousands of such systems discovered by TESS.

We explore the Baryonic Tully-Fisher Relation in the Local Group. Rotationally supported Local Group galaxies adhere precisely to the relation defined by more distant galaxies. For pressure supported dwarf galaxies, we determine the scaling factor $\beta_c$ that relates their observed velocity dispersion to the equivalent circular velocity of rotationally supported galaxies of the same mass such that $V_o = \beta_c \sigma_*$. For a typical mass-to-light ratio $\Upsilon_* = 2\;\mathrm{M}_{\odot}/\mathrm{L}_{\odot}$ in the $V$-band, we find that $\beta_c = 2$. More generally, $\log \beta_c = 0.25 \log \Upsilon_* +0.226$. This provides a common kinematic scale relating pressure and rotationally supported dwarf galaxies.

We present a study of the relation between X-rays and ultraviolet emission in quasars for a sample of broad-line, radio-quiet objects obtained from the cross-match of the Sloan Digital Sky Survey DR14 with the latest Chandra Source Catalog 2.0 (2,332 quasars) and the Chandra COSMOS Legacy survey (273 quasars). The non-linear relation between the ultraviolet (at 2500 A, $L_{O}$) and the X-ray (at 2 keV, $L_{X}$) emission in quasars has been proved to be characterised by a smaller intrinsic dispersion than the observed one, as long as a homogeneous selection, aimed at preventing the inclusion of contaminants in the sample, is fulfilled. By leveraging on the low background of Chandra, we performed a complete spectral analysis of all the data available for the SDSS-CSC2.0 quasar sample (i.e. 3,430 X-ray observations), with the main goal of reducing the uncertainties on the source properties (e.g. flux, spectral slope). We analysed whether any evolution of the $L_{X}-L_{O}$ relation exists by dividing the sample in narrow redshift intervals across the redshift range spanned by our sample, $z \simeq 0.5-4$. We find that the slope of the relation does not evolve with redshift and it is consistent with the literature value of $0.6$ over the explored redshift range, implying that the mechanism underlying the coupling of the accretion disc and hot corona is the same at the different cosmic epochs. We also find that the dispersion decreases when examining the highest redshifts, where only pointed observations are available. These results further confirm that quasars are `standardisable candles', that is we can reliably measure cosmological distances at high redshifts where very few cosmological probes are available.

Dual-Active Galactic Nuclei (AGN) are a natural consequence of the hierarchical structure formation scenario, and can provide an important test of various models for black hole growth. However, due to their rarity and difficulty to find at high redshift, very few confirmed dual-AGN are known at the epoch where galaxy formation peaks. Here we report the discovery of a gravitationally lensed dual-AGN system at redshift 2.37 comprising two optical/IR quasars separated by 6.5+/-0.6 kpc, and a third compact (R_eff = 0.45+/-0.02 kpc) red galaxy that is offset from one of the quasars by 1.7+/-0.1 kpc. From Very Large Array imaging at 3 GHz, we detect 600 and 340 pc-scale radio emission that is associated with both quasars. The 1.4 GHz luminosity densities of the radio sources are about 10^24.35 W / Hz, which is consistent with weak jets. However, the low brightness temperature of the emission is also consistent with star-formation at the level of 850 to 1150 M_sun / yr. Although this supports the scenario where the AGN and/or star-formation is being triggered through an ongoing triple-merger, a post-merger scenario where two black holes are recoiling is also possible, given that neither has a detected host galaxy.

We explore constraints on dark energy and modified gravity with forecast 21cm intensity mapping measurements using the Effective Field Theory approach. We construct a realistic mock data set forecasting a low redshift 21cm signal power spectrum $P_{21}(z,k)$ measurement from the MeerKAT radio-telescope. We compute constraints on cosmological and model parameters through Monte Carlo Markov chain techniques, testing both the constraining power of $P_{21}(k)$ alone and its effect when combined with the latest Planck 2018 CMB data. We complement our analysis by testing the effects of tomography from an ideal mock data set of observations in multiple redshift bins. We conduct our analysis numerically with the codes EFTCAMB/EFTCosmoMC, which we extend by implementing a likelihood module fully integrated with original codes. We find that adding $P_{21}(k)$ to CMB data provides significantly tighter constraints on $\Omega_ch^2$ and $H_0$, with a reduction of the error with respect to Planck results at the level of more than $60\%$. For the parameters describing beyond $\Lambda$CDM theories, we observe a reduction in the error with respect to the Planck constraints at the level of $\lesssim 10\%$. The improvement increases up to $\sim 35\%$ when we constrain the parameters using ideal, tomographic mock observations. We conclude that the power spectrum of the 21cm signal is sensitive to variations of the parameters describing the examined beyond $\Lambda$CDM models and, thus, $P_{21}(k)$ observations could help to constrain dark energy. The constraining power on such theories is improved significantly by tomography.

We use the Cluster-EAGLE (C-EAGLE) hydrodynamical simulations to investigate the effects of self-interacting dark matter (SIDM) on galaxies as they fall into clusters. We find that SIDM galaxies follow similar orbits to their Cold Dark Matter (CDM) counterparts, but end up with ${\sim}$25 per cent less mass by the present day. One in three SIDM galaxies are entirely disrupted, compared to one in five CDM galaxies. However, the excess stripping will be harder to observe than suggested by previous DM-only simulations because the most stripped galaxies form cores and also lose stars: the most discriminating objects become unobservable. The best test will be to measure the stellar-to-halo mass relation (SHMR) for galaxies with stellar mass $10^{10-11}\,\mathrm{M}_{\odot}$. This is 8 times higher in a cluster than in the field for a CDM universe, but 13 times higher for an SIDM universe. Given intrinsic scatter in the SHMR, these models could be distinguished with noise-free galaxy-galaxy strong lensing of ${\sim}32$ cluster galaxies.

Axion-Like Particles (ALPs) are predicted by several Beyond the Standard Model theories, in particular, string theory. In the presence of an external magnetic field perpendicular to the direction of propagation, ALPs can couple to photons. Therefore, if an X-ray source is viewed through a magnetised plasma, such as a luminous quasar in a galaxy cluster, we may expect spectral distortions that are well described by photon-ALP oscillations. We present a $571 \ \mathrm{ks}$ combined High and Low Energy Transmission Grating (HETG/LETG) Chandra observation of the powerful radio-quiet quasar H1821+643, hosted by a cool-core cluster at redshift $0.3$. The spectrum is well described by a double power-law continuum and broad$+$narrow iron line emission typical of type-1 Active Galactic Nuclei (AGN), with remaining spectral features $< 2.5\%$. Using a cell-based approach to describe the turbulent cluster magnetic field, we compare our spectrum with photon-ALP mixing curves for 500 field realisations assuming that the thermal-to-magnetic pressure ratio remains constant up to the virial radius. At $99.7\%$ credibility, we exclude all couplings $g_\mathrm{a\gamma} > 6.3 \times 10^{-13} \ {\mathrm{GeV}}^{-1}$ for most ALP masses $< 10^{-12} \ \mathrm{eV}$. Our results are moderately more sensitive to constraining ALPs than the best previous result from Chandra observations of the Perseus cluster, albeit with a less constrained field model. We provide the best constraints on light ALPs, exceeding the projected sensitivity of next generation axion helioscopes. We reflect on the promising future of ALP studies with bright AGN embedded in rich clusters, especially with the upcoming Athena mission.

The next two decades are expected to open the door to the first coincident detections of electromagnetic (EM) and gravitational wave (GW) signatures associated with massive black hole (MBH) binaries heading for coalescence. These detections will launch a new era of multimessenger astrophysics by expanding this growing field to the low-frequency GW regime and will provide unprecedented understanding of the evolution of MBHs and galaxies. They will also constitute fundamentally new probes of cosmology and would enable unique tests of gravity. The aim of this Living Review is to provide an introduction to this research topic by presenting a summary of key findings, physical processes and ideas pertaining to EM counterparts to MBH mergers as they are known at the time of this writing. We review current observational evidence for close MBH binaries, discuss relevant physical processes and timescales, and summarize the possible EM counterparts to GWs in the precursor, coalescence, and afterglow stages of a MBH merger. We also describe open questions and discuss future prospects in this dynamic and quick paced research area.

We propose a new analysis of small scale CMB data by introducing the cosmological dependency of the foreground signals, focusing first on the thermal Sunyaev-Zel'dovich (tSZ) power spectrum, derived from the halo model. We analyse the latest observations by the South Pole Telescope (SPT) of the high-$\ell$ power (cross) spectra at 90, 150 and 220 GHz, as the sum of CMB and tSZ signals, both depending on cosmological parameters, and remaining contaminants. In order to perform faster analyses, we propose a new tSZ modelling based on machine learning algorithms (namely Random Forest). We show that the additional information contained in the tSZ power spectrum tightens constraints on cosmological and tSZ scaling relation parameters. We combine for the first time the Planck tSZ data with SPT high-$\ell$ to derive even stronger constraints. Finally, we show how the amplitude of the remaining kSZ power spectrum varies depending on the assumptions made on both tSZ and cosmological parameters.

Collisional self-interactions occurring in protostellar jets give rise to strong shocks, the structure of which can be affected by radiative cooling within the flow. To study such colliding flows, we use the AstroBEAR AMR code to conduct hydrodynamic simulations in both one and three dimensions with a power law cooling function. The characteristic length and time scales for cooling are temperature dependent and thus may vary as shocked gas cools. When the cooling length decreases sufficiently rapidly the system becomes unstable to the radiative shock instability, which produces oscillations in the position of the shock front; these oscillations can be seen in both the one and three dimensional cases. Our simulations show no evidence of the density clumping characteristic of a thermal instability, even when the cooling function meets the expected criteria. In the three-dimensional case, the nonlinear thin shell instability (NTSI) is found to dominate when the cooling length is sufficiently small. When the flows are subjected to the radiative shock instability, oscillations in the size of the cooling region allow NTSI to occur at larger cooling lengths, though larger cooling lengths delay the onset of NTSI by increasing the oscillation period.

Over half a century from the discovery of gamma-ray bursts (GRBs), the dominant radiation mechanism responsible for their bright and highly variable prompt emission remains poorly understood. Spectral information alone has proven insufficient for understanding the composition and main energy dissipation mechanism in GRB jets. High-sensitivity polarimetric observations from upcoming instruments in this decade may help answer such key questions in GRB physics. This article reviews the current status of prompt GRB polarization measurements and provides comprehensive predictions from theoretical models. A concise overview of the fundamental questions in prompt GRB physics is provided. Important developments in gamma-ray polarimetry including a critical overview of different past instruments are presented. Theoretical predictions for different radiation mechanisms and jet structures are confronted with time-integrated and time-resolved measurements. The current status and capabilities of upcoming instruments regarding the prompt emission are presented. The very complimentary information that can be obtained from polarimetry of X-ray flares as well as reverse-shock and early to late forward-shock (afterglow) emission is highlighted. Finally, promising directions for overcoming the inherent difficulties in obtaining statistically significant prompt-GRB polarization measurements are discussed, along with prospects for improvements in the theoretical modeling, which may lead to significant advances in the field.

The structure, elemental abundances, physical conditions of the LMC supernova remnant (SNR) 0540-69.3 and its surroundings were investigated using [O III] imaging and spectroscopy. Several new spectral lines are identified, both in central filaments and in interstellar clouds shocked by the supernova blast wave. The central lines are redshifted by $440\pm80$ km s$^{-1}$ with respect to the LMC, and the [O III] emission displays a symmetry axis of ring-like structures which could indicate that the pulsar shares the same general redshift as the central supernova ejecta. [O II], [S II], [Ar III] and H$\beta$ have more compact structures than [O III], and possibly [Ne III]. The average [O III] temperature is $23\,500 \pm 1\,800$ K, and the electron density from [S II] is typically $10^3$ cm$^{-3}$. By mass, the relative elemental abundances of the central shocked ejecta are ${\rm O:Ne:S:Ar} \approx 1:0.07:0.10:0.02$, consistent with explosion models of $13-20$ solar mass progenitors, and similar to that of SN 1987A, as is also the mixing of hydrogen and helium into the center. [O III] is also seen in freely coasting ejecta outside the pulsar-wind nebula out to well above $2\,000$ km s$^{-1}$. From this a pulsar age of $\approx 1\,200$ years is estimated. Four filaments of shocked interstellar medium with a wide range in degree of ionization of iron are identified. One was observed in X-rays, and another has a redshift of $85\pm30$ km~s$^{-1}$ relative to LMC. From this the electron density of the [O III]-emitting gas is estimated to be $10^3$ cm$^{-3}$. The line of the most highly ionized ion, [Fe XIV] $\lambda$5303, likely comes from an evaporation zone in connection with the radiatively cooled gas emitting, e.g., [O III].

Trans-Neptunian Objects (TNOs) provide a window into the history of the Solar System, but they can be challenging to observe due to their distance from the Sun and relatively low brightness. Here we report the detection of 75 moving objects that we could not link to any other known objects, the faintest of which has a VR magnitude of $25.02 \pm 0.93$ using the KBMOD platform. We recover an additional 24 sources with previously-known orbits. We place constraints on the barycentric distance, inclination, and longitude of ascending node of these objects. The unidentified objects have a median barycentric distance of 41.28 au, placing them in the outer Solar System. The observed inclination and magnitude distribution of all detected objects is consistent with previously published KBO distributions. We describe extensions to KBMOD, including a robust percentile-based lightcurve filter, an in-line graphics processing unit (GPU) filter, new coadded stamp generation, and a convolutional neural network (CNN) stamp filter, which allow KBMOD to take advantage of difference images. These enchancements mark a significant improvement in the readiness of KBMOD for deployment on future big data surveys such as LSST.

Using the CLASH-VLT survey, we assembled an unprecedented sample of 1234 spectroscopically confirmed members in Abell~S1063, finding a dynamically complex structure at z_cl=0.3457 with a velocity dispersion \sigma_v=1380 -32 +26 km s^-1. We investigate cluster environmental and dynamical effects by analysing the projected phase-space diagram and the orbits as a function of galaxy spectral properties. We classify cluster galaxies according to the presence and strength of the [OII] emission line, the strength of the H$\delta$ absorption line, and colours. We investigate the relationship between the spectral classes of galaxies and their position in the projected phase-space diagram. We analyse separately red and blue galaxy orbits. By correlating the observed positions and velocities with the projected phase-space constructed from simulations, we constrain the accretion redshift of galaxies with different spectral types. Passive galaxies are mainly located in the virialised region, while emission-line galaxies are outside r_200, and are accreted later into the cluster. Emission-lines and post-starbursts show an asymmetric distribution in projected phase-space within r_200, with the first being prominent at Delta_v/sigma <~-1.5$, and the second at Delta_v/ sigma >~ 1.5, suggesting that backsplash galaxies lie at large positive velocities. We find that low-mass passive galaxies are accreted in the cluster before the high-mass ones. This suggests that we observe as passives only the low-mass galaxies accreted early in the cluster as blue galaxies, that had the time to quench their star formation. We also find that red galaxies move on more radial orbits than blue galaxies. This can be explained if infalling galaxies can remain blue moving on tangential orbits.

The effect of radiation pressure compression (RPC) on ionized gas in Active Galactic Nuclei (AGN) likely sets the photoionized gas density structure. The photoionized gas free-free absorption and emission are therefore uniquely set by the incident ionizing flux. We use the photoionization code Cloudy RPC model results to derive the expected relations between the free-free emission and absorption properties and the distance from the AGN centre, for a given AGN luminosity. The free-free absorption frequency of RPC gas is predicted to increases from $\sim$100 MHz on the kpc scale, to $\sim$100 GHz on the sub-pc scale, consistent with observations of spatially resolved free-free absorption. The free-free emission at 5 GHz is predicted to yield a radio loudness of $R\sim 0.03$, below the typical observed values of $R\sim 0.1-1$ in radio-quiet AGN. However, the flat free-free radio continuum may become dominant above 100 GHz. The suggested detection of optically thin free-free emission in NGC 1068, on the sub pc torus scale, is excluded as the brightness temperature is too high for optically thin free-free emission. However, excess emission observed with ALMA above 150 GHz in NGC 1068, is consistent with the predicted free-free emission from gas just outside the broad line region, a region which overlaps the hot dust disc resolved with GRAVITY. Extended $\sim$100 pc scale free-free emission is also likely present in NGC 1068. Future sub mm observation of radio quiet AGN with ALMA may allow to image the free-free emission of warm photoionized gas in AGN down to the 30 mas scale, including highly absorbed AGN.

The Santa Cruz Extreme AO Lab (SEAL) is a new visible-wavelength testbed designed to advance the state of the art in wavefront control for high contrast imaging on large, segmented, ground-based telescopes. SEAL provides multiple options for simulating atmospheric turbulence, including rotating phase plates and a custom Meadowlark spatial light modulator that delivers phase offsets of up to 6pi at 635nm. A 37-segment IrisAO deformable mirror (DM) simulates the W. M. Keck Observatory segmented primary mirror. The adaptive optics system consists of a woofer/tweeter deformable mirror system (a 97-actuator ALPAO DM and 1024-actuator Boston Micromachines MEMs DM, respectively), and four wavefront sensor arms: 1) a high-speed Shack-Hartmann WFS, 2) a reflective pyramid WFS, designed as a prototype for the ShaneAO system at Lick Observatory, 3) a vector-Zernike WFS, and 4) a Fast Atmospheric Self Coherent Camera Technique (FAST) demonstration arm, consisting of a custom focal plane mask and high-speed sCMOS detector. Finally, science arms preliminarily include a classical Lyot-style coronagraph as well as FAST (which doubles as a WFS and science camera). SEAL's real time control system is based on the Compute and Control for Adaptive optics (CACAO) package, and is designed to support the efficient transfer of software between SEAL and the Keck II AO system. In this paper, we present an overview of the design and first light performance of SEAL.

The Planetary Systems Imager (PSI) is a proposed instrument for the Thirty Meter Telescope (TMT) that provides an extreme adaptive optics (AO) correction to a multi-wavelength instrument suite optimized for high contrast science. PSI's broad range of capabilities, spanning imaging, polarimetry, integral field spectroscopy, and high resolution spectroscopy from 0.6-5 microns, with a potential channel at 10 microns, will enable breakthrough science in the areas of exoplanet formation and evolution. Here, we present a preliminary optical design and performance analysis toolset for the 2-5 microns component of the PSI AO system, which must deliver the wavefront quality necessary to support infrared high contrast science cases. PSI-AO is a two-stage system, with an initial deformable mirror and infrared wavefront sensor providing a common wavefront correction to all PSI science instruments followed by a dichroic that separates "PSI-Red" (2-5 microns) from "PSI-Blue" (0.5-1.8 microns). To meet the demands of visible-wavelength high contrast science, the PSI-Blue arm will include a second deformable mirror and a visible-wavelength wavefront sensor. In addition to an initial optical design of the PSI-Red AO system, we present a preliminary set of tools for an end-to-end AO simulation that in future work will be used to demonstrate the planet-to-star contrast ratios achievable with PSI-Red.

We present 226 large ultra-diffuse galaxy (UDG) candidates ($r_e > 5.3$\arcsec, $\mu_{0,g} > 24$ mag arcsec$^{-2}$) in the SDSS Stripe 82 region recovered using our improved procedure developed in anticipation of processing the entire Legacy Surveys footprint. The advancements include less constrained structural parameter fitting, expanded wavelet filtering criteria, consideration of Galactic dust, estimates of parameter uncertainties and completeness based on simulated sources, and refinements of our automated candidate classification. We have a sensitivity $\sim$1 mag fainter in $\mu_{0,g}$ than the largest published catalog of this region. Using our completeness-corrected sample, we find that (1) there is no significant decline in the number of UDG candidates as a function of $\mu_{0,g}$ to the limit of our survey ($\sim$ 26.5 mag arcsec$^{-2}$); (2) bluer candidates have smaller S\'ersic $n$; (3) most blue ($g-r < 0.45$ mag) candidates have $\mu_{0,g} \lesssim 25$ mag arcsec$^{-2}$ and will fade to populate the UDG red sequence we observe to $\sim 26.5$ mag arcsec$^{-2}$; (4) any red UDGs that exist significantly below our $\mu_{0,g}$ sensitivity limit are not descended from blue UDGs in our sample; and (5) candidates with lower $\mu_{0,g}$ tend to smaller $n$. We anticipate that the final SMUDGes sample will contain $\sim$ 30$\times$ as many candidates.

We present the first results of K2-OjOS, a collaborative project between professional and amateur astronomers primarily aimed to detect, characterize and validate new extrasolar planets. For this work, 10 amateur astronomers looked for planetary signals by visually inspecting the 20,427 light curves of K2 campaign 18 (C18). They found 42 planet candidates, of which 18 are new detections and 24 had been detected in the overlapping C5 by previous works. We used archival photometric and spectroscopic observations, as well as new high-spatial resolution images in order to carry out a complete analysis of the candidates found, including a homogeneous characterization of the host stars, transit modeling, search for transit timing variations and statistical validation. As a result, we report four new planets (K2-XXX b, K2-XXX b, K2-XXX b, and K2-XXX b) and 14 planet candidates. Besides, we refine the transit ephemeris of the previously published planets and candidates by modeling C5, C16 (when available) and C18 photometric data jointly, largely improving the period and mid-transit time precision. Regarding individual systems, we highlight the new planet K2-XXX b and candidate EPIC211537087.02 being near a 2:1 period commensurability, the detection of significant TTVs in the bright star K2-184 (V = 10.35), the location of K2-103 b inside the habitable zone according to optimistic models, the detection of a new single transit in the known system K2-274, and the disposition reassignment of K2-120 b, which we consider as a planet candidate as the origin of the signal can not be ascertained.

We introduce our new code, SMERCURY-T, which is based on existing codes SMERCURY (Lissauer et al. 2012) and Mercury-T (Bolmont et al. 2015). The result is a mixed-variable symplectic N-body integrator that can compute the orbital and spin evolution of a planet within a multi-planet system under the influence of tidal spin torques from its star. We validate our implementation by comparing our experimental results to that of a secular model. As we demonstrate in a series of experiments, SMERCURY-T allows for the study of secular spin-orbit resonance crossings and captures for planets within complex multi-planet systems. These processes can drive a planet's spin state to evolve along vastly different pathways on its road toward tidal equilibrium, as tidal spin torques dampen the planet's spin rate and evolve its obliquity. Additionally, we show the results of a scenario that exemplifies the crossing of a chaotic region that exists as the overlap of two spin-orbit resonances. The test planet experiences violent and chaotic swings in its obliquity until its eventual escape from resonance as it tidally evolves. All of these processes are and have been important over the obliquity evolution of many bodies within the Solar System and beyond, and have implications for planetary climate and habitability. SMERCURY-T is a powerful and versatile tool that allows for further study of these phenomena.

There is a wide consensus that the ubiquitous presence of magnetic reconnection events and the associated impulsive heating (nanoflares) is a strong candidate for solving the solar coronal heating problem. Whether nanoflares accelerate particles to high energies like full-sized flares is unknown. We investigate this question by studying the type III radio bursts that the nanoflares may produce on closed loops. The characteristic frequency-drifts that type III bursts exhibit can be detected using a novel application of the time-lag technique developed by Viall & Klimchuk (2012) even when there are multiple overlapping bursts. We present a simple numerical model that simulates the expected radio emission from nanoflares in an active region (AR), which we use to test and calibrate the technique. We find that in the case of closed loops the frequency spectrum of type III bursts is expected to be extremely steep such that significant emission is produced at a given frequency only for a rather narrow range of loop lengths. We also find that the signature of bursts in the time-lag signal diminishes as: (1)the variety of participating loops within that range increases; (2)the occurrence rate of bursts increases; (3) the duration of bursts increases; and (4) the brightness of the bursts decreases relative to noise. In addition, our model suggests a possible origin of type I bursts as a natural consequence of type III emission in a closed-loop geometry.

We present and analyze a near infrared(NIR) spectrum of the under-luminous Type Ia supernova SN~2020qxp/ASASSN-20jq obtained with NIRES at the Keck Observatory 191 days after B-band maximum. The spectrum is dominated by a number of broad emission features including the [FeII] at 1.644mu which is highly asymmetric with a tilted top and a peak red-shifted by ~2,000km/s. In comparison with 2-D non-LTE synthetic spectra computed from 3-D simulations of off-center delayed-detonation Chandrasekhar-mass white-dwarf(WD) models, we find good agreement between the observed lines and the synthetic profiles, and are able to unravel the structure of the progenitor's envelope. We find that the size and tilt of the [Fe II] 1.644mu-profile (in velocity space) is an effective way to determine the location of an off-center delayed-detonation transition (DDT) and the viewing angle, and it requires a WD with a high central density of ~4E9$g/cm^3$. We also tentatively identify a stable Ni feature around 1.9mu characterized by a `pot-belly' profile that is slightly offset with respect to the kinematic center. In the case of SN~2020qxp/ASASSN-20jq, we estimate that the location of the DDT is ~0.3M(WD) off-center, which gives rise to an asymmetric distribution of the underlying ejecta. We also demonstrate that low-luminosity and high-density WD SNIa progenitors exhibit a very strong overlap of Ca and 56Ni in physical space. This results in the formation of a prevalent [Ca II] 0.73mu emission feature, which is sensitive to asymmetry effects. Our findings are discussed within the context of alternative scenarios, including off-center C/O detonations in He-triggered sub-M(Ch)-WDs and the direct collision of two WDs. Snapshot programs with Gemini/Keck/VLT/ELT class instruments and our spectropolarimetry program are complementary to mid-IR spectra by JWST.

Silicon monosulfide is an important silicon bearing molecule detected in circumstellar envelopes and star forming regions. Its formation and destruction routes are not well understood, partially due to the lack of a detailed knowledge on the involved reactions and their rate coefficients. In this work we have calculated and modeled the potential energy surface (PES) of the HSiS system employing highly accurate multireference electronic structure methods. After obtaining an accurate analytic representation of the PES, which includes long-range energy terms in a realistic way via the DMBE method, we have calculated rate coefficients for the Si+SH$\rightarrow$SiS+H reaction over the temperature range of 25-1000K. This reaction is predicted to be fast, with a rate coefficient of $\sim 1\times 10^{-10}\rm cm^3\, s^{-1}$ at 200K, which substantially increases for lower temperatures (the temperature dependence can be described by a modified Arrhenius equation with $\alpha=0.770\times 10^{-10}\rm cm^3\,s^{-1}$, $\beta=-0.756$ and $\gamma=9.873\, \rm K$). An astrochemical gas-grain model of a shock region similar to L1157-B1 shows that the inclusion of the Si+SH reaction increases the SiS gas-phase abundance relative to \ce{H2} from $5\times 10^{-10}$ to $1.4\times 10^{-8}$, which perfectly matches the observed abundance of $\sim 2\times 10^{-8}$.

Thanks to recent advances in integral field spectroscopy (IFS), modern surveys of nearby galaxies are capable of resolving metallicity maps of Hii regions down to scales of ~50pc. However, statistical analysis of these metallicity maps has seldom gone beyond fitting basic linear regressions and comparing parameters to global galaxy properties. In this paper (the first of a series), we introduce techniques from spatial statistics that are well suited for detailed analysis of both small- and large-scale metallicity variations within the interstellar media (ISMs) of local galaxies. As a first application, we compare the observed structure of small-scale metallicity fluctuations within 7 local galaxies observed by the PHANGS collaboration to predictions from a stochastic, physically motivated, analytical model developed by Krumholz & Ting. We show that while the theoretical model underestimates the amount of correlated scatter in the galactic metallicity distributions by 3-4 orders of magnitude, it provides good estimates of the physical scale of metallicity correlations. We conclude that the ISM of local spiral galaxies is far from homogeneous, with regions of size ~1 kpc showing significant departures from the mean metallicity at each galactocentric radius.

The advantage of having a high-fidelity instrument simulation tool developed in tandem with novel instrumentation is having the ability to investigate, in isolation and in combination, the wide parameter space set by the instrument design. SCALES, the third generation thermal-infrared diffraction limited imager and low/med-resolution integral field spectrograph being designed for Keck, is an instrument unique in design in order to optimize for its driving science case of direct detection and characterization of thermal emission from cold exoplanets. This warranted an end-to-end simulation tool that systematically produces realistic mock data from SCALES to probe the recovery of injected signals under changes in instrument design parameters. In this paper, we quantify optomechanical tolerance and detector electronic requirements set by the fiducial science cases using information content analysis, and test the consequences of updates to the design of the instrument on meeting these requirements.

It has been implicitly assumed that ices on grains in molecular clouds and proto planetary disks are formed by homogeneous layers regardless of their composition or crystallinity. To verify this assumption, we observed the H2O deposition onto refractory substrates and the crystallization of amorphous ices (H2O, CO2, and CO) using an ultra-high-vacuum transmission electron microscope. In the H2O-deposition experiments, we found that three-dimensional islands of crystalline ice (Ic) were formed at temperatures above 130 K. The crystallization experiments showed that uniform thin films of amorphous CO and H2O became three-dimensional islands of polyhedral crystals; amorphous CO2, on the other hand, became a thin film of nano crystalline CO2 covering the amorphous H2O. Our observations show that crystal morphologies strongly depend not only on the ice composition, but also on the substrate. Using experimental data concerning the crystallinity of deposited ices and the crystallization timescale of amorphous ices, we illustrated the criteria for ice crystallinity in space and outlined the macroscopic morphology of icy grains in molecular clouds as follows: amorphous H2O covered the refractory grain uniformly, CO2 nano-crystals were embedded in the amorphous H2O, and a polyhedral CO crystal was attached to the amorphous H2O. Furthermore, a change in the grain morphology in a proto-planetary disk is shown. These results have important implications for the chemical evolution of molecules, non-thermal desorption, collision of icy grains, and sintering.

We explore the cosmic evolution of the bar length, strength, and light deficit around the bar for 379 barred galaxies at 0.2 < z $\leq$ 0.835 using F814W images from the COSMOS survey. Our sample covers galaxies with stellar mass 10.0 $\leq$ log(M*/Msun) $\leq$ 11.4 and various Hubble types. The bar length is strongly related to the galaxy mass, the disk scale length (h), R50, and R90, where the last two are the radii containing 50 and 90% of total stellar mass, respectively. Bar length remains almost constant, suggesting little or no evolution in bar length over the last 7 Gyrs. The normalized bar lengths (Rbar/h, Rbar/R50, and Rbar/R90) do not show any clear cosmic evolution. Also, the bar strength (A2 and Qb) and the light deficit around the bar reveal little or no cosmic evolution. The constancy of the normalized bar lengths over cosmic time implies that the evolution of bars and of disks is strongly linked over all times. We discuss our results in the framework of predictions from numerical simulations. We conclude there is no strong disagreement between our results and up-to-date simulations.

High-redshift dusty star-forming galaxies with very high star formation rates (500 -- 3000 M$_\odot$ yr$^{-1}$) are key to understanding the formation of the most extreme galaxies in the early Universe. Characterising the gas reservoir of these systems can reveal the driving factor behind the high star formation. Using molecular gas tracers like high-J CO lines, neutral carbon lines and the dust continuum, we can estimate the gas density and radiation field intensity in their interstellar medium. In this paper, we present high resolution($\sim$0.4") observations of CO(7-6), [CI](2-1) and dust continuum of 3 lensed galaxies from the SPT-SMG sample at z$\sim$3 with the Atacama Large Millimeter/submillimeter Array. Our sources have high intrinsic star-formation rates (>850 M$_\odot$ yr$^{-1}$) and rather short depletion timescales (<100 Myr). Based on the $\rm L_{[\rm CI](2-1)}/ \rm L_{\rm CO(7-6)}$ and $\rm L_{[\rm CI](2-1)}/\rm L_{\rm IR}$ ratios, our sample galaxies exhibit higher radiation field intensity compared to other submillimetre galaxies but have similar gas densities. We perform visibility-based lens modelling on these objects to reconstruct the kinematics in the source plane. We find that the cold gas masses of the sources are compatible with simple dynamical mass estimates using ULIRG-like values of the CO-H$_2$ conversion factor $\alpha_{\rm CO}$ but not Milky Way-like values. We find diverse source kinematics in our sample: SPT0103-45 and SPT2147-50 are likely rotating disks while SPT2357-51 is a probable major merger. The analysis presented in the paper could be extended to a larger sample to determine better statistics of morphologies and interstellar medium properties of high-z dusty star-forming galaxies.

Quasi-periodic eruptions (QPEs), which are a new kind of X-ray bursts with the recurrence time of several hours, have been detected from the central supermassive black holes (SMBHs) of galactic nuclei, both active and quiescent. Recently, the two newly QPEs discovered by the eROSITA show asymmetric light curves with a fast rise and a slow decline. Current models cannot explain the observational characteristics of QPEs. Here we show that QPEs can be generated from the Roche lobe overflows at each pericentre passage of an evolved star orbiting an SMBH. The evolved stars with masses of $1-10~M_\odot$, which have lost Hydrogen envelopes in the post asymptotic giant branch phase, can fulfill the requirement to produce the properties of QPEs, including the fast rise and slow decay light curves, periods, energetics, and rates. Furthermore, the extreme mass ratio $\sim 10^5$ between the SMBH and the companion star will produce millihertz gravitational waves, called extreme mass-ratio inspirals (EMRIs). These QPEs would be detected as EMRI sources with electromagnetic counterparts for space-based GW detectors, such as Laser Interferometer Space Antenna (LISA) and Tianqin. They would provide a new way to measure the Hubble constant and further test the so-called Hubble constant tension.

X-ray reverberation is a powerful technique which maps out the structure of the inner regions of accretion disks around black holes using the echoes of the coronal emission reflected by the disk. While the theory of X-ray reverberation has been developed almost exclusively for standard thin disks, recently reverberation lags have been observed from likely super-Eddington accretion sources such as the jetted tidal disruption event Swift J1644+57. In this paper, we extend X-ray reverberation studies into the super-Eddington accretion regime, focusing on investigating the lags in the Fe K{\alpha} line region. We find that the coronal photons are mostly reflected by the fast and optically thick winds launched from super-Eddington accretion flow, and this funnel-like reflection geometry produces lag-frequency and lag-energy spectra with unique characteristics. The lag-frequency spectra exhibits a step-function like decline near the first zero-crossing point. As a result, the shape of the lag-energy spectra remains almost independent of the choice of frequency bands and linearly scales with the black hole mass for a large range of parameter spaces. Not only can these morphological differences be used to distinguish super-Eddington accretion systems from sub-Eddington systems, they are also key for constraining the reflection geometry and extracting parameters from the observed lags. When explaining the X-ray reverberation lags of Swift J1644+57, we find that the super-Eddington disk geometry is preferred over the thin disk, for which we obtain a black hole mass of 5-6 million solar masses and a coronal height around 10 gravitational radii by fitting the lag spectra to our modeling.

The jets emanating from the centers of active galactic nuclei (AGN) are among the most energetic objects in the universe. Investigating how the morphology of the jet's synchrotron emission depends on the magnetic nature of the jet's relativistic plasma is fundamental to the comparison between numerical simulations and the observed polarization of relativistic jets. Through the use of 3D relativistic magnetohydrodynamic (RMHD) jet simulations (computed using the PLUTO code) we study how the jet's synchrotron emission depends upon the morphology of the jet's magnetic field structure. Through the application of polarized radiative transfer and ray-tracing (via the RADMC-3D code) we create synthetic radio maps of the jet's total intensity as well as the linearly and circularly polarized intensity for each jet simulation. In particular, we create synthetic ray-traced images of the jet's polarized synchrotron emission when the jet carries a predominantly poloidal, helical, and toroidal magnetic field. We also explore several scaling relations in which the underlying electron power-law distribution is set proportional to: (i) the jet's thermal plasma density, (ii) the jet's internal energy density, and (iii) the jet's magnetic energy density. We find that: (i) the jet emission is edge brightened when the magnetic field is toroidal in nature and spine brightened when the magnetic field is poloidal in nature, (ii) the circularly polarized emission exhibits both negative and positive signs for the toroidal magnetic field morphology at an inclination of 45{\deg} as well as 5{\deg}, and (iii) the relativistic jet's emission is largely independent of different emission scaling relations when the ambient medium is excluded.

The Large-Sized Telescope (LST) prototype of the future Cherenkov Telescope Array (CTA) is located at the Northern site of CTA, on the Canary Island of La Palma. It is designed to provide optimal performance in the lowest part of the energy range covered by CTA, observing gamma rays down to energies of tens of GeV. The LST prototype started performing astronomical observations in November 2019 during the commissioning of the telescope and it has been taking data since then. In this contribution, we will present the tuning of the characteristics of the telescope in the Monte Carlo (MC) simulations to describe the data obtained, the estimation of its angular and energy resolution, and an evaluation of its sensitivity, both with simulations and with observations of the Crab Nebula.

The Southern Wide-field Gamma-ray Observatory (SWGO) is a proposed ground-based gamma-ray detector that will be located in the Southern Hemisphere and is currently in its design phase. In this contribution, we will outline the prospects for Galactic science with this Observatory. Particular focus will be given to the detectability of extended sources, such as gamma-ray halos around pulsars; optimisation of the angular resolution to mitigate source confusion between known TeV sources; and studies of the energy resolution and sensitivity required to study the spectral features of PeVatrons at the highest energies. Such a facility will ideally complement contemporaneous observatories in studies of high energy astrophysical processes in our Galaxy.

Fast radio bursts (FRBs) are one of the most exciting new mysteries of astrophysics. Their origin is still unknown, but recent observations seems to link them to Soft Gamma Repeaters and, in particular, to magnetar giant flares (MGFs). The recent detection of a MGF at GeV energies by the \textit{Fermi} Large Area Telescope (LAT) motivated the search for GeV counterparts to the >100 currently known FRBs. Taking advantage of more than 12 years of \textit{Fermi}-LAT data, we perform a search for gamma-ray emission from all the reported repeating and non-repeating FRBs. We analyse on different-time scales the \textit{Fermi}-LAT data of each individual source separately, including a cumulative analysis on the repeating ones. In addition, we perform the first stacking analysis at GeV energies of this class of sources in order to constrain the gamma-ray properties of the FRBs that are undetected at high energies. The stacking analysis is a powerful method that allow a possible detection from below-threshold FRBs providing important information on these objects. In this talk we present the preliminary results of our study and we discuss their implications for the predictions of gamma-ray emission from this class of sources

According to radiative models, radio galaxies are predicted to produce gamma rays from the earliest stages of their evolution onwards.The study of the high-energy emission from young radio sources is crucial for providing information on the most energetic processes associated with these sources, the actual region responsible for this emission, as well as the structure of the newly born radio jets. Despite systematic searches for young radio sources at gamma-ray energies, only a handful of detections have been reported so far. Taking advantage of more than 11 years of \textit{Fermi}-LAT data, we investigate the gamma-ray emission of 162 young radio sources (103 galaxies and 59 quasars), the largest sample of young radio sources used so far for a gamma-ray study. We analyse the \textit{Fermi}-LAT data of each individual source separately to search for a significant detection. In addition, we perform the first stacking analysis of this class of sources in order to investigate the gamma-ray emission of the young radio sources that are undetected at high energies. We report the detection of significant gamma-ray emission from 11 young radio sources, including the discovery of significant gamma-ray emission from the compact radio galaxy PKS 1007+142. Although the stacking analysis of below-threshold young radio sources does not result in a significant detection, it provides stringent upper limits to constrain the gamma-ray emission from these objects. In this talk we present the results of our study and we discuss their implications for the predictions of gamma-ray emission from this class of sources.

Upcoming studies of extrasolar gas giants will give precise insights into the composition of planetary atmospheres with the ultimate goal to link it to the formation history of the planet. Here, we investigate how drifting and evaporating pebbles that enrich the gas phase of the disk influence the chemical composition of growing and migrating gas giants. To achieve this goal, we perform semi analytical 1D models of protoplanetary disks including viscous evolution, pebble drift and evaporation to simulate the growth of planets from planetary embryos to Jupiter mass objects by the accretion of pebbles and gas while they migrate through the disk. The gas phase of the protoplanetary disk is enriched due to the evaporation of inward drifting pebbles crossing evaporation lines, leading to the accretion of large amounts of volatiles into the planetary atmosphere. As a consequence, gas accreting planets are enriched in volatiles (C, O, N) compared to refractories (e.g., Mg, Si, Fe) by up to a factor of 100, depending on the chemical species, its exact abundance and volatility as well as the disk's viscosity. A simplified model for the formation of Jupiter reveals that its nitrogen content can be explained by inward diffusing nitrogen rich vapor, implying that Jupiter does not need to form close to the N2 evaporation front as indicated by previous simulations. However, our model predicts a too low oxygen abundance for Jupiter, implying either Jupiter's migration across the water ice line or an additional accretion of solids into the atmosphere, which can also increase Jupiter's carbon abundance. The accretion of solids will increase the refractory to volatile ratio in planetary atmospheres substantially. We thus conclude that the volatile to refractory ratio can place a strong constraint on planet formation theories making it an important target for future observations.

Code comparisons are a valuable tool for the verification of supernova simulation codes and the quantification of model uncertainties. Here we present a first comparison of axisymmetric magnetohydrodynamic (MHD) supernova simulations with the CoCoNuT-FMT and Aenus-Alcar codes, which use distinct methods for treating the MHD induction equation and the neutrino transport. We run two sets of simulations of a rapidly rotating 35M gamma-ray burst progenitor model with different choices for the initial field strength, namely 10^12 G for the maximum poloidal and toroidal field in the strong-field case and 10^10 G in the weak-field case. We also investigate the influence of the Riemann solver and the resolution in CoCoNuT-FMT. The dynamics is qualitatively similar for both codes and robust with respect to these numerical details, with a rapid magnetorotational explosion in the strong-field case and a delayed neutrino-driven explosion in the weak-field case. Despite relatively similar shock trajectories, we find sizeable differences in many other global metrics of the dynamics, like the explosion energy and the magnetic energy of the proto-neutron star. Further differences emerge upon closer inspection, for example, the disk-like surface structure of the proto-neutron star proves highly sensitives to numerical details. The electron fraction distribution in the ejecta as a crucial determinant for the nucleosynthesis is qualitatively robust, but the extent of neutron- or proton-rich tails is sensitive to numerical details. Due to the complexity of the dynamics, the ultimate cause of model differences can rarely be uniquely identified, but our comparison helps gauge uncertainties inherent in current MHD supernova simulations.

Understanding the evolution of a supernova remnant shell in time is fundamental. Such understanding is critical to build reliable models of the dynamics of the supernova remnant shell interaction with any pulsar wind nebula it might contain. Here, we perform a large study of the parameter space for the one-dimensional spherically symmetric evolution of a supernova remnant, accompanying it by analytical analysis. Assuming, as is usual, an ejecta density profile with a power-law core and an envelope, and a uniform ambient medium, we provide a set of highly-accurate approximations for the evolution of the main structural features of supernova remnants, such as the reverse and forward shocks and the contact discontinuity. We compare our results with previously adopted approximations, showing that existing simplified prescriptions can easily lead to large errors. In particular, in the context of pulsar wind nebulae modelling, an accurate description for the supernova remnant reverse shock is required. We also study in depth the self-similar solutions for the initial phase of evolution, when the reverse shock propagates through the envelope of the ejecta. Since these self-similar solutions are exact, but not fully analytical, we here provide highly-accurate approximations as well.

We present a novel Probabilistic Flux Variation Gradient (PFVG) approach to disentangle the active galactic nuclei (AGN) and host-galaxy contributions in the context of photometric reverberation mapping (PRM) of AGN. We explored the ability to recover the fractional contribution, in a model-independent way, by using simultaneously the entire set of light curves obtained through different filters and photometric apertures. The method is based on the observed bluer when brighter phenomena attributed to the superimposition of a two component structure; the red-host galaxy, which is constant in time, and the varying blue AGN. We describe the PFVG mathematical formalism and demonstrate its performance using simulated light curves and available PRM observations. The new probabilistic approach is able to recover host-galaxy fluxes within 1\% precision, as long as the light curves do not show a significant contribution from time-delays. This represent a significant improvement with respect to previous applications of the traditional FVG method on PRM data. The proposed PFVG provides an efficient and accurate way to disentangling the AGN and host-galaxy luminosities in PRM monitoring data. The method will be specially helpful in the case of large upcoming photometric survey telescopes like the LSST. Finally, we have made the algorithms freely available as part of our Julia PFVG package.

(Abridged) Complex organic molecules (COMs) are commonly detected in and near star-forming regions. However, the dominant process in the release of these COMs from the icy grains - where they predominately form - to the gas phase is still an open question. We investigate the origin of COM emission in a protostellar source, CygX-N30, through high-angular-resolution interferometric observations over a continuous broad frequency range. We used 32 GHz Submillimeter Array observations with continuous frequency coverage from 329 to 361 GHz at an angular resolution of ~1" to do a line survey and obtain a chemical inventory of the source. The line emission was used to determine column densities and excitation temperatures for the COMs. We mapped out the intensity distribution of the different species and identified approximately 400 lines that can be attributed to 29 different molecular species and their isotopologues. We find that the molecular peak emission is along a linear gradient, coinciding with the axis of red- and blueshifted H2CO and CS emission. Chemical differentiation is detected along this gradient, with the O-bearing molecular species peaking towards one component of the system and the N- and S-bearing species peaking towards the other. The inferred column densities and excitation temperatures are compared to other sources where COMs are abundant. The origin of the observed COM emission is probably a combination of the young stellar sources along with accretion of infalling material onto a disc-like structure surrounding a young protostar. The low D/H ratio observed (<0.1%) likely reflects a pre-stellar phase where COMs formed on the ices at warm temperatures (~ 30 K), with inefficient deuterium fractionation. The observations and results presented here demonstrate the importance of good frequency coverage and high angular resolution when disentangling the origin of COM emission.

We study the evolution of nitrogen in the Galactic halo, thick disc, thin disc and bulge by comparing detailed chemical evolution models with recent observations. The models used in this work have already been constrained to explain the abundance patterns of $\alpha$-elements and the metallicity distribution functions of halo, disc and bulge stars; here, we adopt them to investigate the origin and evolution of N in the different Galactic components. First, we consider different sets of yields and study the importance of the various channels proposed for N production. Secondly, we apply the reference models to study the evolution of both the Galactic discs and bulge. We conclude that: i) primary N produced by rotating massive stars is required to reproduce the plateau in log(N/O) and [N/Fe] ratios at low metallicity, as well as the secondary and primary production from low- and intermediate-mass stars to reproduce the data of the thin disc. ii) The parallel model can provide a good explanation of the evolution of N abundance in the thick and thin discs; we confirm that the thick disc has evolved much faster than the thin disc, in agreement with the results from the abundance patterns of other chemical elements. iii) Finally, we present new model predictions for N evolution in the Galactic bulge, and we show that the observations in bulge stars can be explained if massive stars rotate fast during the earliest phases of Galactic evolution, in agreement with findings from the abundance pattern of carbon.

Recent improvements to GPU hardware and the symplectic N-body code GENGA allow for unprecedented resolution in simulations of planet formation. In this paper, we report results from high-resolution N-body simulations of terrestrial planet formation that are mostly direct continuation of our previous 10 Myr simulations (Woo et al. 2021a) until 150 Myr. By assuming that Jupiter and Saturn have always maintained their current eccentric orbits (EJS), we are able to achieve a reasonably good match to the current inner solar system architecture. However, due to the strong radial mixing that occurs in the EJS scenario, it has difficulties in explaining the known isotopic differences between bodies in the inner solar system, most notably between Earth and Mars. On the other hand, assuming initially circular orbits for Jupiter and Saturn (CJS) can reproduce the observed low degree of radial mixing in the inner solar system, while failing to reproduce the current architecture of the inner solar system. These outcomes suggest a possible paradox between dynamical structure and cosmochemical data for the terrestrial planets within the classical formation scenario.

We present Lyman continuum (LyC) radiation escape fraction $f_{\rm{esc}}$ measurements for 183 spectroscopically confirmed star-forming galaxies in the redshift range $3.11 < z < 3.53$ in the \textit{Chandra} Deep Field South. We use ground-based imaging to measure $f_{\rm{esc}}$, and use ground- and space-based photometry to derive galaxy physical properties using spectral energy distribution (SED) fitting. We additionally derive [O\,\textsc{iii}]\,+\,H$\beta$ equivalent widths (that fall in the observed $K$ band) by including nebular emission in the SED fitting. After removing foreground contaminants, we report the discovery of 11 new candidate LyC leakers, with absolute LyC escape fractions, $f_{\rm{esc}}$ in the range $0.07-0.52$. Most galaxies in our sample ($\approx94\%$) do not show any LyC leakage, and we place $1\sigma$ upper limits of $f_{\rm{esc}} < 0.07$ through weighted averaging, where the Lyman-break selected galaxies have $f_{\rm{esc}} < 0.07$ and `blindly' discovered galaxies with no prior photometric selection have $f_{\rm{esc}} < 0.10$. We additionally measure $f_{\rm{esc}} < 0.09$ for extreme emission line galaxies in our sample with rest-frame [O\,\textsc{iii}]\,+\,H$\beta$ equivalent widths $>300$\,\AA. For the candidate LyC leakers, we do not find a strong dependence of $f_{\rm{esc}}$ on their stellar masses and/or specific star-formation rates, and no correlation between $f_{\rm{esc}}$ and EW$_0$([O\,\textsc{iii}]\,+\,H$\beta$). We suggest that this lack of correlations may be explained by viewing angle and/or non-coincident timescales of starburst activity and periods of high $f_{\rm{esc}}$. Alternatively, escaping radiation may predominantly occur in highly localised star-forming regions, thereby obscuring any global trends with galaxy properties. Both hypotheses have important consequences for models of reionisation.

Sudden changes in the internal structure of stars, placed at the interface between convective and radiative regions, regions of partial ionisation, or between layers that have acquired different chemical composition as a result of nuclear burning, often produce specific signatures in the stars oscillation spectra. Through the study of these signatures one may gain information on the physical processes that shape the regions that produce them, including diffusion and chemical mixing beyond the convectively unstable regions, as well as information about the helium content of stars. In this talk, I will review important theoretical and observational efforts conducted over the years towards this goal. I will emphasise the potential offered by the study of acoustic, gravity, and mixed modes observed in stars of different mass and evolutionary stages, at a time when space-based data is allowing us to build on the knowledge gained from the study of the sun and white dwarfs, where these efforts have long been undertaken, extending the methods developed to stars across the HR diagramme.

The distribution of ejecta in young supernova remnants offers a powerful observational probe of their explosions and progenitors. Here we present a 3D reconstruction of the ejecta in SNR 0540-69.3, which is an O-rich remnant with a pulsar wind nebula located in the LMC. We use observations from VLT/MUSE to study H\beta, [O III] \lambda \lambda 4959, 5007, H\alpha, [S II] \lambda \lambda 6717, 6731, [Ar III] \lambda 7136 and [S III] \lambda 9069 emission lines. This is complemented by 2D spectra from VLT/X-shooter, which also cover [O II] \lambda \lambda 3726, 3729 and [Fe II] \lambda 12567. We identify three main emission components: (i) Clumpy rings in the inner nebula (<1000 km/s) with similar morphologies in all lines; (ii) Faint extended [O III] emission dominated by an irregular ring-like structure with radius ~1600 km/s and inclination ~40 \dg, but with maximal velocities reaching ~3000 km/s; and (iii) A blob of H\alpha and H\beta located southeast of the pulsar at velocities ~1500-3500 km/s. We analyze the geometry using a clump-finding algorithm and use the clumps in the [O III] ring to estimate an age of 1146 \pm 116 years. The observations favor an interpretation of the [O III] ring as ejecta, while the origin of the H-blob is more uncertain. An alternative explanation is that it is the blown-off envelope of a binary companion. From the detection of Balmer lines in the innermost ejecta we confirm that SNR 0540 was a Type II supernova and that hydrogen was mixed down to low velocities in the explosion.

The High Energy Stereoscopic System (H.E.S.S.) is one of the currently operating Imaging Atmospheric Cherenkov Telescopes. H.E.S.S. operates in the broad energy range from a few tens of GeV to more than 50 TeV reaching its best sensitivity around 1 TeV. In this contribution, we present an analysis technique, which is optimised for the detection at the highest energies accessible to H.E.S.S. and aimed to improve the sensitivity above 10 TeV. It includes the employment of improved event direction reconstruction and gamma-hadron separation. For the first time, also extensive air showers with event offsets up to 4.5$^{\circ}$ from the camera centre are considered in the analysis, thereby increasing the effective Field-of-View of H.E.S.S. from 5$^{\circ}$ to 9$^{\circ}$. Key performance parameters of the new high-energy analysis are presented and its applicability demonstrated for representative hard-spectrum sources in the Milky Way.

We use 15 years of $\gamma$-ray data from INTEGRAL/SPI in a refined investigation of the morphology of the Galactic bulge positron annihilation signal. Our spatial analysis confirms that the signal traces the old stellar population in the bulge and reveals for the first time that it traces the boxy bulge and nuclear stellar bulge. Using a 3D smoothing kernel, we find that the signal is smeared out over a characteristic length scale of $150 \pm 50\,$pc, suggesting either annihilation in situ at astrophysical sources kicked at formation or positron propagation away from sources. The former is disfavoured by its requiring kick velocities different between the Galactic nucleus ($\gtrsim 50\,\mathrm{km\,s^{-1}}$) and wider bulge ($\lesssim 15\,\mathrm{km\,s^{-1}}$) source. Positron propagation prior to annihilation can explain the overall phenomenology of the 511 keV signal for positrons injection energies $\lesssim 1.4\,$MeV, suggesting a nucleosynthesis origin.

Observations with the current generation of very-high-energy gamma-ray telescopes have revealed an astonishing variety of particle accelerators in the Milky Way, such as supernova remnants, pulsar wind nebulae, and binary systems. The upcoming Cherenkov Telescope Array (CTA) will be the first instrument to enable a survey of the entire Galactic plane in the energy range from a few tens of GeV to 300 TeV with unprecedented sensitivity and improved angular resolution. In this contribution we will revisit the scientific motivations for the survey, proposed as a Key ScienceProject for CTA. We will highlight recent progress, including improved physically-motivated models for Galactic source populations and interstellar emission, advance on the optimization of the survey strategy, and the development of pipelines to derive source catalogues tested on simulated data. Based on this, we will provide a new forecast on the properties of the sources thatCTA will detect and discuss the expected scientific return from the study of gamma-ray source populations.

The structure and dynamics of the magnetospheres of the outer planets, particularly Saturn and Jupiter, have been explored both through remote and in-situ observations. Interpreting these observations often necessitates simultaneous knowledge of the solar-wind conditions impinging on the magnetosphere. Without an available upstream monitor, solar-wind context is typically provided using models initiated with either the output of magnetogram-constrained coronal models or, more commonly, in-situ observations from 1 AU. While 1-AU observations provide a direct measure of solar wind conditions, they are single-point observations and thus require interpolation to provide inputs to outer-heliosphere solar-wind models. In this study we test the different interpolation methods using synthetic 1-AU observations of time-evolving solar-wind structure. The simplest method is "corotation", which assumes solar-wind structure is steady state and rotates with the Sun. This method of reconstruction produces discontinuities in the solar-wind inputs as new observations become available. This can be reduced by corotating both back and forward in time, but this still introduces large errors in the magnitude and timing of solar wind streams. We show how the dynamic time warping (DTW) algorithm can provide around an order-of-magnitude improvement in solar-wind inputs to outer-heliosphere model from in-situ observations near 1 AU. This is intended to build the foundation for further work demonstrating and validating methods to improve inner-boundary conditions to outer-heliosphere solar-wind models, including dealing with solar wind transients and quantifying the improvements at Saturn and Jupiter.

The presence of current sheet instabilities, such as the tearing mode instability, are needed to account for the observed rate of energy release in solar flares. Insights into these current sheet dynamics can be revealed by the behaviour of flare ribbon substructure, as magnetic reconnection accelerates particles down newly reconnected field lines into the chromosphere to mark the flare footpoints. Behaviour in the ribbons can therefore be used to probe processes occurring in the current sheet. In this study, we use high-cadence (1.7 s) IRIS Slit Jaw Imager observations to probe for the growth and evolution of key spatial scales along the flare ribbons - resulting from dynamics across the current sheet of a small solar flare on December 6th 2016. Combining analysis of spatial scale growth with Si IV non-thermal velocities, we piece together a timeline of flare onset for this confined event, and provide evidence of the tearing-mode instability triggering a cascade and inverse cascade towards a power spectrum consistent with plasma turbulence.

We present the results of a search for outer Solar System objects in the full six years of data (Y6) from the Dark Energy Survey (DES). The DES covered a contiguous $5000$ deg$^2$ of the southern sky with $\approx 80,000$ $3$ deg$^2$ exposures in the $grizY$ optical/IR filters between 2013 and 2019. This search yielded 815 trans-Neptunian objects (TNOs), one Centaur and one Oort cloud comet, with 461 objects reported for the first time in this paper. We present methodology that builds upon our previous search carried out on the first four years of data. Here, all DES images were reprocessed with an improved detection pipeline that leads to an average completeness gain of 0.47 mag per exposure, as well as an improved transient catalog production and optimized algorithms for linkage of detections into orbits. All objects were verified by visual inspection and by computing the sub-threshold significance, the total signal-to-noise ratio in the stack of images in which the object's presence is indicated by the orbit fit, but no detection was reported. This yields a highly pure catalog of TNOs complete to $r \approx 23.8$ mag and distances $29<d<2500$ au. The Y6 TNOs have minimum (median) of 7 (12) distinct nights' detections and arcs of 1.1 (4.2) years, and will have $grizY$ magnitudes available in a further publication. We present software for simulating our observational biases that enable comparisons of population models to our detections. Initial inferences demonstrating the statistical power of the DES catalog are: the data are inconsistent with the CFEPS-L7 model for the classical Kuiper Belt; the 16 ``extreme'' TNOs ($a>150$ au, $q>30$ au) are consistent with the null hypothesis of azimuthal isotropy; and non-resonant TNOs with $q>38$ au, $a>50$ au show a highly significant tendency to be sunward of the major mean motion resonances, whereas this tendency is not present for $q<38$ au.

Measuring the total neutrino mass is one of the most exciting opportunities available with next-generation cosmological data sets. We study the possibility of detecting the total neutrino mass using large-scale clustering in 21cm intensity mapping and photometric galaxy surveys, together with CMB information. We include the scale-dependent halo bias contribution due to the presence of massive neutrinos, and use a multi-tracer analysis in order to reduce cosmic variance. The multi-tracer combination of an SKAO-MID 21cm intensity map with Stage~4 CMB lensing dramatically shrinks the uncertainty on total neutrino mass to $\sigma(M_\nu) \simeq 45\,$meV, using only linear clustering information ($k_{\rm max} = 0.1\, h/$Mpc) and without a prior on optical depth. When we add to the multi-tracer the clustering information expected from LSST, the forecast is $\sigma(M_\nu) \simeq 12\,$meV.

Binary stars play a vital role in astrophysical research, as a good fraction of stars are in binaries. Binary fraction (BF) is known to change with stellar mass in the Galactic field, but such studies in clusters require binary identification and membership information. Here, we estimate the total and spectral-type-wise high mass-ratio (HMR) BF ($f^{0.6}$) in 23 open clusters using unresolved binaries in color-magnitude diagrams using \textit{Gaia} DR2 data. We introduce the segregation index (SI) parameter to trace mass segregation of HMR (total and mass-wise) binaries and the reference population. This study finds that in open clusters, (1) HMR BF for the mass range 0.4--3.6 Msun (early M to late B type) has a range of 0.12 to 0.38 with a peak at 0.12--0.20, (2) older clusters have a relatively higher HMR BF, (3) the mass-ratio distribution is unlikely to be a flat distribution and BF(total) $\sim$ (1.5 to 2.5) $\times f^{0.6}$, (4) a decreasing BF(total) from late B-type to K-type, in agreement with the Galactic field stars, (5) older clusters show radial segregation of HMR binaries, (6) B and A/F type HMR binaries show radial segregation in some young clusters suggesting a primordial origin. This study will constrain the initial conditions and identify the major mechanisms that regulate binary formation in clusters. Primordial segregation of HMR binaries could result from massive clumps spatially segregated in the collapse phase of the molecular cloud.

Reticulum II (Ret II) is a satellite galaxy of the Milky Way and presents a prime target to investigate the nature of dark matter (DM) because of its high mass-to-light ratio. We evaluate a dedicated INTEGRAL observation campaign data set to obtain $\gamma$-ray fluxes from Ret II and compare those with expectations from DM. Ret II is not detected in the $\gamma$-ray band 25--8000 keV, and we derive a flux limit of $\lesssim 10^{-8}\,\mathrm{erg\,cm^{-2}\,s^{-1}}$. The previously reported 511 keV line is not seen, and we find a flux limit of $\lesssim 1.7 \times 10^{-4}\,\mathrm{ph\,cm^{-2}\,s^{-1}}$. We construct spectral models for primordial black hole (PBH) evaporation and annihilation/decay of particle DM, and subsequent annihilation of positrons produced in these processes. We exclude that the totality of DM in Ret II is made of a monochromatic distribution of PBHs of masses $\lesssim 8 \times 10^{15}\,\mathrm{g}$. Our limits on the velocity-averaged DM annihilation cross section into $e^+e^-$ are $\langle \sigma v \rangle \lesssim 5 \times 10^{-28} \left(m_{\rm DM} / \mathrm{MeV} \right)^{2.5}\,\mathrm{cm^3\,s^{-1}}$. We conclude that analysing isolated targets in the MeV $\gamma$-ray band can set strong bounds on DM properties without multi-year data sets of the entire Milky Way, and encourage follow-up observations of Ret II and other dwarf galaxies.

We study an impact of self-interacting bosonic dark matter (DM) on various observable properties of neutron stars (NSs). The analysis is performed for asymmetric DM with masses from few MeV to GeV, the self-coupling constant of order $\mathcal{O}(1)$ and various DM fractions. Allowing a mixture between DM and baryonic matter, the formation of a dense DM core or an extended dark halo have been explored. We find that both distribution regimes crucially depend on the mass and fraction of DM for sub-GeV boson masses in the strong coupling regime. From the combined analysis of the mass-radius relation and the tidal deformability of compact stars including bosonic DM, we set a stringent constraint on DM fraction. We conclude that observations of 2$M_{\odot}$ NSs together with $\Lambda_{1.4}\leq580$ constraint, set by LIGO/Virgo Collaboration, favour sub-GeV DM particles with low fractions below $\sim 5 \%$.

Primordial magnetic fields can change the recombination history of the universe by inducing clumping in the baryon density at small scales. They were recently proposed as a candidate model to relieve the Hubble tension. We investigate the consistency of the constraints on a clumping factor parameter $b$ in a simplistic model, using the latest CMB data from Planck, ACT DR4 and SPT-3G 2018. For the combined CMB data alone, we find no evidence for clumping being different from zero, though when adding a prior on $H_0$ based on the latest distance-ladder analysis of the SH0ES team, we report a weak detection of $b$. Our analysis of simulated datasets shows that ACT DR4 has more constraining power with respect to SPT-3G 2018 due to the degeneracy breaking power of the TT band powers (not included in SPT). Simulations also suggest that the TE,EE power spectra of the two datasets should have the same constraining power. However, the ACT DR4 TE,EE constraint is tighter than expectations, while the SPT-3G 2018 one is looser. While this is compatible with statistical fluctuations, we explore systematic effects which may account for such deviations. Overall, the ACT results are only marginally consistent with Planck or SPT-3G, at the $2-3\sigma$ level within $\Lambda$CDM+$b$ and $\Lambda$CDM, while Planck and SPT-3G are in good agreement. Combining the CMB data together with BAO and SNIa provides an upper limit of b<0.4 at 95\% c.l. (b<0.5 without ACT). Adding a SH0ES-based prior on the Hubble constant gives $b = 0.31^{+0.11}_{-0.15}$ and $H_0=69.28 \pm 0.56$ km/s/Mpc ($b = 0.41^{+0.14}_{-018}$ and $H_0=69.70 \pm 0.63$ km/s/Mpc without ACT). Finally, we forecast constraints on $b$ for the full SPT-3G survey, Simons Observatory, and CMB-S4, finding improvements by factors of 1.5 (2.7 with Planck), 5.9 and 7.8, respectively, over Planck alone.

We study for the first time the possibility of probing long-range fifth forces utilizing asteroid astrometric data, via the fifth force-induced orbital precession. We examine nine Near-Earth Object (NEO) asteroids whose orbital trajectories are accurately determined via optical and radar astrometry. Focusing on a Yukawa-type potential mediated by a new gauge field (dark photon) or a baryon-coupled scalar, we estimate the sensitivity reach for the fifth-force coupling strength and mediator mass in the mass range $m \simeq 10^{-21}-10^{-15}\,{\rm eV}$. Our estimated sensitivity is comparable to leading limits from torsion balance experiments, potentially exceeding these in a specific mass range. The fifth forced-induced precession increases with the orbital semi-major axis in the small $m$ limit, motivating the study of objects further away from the Sun. We discuss future exciting prospects for extending our study to more than a million asteroids (including NEOs, main-belt asteroids, Hildas, and Jupiter Trojans), as well as trans-Neptunian objects and exoplanets.

Nature's most powerful high-energy sources are capable of accelerating particles to high energy and radiate it away on extremely short timescales, even shorter than the light crossing time of the system. It is yet unclear what physical processes can produce such an efficient acceleration, despite the copious radiative losses. By means of radiative particle-in-cell simulations, we show that magnetically dominated turbulence in pair plasmas subject to strong synchrotron cooling generates a nonthermal particle spectrum with a hard power-law range (slope $p \sim 1$) within a few eddy turnover times. Low pitch-angle particles can significantly exceed the nominal radiation-reaction limit, before abruptly cooling down. The particle spectrum becomes even harder ($p < 1$) over time owing to particle cooling with an energy-dependent pitch-angle anisotropy. The resulting synchrotron spectrum is hard ($\nu F_\nu \propto \nu^s$ with $s \sim 1$). Our findings have important implications for understanding the nonthermal emission from high-energy astrophysical sources, most notably the prompt phase of gamma-ray bursts and gamma-ray flares from the Crab nebula.

We present an explicit string realisation of a cosmological inflationary scenario we proposed recently within the framework of type IIB flux compactifications in the presence of three magnetised D7-brane stacks. Inflation takes place around a metastable de Sitter vacuum. The inflaton is identified with the volume modulus and has a potential with a very shallow minimum near the maximum. Inflation ends due to the presence of "waterfall" fields that drive the evolution of the Universe from a nearby saddle point towards a global minimum with tuneable vacuum energy describing the present state of our Universe.

We derive supernova (SN) bounds on muon-philic bosons, taking advantage of the recent emergence of muonic SN models. Our main innovations are to consider scalars $\phi$ in addition to pseudoscalars $a$ and to include systematically the generic two-photon coupling $G_{\gamma\gamma}$ implied by a muon triangle loop. This interaction allows for Primakoff scattering and radiative boson decays. The globular-cluster bound $G_{\gamma\gamma}<0.67\times10^{-10}~{\rm GeV}^{-1}$ derived for axion-like particles carries over to the muonic Yukawa couplings as $g_a<3.1\times10^{-9}$ and $g_\phi< 4.6\times10^{-9}$ for $m_{a,\phi}\lesssim 100$ keV, so SN arguments become interesting mainly for larger masses. If bosons escape freely from the SN core the main constraints originate from SN1987A $\gamma$ rays and the diffuse cosmic $\gamma$-ray background. The latter allows at most $10^{-4}$ of a typical total SN energy of $E_{\rm SN}\simeq3\times10^{53}$erg to show up as $\gamma$ rays, for $m_{a,\phi}\gtrsim 100$keV implying $g_a \lesssim 0.9\times10^{-10}$ and $g_\phi \lesssim 0.4\times10^{-10}$. In the trapping regime the bosons emerge as quasi-thermal radiation from a region near the neutrino sphere and match $L_\nu$ for $g_{a,\phi}\simeq 10^{-4}$. However, the $2\gamma$ decay is so fast that all the energy is dumped into the surrounding progenitor-star matter, whereas at most $10^{-2}E_{\rm SN}$ may show up in the explosion. To suppress boson emission below this level we need yet larger couplings, $g_{a}\gtrsim 2\times10^{-3}$ and $g_{\phi}\gtrsim 4\times10^{-3}$. Muonic scalars can explain the muon magnetic-moment anomaly for $g_{\phi}\simeq 0.4\times10^{-3}$, a value hard to reconcile with SN physics despite the uncertainty of the explosion-energy bound. For generic axion-like particles, this argument covers the "cosmological triangle" in the $G_{a\gamma\gamma}$--$m_a$ parameter space.

Self-interaction among the neutrinos in the early Universe has been proposed as a solution to the Hubble tension, a discrepancy between the measured values of the Hubble constant from CMB and low-redshift data. However, flavor-universal neutrino self-interaction is highly constrained by Big Bang Nucleosynthesis and other laboratory experiments such as, tau and K-meson decay, double-neutrino beta decay etc. We study the cosmology when only one or two neutrino states are self-interacting. Such flavor-specific interactions are less constrained by the laboratory experiments. Lastly, we address the feasibility of resolving the Hubble tension within the framework of such flavor-specific neutrino self-interaction.

Production of domain walls and string-like solitons in the model with two real scalar fields and potential with at least one saddle point and a local maximum is considered. The model is regarded as 2-dimensional spatial slices of 3-dimensional entire structures. It is shown that, in the early Universe, both types of solitons may appear. In addition, the qualitative estimate of the domain walls and strings formation probabilities is presented. It is found that the probability of the formation of string-like solitons is suppressed compared to that of domain walls.

Ultralight bosons, as important candidates of dark matter, can condense around spinning black holes (BHs) to form long-lived ``boson clouds'' due to superradiance instability. The boson-BH system can be observed through gravitational wave detection and may become a new window to find traces of ultralight bosons. In this letter we explore the effects on the superradiant instability of BHs from the near-horizon microstructure. By introducing the reflection parameter near a BH horizon, we derived analytical results on the corrections to both energy levels of bosonic cloud and its characteristic frequencies of superradiance instability. Our results imply that the evolution of a boson-BH system and gravitational waves it emits would be influenced by the near-horizon physics of a BH.

A correspondence between fluctuations of non-minimally coupled scalar fields and that of an effective fluid with heat flux and anisotropic stresses, is shown. Though the correspondence between respective stress tensors of scalar fields and fluids is known and widely used in literature, the fluctuations in the two cases still await a formal correspondence and are open to investigation in all details. Using results obtained in the newly established theory of semiclassical stochastic gravity which focuses on the fluctuations of the quantum stress tensor, we show new relations in this regard. This development is expected to give insight to the mesoscopic phenomena for gravitating systems, and enable backreaction studies of the fluctuations on the perturbations of astrophysical objects. Such a development is aimed to enhance the perturbative analysis for cosmological spacetimes and astrophysical objects specifically in the decoherence limit. A kinetic theory, which can be based on stochastic fluctuations vs particle picture in curved spacetime may find useful insights from such correspondences in future work.

The computational cost of inspiral and merger simulations for black-hole binaries increases in inverse proportion to the square of the mass ratio $q:=m_2/m_1\leq 1$. One factor of $q$ comes from the number of orbital cycles, which is proportional to $1/q$, and another is associated with the required number of time steps per orbit, constrained (via the Courant-Friedrich-Lewy condition) by the need to resolve the two disparate length scales. This problematic scaling makes simulations progressively less tractable at smaller $q$. Here we propose and explore a method for alleviating the scale disparity in simulations with mass ratios in the intermediate astrophysical range ($10^{-4} \lesssim q\lesssim 10^{-2}$), where purely perturbative methods may not be adequate. A region of radius much larger than $m_2$ around the smaller object is excised from the numerical domain, and replaced with an analytical model approximating a tidally deformed black hole. The analytical model involves certain a priori unknown parameters, associated with unknown bits of physics together with gauge-adjustment terms; these are dynamically determined by matching to the numerical solution outside the excision region. In this paper we develop the basic idea and apply it to a toy model of a scalar charge in a circular geodesic orbit around a Schwarzschild black hole, solving for the massless Klein-Gordon field in a 1+1D framework. Our main goal here is to explore the utility and properties of different matching strategies, and to this end we develop two independent implementations, a finite-difference one and a spectral one. We discuss the extension of our method to a full 3D numerical evolution and to gravity.

We study in some detail an interacting cosmological model based on two canonical scalar fields starting from a Lagrangian description. Contrary to other more phenomenological approaches where non-relativistic matter and dark energy are cosmological perfect fluids, and where a source term is added by hand at the level of the continuity equations, here within Einstein's theory we model the dark sector, which dominates the evolution of the universe, as two minimally coupled scalar fields, out of which the first play the role of dark matter and the second plays the role of dark energy. We compute both the deceleration parameter and the distance modulus versus red-shift, and we demonstrate that the model is capable of explaining the current cosmic acceleration. We find that a negative coupling constant implies two distinctive features, which comprise two robust predictions of the model studied in the present work, and which are the following: a) a transient acceleration phase, and b) a collapsing Universe, or in other words an initial expansion which is followed by a contraction leading eventually to a Big-Crunch.

Kaluza-Klein (KK) axions appear in theories with extra dimensions as higher mass, significantly shorter lifetime, excitations of the Peccei-Quinn axion. When produced in the Sun, they would remain gravitationally trapped in the solar system, and their decay to a pair of photons could provide an explanation of the solar corona heating problem. A low-density detector would discriminate such a signal from the background, by identifying the separation of the interaction point of the two photons. The NEWS-G collaboration uses large volume Spherical Proportional Counters, gas-filled metallic spheres with a spherical anode in their centre. After observation of a single axion-like event in a 42 day long run with the SEDINE detector, a $90\%$ C.L. upper limit of $g_{a\gamma\gamma}<7.76\cdot10^{-13}\,GeV^{-1}$ is set on the axion-photon coupling for a KK axion density on Earth of $n_{a}=4.07\cdot10^{13}\,m^{-3}$ and two extra dimensions of size $R = 1\,eV^{-1}$.

In this paper we infer the mass composition of the ultra high energy cosmic rays (UHECRs) from measurements of $X_{\rm max}$ distributions recorded at the Pierre Auger (2014) and Telescope Array (TA) (2016) Observatories, by fitting them with all possible combinations of Monte Carlo (MC) templates from a large set of primary species (p, He, C, N, O, Ne, Si and Fe), as predicted by EPOS-LHC, QGSJETII-04 and Sibyll 2.1 hadronic interaction models. We use the individual fractions of nuclei reconstructed from one experiment in each energy interval to build equivalent MC $X_{\rm max}$ distributions, which we compare with the experimental $X_{\rm max}$ distributions of the other experiment, applying different statistical tests of compatibility. The results obtained from both experiments confirm that the mass composition of the UHECRs is dominated ($\gtrsim$$70\%$) by protons and He nuclei {in the energy range investigated $\lg E (\rm eV)$ = [17.8--19.3] (Auger) and $\lg E \rm (eV)$ = [18.2--19.0] (TA).} % We changed the minus sign into endash in the paper, please confirm. OK The indirect comparisons between the $X_{\rm max}$ distributions recorded by the two experiments show that the degree of compatibility of the two datasets is good, even excellent in some high energy intervals, especially above the ankle ($\lg E (\rm eV) \sim 18.7$). However, our study reveals that, at low energies, further effort in data analysis is required in order to harmonize the results of the two experiments.