Papers for Friday, Dec 17 2021

Galaxy evolution is driven by a variety of physical processes which are predicted to proceed at different rates for different dark matter haloes and environments across cosmic times. A record of this evolution is preserved in galaxy stellar populations, which we can access using absorption-line spectroscopy. Here we explore the large LEGA-C survey (DR3) to investigate the role of the environment and stellar mass on stellar populations at z~0.6-1.0 in the COSMOS field. Leveraging the statistical power and depth of LEGA-C, we reveal significant gradients in D4000 and H-delta equivalent widths (EWs) distributions over the stellar mass vs environment 2D spaces for the massive galaxy population (M>10^10 M$_{\odot}$) at z~0.6-1.0. D4000 and H-delta EWs primarily depend on stellar mass, but they also depend on environment at fixed stellar mass. By splitting the sample into centrals and satellites, and in terms of star-forming galaxies and quiescent galaxies, we reveal that the significant environmental trends of D4000 and H-delta EW when controlling for stellar mass are driven by quiescent galaxies. Regardless of being centrals or satellites, star-forming galaxies reveal D4000 and H-delta EWs which depend strongly on their stellar mass and are completely independent of the environment at 0.6<z<1.0. The environmental trends seen for satellite galaxies are fully driven by the trends that hold only for quiescent galaxies, combined with the strong environmental dependency of the quiescent fraction at fixed stellar mass. Our results are consistent with recent predictions from simulations that point towards massive galaxies forming first in over-densities or the most compact dark matter haloes.

The AMS-02 experiment has provided high-precision measurements of several cosmic-ray (CR) species. The achieved percent-level accuracy gives access to small spectral differences among the different species and, in turn, this allows scrutinizing the universality of CR acceleration, which is expected in the standard scenario of CR shock acceleration. While pre-AMS-02 data already indicated a violation of the universality between protons and helium, it is still an open question if at least helium and heavier nuclei can be reconciled. To address this issue, we performed a joint analysis using the AMS-02 CR measurements of antiprotons, protons, helium, helium 3, boron, carbon, nitrogen, and oxygen. We explore two competing propagation scenarios, one with a break in the diffusion coefficient at a few GVs and no reacceleration, and another one with reacceleration and with a break in the injection spectra of primaries. Furthermore, we explicitly consider the impact of the uncertainties in the nuclear production cross-sections of secondaries by including nuisance parameters in the fit. The resulting parameter space is explored with the help of Monte Carlo methods. We find that, contrary to the naive expectation, in the standard propagation scenarios CR universality is violated also for He, on the one hand, and C, N, and O, on the other hand, i.e., different injection slopes are required to explain the observed spectra. As an alternative, we explore further propagation scenarios, inspired by non-homogeneous diffusion, which might save universality. Finally, we also investigate the universality of CR propagation, i.e., we compare the propagation properties inferred using only light nuclei ($\bar{p}$, p, He, $^3$He) with the ones inferred using only heavier nuclei (B, C, N, O).

Ill-posed inverse problems are common in astronomy, and their solutions are unstable with respect to noise in the data. Solutions of such problems are typically found using two classes of methods: parametrization and fitting the data against some predefined function or a solution with a non-parametrical function using regularization. Here we are focusing on the latter non-parametric approach applied for the recovery of complex stellar line-of-sight velocity distribution (LOSVD) from the observed galaxy spectra. Development of such an approach is crucial for galaxies hosting multiple kinematically misaligned stellar components, such as 2 stellar counter-rotating disks, thin and thick disks, kinematically decoupled cores, and others. Stellar LOSVD recovery from the observed galaxy spectra is equivalent to a deconvolution and can be solved as a linear inverse problem. To overcome its ill-posed nature we apply smoothing regularization. Searching for an optimal degree of smoothing regularization is a challenging part of this approach. Here we present a non-parametric fitting technique, discuss its potential caveats, perform numerous tests based on synthetic mock spectra, and show real-world application to MaNGA spectral data cubes and some long-slit spectra of stellar counter-rotating galaxies. GitHub repository: https://github.com/gasymovdf/sla

The presence of scalar fields with non-minimal gravitational interactions of the form $\xi |\phi|^2 R$ may have important implications for the physics of the early universe. While many studies solve the dynamics of non-minimally coupled scalars in the Einstein frame, where gravity is simply described by the Einstein-Hilbert action, we instead propose a procedure to solve the dynamics directly in the original Jordan frame where the non-minimal couplings are maintained explicitly. Our algorithm can be applied to scenarios that include minimally coupled fields and an arbitrary number of non-minimally coupled scalars, with the expansion of the universe sourced by all fields present. This includes situations when the dynamics become fully inhomogeneous, fully non-linear (due to e.g.~backreaction or mode rescattering effects), and/or when the expansion of the universe is dominated by non-minimally coupled species. As an example, we study geometric preheating with a non-minimally coupled scalar spectator field when the inflaton oscillates following the end of inflation. In the future, our technique may be used to shed light on aspects of the equivalence of the Jordan and Einstein frames at the quantum level.

The analysis of cosmological galaxy surveys requires realistic simulations for their interpretation. Forward modelling is a powerful method to simulate galaxy clustering without the need for an underlying complex model. This approach requires fast cosmological simulations with a high resolution and large volume, to resolve small dark matter halos associated to single galaxies. In this work, we present fast halo and subhalo clustering simulations based on the Lagrangian perturbation theory code PINOCCHIO, which generates halos and merger trees. The subhalo progenitors are extracted from the merger history and the survival of subhalos is modelled. We introduce a new fitting function for the subhalo merger time, which includes an additional dependence on subhalo mass. The spatial distribution of subhalos within their hosts is modelled using a number density profile. We compare our simulations with the halo finder Rockstar applied to the full N-body code GADGET-2. We find a good agreement for the number ratio of subhalos, for halo masses down to $5.7 \cdot 10^9$ M$_\odot$/h. The subhalo mass function and the correlation function of halos and subhalos are also in good agreement. We investigate the effect of the chosen number density profile on the resulting subhalo clustering. Our simulation is approximate yet realistic and significantly faster compared to a full N-body simulation combined with a halo finder. The fast halo and subhalo clustering simulations offer good prospects for galaxy forward models using subhalo abundance matching.

The Radcliffe Wave (RW) is a recently discovered sinusoidal vertical feature of dense gas in the proximity of the Sun. In the disk plane, it is aligned with the Local Arm. However, the origin of its vertical undulation is still unknown. This study constrains the kinematics of the RW, using young stars and open clusters as tracers, and explores the possibility of this oscillation being part of a more extended vertical mode. We study the median vertical velocity trends of the young stars and clusters along with the RW and extend it further to the region beyond it. We discover a kinematic wave in the Galaxy, distinct from the warp, with the amplitude of oscillation depending on the age of the stellar population. We perform a similar analysis in the N-body simulation of a satellite as massive as the Sagittarius dwarf galaxy impacting the galactic disk. When projected in the plane, the spiral density wave induced by the satellite impact is aligned with the RW, suggesting that both may be the response of the disk to an external perturbation. However, the observed kinematic wave is misaligned. It appears as a kinematic wave travelling radially, winding up faster than the density wave matched by the RW, questioning its origin. If a satellite galaxy is responsible for this kinematic wave, we predict the existence of a vertical velocity dipole that should form across the disk and this may be measurable with the upcoming Gaia DR3 and DR4.

$\texttt{PyCosmo}$ is a Python-based framework for the fast computation of cosmological model predictions. One of its core features is the symbolic representation of the Einstein-Boltzmann system of equations. Efficient $\texttt{C/C++}$ code is generated from the $\texttt{SymPy}$ symbolic expressions making use of the $\texttt{sympy2c}$ package. This enables easy extensions of the equation system for the implementation of new cosmological models. We illustrate this with three extensions of the $\texttt{PyCosmo}$ Boltzmann solver to include a dark energy component with a constant equation of state, massive neutrinos and a radiation streaming approximation. We describe the $\texttt{PyCosmo}$ framework, highlighting new features, and the symbolic implementation of the new models. We compare the $\texttt{PyCosmo}$ predictions for the $\Lambda$CDM model extensions with $\texttt{CLASS}$, both in terms of accuracy and computational speed. We find a good agreement, to better than 0.1% when using high-precision settings and a comparable computational speed. Links to the Python Package Index (PyPI) page of the code release and to the PyCosmo Hub, an online platform where the package is installed, are available at: https://cosmology.ethz.ch/research/software-lab/PyCosmo.html.

We employ high-resolution zoom-in cosmological simulations to analyze the emerging morphology of galaxies in dark matter halos at redshifts z > 2. We choose DM halos of similar masses of log (Mvir/Mo) ~11.65 +- 0.05 at the target redshifts of z_f = 6, 4 and 2. The rationale for this choice, among others, allows us to analyze how the different growth rate in these halos propagates down to galaxy scales. Halos were embedded in high or low overdensity regions, and two different versions of a galactic wind feedback have been employed. Our main results are: (1) Although our galaxies evolve in different epochs, their global parameters remain within a narrow range. Their morphology, kinematics and stellar populations differ substantially, yet all of them host sub-kpc stellar bars; (2) The SFRs appear higher for larger z_f, in tandem with their energy and momentum feedback; (3) The stellar kinematics allowed separation of bulge from the stellar spheroid. The existence of disk-like bulges has been revealed based on stellar surface density and photometry, but displayed a mixed disk-like and classical bulges based on their kinematics. The bulge-to-total mass ratios appear independent of the last merger time for all z_f. The stellar spheroid-to-total mass ratios of these galaxies lie in the range of ~0.5-0.8; (4) The synthetic redshifted, pixelized and PSF-degraded JWST images allow to detect stellar disks at all z_f. Some bars disappear in degraded images, but others remain visible; (5) Based on the kinematic decomposition, for stellar disks separated from bulges and spheroids. we observe that rotational support in disks depends on the feedback type, but increases with decreasing z_f; (6) Finally, the ALMA images detect disks at all z_f, but their spiral structure is only detectable in z_f=2 galaxies.

Context. The Perseus complex in the outer disk of the Galaxy hosts a number of clusters and associations of young stars. Gaia is providing a detailed characterization of their kinematic structure and evolutionary properties. Aims. Within the SPA Large Programme at the TNG, we secured HARPS-N and GIANO-B high-resolution optical and near-infrared (NIR) spectra of the young red supergiant (RSG) stars in the Perseus complex, in order to obtain accurate radial velocities, stellar parameters and detailed chemical abundances. Methods. We used spectral synthesis to best-fit hundreds of atomic and molecular lines in the spectra of the observed 27 RSGs. We obtained accurate estimates of the stellar temperature, gravity, micro and macro turbulence velocities and chemical abundances for 25 different elements. We also measured the $^{12}$C/$^{13}$C abundance ratio. Results. Our combined optical and NIR chemical study provides homogeneous half-solar iron with a small dispersion, about solar-scaled abundance ratios for the iron-peak, alpha and other light elements and a small enhancement of Na, K and neutron-capture elements, consistent with the thin disk chemistry traced by older stellar populations at a similar Galactocentric distance of about 10 kpc. We inferred enhancement of N, depletion of C and of the $^{12}$C/$^{13}$C isotopic abundance ratio, consistent with mixing processes in the stellar interiors during the RSG evolution.

In simulation-based models of the galaxy-halo connection, theoretical predictions for galaxy clustering and lensing are typically made based on Monte Carlo realizations of a mock universe. In this paper, we use Subhalo Abundance Matching (SHAM) as a toy model to introduce an alternative to stochastic predictions based on mock population, demonstrating how to make simulation-based predictions for clustering and lensing that are both exact and differentiable with respect to the parameters of the model. Conventional implementations of SHAM are based on iterative algorithms such as Richardson-Lucy deconvolution; here we use the JAX library for automatic differentiation to train SHAMNet, a neural network that accurately approximates the stellar-to-halo mass relation (SMHM) defined by abundance matching. In our approach to making differentiable predictions of large-scale structure, we map parameterized PDFs onto each simulated halo, and calculate gradients of summary statistics of the galaxy distribution by using autodiff to propagate the gradients of the SMHM through the statistical estimators used to measure one- and two-point functions. Our techniques are quite general, and we conclude with an overview of how they can be applied in tandem with more complex, higher-dimensional models, creating the capability to make differentiable predictions for the multi-wavelength universe of galaxies.

We present a novel machine learning method for predicting the baryonic properties of dark matter only subhalos from N-body simulations. Our model is built using the extremely randomized tree (ERT) algorithm and takes subhalo properties over a wide range of redshifts as its input features. We train our model using the IllustrisTNG simulations to predict blackhole mass, gas mass, magnitudes, star formation rate, stellar mass, and metallicity. We compare the results of our method with a baseline model from previous works, and against a model that only considers the mass history of the subhalo. We find that our new model significantly outperforms both of the other models. We then investigate the predictive power of each input by looking at feature importance scores from the ERT algorithm. We produce feature importance plots for each baryonic property, and find that they differ significantly. We identify low redshifts as being most important for predicting star formation rate and gas mass, with high redshifts being most important for predicting stellar mass and metallicity, and consider what this implies for nature vs nurture. We find that the physical properties of galaxies investigated in this study are all driven by nurture and not nature. The only property showing a somewhat stronger impact of nature is the present-day star formation rate of galaxies. Finally we verify that the feature importance plots are discovering physical patterns, and that the trends shown are not an artefact of the ERT algorithm.

We present measurements of the expansion of Kepler's Supernova Remnant (SNR) over three epochs of Chandra X-ray observations from 2000, 2006, and 2014. As the remnant of a historical supernova (observed in 1604 CE), Kepler's SNR presents the rare opportunity to study the dynamical evolution of such an object in real time. Measurements of the asymmetry in forward shock velocity can also provide insight into the nature of the explosion and density of the circumstellar material. Combining data from 2014 with previous epochs in 2000 and 2006, we can observe the proper motion of filaments along the outer rim of the SNR. Prior studies of Kepler's SNR have shown proper motion differences up to a factor of 3 between northern and southern regions around the remnant. With the longer time baseline we use here, we find results that are consistent with previous studies, but with smaller uncertainties. Additionally, by adding a third epoch of observations, we search for any systemic change in the velocity in the form of a deceleration of the blast wave, as was recently reported in Tycho's SNR. We find little to no conclusive evidence of such deceleration, and conclude that Kepler's SNR is encountering circumstellar material that is roughly constant in density, though substantially varied around the periphery.

We report 350 pulsating variable stars found in four DECam fields ($\sim 12$ sq. deg.) covering the Antlia 2 satellite galaxy. The sample of variables includes 318 RR Lyrae stars and eight anomalous Cepheids in the galaxy. Reclassification of several objects designated previously to be RR Lyrae as Anomalous Cepheids gets rid of the satellite's stars intervening along the line of sight. This in turn removes the need for prolific tidal disruption of the dwarf, in agreement with the recently updated proper motion and peri-centre measurements based on Gaia EDR3. There are also several bright foreground RR Lyrae stars in the field, and two distant background variables located $\sim 45$ kpc behind Antlia 2. We found RR Lyrae stars over the full search area, suggesting that the galaxy is very large and likely extends beyond our observed area. The mean period of the RRab in Antlia 2 is 0.599 days, while the RRc have a mean period of 0.368 days, indicating the galaxy is an Oosterhoff-intermediate system. The distance to Antlia 2 based on the RR Lyrae stars is $124.1$ kpc ($\mu_0=20.47$) with a dispersion of $5.4$ kpc. We measured a clear distance gradient along the semi-major axis of the galaxy, with the South-East side of Antlia 2 being $\sim13$ kpc farther away from the North-West side. This elongation along the line of sight is likely due to the ongoing tidal disruption of Ant 2.

We investigate the characteristic modifications in the evolving cosmological perturbations when dark energy interacts with dust-like matter, causing the latter's background energy density fall off with time faster than usual. Focusing in particular to the late-time cosmic evolution, we show that such an interaction (of a specific form, arising naturally in a scalar-tensor formulation, or a wide range of modified gravity equivalents thereof), can have a rather significant effect on the perturbative spectrum, than on the background configuration which is not expected to get distorted much from $\L$CDM. Specifically, the matter density contrast, which is by and large scale-invariant in the deep sub-horizon limit, not only gets dragged as the interaction affects the background Hubble expansion rate, but also receives a contribution from the perturbation in the (scalar field induced) dark energy, which oscillates about a non-zero mean value. As such, the standard parametrization ansatz for the the matter density growth factor becomes inadequate. So we modify it suitably, and also find a numerical fit of the growth index in terms of the background parameters, in order to alleviate the problems that arise otherwise. Such a fit enables direct estimations of the background parameters, as well as the growth parameter and the reduced Hubble parameter, which we duly carry out using a redshift space distortion (RSD) subsample and its combination with the observational Hubble data. On the whole, the parametric estimates show consistency with the general observational constraints on the background level cosmology, as well as the constraints on scalar-tensor gravity from astrophysical observations, apart from having significance in the domain of cosmological perturbations.

We study the structural and environmental dependence of the star formation on the plane of stellar mass versus central core density ($\Sigma_{\rm 1\ kpc}$) in the nearby universe. We study the central galaxies in the sparse environment and find a characteristic population-averaged $\rm \Sigma_{1\ kpc} \sim 10^9-10^{9.2}\ M_{\odot}\ kpc^{-2}$, above which quenching is operating. This $\rm \Sigma^{crit}_{1\ kpc}$ only weakly depends on the stellar mass, suggesting that the mass-quenching of the central galaxies is closely related to the processes that operate in the central region of galaxies. For satellites, at a given stellar mass, environment-quenching appears to operate in a similar fashion as mass-quenching in centrals, also starting from galaxies with high $\rm \Sigma_{1\ kpc}$ to low $\rm \Sigma_{1\ kpc}$, and $\rm \Sigma^{crit}_{1\ kpc}$ becomes strongly mass-dependent, in particular in dense regions. This is because (1) more low-mass satellites are quenched by the environmental effects in denser regions and (2) at fixed stellar mass and environment, the environment-quenched satellites on average have, on average, larger $\Sigma_{\rm 1\ kpc}$, $\rm M_{1\ kpc}/M_{\star}$, Sersic index $n$ and as well as smaller size. These results imply that either some dynamical processes change the structure of the satellites during quenching or the satellites with higher $\Sigma_{\rm 1\ kpc}$ are more susceptible to the environmental effects.

Using multi-band data we examine the star formation activity of the nearby group-dominant early-type galaxies of the Complete Local-volume Groups Sample (CLoGS), and the relation between star formation, gas content, and local environment. Only a small fraction of the galaxies (13%; 6/47) are found to be Far-Ultraviolet (FUV) bright, with FUV to near-infrared colours indicative of recent active star formation (NGC 252, NGC 924, NGC 940, NGC 1106, NGC 7252, and ESO 507-25). These systems are lenticulars presenting the highest FUV specific star-formation rates in the sample (sSFR FUV > 5$\times$10$^{13}$ yr$^{-1}$), significant cold gas reservoirs (M(H2)=0.5 - 61$\times$10$^8$ M$_\odot$), reside in X-ray faint groups, and none of them hosts a powerful radio AGN (P$_{1.4GHz}$<$10^{23}$ W Hz$^{-1}$). The majority of the group-dominant galaxies (87%; 41/47) are FUV faint, with no significant star formation, classified in most cases as spheroids based on their position on the infrared star-forming main sequence (87%; 46/53). Examining the relationships between radio power, SFR FUV, and stellar mass we find a lack of correlation that suggests a combination of origins for the cool gas in these galaxies, including stellar mass loss, cooling from the intra-group medium (IGrM) or galaxy halo, and acquisition through mergers or tidal interactions. X-ray bright systems, in addition to hosting radio powerful AGN, have a range of SFRs but, with the exception of NGC 315, do not rise to the highest rates seen in the FUV bright systems. We suggest that central group galaxy evolution is linked to gas mass availability, with star formation favoured in the absence of a group-scale X-ray halo, but AGN jet launching is more likely in systems with a cooling IGrM.

The Extreme Universe Space Observatory on a Super Pressure Balloon II (EUSO-SPB2) is a second generation stratospheric balloon instrument for the detection of Ultra High Energy Cosmic Rays (UHECRs, E > 1 EeV) via the fluorescence technique and of Very High Energy (VHE, E > 10 PeV) neutrinos via Cherenkov emission. EUSO-SPB2 is a pathfinder mission for instruments like the proposed Probe Of Extreme Multi-Messenger Astrophysics (POEMMA). The purpose of such a space-based observatory is to measure UHECRs and UHE neutrinos with high statistics and uniform exposure. EUSO-SPB2 is designed with two Schmidt telescopes, each optimized for their respective observational goals. The Fluorescence Telescope looks at the nadir to measure the fluorescence emission from UHECR-induced extensive air shower (EAS), while the Cherenkov Telescope is optimized for fast signals ($\sim$10 ns) and points near the Earth's limb. This allows for the measurement of Cherenkov light from EAS caused by Earth skimming VHE neutrinos if pointed slightly below the limb or from UHECRs if observing slightly above. The expected launch date of EUSO-SPB2 is Spring 2023 from Wanaka, NZ with target duration of up to 100 days. Such a flight would provide thousands of VHECR Cherenkov signals in addition to tens of UHECR fluorescence tracks. Neither of these kinds of events have been observed from either orbital or suborbital altitudes before, making EUSO-SPB2 crucial to move forward towards a space-based instrument. It will also enhance the understanding of potential background signals for both detection techniques. This contribution will provide a short overview of the detector and the current status of the mission as well as its scientific goals.

3A 1954+319 has been classified for a long time as a symbiotic X-ray binary, hosting a slowly rotating neutron star and an aged M red giant. Recently, this classification has been revised thanks to the discovery that the donor star is an M supergiant. This makes 3A 1954+319 a rare type of high mass X-ray binary consisting of a neutron star and a red supergiant donor. In this paper, we analyse two archival and still unpublished XMM-Newton and NuSTAR observations of the source. We perform a detailed hardness ratio-resolved spectral analysis to search for spectral variability that could help investigating the structures of the inhomogeneous M supergiant wind from which the neutron star is accreting. We discuss our results in the context of wind-fed supergiant X-ray binaries and show that the newest findings on 3A 1954+319 reinforce the hypothesis that the neutron star in this system is endowed with a magnetar-like magnetic field strength ($\gtrsim10^{14}$ G).

Titan's thick atmosphere is primarily composed of nitrogen and methane. Complex chemistry happening in Titan's atmosphere produces optically thick organic hazes. These hazes play significant roles in Titan's atmosphere and on its surface, and their optical properties are crucial for understanding many processes happening on Titan. Due to the lack of such information, the optical constants of laboratory prepared Titan haze analogues are essential inputs for atmospheric modeling and data analysis of remote sensing observations of Titan. Here, we perform laboratory simulations in a Titan relevant environment, analyze the resulting Titan haze analogues using vacuum Fourier transform infrared spectroscopy, and calculate the optical constants from the measured transmittance and reflectance spectra. We provide a reliable set of optical constants of Titan haze analogue in the wavelength range from 0.4 to 3.5 micron and will extend to 28.5 micron in the near future, which can both be used for analyzing existing and future observational data of Titan. This study establishes a feasible method to determine optical constants of haze analogues of (exo)planetary bodies.

We analyzed TESS data of the IW And-type dwarf nova ST Cha during ordinary dwarf nova states. We have identified an orbital period of 0.285360(1) d. The object was reported to be eclipsing in the past using the data obtained in the 1970s, which is not in agreement with the present data. Despite the constant mean brightness, the strength of the orbital signal varied significantly, suggesting that the strength of the orbital signal does not always reflect the mass-transfer rate. During an outburst with a shoulder, we did not find evidence of humps recurring with a period longer than the orbital period which were recorded in V363 Lyr. This finding strengthens the idea that V363 Lyr is an unusual object. We found that the strength of the orbital signal increased after an outburst with a shoulder. This outburst may have changed the state of the disk and the hot spot became more apparent. Such a change in the disk may have triggered a transition from an ordinary dwarf nova-type state to an IW And-type state and this possibility would require further examination.

We examine massive black hole (MBH) mergers and their associated gravitational wave signals from the large-volume cosmological simulation Astrid. Astrid includes galaxy formation and black hole models recently updated with a MBH seed population between $3\times 10^4M_{\odot}/h$ and $3\times 10^5M_{\odot}/h$ and a sub-grid dynamical friction (DF) model to follow the MBH dynamics down to $1.5\;\text{ckpc}/h$. We calculate initial eccentricities of MBH orbits directly from the simulation at kpc-scales, and find orbital eccentricities above $0.7$ for most MBH pairs before the numerical merger. After approximating unresolved evolution on scales below ${\sim 200\,\text{pc}}$, we find that the in-simulation DF on large scales accounts for more than half of the total orbital decay time ($\sim 500\,\text{Myrs}$) due to DF. The binary hardening time is an order of magnitude longer than the DF time, especially for the seed-mass binaries ($M_\text{BH}<2M_\text{seed}$). As a result, only $\lesssim20\%$ of seed MBH pairs merge at $z>3$ after considering both unresolved DF evolution and binary hardening. These $z>3$ seed-mass mergers are hosted in a biased population of galaxies with the highest stellar masses of $>10^9\,M_\odot$. With the higher initial eccentricity prediction from Astrid, we estimate an expected merger rate of $0.3-0.7$ per year from the $z>3$ MBH population. This is a factor of $\sim 7$ higher than the prediction using the circular orbit assumption. The LISA events are expected at a similar rate, and comprise $\gtrsim 60\%$ seed-seed mergers, $\sim 30\%$ involving only one seed-mass MBH, and $\sim 10\%$ mergers of non-seed MBHs.

Modern cosmological data demand modern data analysis techniques. We introduce BayOp, a new likelihood sampling and maximisation method which is based on the Bayesian Optimisation algorithm and learns a function instead of randomly sampling from it. We apply BayOp to analyse Planck data for traces of inflationary features models with global periodic modulations of the primordial power spectrum. While we do not find any new evidence for features, we demonstrate that BayOp provides an extremely efficient way of sampling likelihoods over low-to-moderate-dimensional parameter spaces, even for very complex likelihood landscapes.

Building upon three late-type galaxies in the Virgo cluster with both a predicted black hole mass of less than $\sim$10$^5$ M$_{\odot}$ and a centrally-located X-ray point-source, we reveal 11 more such galaxies, more than tripling the number of active intermediate-mass black hole candidates among this population. Moreover, this amounts to a 36$\pm$8% X-ray detection rate (despite the sometimes high, X-ray-absorbing, HI column densities), compared to just 10$\pm$5% for (the largely HI-free) dwarf early-type galaxies in the Virgo cluster. The expected contribution of X-ray binaries from the galaxies' inner field stars is negligible. Moreover, given that both the spiral and dwarf galaxies contain nuclear star clusters, the above inequality appears to disfavor X-ray binaries in nuclear star clusters. The higher occupation, or rather detection, fraction among the spiral galaxies may instead reflect an enhanced cool gas/fuel supply and Eddington ratio. Indeed, four of the 11 new X-ray detections are associated with known LINERs or LINER/HII composites. For all (four) of the new detections for which the X-ray flux was strong enough to establish the spectral energy distribution in the Chandra band, it is consistent with power-law spectra. Furthermore, the X-ray emission from the source with the highest flux (NGC 4197: $L_X \approx 10^{40}$ erg s$^{-1}$) suggests a non-stellar-mass black hole if the X-ray spectrum corresponds to the `low/hard state'. Follow-up observations to further probe the black hole masses, and prospects for spatially resolving the gravitational spheres-of-influence around intermediate-mass black holes, are reviewed in some detail.

The study of secondary particles produced by the cosmic-ray interaction in the Earth's atmosphere is very crucial as these particles mainly constitute the background counts produced in the high-energy detectors at balloon and satellite altitudes. In the present work, we calculate the abundance of cosmic-ray generated secondary particles at various heights of the atmosphere by means of a Monte Carlo simulation and use this result to understand the background counts in our X-ray observations using balloon-borne instruments operating near the tropical latitude (geomagnetic latitude: $\sim 14.50^{\circ}$ N). For this purpose, we consider a 3D description of the atmospheric and geomagnetic field configurations surrounding the Earth, as well as the electromagnetic and nuclear interaction processes using Geant4 simulation toolkit. Subsequently, we use a realistic mass model description of the detector under consideration, to simulate the counts produced in the detector due to secondary cosmic-ray particles.

The Cosmic Ray Albedo Neutron Decay (CRAND) is believed to be the principal mechanism for the formation of inner proton radiation belt -- at least for relatively higher energy particles. We implement this mechanism in a Monte Carlo simulation procedure to calculate the trapped proton radiation at the low Earth orbits, through event-by-event interaction of the cosmic ray particles in the Earth's atmosphere and their transportation in the magnetosphere. We consider the generation of protons from subsequent decay of the secondary neutrons from the cosmic ray interaction in the atmosphere and their transport (and/or trapping) in the geomagnetic field. We address the computational challenges for this type of calculations and develop an optimized algorithm to minimize the computation time. We consider a full 3D description of the Earth's atmospheric and magnetic-field configurations using the latest available models. We present the spatial and phase-space distribution of the trapped protons considering the adiabatic invariants and other parameters at the low Earth orbits. We compare the simulation results with the trapped proton flux measurements made by PAMELA experiment at low Earth orbit and explain certain features observed by the measurement.

The accuracy of galaxy photometric redshift (photo-z) can significantly affect the analysis of weak gravitational lensing measurements, especially for future high-precision surveys. In this work, we try to extract photo-z information from both galaxy flux and image data obtained by China Space Station Telescope (CSST) using neural networks. We generate mock galaxy images based on the observational images from the Advanced Camera for Surveys of Hubble Space Telescope (HST-ACS) and COSMOS catalogs, considering the CSST instrumental effects. Galaxy flux data are then measured directly from these images by aperture photometry. The Multi-Layer Perceptrons (MLP) and Convolutional Neural Network (CNN) are constructed to predict photo-$z$ from fluxes and images, respectively. We also propose to use an efficient hybrid network, which combines MLP and CNN, by employing transfer learning techniques to investigate the improvement of the result with both flux and image data included. We find that the photo-z accuracy and outlier fraction can achieve sigma_NMAD = 0.023 and eta = 1.43 for the MLP using flux data only, and sigma_NMAD = 0.025 and eta = 1.21% for the CNN using image data only. The result can be further improved in high efficiency as sigma_NMAD = 0.020 and eta = 0.90% for the hybrid transfer network. This indicates that our networks can effectively and properly extract photo-z information from the galaxy flux and image data, which can achieve the scientific requirement of the CSST photometric surveys.

AR Sco is a binary system that contains a white and red dwarf. The rotation rate of the white dwarf has been observed to slow down, analogous to rotation-powered radio pulsars; it has thus been dubbed a "white dwarf pulsar". We previously fit the traditional radio pulsar rotating vector model to the linearly polarised optical data from this source, constraining the system geometry as well as the white dwarf mass. Using a much more extensive dataset, we now explore the application of the same model to binary phase-resolved optical polarimetric data, thought to be the result of non-thermal synchrotron radiation, and derive the magnetic inclination angle $\alpha$ and the observer angle $\zeta$ at different orbital phases. We obtain a $\sim 10^{\circ}$ variation in $\alpha$ and $\sim 30^{\circ}$ variation in $\zeta$ over the orbital period. The variation patterns in these two parameters is robust, regardless of the binning and epoch of data used. We speculate that the observer is detecting radiation from an asymmetric emission region that is a stable structure over several orbital periods. The success of this simple model lastly implies that the pitch angles of the particles are small and the pulsed, non-thermal emission originates relatively close to the white dwarf surface.

In this paper we report the first identification of jellyfish galaxies in the Perseus cluster (Abell 426). We identified four jellyfish galaxies (LEDA 2191078, MCG +07-07-070, UGC 2654, UGC 2665) within the central $2^\circ \times 2^\circ$ ($2.6\,\mathrm{Mpc} \times 2.6\,\mathrm{Mpc}$) of Perseus based on the presence of one-sided radio continuum tails that were detected at $144\,\mathrm{MHz}$ by the LOw Frequency ARray (LOFAR). The observed radio tails, as well as the orientation of morphological features in the rest-frame optical, are consistent with these four galaxies being impacted by ram pressure stripping as they orbit through the Perseus intracluster medium. By combining the LOFAR imaging at 144 MHz with 344 MHz imaging from the Karl G. Jansky Very Large Array, we derived spectral indices for the disks and the stripped tails of these jellyfish galaxies. We show that the spectral indices over the galaxy disks are quite flat, while the indices of the stripped tails are substantially steeper. We also identified a number of compact $\mathrm{H\alpha + [NII]}$ sources with narrowband imaging from the Isaac Newton Telescope. These sources are brighter along the leading side of the galaxy (i.e., opposite to the direction of the stripped tail), which is consistent with ram pressure induced star formation. Lastly, consistent with previous works in other clusters, we find that these jellyfish galaxies show enhanced radio luminosities for their observed star formation rates. Given the small distance to the Perseus cluster ($D \sim 70\,\mathrm{Mpc}$, $1'' \simeq 340\,\mathrm{pc}$), these galaxies are excellent candidates for multiwavelength follow-up observations to probe the impact of ram pressure stripping on galaxy star formation at subkiloparsec scales.

HB9 (G160.9+2.6) is a mixed-morphology Galactic supernova remnant (SNR) at a distance of $\sim$0.6 kpc. Previous analyses revealed recombining plasma emission in X-rays and an expanding shell structure in HI and CO emission, which were correlating with the spatial extent of HB9. In GeV energies, HB9 was found to show extended gamma-ray emission with a morphology that is consistent with the radio continuum emission showing a log-parabola-type spectrum. The overlap reported between the gas data and the excess gamma-ray emission at the southern region of the SNR's shell could indicate a possible interaction between them. We searched for hadronic gamma-ray emission signature in the spectrum to uncover possible interaction between the molecular environment and the SNR. Here we report the results of the gamma-ray spectral modelling studies of HB9.

The discovery of the sources of ultra-high energy photons by High-Altitude Water Cerenkov Gamma ray Observatory and Large High Altitude Air Shower Observatory in our Galaxy has revolutionised the field of gamma ray astronomy in the last few years. These emissions are sometimes found in the vicinity of powerful pulsars or supernova remnants associated with giant molecular clouds. Inverse Compton emission by shock accelerated electrons emitted by pulsars and proton-proton interactions of shock accelerated protons emitted by supernova remnants with cold protons in molecular clouds are often identified as the causes of these emissions. In this paper we have selected two ultra-high energy photon sources LHAASO J2108+5157 and LHAASO J0341+5258 which are associated with giant molecular clouds, but no powerful pulsar or supernova remnant has been detected in their vicinity.We have proposed a scenario where shock accelerated electrons and protons are injected in the local environment of these sources from past explosions, which happened thousands of years ago. We show that the observed ultra-high energy photon flux can be explained with the secondary gamma rays produced by the time evolved relativistic electron and proton spectra.

Mini-EUSO is the first detector of the JEM-EUSO program deployed on the ISS. It is a wide field of view telescope currently operating from a nadir-facing UV-transparent window on the ISS. It is based on an array of MAPMTs working in photon counting mode with a 2.5 $\mu$s time resolution. Among the different scientific objectives it searches for light signals with time duration compatible to those expected from Extensive Air Showers (EAS) generated by EECRs interacting in the atmosphere. Although the energy threshold for cosmic ray showers is above $E>10^{21}$ eV, due the constraints given by the size of the UV-transparent window, the dedicated trigger logic has been capable of the detection of other interesting classes of events, like elves and ground flashers. An overview of the general performance of the trigger system is provided, with a particular focus on the identification of classes of events responsible for the triggers.

The physical nature of the mechanism responsible for the emission of neutrinos in active galactic nuclei (AGN) has been matter of debate in the literature, with relativistic jets of radio-loud AGNs as possible candidates to be the sources of high energy neutrinos. The most prominent candidate so far is the blazar TXS 0506+056, which is found to be associated with the neutrino event IceCube-170922A. Furthermore, the IceCube reported an excess of neutrinos towards TXS 0506+056 between September 2014 and March 2015, even though this association needs additional investigation, considering the presence of a nearby gamma-ray source, the quasar PKS 0502+049. Motivated by this, we studied the parsec-scale structures of TXS 0506+056 and PKS 0502+049 through radio interferometry at 8 and 15\,GHz. We identified twelve jet components in TXS 0506+056 and seven components in PKS 0502+049. The most reliable jet components show superluminal speeds ranging from 9.5c to 66c in the case of TXS 0506+056, and from 14.3c to 59c for PKS 0502+049, which were used to estimate a lower (upper) limit for the Lorentz factor (jet viewing angle) for both sources. A novel approach using simultaneously the brightness temperature of the core region and the apparent speeds of the jet components allowed us to infer basic jet parameters for TXS 0506+056 at distinct epochs. We also found that the emergence of new jet components coincides with the occurrence of gamma-ray flares. Interestingly, two of these coincidences in the case of PKS 0502+049 and one for TXS 0506+056 seems to be correlated with neutrino events detected by the IceCube Observatory.

The optical radiation emitted by blazars contains contributions from synchrotron radiation by relativistic electrons in the jets, as well as thermal radiation emitted mainly by the Accretion Disk (AD), the Broad Line Region (BLR) and the host galaxy. The unpolarized radiation components from the AD, BLR and host galaxy present themselves by decreasing the total polarization in the optical/ultraviolet(UV) spectrum. A combined model for the Spectral Energy Distribution (SED) and degree of optical/UV polarization is constructed, enabling the disentanglement of the synchrotron and AD components. Our model is applied to the multi-wavelength SED and spectropolarimetry observations of the Flat Spectrum Radio Quasar 4C+01.02 ($z = 2.1$) in its 2016 July-August flaring state and July-August 2017 quiescent state, using data from the Fermi Large Area Telescope, the Southern African Large Telescope and the Las Cumbres Observatory network of telescopes. By constraining the AD component, the mass of the super massive black hole is obtained as $3 \times 10^9 \rm M_{\odot}$. Furthermore, the model retrieves the characteristics of the relativistic electron distribution in the jet and the degree of ordering of the magnetic field. Our results highlight the potential of spectropolarimetry observations for disentangling thermal from non-thermal (jet) emission components and thus revealing the physics of particle acceleration and high-energy emission in active galactic nuclei jets.

The accretion of planetary debris into the atmospheres of white dwarfs leads to the presence of metal lines in their spectra. Cool metal-rich white dwarfs, which left the main-sequence many Gyr ago, allow the study of the remnants of the oldest planetary systems. Despite their low effective temperatures ($T_\mathrm{eff}$), a non-neglible amount of their flux is emitted in the near ultraviolet (NUV), where many overlapping metal lines can potentially be detected. We have observed three metal-rich cool white dwarfs with the Space Telescope Imaging Spectrograph (STIS) onboard the Hubble Space Telescope (HST), and compare the results determined from the NUV data with those previously derived from the analysis of optical spectroscopy. For two of the white dwarfs, SDSSJ1038-0036 and SDSSJ1535+1247, we find reasonable agreement with our previous analysis and the new combined fit of optical and NUV data. For the third object, SDSSJ0956+5912, including the STIS data leads to a ten percent lower $T_\mathrm{eff}$, though we do not identify a convincing explanation for this discrepancy. The unusual abundances found for SDSSJ0956+5912 suggest that the accreted parent-body was composed largely of water ice and magnesium silicates, and with a mass of up to $\simeq 2\times 10^{25}$g. Furthermore SDSSJ0956+5912 shows likely traces of atomic carbon in the NUV. While molecular carbon is not observed in the optical, we demonstrate that the large quantity of metals accreted by SDSSJ0956+5912 can suppress the C$_2$ molecular bands, indicating that planetary accretion can convert DQ stars into DZs (and not DQZs/DZQs).

The first set of theoretical cross sections for propylene oxide (CH3CHCH2O) colliding with cold He atoms has been obtained at the full quantum level using a high-accuracy potential energy surface. By scaling the collision reduced mass, rotational rate coefficients for collisions with para-H2 are deduced in the temperature range 5-30 K. These collisional coefficients are combined with radiative data in a non-LTE radiative transfer model in order to reproduce observations of propylene oxide made towards the Sagittarius B2(N) molecular cloud with the Green Bank and Parkes radio telescopes. The three detected absorption lines are found to probe the cold (~ 10 K) and translucent (nH ~ 2000 cm-3) gas in the outer edges of the extended Sgr B2(N) envelope. The derived column density for propylene oxide is Ntot ~ 3e12 cm-2, corresponding to a fractional abundance relative to total hydrogen of ~ 2.5e-11. The present results are expected to help our understanding of the chemistry of propylene oxide, including a potential enantiomeric excess, in the cold interstellar medium.

We present periodicity search analyses on the $\gamma$-ray lightcurve of the TeV blazar PKS 1510-089 observed by the \textit{Fermi} Large Area Telescope. We report the detection of two transient quasi-periodic oscillations: a 3.6-day QPO during the outburst in 2009 that lasted five cycles (MJD 54906--54923); and a periodicity of 92 days spanning over 650 days from 2018 to 2020 (MJD 58200--58850), which lasted for seven cycles. We employed the Lomb-Scargle periodogram, Weighted Wavelet Z-transform, REDFIT, and the Monte Carlo lightcurve simulation techniques to find any periodicity and the corresponding significance. The 3.6-day QPO was detected at a moderate significance of $\sim$3.5$\sigma$, while the detection significance of the 92-day QPO was $\sim$7.0$\sigma$. We explore a few physical models for such transient QPOs including a binary black hole system, precession of the jet, a non-axisymmetric instability rotating around the central black hole near the innermost stable circular orbit, the presence of quasi-equidistant magnetic islands inside the jet, and a geometric model involving a plasma blob moving helically inside a curved jet.

Since the first discovery of an extra-solar planet around a main-sequence star, in 1995, the number of detected exoplanets has increased enormously. Over the past two decades, observational instruments (both onboard and on ground-based facilities) have revealed an astonishing diversity in planetary physical features (i. e. mass and radius), and orbital parameters (e.g. period, semi-major axis, inclination). Exoplanetary atmospheres provide direct clues to understand the origin of these differences through their observable spectral imprints. In the near future, upcoming ground and space-based telescopes will shift the focus of exoplanetary science from an era of 'species discovery' to one of 'atmospheric characterization'. In this context, the Atmospheric Remote-sensing Infrared Exoplanet Large (Ariel) survey, will play a key role. As it is designed to observe and characterize a large and diverse sample of exoplanets, Ariel will provide constraints on a wide gamut of atmospheric properties allowing us to extract much more information than has been possible so far (e.g. insights into the planetary formation and evolution processes). The low resolution spectra obtained with Ariel will probe layers different from those observed by ground-based high resolution spectroscopy, therefore the synergy between these two techniques offers a unique opportunity to understanding the physics of planetary atmospheres. In this paper, we set the basis for building up a framework to effectively utilise, at near-infrared wavelengths, high-resolution datasets (analyzed via the cross-correlation technique) with spectral retrieval analyses based on Ariel low-resolution spectroscopy. We show preliminary results, using a benchmark object, namely HD 209458 b, addressing the possibility of providing improved constraints on the temperature structure and molecular/atomic abundances.

Results. Gas phase silicon was present throughout the Rosetta mission. Furthermore, the presence of sodium and iron atoms near the comet's perihelion confirms that sputtering cannot be the sole release process for refractory elements into the gas phase. Nickel was below detection limit. The search for parent species of any of the identified gas phase refractories has not been successful. Upper limits for a suite of possible fragment species (SiH, SiC, NaH, ...) of larger parent and daughter species have been obtained. Furthermore, Si did not exhibit the same drop in signal like common cometary gases when the spacecraft was pointed away from the nucleus. The combined results suggest that direct release of elemental species from small grains on the surface of the nucleus and/or from small grains in the surrounding coma are more likely than release through dissociation of gaseous parent molecules.

A model based on celestial geometry and atmospheric physics predicts the dimming and the color of lunar eclipses. Visual magnitudes and color indices for eclipses from year 2000 through 2050 are listed. The enlargement of the Earth's umbral shadow reported by observers for over 300 years is explained. The geometrical aspects of the model are the sizes and separations of the Sun, Moon and Earth. Atmospheric effects include refraction, absorption and focusing of sunlight.

The Thesis is an attempt to combine data from three previously independent areas - the structure and kinematics of stellar populations of the Galaxy, photometric and dynamical models of galaxies, and models of the dynamical and physical evolution of galaxies. This synthesis was made with the goal to understand better the vast topic of the structure and evolution of galaxies. The main results of the study can be divided into methodical and astronomical. The methodical results are 1) extrapolation of the mass distribution function beyond the Sun distance, and the determination of the circular velocity at the Sun distance from the Galactic centre; 2) conditions of physical correctness of models are developed, and the generalised exponential model is suggested; 3) a method has been developed to construct spatial and hydrodynamical models of stellar systems, using a combination of observational data on populations, and data on physical evolution of populations. The main astronomical results are - 1) spatial and hydrodynamical models of the Galaxy and the Andromeda galaxy M31 are elaborated in several approximations; 2) the dynamical evolution of the Galaxy is reconstructed using kinematical characteristics of stellar populations of different ages; 3) on the basis of stellar evolutionary tracks and star formation function, a theory of the evolution of galaxies is elaborated. The basic conclusion of the study is: it is impossible to reproduce the observed rotation curves of galaxies with known stellar populations.

The gravitational wave/gamma-ray burst GW/GRB170817 event marked the beginning of the era of multi-messenger astrophysics, in which new observations of Gravitational Waves (GW) are combined with traditional electromagnetic observations from the very same astrophysical source. In the next few years, Advanced LIGO/VIRGO and KAGRA in Japan and LIGO-India will reach their nominal/ultimate sensitivity. In the electromagnetic domain, the Vera C. Rubin Observatory and the Cherenkov Telescope Array (CTA) will come online in the next few years, and they will revolutionize the investigation of transient and variable cosmic sources in the optical and TeV bands. The operation of an efficient X-ray/gamma-ray all-sky monitor with good localisation capabilities will play a pivotal role in providing the high-energy counterparts of the GW interferometers and Rubin Observatory, bringing multi-messenger astrophysics to maturity. To reach the required precision in localisation and timeliness for an unpredictable physical event in time and space requires a sensor distribution covering the whole sky. We discuss the potential of large-scale, small-platform-distributed architectures and constellations to build a sensitive X-ray/gamma-ray all-sky monitor and the programmatic implications of this, including the set-up of an efficient assembly line for both hardware development and data analysis. We also discuss the potential of a constellation of small platforms operating at other wavelengths (UV/IR) that are capable of repointing quickly to follow-up high-energy transients.

We present a numerical study, based on Monte Carlo simulations, aimed at defining new empirical parameters measurable from observations and able to trace the different phases of star cluster dynamical evolution. As expected, a central density cusp, deviating from the King model profile, develops during the core collapse (CC) event. Although the slope varies during the post-CC oscillations, the cusp remains a stable feature characterizing the central portion of the density profile in all post-CC stages. We then investigate the normalized cumulative radial distribution (nCRD) drawn by all the cluster stars included within one half the tridimensional half-mass radius (R<0.5 rh), finding that its morphology varies in time according to the cluster's dynamical stage. To quantify these changes we defined three parameters: A5, the area subtended by the nCRD within 5% of the half-mass radius, P5, the value of the nCRD measured at the same distance, and S2.5, the slope of the straight line tangent to the nCRD measured at R=2.5% rh. The three parameters evolve similarly during the cluster's dynamical evolution: after an early phase in which they are essentially constant, their values rapidly increase, reaching their maximum at the CC epoch and slightly decreasing in the post-CC phase, when their average value remains significantly larger than the initial one, in spite of some fluctuations. The results presented in the paper suggest that these three observable parameters are very promising empirical tools to identify the star cluster's dynamical stage from observational data.

We investigate the resolution dependence of HII regions expanding past their Str\"{o}mgren spheres. We find that their structure and size, and the radial momentum that they attain at a given time, is in good agreement with analytical expectations if the Str\"{o}mgren radius is resolved with $dr \leq 0.3\,R_{\rm st}$. If this is not satisfied, the radial momentum may be over- or under-estimated by factors up to 10 or more. Our work has significance for the amount of radial momentum that a HII region can impart to the ambient medium in numerical simulations, and thus on the relative importance of ionizing feedback from massive stars.

Our aim in this note is to compare a recent explanation of the galactic X pattern in Faraday Rotation (XRM), to that produced by the advection part of the classical dynamo.We find that the characteristic X magnetic field polarization in the plane of the sky, found in edge-on spiral galaxies, can develop magnetohydromagnetically from an initial disc magnetic field combined with wind and rotation. The Rotation Measure develops a corresponding X distribution in sign, but this distribution is not a `universal' behaviour because it depends primarily on the velocity field. We use Cauchy evolution of an initial magnetic field to find the field at some later time and place. A `battery' mechanism that requires current to always flow out of a galaxy has been recently suggested, which contrasts with the conclusions of this paper. Either explanation has significant consequences for the structure of a spiral galaxy. If the battery mechanism applies then we have a new method of producing magnetic field independent of the traditional dynamo. If however the flow mechanism applies; then constraints concerning the presence and nature of a galactic wind, together with the signature of the mean radial magnetic field, can be inferred.

We have obtained 1 and 3 mm spectral scans of ASXDF1100.053.1 using NOEMA. ASXDF1100.053.1 is an unlensed optically dark millimetre-bright SMG with $K_{\rm AB}>25.7$ ($2\sigma$), which was expected to lie at $z=$5-7 based on its radio-submm photo-$z$. Our data detected line emission due to $^{12}$CO($J=$5-4) and ($J=$6-5), providing a $z_{\rm CO}= 5.2383\pm0.0005$. Energy-coupled SED modelling indicates properties of $L_{\rm IR}=8.3^{+1.5}_{-1.4}\times10^{12}$ L$_{\odot}$, SFR $=630^{+260}_{-380}$ M$_{\odot}$ yr$^{-1}$, $M_{\rm dust}=4.4^{+0.4}_{-0.3}\times10^{8}$ M$_{\odot}$, $M_{\rm stellar}=3.5^{+3.6}_{-1.4}\times10^{11}$ M$_{\odot}$, and $T_{\rm dust}=37.4^{+2.3}_{-1.8}$ K. The CO luminosity allows us to estimate a gas mass $M_{\rm gas}=3.1\pm0.3\times10^{10}$ M$_{\odot}$, suggesting a gas-to-dust mass ratio of around 70, fairly typical for $z\sim2$ SMGs. ASXDF1100.053.1 has $R_{\rm e, mm}=1.0^{+0.2}_{-0.1}$ kpc, so its surface $L_{\rm IR}$ density $\Sigma_{\rm IR}$ is $1.2^{+0.1}_{-0.2}\times10^{12}$ L$_{\odot}$ kpc$^{-2}$. These properties indicate that ASXDF1100.053.1 is a massive dusty star-forming (SF) galaxy with an unusually compact starburst. It lies close to the SF main sequence at $z\sim5$, with low $M_{\rm gas}$/$M_{\rm stellar}=0.09$, SFR/SFR$_{\rm MS} (R_{\rm SB})=0.6$, and a gas-depletion time $\tau_{\rm dep}$ of $\approx 50$ Myr, modulo assumptions about the stellar initial mass function in such objects. ASXDF1100.053.1 has extreme values of $M_{\rm gas}/M_{\rm stellar}$, $R_{\rm SB}$, and $\tau_{\rm dep}$ compared to SMGs at $z\sim$2-4, and those of ASXDF1100.053.1 are the smallest among SMGs at $z>5$. ASXDF1100.053.1 is likely a late-stage dusty starburst prior to passivisation. The number of $z=$5.1-5.3 unlensed SMGs now suggests a number density $dN/dz=30.4\pm19.0$ deg$^{-2}$, barely consistent with the latest cosmological simulations.

Aims. We use the recently published database (Trifonov et al. 2020) of radial velocities (RVs) that were derived from fifteen years of HARPS/ESO observations to search for planet candidates. Methods. For targets with sufficient RV data, we apply an automated algorithm to identify significant periodic signals and fit a Keplerian model for orbital estimates. We also search the auxiliary data of stellar-activity indices and compare our findings with existing literature, to detect periodic RV signals that have no counterpart in the activity timeseries. The most convincing signals are then manually inspected to designate additional false planet detection, focusing the search on long-period (P > 1 000 d) massive candidates around FGK dwarf stars. Results. We identify two Jupiter analogs, in orbit around the slightly evolved F8V star HD 103891 and the Solar-like star HD 105779. We use nested sampling to derive their orbital parameters, and find their orbital periods to be 1919 +/- 16 d and 2412 +/- 54 d, while their minimum masses are 1.44 +/- 0.02 M Jup and 0.64 +/- 0.06 M Jup , respectively. While the orbit of HD 103891 b is slightly eccentric (e = 0.31 +/- 0.03), that of HD 105779 b is likely circular (e < 0.16). Conclusions. With minimum astrometric signatures of 59 and 42 $\mu$as, HD 103891 b and HD 105779 b join the growing sample of planets whose exact masses may soon be derived with Gaia astrometry. This finding also highlights the importance of long-term RV surveys to study planetary occurrence beyond the snow line of Solar-like stars.

Star-forming galaxies (SFGs) are expected to harbour an abundant reservoir of cosmic rays (CRs). At GeV energies, these CRs can undergo hadronic interactions with interstellar gases to produce $\gamma$-rays, and the unresolved $\gamma$-ray emission from populations of SFGs form a component of the extragalactic $\gamma$-ray background (EGB). In this work, we investigate the contribution to the 0.01 - 50 GeV EGB from SFG populations located up to redshift $z=3$. We consider their redshift evolution and variations in their physical properties, and model how this affects their contribution to the EGB. We find this is dominated by starbursts, while the contribution from main sequence SFGs is marginal at all energies. We also demonstrate that most of the $\gamma$-ray contribution from SFGs emanates from low mass galaxies, with over 80 per cent of the emission originating from galaxies with stellar masses below $10^8 \;\!{\rm M}_{\odot}$. We show that the EGB at different energies captures different stages of the evolution of the source galaxies. At 0.01 GeV, the emission is dominated by galaxies at the noon of cosmic star-formation, around $z\sim 2$, however higher energy $\gamma$-rays are instead mainly contributed from low mass starburst populations at higher redshifts, $\sim$ 700 Myr earlier. The redshift distributions of the EGB sources at different energies imprint intensity signatures at different angular scales, allowing their contribution to be distinguished using analyses of small-scale EGB intensity anisotropies. We show that the EGB is sensitive to the evolution of low mass populations of galaxies, particularly around $z\sim2-2.5$, and that it provides a new means to probe the engagement of CRs in these galaxies before and during the high noon of cosmic star-formation.

Solar radio flux along with geomagnetic indices are important indicators of solar activity and its effects. Extreme solar events such as flares and geomagnetic storms can negatively affect the space environment including satellites in low-Earth orbit. Therefore, forecasting these space weather indices is of great importance in space operations and science. In this study, we propose a model based on long short-term memory neural networks to learn the distribution of time series data with the capability to provide a simultaneous multivariate 27-day forecast of the space weather indices using time series as well as solar image data. We show a 30-40\% improvement of the root mean-square error while including solar image data with time series data compared to using time series data alone. Simple baselines such as a persistence and running average forecasts are also compared with the trained deep neural network models. We also quantify the uncertainty in our prediction using a model ensemble.

High-energy astrophysical neutrinos have been observed by several telescopes in the last decade, but their sources still remained unknown. We address the problem of locating astrophysical neutrinos' sources in a statistical manner. We show that blazars positionally associated with IceCube neutrino detections have stronger parsec-scale radio cores than the rest of the sample. The probability of a chance coincidence is only 4*10^-5 corresponding to a significance of 4.1 sigma. We explicitly list four strong radio blazars as highly probable sources of neutrinos above 200 TeV: 3C 279, NRAO 530, PKS 1741-038, and PKS 2145+067. There are at least 70 more radio-bright blazars that emit neutrinos of lower energies, starting from TeVs. Using continuous RATAN-600 monitoring of VLBI-selected blazars, we find that radio flares at frequencies above 10 GHz coincide with neutrino arrival dates. The most pronounced example of such behavior is PKS 1502+106 that experienced a major flare in 2019. We demonstrate that the majority of IceCube astrophysical neutrino flux derived from muon-track analyses may be explained by blazars, that is, AGNs with bright Doppler-boosted jets. High-energy neutrinos can be produced in photohadronic interactions within parsec-scale relativistic jets. Radio-bright blazars associated with neutrino detections have very diverse gamma-ray properties, which suggests that gamma-rays and neutrinos may be produced in different regions of blazars and not directly related. A narrow jet viewing angle is, however, required to detect either of them. We conclude with discussion of recent independent tests and extensions of our findings.

Clusters of galaxies evolve and accrete mass, mostly from small galaxy systems. Our aim is to study the velocity field of the galaxy cluster Abell 780, which is known for the powerful radio source Hydra A at its center and where a spectacular X-ray tail associated with the galaxy LEDA 87445 has been discovered. Our analysis is based on the new spectroscopic data for hundreds of galaxies obtained with the Italian Telescopio Nazionale {\em Galileo} and the Very Large Telescope. We have constructed a redshift catalog of 623 galaxies and selected a sample of 126 cluster members. We analyze the internal structure of the cluster using a number of techniques. We estimate the mean redshift z=0.0545, the line-of-sight velocity dispersion sigmav about 800 km/s, and the dynamical mass M200 about 5.4 10E14 solar masses. The global properties of Abell 780 are typical of relaxed clusters. On a smaller scale, we can detect the presence of a galaxy group associated with LEDA 87445 in projected phase space. The mean velocity and position of the center of the group agree well with the velocity and position of LEDA 87445. We estimate the following parameters of the collision. The group is characterized by a higher velocity relative to the main system. It is infalling at a rest frame velocity of Vrf=+870 km/s and lies at a projected distance D=1.1 Mpc to the south, slightly southeast of the cluster center. The mass ratio between the group and the cluster is about 1:5. We also find evidence of an asymmetry in the velocity distribution of galaxies in the inner cluster region, which might be related to a small low-velocity group detected as a substructure at Vrf=-750 km/s. We conclude that A780, although dynamically relaxed at first sight, contains small substructures that may have some impact on the energetics of the core region.

We present a novel analysis for cluster cosmology that fully forward models the abundances, weak lensing, and the clustering of galaxy clusters. Our analysis notably includes an empirical model for the anisotropic boosts impacting the lensing and clustering signals of optical clusters. These boosts arise from a preferential selection of clusters surrounded by anisotropic large scale structure, a consequence of the limited discrimination between line-of-sight interlopers and true cluster members offered by photometric surveys. We validate our analysis via a blind cosmology challenge on mocks, and find that we can obtain tight and unbiased cosmological constraints without informative priors or external calibrations on any of our model parameters. We then apply our analysis on the SDSS redMaPPer clusters, and find results favoring low $\Omega_\mathrm{m}$ and high $\sigma_8$, combining to yield the lensing strength constraint $S_8 = 0.718_{-0.015}^{+0.016}$. We investigate potential drivers behind these results through a series of post-unblinding tests, noting that our results are consistent with existing cluster cosmology constraints but clearly inconsistent with other CMB/LSS based cosmology results. From these tests, we find hints that a suppression in the cluster lensing signal may be driving our results.

We report on the X-ray time resolved spectral analysis of XMM-Newton observations of NGC 3783. The main goal is to detect transient features in the Fe K line complex, in order to study the dynamics of the innermost accretion flow. We reanalize archival observations of NGC 3783, a bright local AGN, for which a transient Fe line was reported, complementing this data set with new available observations. This results in a long set of observations which can allow us to better assess the significance of transient features and possibly test their recurrence time. Moreover, since the new data catch the source in an obscured state, this analysis allows also to test whether the appearance/disappearance of transient features is linked to the presence of obscuring gas. We detect discrete features at the >=90% significance level both in emission and in absorption at different times of the observations, split into 5ks time-resolved spectra. The overall significance of individual features is higher in the obscured dataset. The energy distribution of the detections changes between the two states of the source, and the features appear to cluster at different energies. Counting the occurrences of emission/absorption lines at the same energies, we identify several groups of $\geq3\sigma$ detections: emission features in the 4-6 keV band are present in all observations and are most likely due to effects of the absorber present in the source; an emission line blend of neutral Fe K$\beta$/ionized Fe Ka is present in the unobscured dataset; absorption lines produced by gas at different ouflowing velocities and ionization states show an increase in energy between the two epochs, shifting from ~6 keV to ~6.7-6.9 keV. The representation of the features in a time-energy plane via residual maps highlighted a possible modulation of the Fe Ka line intensity, linked to the clumpiness of the absorbing medium.

Sub-subgiant stars (SSGs) fall below the subgiant branch and/or red of the giant branch in open and globular clusters, an area of the color-magnitude diagram (CMD) not populated by standard stellar evolution tracks. One hypothesis is that SSGs result from rapid rotation in subgiants or giants due to tidal synchronization in a close binary. The strong magnetic fields generated inhibit convection, which in turn produces large starspots, radius inflation, and lower-than-expected average surface temperatures and luminosities. Here we cross-reference a catalog of active giant binaries (RS CVns) in the field with Gaia EDR3. Using the Gaia photometry and parallaxes we precisely position the RS CVns in a CMD. We identify stars that fall below a 14 Gyr, metal-rich isochrone as candidate field SSGs. Out of a sample of 1723 RS CVn, we find 448 SSG candidates, a dramatic expansion from the 65 SSGs previously known. Most SSGs have rotation periods of 2-20 days, with the highest SSG fraction found among RS CVn with the shortest periods. The ubiquity of SSGs among this population indicates SSGs are a normal phase in evolution for RS CVn-type systems, not rare by-products of dynamical encounters found only in dense star clusters as some have suggested. We present our catalog of 1723 active giants, including Gaia photometry and astrometry, and rotation periods from TESS and VSX. This catalog can serve as an important sample to study the impacts of magnetic fields in evolved stars.

The Lyman-$\alpha$ absorption spectrum associated with photons traversing the intergalactic medium allows us to probe the linear matter power spectrum down to relatively small distance scales. Finding ways of accurately evaluating Lyman-$\alpha$ constraints across large classes of candidate models of dark-matter physics is thus of paramount importance. While such constraints have been evaluated for dark-matter models with relatively simple dark-matter velocity distributions, more complex models -- particularly those with dark-matter velocity distributions stretching across multiple scales -- are receiving increasing attention. In this paper, we undertake a study of the Lyman-$\alpha$ constraints associated with general dark-matter velocity distributions. Although these constraints are difficult to evaluate in principle, in practice there exist two ways of recasting them into forms which are easier to evaluate and which therefore allow a more rapid determination of whether a given dark-matter model is ruled in or out. We utilize both of these recasts in order to determine the Lyman-$\alpha$ bounds on different classes of dark-matter velocity distributions. We also develop a general method by which the results of these different recasts can be compared. For relatively simple dark-matter velocity distributions, we find that these two classes of recasts tend to align and give similar results. However, the situation is far more complex for distributions involving multiple velocity scales: while these two recasts continue to yield similar results within certain regions of parameter space, they nevertheless yield dramatically different results within precisely those regions of parameter space which are likely to be phenomenologically relevant. This, then, serves as a cautionary tale regarding the use of such recasts for complex dark-matter velocity distributions.

The complexity of atmospheric retrieval models is largely data-driven and one-dimensional models have generally been considered adequate with current data quality. However, recent studies have suggested that using 1D models in retrievals can result in anomalously cool terminator temperatures and biased abundance estimates even with existing transmission spectra of hot Jupiters. Motivated by these claims and upcoming high-quality transmission spectra we systematically explore the limitations of 1D models using synthetic and current observations. We use 1D models of varying complexity, both analytic and numerical, to revisit claims of biases when interpreting transmission spectra of hot Jupiters with inhomogeneous terminator compositions. Overall, we find the reported biases to be resulting from specific model assumptions rather than intrinsic limitations of 1D atmospheric models in retrieving current observations of asymmetric terminators. Additionally, we revise atmospheric retrievals of the hot Jupiter WASP-43b ($T_{\rm eq}=1440$ K) and the ultra-hot Jupiter WASP-103b ($T_{\rm eq}=2484$ K ) for which previous studies inferred abnormally cool atmospheric temperatures. We retrieve temperatures consistent with expectations. We note, however, that in the limit of extreme terminator inhomogeneities and high data quality some atmospheric inferences may conceivably be biased, although to a lesser extent than previously claimed. To address such cases, we implement a 2D retrieval framework for transmission spectra which allows accurate constraints on average atmospheric properties and provides insights into the spectral ranges where the imprints of atmospheric inhomogeneities are strongest. Our study highlights the need for careful considerations of model assumptions and data quality before attributing biases in retrieved estimates to unaccounted atmospheric inhomogeneities.

We consider the interplay of the Early Dark Energy (EDE) model, the Swampland Distance Conjecture (SDC), and cosmological parameter tensions. EDE is a proposed resolution of the Hubble tension relying upon a near-Planckian scalar field excursion, while the SDC predicts an exponential sensitivity of masses of other fields to such an excursion, $m\propto e^{-c|\Delta \phi|/M_{\rm pl}}$ with $c\sim{\cal O}(1)$. Meanwhile, EDE is in tension with large-scale structure (LSS) data, due to shifts in the standard $\Lambda$CDM parameters necessary to fit the cosmic microwave background (CMB). One might hope that a proper treatment of the model, e.g., accounting for the SDC, may ameliorate the tension with LSS. Motivated by these considerations, we introduce the Early Dark Sector (EDS) model, wherein the mass of dark matter is exponentially sensitive to super-Planckian field excursions of the EDE scalar. The EDS model exhibits new phenomenology in both the early and late universe, the latter due to an EDE-mediated dark matter self-interaction. This dark matter-philic "fifth force", while constrained to be small, remains active in the late universe and is not screened in virialized halos. We find that the new interaction with dark matter partially resolves the LSS tension. However, the marginalized posteriors are nonetheless consistent with $f_{\rm EDE}=0$ at 95$\%$ CL once the Dark Energy Survey Year 3 measurement of $S_8$ is included. We study constraints on the model from Atacama Cosmology Telescope data, and find a factor of two improvement on the error bar on the SDC parameter $c$, along with an increased preference for the EDE component. We discuss the implications of these constraints for the SDC, and find the tightest observational constraints to date on a swampland parameter, suggesting that an EDE description of cosmological data is in tension with the SDC.

Feebly interacting particles with masses with O(10-100) MeV can be copiously produced by core-collapse supernovae (SNe). In this paper we consider the case of MeV-ish sterile neutrinos and dark photons mixed with ordinary neutrinos and photons, respectively. Furthermore, both sterile neutrinos and dark photons may decay into positrons on their route to Earth. Such positrons would annihilate with electrons in the Galactic medium and contribute to the photon flux in the 511 keV line. Using the SPI (SPectrometer on INTEGRAL) observation of this line improves the bounds on the mixing parameters for these particles by several orders of magnitude below what is already excluded by the SN 1987A energy-loss argument.

The approach of dynamical systems is a useful tool to investigate the cosmological history that follows from modified theories of gravity. It provides qualitative information on the typical background solutions in a parametrized family of models, through the computation of the fixed points and their characters (attractor, repeller or saddle), allowing, for instance, the knowledge of which regions on the parameter space of the models generate the desired radiation, matter and dark energy dominated eras. However, the traditional proposal for building dynamical systems for an $f(R)$ theory in the Palatini formalism assumes the invertibility of a function that depends on the specific Lagrangian functional form, which is not true, for example, for the particular theory of exponential gravity ($f(R)=R-\alpha R_*(1-e^{-R/R_*})$). In this work, we propose an alternative choice of variables to treat $f(R)$ models in their Palatini formulation, which do include exponential gravity. We derive some general results that can be applied to a given model of interest and present a complete description of the phase space for exponential gravity. We show that Palatini exponential gravity theories have a final attractor critical point with an effective equation of state parameter $w_{\text{eff}} = -1$ (for $\alpha>1$), $w_{\text{eff}} = -2/3$ (for $\alpha=1$) and $w_{\text{eff}} = 0$ (for $\alpha<1$). Finally, our analytical results are compared with numerical solutions of the field equations.

The detection of a binary neutron star merger in 2017 through both gravitational waves and electromagnetic emission opened a new era of multimessenger astronomy. The understanding of the magnetic field amplification triggered by the Kelvin-Helmholtz instability during the merger is still a numerically unresolved problem because of the relevant small scales involved. One of the uncertainties comes from the simplifications usually assumed in the initial magnetic topology of merging neutron stars. We perform high-resolution, convergent large-eddy simulations of binary neutron star mergers, following the newly formed remnant for up to $30$ milliseconds. Here we specifically focus on the comparison between simulations with different initial magnetic configurations, going beyond the widespread-used aligned dipole confined within each star. The results obtained show that the initial topology is quickly forgotten, in a timescale of few miliseconds after the merger. Moreover, at the end of the simulations, the average intensity ($B\sim 10^{16}$ G) and the spectral distribution of magnetic energy over spatial scales barely depend on the initial configuration. This is expected due to the small-scale efficient dynamo involved, and thus it holds as long as: (i) the initial large-scale magnetic field is not unrealistically high (as often imposed in mergers studies); (ii) the turbulent instability is numerically (at least partially) resolved, so that the amplified magnetic energy is distributed across a wide range of scales and becomes orders-of-magnitude larger than the initial one.

Magnetic fields are expected to play a key role in the dynamics and the ejection mechanisms that accompany the merger of two neutron stars. General relativistic magnetohydrodynamic (MHD) simulations offer a unique opportunity to unravel the details of the ongoing physical processes. Nevertheless, current numerical studies are severely limited by the fact that any affordable resolution remains insufficient to fully capture the small-scale dynamo, initially triggered by the Kelvin-Helmholtz instability, and later sourced by several MHD processes involving differential rotation. Here, we alleviate this limitation by using explicit large-eddy simulations, a technique where the unresolved dynamics occurring at the sub-grid scales (SGS) is modeled by extra terms, which are functions of the resolved fields and their derivatives. The combination of high-order numerical schemes, high resolutions, and the gradient SGS model allow us to capture the small-scale dynamos produced during the binary neutron star mergers. Here we follow the first 50 milliseconds after the merger and, for the first time, we find numerical convergence on the magnetic field amplification, in terms of integrated energy and spectral distribution over spatial scales. We also find that the average intensity of the magnetic field in the remnant saturates at $\sim 10^{16}$~G around $5$~ms after the merger. After $20-30$~ms, both toroidal and poloidal magnetic field components grow continuously, fed by the winding mechanism that provides a slow inverse cascade. We find no clear hints for magneto-rotational instabilities, and no significant impact of the magnetic field on the redistribution of angular momentum in the remnant in our simulations, probably due to the very turbulent and dynamical topology of the magnetic field at all stages, with small-scale components largely dominating over the large-scale ones.

Astronomical transients are stellar objects that become temporarily brighter on various timescales and have led to some of the most significant discoveries in cosmology and astronomy. Some of these transients are the explosive deaths of stars known as supernovae while others are rare, exotic, or entirely new kinds of exciting stellar explosions. New astronomical sky surveys are observing unprecedented numbers of multi-wavelength transients, making standard approaches of visually identifying new and interesting transients infeasible. To meet this demand, we present two novel methods that aim to quickly and automatically detect anomalous transient light curves in real-time. Both methods are based on the simple idea that if the light curves from a known population of transients can be accurately modelled, any deviations from model predictions are likely anomalies. The first approach is a probabilistic neural network built using Temporal Convolutional Networks (TCNs) and the second is an interpretable Bayesian parametric model of a transient. We show that the flexibility of neural networks, the attribute that makes them such a powerful tool for many regression tasks, is what makes them less suitable for anomaly detection when compared with our parametric model.

In this review article, we describe the role of holography in deciphering the physics of dense QCD matter, relevant for the description of compact stars and their binary mergers. We review the strengths and limitations of the holographic duality in describing strongly interacting matter at large baryon density, walk the reader through the most important results derived using the holographic approach so far, and highlight a number of outstanding open problems in the field. Finally, we discuss how we foresee holography contributing to compact-star physics in the coming years.

Non-thermalized dark matter is a cosmologically viable alternative to the widely studied weakly interacting massive particle. We study the evolution of the dark matter phase-space distributions arising from freeze-in and superWIMP production as well as the combination of both. Utilizing our implementation in CLASS, we investigate the cosmological imprints on the matter power spectrum, constrained by Lyman-$\alpha$ forest observations. For the explicit example of a colored $t$-channel mediator model, we explore the cosmologically allowed parameter space highlighting the interplay of Lyman-$\alpha$ constraints with those from Big Bang Nucleosynthesis and the LHC.

The extensive focus on performance indicators in research evaluation has been facing critique in science studies. Stemming from a neoliberalist paradigm, metrics allegedly objectify and create certainty about researchers' performance. This has created a publish-or-perish culture where deviant behaviour, such as research misconduct, may have become the rule of the game. Not only does this culture foster a decrease of scientists' well-being, but also a decrease in research quality. In recent years, calls for a culture change have accumulated, from discussing detrimental cultural aspects under IchbinHannah to studies that demonstrate the connection between research culture and research integrity. This study is a qualitative analysis of how astronomers reimagine their research culture. This includes alternative output formats, alternative evaluation criteria and how they aspire to do research differently. In summary, we find that the time is ripe for a transformation in research culture towards a more participative and diverse work environment. This may include the use of an open knowledge management infrastructure, where all sorts of research output of various stages can be stored and shared. Moreover, through a more reflexive evaluation, which is continuously adapted to the needs of the researchers, scientific quality may be encouraged, rather than producing more and more publications. This study sets the basis for future action research with the aim of transforming academic cultures towards more participative ones, where scientists can let their creative minds thrive and collaborate beyond disciplinary boundaries.

We propose a new mechanism to simultaneously explain the observed dark matter abundance and the baryon asymmetry of the Universe. The mechanism is based on the Filtered Dark Matter scenario, where dark matter particles acquire a large mass during a first-order phase transition. This implies that only a small fraction of them are energetic enough to enter the advancing true vacuum bubbles and survive until today, while the rest are reflected and annihilate away quickly. We supplement this scenario with a CP-violating interaction, which creates a chiral asymmetry in the population of dark matter particles. In the false vacuum phase, a portal interaction quickly converts the dark sector chiral asymmetry into a Standard Model lepton asymmetry. The lepton asymmetry is then partially converted to a baryon asymmetry by standard electroweak sphaleron processes. We discuss the dependence of the generated asymmetry on the parameters of the model for two different portal interactions and demonstrate successful baryogenesis for both. For one of the portals, it is also possible to simultaneously explain the observed dark matter abundance, over many orders of magnitude in the dark matter mass.

The E989 experiment at the Fermi National Laboratory reported a 4.2$\sigma$ discrepancy between the measured magnetic dipole moment of the muon, and its prediction in the Standard Model (SM). In this study, we address the anomaly by considering a minimal and generic extension to the SM which also provides for a dark matter (DM) candidate. The extra states in this framework are: a SM singlet Majorana fermion, referred to as the Bino, playing the role of DM; and muonic scalars, referred to as sleptons. The couplings between the sleptons, SM muons and the Bino can account for the muon $g-2$ anomaly if the scalar muon partners, or smuons, mix chirality. On the other hand, the DM relic density is satisfied primarily through coannihilation effects involving the Bino and the lighter sleptons. The viable parameter space of our model includes regions with relatively light coannihilating particles, similar to what has been found in previous scans of the Minimal Supersymmetric Standard Model (MSSM). Relaxing the assumption of minimal flavor violation typically assumed in the MSSM, we see that scenarios with sizable smuon mixing and large mass splittings between the smuons can satisfy both the muon $g-2$ anomaly and the DM relic density for coannihilating particle masses up to and beyond the TeV scale. When we specify the origin of the left-right smuon mixing to be trilinear couplings between the smuons and the SM Higgs boson, the constraints on these scenarios arising from perturbative unitarity and electroweak vacuum stability confine the coannihilating particle masses to be $\lesssim$ 1 TeV. We demonstrate that next generation direct detection experiments are only marginally sensitive to the viable parameter space of our model and, thus, a future lepton collider could be the essential probe necessary to distinguish our model from other BSM solutions to the muon $g-2$ anomaly.

The Weakly Interacting Massive Particle (WIMP) paradigm is one of the most popular scenarios for Dark Matter (DM) theories that however is strongly constrained, in particular by direct detection experiments. We stick with the WIMP hypothesis and consider a Dirac fermion candidate for DM that interacts with the Standard Model (SM) via a spin-1 $Z'$, arising from the spontaneous breaking of an Abelian $U(1)'_{\mu}$ gauge symmetry, under which only second generation leptons and the DM are appropriately charged. Due to the charge assignment, the model is gauge anomalous and can only be interpreted as an effective field theory (EFT) at low energy. The $Z'$ couples at tree level only to the vector DM current, to the axial muon current and to left-handed muonic neutrinos, so the WIMP-nucleon cross section is beyond the experimental reach of spin-independent (SI) direct detection searches. We study the current bounds on this model coming from direct and indirect detection of DM, collider searches, contributions to $(g-2)_{\mu}$ and to neutrino trident production. We find that large regions of the parameter space remains to be explored. In the context of LHC searches, we study the impact of a muon-exclusive signal region for the $3\mu$ + ${E}^{{\rm miss}}_T$ channel with an invariant mass window around $m_{Z'}$. We show that this search can significantly improve the current collider bounds. Finally, from the anomalous nature of our EFT, there remain at low energy triboson anomalous interactions between the $Z'$ and the electroweak (EW) SM gauge bosons. We explore the possibilities of probing these interactions at the LHC and at a 100 TeV proton collider finding it extremely challenging. On the other hand, for a muon collider the resonant channel $\mu^{+}\mu^{-}\to Z'\to ZZ$ could be discovered in the most promising scenario with luminosity of $\mathcal{O}({\rm few}\; 10)$ ${\rm fb}^{-1}$.

We add to the Standard Model a new fermion $\chi$ with minimal baryon number 1/3. Neutron decay $n \to \chi\chi\chi$ into non-relativistic $\chi$ can account for the neutron decay anomaly, compatibly with bounds from neutron stars. $\chi$ can be Dark Matter, and its cosmological abundance can be generated by freeze-in dominated at $T \sim m_n$. The associated processes $n \to \chi\chi\chi \gamma$, hydrogen decay ${\rm H}\to \chi\chi \chi\nu(\gamma)$ and DM-induced neutron disappearance $\bar\chi n \to \chi \chi (\gamma)$ have rates below experimental bounds and can be of interest for future experiments.

We demonstrate the existence of Q-balls in non-minimally coupled inflation models with a complex inflaton in the Palatini formulation of gravity. We show that there exist Q-ball solutions which are compatible with inflation and we derive a window in the inflaton mass squared for which this is the case. In particular, we confirm the existence of Q-ball solutions with $\phi \sim 10^{17}-10^{18} \GeV $, consistent with the range of field values following the end of slow-roll Palatini inflation. We study the Q-balls and their properties both numerically and in an analytical approximation. The existence of such Q-balls suggests that the complex inflaton condensate can fragment into Q-balls, and that there may be an analogous process for the case of a real inflaton with fragmentation to neutral oscillons. We discuss the possible post-inflationary cosmology following the formation of Q-balls, including an early Q-ball matter domination (eMD) period and the effects of this on the reheating dynamics of the model, gravitational wave signatures which may be detectable in future experiments, and the possibility that Q-balls could lead to the formation of primordial black holes (PBHs). In particular, we show that Palatini Q-balls with field strengths typical of inflaton condensate fragmentation can directly form black holes with masses around 500 kg or more when the self-coupling is $\lambda = 0.1$, resulting in very low (less than 100 GeV) reheating temperatures from black hole decay, with smaller black hole masses and larger reheating temperatures possible for smaller values of $\lambda$. Q-ball dark matter from non-minimally coupled Palatini inflation may also be a direction for future work.