Papers for Wednesday, Sep 29 2021

Active galactic nuclei (AGN) forming in the early universe are thought to be the primary source of hard ionizing photons contributing to the reionization of intergalactic helium. However, the number density and spectral properties of high-redshift AGN remain largely unconstrained. In this work, we make use of physically-informed models calibrated with a wide variety of available observations to provide estimates for the role of AGN throughout the Epoch of Reionization. We present AGN luminosity functions in various bands between z = 2 to 7 predicted by the well-established Santa Cruz semi-analytic model, which includes modelling of black hole accretion and AGN feedback. We then combine the predicted AGN populations with a physical spectral model for self-consistent estimates of ionizing photon production rates, which depend on the mass and accretion rate of the accreting supermassive black hole. We then couple the predicted comoving ionizing emissivity with an analytic model to compute the subsequent reionization history of intergalactic helium and hydrogen. This work demonstrates the potential of coupling physically motivated analytic or semi-analytic techniques to capture multi-scale physical processes across a vast range of scales (here, from AGN accretion disks to cosmological scales). Our physical model predicts an intrinsic ionizing photon budget well above many of the estimates in the literature, meaning that helium reionization can comfortably be accomplished even with a relatively low escape fraction. We also make predictions for the AGN populations that are expected to be detected in future \emph{James Webb Space Telescope} surveys.

During a failed core-collapse supernova, the protoneutron star eventually collapses under its own gravitational field and forms a black hole. This collapse happens quickly, on the dynamical time of the protoneutron star, $\lesssim$0.5 ms. During this collapse, barring any excessive rotation, the entire protoneutron star is accreted into the newly formed black hole. The main source of neutrinos is now removed and the signal abruptly shuts off over this formation timescale. However, while the source of neutrinos is turned off, the arrival times at an Earth-based detector will depend on the neutrino path. We show here that a modest amount of neutrinos, emitted just prior to the black hole forming, scatter on the infalling material into our line of sight and arrive after the formation of the black hole, up to 15 ms in our model. This neutrino echo, which we characterize with Monte Carlo simulations and analytic models, has a significantly higher average energy (upwards of $\sim$ 50 MeV) compared to the main neutrino signal, and for the canonical failed supernova explored here, is likely detectable in $\mathcal{O}$(10 kT) supernova neutrino detectors for Galactic failed supernovae. The presence of this signal is important to consider if using black hole formation as a time post for triangulation or the post black hole timing profile for neutrino mass measurements. On its own, it can also be used to characterize or constrain the structure and nature of the accretion flow.

Using the ARTEMIS set of 45 high-resolution cosmological simulations, we investigate a range of merger-induced dynamical transformations of Milky Way-like galaxies. We first identify populations of accreted stars on highly radial orbits, similar to the 'Gaia Sausage' in the Milky Way. We show that $\approx1/3$ of the ARTEMIS galaxies contain a similar feature, and confirm that they usually comprise stellar debris from the most massive accreted satellite. Selecting these 15 galaxies, we study their changes around the times of the GS-like mergers. Dark matter haloes of many of these exhibit global changes in shape and orientation. Focusing on the galaxies themselves, we find multiple examples of stellar discs whose angular momentum (AM) axes change orientation at rapid rates of $\sim60$ degrees Gyr$^{-1}$. By calculating the orbital angular momentum axes of the satellites before they are accreted, we show that there is a tendency for the disc's AM to become more aligned with this axis after the merger. We also investigate the origin of in situ retrograde stars, analogous to the 'Splash' in the Milky Way. Tracing them back to earlier snapshots, we demonstrate that they were often disrupted onto their extreme orbits by multiple early mergers. We also find that the total mass of these stars outside the central regions positively correlates with the total accreted stellar mass. Finally, we conduct a brief investigation into whether bars form soon after the mergers. In a few galaxies we find a bar-like feature whose emergence coincides with a significant merger.

Context. Existing samples of strong lenses have been assembled by giving priority to sample size, at the cost of having a complex selection function. With the advent of the next generation of wide-field photometric surveys, however, it might become possible to identify subsets of the lens population with well-defined selection criteria, trading sample size for completeness. Aims. There are two main advantages of working with a complete sample of lenses. First, it is possible to recover the properties of the general population of galaxies, of which strong lenses are a biased subset. Second, the relative number of lenses and non-detections can be used to further constrain models of galaxy structure. This work illustrates how to carry out a statistical strong lensing analysis that takes advantage of these features. Methods. I introduced a general formalism for the statistical analysis of a sample of strong lenses with known selection function, then tested it on simulated data. The simulation consists of a population of $10^5$ galaxies with an axisymmetric power-law density profile, a population of background point sources, and a subset of $\sim10^3$ strong lenses, complete above an observational cut. Results. The method allows to recover the distribution in Einstein radius and mass density slope of the galaxy population in an unbiased way. The number of non-lenses helps to constrain the model when magnification data are not available. Conclusions. Complete samples of lenses are a powerful asset to turn precise strong lensing measurements into accurate statements on the properties of the general galaxy population.

Exoplanetary observations reveal that the occurrence rate of hot Jupiters is correlated with star clustering. In star clusters, interactions between planetary systems and close fly-by stars can significantly change the architecture of primordially coplanar, circular planetary systems. Flybys in dense clusters have a significant impact on hot Jupiter formation via activation of high eccentricity excitation mechanisms such as the Zeipel-Lidov-Kozai (ZLK) effect and planet-planet scattering. Previous studies have shown that if there are two giant planets in the planetary system, close flybys can efficiently activate the ZLK mechanism, thus triggering high eccentricity tidal migration and ultimately form hot Jupiters in star clusters. Here we extend our previous study with a multi-planet (triple) system. We perform high precision, high-accuracy few-body simulations of stellar flybys and subsequent planetary migration within the perturbed planetary systems using the code {\tt SpaceHub}. Our simulations demonstrate that a single close flyby on a multi-planet system in a cluster can activate secular chaos and ultimately lead to hot Jupiter formation via high eccentricity migration. We find that the hot Jupiter formation rate per system increases with both the size of the planetary system as well as with the mass of the outer planet, and we quantify the relative formation fractions for a range of parameters. Hot Jupiters formed via secular chaos are expected to be accompanied by massive companions with very long periods. Our study further shows that this flyby-induced secular chaos is preferred in low-density clusters where multi-planet systems are more likely to survive, and that it contributes a significant fraction of the hot Jupiter formation in star clusters compared to the flyby-induced ZLK mechanism.

Recently, cosmic ray (CR) feedback has been identified as a critical process in galaxy formation but most previous simulations have integrated out the energy-dependence of the CR distribution, despite its large extent over more than twelve decades in particle energy. To improve upon this simplification, we present the implementation and first application of spectrally resolved CRs which are coupled to the magneto-hydrodynamics in simulations of galaxy formation. The spectral model for the CRs enables more accurate cooling of CRs and allows for an energy-dependent spatial diffusion, for which we introduce a new stable numerical algorithm that proves essential in highly dynamical systems. We perform galaxy formation simulations with this new model and compare the results to a grey CR approach with a simplified diffusive transport and effective cooling that assumes steady-state spectra. We find that the galaxies with spectrally resolved CRs differ in morphology, star formation rate, and strength and structure of the outflows. Interestingly, the first outflow front is driven by CRs with average momenta of $\sim200-600\,\mathrm{Gev}~c^{-1}$. The subsequent formation of outflows, which reach mass loading factors of order unity, are primarily launched by CRs of progressively smaller average momenta of $\sim8-15\,\mathrm{GeV}~c^{-1}$. The CR spectra in the galactic centre quickly approach a steady state, which does not significantly vary over time. In the outer disc and outflow regions, the spectral shape approaches steady state only after $\sim2\,\mathrm{Gyr}$ of evolution. Furthermore, the shapes of the approximate steady state spectra differ for individual regions of the galaxy, which highlights the importance of actively including the full CR spectrum.

We explore the evolution of the flux power spectrum in the Cosmic Reionization On Computers (CROC) simulations. We find that, contrary to some previous studies, the shape of the flux power spectrum is rather insensitive to the timing of reionization. However, the amplitude of the flux power spectrum does strongly evolve with time, and that evolution is almost perfectly correlated with the timing of reionization. We show how such correlation can be used in a (futuristic) measurement to determine the redshift of overlap of ionized bubbles.

Half of the satellite galaxies of Andromeda form a narrow plane termed the Great Plane of Andromeda (GPoA), and their line-of-sight velocities display correlation reminiscent of a rotating structure. Recently reported first proper motion measurements for the on-plane satellites NGC 147 and NGC 185 indicate that they indeed co-orbit along the GPoA. This provides a novel opportunity to compare the M31 satellite system to $\Lambda$CDM expectations. We perform the first detailed comparison of the orbital alignment of two satellite galaxies beyond the Milky Way with several hydrodynamical and dark-matter-only cosmological simulations (Illustris TNG-50, TNG-100, ELVIS, PhatELVIS), in the context of the Planes of Satellite Galaxies Problem. In line with previous works, we find that the spatial flattening and line-of-sight velocity correlation alone is already in substantial tension with $\Lambda$CDM, with none of the simulated analogs simultaneously reproducing both parameters. Almost none (3 to 4\%) of the simulated systems contain two satellites with orbital poles as well aligned with their satellite plane as indicated by the most-likely proper motions of NGC 147 and NGC 185. However, within current measurement uncertainties it is common (~70%) that the two best-aligned satellites of simulated systems are consistent with the orbital alignment. Yet, the chance that any two simulated on-plane satellites have as well aligned orbital poles as observed is low (~4%). We conclude that confirmation of the tight orbital alignment for these two objects via improved measurements, or the discovery of similar alignments for additional GPoA members, hold the potential to further raise the tension with $\Lambda$CDM expectations.

Gas and dust in inclined orbits around binaries experience precession induced by the binary gravitational torque. The difference in precession between gas and dust alters the radial drift of weakly coupled dust and leads to density enhancements where the radial drift is minimised. We explore this phenomenon using 3D hydrodynamical simulations to investigate the prominence of these `dust traffic jams' and the evolution of the resulting dust sub-structures at different disc inclinations and binary eccentricities. We then derive evolution equations for the angular momentum of warped dust discs and implement them in a 1D code and present calculations to further explain these traffic jams. We find that dust traffic jams in inclined circumbinary discs provide significant dust density enhancements that are long lived and can have important consequences for planetesimal formation.

Large spectroscopic surveys of the Milky Way need to be calibrated against a sample of benchmark stars to ensure the reliable determination of atmospheric parameters. We present new fundamental stellar parameters of seven giant and subgiant stars that will serve as benchmarks. The aim is to reach a precision of 1% in the effective temperature. This precision is essential for accurate determinations of the full set of fundamental parameters and abundances of stars observed by the surveys. We observed HD121370 (etaBoo), HD161797 (muHer), HD175955, HD182736, HD185351, HD188512 (betaAql), and HD189349 using the high angular resolution optical interferometric instrument PAVO/CHARA. The limb-darkening corrections were determined from 3D model atmospheres based on the STAGGER grid. The Teff were determined directly from the Stefan-Boltzmann relation, with an iterative procedure to interpolate over tables of bolometric corrections. We estimated surface gravities from comparisons to Dartmouth stellar evolution model tracks. The spectroscopic observations were collected from the ELODIE and FIES spectrographs. We estimated metallicities ([Fe/H]) from a 1D non-local thermodynamic equilibrium (NLTE) abundance analysis of unblended lines of neutral and singly ionised iron. For six of the seven stars we measure Teff to better than 1%. For one star, HD189349, the uncertainty in Teff is 2% due to an uncertain bolometric flux. We do not recommend this star as a benchmark until this measurement can be improved. Median uncertainties for all stars in logg and [Fe/H]} are 0.034dex and 0.07dex, respectively. All the fundamental stellar parameters were based on consistently combining interferometric observations, 3D limb-darkening modelling and spectroscopic analysis. This paper follows our previous papers including dwarfs and metal-poor stars.

The shock breakout emission is the first light that emerges from a supernova. In the spherical case it is characterized by a brief UV flash. In an axisymmetric, non-spherical prolate explosion, the shock first breaches the surface along the symmetry axis, then peels around to larger angles, producing a breakout light curve which may differ substantially from the spherically symmetric case. We study the emergence of a non-relativistic, bipolar shock from a spherical star, and estimate the basic properties of the associated bolometric shock breakout signal. We identify four possible classes of breakout light curves, depending on the degree of asphericity. Compared to spherical breakouts, we find that the main distinguishing features of significantly aspherical breakouts are 1) a longer and fainter initial breakout flash and 2) an extended phase of slowly-declining, or even rising, emission which is produced as ejecta flung out by the oblique breakout expand and cool. We find that the breakout flash has a maximum duration of roughly $\sim R_*/v_{\rm bo}$ where $R_*$ is the stellar radius and $v_{\rm bo}$ is the velocity of the fastest-moving ejecta. For a standard Wolf--Rayet progenitor, the duration of the X-ray flash seen in SN 2008D exceeds this limit, and the same holds true for the prompt X-ray emission of low-luminosity GRBs such as GRB 060218. This suggests that these events cannot be explained by an aspherical explosion within a typical Wolf--Rayet star, implying that they originate from non-standard progenitors with larger breakout radii.

Fitting parameterized models to images of galaxies has become the standard for measuring galaxy morphology. This forward modelling technique allows one to account for the PSF to effectively study semi-resolved galaxies. However, using a specific parameterization for a galaxy's surface brightness profile can bias measurements if it is not an accurate representation. Furthermore, it can be difficult to assess systematic errors in parameterized profiles. To overcome these issues we employ the Multi-Gaussian expansion (MGE) method of representing a galaxy's profile together with a Bayesian framework for fitting images. MGE flexibly represents a galaxy's profile using a series of Gaussians. We introduce a novel Bayesian inference approach which uses pre-rendered Gaussian components, which greatly speeds up computation time and makes it feasible to run the fitting code on large samples of galaxies. We demonstrate our method with a series of validation tests. By injecting galaxies, with properties similar to those observed at $z\sim1.5$, into deep HST observations we show that it can accurately recover total fluxes and effective radii of realistic galaxies. Additionally we use degraded images of local galaxies to show that our method can recover realistic galaxy surface brightness and color profiles. Our implementation is available in an open source python package $\texttt{imcascade}$, which contains all methods needed for the preparation of images, fitting and analysis of results.

We image the spatial extent of a cool galactic outflow with resonant Mg II emission in a gravitationally lensed star-forming galaxy at $z \sim 1.7$ using VLT/MUSE observations. We observe Mg II residual (continuum-subtracted) emission out to an observed radial distance of $26.5_{-0.4}^{+0.5}$ kpc from the galaxy, with an observed maximum spatial extent of $\approx39_{-0.6}^{+0.8}$ kpc ($30_{-0.5}^{+0.7}$ kpc after correcting for seeing). Mg II residual emission is patchy and covers a total area of ${\rm 184_{-10}^{+5} kpc^2}$, constraining the minimum area covered by the outflowing gas to be $13.27_{-1.02}^{+0.55}$ % of the total area. The spatial extent of the Mg II emission is asymmetric and shows an observed $27.6_{-0.7}^{+0.8}$ % ($20.9_{-0.6}^{+0.7}$ % after seeing correction) larger extent along the declination direction. We constrain the covering fraction of the Mg II emission as a function of radial distance and characterize it with a power-law model convolved with the seeing with an index $\gamma= -1.25_{-0.02}^{+0.02}$. We find two kinematically distinct Mg II emission components ($\Delta v \approx400\ {\rm km\ s^{-1}}$) which extend out to similar distances, and may correspond to two distinct shells of outflowing gas. By using multiple images with different magnifications of the galaxy in the image plane, we trace the Mg II residual emission in three individual star-forming regions inside the galaxy out to $6.0_{- 0.2 }^{+0.2}$ , $7.0_{-0.2}^{+0.3}$ , and $8.5_{-0.1}^{+0.1}$ kpc. Both the Fe II* fine structure emission, and the nebular [O II] emission are not spatially extended relative to the stellar continuum. These findings provide robust constraints on the spatial extent of the outflowing gas and combined with outflow velocity and column density measurements will give stringent constraints on mass outflow rates of the galaxy.

Measurements of secondary cosmic microwave background (CMB) anisotropies, such as the Sunyaev-Zel'dovich (SZ) effect, will enable new tests of neutrino and dark sector properties. The kinetic SZ (kSZ) effect is produced by cosmological flows, probing structure growth. Ultra-light axions (ULAs) are a well-motivated dark-matter candidate. Here the impact of ULA dark matter (with mass $10^{-27}~{\rm eV}$ to $10^{-23}~{\rm eV}$) on kSZ observables is determined, applying new analytic expressions for pairwise cluster velocities and Ostriker-Vishniac signatures in structure-suppressing models. For the future CMB-S4 and ongoing DESI galaxy surveys, the kSZ effect (along with primary anisotropies) will probe ULA fractions $\eta_a = \Omega_{\rm{axion}}/\Omega_{\rm DM}$ as low as $\sim 5\%$ if $m_{a}\simeq 10^{-27}~{\rm eV}$ (at 95\% C.L.), with sensitivity extending up to $m_{a}\simeq 10^{-25}~{\rm eV}$. If reionization and the primary CMB can be adequately modeled, Ostriker-Vishniac measurements could probe values $\eta_{a}\simeq 10^{-3}$ if $10^{-27}~{\rm eV}\lesssim m_{a}\lesssim 10^{-24}~{\rm eV}$, or $\eta_{a}\simeq 1$ if $m_{a}\simeq 10^{-22}~{\rm eV}$, within the fuzzy dark matter window.

Using a suite of 3D hydrodynamical simulations of star-forming molecular clouds, we investigate how the density probability distribution function (PDF) changes when including gravity, turbulence, magnetic fields, and protostellar outflows and heating. We find that the density PDF is not lognormal when outflows and self-gravity are considered. Self-gravity produces a power-law tail at high densities and the inclusion of stellar feedback from protostellar outflows and heating produces significant time-varying deviations from a lognormal distribution at the low densities. The simulation with outflows has an excess of diffuse gas compared to the simulations without outflows, exhibits increased average sonic Mach number, and maintains a slower star formation rate over the entire duration of the run. We study the mass transfer between the diffuse gas in the lognormal peak of the PDF, the collapsing gas in the power-law tail, and the stars. We find that the mass fraction in the power-law tail is constant, such that the stars form out of the power-law gas at the same rate at which the gas from the lognormal part replenishes the power-law. We find that turbulence does not provide significant support in the dense gas associated with the power-law tail. When including outflows and magnetic fields in addition to driven turbulence, the rate of mass transfer from the lognormal to the power-law, and then to the stars, becomes significantly slower, resulting in slower star formation rates and longer depletion times.

We reanalyze the multi-epoch direct detections of HD 88133 b and ups And b that were published in Piskorz et al. 2016 and Piskorz et al. 2017, respectively. Using simulations to attempt to reproduce the detections, we find that with the 6 and 7 $L$ band Keck/NIRSPEC epochs analyzed in the original works, the planets would not have been detectable unless they had unreasonably large radii. HD88133 and ups And both have fairly large stellar radii, which contributed to the difficulty in detecting the planets. We take this opportunity to consider how these planets may have been detectable with the small number of epochs originally presented by running simulations both with the upgraded NIRSPEC instrument and with near-zero primary velocities, as recommended by Buzard et al. 2021. While 7 $L$ band NIRSPEC2.0 epochs with near-zero primary velocities would have allowed a strong ($10.8\sigma$) detection of ups And b, many more than 6 $L$ band epochs would have been required for a strong detection of HD88133b, which could be due in part to both this system's large stellar radius and low stellar temperature. This work stresses the importance of careful analytic procedures and the usefulness of simulations in understanding the expected sensitivity of high-resolution spectroscopic data.

Filamentary structures identified in far-infrared observations of molecular clouds are typically found to have full-widths at half-maximum $\sim\!0.1$ pc. However, the physical explanation for this phenomenon is currently uncertain. We use hydrodynamic simulations of cylindrically-symmetric converging flows to show that the full-width at half-maximum of the resulting filament's surface density profile, FWHM$_{\Sigma}$, is closely related to the location of the accretion shock, where the inflow meets the boundary of the filament. For inflow Mach Number, ${\cal M}$, between 1 and 5, filament FWHM$_{\Sigma}$s fall in the range $0.03$ pc $\lesssim$ FWHM$_{\Sigma} \lesssim 0.3$ pc, with higher ${\cal M}$ resulting in narrower filaments. A large sample of filaments, seen at different evolutionary stages and with different values of ${\cal M}$, naturally results in a peaked distribution of FWHM$_{\Sigma}$s similar in shape to that obtained from far-infrared observations of molecular clouds. However, unless the converging flows are limited to ${\cal M} \lesssim 3$, the peak of the distribution of FWHM$_{\Sigma}$s is below the observed $\sim 0.1$ pc.

Identifying stars formed in pristine environments (Pop III) within the first billion years is vital to uncovering the earliest growth and chemical evolution of galaxies. Pop III galaxies, however, are typically expected to be too faint and too few in number to be detectable by forthcoming instruments without extremely long integration times and/or extreme lensing. In an environment, however, where star formation is suppressed until a halo crosses the atomic cooling limit (e.g., by a modest Lyman-Werner flux, high baryonic streaming velocities, and/or dynamical heating effects),primordial halos can form substantially more numerous and more massive stars. Some of these stars will in-turn be accreting more rapidly than they can thermally relax at any given time. Using high resolution cosmological zoom-in simulations of massive star formation in high-z halos, we find that such rapidly accreting stars produce prominent spectral features which would be detectable by {\it JWST}. The rapid accretion episodes within the halo lead to stochastic reprocessing of 0--20\% of the total stellar emission into the rest-frame optical over long timescales, a unique signature which may allow deep observations to identify such objects out to $z \sim 10-13$ using mid- and wide-band NIRCam colors alone.

We present the results of simultaneous observations of the transitional millisecond pulsar (tMSP) candidate CXOU J110926.4-650224 with the XMM-Newton satellite and the MeerKAT telescope. The source was found at an average X-ray luminosity of L_{\rm X}\simeq7\times10^{33} erg s^{-1} over the 0.3-10 keV band (assuming a distance of 4 kpc) and displayed a peculiar variability pattern in the X-ray emission, switching between high, low and flaring modes on timescales of tens of seconds. A radio counterpart was detected at a significance of 7.9\sigma with an average flux density of \simeq33\muJy at 1.28 GHz. It showed variability over the course of hours and emitted a \simeq10-min long flare just a few minutes after a brief sequence of multiple X-ray flares. No clear evidence for a significant correlated or anticorrelated variability pattern was found between the X-ray and radio emissions over timescales of tens of minutes and longer. CXOU J110926.4-650224 was undetected at higher radio frequencies in subsequent observations performed with the Australia Telescope Compact Array, when the source was still in the same X-ray sub-luminous state observed before, down to a flux density upper limit of 15\muJy at 7.25 GHz (at 3\sigma). We compare the radio emission properties of CXOU J110926.4-650224 with those observed in known and candidate tMSPs and discuss physical scenarios that may account for its persistent and flaring radio emissions.

Gravitationally lensed Type Ia supernovae are an emerging probe with great potential for constraining dark energy, spatial curvature, and the Hubble constant. The multiple images and their time delayed and magnified fluxes may be unresolved, however, blended into a single lightcurve. We demonstrate methods without a fixed source template matching for extracting the individual images, determining whether there are one (no lensing) or two or four (lensed) images, and measuring the time delays between them that are valuable cosmological probes. We find 100% success for determining the number of images for time delays greater than $\sim10$ days.

We first stress that any spherically symmetric galactic model whose integrated mass profile $M (r) \to 0$ as $r \to 0$ is physically consistent close to the centre only provided that the circular velocity $v_c (r) \to 0$ and the gravitational field $g (r) \to 0$ as $r \to 0$. Next, we apply such a statement to a broad class of five-parameter spherical galactic models, which includes most of those used in astrophysics and cosmology, In particular, we discover that the Jaffe and Hernquist models can only be trusted for $r \gtrsim 0.2 \, R_e$, while the NFW model cannot describe the central region either of regular galaxy clusters or of pure dark matter halos, thereby failing to predict any central cusp. As a consequence, we show that the observed gamma-ray excess from the Galactic Centre cannot be explained in terms of a cusp-enhanced WIMP annihilation into photons.

A distant, massive planet in the outer solar system has recently been proposed to explain some observed features of extreme trans-Neptunian objects (TNOs). Here we use N-body simulations of the formation of the Kuiper belt and Oort cloud as well as a survey simulator to compare models of the solar system with and without a 9th planet to one another as well as to observations. The main mechanism for TNOs to be deposited into the distant ($a$ > au), detached ($q$ > au) region of the Kuiper Belt in the 8-planet model is Kozai-Lidov oscillation of objects in mean motion resonances (MMR) with Neptune. This effect does not deposit low-inclination ($i \lesssim$ 20{\deg}) objects into this region. However, we find that the 9th planet generates a group of distant, detached TNOs at low inclinations that are not present in the 8-planet model. This disparity between the 8-planet and the 9-planet models could provide a strong constraint on a possible planet 9 with further detections of TNOs in the distant, detached region of the Kuiper Belt.

On 2019 August 14 at 21:10:39 UTC, the LIGO/Virgo Collaboration (LVC) detected a possible neutron star-black hole merger (NSBH), the first ever identified. An extensive search for an optical counterpart of this event, designated GW190814, was undertaken using DECam on the CTIO Blanco 4-m telescope. Target of opportunity interrupts were issued on 8 separate nights to observe 11 candidates using the SOAR Goodman Spectrograph in order to assess whether any of these transients was likely to be an optical counterpart of the possible NSBH merger. Here, we describe the process of observing with the SOAR Goodman spectrograph, the analysis of our spectra, our spectroscopic typing methodology, and our resultant conclusion that none of the candidates corresponded to the black hole-neutron star merger but were all instead other transient events. Finally, we describe the lessons learned from this effort. Application of these lessons will be critical for a successful community spectroscopic follow-up program for LVC season 4 (O4) and beyond.

Obtaining constraints from the largest scales of a galaxy survey is challenging due to the survey mask allowing only partial measurement of large angular modes. This scatters information from the harmonic-space 2-point function away from the diagonal and introduces coupling between modes. In this paper, we derive a custom eigenbasis adapted to any particular survey geometry so that all information is retained on the diagonal. At the expense of a somewhat complex pixel- and selection-function-window, the result is a diagonal 2-point function with a simple shot noise, and a diagonal covariance matrix in the case of a Gaussian random field. We derive the basis on the surface of a sphere, and we use it to construct a 3D spherical Fourier-Bessel power spectrum estimator assuming a survey geometry that is separable in the angular and radial directions.

To directly detect exoplanets and protoplanetary disks, the development of high accuracy wavefront sensing and control (WFS&C) technologies is essential, especially for ground-based Extremely Large Telescopes (ELTs). The Subaru Coronagraphic Extreme Adaptive Optics (SCExAO) instrument is a high-contrast imaging platform to discover and characterize exoplanets and protoplanetary disks. It also serves as a testbed to validate and deploy new concepts or algorithms for high-contrast imaging approaches for ELTs, using the latest hardware and software technologies on an 8-meter class telescope. SCExAO is a multi-band instrument, using light from 600 to 2500 nm, and delivering a high Strehl ratio (>80% in median seeing in H-band) downstream of a low-order correction provided by the facility AO188. Science observations are performed with coronagraphs, an integral field spectrograph, or single aperture interferometers. The SCExAO project continuously reaches out to the community for development and upgrades. Existing operating testbeds such as the SCExAO are also unique opportunities to test and deploy the new technologies for future ELTs. We present and show a live demonstration of the SCExAO capabilities (Real-time predictive AO control, Focal plane WFS&C, etc) as a host testbed for the remote collaborators to test and deploy the new WFS&C concepts or algorithms. We also present several high-contrast imaging technologies that are under development or that have already been demonstrated on-sky.

We present a method to implement relativistic corrections to the evolution of dark matter structures in cosmological simulations of a {\Lambda}CDM universe. We take the non-linear correspondence between the Lagrangian (Newtonian) evolution of dark matter inhomogeneities and the synchronous-comoving (relativistic) matter density description, and use it to promote the relativistic constraint as the initial condition for numerical simulations of structure formation. In this method, the incorporation of Primordial non-Gaussianity (PNG) contributions as initial conditions is straightforward. We implement the relativistic, fNL and gNL contributions as initial conditions for the L-PICOLA code, and derive the power spectrum and bispectrum of the evolved matter field. We focus specifically in the case of largest values of non-Gaussianity allowed at 1{\sigma}by Planck (fNL=-4.2 and gNL=-7000). As a checkup, we show consistency with the one-loop perturbative prescription in the adequate scales. Our results confirm that both relativistic and PNG features are most prominent at very large scales and for squeezed triangulations, but there is also a small relativistic effect at scales beyond the perturbative regime in the bispectrum. We discuss future prospects to probe these two contributions in the bispectrum of the matter density distribution.

Ultra-hot Jupiters represent an exciting avenue for testing extreme physics and observing atmospheric circulation regimes not found in our solar system. Their high temperatures result in thermally ionized particles embedded in atmospheric winds interacting with the planet's interior magnetic field by generating current and experiencing bulk Lorentz force drag. Previous treatments of magnetic drag in 3D General Circulation Models (GCMs) of ultra-hot Jupiters have mostly been uniform drag timescales applied evenly throughout the planet, which neglects the strong spatial dependence of these magnetic effects. In this work, we apply our locally calculated active magnetic drag treatment in a GCM of the planet WASP-76b. We find the effects of this treatment to be most pronounced in the planet's upper atmosphere, where strong differences between the day and night side circulation are present. These circulation effects alter the resulting phase curves by reducing the hotspot offset and increasing the day-night flux contrast. We compare our models to Spitzer phase curves which imply a magnetic field of at least 3 G for the planet. We additionally contrast our results to uniform drag timescale models. This work highlights the need for more careful treatment of magnetic effects in atmospheric models of hot gas giants.

Cosmic voids provide a unique environment to study galaxy formation and evolution. In this paper, we analyse a set of hydrodynamic zoom-in simulations of voids, to analyse in detail their inner structures. These voids were identified in a cosmological simulation and classified according to their surrounding dynamics at very large scales: whether they are in expansion or contraction at their outskirts. We study how these environments and the dynamics of voids impact the baryonic processes inside haloes and their mechanisms of formation and evolution. We find an under-abundance of processed gas within the voids that can be associated with the lack of massive haloes. By studying the dynamical phase-space diagram of haloes and the halo-particle correlation function, we find that haloes inside of contracting voids are slightly affected by the presence of bigger structures, in comparison to haloes in the inner regions of expanding voids. Consistent signals are obtained both when using dark matter and gas particles. We show that the halo assembly depends on the void dynamical state: haloes in expanding voids assemble slowly than those in contracting voids and in the general universe. This difference in the assembly impacts the baryonic evolution of haloes. Overall the redshift range analysed, haloes in voids have less baryon content than haloes in the general universe and particularly at z = 0 less stellar content. Our results suggest that the large scale void environment modulate the baryonic process occurring inside haloes according to the void dynamical state.

Using quiet-time measurements of element oxygen within the energy range 7.3--237.9 MeV nuc$^{-1}$ from the ACE spacecraft at 1 au, we compare the energy spectra and intensities of anomalous and Galactic cosmic rays (ACRs and GCRs, respectively) during 1997--2020. Our analysis shows that the transition from ACR-dominated spectrum to GCR-dominated spectrum occurs at energies $\sim$15 to $\sim$35 MeV nuc$^{-1}$, and the transition energy $E_t$ is found to be well anticorrelated with varying solar activity. This is the first study of ACR-GCR transition energy dependence on the solar cycle variation. At energies below $E_t$, the index of the power-law ACR-dominated spectrum ($\gamma_1$) ranges from -2.0 to -0.5, whereas the GCR-dominated spectrum has a power-law index ($\gamma_2$) changing from 0.3 to 0.8 at energies ranging from $E_t$ to 237.9 MeV nuc$^{-1}$. Both $\gamma_1$ and $\gamma_2$ are positively correlated with solar activity. In addition, during the solar cycle 24/25 minimum period, the peak GCR intensity observed by ACE spacecraft is about 8\% above its 2009 value, setting a new record since the space age, while the peak ACR intensity is almost similar to that of the previous two solar cycles with the same pattern of solar magnetic polarity, indicating a different modulation mechanism between ACRs and GCRs.

The magnetosphere of, and electromagnetic (EM) radiation from pulsars are usually described in the framework of classical electrodynamics. For some pulsars, however, whose emission heights are relatively close to the surface of the neutron star, general relativistic effects might modify the emission from the pulsar. We consider a magnetic dipole model of a pulsar to investigate general relativistic effects on EM radiation from it. Our study includes general relativistic modifications applicable to some significant issues in pulsar astronomy, such as the magnetosphere structure and pulse profiles. We implement computation of the magnetic field in the pulsar magnetosphere from a solution to Maxwell's equations defined in the strongly curved spacetime around a pulsar and find that the field exhibits a strong gravitational effect. The effect modifies curvature radiation of a pulsar, which then leads to modifications of the pulse profiles of radio emission. We take the pulsar PSR J1828-1101 as an example and work out Stokes parameters to simulate the pulse profiles for its main and interpulse emissions theoretically, which exhibit the gravitational effects clearly; however, their testability is beyond the current detection capabilities, with the absolute magnitude of the pulse profiles not being precisely predictable.

We present a measurement of the cosmic-ray spectrum above 100\,PeV using the part of the surface detector of the Pierre Auger Observatory that has a spacing of 750~m. An inflection of the spectrum is observed, confirming the presence of the so-called \emph{second-knee} feature. The spectrum is then combined with that of the 1500\,m array to produce a single measurement of the flux, linking this spectral feature with the three additional breaks at the highest energies. The combined spectrum, with an energy scale set calorimetrically via fluorescence telescopes and using a single detector type, results in the most statistically and systematically precise measurement of spectral breaks yet obtained. These measurements are critical for furthering our understanding of the highest energy cosmic rays.

We calculate the stellar evolution of both white dwarfs (WDs) in AM CVn binaries with orbital periods of $P_{\mathrm{orb}} \approx 5-70$ minutes. We focus on the cases where the donor starts as a $M_{\mathrm{He}} < 0.2 \, M_{\odot}$ Helium WD and the accretor is a $M_{\mathrm{WD}} > 0.6 \, M_{\odot}$ WD. Using Modules for Experiments in Stellar Astrophysics (MESA), we simultaneously evolve both WDs assuming conservative mass transfer and angular momentum loss from gravitational radiation. This self-consistent evolution yields the important feedback of the properties of the donor on the mass transfer rate, $\dot{M}$, as well as the thermal evolution of the accreting WD. Consistent with earlier work, we find that the high $\dot{M}$'s at early times forces an adiabatic evolution of the donor for $P_{\mathrm{orb}} < 30$ minutes so that its mass-radius relation depends primarily on its initial entropy. As the donor reaches $ M_{\mathrm{He}} \approx 0.02-0.03 \, M_{\odot}$ at $P_{\mathrm{orb}} \simeq 30 $ minutes, it becomes fully convective and could lose entropy and expand much less than expected under further mass loss. However, we show that the lack of reliable opacities for the donor's surface inhibit a secure prediction for this possible cooling. Our calculations capture the core heating that occurs during the first $\approx 10^7$ years of accretion and continue the evolution into the phase of WD cooling that follows. When compared to existing data for accreting WDs, as seen by Cheng and collaborators for isolated WDs, we also find that the accreting WDs are not as cool as we would expect given the amount of time they have had to cool.

This paper explores the potential of AstroPix, a project to develop Complementary Metal Oxide Semiconductor (CMOS) pixels for the next generation of space-based high-energy astrophysics experiments. Multimessenger astrophysics is a rapidly developing field whose upcoming missions need support from new detector technology such as AstroPix. ATLASPix, a monolithic silicon detector optimized for the ATLAS particle detector at CERN, is the beginning of the larger AstroPix project. Energy resolution is a driving parameter in the reconstruction of gamma-ray events, and therefore the characterization of ATLASPix energy resolution is the focus of this paper. The intrinsic energy resolution of the detector exceeded our baseline requirements of <10% at 60 keV. The digital output of ATLASPix results in energy resolutions insufficient to advance gamma-ray astronomy. However, the results from the intrinsic energy resolution indicate the digital capability of the detector can be redesigned, and the next generation of pixels for the larger AstroPix project have already been constructed. Iterations of AstroPix-type pixels are an exciting technology candidate to support new space-based missions.

Luminous blue variable stars (LBVs) are of great interest in massive-star evolution as they experience very high mass-loss episodes within short periods of time. HR Car is a famous member of this class in the Galaxy. It has a large circumstellar nebula and has also been confirmed as being in a binary system. One means of gaining information about the evolutionary status and physical nature of LBVs is studying their environments. We investigated the stellar content within ~100 pc of HR Car and also its circumstellar nebula. Very Large Telescope (VLT) Multi Unit Spectroscopic Explorer (MUSE) observations of a 2'x2' region around the star highlight the incompleteness of stellar classification for stars with magnitudes of V > 13 mag. Eight B0 to B9 stars have been identified which may lie in close spatial vicinity to HR Car. For a region with a radius of r =1.2 degree (~100 pc at a distance of 4.8 kpc) around HR Car, existing catalogs list several late O-type and early B-type stars, but only one early O-type star. Given the relatively low stellar and nebular masses in the HR Car system, no early O-type stars and only a few late O-type stars would be expected in association with HR Car. Instead, HR Car's location in a point vector diagram suggests that HR Car is not isolated, but is part of a moving group with a population of B-type stars in a spiral arm, and it has not received a strong kick from a supernova explosion of a companion star or a merger event. Potential binary evolution pathways for the HR Car system cannot be fully explored because of the unknown nature of the companion star. Furthermore, the MUSE observations reveal the presence of a fast outflow and "bullets" that have been ejected at intervals of about 400 years. These features may have been caused by recurrent mass transfer in the system.

Improving the performance of solar flare forecasting is a hot topic in solar physics research field. Deep learning has been considered a promising approach to perform solar flare forecasting in recent years. We first used the Generative Adversarial Networks (GAN) technique augmenting sample data to balance samples with different flare classes. We then proposed a hybrid convolutional neural network (CNN) model M for forecasting flare eruption in a solar cycle. Based on this model, we further investigated the effects of the rising and declining phases for flare forecasting. Two CNN models, i.e., Mrp and Mdp, were presented to forecast solar flare eruptions in the rising phase and declining phase of solar cycle 24, respectively. A series of testing results proved: 1) Sample balance is critical for the stability of the CNN model. The augmented data generated by GAN effectively improved the stability of the forecast model. 2) For C-class, M-class, and X-class flare forecasting using Solar Dynamics Observatory (SDO) line-of-sight (LOS) magnetograms, the means of true skill statistics (TSS) score of M are 0.646, 0.653 and 0.762, which improved by 20.1%, 22.3%, 38.0% compared with previous studies. 3) It is valuable to separately model the flare forecasts in the rising and declining phases of a solar cycle. Compared with model M, the means of TSS score for No-flare, C-class, M-class, X-class flare forecasting of the Mrp improved by 5.9%, 9.4%, 17.9% and 13.1%, and the Mdp improved by 1.5%, 2.6%, 11.5% and 12.2%.

In Mbarek & Caprioli (2019), we laid the groundwork for studying the espresso paradigm Caprioli (2015), a reacceleration mechanism to boost galactic CRs to UHECR levels. Our bottom-up approach uses realistic 3D MHD simulations of relativistic AGN jets and accounts for all of the crucial ingredients of a universal acceleration theory: injection, acceleration, and escape in realistic environments. Our results are consistent with the main features of UHECR spectra, i.e., power-law slopes, chemical composition, and anisotropy. In Mbarek & Caprioli (2021), we refine our model by including sub-grid particle scattering to model small-scale magnetic turbulence that cannot be resolved by MHD simulations, constraining for the first time one crucial but hard-to-model ingredient, and allowing us to establish the relative importance of espresso and stochastic shear acceleration in relativistic jets. Our framework also enables us to analyze high-energy neutrinos produced from our accelerated UHECRs considering the effects of external photon fields, and to incorporate nucleus photodisintegration. The spectra we obtain are consistent with the picture drawn by observations with Auger, Telescope Array, and IceCube observatory.

Fast Radio Bursts (FRBs) are a promising tool for studying the low-density universe as their dispersion measures (DM) are extremely sensitive probes of electron column density. Active Galactic Nuclei (AGN) inject energy into the intergalactic medium, affecting the DM and its scatter. To determine the effectiveness of FRBs as a probe of AGN feedback, we analysed three different AGN models from the EAGLE simulation series. We measured the mean DM-redshift relation, and the scatter around it, using $2.56 \times 10^8$ sightlines at 131 redshift ($z$) bins between $0 \leq z \leq 3$. While the DM-redshift relation itself is highly robust against different AGN feedback models, significant differences are detected in the scatter around the mean: weaker feedback leads to more scatter. We find $\sim 10^4$ localised FRBs are needed to discriminate between the scatter in standard feedback and stronger, more intermittent feedback models. The number of FRBs required is dependent on the redshift distribution of the detected population. A log-normal redshift distribution at $z=0.5$ requires approximately 50% fewer localised FRBs than a distribution centred at $z=1$. With the Square Kilometre Array expected to detect $>10^3$ FRBs per day, in the future, FRBs will be able to provide constraints on AGN feedback.

GRBs are the most energetic electromagnetic events in the Universe. Those whose typical duration is longer than a few seconds are known as long GRBs and shorter than a few seconds are short GRBs. It is widely accepted that these events are associated with the collapse of a very massive star and the neutron star (NS) binary merger, respectively. A fast-spinning, strongly magnetized NS could be expected before a black hole (BH) in both scenarios. We allude to the thermal neutrinos' particular properties propagating inside the fireball for differentiating both scenarios in this work. We first derive the neutrino effective potential associated with each medium in a strong and weak magnetic field. We calculate the three-flavor oscillation probabilities, and finally, we get the expected neutrino rate in both scenarios. Given these observables' evolution, we can determine whether the progenitor could be associated with a strongly magnetized NS or a BH.

Hayabusa2 collected 5.4 g of samples from the asteroid Ryugu and brought them back to Earth. Measuring these samples using a camera system comparable to that used for remote-sensing observations is important for characterizing the collected samples, examining sample representativeness, and identifying materials collected from the asteroid. In this study, we constructed an instrument that enables both visual multispectral imaging and 3D shape reconstruction of samples based on stereo imaging. The imaging system has the pixel resolution of 1.93 {\mu}m/pix and field of view of 7.9 mm x 4.2 mm. Our validation measurements demonstrate that our multispectral imaging system, comparable to the telescopic optical navigation camera (ONC-T) on Hayabusa2, yields reflectance spectra with an accuracy of 3% and a 3D model with a precision of 5%. Using this instrument, we conducted multispectral measurements of two Ryugu samples acquired from two locations on the asteroid. The average spectra of the samples are consistently flat with v-band reflectance of 2.4%, higher than the overall average of Ryugu as observed with ONC-T. This apparent difference could be attributed to the potentially reflective surface of the returned grain samples.

The Kepler transit survey with follow-up spectroscopic observations has discovered numerous super-Earth sized planets and revealed intriguing features of their sizes, orbital periods, and their relations between adjacent planets. For the first time, we investigate the size evolution of planets via both giant impacts and photoevaporation to compare with these observed features. We calculate the size of a protoplanet, which is the sum of its core and envelope sizes, by analytical models. $N$-body simulations are performed to evolve planet sizes during the giant impact phase with envelope stripping via impact shocks. We consider the initial radial profile of the core mass and the initial envelope mass fractions as parameters. Inner planets can lose their whole envelopes via giant impacts, while outer planets can keep their initial envelopes since they do not experience giant impacts. Photoevaporation is simulated to evolve planet sizes afterward. Our results suggest that the period-radius distribution of the observed planets would be reproduced if we perform simulations in which the initial radial profile of the core mass follows a wide range of power-law distributions and the initial envelope mass fractions are $\sim0.1$. Moreover, our model shows that the adjacent planetary pairs have similar sizes and regular spacings, with slight differences from detailed observational results such as the radius gap.

Andes Large-area PArticle detector for Cosmic-ray physics and Astronomy (ALPACA) is an international experiment that applies southern very-high-energy (VHE) gamma-ray astronomy to determine the origin of cosmic rays around the knee energy region ($10^{15}\, {\rm eV} - 10^{16}\, {\rm eV}$). The experiment consists of an air shower (AS) array with a surface of $83, 000\, {\rm m}^2$ and an underground water Cherenkov muon detector (MD) array covering $5, 400\, {\rm m}^2$. The experimental site is at the Mt. Chacaltaya plateau in La Paz, Bolivia, with an altitude of $4, 740\, {\rm m}$ corresponding to $572\, {\rm g}/{\rm cm}^2$ atmospheric thickness. As the prototype experiment of ALPACA, the ALPAQUITA experiment aims to begin data acquisition in late 2021. The ALPAQUITA array consists of a smaller AS array ($18, 450\, {\rm m}^2$) and underground MD ($900\, {\rm m}^2$), which are now under construction. ALPAQUITA's sensitivity to gamma-ray sources is evaluated with Monte Carlo simulations. The simulation finds that five gamma-ray sources observed by H.E.S.S. and HAWC experiments will be detected by ALPAQUITA beyond 10 TeV and one out of these five - HESS J1702-420A - above 300 TeV in one calendar year observation. The latter finding means that scientific discussions can be made on the emission mechanism of gamma rays beyond 100 TeV from southern sources on the basis of the observational results of this prototype experiment.

We present high angular resolution ALMA CO(2-1) and 1.7 mm continuum observations of the far-infrared-bright galaxy PKS0023-26, which hosts a young radio source and a luminous optical AGN. Although young, the powerful radio source has grown to a size of a few kpc, making it capable of affecting the ISM of the host galaxy. We detect an extended distribution of molecular gas with a mass between 0.3 and 3x10^10 Msun, depending on the CO conversion factor. The gas has a maximum extent of ~24 kpc and is distributed asymmetrically wrt the radio galaxy. Overall, the observed properties are reminiscent of molecular structures observed in some galaxy clusters. However, in the inner few kpc, the kinematics of the gas appears to be affected by the radio source. In the central regions we observe the brightest emission and the broadest profiles (FWZI ~ 500 km/s), which indicate a direct interaction of the jet with dense clouds. On larger scales, the molecular gas appears to avoid the radio lobes and gas with smaller velocity dispersion (FWZI ~350 km/s) is observed around the them. There, the gas appears to be affected by the expanding cocoon surrounding the radio source, dispersing and heating preexisting molecular clouds. The data suggest that the mode of coupling between radio jets and ISM changes from an outflowing phase in the inner regions to a maintenance phase at larger radii. This reveals that on galaxy scales the impact of the AGN is not limited to outflows. With a star-formation rate of 25 Msun/yr, PKS0023-26 is located on the SFR-M* relation for star forming galaxies. The AGN does not appear to have, at present, a major impact on the overall level of star formation of the host galaxy. However, as the jet and lobes expand throughout the galaxy, they will carry enough energy to prevent further gas cooling and/or to inject turbulence and affect future star formation.

Primordial magnetic fields could explain the large-scale magnetic fields present in the Universe. Inflation and phase transitions in the early Universe could give rise to such fields with unique characteristics. We investigate the magneto-hydrodynamic evolution of these magnetogenesis scenarios with cosmological simulations. We evolve inflation-generated magnetic fields either as (i) uniform (homogeneous) or as (ii) scale-invariant stochastic fields, and phase transition-generated ones either as (iii) helical or as (iv) non-helical fields from the radiation-dominated epoch. We find that the final distribution of magnetic fields in the simulated cosmic web shows a dependence on the initial strength and the topology of the seed field. Thus, the observed field configuration retains information on the initial conditions at the moment of the field generation. If detected, primordial magnetic field observations would open a new window for indirect probes of the early universe. The differences between the competing models are revealed on the scale of galaxy clusters, bridges, as well as filaments and voids. The distinctive spectral evolution of different seed fields produces imprints on the correlation length today. We discuss how the differences between rotation measures from highly ionized regions can potentially be probed with forthcoming surveys.

The existence of flat or bulgeless galaxies poses a challenge to the hierarchical structure formation scenario advocated by modern cosmology. We determine the geometrical environment of a sample of $315$ flat galaxies from the Revised Flat Galaxy Catalog (RFGC) using `local dimension' $D$, which, on a given length scale, quantifies the dimension of the cosmic structure in which a galaxy is embedded. For galaxies residing in filaments, nodes and sheets, $D \sim 1$, $D \sim 1.5$ and $D \sim 2$ respectively; $D \sim 3$ represents field galaxies. We also determine the local dimensions of a sample of 15,622 non-flat galaxies identified in the Galaxy Zoo project from the Sloan Digital Sky Survey (SDSS). We find that the median values of $D$ for the flat and the non-flat galaxies are $2.2$ and $1.8$ respectively, implying that flat galaxies are located in a relatively sparser environment compared to non-flat galaxies; a Kolmogorov-Smirnov (KS) test indicates that their geometrical environments are different at $>99$\% confidence level. Further, using a group finding algorithm, we study the local environment of a subset of $779$ flat galaxies with major-to-minor axes ratio $a/b >10$ identified as superthin galaxies. We find that the median clusterization index $k_{\rm{min}}$ for superthin flat galaxies $\sim$ $2.3$ while $\sim$ $1.7$ for other flat galaxies, confirming that the superthins reside in an under-dense environment compared to other flat galaxies at $> 98 \%$ confidence level. Our results may therefore have important implications for the formation and evolution models of flat galaxies in the universe.

High-energy solar flares have shown to have at least two distinct phases: prompt-impulsive and delayed-gradual. Identifying the mechanism responsible for accelerating the electrons and ions and the site at which it occurs during these two phases is one of the outstanding questions in solar physics. Many advances have been made over the past decade thanks to new observational data and refined simulations that together help to shed light on this topic. For example, the detection by Fermi Large Area Telescope (LAT) of GeV emission from solar flares originating from behind the visible solar limb and >100 MeV emission lasting for more than 20 hours have suggested the need for a spatially extended source of acceleration during the delayed emission phase. In this work we will review some of the major results from Fermi-LAT observations of the 24th solar cycle and how this new observational channel combined with observations from across the electromagnetic spectrum can provide a unique opportunity to diagnose the mechanisms of high-energy emission and particle acceleration in solar flares.

We present the results of simultaneous 450 $\mu$m and 850 $\mu$m polarization observations toward the massive star forming region NGC 2071IR, a target of the BISTRO (B-fields in Star-Forming Region Observations) Survey, using the POL-2 polarimeter and SCUBA-2 camera mounted on the James Clerk Maxwell Telescope. We find a pinched magnetic field morphology in the central dense core region, which could be due to a rotating toroidal disk-like structure and a bipolar outflow originating from the central young stellar object, IRS 3. Using the modified Davis-Chandrasekhar-Fermi method, we obtain a plane-of-sky magnetic field strength of 563$\pm$421 $\mu$G in the central $\sim$0.12 pc region from 850 $\mu$m polarization data. The corresponding magnetic energy density of 2.04$\times$10$^{-8}$ erg cm$^{-3}$ is comparable to the turbulent and gravitational energy densities in the region. We find that the magnetic field direction is very well aligned with the whole of the IRS 3 bipolar outflow structure. We find that the median value of polarization fractions, 3.0 \%, at 450 $\mu$m in the central 3 arcminute region, which is larger than the median value of 1.2 \% at 850 $\mu$m. The trend could be due to the better alignment of warmer dust in the strong radiation environment. We also find that polarization fractions decrease with intensity at both wavelengths, with slopes, determined by fitting a Rician noise model, of $0.59 \pm 0.03$ at 450 $\mu$m and $0.36 \pm 0.04$ at 850 $\mu$m, respectively. We think that the shallow slope at 850 $\mu$m is due to grain alignment at the center being assisted by strong radiation from the central young stellar objects.

The emission of very-high-energy photons (VHE, E$>$100 GeV) in active galactic nuclei (AGN) is closely connected with the production of ultra-relativistic particles. Among AGN, the subclass of extreme BL Lacertae are of particular interest because they challenge state-of-art models on how these cosmic particle accelerators operate. By cross-matching two gamma-ray catalogs (this is, 4FGL-DR2 and 2BIGB), we identified 23 high-synchrotron-peaked (HSP) blazar candidates with photometric or spectroscopic redshifts, good multi-wavelength coverage, that are possibly detectable by VHE instruments. We performed a new analysis of Fermi Large Area Telescope data including the effects of attenuation from the extragalactic background light and complemented these results by collecting multiwavelength data from optical, radio and X-ray archival observations. Their broadband spectral energy distributions were interpreted in terms of synchrotron-self-Compton models with external-Compton components and compared with the properties of prototypical extreme HSP blazars. Finally, we test their detectability with imaging atmospheric Cherenkov telescopes (IACTs) and propose a new method for selecting these extreme targets for these ground-based telescopes.

ScopeSim is a flexible multipurpose instrument data simulation framework built in Python. It enables both raw and reduced observation data to be simulated for a wide range of telescopes and instruments quickly and efficiently on a personal computer. The software is currently being used to generate simulated raw input data for developing the data reduction pipelines for the MICADO and METIS instruments at the ELT. The ScopeSim environment consists of three main packages which are responsible for providing on-sky target templates (ScopeSim_templates), the data to build the optical models of various telescopes and instruments (instrument reference database), and the simulation engine (ScopeSim). This strict division of responsibilities allows ScopeSim to be used to simulate observation data for many different instrument and telescope configurations for both imaging and spectroscopic instruments. ScopeSim has been built to avoid redundant calculations wherever possible. As such it is able to deliver simulated observations on time scales of seconds to minutes. All the code and data is open source and hosted on Github. The community is also most welcome, and indeed encouraged to contribute to code ideas, target templates, and instrument packages.

The spectral evolution of transient X-ray binaries can be reproduced by an interplay between two flows separated at a radius $R_J$: a standard accretion disk (SAD) in the outer parts and a jet-emitting disk (JED) in the inner parts. In the previous papers of this series, we recover the spectral evolution in both X-rays and radio for four outbursts of GX339-4 by playing independently with the two parameters: $R_J$ and the accretion rate $\dot{M}_{in}$. In this paper, we compare the time evolution of $R_J$ and $\dot{M}_{in}$ for the four outbursts. We show that despite the undeniable differences between the time evolution of each outburst, a unique pattern in the $\dot{M}_{in}-R_J$ plane seems to be followed by all cycles within the JED-SAD model. We call this pattern a fingerprint, and show that even the 'failed' outburst considered follows it. We also compute the radiative efficiency in X-rays during the cycles and consider its impact on the radio--X-ray correlation. Within the JED-SAD paradigm, we find that the accretion flow is always radiatively efficient in the hard states, with between $15\%$ and $40\%$ of the accretion power being radiated away at any given time. Moreover, we show that the radiative efficiency evolves with the accretion rate because of key changes in the JED thermal structure. These changes give birth to two different regimes with different radiative efficiencies: the thick disk and the slim disk. While the existence of these two regimes is intrinsically linked to the JED-SAD model, we show direct observational evidence of the presence of two different regimes using the evolution of the X-ray power-law spectral index, a model-independent estimate. We then argue that these two regimes could be the origin of the gap in X-ray luminosity in the hard state, the wiggles and different slopes seen in the radio--X-ray correlation, and even the existence of outliers.

Galaxies living in rich environments are suffering different perturbations able to drastically affect their evolution. Among these, ram pressure stripping, i.e. the pressure exerted by the hot and dense intracluster medium (ICM) on galaxies moving at high velocity within the cluster gravitational potential well, is a key process able to remove their interstellar medium (ISM) and quench their activity of star formation. This review is aimed at describing this physical mechanism in different environments, from rich clusters of galaxies to loose and compact groups. We summarise the effects of this perturbing process on the baryonic components of galaxies, from the different gas phases (cold atomic and molecular, ionised, hot) to magnetic fields and cosmic rays, and describe their induced effects on the different stellar populations, with a particular attention to its role in the quenching episode generally observed in high density environments. We also discuss on the possible fate of the stripped material once removed from the perturbed galaxies and mixed with the ICM, and we try to estimate its contribution to the pollution of the surrounding environment. Finally, combining the results of local and high redshift observations with the prediction of tuned models and simulations, we try to quantify the importance of this process on the evolution of galaxies of different mass, from dwarfs to giants, in various environments and at different epochs.

We analyse properties of the mass outflow from the Roche-lobe filling component of a semi-detached binary system. We follow the approaches published by Paczy\'nski \& Sienkiewicz and by Lubow \& Shu, which we compare with other simplified approaches. We find that the density of the flow at $L_1$ is orders of magnitude lower than the density on the same equipotential but away from $L_1$. Furthermore, the effective cross section of the flow, after averaging over its profile of the momentum density, is much lower than some published estimates done without accounting for the averaging. Thus, the use of some simplified formulae for the density and the flow cross section can lead to overestimates of the accretion rate and of the mass contained in the $L_1$ regions by very large factors unless they are supported by simultaneous integrations of the equations of stellar structure for the outer layers of the donor.

We apply a turbulence-regulated model of star formation to calculate the star formation rate (SFR) of dense star-forming clouds in simulations of jet-ISM interactions. The method isolates individual clumps and accounts for the impact of virial parameter and Mach number of the clumps on the star formation activity. This improves upon other estimates of the SFR in simulations of jet--ISM interactions, which are often solely based on local gas density, neglecting the impact of turbulence. We apply this framework to the results of a suite of jet-ISM interaction simulations to study how the jet regulates the SFR both globally and on the scale of individual star-forming clouds. We find that the jet strongly affects the multi-phase ISM in the galaxy, inducing turbulence and increasing the velocity dispersion within the clouds. This causes a global reduction in the SFR compared to a simulation without a jet. The shocks driven into clouds by the jet also compress the gas to higher densities, resulting in local enhancements of the SFR. However, the velocity dispersion in such clouds is also comparably high, which results in a lower SFR than would be observed in galaxies with similar gas mass surface densities and without powerful radio jets. We thus show that both local negative and positive jet feedback can occur in a single system during a single jet event, and that the star-formation rate in the ISM varies in a complicated manner that depends on the strength of the jet-ISM coupling and the jet break-out time-scale.

As a novel approach for tracing interstellar magnetic fields, the Velocity Gradient Technique (VGT) has been proven to be effective for probing magnetic fields in the diffuse interstellar medium (ISM). In this work, we verify the VGT in a broader context by applying the technique to a molecular cloud interacting with the supernovae remnant (SNR) W44. We probe the magnetic fields with the VGT using CO, $\rm HCO^+$, and H I emission lines and make a comparison with the Planck 353 GHZ dust polarization. We show that the VGT gives an accurate measurement that coheres with the Planck polarization especially in intense molecular gas emission regions. We further study the foreground's contribution to the polarization that results in a misalignment between the VGT and the Planck measurements in low-intensity molecular gas areas. We advance the VGT to achieve magnetic field tomography by decomposing the W44 into various velocity components. We show that W44's velocity component at $v\sim45$ km s$^{-1}$ exhibits the largest coverage and gives the best agreement with Planck polarization in terms of magnetic field orientation.

We use multi-spacecraft observations of invididual type III radio bursts in order to calculate the directivity of the radio emission, to be compared to the results of ray-tracing simulations of the radio-wave propagation and probe the plasma properties of the inner heliosphere. Ray-tracing simulations of radio-wave propagation with anisotropic scattering on density inhomogeneities are used to study the directivity of radio emissions. Simultaneous observations of type III radio bursts by four widely-separated spacecraft are used to calculate the directivity and position of the radio sources. The shape of the directivity pattern deduced for individual events is compared to the directivity pattern resulting from the ray-tracing simulations. We show that simultaneous observations of type radio III bursts by 4 different probes provide the opportunity to estimate the radio source positions and the directivity of the radio emission. The shape of the directivity varies from one event to another, and is consistent with anisotropic scattering of the radio-waves.

We present an all-sky 90\% confidence level upper limit on the cosmic flux of relativistic magnetic monopoles using 2886 days of IceCube data. The analysis was optimized for monopole speeds between 0.750$c$ and 0.995$c$, without any explicit restriction on the monopole mass. We constrain the flux of relativistic cosmic magnetic monopoles to a level below $2.0\times 10^{-19} {\textrm{cm}}^{-2} {\textrm{s}}^{-1} {\textrm{sr}}^{-1}$ over the majority of the targeted speed range. This result constitutes the most strict upper limit to date for magnetic monopoles above the Cherenkov threshold and up to $\beta \sim 0.995$ and fills the gap between existing limits on the cosmic flux of non-relativistic and ultrarelativistic magnetic monopoles

This conference proceedings is a write-up of the Gamma-ray Indirect rapporteur talk given at the 37th International Cosmic Ray Conference (ICRC 2021). In contrast to previous ICRCs, this years edition was held in a fully virtual format, with dedicated discussion sessions organised around specific scientific themes. Many of these topics span the two categories of Gamma-ray Indirect (GAI) and Gamma-ray Direct (GAD), observations of gamma-rays by ground-based and space-based facilities respectively. To cover this organisation by topic in a coherent manner, this GAI rapporteur contribution focuses predominantly (but not exclusively) on Galactic gamma-ray astronomy, whereas the GAD rapporteur contribution focuses predominantly (but not exclusively) on Extra-galactic gamma-ray astronomy. In recent years, the field has seen enormous progress in both theory and observation, particularly in identifying PeVatrons (accelerators of Cosmic Rays to PeV energies), studies of particle escape from the accelerator, and detection of gamma-ray transients, especially gamma-ray bursts.

The Hydrogen Intensity and Real-time Analysis eXperiment (HIRAX) is a radio interferometer array currently in development, with an initial 256-element array to be deployed at the South African Radio Astronomy Observatory (SARAO) Square Kilometer Array (SKA) site in South Africa. Each of the 6~m, $f/0.23$ dishes will be instrumented with dual-polarisation feeds operating over a frequency range of 400-800 MHz. Through intensity mapping of the 21 cm emission line of neutral hydrogen, HIRAX will provide a cosmological survey of the distribution of large-scale structure over the redshift range of $0.775 < z < 2.55$ over $\sim$15,000 square degrees of the southern sky. The statistical power of such a survey is sufficient to produce $\sim$7 percent constraints on the dark energy equation of state parameter when combined with measurements from the Planck satellite. Additionally, HIRAX will provide a highly competitive platform for radio transient and HI absorber science while enabling a multitude of cross-correlation studies. In these proceedings, we describe the science goals of the experiment, overview of the design and status of the sub-components of the telescope system, and describe the expected performance of the initial 256-element array as well as the planned future expansion to the final, 1024-element array.

We present Northern Extended Millimeter Array (NOEMA) CO(2-1) maps of the z=0.4418 cluster-central QSO IRAS 09104+4109, which trace ~4.5x10^10 Msol of molecular gas in and around the galaxy. As in many low redshift cool core clusters, the molecular gas is located in a series of clumps extending along the old radio jets and lobes. It has a relatively low velocity dispersion (336 [+39,-35] km/s FWHM) and shows no velocity gradients indicative of outflow or infall. Roughly half the gas is located in a central clump on the northeast side of the galaxy, overlapping a bright ionized gas filament and a spur of excess X-ray emission, suggesting that this is a location of rapid cooling. The molecular gas is unusually extended, out to ~55 kpc radius, comparable to the scale of the filamentary nebula in the Perseus cluster, despite the much higher redshift of this system. The extent falls within the thermal instability radius of the intracluster medium (ICM), with t_cool/t_ff<25 and t_cool}/t_eddy~1 within ~70 kpc. Continuum measurements at 159.9 GHz from NOEMA and 850 micron from the JCMT SCUBA-2 show excess far infrared emission, which we interpret as free-free emission arising from the ongoing starburst. These observations suggest that ICM cooling is not strongly affected by the buried QSO, and that cooling from the ICM can build gas reservoirs sufficient to fuel quasar-mode activity and drive the reorientation of the central AGN.

Stephan's Quintet (SQ) is a nearby compact galaxy group and a perfect laboratory for studying the process of galaxy evolution through galaxy harassment and interaction. By analysing the kinematics of SQ we aim to provide an increased understanding of the group, the history of the interactions, their cause and effect, and the details regarding the physical processes occurring as galaxies interact. We have studied the ionised gas and stellar kinematics using the Large Binocular Telescope, and the molecular gas kinematics via CO using the IRAM 30m. Large areas of the group have been mapped and analysed. We obtain a total ionised gas mass in the regions chosen for closer analysis of 20.1$\pm$0.2x10^10 Msun and a total H2 gas mass of 21$\pm$2x10^9 Msun in the observed area (spectra integrated over the velocity range of SQ), while the star-forming (SF) clouds show an impressive complexity, with gas congregations at multiple velocities throughout the group. We map the large-scale nuclear wind in NGC7319 and its decoupled gas and stellar disk. With our high resolution data we can, for the first time, reveal the Seyfert 1 nature of NGC7319 and fit the narrow-line and broad-line regions. While the 12CO(1-0) map shows significant emission in the area of NGC7319, the bridge, and the SF ridge, the 12CO(2-1) emission shows a prevalence to the SF ridge, an area south of the NGC7318 pair, and an extension towards NGC7317 - connecting NGC7317 to the centre of the group, indicating a previous interaction. NGC7317 may also be a prime candidate for studies of the process of galaxy harassment. Furthermore, we connect the kinematical structures in SQ to the history of the group and the ongoing interaction with NGC7318B. Through our extensive observations of SQ we trace the kinematics and evolution of the complex processes and structures occurring in this nearby interactive group. [Abstract abridged]

We report here the first results from a 15-yr long variability monitoring of the z=2.0 quasar QSO B1312+7837. It shows luminosity changes with a period P~6.13 yr (P~2.04 yr at rest frame) and an amplitude of ~0.2mag, superimposed on a gradual dimming at a rate of ~0.55mag per 100 yrs. Two false periods associated with power peaks in the data windowing function were discarded. The measured period is confirmed with a bootstrapping Monte-Carlo simulation. A damped random walk model yields a better fit to the data than a sine-function model, but at the cost of employing some high frequency variations which are typically not seen in quasars. We consider the possible mechanisms driving this variability, and conclude that orbital motion of two supermassive black holes - result from a recent galaxy merger - is a possible explanation.

Interferometric observations of the Sun with ALMA provide valuable diagnostic tools for studying the small-scale dynamics of the solar atmosphere. Estimations of the observability of the small-scale dynamics as a function of spatial resolution for regions with different characteristic magnetic field topology facilitate a more robust analysis of ALMA observations of the Sun. A 3D model of the solar atmosphere from the MHD code Bifrost was used to produce high-cadence observables at mm and submm wavelengths. The synthetic observables for receiver bands 3-10 were degraded to the angular resolution corresponding to ALMA observations with different configurations of the interferometric array from the most compact C1, to the more extended C7. The observability of the small-scale dynamics was analyzed in each case. The analysis was thus also performed for predicting the potential of future capabilities. The minimum resolution required to study the typical small spatial scales in the solar chromosphere depends on the characteristic properties of the target region. Here, a range from quiet Sun to enhanced network loops is considered. Limited spatial resolution affects the observable signatures of dynamic small-scale brightening events in the form of reduced brightness temperature amplitudes, potentially leaving them undetectable, and even shifts in the times at which the peaks occur of up to tens of seconds. Conversion factors between the observable brightness amplitude and the original amplitude in the fully resolved simulation are provided that can be applied to observational data in principle, but are subject to wavelength-dependent uncertainties. Predictions of the typical appearance at the different combinations of receiver band, array configuration, and properties of the target region are conducted.

Short Gamma-Ray Bursts (SGRBs) are produced by the coalescence of compact binary systems which are remnants of massive stars. GRB 160410A is classified as a short-duration GRB with extended emission and is currently the farthest SGRB with a redshift determined from an afterglow spectrum and also one of the brightest SGRBs to date. The fast reaction to the Neil Gehrels Swift Observatory alert allowed us to obtain a spectrum of the afterglow using the X-shooter spectrograph at the Very Large Telescope (VLT). The spectrum shows a number of absorption features at a redshift of z=1.7177, in addition, we detect two intervening systems at z=1.581 and z=1.444. The spectrum shows ly-alpha in absorption with a column density of log N(HI)=21.3+/-0.3 cm$^{-2}$ which, together with FeII, CII, SiII and OI, allow us to perform the first study of chemical abundances in a SGRB host galaxy. We determine a metallicity of [Fe/H]=-2.7+/-0.3, significantly lower than observed for any long GRB host. We find no evidence for extinction in the afterglow spectral energy distribution (SED). The environment has a low degree of ionisation and the CIV and SiIV lines are completely absent. We do not detect an underlying host galaxy down to deep limits. Additionally, we present the spectrum of GRB 201221D, another high-z short GRB that shows absorption lines at z=1.045 but whose environment seems to be more similar to the one of short GRBs as derived from the SED fitting to the host photometry.

We present the $NuSTAR$ extragalactic survey of the $James Webb Space Telescope$ ($JWST$) North Ecliptic Pole (NEP) Time-Domain Field. The survey covers a $\sim$0.16 deg$^2$ area with a total exposure of 681 ks acquired in a total of nine observations from three epochs. The survey sensitivities at 20% of the area are 2.39, 1.14, 2.76, 1.52, and 5.20 $\times$ 10$^{-14}$ erg cm$^{-2}$ s$^{-1}$ in the 3-24, 3-8, 8-24, 8-16, and 16-24 keV bands, respectively. The NEP survey is one of the most sensitive extragalactic surveys with $NuSTAR$ so far. A total of 33 sources were detected above 95% reliability in at least one of the five bands. We present the number counts, log$N$-log$S$, measured in the hard X-ray 8-24 and 8-16 keV bands, uniquely accessible by $NuSTAR$ down to such faint fluxes. We performed source detection on the XMM-$Newton$ and $Chandra$ observations of the same field to search for soft X-ray counterparts of each $NuSTAR$ detection. The soft band positions were used to identify optical and infrared associations. We present the X-ray properties (hardness ratio and luminosity) and optical-to-X-ray properties of the detected sources. The measured fraction of candidate Compton-thick (N$\rm _H\ge10^{24} cm^{-2}$) active galactic nuclei, derived from the hardness ratio, is between 3% to 27%. As this survey was designed to have variability as its primary focus, we present preliminary results on multi-epoch flux variability in the 3-24 keV band.

There are now strong indications that white dwarf (WD) stars with masses well below the Chandrasekhar limit (MCh ~ 1.4 Msun) contribute a significant fraction of SN Ia progenitors. The relative fraction of stable iron-group elements synthesized in the explosion has been suggested as a possible discriminant between MCh and sub-MCh events. In particular, it is thought that the higher-density ejecta of MCh WDs, which favours the synthesis of stable isotopes of nickel, results in prominent [Ni II] lines in late-time spectra. We study the explosive nucleosynthesis of stable nickel in SNe Ia resulting from MCh and sub-MCh progenitors, and explore the potential for lines of [Ni II] at 7378 A and 1.94 microns in late-time spectra to serve as a diagnostic of the exploding WD mass, using nonlocal thermodynamic equilibrium radiative-transfer simulations with the CMFGEN code. We find that the radiative proton-capture reaction 57Co(p,gamma)58Ni is the dominant production mode for 58Ni in both MCh and sub-MCh models, while the alpha-capture reaction on 54Fe has a negligible impact on the final 58Ni yield. More importantly, we demonstrate that the lack of [Ni II] lines in late-time spectra of sub-MCh events is not always due to an under-abundance of stable Ni, but results from the higher ionization of Ni in the inner ejecta. Conversely, the strong [Ni II] lines predicted in our 1D MCh models are completely suppressed when 56Ni is sufficiently mixed with the innermost layers rich in stable iron-group elements. [Ni II] lines in late-time SN Ia spectra have a complex dependency on the abundance of stable Ni, which limits their use alone in distinguishing between MCh and sub-MCh progenitors. However, we argue that a low-luminosity SN Ia displaying strong [Ni II] lines would most likely result from a Chandrasekhar-mass progenitor. [Abridged]

To date, 240 individual molecular species, comprised of 19 different elements, have been detected in the interstellar and circumstellar medium by astronomical observations. These molecules range in size from two atoms to seventy, and have been detected across the electromagnetic spectrum from cm-wavelengths to the ultraviolet. This census presents a summary of the first detection of each molecular species, including the observational facility, wavelength range, transitions, and enabling laboratory spectroscopic work, as well as listing tentative and disputed detections. Tables of molecules detected in interstellar ices, external galaxies, protoplanetary disks, and exoplanetary atmospheres are provided. A number of visual representations of this aggregate data are presented and briefly discussed in context.

Observed Cosmic Microwave Background (CMB) maps are contaminated by foregrounds, some of which are usually masked to perform cosmological analyses. If masks are correlated to the lensing signal, such as those removing extragalactic emissions located in matter overdensities, measurements over the unmasked sky may give biased estimates. We quantify the impact of these mask-induced biases for the reconstructed CMB lensing auto- and cross-correlation power spectra with external matter tracers. We show that they arise both from changes in the lensing power, and via modifications to the reconstruction power spectrum corrections, $N_L^{(0)}$, $N_L^{(1)}$ and $N_L^{(3/2)}$). For direct masking of the CMB lensing field, we derive simple analytic models of the masking effect and show that it is potentially large. We show that mask-induced biases are significantly reduced by optimal filtering of the CMB maps in the lensing reconstruction. We test the resulting lensing power spectrum biases on numerical simulations, masking radio sources, and peaks of thermal Sunyaev-Zeldovich (tSZ) and cosmic infrared background (CIB) emission. For the lensing auto spectrum, masking biases can only be measured with a statistical significance $\lesssim 3\sigma$ for future data sets. The same applies to the cross-correlation power spectra between CMB lensing and tSZ and CIB even though biases are larger (up to ~30%). We find that masking tSZ-selected galaxy clusters leads to the largest mask biases, potentially detectable with high significance. We find that the calibration of cluster masses using CMB lensing, in particular for objects at $z\lesssim 0.6$, might be significantly affected by mask biases for near-future observations if the lensing signal recovered inside the mask holes is used without further corrections. Conversely, mass calibration of high redshift objects will still deliver unbiased results.

We report on the extinction properties in the fields around the clusters NGC 1854, NGC 1856, and NGC 1858 in the bar of the Large Magellanic Cloud. The colour-magnitude diagrams of the stars in all these regions show an elongated red giant clump that reveals a variable amount of extinction across these fields, ranging from Av~0.2 to Av~1.9, including Galactic foreground extinction. The extinction properties nonetheless are remarkably uniform. The slope of the reddening vectors measured in the (V-I,V) and (B-I,B) colour-magnitude planes is fully in line with the Av/E(B-V)~5.5 value found in the outskirts of 30 Dor. This indicates the presence of an additional grey extinction component in the optical requiring big grains to be about twice as abundant as in the diffuse Galactic interstellar medium (ISM). Areas of higher extinction appear to be systematically associated with regions of more intense star formation, as measured by the larger number of stars more massive than 8 Msun, thus making injection of big grains into the ISM by SNII explosion the likely mechanism at the origin of the observed grey extinction component.

The Forward Physics Facility (FPF) is a proposal to create a cavern with the space and infrastructure to support a suite of far-forward experiments at the Large Hadron Collider during the High Luminosity era. Located along the beam collision axis and shielded from the interaction point by at least 100 m of concrete and rock, the FPF will house experiments that will detect particles outside the acceptance of the existing large LHC experiments and will observe rare and exotic processes in an extremely low-background environment. In this work, we summarize the current status of plans for the FPF, including recent progress in civil engineering in identifying promising sites for the FPF and the experiments currently envisioned to realize the FPF's physics potential. We then review the many Standard Model and new physics topics that will be advanced by the FPF, including searches for long-lived particles, probes of dark matter and dark sectors, high-statistics studies of TeV neutrinos of all three flavors, aspects of perturbative and non-perturbative QCD, and high-energy astroparticle physics.

The strong interactions at low energy scales determine the state of the supranuclear matter in the pulsar-like compact objects. It is proposed that the bulk strong matter could be composed of strangeons, which are quark clusters with a nearly equal number of three light-flavor quarks. In this work, to characterize the strong-repulsive interactions at short distances and the non-relativistic nature of the strangeons, the Lennard-Jones model is used to describe the equation of state (EoS) of strangeon stars (SSs). We investigate the static, the slowly rotating, and the tidally deformed SSs in detail. The corrections resulted from the finite surface densities are considered crucially in the perturbative approaches. We also study the universal relations between the moments of inertia, the tidal deformabilities, and the quadrupole moments. Those results are ready to be used for various purposes in astrophysics, and possible constraints from contemporary observations on the parameter space of the Lennard-Jones model are discussed. Future observations of the pulsars' radio signals, the X-ray emissions from the hot spots on the surface of the stars, and the gravitational waves (GWs) from binary mergers can give tighter constraints or even verify or falsify the existence of SSs.

The AstroCamp is an academic excellence program in the field of astronomy and physics for students in the last 3 years of pre-university education, which often includes a course (or a significant part thereof) on Relativity. After an introduction to the principles, goals and structure of the camp, I describe the approach followed by camp lecturers (myself and others) for teaching Special and General Relativity, and some lessons learned and feedback from the students. I also provide some thoughts on the differences between the physics and mathematics secondary school curricula in Portugal and in other countries, and on how these curricula could be modernized.

We consider the production of dark matter during the process of reheating after inflation. The relic density of dark matter from freeze-in depends on both the energy density and energy distribution of the inflaton scattering or decay products composing the radiation bath. We compare the perturbative and non-perturbative calculations of the energy density in radiation. We also consider the (likely) possibility that the final state scalar products are unstable. Assuming either thermal or non-thermal energy distribution functions, we compare the resulting relic density based on these different approaches. We show that the present-day cold dark matter density can be obtained through freeze-in from preheating for a large range of dark matter masses.

We have recently proposed a simple relativistic theory which reduces to Modified Newtonian Dynamics (MOND) for the weak-field quasistatic situations applied to galaxies, and to cosmological behaviour as in the $\Lambda$CDM model, yielding a realistic cosmology in line with observations. A key requirement of any such model is that Minkowski space is stable against linear perturbations. We expand the theory action to 2nd order in perturbations on a Minkowski background and show that it leads to healthy dispersion relations involving propagating massive modes in the vector and the scalar sector. We use Hamiltonian methods to eliminate constraints present, demonstrate that the massive modes have positive Hamiltonian and show that a non-propagating mode with a linear time dependence may have negative Hamiltonian for wavenumbers $k< k_*$ and positive otherwise. The scale $k_*$ is estimated to be around $\lesssim \mathrm{Mpc}^{-1}$ so that the low momenta instability may only play a role on cosmological scales.

Neutron stars (NSs) in scalar-tensor theories of gravitation with the phenomenon of spontaneous scalarization can develop significant deviations from general relativity. Cases with a massless scalar were studied widely. Here we compare the NS scalarizations in the Damour--Esposito-Far{\`e}se theory, the Mendes-Ortiz theory, and the $\xi$-theory with a massive scalar field. Numerical solutions for slowly rotating NSs are obtained. They are used to construct the X-ray pulse profiles of a pair of extended hot spots on the surface of NSs. We also calculate the tidal deformability for NSs with spontaneous scalarization which is done for the first time with a massive scalar field. We show the universal relation between the moment of inertia and the tidal deformability. The X-ray pulse profiles, the tidal deformability, and the universal relation may help to constrain the massive scalar-tensor theories in X-ray and gravitational-wave observations of NSs, including the Neutron star Interior Composition Explorer (NICER) satellite, Square Kilometre Array (SKA) telescope, and LIGO/Virgo/KAGRA laser interferometers.

The interior of a neutron star is expected to exhibit different states of matter. In particular, complex non-spherical configurations known as `pasta' phases may exist at the highest densities in the inner crust, potentially having an impact on different neutron-star phenomena. We study the properties of the pasta phase and the uncertainties in the pasta observables which are due to our incomplete knowledge of the nuclear energy functional. To this aim, we employed a compressible liquid-drop model approach with surface parameters optimised either on experimental nuclear masses or theoretical calculations. To assess the model uncertainties, we performed a Bayesian analysis by largely varying the model parameters using uniform priors, and generating posterior distributions with filters accounting for both our present low-density nuclear physics knowledge and high-density neutron-star physics constraints. Our results show that the nuclear physics constraints, such as the neutron-matter equation of state at very low density and the experimental mass measurements, are crucial in determining the crustal and pasta observables. Accounting for all constraints, we demonstrate that the presence of pasta phases is robustly predicted in an important fraction of the inner crust. We estimate the relative crustal thickness associated with pasta phases as $R_{\rm pasta}/R_{\rm crust}=0.128\pm 0.047$ and the relative moment of inertia as $I_{\rm pasta}/I_{\rm crust}=0.480\pm 0.137$. Our findings indicate that the surface and curvature parameters are more influential than the bulk parameters for the description of pasta observables. We also show that using a surface tension that is inconsistent with the bulk functional leads to an underestimation of both the average values and the uncertainties in the pasta properties, thus highlighting the importance of a consistent calculation of the nuclear functional.

We discuss the computation of the quantum effective action of strongly interacting field theories using holographic duality, and its use to determine quasi-equilibrium parameters of first order phase transitions relevant for gravitational wave production. A particularly simple holographic model is introduced, containing only the metric and a free massive scalar field. Despite the simplicity, the model contains a rich phase diagram, including first order phase transitions at non-zero temperature, due to various multi-trace deformations. We obtain the leading terms in the effective action from homogeneous black brane solutions in the gravity dual, and linearised perturbations around them. We then employ the effective action to construct bubble and domain wall solutions in the field theory side and study their properties. In particular, we show how the scaling of the effective action with the effective number of degrees of freedom of the quantum field theory determines the corresponding scaling of gravitational wave parameters.