Papers for Thursday, Nov 03 2022

We present the analysis of the atmospheres of 70 gaseous extrasolar planets via transit spectroscopy with Hubble's Wide Field Camera 3 (WFC3). For over half of these, we statistically detect spectral modulation which our retrievals attribute to molecular species. Among these, we use Bayesian Hierarchical Modelling to search for chemical trends with bulk parameters. We use the extracted water abundance to infer the atmospheric metallicity and compare it to the planet's mass. We also run chemical equilibrium retrievals, fitting for the atmospheric metallicity directly. However, although previous studies have found evidence of a mass-metallicity trend, we find no such relation within our data. For the hotter planets within our sample, we find evidence for thermal dissociation of dihydrogen and water via the H$^-$ opacity. We suggest that the general lack of trends seen across this population study could be due to i) the insufficient spectral coverage offered by HST WFC3 G141, ii) the lack of a simple trend across the whole population, iii) the essentially random nature of the target selection for this study or iv) a combination of all the above. We set out how we can learn from this vast dataset going forward in an attempt to ensure comparative planetology can be undertaken in the future with facilities such as JWST, Twinkle and Ariel. We conclude that a wider simultaneous spectral coverage is required as well as a more structured approach to target selection.

Hubble Space Telescope (HST) photometry is providing an extensive analysis of globular clusters (GCs). In particular, the pseudo two-colour diagram dubbed 'chromosome map (ChM)' allowed to detect and characterize their multiple populations with unprecedented detail. The main limitation of these studies is the small field of view of HST, which makes it challenging to investigate some important aspects of the multiple populations, such as their spatial distributions and the internal kinematics in the outermost cluster regions. To overcome this limitation, we analyse state-of-art wide-field photometry of 43 GCs obtained from ground-based facilities. We derived high-resolution reddening maps and corrected the photometry for differential reddening when needed. We use photometry in the U, B, and I bands to introduce the $\Delta c_{\rm U,B,I}$ vs. $\Delta_{\rm B,I}$ ChM of red-giant branch (RGB) and asymptotic-giant branch (AGB) stars. We demonstrate that this ChM, which is built with wide-band ground-based photometry, is an efficient tool to identify first- and second-generation stars (1G and 2G) over a wide field of view. To illustrate its potential, we derive the radial distribution of multiple populations in NGC 288 and infer their chemical composition. We present the ChMs of RGB stars in 29 GCs and detect a significant degree of variety. The fraction of 1G and 2G stars, the number of subpopulations, and the extension of the ChMs significantly change from one cluster to another. Moreover, the metal-poor and metal-rich stars of Type II GCs define distinct sequences in the ChM. We confirm the presence of extended 1G sequences.

We present Hubble Space Telescope imaging of 14 gas-rich, low surface brightness and ultra-diffuse galaxies (UDGs) in the field at distances of 25-36 Mpc. An inspection of point-like sources brighter than the turnover magnitude of the globular cluster luminosity function and within twice the half-light radii of each galaxy reveals that, unlike those in denser environments, gas-rich, field UDGs host very few old globular clusters (GCs). Most of the targets (nine) have zero candidate GCs, with the remainder having one or two candidates each. These findings are broadly consistent with expectations for normal dwarf galaxies of similar stellar mass. This rules out gas-rich, field UDGs as potential progenitors of the GC-rich UDGs that are typically found in galaxy clusters. However, some in galaxy groups may be directly accreted from the field. In line with other recent results, this strongly suggests that there must be at least two distinct formation pathways for UDGs, and that this sub-population is simply an extreme low surface brightness extension of the underlying dwarf galaxy population. The root cause of their diffuse stellar distributions remains unclear, but the formation mechanism appears to only impact the distribution of stars (and potentially dark matter), without strongly impacting the distribution of neutral gas, the overall stellar mass, or the number of GCs.

The interaction between Active Galactic Nuclei (AGN) and their host galaxies is scarcely resolved. Narrow-line Seyfert 1 (NLS1) galaxies are believed to represent AGN at early stages of their evolution and allow to observe AGN feeding and feedback processes at high accretion rates. We apply a spectroastrometric analysis to VLT MUSE NFM-AO observations of Mrk 1044, a nearby super-Eddington accreting NLS1. This allows us to map two ionised gas outflows traced by [O$\,$III] which have velocities of $-560\pm20\,{\rm km\:s}^{-1}$ and $-144 \pm 5 \,{\rm km\:s}^{-1}$. Both outflows are spatially unresolved and located close to the galaxy nucleus ($<1\,{\rm pc}$). They have gas densities higher than $10^5\,{\rm cm}^{-3}$, which implies that the BPT diagnostic cannot be used to constrain the underlying ionisation mechanism. We explore whether an expanding shell model can describe the velocity structure of Mrk 1044's unresolved multi-phase outflow. A kinematic analysis suggests that significant turbulence may be present in the ISM around the nucleus, which may lead to a condensation rain, potentially explaining the efficient feeding of Mrk 1044's AGN. We identify an additional ionised gas outflowing component that is spatially resolved. It has a velocity of $-211 \pm 22 \,{\rm km\:s}^{-1}$ and a projected size of $4.6 \pm 0.6 \,{\rm pc}$. Within the innermost 0.5" (160$\,{\rm pc}$) around the nucleus we detect modest star formation hidden by the beam-smeared emission from the outflow, which suggests that Mrk 1044's AGN phase set on recently. We estimate that the multi-phase outflow has been launched $<10^4 \,{\rm yrs}$ ago. It carries enough mass and energy to impact the host galaxy star formation on different spatial scales, highlighting the complexity of the AGN feeding and feedback cycle in its early stages.

A large fraction of gamma-ray bursts (GRBs) shows a plateau phase during the X-ray afterglow emission, whose physical origin is still debated. In this work we define a sample of 30 GRBs with simultaneous X-ray and optical data during and after the plateau phase. Through a time-resolved spectral analysis of the X-ray plateaus, we test the consistency of the unabsorbed optical fluxes with those obtained via X-ray-to-optical spectral extrapolation by assuming a synchrotron spectrum. Combining X-ray with optical data, we find that 63% (19/30) GRBs are compatible with a single synchrotron spectrum thus suggesting that both the optical and X-ray radiations are produced from a single emitting region. For these GRBs we derive the temporal evolution of the break frequency and we compare it with the expectations predicted by several models. For 11/30 GRBs the optical emission is above the predicted range of values extrapolated from the X-rays in at least one temporal bin of the light curve. These GRBs may not be explained with a single zone emission, indicating the necessity of invoking two cooperating processes in order to explain the broad band spectral behaviour during X-ray plateaus. We discuss our findings in the framework of different scenarios invoked to explain the plateau feature, including the energy injection from a spinning-down magnetar and the high latitude emission from a structured jet.

We report the detection of a 219 $^{+10}_{-28}$ pc-sized dark matter core structure in the center of Milky Way at one sigma level by using the micro-lensing event rate sky map data from the Optical Gravitational Lensing Experiment (OGLE) survey. For the first time, we apply the spacial information of the micro-lensing sky map and model it with the detailed Milky Way structure, the Mini Dark Matter Structure (MDMS) fraction ($f_{\rm MDMS}=\Omega_{\rm MDMS}/\Omega_{\rm DM}$) and the core size.We find that this sky map can constrain both $f_{\rm MDMS}$ and the core size simultaneously without strong degeneracy.This discovery provides not only guidance for dark matter particle models, such as self-interacting dark matter (SIDM), but also the baryonic physics of Milky Way.

HD 191939 (TOI-1339) is a nearby (d=54pc), bright (V=9mag), and inactive Sun-like star (G9 V) known to host a multi-planet transiting system. Ground-based spectroscopic observations confirmed the planetary nature of the three transiting sub-Neptunes (HD 191939 b, c, and d) originally detected by TESS and were used to measure the masses for planets b and c with 3$\sigma$ precision. These previous observations also reported the discovery of an additional Saturn-mass planet (HD 191939 e) and evidence for a further, very long-period companion (HD 191939 f). Here, we report the discovery of a new non-transiting planet in the system and a refined mass determination of HD 191939 d. The new planet, HD 191939 g, has a minimum mass of 13.5$\pm$2.0 M$_\oplus$ and a period of about 280 d. This period places the planet within the conservative habitable zone of the host star, and near a 1:3 resonance with HD 191939 e. The compilation of 362 radial velocity measurements with a baseline of 677 days from four different high-resolution spectrographs also allowed us to refine the properties of the previously known planets, including a 4.6$\sigma$ mass determination for planet d, for which only a 2$\sigma$ upper limit had been set until now. We confirm the previously suspected low density of HD 191939 d, which makes it an attractive target for attempting atmospheric characterisation. Overall, the planetary system consists of three sub-Neptunes interior to a Saturn-mass and a Uranus-mass planet plus a high-mass long-period companion. This particular configuration has no counterpart in the literature and makes HD 191939 an exceptional multi-planet transiting system with an unusual planet demographic worthy of future observation.

Temperature measurements of galaxy clusters are used to determine their masses, which in turn are used to determine cosmological parameters. However, systematic differences between the temperatures measured by different telescopes imply a significant source of systematic uncertainty on such mass estimates. We perform the first systematic comparison between cluster temperatures measured with Chandra and NuSTAR. This provides a useful contribution to the effort of cross-calibrating cluster temperatures due to the harder response of NuSTAR compared with most other observatories. We measure average temperatures for 8 clusters observed with NuSTAR and Chandra. We fit the NuSTAR spectra in a hard (3-10 keV) energy band, and the Chandra spectra in both the hard and a broad (0.6-9 keV) band. We fit a power-law cross-calibration model to the resulting temperatures. At a Chandra temperature of 10 keV, the average NuSTAR temperature was $(10.5 \pm 3.7)\%$ and $(15.7 \pm 4.6)\%$ lower than Chandra for the broad and hard band fits respectively. We explored the impact of systematics from background modelling and multiphase temperature structure of the clusters, and found that these did not affect our results. Our sample are primarily merging clusters with complex thermal structures so are not ideal calibration targets. However, given the harder response of NuSTAR it would be expected to measure a higher average temperature than Chandra for a non-isothermal cluster, so we interpret our measurement as a lower limit on the difference in temperatures between NuSTAR and Chandra.

We describe methods to measure simultaneously the orientation angle $\psi$ and pattern speed $\Omega=\mathrm{d}\psi/\mathrm{d}t$ from single snapshots of simulated barred galaxies. Unlike previous attempts, which suffer from systematic errors of 5-25%, our approach is unbiased. It can be extended to obtain the rate and axis of rotation, i.e. the vector $\boldsymbol{\Omega}$. We provide computer code implementing our method.

In the era of big astronomical surveys, our ability to leverage artificial intelligence algorithms simultaneously for multiple datasets will open new avenues for scientific discovery. Unfortunately, simply training a deep neural network on images from one data domain often leads to very poor performance on any other dataset. Here we develop a Universal Domain Adaptation method DeepAstroUDA, capable of performing semi-supervised domain alignment that can be applied to datasets with different types of class overlap. Extra classes can be present in any of the two datasets, and the method can even be used in the presence of unknown classes. For the first time, we demonstrate the successful use of domain adaptation on two very different observational datasets (from SDSS and DECaLS). We show that our method is capable of bridging the gap between two astronomical surveys, and also performs well for anomaly detection and clustering of unknown data in the unlabeled dataset. We apply our model to two examples of galaxy morphology classification tasks with anomaly detection: 1) classifying spiral and elliptical galaxies with detection of merging galaxies (three classes including one unknown anomaly class); 2) a more granular problem where the classes describe more detailed morphological properties of galaxies, with the detection of gravitational lenses (ten classes including one unknown anomaly class).

Studies of gamma-ray bursts (GRBs) and their multi-wavelength afterglows have led to insights in electron acceleration and emission properties from relativistic, high-energy astrophysical sources. Broadband modeling across the electromagnetic spectrum has been the primary means of investigating the physics behind these sources, although independent diagnostic tools have been developed to inform and corroborate assumptions made in particle acceleration simulations and broadband studies. We present a methodology to constrain three physical parameters related to electron acceleration in GRB blast waves: the fraction of shock energy in electrons, $\epsilon_e$; the fraction of electrons that gets accelerated into a power-law distribution of energies, $\xi_e$; and the minimum Lorentz factor of the accelerated electrons, $\gamma_m$. These parameters are constrained by observations of the peaks in radio afterglow light curves and spectral energy distributions. From a sample of 49 radio afterglows, we are able to find narrow distributions for these parameters, hinting at possible universality of the blast wave microphysics, although observational bias could play a role in this. Using radio peaks and considerations related to the prompt gamma-ray emission efficiency, we constrain the allowed parameter ranges for both $\epsilon_e$ and $\xi_e$ to within about one order of magnitude, $0.01\lesssim\epsilon_e\lesssim0.2$ and $0.1\lesssim\xi_e\lesssim1$. Such stringent constraints are inaccessible for $\xi_e$ from broadband studies due to model degeneracies.

How ions are energized and heated is a fundamental problem in the study of energy dissipation in magnetized plasmas. In particular, the heating of heavy ions (including ${}^{4}\mathrm{He}^{2+}$, ${}^{3}\mathrm{He}^{2+}$ and others) has been a constant concern for understanding the microphysics of impulsive solar flares. In this article, via two-dimensional hybrid-kinetic Particle-in-Cell simulations, we study the heating of Helium ions (${}^{4}\mathrm{He}^{2+}$) by turbulence driven by cascading waves launched at large scales from the left-handed polarized Helium ion cyclotron wave branch of a multi-ion plasma composed of electrons, protons, and Helium ions. We find significant parallel (to the background magnetic field) heating for both Helium ions and protons due to the formation of beams and plateaus in their velocity distribution functions along the background magnetic field. The heating of Helium ions in the direction perpendicular to the magnetic field starts with a lower rate than that in the parallel direction, but overtakes the parallel heating after a few hundreds of the proton gyro-periods due to cyclotron resonances with mainly obliquely propagating waves induced by the cascade of injected Helium ion cyclotron waves at large scales. There is however little evidence for proton heating in the perpendicular direction due to the absence of left-handed polarized cyclotron waves near the proton cyclotron frequency. Our results are useful for understanding the preferential heating of ${}^{3}\mathrm{He}$ and other heavy ions in the ${}^{3}\mathrm{He}$-rich solar energetic particle events, in which Helium ions play a crucial role as a species of background ions regulating the kinetic plasma behavior.

Young star clusters enable us to study the effects of stellar rotation on an ensemble of stars of the same age and across a wide range in stellar mass and are therefore ideal targets for understanding the consequences of rotation on stellar evolution. We combine MUSE spectroscopy with HST photometry to measure the projected rotational velocities (Vsini) of 2,184 stars along the split main sequence and on the main sequence turn-off (MSTO) of the 100 Myr-old massive (10^5 M_sun) star cluster NGC 1850 in the Large Magellanic Cloud. At fixed magnitude, we observe a clear correlation between Vsini and colour, in the sense that fast rotators appear redder. The average Vsini values for stars on the blue and red branches of the split main sequence are ~100 km/s and ~200 km/s, respectively. The values correspond to about 25-30% and 50-60% of the critical rotation velocity and imply that rotation rates comparable to those observed in field stars of similar masses can explain the split main sequence. Our spectroscopic sample contains a rich population of ~200 fast rotating Be stars. The presence of shell features suggests that 23% of them are observed through their decretion disks, corresponding to a disk opening angle of 15 degrees. These shell stars can significantly alter the shape of the MSTO, hence care should be taken when interpreting this photometric feature. Overall, our findings impact our understanding of the evolution of young massive clusters and provide new observational constraints for testing stellar evolutionary models.

Understanding planet formation requires robust population studies, which are designed to reveal trends in planet properties. In this work, we aim to determine if different methods for selecting populations of exoplanets for atmospheric characterization with JWST could influence population-level inferences. We generate three hypothetical surveys of super-Earths/sub-Neptunes, each spanning a similar radius-insolation flux space. The survey samples are constructed based on three different selection criteria (evenly-spaced-by-eye, binned, and a quantitative selection function). Using an injection-recovery technique, we test how robustly individual-planet atmospheric parameters and population-level parameters can be retrieved. We find that all three survey designs result in equally suitable targets for individual atmospheric characterization, but not equally suitable targets for constraining population parameters. Only samples constructed with a quantitative method or that are sufficiently evenly-spaced-by-eye result in robust population parameter constraints. Furthermore, we find that the sample with the best targets for individual atmospheric study does not necessarily result in the best constrained population parameters. The method of sample selection must be considered. We also find that there may be large variability in population-level results with a sample that is small enough to fit in a single JWST cycle ($\sim$12 planets), suggesting that the most successful population-level analyses will be multi-cycle. Lastly, we infer that our exploration of sample selection is limited by the small number of transiting planets with measured masses around bright stars. Our results can guide future development of programs that aim to determine underlying trends in exoplanet atmospheric properties and, by extension, formation and evolution processes.

Kepler-33 hosts five validated transiting planets ranging in period from 5 to 41 days. The planets are in nearly co-planar orbits and exhibit remarkably similar (appropriately scaled) transit durations indicative of similar impact parameters. The outer three planets have radii of $3.5\lesssim R_{\rm p}/R_\oplus\lesssim4.7$ and are closely-packed dynamically, and thus transit timing variations can be observed. Photodynamical analysis of transit timing variations provide $2\sigma$ upper bounds on the eccentricity of the orbiting planets (ranging from $<0.02$ to $<0.2$) and the mean density of the host-star ($0.39_{-0.02}^{+0.01}\,{\rm g/cm^3}$). We combine \emph{Gaia} Early Data Release 3 parallax observations, the previously reported host-star effective temperature and metallicity, and our photodynamical model to refine properties of the host-star and the transiting planets. Our analysis yields well-constrained masses for Kepler-33~e ($6.6_{-1.0}^{+1.1}\,M_\oplus$) and f ($8.2_{-1.2}^{+1.6}\,M_\oplus$) along with $2\sigma$ upper limits for planets c ($<19\,M_\oplus$) and d ($<8.2\,M_\oplus$). We confirm the reported low bulk densities of planet d ($<0.4\,{\rm g/cm^3}$), e ($0.8\pm0.1\,{\rm g/cm^3}$), and f ($0.7\pm0.1\,{\rm g/cm^3}$). Based on comparisons with planetary evolution models, we find that Kepler-33~e and f exhibit relatively high envelope mass fractions of $f_{\rm env}=7.0_{-0.5}^{+0.6}\%$ and $f_{\rm env}=10.3\pm0.6\%$, respectively. Assuming a mass for planet d $\sim4\,M_\oplus$ suggests that it has $f_{\rm env}\gtrsim12\%$.

The classical radiation pressure instability has been a persistent theoretical feature of thin, radiatively efficient accretion disks with accretion rates 1 to 100 per cent of the Eddington rate. But there is only limited evidence of its occurrence in nature: rapid heartbeat oscillations of a few X-ray binaries and now, perhaps, the new class of hourly X-ray transients called quasi-periodic eruptions (QPEs). The accretion disks formed in tidal disruption events (TDEs) have been observed to peacefully trespass through the range of unstable accretion rates without exhibiting any clear sign of the instability. We try to explain the occurrence or otherwise of this instability in these systems, by constructing steady state 1D models of thin magnetic accretion disks. The local magnetic pressure in the disk is assumed to be dominated by toroidal fields arising from a dynamo sourced by magneto-rotational instability (MRI). We choose a physically motivated criterion of MRI saturation, validated by recent magnetohydrodynamic simulations, to determine the strength of magnetic pressure in the disk. The resulting magnetic pressure support efficiently shrinks: (1) the parameter space of unstable mass accretion rates, explaining the absence of instability in systems such as TDEs and (2) the range of unstable radii in the inner accretion disk, which can shorten the quasi-periods of instability limit-cycles by more than three orders of magnitude, explaining the observed periods ( a few hrs) of QPEs. In addition to examining stability properties of strongly magnetized disks, we predict other observational signatures such as spectral hardening factors and jet luminosities to test the compatibility of our disk models with observations of apparently stable TDE disks.

Enceladus, an icy moon of Saturn, possesses an internal water ocean and jets expelling ocean material into space. Cassini investigations indicated that the subsurface ocean could be a habitable environment having a complex interaction with the rocky core. Further investigation of the composition of the plume formed by the jets is necessary to fully understand the ocean, its potential habitability, and what it tells us about Enceladus' origin. Moonraker has been proposed as an ESA M-class mission designed to orbit Saturn and perform multiple flybys of Enceladus, focusing on traversals of the plume. The proposed Moonraker mission consists of an ESA-provided platform, with strong heritage from JUICE and Mars Sample Return, and carrying a suite of instruments dedicated to plume and surface analysis. The nominal Moonraker mission has a duration of 13.5 years. It includes a 23-flyby segment with 189 days allocated for the science phase, and can be expanded with additional segments if resources allow. The mission concept consists in investigating: i) the habitability conditions of present-day Enceladus and its internal ocean, ii) the mechanisms at play for the communication between the internal ocean and the surface of the South Polar Terrain, and iii) the formation conditions of the moon. Moonraker, thanks to state-of-the-art instruments representing a significant improvement over Cassini's payload, would quantify the abundance of key species in the plume, isotopic ratios, and physical parameters of the plume and the surface. Such a mission would pave the way for a possible future landed mission.

At a global scale, the universe is generally fitted by an idealized manifold described by the FLRW metric. This is in particular the case when probing the universe to determine its dynamics. The process that fits the idealized manifold to the real universe is however not uniquely defined. This process may depend on the cosmic probe that has been used for the measurements, and could hence lead to different observed temporal evolutions of the scale factor. A correct interpretation of the observed accelerated expansion of the universe requires therefore first a thorough understanding of the fitting process that has been implicitly applied. In this article we establish the fitting processes for the SNIa, BAO and CMB cosmic probes, and deduce the related averaged Einstein equations. We demonstrate that the way these fittings have been applied in practice lead to an apparent dark energy effect. We also highlight the conceptual differences in the fitting processes between the SNIa and BAO probes on the one hand, and the CMB probe on the other hand. Considering those differences, we then show how the so-called Hubble tension can be explained.

Active galactic nuclei (AGN) are accreting supermassive black holes co-evolving with their host galaxies through a complex interplay of feeding and feedback. In this chapter, we first discuss AGN fuelling in galaxies, both in interacting and isolated systems, focusing on the role that instabilities have on the angular momentum budget of the gas. We then review observations and models of feedback through AGN-driven winds from nuclear, sub-pc scales out to galactic and circumgalactic medium scales. We continue with an overview of surveys and statistical properties of the AGN population, before concluding with a discussion on the prospects for the future facilities, focusing in particular on Athena.

The scarce optical variability studies in spectrally classified Type 2 active galactic nuclei (AGNs) have led to the discovery of anomalous objects that are incompatible with the simplest unified models (UM). This paper focuses on the exploration of different variability features that allows to separate between obscured, Type 2 AGNs, and the variable, unobscured Type 1s. We analyse systematically the Zwicky Transient Facility, 2.5 years long light curves of ~ 15000 AGNs from the Sloan Digital Sky Survey Data Release 16, which are generally considered Type 2s due to the absence of strong broad emission lines (BELs). Consistently with the expectations from the UM, the variability features are distributed differently for distinct populations, with spectrally classified weak Type 1s showing 1 order of magnitude larger variances than the Type 2s. We find that the parameters given by the damped random walk model leads to broader H{\alpha} equivalent width for objects with {\tau}_g > 16 d and long term structure function SF{\infty},g> 0.07 mag. By limiting the variability features, we find that ~ 11 per cent of Type 2 sources show evidence for optical variations. A detailed spectral analysis of the most variable sources (~1 per cent of the Type 2 sample) leads to the discovery of misclassified Type 1s with weak BELs and changing-state candidates. This work presents one of the largest systematic investigations of Type 2 AGN optical variability to date, in preparation for future large photometric surveys.

The primary goal of this thesis was the correction of Non-Common-Path-Aberrations in the SPHERE instrument for helping it meeting its contrast requirements. SPHERE's purpose is the search and characterization of giant exo-planets around nearby stars. The author implemented a method called Electric Field Conjugation that he tested in simulation as well as on the real system. A full week was booked in SPHERE schedule a few days before the second commissioning in June 2014. It gave the opportunity to the author to travel to the VLT in Chile and experiment directly on the system. The contrast gain objective of another order of magnitude in a medium-sized area has successfully been reached bringing SPHERE raw speckle contrast from about $10^{-6}$ to $10^{-7}$. The algorithm has therefore proven its value and will be further investigated and hopefully automated by the SPHERE team based on the codes developed by the author. However it is important to keep in mind that Electric Field Conjugation is more effective for follow-up studies in order to improve the quality of the observations. Indeed the area for a good correction is very limited. It can't be used for exo-planets discovery unless the corrected area is made big enough but the performance will be less.

We present our new Atacama Large Millimeter/submillimeter Array (ALMA) observations targeting CO(6-5) emission from three luminous Lyman break galaxies (LBGs) at $z_{\rm spec} = 6.0293$-$6.2037$ found in the Subaru/Hyper Suprime-Cam survey, whose [OIII]$88\mu$m and [CII]$158\mu$m emission have been detected with ALMA. We find a marginal detection of the CO(6-5) line from one of our LBGs, J0235-0532, at the $\simeq 4 \sigma$ significance level and obtain upper limits for the other two LBGs, J1211-0118 and J0217-0208. Our $z=6$ luminous LBGs are consistent with the previously found correlation between the CO luminosity and the infrared luminosity. The unique ensemble of the multiple far-infrared emission lines and underlying continuum fed to a photodissociation region model reveal that J0235-0532 has a relatively high hydrogen nucleus density that is comparable to those of low-$z$ (U)LIRGs, quasars, and Galactic star-forming regions with high $n_{\rm H}$ values, while the other two LBGs have lower $n_{\rm H}$ consistent with local star-forming galaxies. By carefully taking account of various uncertainties, we obtain total gas mass and gas surface density constraints from their CO luminosity measurements. We find that J0235-0532 locates below the Kennicutt-Schmidt (KS) relation, comparable to the previously CO(2-1) detected $z=5.7$ LBG, HZ10. Combined with previous results for dusty starbursts at similar redshifts, the KS relation at $z=5$-$6$ is on average consistent with the local one.

We report $I$-band photometric observations of the radio-detected M9.5 dwarf BRI 0021-0214, obtained with the Galway Ultra Fast Imager (GUFI) on the 1.8m Vatican Advanced Technology Telescope VATT at Mt. Graham International Observatory, Arizona. In total, 19 hours of observations over a 73 day baseline were obtained. BRI 0021-0214 was shown to exhibit modulated emission with a period of $ 3.052 \pm 0.004$ hours with a mean amplitude variability of 0.0044 mag. When combined with rotational velocity data obtained from previous work, our newly discovered rotation period gives an inclination angle of 51.7$^{+5.0}_{-4.5}$ degrees for the rotation axis of BRI 0021-0214 relative to our line of sight. Previous studies have reported that the most plausible cause for optical variability from this dwarf is a consequence of suspended co-rotating dust clouds in its atmosphere. However reports of enhanced H$_{\alpha}$ and intermittent coherent radio emission suggest the possibility of auroral activity in its magnetosphere. Further, more coordinated multiwavlength observations of this dwarf could fully resolve the nature of this elusive rapid-rotator object's observational properties.

The Atacama Large Millimetre/Sub-millimetre Array (ALMA) is obtaining the deepest observations of early galaxies ever achieved at (sub-)millimetre wavelengths, and detecting the dust emission of young galaxies in the first billion years of cosmic history, well in the epoch of reionization. Here I review some of the latest results from these observations, with special focus on the REBELS large programme, which targets a sample of 40 star-forming galaxies at z~7. ALMA detects significant amounts of dust in very young galaxies, and this dust might have different properties to dust in lower-redshift galaxies. I describe the evidence for this, and discuss theoretical/modelling efforts to explain the dust properties of these young galaxies. Finally, I describe two additional surprising results to come out of the REBELS survey: (i) a new population of completely dust-obscured galaxies at z~7, and (ii) the prevalence of spatial offsets between the ultraviolet and infrared emission of UV-bright, high-redshift star-forming galaxies.

We explore the effects of electromagnetic (EM) fluctuations in plasmas on solar neutrino fluxes exploiting the fluctuation-dissipation theorem. We find that the EM spectrum in the solar core is enhanced by the EM fluctuations due to the high density of the Sun, which increases the radiation energy density and pressure. By the EM fluctuations involving the modified radiation formula, the central temperature decreases when the central pressure of the Sun is fixed. With a help of the empirical relation between central temperature and neutrino fluxes deduced from the numerical solar models, we present the change in each of the solar neutrino fluxes by the EM fluctuations. We also discuss the enhanced radiation pressure and energy density by the EM fluctuations for other astronomical objects.

We present the first Gaia catalogue of eclipsing binary candidates released in Gaia DR3, describe its content, provide tips for its usage, estimate its quality, and show illustrative samples. The catalogue contains 2,184,477 sources with G magnitudes up to 20 mag. Candidate selection is based on the results of variable object classification performed within the Gaia Data Processing and Analysis Consortium, further filtered using eclipsing binary-tailored criteria based on the G light curves. To find the orbital period, a large ensemble of trial periods is first acquired using three distinct period search methods applied to the cleaned G light curve. The G light curve is then modelled with up-to two Gaussians and a cosine for each trial period. The best combination of orbital period and geometric model is finally selected using Bayesian model comparison based on the BIC. A global ranking metric is provided to rank the quality of the chosen model between sources. The catalogue is restricted to orbital periods larger than 0.2 days. About 530,000 of the candidates are classified as eclipsing binaries in the literature as well, out of ~600,000 available crossmatches, and 93% of them have published periods compatible with the Gaia periods. Catalogue completeness is estimated to be between 25% and 50%, depending on the sky region, relative to the OGLE4 catalogues of eclipsing binaries towards the Galactic Bulge and the Magellanic Clouds. The analysis of an illustrative sample of ~400,000 candidates with significant parallaxes shows properties in the observational HR diagram as expected for eclipsing binaries. The subsequent analysis of a sub-sample of detached bright candidates provides further hints for the exploitation of the catalogue. The orbital periods, light curve model parameters, and global rankings are all published in the catalogue with their related uncertainties where applicable.

Based on the recent advancements in the numerical simulations of galaxy formation, we anticipate the achievement of realistic models of galaxies in the near future. Morphology is the most basic and fundamental property of galaxies, yet observations and simulations still use different methods to determine galaxy morphology, making it difficult to compare them. We hereby perform a test on the recent NewHorizon simulation which has spatial and mass resolutions that are remarkably high for a large-volume simulation, to resolve the situation. We generate mock images for the simulated galaxies using SKIRT that calculates complex radiative transfer processes in each galaxy. We measure morphological indicators using photometric and spectroscopic methods following observer's techniques. We also measure the kinematic disk-to-total ratios using the Gaussian mixture model and assume that they represent the true structural composition of galaxies. We found that spectroscopic indicators such as $V/{\sigma}$ and ${\lambda}_{R}$ closely trace the kinematic disk-to-total ratios. In contrast, photometric disk-to-total ratios based on the radial profile fitting method often fail to recover the true kinematic structure of galaxies, especially for small galaxies. We provide translating equations between various morphological indicators.

Since the unexpected discovery of fast radio bursts (FRBs), researchers have proposed varied theories and models to explain their phenomena. One such model that has recently been developed incorporates the so-called Gertsenshtein-Zel'dovich (GZ) effect, which states that when gravitational waves traverse the pulsar magnetosphere, a portion of their gravitational radiation is transformed into electromagnetic (EM) radiation. The observed properties of FRBs are consistent with the properties of this EM radiation, implying, remarkably, that the GZ effect can account for both repeating and non-repeating FRBs. If this model is correct, the pulsar's properties should not change over time, and it would continue to emit both EM dipole and gravitational quadrupole radiation for a long period of time. This article targets the gravitational radiation produced by the pulsar mechanism and shows that several proposed gravitational wave detectors can detect these gravitational waves. If such detections are performed in the future from the positions of FRBs, it will validate the GZ process for FRB production and can rule out several other FRB theories.

While the temperature of the X-ray corona ($\rm{kT_e}$) in active galactic nuclei (AGN) are known for many sources, its variation, if any, is limited to a handful of objects. This is in part due to the requirement of good signal-to-noise X-ray spectra covering a wide range of energies. We present here results on the X-ray spectral analysis of 18 Seyferts, having more than one epoch of observations to look for variation in $\rm{kT_e}$. The data for a total of 52 epochs on these 18 AGN were taken from observations carried out by NuSTAR in the 3$-$79 keV energy band. From phenomenological and physical model fits to the multi-epoch data on these 18 sources from {\it NuSTAR}, we could constrain the cut-off energy ($E_{cut}$) in a large fraction of the sources. Also, from Comptonized model fits, we could obtain $\rm{kT_e}$ for our sample. Of the 18 sources, at the 90 per cent confidence level, evidence for variation in $\rm{kT_e}$ was found for only one source, namely MCG+08-11-011. For this source, between two epochs, separated by about five years, we found $\rm{kT_e}$ to decrease from 57$^{+29}_{-16}$ keV to 30$^{+11}_{-7}$ keV. During the same period, the flux decreased from (12.60 to 14.02) $\times$ 10$^{-11}$ erg cm$^{-2}$ s$^{-1}$ and the optical depth increased from 1.68 to 2.73. We thus found a positive correlation between flux and coronal temperature with a reduction of about 40 per cent in optical depth. Our observations tend to favour the vertically outflowing corona scenario for the observed variation in $\rm{kT_e}$ in MCG+08-11-011.

Planet(s) in binaries are unique architectures for testing predictions of planetary formation and evolution theories in very hostile environments. We used the IRDIS dual-band imager of SPHERE at VLT, and the speckle interferometric camera HRCAM of SOAR, to acquire high-angular resolution images of HD 196885 AB between 2015 and 2020. Radial velocity observations have been extended over almost 40 yr extending the radial velocity measurements HD 196885 A and resolving both the binary companion and the inner giant planet HD 196885 Ab. Finally, we took advantage of the exquisite astrometric precision of the dual-field mode of VLTI/GRAVITY (down to 30 {\mu}as) to monitor the relative position of HD 196885 A and B to search for the 3.6 yr astrometric wobble of the circumprimary planet Ab imprinted on the binary separation. Our observations enable to accurately constrain the orbital properties of the binary HD 196885 AB, seen on an inclined and retrograde orbit (iAB = 120.43 deg) with a semi-major axis of 19.78 au, and an eccentricity of 0.417. The GRAVITY measurements confirm for the first time the nature of the inner planet HD 196885 Ab by rejecting all families of pole-on solutions in the stellar or brown dwarf masses. The most favored island of solutions is associated with a Jupiter-like planet (MAb = 3.39 MJup), with moderate eccentricity (eAaAb = 0.44), and inclination close to 143.04 deg. This results points toward a significant mutual inclination (Phi = 24.36 deg) between the orbital planes (relative to the star) of the binary companion B and the planet Ab. Our dynamical simulations indicate that the system is dynamically stable over time. Eccentricity and mutual inclination variations could be expected for moderate von Zipele Kozai Lidov cycles that may affect the inner planet.

We employ our Bayesian Machine Learning framework BINGO (Bayesian INference for Galactic archaeOlogy) to obtain high-quality stellar age estimates for 68,360 red giant and red clump stars present in the 17th data release of the Sloan Digital Sky Survey, the APOGEE-2 high-resolution spectroscopic survey. By examining the denoised age-metallicity relationship of the Galactic disc stars, we identify a drop in metallicity with an increase in [Mg/Fe] at an early epoch, followed by a chemical enrichment episode with increasing [Fe/H] and decreasing [Mg/Fe]. This result is congruent with the chemical evolution induced by an early-epoch gas-rich merger identified in the Milky Way-like zoom-in cosmological simulation Auriga. In the initial phase of the merger of Auriga 18 there is a drop in metallicity due to the merger diluting the metal content and an increase in the [Mg/Fe] of the primary galaxy. Our findings suggest that the last massive merger of our Galaxy, the Gaia-Sausage-Enceladus, was likely a significant gas-rich merger and induced a starburst, contributing to the chemical enrichment and building of the metal-rich part of the thick disc at an early epoch.

We carry out detailed spectral and timing analyses of the $Chandra$ X-ray data of HD 179949, a prototypical example of a star with a close-in giant planet with possible star-planet interaction (SPI) effects. We find a low coronal abundance Fe/H$\approx$0.2 relative to the solar photosphere, as well as lower abundances of high FIP elements O/Fe $\lesssim$1, Ne/Fe $\lesssim$ 0.1, but with indications of higher abundances of N and Al. This star also has an anomalous FIP bias of $\approx 0.03 \pm 0.03$, larger than expected for stars of this type. We detect significant intensity variability over time scales ranging from 100~s - 10~ks, and also evidence for spectral variability over time scales of 1-10~ks. We combine the $Chandra$ flux measurements with $Swift$ and $XMM-Newton$ measurements to detect periodicities, and determine that the dominant signal is tied to the stellar polar rotational period, consistent with expectations that the corona is rotational-pole dominated. We also find evidence for periodicity at both the planetary orbital frequency and at its beat frequency with the stellar polar rotational period, suggesting the presence of a magnetic connection between the planet and the stellar pole. If these periodicities represent an SPI signal, the lack of phase dependence in coronal temperature or flaring suggests that the SPI in this system is driven by a quasi-continuous form of heating (e.g., magnetic field stretching) rather than a highly sporadic, hot, impulsive form (e.g., flare-like reconnection).

MAGIC is a system of two imaging atmospheric Cherenkov telescopes (IACTs) located on the Canary island of La Palma. Each telescope's imaging camera consists of 1039 photomultiplier tubes (PMTs). We developed three detector modules based on silicon photomultipliers (SiPMs) of seven pixels each that are mechanically and electronically compatible with those used in the MAGIC camera. These prototype modules are installed next to the PMTs in the imaging camera and are operated in parallel. To achieve a similar active area per pixel we used seven to nine SiPMs for producing a composite pixel. The SiPM signals within one such pixel are actively summed up for retaining the fast signal pulse shapes. Two different PCB designs are tested for thermal performance. We present our simulations of Cherenkov and light of the night sky (LoNS) responses. Based on those we calculate the signal-to-noise ratio (SNR) for this imaging application. We compare our expectations with the measurements of one of the SiPM-based detector modules.

According to dynamo theory, stars with convective envelopes efficiently generate surface magnetic fields, which manifest as magnetic activity in the form of starspots, faculae, flares, when their rotation period is shorter than their convective turnover time. Most red giants, having undergone significant spin down while expanding, have slow rotation and no spots. However, based on a sample of about 4500 red giants observed by the NASA Kepler mission, a previous study showed that about 8 % display spots, including about 15 % that belong to close binary systems. Here, we shed light on a puzzling fact: for rotation periods less than 80 days, a red giant that belongs to a close binary system displays a photometric modulation about an order of magnitude larger than that of a single red giant with similar rotational period and physical properties. We investigate whether binarity leads to larger magnetic fields when tides lock systems, or if a different spot distribution on single versus close binary stars can explain this fact. For this, we measure the chromospheric emission in the CaII H & K lines of 3130 of the 4465 stars studied in a previous work thanks to the LAMOST survey. We show that red giants in a close-binary configuration with spin-orbit resonance display significantly larger chromospheric emission than single stars, suggesting that tidal locking leads to larger magnetic fields at a fixed rotational period. Beyond bringing interesting new observables to study the evolution of binary systems, this result could be used to distinguish single versus binary red giants in automatic pipelines based on machine learning.

We report the discovery of three new hot Jupiters with the Next Generation Transit Survey (NGTS) as well as updated parameters for HATS-54b, which was independently discovered by NGTS. NGTS-23b, NGTS-24b and NGTS-25b have orbital periods of 4.076, 3.468, and 2.823 days and orbit G-, F- and K-type stars, respectively. NGTS-24 and HATS-54 appear close to transitioning off the main-sequence (if they are not already doing so), and therefore are interesting targets given the observed lack of Hot Jupiters around sub-giant stars. By considering the host star luminosities and the planets' small orbital separations (0.037 - 0.050 au), we find that all four hot Jupiters are above the minimum irradiance threshold for inflation mechanisms to be effective. NGTS-23b has a mass of 0.61 $M_{J}$ and radius of 1.27 $R_{J}$ and is likely inflated. With a radius of 1.21 $R_{J}$ and mass of 0.52 $M_{J}$, NGTS-24b has a radius larger than expected from non-inflated models but its radius is smaller than the predicted radius from current Bayesian inflationary models. Finally, NGTS-25b is intermediate between the inflated and non-inflated cases, having a mass of 0.64 $M_{J}$ and a radius of 1.02 $R_{J}$. The physical processes driving radius inflation remain poorly understood, and by building the sample of hot Jupiters we can aim to identify the additional controlling parameters, such as metallicity and stellar age.

We present the results of a deep study of the neutron star (NS) population in the globular cluster M28 (NGC 6626), using the full 330-ks 2002-2015 ACIS dataset from the Chandra X-ray Observatory and coordinated radio observations taken with the Green Bank Telescope (GBT) in 2015. We investigate the X-ray luminosity (Lx), spectrum, and orbital modulation of the 7 known compact binary millisecond pulsars (MSPs) in the cluster. We report two simultaneous detections of the redback PSR J1824-2452I (M28I) and its X-ray counterpart. We discover a double-peaked X-ray orbital flux modulation in M28I during its pulsar state, centered around pulsar inferior conjunction. We analyze the spectrum of the quiescent neutron star low-mass X-ray binary to constrain its mass and radius. Using both hydrogen and helium NS atmosphere models, we find a NS radius of R = 9.5-11.5 km and R = 13.5 - 16.7 km, respectively, for a neutron star mass of 1.4 Msun. We also search for long-term variability in the 46 brightest X-ray sources and report the discovery of six new variable low luminosity X-ray sources in M28.

We perform three-dimensional magnetohydrodynamic simulations of magnetospheric accretion in a T Tauri star to study the accretion and wind structures in the close vicinity of the star. The gas accreting onto the star consists of the gas from the magnetospheric boundary and the failed disk winds. The accreting gas is commonly found as a multi-column accretion, which is consistent with observations. A significant fraction of the angular momentum of the accreting flows is removed by the magnetic fields of conical disk winds and turbulent failed winds inside and near the magnetosphere. As a result, the accretion torque is significantly reduced compared to the simple estimation based on the mass accretion rate. The stellar spin affects the time variability of the conical disk wind by changing the stability condition of the magnetospheric boundary. However, the time-averaged magnetospheric radius only weakly depends on the stellar spin, which is unlike the prediction of classical theories that the stellar spin controls the magnetospheric radius through the magnetic torque. The ratio of the toroidal to the poloidal field strengths at the magnetospheric boundary, which is a key parameter for the magnetic torque, is also insensitive to the spin; it is rather determined by the disk dynamics. Considering newly found three-dimensional effects, we obtain a scaling relation of the magnetospheric radius very similar to the Ghosh & Lamb relation from the steady angular momentum transport equation.

The distribution of magnetic fields of positive and negative polarities over the surface of the Sun was studied on the basis of synoptic maps NSO Kitt Peak (1978-2016). To emphasize the contribution of weak fields the following transformation of synoptic maps was made: for each synoptic map only magnetic fields with modulus less than 5 G (|B| < 5 G) were left unchanged on each synoptic map while larger or smaller fields were replaced by the corresponding limiting values +5 G or -5 G. Cyclic variations of the magnetic field polarity have been observed associated with two types of magnetic field flows in the photosphere. Rush-to-the-Poles (RTTP) form near the maximum of solar activity and have the same sign as the following sunspots. The lifetime of RTTP is 3 yrs, during which time they drift from latitudes 30 - 40 deg. to the pole, causing the polarity change of the Sun's polar field. We studied another type of variations which has the form of series of flows with individual flows of 0.5-1 yr and with alternating polarity (ripples). Ripples are located in time between two RTTP and drift from the equator to the latitudes of 50 deg. Magnetic field variations were considered in 6 time intervals along the latitudes +33 deg. in the northern and -33 deg. in the southern hemispheres. The time change of the field strength was approximated by the sinusoidal function. The period of variation of ripples was 1.1 yr for the N-hemisphere and 1.3 yr for the S-hemisphere. The amplitude of variation was higher for the time intervals where the polar field had a positive sign. Within the same flow, fields of positive and negative signs developed in anti-phase.

Recent results obtained from leading cosmic ray experiments indicate that simulations using LHC-tuned hadronic interaction models underestimate the number of muons in extensive air showers compared to experimental data. This is the so-called muon deficit problem. Determination of the muon component in the air shower is crucial for inferring the mass of the primary particle, which is a key ingredient in the efforts to pinpoint the sources of ultra-high energy cosmic rays.In this paper, we present a new method to derive the muon signal in detectors, which uses the difference between the total reconstructed (data) and simulated signals is roughly independent of the zenith angle, but depends on the mass of the primary cosmic ray. Such a method offers an opportunity not only to test/calibrate the hadronic interaction models, but also to derive the $\beta$ exponent, which describes an increase of the number of muons in a shower as a function of the energy and mass of the primary cosmic ray. Detailed simulations show a dependence of the $\beta$ exponent on hadronic interaction properties, thus the determination of this parameter is important for understanding the muon deficit problem. We validate the method by using Monte Carlo simulations for the EPOS-LHC and QGSJetII-04 hadronic interaction models, and showing that this method allows us to recover the ratio of the muon signal between EPOS-LHC and QGSJetII-04 and the average $\beta$ exponent for the studied system, within less than a few percent. This is a consequence of the good recovery of the muon signal for each primary included in the analysis.

Non-equilibrium ionization (NEI) is essentially required for astrophysical plasma diagnostics once the plasma status departs from ionization equilibrium assumptions. In this work, we perform fast NEI calculations combined with magnetohydrodynamic (MHD) simulations and analyze the ionization properties of a Petschek-type magnetic reconnection current sheet during solar eruptions. Our simulation reveals Petschek-type slow-mode shocks in the classical Spitzer thermal conduction models and conduction flux-limitation situations. The results show that under-ionized features can be commonly found in shocked reconnection outflows and thermal halo regions outside the shocks. The departure from equilibrium ionization strongly depends on plasma density. In addition, this departure is sensitive to the observable target temperature: the high-temperature iron ions are strongly affected by NEI effects. The under-ionization also affects the synthetic SDO/AIA intensities, which indicates that the reconstructed hot reconnection current sheet structure may be significantly under-estimated either for temperature or apparent width. We also perform the MHD-NEI analysis on the reconnection current sheet in the classical solar flare geometry. Finally, we show the potential reversal between the under-ionized and over-ionized state at the lower tip of reconnection current sheets where the downward outflow collides with closed magnetic loops, which can strongly affect multiple SDO/AIA band ratios along the reconnection current sheet.

The $Euclid$ mission will provide first-of-its-kind coverage in the near-infrared over deep (three fields, $\sim$10-20 square degrees each) and wide ($\sim$10000 square degrees) fields. While the survey is not designed to discover transients, the deep fields will have repeated observations over a two-week span, followed by a gap of roughly six months. In this analysis, we explore how useful the deep field observations will be for measuring properties of Type Ia supernovae (SNe Ia). Using simulations that include $Euclid$'s planned depth, area and cadence in the deep fields, we calculate that more than 3700 SNe between $0.0<z<1.5$ will have at least five $Euclid$ detections around peak with signal-to-noise ratio larger than 3. While on their own, $Euclid$ light curves are not good enough to directly constrain distances, when combined with LSST deep field observations, we find that uncertainties on SN distances are reduced by 20-30% for $z<0.8$ and by 40-50% for $z>0.8$. Furthermore, we predict how well additional $Euclid$ mock data can be used to constrain a key systematic in SN Ia studies - the size of the luminosity 'step' found between SNe hosted in high mass ($>10^{10} M_{\odot}$) and low mass ($>10^{10} M_{\odot}$) galaxies. This measurement has unique information in the rest-frame NIR. We predict that if the step is caused by dust, we will be able to measure its reduction in the NIR compared to optical at the 4$\sigma$ level. We highlight that the LSST and $Euclid$ observing strategies used in this work are still provisional and some level of joint processing is required. Still, these first results are promising, and assuming $Euclid$ begins observations well before the Nancy Roman Space Telescope (Roman), we expect this dataset to be extremely helpful for preparation for Roman itself.

We present a detailed compilation and analysis of the X-ray phase space of low- to intermediate-redshift transients that consolidates observed light curves (and theory where necessary) for a large variety of classes of transient/variable phenomena in the 0.3--10 keV band including gamma-ray burst afterglows, supernovae, supernova shocks interacting with the environment, tidal disruption events and active galactic nuclei, fast blue optical transients, cataclysmic variables, magnetar flares/outbursts and fast radio bursts, cool stellar flares, X-ray binary outbursts, and ultraluminous X-ray sources. Our over-arching goal is to offer a comprehensive resource for the examination of these ephemeral events, extending the X-ray duration-luminosity phase space (DLPS) to show luminosity evolution. We use existing observations (both targeted and serendipitous) to characterize the behavior of various transient/variable populations. Contextualizing transient signals in the larger DLPS serves two primary purposes: to identify areas of interest (i.e., regions in the parameter space where one would expect detections, but in which observations have historically been lacking) and to provide initial qualitative guidance in classifying newly discovered transient signals. We find that while the most luminous (largely extragalactic) and least luminous (largely Galactic) part of the phase space is well-populated at $t > 0.1$ days, intermediate luminosity phenomena (L$_x = 10^{34} - 10^{42}$ erg s$^{-1}$) represent a gap in the phase space. We thus identify L$_x = 10^{34} - 10^{42}$ erg s$^{-1}$ and $t = 10^{-4} - 0.1$ days as a key discovery phase space in transient X-ray astronomy.

We use the MACSIS hydrodynamical simulations to estimate the extent of gas clumping in the intracluster medium of massive galaxy clusters and how it affects the hydrostatic mass bias. By comparing the clumping to the azimuthal scatter in the emission measure, an observational proxy, we find that they both increase with radius and are larger in higher-mass and dynamically perturbed systems. Similar trends are also seen for the azimuthal temperature scatter and non-thermal pressure fraction, both of which correlate with density fluctuations, with these values also increasing with redshift. However, in agreement with recent work, we find only a weak correlation between the clumping, or its proxies, and the hydrostatic mass bias. To reduce the effect of clumping in the projected profiles, we compute the azimuthal median following recent observational studies, and find this reduces the scatter in the bias. We also attempt to correct the cluster masses by using a non-thermal pressure term and find over-corrected mass estimates ($1-b=0.86$ to $1-b=1.15$) from 3D gas profiles but improved mass estimates ($1-b=0.75$ to $1-b=0.85$) from projected gas profiles. We conclude that the cluster-averaged mass bias is minimised from applying a non-thermal pressure correction ($1-b=0.85$) with more modest reductions from selecting clusters that have low clumping ($1-b=0.79$) or are dynamically relaxed ($1-b=0.80$). However, the latter selection is most effective at minimising the scatter for individual objects. Such results can be tested with next generation X-ray missions equipped with high-resolution spectrometers such as Athena.

The Shower Database (SD) of the Meteor Data Center (MDC) had been operating on the basis of stream-naming rules which were too complex and insufficiently precise for 15 years. With a gradual increase in the number of discovered meteor showers, the procedure for submitting new showers to the database and naming them lead to situations that were inconsistent with the fundamental role of the SD - the disambiguation of stream names in the scientific literature. Our aim is to simplify the meteor shower nomenclature rules. We propose a much simpler set of meteor shower nomenclature rules, based on a two-stage approach, similar to those used in the case of asteroids. The first stage applies to a new shower just after its discovery. The second stage concerns the repeatedly observed shower, the existence of which no longer raises any doubts. Our proposed new procedure was approved by a vote of the commission F1 of the IAU in July 2022.

To better understand the ionizing properties of galaxies in the EoR, we investigate deep, rest-frame ultraviolet (UV) spectra of $\simeq 500$ star-forming galaxies at $3 \leq z \leq 5$ selected from the public ESO-VANDELS spectroscopic survey. The absolute ionizing photon escape fraction ($f_{\rm esc}^{\rm abs}$) is derived by combining absorption line measurements with estimates of the UV attenuation. The ionizing production efficiency ($\xi_{ion}$) is calculated by fitting the far-UV (FUV) stellar continuum of the VANDELS galaxies. We find that the $f_{\rm esc}^{\rm abs}$ and $\xi_{ion}$ parameters increase towards low-mass, blue UV-continuum slopes and strong Ly$\alpha$ emitting galaxies, and both are just slightly higher-than-average for the UV-faintest galaxies in the sample. Potential Lyman Continuum Emitters (LCEs) and selected Lyman Alpha Emitters (LAEs) show systematically higher $\xi_{ion}$ ($\log \xi_{ion}$ (Hz\erg) $\approx 25.38, 25.41$) than non-LCEs and non-LAEs ($\log \xi_{ion}$ (Hz\erg) $\approx 25.18, 25.14$) at similar UV magnitudes. This indicates very young underlying stellar populations ($\approx 10~{\rm Myr}$) at relatively low metallicities ($\approx 0.2~{\rm Z_{\odot}}$). The FUV non-ionizing spectra of potential LCEs is characterized by very blue UV slopes ($\leq -2$), enhanced Ly$\alpha$ emission ($\leq -25$A), strong UV nebular lines (e.g., high CIV1550/CIII]1908 $\geq 0.75$ ratios), and weak absorption lines ($\leq 1$A). The latter suggests the existence of low gas-column-density channels in the interstellar medium which enables the escape of ionizing photons. By comparing our VANDELS results against other surveys in the literature, our findings imply that the ionizing budget in the EoR was likely dominated by UV-faint, low-mass and dustless galaxies.

Magnetic helicity is a fundamental constraint in both ideal and resistive magnetohydrodynamics. Measurements of magnetic helicity density on the Sun and other stars are used to interpret the internal behaviour of the dynamo generating the global magnetic field. In this note, we study the behaviour of the global relative magnetic helicity in three self-consistent spherical dynamo solutions of increasing complexity. Magnetic helicity describes the global linkage of the poloidal and toroidal magnetic fields (weighted by magnetic flux), and our results indicate that there are preferred states of this linkage. This leads us to propose that global magnetic reversals are, perhaps, a means of preserving this linkage, since, when only one of the poloidal or toroidal fields reverses, the preferred state of linkage is lost. It is shown that magnetic helicity indicates the onset of reversals and that this signature may be observed at the outer surface.

With the maturation of near-infrared high-resolution spectroscopy, especially when used for precision radial velocity, data reduction has faced unprecedented challenges in terms of how one goes from raw data to calibrated, extracted, and corrected data with required precisions of thousandths of a pixel. Here we present APERO (A PipelinE to Reduce Observations), specifically focused on SPIRou, the near-infrared spectropolarimeter on the Canada--France--Hawaii Telescope (SPectropolarim\`etre InfraROUge, CFHT). In this paper, we give an overview of APERO and detail the reduction procedure for SPIRou. APERO delivers telluric-corrected 2D and 1D spectra as well as polarimetry products. APERO enables precise stable radial velocity measurements on sky (via the LBL algorithm), good to at least ~2 m/s over the current 5-year lifetime of SPIRou.

Standard Model extensions with a strongly coupled dark sector can induce high-multiplicity states of soft quarks. Such final states trigger extremely efficient antinucleus formation. We show that dark matter annihilation or decay into a strongly coupled sector can dramatically enhance the cosmic-ray antinuclei flux -- by six orders of magnitude in the case of ${^4\overline{\text{He}}}$. In this work, we argue that the tentative ${^3\overline{\text{He}}}$ and ${^4\overline{\text{He}}}$ events reported by the AMS-02 collaboration could be the first sign of a strongly coupled dark sector observed in nature.

The paper is a brief review on the existence and basic properties of static, spherically symmetric regular black hole solutions of general relativity, where the source of gravity is represented by nonlinear electromagnetic fields with the Lagrangian function $L$ depending on the single invariant $f = F_{\mu\nu}F^{\mu\nu}$ or on two variables: either $L(f, h)$, where $h = {^*}F_{\mu\nu} F^{\mu\nu}$, where ${^*}F_{\mu\nu}$ is the Hodge dual of $F_{\mu\nu}$, or $L(f, J)$, where $J = F_{\mu\nu}F^{\nu\rho} F_{\rho\sigma} F^{\sigma\mu}$. A number of no-go theorems are discussed, revealing the conditions under which the space-time cannot have a regular center, among which the theorems concerning $L(f,J)$ theories are probably new. These results concern both regular black holes and regular particlelike or starlike objects (solitons) without horizons. Thus, a regular center in solutions with an electric charge $q_e\ne 0$ is only possible with nonlinear electrodynamics (NED) having no Maxwell weak field limit. Regular solutions with $L(f)$ and $L(f, J)$ NED, possessing a correct (Maxwell) weak-field limit, are possible if the system contains only a magnetic charge $q_m \ne 0$. It is shown, however, that in such solutions the causality and unitarity as well as dynamic stability conditions are inevitably violated in a neighborhood of the center. Some particular examples are discussed.

We show that the space-time foam model from string/D-brane theory predicts a scenario in which neutrinos can possess linearly energy dependent speed variation, together with an asymmetry between neutrinos and antineutrinos, indicating the possibility of Lorentz and CPT symmetry violation for neutrinos. Such scenario is supported by a phenomenological conjecture from the possible associations of IceCube ultrahigh-energy neutrino events with the gamma-ray bursts. It is also consistent with the constraints set by the energy-losing decay channels~(e.g., $e^{+}e^{-}$ pair emission, or neutrino splitting) upon superluminal neutrino velocities. We argue that the plausible violations of energy-momentum conservation during decay may be responsible for the stable propagation of these neutrinos, and hence for the evasion of relevant constraints.

Herein, we review the nuclear equations of state (EOSs) %for core-collapse supernova simulations and the constituent nuclei of core-collapse supernovae (CCSNe) and their roles in CCSN simulations. Various nuclei such as deuterons, iron, and extremely neutron-rich nuclei compose in the central engines of CCSNe. The center of a collapsing core is dominated by neutron-rich heavy nuclei prior to the occurrence of core bounce. Their weak interactions significantly affect the neutrino emission and the size of the produced proto-neutron star. After a core bounce, heavy nuclei are dissolved to protons, neutrons, and light nuclei between the expanding shock wave and the newly formed neutron star (NS). Some of the key components in determining the shock-wave dynamics and supernova explosion of outer envelopes are neutrino interactions of nucleons and light nuclei such as deuterons. An EOS provides the relations between thermodynamical properties and the nuclear composition, and is needed to simulate this explosion. Further investigations on uniform and non-uniform nuclear matter are needed to improve the understanding of the mechanism of CCSNe and the properties of supernova nuclei. The knowledge of the EOS for uniform nuclear matter is being continually improved by a combination of microscopic calculations, terrestrial experiments, and NS observations. With reference to various nuclear experiments and current theories, the finite temperature effects on heavy nuclei, formation of light nuclei in dilute nuclear matter, and transition to uniform nuclear matter should be improved in the model of the EOS for non-uniform nuclear matter.

The recent first detection of gravitational waves (GWs) from binary black hole mergers has spurred a renewed interest in possible deviations from General Relativity (GR), since they could be detected in the GWs emitted by such systems. Of particular interest is the ringdown phase of a binary black hole merger, which can be described by linear perturbations about a background stationary black hole solution. These perturbations mainly correspond to a superposition of 'quasi-normal modes' (QNMs), whose frequencies form a discrete set. One expects that modified gravity models could predict QNMs that differ from their GR counterpart: the detailed analysis of the GW signal represents an invaluable window to test GR and to look for specific signatures of modified gravity. The work done in this thesis takes place in the context of scalar-tensor theories of gravity, and more particularly the Degenerate Higher-Order Scalar-Tensor theories. We start by a review of these theories and their properties, and describe a way to reformulate them in a framework with a clear geometrical interpretation. We then study linear perturbations about several existing nonrotating black hole solutions of such theories, and show why the perturbation equations obtained are very hard to decouple in general. When it is possible, in the case of odd parity perturbations, we describe the propagation of waves and relate it to the stability of the underlying spacetime. When it is not, we circumvent the difficulty by making use of an algorithm proposed recently in the mathematical literature that allows us to decouple the equations both at the black hole horizon and at infinity. This allows us to get the asymptotic behaviour of waves on such spacetimes, yielding valuable information that can allow us to rule some of them out. Finally, we use the asymptotic behaviours obtained to compute QNMs numerically.

Alternative theories of gravity predict up to six distinct polarization modes for gravitational waves. Strong gravitational lensing of gravitational waves allows us to probe the polarization content of these signals by effectively increasing the number of observations from the same astrophysical source. The lensing time delays due to the multiple observed lensed images combined with the rotation of the Earth allows for effective non-collocated interferometers to be defined with respect to the source location and hence probe the alternative polarization amplitudes with more observations. To measure these amplitudes, we jointly fit the image observations to a single gravitational wave signal model that takes into account the image magnifications, time delays, and polarization mode amplitudes. We show that for certain systems, we can make a measurement of the relative mode amplitudes for lensed events with two detectable images.

Future gravitational wave interferometers such as LISA, Taiji, DECIGO, and TianQin, will enable precision studies of the environment surrounding black holes. In this paper, we study intermediate and extreme mass ratio binary black hole inspirals, and consider three possible environments surrounding the primary black hole: accretion disks, dark matter spikes, and clouds of ultra-light scalar fields, also known as gravitational atoms. We present a Bayesian analysis of the detectability and measurability of these three environments. Focusing for concreteness on the case of a detection with LISA, we show that the characteristic imprint they leave on the gravitational waveform would allow us to identify the environment that generated the signal, and to accurately reconstruct its model parameters.