Papers for Friday, Mar 17 2023

We present a one-dimensional radiation-hydrodynamic model of a spherically symmetric cloud evolving under the influence of the self-gravity and the feedback from a star cluster forming in its centre. On one hand, the model is simple due to its 1D geometry, on the other hand, the feedback includes the ionising radiation, stellar winds and the radiation pressure acting on gas and dust. The star cluster is formed from the gas flowing into the cloud centre and the feedback parameters are determined from stellar evolution models and the cluster star forming history. The model is compared to the semi-analytic code WARPFIELD implementing similar physical processes and exploring the scenario that the young cluster R136 in the Large Magellanic Cloud was formed due to re-collapse of the shell formed by the previous generation star cluster. A good qualitative agreement is found, however, $3 - 4$ times higher stellar mass is needed to disrupt the cloud in our model, because it takes into account (contrary to WARPFIELD) self-gravity of the cloud surrounding the shell. We use the model to explore star formation in clouds with different mass, radius and density profile measuring their star formation efficiency (SFE), i.e. the fraction of the cloud mass converted to stars. We found that SFE is a function of a single parameter, $\mathrm{log(SFE)} \propto -n_{hm}^{-0.46}$, with $n_{hm}$ being the cloud mean particle density within its half-mass radius. Furthermore, we found that the feedback efficiency, i.e. a fraction of the feedback energy retained by gas, has a nearly constant value $\sim 10^{-3}$.

We derive optimal estimators for the two-, three-, and four-point correlators of statistically isotropic scalar fields defined on the sphere, such as the Cosmic Microwave Background temperature fluctuations, allowing for arbitrary (linear) masking and inpainting schemes. In each case, we give the optimal unwindowed estimator (obtained via a maximum-likelihood prescription, with an associated Fisher deconvolution matrix), and an idealized form, and pay close attention to their efficient computation. For the trispectrum, we include both parity-even and parity-odd contributions, as allowed by symmetry. The estimators can include arbitrary weighting of the data (and remain unbiased), but are shown to be optimal in the limit of inverse-covariance weighting and Gaussian statistics. The normalization of the estimators is computed via Monte Carlo methods, with the rate-limiting steps (involving spherical harmonic transforms) scaling linearly with the number of bins. An accompanying code package, PolyBin, implements these estimators in Python, and we demonstrate the estimators' efficacy via a suite of validation tests.

The promise of gyrochronology is that given a star's rotation period and mass, its age can be inferred. The reality of gyrochronology is complicated by effects other than ordinary magnetized braking that alter stellar rotation periods. In this work, we present an interpolation-based gyrochronology framework that reproduces the time- and mass-dependent spin-down rates implied by the latest open cluster data, while also matching the rate at which the dispersion in initial stellar rotation periods decreases as stars age. We validate our technique for stars with temperatures of 3800-6200 K and ages of 0.08-2.6 gigayears (Gyr), and use it to reexamine the empirical limits of gyrochronology. In line with previous work, we find that the uncertainty floor varies strongly with both stellar mass and age. For Sun-like stars (5800 K), the statistical age uncertainties improve monotonically from $\pm$38% at 0.2 Gyr to $\pm12$% at 2 Gyr, and are caused by the empirical scatter of the cluster rotation sequences combined with the rate of stellar spin-down. For low-mass K-dwarfs (4200 K), the posteriors are highly asymmetric due to stalled spin-down, and $\pm$1$\sigma$ age uncertainties vary non-monotonically between 10% and 50% over the first few gigayears. High-mass K-dwarfs (5000 K) older than 1.5 Gyr yield the most precise ages, with limiting uncertainties currently set by possible changes in the spin-down rate (12% systematic), the calibration of the absolute age scale (8% systematic), and the width of the slow sequence (4% statistical). An open-source implementation, called gyro-interp, is available online at https://github.com/lgbouma/gyro-interp

VISIONS is an ESO public survey of five nearby (d < 500 pc) star-forming molecular cloud complexes that are canonically associated with the constellations of Chamaeleon, Corona Australis, Lupus, Ophiuchus, and Orion. The survey was carried out with VISTA, using VIRCAM, and collected data in the near-infrared passbands J, H, and Ks. With a total on-sky exposure time of 49.4 h VISIONS covers an area of 650 deg$^2$, and it was designed to build an infrared legacy archive similar to that of 2MASS. Taking place between April 2017 and March 2022, the observations yielded approximately 1.15 million images, which comprise 19 TB of raw data. The observations are grouped into three different subsurveys: The wide subsurvey comprises shallow, large-scale observations and has visited the star-forming complexes six times over the course of its execution. The deep subsurvey of dedicated high-sensitivity observations has collected data on the areas with the largest amounts of dust extinction. The control subsurvey includes observations of areas of low-to-negligible dust extinction. Using this strategy, the VISIONS survey offers multi-epoch position measurements, is able to access deeply embedded objects, and provides a baseline for statistical comparisons and sample completeness. In particular, VISIONS is designed to measure the proper motions of point sources with a precision of 1 mas/yr or better, when complemented with data from VHS. Hence, VISIONS can provide proper motions for sources inaccessible to Gaia. VISIONS will enable addressing a range of topics, including the 3D distribution and motion of embedded stars and the nearby interstellar medium, the identification and characterization of young stellar objects, the formation and evolution of embedded stellar clusters and their initial mass function, as well as the characteristics of interstellar dust and the reddening law.

The galaxy group NGC6338 is one of the most violent group-group mergers known to date. While the central dominant galaxies rush at each other at 1400km/s along the line of sight, with dramatic gas heating and shock fronts detected, the central gas in the BCGs remains cool. There are also indications of feedback from active galactic nuclei (AGNs), and neither subcluster core has been disrupted. With our deep radio uGMRT data at 383MHz and 650MHz we clearly detect a set of large, old lobes in the southern BCG coinciding with the X-ray cavities, while the northern, and smaller BCG appears slightly extended in the radio. The southern BCG also hosts a smaller younger set of lobes, perpendicular to the larger lobes, but also coinciding with the inner X-ray cavities, and matching the jet direction in the parsec-resolution VLBA image. Our spectral analysis confirms the history of two feedback cycles. The high radio frequency analysis classifies the compact source in the southern BCG with a powerlaw, while ruling out a significant contribution from accretion. The radio lightcurve over 3 decades shows a change about 10 years ago, which might be related to ongoing feedback in the core. The southern BCG in the NGC6338 merger remains another prominent case where the direction of jet-mode feedback between two cycles changed dramatically.

The conditions under which accreting neutron stars launch radio-emitting jets and/or outflows are still poorly understood. The ultracompact X-ray binary X1850--087, located in the globular cluster NGC 6712, is a persistent atoll-type X-ray source that has previously shown unusual radio continuum variability. Here we present the results of a pilot radio monitoring program of X1850--087 undertaken with the Karl G. Jansky Very Large Array, with simultaneous or quasi-simultaneous Swift/XRT data obtained at each epoch. The binary is clearly detected in the radio in two of the six new epochs. When combined with previous data, these results suggest that X1850--087 shows radio emission at a slightly elevated hard state X-ray luminosity of L_X > 2x10^36 erg/s, but no radio emission in its baseline hard state L_X ~10^36 erg/s. No clear X-ray spectral changes are associated with this factor of > 10 radio variability. At all detected epochs X1850--087 has a flat-to-inverted radio spectral index, more consistent with the partially absorbed optically thick synchrotron of a compact jet rather than the evolving optically thick to thin emission associated with transient expanding synchrotron-emitting ejecta. If the radio emission in X1850--087 is indeed due to a compact jet, then it is plausibly being launched and quenched in the hard state on timescales as short as a few days. Future radio monitoring of X1850--087 could help elucidate the conditions under which compact jets are produced around hard state accreting neutron stars.

We present simultaneous optical and near-infrared spectra obtained during the 2021 outburst of the black hole transient 4U 1543-47. The X-ray hardness-intensity diagram and the comparison with similar systems reveal a luminous outburst, probably reaching the Eddington luminosity, as well as a long-lasting excursion to the so-called ultra-luminous state. VLT/X-shooter spectra were taken in two epochs 14 days apart during the early and brightest part of the outburst, while the source was in this ultra-luminous accretion state. The data show strong H and HeI emission lines, as well as high-excitation HeII and OIII transitions. Most lines are single-peaked in both spectra, except for the OIII lines that exhibit evident double-peaked profiles during the second epoch. The Balmer lines are embedded in broad absorption wings that we believe are mainly produced by the contribution of the A2V donor to the optical flux, which we estimate to be in the range of 11 to 14 per cent in the $r$ band during our observations. Although no conspicuous outflow features are found, we observe some wind-related line profiles, particularly in the near-infrared. Such lines include broad emission line wings and skewed red profiles, suggesting the presence of a cold (i.e. low ionisation) outflow with similar observational properties to those found in other low-inclination black hole transients.

We use spectroscopic data for ${\sim}6,000$ Red Giant Branch (RGB) stars in the Small Magellanic Cloud (SMC), together with proper motion data from \textit{Gaia} Early Data Release 3 (EDR3), to build a mass model of the SMC. We test our Jeans mass modelling method (\textsc{Binulator}+\textsc{GravSphere}) on mock data for an SMC-like dwarf undergoing severe tidal disruption, showing that we are able to successfully remove tidally unbound interlopers, recovering the Dark Matter density and stellar velocity anisotropy profiles within our 95\% confidence intervals. We then apply our method to real SMC data, finding that the stars of the cleaned sample are isotropic at all radii (at 95\% confidence), and that the inner Dark Matter density profile is dense, $\rho_{\rm DM}(150\,{\rm pc}) = 2.81_{-1.07}^{+0.72}\times 10^8 M_{\odot} \rm kpc^{-3} $, consistent with a $\Lambda$ Cold Dark Matter ($\Lambda$CDM) cusp at least down to 400\,pc from the SMC's centre. Our model gives a new estimate of the SMC's total mass within 3\,kpc ($M_{\rm tot} \leq 3\,{\rm kpc})$ of $2.34\pm0.46 \times 10^9 M_{\odot}$. We also derive an astrophysical \textquote{$J$-factor} of $19.22\pm0.14$\, GeV$^2$\,cm$^{-5}$ and a \textquote{$D$-factor} of $18.80\pm0.03$\, GeV$^2$\,cm$^{-5}$, making the SMC a promising target for Dark Matter annihilation and decay searches. Finally, we combine our findings with literature measurements to test models in which Dark Matter is \textquote{heated up} by baryonic effects. We find good qualitative agreement with the Di Cintio et al. 2014 model but we deviate from the Lazar et al. 2020 model at high $M_*/M_{200} > 10^{-2}$. We provide a new, analytic, density profile that reproduces Dark Matter heating behaviour over the range $10^{-5} < M_*/M_{200} < 10^{-1}$.

The VISIONS public survey provides large-scale, multiepoch imaging of five nearby star-forming regions at subarcsecond resolution in the near-infrared. All data collected within the program and provided by the European Southern Observatory (ESO) science archive are processed with a custom end-to-end pipeline infrastructure to provide science-ready images and source catalogs. The data reduction environment has been specifically developed for the purpose of mitigating several shortcomings of the bona fide data products processed with software provided by the Cambridge Astronomical Survey Unit (CASU), such as spatially variable astrometric and photometric biases of up to 100 mas and 0.1 mag, respectively. At the same time, the resolution of coadded images is up to 20% higher compared to the same products from the CASU processing environment. Most pipeline modules are written in Python and make extensive use of C extension libraries for numeric computations, thereby simultaneously providing accessibility, robustness, and high performance. The astrometric calibration is performed relative to the Gaia reference frame, and fluxes are calibrated with respect to the source magnitudes provided in the Two Micron All Sky Survey (2MASS). For bright sources, absolute astrometric errors are typically on the order of 10 to 15 mas and fluxes are determined with subpercent precision. Moreover, the calibration with respect to 2MASS photometry is largely free of color terms. The pipeline produces data that are compliant with the ESO Phase 3 regulations and furthermore provides curated source catalogs that are structured similarly to those provided by the 2MASS survey.

The bulk of X-ray spectroscopic studies of active galactic nuclei (AGN) are focused on local ($z < 0.1$) sources with low-to-moderate ($< 0.3$) Eddington ratio ($\lambda_\mathrm{Edd}$). It is then mandatory to overcome this limitation and improve our understanding of highly accreting AGN. In this work we present the preliminary results from the analysis of a sample of $\sim70$ high-$\lambda_\mathrm{Edd}$ radio-quiet AGN at $0.06 \leq z \leq 3.3$, based on the 10th release of the XMM-Newton serendipitous source catalogue, that we named as XMM-Newton High-Eddington Serendipitous AGN Sample (X-HESS). Almost $\sim35\%$ of the X-HESS AGN have multi-epoch archival observations and $\sim70\%$ of the sources can rely on simultaneous OM optical data. First results reveal sources showing signatures of ultra-fast outflows and remarkable long- and short-term X-ray flux variations. Indeed in J095847.88+690532.7 ($z \sim 1.3$), one of the most densely monitored objects hosting a $\sim$$10^9\,M_\odot$ supermassive black hole, we discovered a variation of the soft X-ray flux by a factor of > 2 over approximately one week (rest-frame). Large variations in the power-law continuum photon index $\Gamma$ are also observed, questioning expectations from previously reported $\Gamma - \lambda_\mathrm{Edd}$ relations, for which $\Gamma \geq 2$ would be a ubiquitous hallmark of AGN with $\lambda_\mathrm{Edd} \sim 1$.

PBC J2333.9-2343 is a giant radio galaxy at z = 0.047 with a bright central core associated to a blazar nucleus. If the nuclear blazar jet is a new phase of the jet activity, then the small orientation angle suggest a dramatic change of the jet direction. We present observations obtained between September 2018 and January 2019 (cadence larger than three days) with Effeslberg, SMARTS-1.3m, ZTF, ATLAS, Swift, and Fermi-LAT, and between April-July 2019 (daily cadence) with SMARTS-1.3m and ATLAS. Large (>2x) flux increases are observed on timescales shorter than a month, which are interpreted as flaring events. The cross correlation between the SMARTS-1.3m monitoring in the NIR and optical shows that these data do not show significant time lag within the measured errors. A comparison of the optical variability properties between non-blazars and blazars AGN shows that PBC J2333.9-2343 has properties more comparable to the latter. The SED of the nucleus shows two peaks, that were fitted with a one zone leptonic model. Our data and modelling shows that the high energy peak is dominated by External Compton from the dusty torus with mild contribution from Inverse Compton from the jet. The derived jet angle of 3 degrees is also typical of a blazar. Therefore, we confirm the presence of a blazar-like core in the center of this giant radio galaxy, likely a Flat Spectrum Radio Quasar with peculiar properties.

We revisit the cosmology of neutrino self-interaction and use the latest cosmic microwave background data from the Atacama Cosmology Telescope (ACT) and the Planck experiment to constrain the interaction strength. In both flavor-universal and nonuniversal coupling scenarios, we find that the ACT data prefers strong neutrino self-interaction that delays neutrino free streaming until just before the matter-radiation equality. When combined with the Planck 2018 data, the preference for strong interaction decreases due to the Planck polarization data. For the combined dataset, the flavor-specific interaction still provides a better fit to the CMB data than $\Lambda$CDM. This trend persists even when neutrino mass is taken into account and extra radiation is added. We also study the prospect of constraining such strong interaction by future terrestrial and space telescopes, and find that the upcoming CMB-S4 experiment will improve the upper limit on neutrino self-interaction by about a factor of three.

Fuzzy Dark Matter (FDM), consisting of ultralight bosons, is an intriguing alternative to Cold Dark Matter. Numerical simulations solving the Schr\"odinger-Poisson (SP) equation, which governs FDM dynamics, show that FDM halos consist of a central solitonic core (representing the ground state of the SP equation), surrounded by a large envelope of excited states. Wave interference gives rise to order unity density fluctuations throughout the envelope and causes the soliton to undergo density oscillations and execute a confined random walk in the central region of the halo. The resulting gravitational potential perturbations are an efficient source of dynamical heating. Using high-resolution numerical simulations of a $6.6 \times 10^{9} \rm M_{\odot}$ FDM halo with boson mass, $m_{\rm b}=8 \times 10^{-23} \ \rm eV$, we investigate the impact of this dynamical heating on the structure and kinematics of spheroidal dwarf galaxies of a fixed mass but different initial sizes and ellipticities. The galaxies are set up in equilibrium in the time-and-azimuthally averaged halo potential and evolved for $10 \ \rm Gyr$ in the live FDM halo. We find that they continuously increase their sizes and central velocity dispersions. In addition, their kinematic structures become strongly radially anisotropic, especially in the outskirts. Dynamical heating also causes initially ellipsoidal galaxies to become more spherical over time from the inside out and gives rise to distorted, non-concentric isodensity contours. These tell-tale characteristics of dynamical heating of dwarf galaxies in FDM halos can potentially be used to constrain the boson mass.

Collisionless shock waves in supernova remnants and the solar wind heat electrons less effectively than they heat ions, as is predicted by kinetic simulations. However, the values of T$_e$/T$_p$ inferred from the H alpha profiles of supernova remnant shocks behave differently as a function of Mach number or Alfv\'{e}n Mach number than what is measured in the solar wind or predicted by simulations. Here we determine T$_e$/T$_p$ for supernova remnant shocks using H alpha profiles, shock speeds from proper motions, and electron temperatures from X-ray spectra. We also improve the estimates of sound speed and Alfv\'{e}n speed used to determine Mach numbers. We find that the H alpha determinations are robust and that the discrepancies among supernova remnant shocks, solar wind shocks and computer-simulated shocks remain. We discuss some possible contributing factors, including shock precursors, turbulence and varying preshock conditions.

Low-ionization Broad Absorption Line QSOs (LoBALs) are suspected to be merging systems in which extreme, AGN-driven outflows have been triggered. Whether or not LoBALs are uniquely associated with mergers, however, has yet to be established. To characterize the morphologies of LoBALs, we present the first high-resolution morphological analysis of a volume-limited sample of 22 SDSS-selected LoBALs at 0.5 < z < 0.6 from Hubble Space Telescope Wide Field Camera 3 observations. Host galaxies are resolved in 86% of the systems in F125W, which is sensitive to old stellar populations, while only 18% are detected in F475W, which traces young, unobscured stellar populations. Signs of recent or ongoing tidal interaction are present in 45-64% of the hosts, including double nuclei, tidal tails, bridges, plumes, shells, and extended debris. Ongoing interaction with a companion is apparent in 27-41% of the LoBALs, with as much as 1/3 of the sample representing late-stage mergers at projected nuclear separations <10 kpc. Detailed surface brightness modeling indicates that 41% of the hosts are bulge-dominated while only 18% are disks. We discuss trends in various properties as a function of merger stage and parametric morphology. Notably, mergers are associated with slower, dustier winds than those seen in undisturbed/unresolved hosts. Our results favor an evolutionary scenario in which quasar-level accretion during various merger stages is associated with the observed outflows in low-z LoBALs. We discuss differences between LoBALs and FeLoBALs and show that selection via the traditional Balnicity index would have excluded all but one of the mergers.

We systematically investigate the claim by Vandorou et al. (2023) to have detected the host star of the low mass-ratio ($q<10^{-4}$) microlensing planet OGLE-2016-BLG-1195Lb, via Keck adaptive optics (AO) measurements $\Delta t=4.12\,$yr after the peak of the event ($t_0$). If correct, this measurement would contradict the microlens parallax measurement derived from Spitzer observations in solar orbit taken near $t_0$. We show that this host identification would be in $4\,\sigma$ conflict with the original ground-based lens-source relative proper-motion measurements. By contrast, Gould (2022) estimated a probability $p=10\%$ that the ``other star'' resolved by single-epoch late-time AO would be a companion to the host or the microlensed source, which is much more probable than a 4$\,\sigma$ statistical fluctuation. In addition, independent of this proper-motion discrepancy, the kinematics of this host-identification are substantially less probable than those of the Spitzer solution. Hence, this identification should not be accepted, pending additional observations that would either confirm or contradict it, which could be taken in 2023. Motivated by this tension, we present two additional investigations. We explore the possibility that Vandorou et al. (2023) identified the wrong ``star'' (or stellar asterism) on which to conduct their analysis. We find that astrometry of KMT and Keck images favors a star (or asterism) lying about 175 mas northwest of the one that they chose. We also present event parameters from a combined fit to all survey data, which yields, in particular, a more precise mass ratio, $q=(4.6\pm 0.4)\times 10^{-5}$. Finally, we discuss the broader implications of minimizing such false positives for the first measurement of the planet mass function, which will become possible when AO on next-generation telescopes are applied to microlensing planets.

All single stars that are born with masses up to 8.5 - 10 $M_\odot$ will end their lives as a white dwarf (WD) star. In this evolutionary stage, WDs enter the cooling sequence, where the stars radiate away their thermal energy, and are basically cooling. As these stars cool, they reach temperatures and conditions that cause the stars to pulsate. Using differential photometry to produce light curves, we can determine the observed periods of pulsation from the WD. We used the White Dwarf Evolution Code (WDEC) to calculate a grid of over one million models with various temperature, stellar mass and mass of helium and hydrogen layers, and calculated their theoretical pulsation periods. In this paper, we describe our approach to WD asteroseismology using WDEC models and we present seismological studies for 29 observed DAVs in the Kepler and K2 datasets, 25 of which have never been analyzed using these observations, and 19 of which have never been seismically analyzed in any capacity before. Learning about the internal structure of WDs place important constraints on the WD cooling sequence and our overall understanding of stellar evolution for low mass stars.

The 'Great Dimming' of the prototypical red supergiant Betelgeuse, which occurred between December 2019 and April 2020, gives us unprecedented insight into the processes occurring on the stellar surface and in the inner wind of this type of star. In particular it may bring further understanding of their dust nucleation and mass loss processes. Here, we present and analyse VLTI/MATISSE observations in the N-band (8 - 13 $\mu$m) taken near the brightness minimum in order to assess the status of the dusty circumstellar environment. We explore the compatibility of a dust clump obscuring the star with our mid-infrared interferometric observations using continuum 3D radiative transfer modelling, and probe the effect of adding multiple clumps close to the star on the observables. We also test the viability of a large cool spot on the stellar surface without dust present in the ambient medium. Using the visibility data, we derive a uniform disk diameter of 59.02 $\pm$ 0.64 mas in the spectral range 8 to 8.75 $\mu$m. We find that both the dust clump and the cool spot models are compatible with the data. Further to this, we note that the extinction and emission of our localised dust clump in the line of sight of the star, directly compensate each other making the clump undetectable in the spectral energy distribution and visibilities. The lack of infrared brightening during the 'Great Dimming' therefore does not exclude extinction due to a dust clump as one of the possible mechanisms. The visibilities can be reproduced by a spherical wind with dust condensing at 13 stellar radii and a dust mass-loss rate of (2.1 - 4.9) $\times$ 10$^{-10}$ $\mathit{M}_{\odot} {\rm yr}^{-1}$, however, in order to reproduce the complexity of the observed closure phases, additional surface features or dust clumps would be needed.

Despite the fact that jets from black holes were first understood to exist over 40 years ago, we are still in ignorance about many primary aspects of these systems -- including the radiation mechanism at high energies, the particle makeup of the jets, and how particles are accelerated, possibly to energies as high as 100 TeV and hundreds of kpc from the central engine. We focus in particular on the discovery (and mystery) of strong X-ray emission from radio jets on kpc-scales, enabled by the unequaled high resolution of the \emph{Chandra} X-ray observatory. We review the main evidence for and against the viable models to explain this X-ray emission over the last 20 years. Finally, we present results of a recent study on the X-ray variability of kpc-scale jets, where we find evidence that between 30-100\% of the X-ray jet population is variable at the tens-of-percent level. The short ($\sim$years) variability timescale is incompatible with the IC/CMB model for the X-rays and implies extremely small structures embedded within the kpc-scale jet, and thus requires a reconsideration of many assumptions about jet structure and dynamics.

In order to probe modifications of gravity at cosmological scales, one needs accurate theoretical predictions. N-body simulations are required to explore the non-linear regime of structure formation but are very time consuming. In this work, we build an emulator, dubbed e-MANTIS, that performs an accurate and fast interpolation between the predictions of a given set of cosmological simulations, in $f(R)$ modified gravity, run with ECOSMOG. We sample a wide 3D parameter space given by the current background scalar field value $10^{-7} < \left|f_{R_0} \right| < 10^{-4}$, matter density $0.24<\Omega_\mathrm{m}<0.39$, and primordial power spectrum normalisation $0.6<\sigma_8<1.0$, with 110 points sampled from a Latin Hypercube. For each model we perform pairs of $f(R)$CDM and $\Lambda$CDM simulations covering an effective volume of $\left(560 \, h^{-1}\mathrm{Mpc}\right)^3$ with a mass resolution of $\sim 2 \times 10^{10} h^{-1} M_\odot$. We compute the matter power spectrum boost due to $f(R)$ gravity $B(k)=P_{f(R)}(k)/P_{\Lambda\mathrm{CDM}}(k)$ and build an emulator using a Gaussian Process Regression method. The boost is mostly independent of $h$, $n_{s}$, and $\Omega_{b}$, which reduces the dimensionality of the relevant cosmological parameter space. Additionally, it is much more robust against statistical and systematic errors than the raw power spectrum, thus strongly reducing our computational needs. The resulting emulator has a maximum error of $3\%$ across the whole cosmological parameter space, for scales $0.03 \ h\mathrm{Mpc}^{-1} < k < 7 \ h\mathrm{Mpc}^{-1}$, and redshifts $0 < z < 2$, while in most cases the accuracy is better than $1\%$. Such an emulator could be used to constrain $f(R)$ gravity with weak lensing analyses.

Local primordial non-Gaussianity (LPNG) is predicted by many non-minimal models of inflation, and creates a scale-dependent contribution to the power spectrum of large-scale structure (LSS) tracers, whose amplitude is characterized by $b_{\phi}$. Knowledge of $b_{\phi}$ for the observed tracer population is therefore crucial for learning about inflation from LSS. Recently, it has been shown that the relationship between linear bias $b_1$ and $b_{\phi}$ for simulated halos exhibits significant secondary dependence on halo concentration. We leverage this fact to forecast multi-tracer constraints on $f_{NL}^{\mathrm{loc}}$. We train a machine learning model on observable properties of simulated Illustris-TNG galaxies to predict $b_{\phi}$ for samples constructed to approximate DESI emission line galaxies (ELGs) and luminous red galaxies (LRGs). We find $\sigma(f_{NL}^{\mathrm{loc}}) = 2.3$, and $\sigma(f_{NL}^{\mathrm{loc}}) = 3.7$, respectively. These forecasted errors are roughly factors of 3, and 35\% improvements over the single-tracer case for each sample, respectively. When considering both ELGs and LRGs in their overlap region, we forecast $\sigma(f_{NL}^{\mathrm{loc}}) = 1.5$ is attainable with our learned model, more than a factor of 3 improvement over the single-tracer case, while the ideal split by $b_{\phi}$ could reach $\sigma(f_{NL}^{\mathrm{loc}}) <1$. We also perform multi-tracer forecasts for upcoming spectroscopic surveys targeting LPNG (MegaMapper, SPHEREx) and show that splitting tracer samples by $b_{\phi}$ can lead to an order-of-magnitude reduction in projected $\sigma(f_{NL}^{\mathrm{loc}})$ for these surveys.

We report the discovery of an accreting supermassive black hole at z=8.679, in CEERS_1019, a galaxy previously discovered via a Ly$\alpha$-break by Hubble and with a Ly$\alpha$ redshift from Keck. As part of the Cosmic Evolution Early Release Science (CEERS) survey, we observed this source with JWST/NIRSpec spectroscopy, MIRI and NIRCam imaging, and NIRCam/WFSS slitless spectroscopy. The NIRSpec spectra uncover many emission lines, and the strong [O III] emission line confirms the ground-based Ly$\alpha$ redshift. We detect a significant broad (FWHM~1200 km/s) component in the H$\beta$ emission line, which we conclude originates in the broad-line region of an active galactic nucleus (AGN), as the lack of a broad component in the forbidden lines rejects an outflow origin. This hypothesis is supported by the presence of high-ionization lines, as well as a spatial point-source component embedded within a smoother surface brightness profile. The mass of the black hole is log($M_{BH}/M_{\odot})=6.95{\pm}0.37$, and we estimate that it is accreting at 1.2 ($\pm$0.5) x the Eddington limit. The 1-8 $\mu$m photometric spectral energy distribution (SED) from NIRCam and MIRI shows a continuum dominated by starlight and constrains the host galaxy to be massive (log M/M$_{\odot}$~9.5) and highly star-forming (SFR~30 M$_{\odot}$ yr$^{-1}$). Ratios of the strong emission lines show that the gas in this galaxy is metal-poor (Z/Z$_{\odot}$~0.1), dense (n$_{e}$~10$^{3}$ cm$^{-3}$), and highly ionized (log U~-2.1), consistent with the general galaxy population observed with JWST at high redshifts. We use this presently highest-redshift AGN discovery to place constraints on black hole seeding models and find that a combination of either super-Eddington accretion from stellar seeds or Eddington accretion from massive black hole seeds is required to form this object by the observed epoch.

The short-lived radionuclides (SLRs) have a half-life $\leq$ 100 Myr. The $\gamma$-ray observations and excess abundance of their daughter nuclides in various meteoritic phases confirm the existence of SLRs in the Galaxy and early solar system (ESS), respectively. In this work, we have developed Galactic Chemical Evolution (GCE) models for SLRs, $^{26}$Al, and $^{60}$Fe along with $^{36}$Cl, $^{41}$Ca, and $^{53}$Mn. These models predict the temporal and spatial evolution of SLR abundance trends in the Galaxy from 2-18 kpc. The abundance of two SLRs, $^{26}$Al, and $^{60}$Fe, are investigated further, as their $\gamma$-ray observations are available for comparison with the model predictions. The predictions for the abundance per unit area for each ring decrease from the inner to outer regions of the Galaxy.The GCE predictions for the total mass of alive $^{26}$Al, and $^{60}$Fe in 2-18 kpc of the Galaxy at present time are 0.2 M$_\odot$ and 0.08 M$_\odot$, respectively.

Abstract. We discuss recent results on the clustering, composition and distribution of Ultra-High Energy Cosmic Rays (UHECR) in the sky; from the energy of several tens of EeV in the dipole anisotropy, up to the highest energy of a few narrow clusters, those of Hot Spots. Following the early UHECR composition records deviations from proton, we noted that the UHECR events above 40 EeV can be made not just by any light or heavy nuclei, but mainly by the lightest ones as He,D, Li,Be. The remarkable Virgo absence and the few localized nearby extragalactic sources, such as CenA, NGC 253 and M82, are naturally understood: lightest UHECR nuclei cannot reach us from the Virgo distance of twenty Mpc, due to their nuclei fragility above a few Mpc distances. Their deflection and smearing in wide hot spots is better tuned to the lighter nuclei than to the preferred proton or heavy nuclei candidate courier. We note that these lightest nuclei still suffer of a partial photodistruction even from such close sources. Therefore, their distruption in fragments, within few tens EeV multiplet chain of events, have been expected and later on observed by Auger collaboration, nearly a decade ago. These multiplet presences, strongly correlate with the same CenA, NGC253 sources. The statistical weight of such correlation is reminded. We conclude that the same role of NGC 253 clustering at lower energies could also feed the Auger dipole anisotropy at lower energy ranges, integrated by nearest Vela, Crab, LMC and Cas A contributes. In our present UHECR model, based on lightest nuclei in local volumes of a few Mpcs, closest AGN, Star-Burst or very close SNR are superimposing their signals, frozen in different epochs, distances and directions, feeding small and wide anisotropy. Possible tests to confirm, or untangle the current model from alternative ones, are suggested and updated.

The elemental carbon-to-oxygen ratio (C/O) in the atmosphere of a giant planet is a promising diagnostic of that planet's formation history in a protoplanetary disk. Alongside efforts in the exoplanet community to measure C/O in planetary atmospheres, observational and theoretical studies of disks are increasingly focused on understanding how the gas-phase C/O varies both with radial location and between disks. This is mostly tied to the icelines of major volatile carriers such as CO and H2O. Using ALMA observations of CS and SO, we have unearthed evidence for an entirely novel type of C/O variation in the protoplanetary disk around HD 100546: an azimuthal variation from a typical, oxygen-dominated ratio (C/O=0.5) to a carbon-dominated ratio (C/O>1.0). We show that the spatial distribution and peculiar line kinematics of both CS and SO molecules can be well-explained by azimuthal variations in the C/O ratio. We propose a shadowing mechanism that could lead to such a chemical dichotomy. Our results imply that tracing the formation history of giant exoplanets using their atmospheric C/O ratios will need to take into account time-dependent azimuthal C/O variations in a planet's accretion zone.

A new golden age in astronomy is upon us, dominated by data. Large astronomical surveys are broadcasting unprecedented rates of information, demanding machine learning as a critical component in modern scientific pipelines to handle the deluge of data. The upcoming Legacy Survey of Space and Time (LSST) of the Vera C. Rubin Observatory will raise the big-data bar for time-domain astronomy, with an expected 10 million alerts per-night, and generating many petabytes of data over the lifetime of the survey. Fast and efficient classification algorithms that can operate in real-time, yet robustly and accurately, are needed for time-critical events where additional resources can be sought for follow-up analyses. In order to handle such data, state-of-the-art deep learning architectures coupled with tools that leverage modern hardware accelerators are essential. We showcase how the use of modern deep compression methods can achieve a $18\times$ reduction in model size, whilst preserving classification performance. We also show that in addition to the deep compression techniques, careful choice of file formats can improve inference latency, and thereby throughput of alerts, on the order of $8\times$ for local processing, and $5\times$ in a live production setting. To test this in a live setting, we deploy this optimised version of the original time-series transformer, t2, into the community alert broking system of FINK on real Zwicky Transient Facility (ZTF) alert data, and compare throughput performance with other science modules that exist in FINK. The results shown herein emphasise the time-series transformer's suitability for real-time classification at LSST scale, and beyond, and introduce deep model compression as a fundamental tool for improving deploy-ability and scalable inference of deep learning models for transient classification.

Accreting compact objects are crucial to understand several important astrophysical phenomena such as Type Ia supernovae, gravitational waves, or X-ray and $\gamma$-ray bursts. In addition, they are natural laboratories to infer fundamental properties of stars, to investigate high-energy phenomena and accretion processes, to test theories of stellar and binary evolution, to explore interactions between high-density plasma and very strong magnetic fields, to examine the interplay between binary evolution and dynamical interactions (in the case they belong to dense star clusters), and they can even be used as a probe for the assembling process of galaxies over cosmic time-scales. Despite the fundamental importance of accreting compact objects for astrophysics and recent progress with the comprehension of these fascinating objects, we still do not fully understand how they form and evolve. In this chapter, we will review the current theoretical status of our knowledge on these objects, and will discuss standing problems and potential solutions to them.

Filaments are ubiquitous in the interstellar medium (ISM), yet their formation and evolution remains the topic of intense debate. In order to obtain a more comprehensive view of the 3D morphology and evolution of the Musca filament, we model the C$^{18}$O(2-1) emission along the filament crest with several large-scale velocity field structures. This indicates that Musca is well described by a 3D curved cylindrical filament with longitudinal mass inflow to the center of the filament unless the filament is a transient structure with a lifetime $\lesssim$~0.1 Myr. Gravitational longitudinal collapse models of filaments appear unable to explain the observed velocity field. To better understand these kinematics, we further analyze a map of the C$^{18}$O(2-1) velocity field at the location of SOFIA HAWC+ dust polarization observations that trace the magnetic field in the filament. This unveils an organized magnetic field that is oriented roughly perpendicular to the filament crest. Although the velocity field is also organized, it progressively changes its orientation by more than 90$^{o}$ when laterally crossing the filament crest and thus appears disconnected from the magnetic field in the filament. This strong lateral change of the velocity field over the filament remains unexplained and might be associated with important longitudinal motion in the filament that can be associated to the large-scale kinematics along the filament.

The properties of the matter density field in the initial conditions have a decisive impact on the features of the large-scale structure of the Universe as observed today. These need to be studied via $N$-body simulations, which are imperative to analyze high density collapsed regions into dark matter halos. In this paper, we train Machine Learning algorithms with information from N -body simulations to infer two properties: dark matter particle halo classification that leads to halo formation prediction with the characteristics of the matter density field traced back to the initial conditions, and dark matter halo formation by calculating the Halo Mass Function (HMF), which offers the number density of dark matter halos with a given threshold. We map the initial conditions of the matter density field into classification labels of dark matter halo structures. The Halo Mass Function of the simulations is calculated and reconstructed with theoretical methods as well as our trained algorithms. We test several Machine Learning techniques where we could find that the Random Forest and Neural Networks proved to be the better performing tools to classify dark matter particles in cosmological simulations. We also show that that it is not compulsory to use a high amount of data to train the algorithms in order to reconstruct the HMF, giving us a very good fitting function for both simulation and theoretical results.

The first multiband photometric solutions of the short-period V Gru eclipsing binary from the southern hemisphere is presented in this study. Light curves of the system were observed through BVI filters at the Congarinni Observatory in Australia for 15 nights. In addition to the new ground-based data, we also used the TESS observations in two sectors. We analyzed the light curves of the system using the PHysics Of Eclipsing BinariEs (PHOEBE) 2.4.7 version code to achieve the best accordance with the photometric observations. The solutions suggest that V Gru is a near-contact binary system with q=1.302(81) mass ratio, f1=0.010(23), f2=-0.0.009(21), and i=73.45(38). We considered the two hot spots on the hotter and cooler components for the light curve analysis. We extracted the minima times from the light curves based on the Markov Chain Monte Carlo (MCMC) approach. Using our new light curves, TESS, and additional literature minima, we computed the ephemeris of V Gru. The system's eclipse timing variation trend was determined using the MCMC method. This system is a good and challenging case for future studies.

LDN1415-IRS, a low-mass young stellar object (YSO) went into an outburst between 2001 and 2006, illuminating a surrounding nebula, LDN1415-Neb. LDN1415-Neb was found to have brightened by I=3.77 mag by April 2006. The optical light curve covering $\sim$ 15.5 years, starting from October 2006 to January 2022, is presented in this study. The initial optical spectrum indicated the presence of winds in the system but the subsequent spectra have no wind indicators. The declining light curve and the absence of the P-Cygni profile in later epoch spectra indicate that the star and nebula system is retrieving back from its outburst state. Two Herbig-Haro objects (HHOs) are positioned linearly with respect to the optical brightness peak of the nebula, probably indicating the circumstellar disk being viewed edge-on. Our recent deep near-infrared (NIR) imaging using TANSPEC has revealed the presence of a nearby star-like source, south of the LDN1415-IRS, at an angular distance of $\sim$ 5.4 arcsec.

We report a trigonometric parallax measurement of 22 GHz water masers in the massive star-forming region G034.43+0.24 as part of the Bar and Spiral Structure Legacy (BeSSeL) Survey using the Very Long Baseline Array. The parallax is 0.330$\pm$50.018 mas, corresponding to a distance of $3.03^{+0.17}_{-0.16}$ kpc. This locates G034.43+0.24 near the inner edge of the Sagittarius spiral arm and at one end of a linear distribution of massive young stars which cross nearly the full width of the arm. The measured 3-dimensional motion of G034.43+0.24 indicates a near-circular Galactic orbit. The water masers display arc-like distributions, possibly bow shocks, associated with winds from one or more massive young stars.

We have shown previously that the Murchison Widefield Array (MWA), can detect hundreds of Interplanetary Scintillation (IPS) sources simultaneously across a field of view $\sim30^\circ$ in extent. To test if we can use this capability to track heliospheric structures, we undertook a search of 88 hours of MWA IPS data, and identified an observation likely to have a significant Coronal Mass Ejection (CME) in the field of view. We demonstrate that in a single 5-minute MWA observation we are able to localise and image a CME $\sim$33 hours after launch at an elongation of $\sim37^\circ$ from the Sun. We use IPS observables to constrain the kinematics of the CME, and describe how MWA IPS observations can be used in the future to make a unique contribution to heliospheric modelling efforts.

Typical accretion disks around massive protostars are hot enough for water ice to sublimate. We here propose to utilize the massive protostellar disks for investigating the collisional evolution of silicate grains with no ice mantle, which is an essential process for the formation of rocky planetesimals in protoplanetary disks around lower-mass stars. We for the first time develop a model of massive protostellar disks that includes the coagulation, fragmentation, and radial drift of dust. We show that the maximum grain size in the disks is limited by collisional fragmentation rather than by radial drift. We derive analytic formulas that produce the radial distribution of the maximum grain size and dust surface density in the steady state. Applying the analytic formulas to the massive protostellar disk of GGD27-MM1, where the grain size is constrained from a millimeter polarimetric observation, we infer that the silicate grains in this disk fragment at collision velocities above ~ 10 m/s. The inferred fragmentation threshold velocity is lower than the maximum grain collision velocity in typical protoplanetary disks around low-mass stars, implying that coagulation alone may not lead to the formation of rocky planetesimals in those disks. With future measurements of grain sizes in massive protostellar disks, our model will provide more robust constraints on the sticking property of silicate grains.

Line profiles from spatially unresolved stellar observations consist of a superposition of local line profiles that result from observing the stellar atmosphere under specific viewing angles. Line profile variability caused by stellar magnetic activity or planetary transit selectively varies the weight and/or shape of profiles at individual surface positions. The effect is usually modeled with radiative transfer calculations because observations of spatially resolved stellar surfaces are not available. For the Sun, we recently obtained a broadband spectroscopic atlas of the solar center-to-limb variation (CLV). We use the atlas to study systematic differences between largely used radiative transfer calculations and solar observations. We concentrate on four strong lines useful for exoplanet transmission analysis, and we investigate the impact of CLV on transmission and Rossiter-McLaughlin (RM) curves. Solar models used to calculate synthetic spectra tend to underestimate line core depths but overestimate the effect of CLV. Our study shows that CLV can lead to significant systematic offsets in transmission curves and particularly in RM curves; transmission curves centered on individual lines are overestimated by up to a factor of two by the models, and simulations of RM curves yield amplitudes that are off by up to 5--10\,m\,s$^{-1}$ depending on the line. For the interpretation of transit observations, it is crucial for model spectra that accurately reproduce the solar CLV to become available which, for now, is the only calibration point available.

Our current capability of space weather prediction in the Earth's radiation belts is limited to only an hour in advance using the real-time solar wind monitoring at the Lagrangian L1 point. To mitigate the impacts of space weather on telecommunication satellites, several frameworks were proposed to advance the lead time of the prediction. We develop a prototype pipeline called "Helio1D" to forecast ambient solar wind conditions (speed, density, temperature, tangential magnetic field) at L1 with a lead time of 4 days. This pipeline predicts Corotating Interaction Regions (CIRs) and high-speed streams that can increase high-energy fluxes in the radiation belts. The Helio1D pipeline connects the Multi-VP model, which provides real-time solar wind emergence at 0.14 AU, and a 1D MHD model of solar wind propagation. We benchmark the Helio1D pipeline for solar wind speed against observations for the intervals in 2004 - 2013 and 2017 - 2018. We developed a framework based on the Fast Dynamic Time Warping technique that allows us to continuously compare time-series outputs containing CIRs to observations to measure the pipeline's performance. In particular, we use this framework to calibrate and improve the pipeline's performance for operational forecasting. To provide timing and magnitude uncertainties, we model several solar wind conditions in parallel, for a total of 21 profiles corresponding to the various virtual targets including the Earth. This pipeline can be used to feed real-time, daily solar wind forecasting that aims to predict the dynamics of the inner magnetosphere and the radiation belts.

Galactic plane radio surveys play a key role in improving our understanding of a wide range of astrophysical phenomena from neutron stars and Galactic magnetic fields to stellar formation and evolution. Performing such a survey on the latest interferometric telescopes produces large data rates necessitating a shift towards fully or quasi-real-time data analysis with data being stored for only the time required to process them. This has instilled a need to re-devise scientific strategies and methods for the effective management of telescope observing time. We describe here the setup for the 3000 hour Max-Planck-Institut f\"ur Radioastronomie (MPIfR) MeerKAT Galactic Plane survey (MMGPS). The survey is unique by operating in a commensal mode. Key science objectives of the survey including the discovery of new pulsars and transients as well as studies of Galactic magnetism, the interstellar medium and star formation rates. We explain the strategy coupled with the necessary hardware and software infrastructure needed for data reduction in the imaging, spectral and time domains. We have so far discovered 78 new pulsars including 17 binary pulsars. We have also developed an imaging pipeline sensitive to the order of a few tens of micro-Jansky and a spatial resolution of a few arcseconds. Further science operations are about to commence with an S-Band receiver system built in-house and operated in collaboration with the South African Radio Astronomy Observatory (SARAO). Spectral line commissioning observations at S-band already illustrate the spectroscopic capabilities of this instrument. These results have not only opened new avenues for Galactic science but also laid a strong foundation for surveys with future telescopes like the Square Kilometre Array (SKA).

The optical emission line spectra of X-ray binaries (XRBs) are thought to be produced in an irradiated atmosphere, possibly the base of a wind, located above the outer accretion disc. However, the physical nature of - and physical conditions in - the line-forming region remain poorly understood. Here, we test the idea that the optical spectrum is formed in the transition region between the cool, geometrically thin part of the disc near the mid-plane and a hot, vertically extended atmosphere or outflow produced by X-ray irradiation. We first present a VLT X-Shooter spectrum of XRB MAXI J1820+070 in the soft state associated with its 2018 outburst, which displays a rich set of double-peaked hydrogen and helium recombination lines. Aided by ancillary X-ray spectra and reddening estimates, we then model this spectrum with the Monte Carlo radiative transfer code Python, using a simple biconical disc wind model inspired by radiation-hydrodynamic simulations of irradiation-driven outflows from XRB discs. Such a model can qualitatively reproduce the observed features; nearly all of the optical emission arising from the transonic 'transition region' near the base of the wind. In this region, characteristic electron densities are on the order of 10$^{12-13}$ cm$^{-3}$, in line with the observed flat Balmer decrement (H$\alpha$/H$\beta \approx 1.3$). We conclude that strong irradiation can naturally give rise to both the optical line-forming layer in XRB discs and an overlying outflow/atmosphere that produces X-ray absorption lines.

We collected the archival data of blazar OJ~287 from heterogeneous very long baseline interferometry (VLBI) monitoring programs at 2.3 GHz, 8.6 GHz, 15 GHz and 43 GHz. The data reduction and observable extraction of those multi-band multi-epoch observations are batch-processed consistently with our automated pipeline. We present the multivariate correlation analysis on the observables at each band. We employ the cross-correlation function to search the correlations and the Monte Carlo (MC) technique to verify the certainty of correlations. Several correlations are found. The foremost findings are the correlations between the core flux density and the jet position angles on different scales, which validated the plausible predictions of the jet with precession characteristics. Meanwhile, there is a variation in the offset between the core EVPA and the inner-jet position angle over time at 15~GHz and 43~GHz.

We explore the possibilities of measuring the longitudinal profile of individual air showers beyond $X_{\rm max}$ when using very dense radio arrays such as SKA. The low-frequency part of the Square Kilometre Array, to be built in Australia, features an enormous antenna density of about $50,000$ antennas in the inner core region of radius 500 m, with a frequency band from 50 to 350 MHz. From CoREAS simulations, a SKA-Low antenna model plus noise contributions, and adapted LOFAR analysis scripts, we obtain a resolution in the shower maximum $X_{\rm max}$ and energy that is considerably better than at LOFAR. Already from this setup, we show that at least one additional parameter of the longitudinal profile can be measured. This would improve mass composition analysis by measuring an additional composition-dependent quantity. Moreover, it would offer an opportunity to discriminate between the different predictions of hadronic interaction models, hence contributing to hadronic physics at energy levels beyond man-made accelerators.

The ancient Large Magellanic Cloud (LMC) globular cluster NGC 2005 has recently been reported to have an ex-situ origin, thus, setting precedents that the LMC could have partially formed from smaller merged dwarf galaxies. We here provide additional arguments from which we conclude that is also fairly plausible an in-situ origin of NGC 2005, based on the abundance spread of a variety of chemical elements measured in dwarf galaxies, their minimum mass in order to form globular clusters, the globular cluster formation imprints kept in their kinematics, and the recent modeling showing that explosions of supernovae are responsible for the observed chemical abundance spread in dwarf galaxies. The present analysis points to the need for further development of numerical simulations and observational indices that can help us to differentiate between two mechanisms of galaxy formation for the LMC, namely, a primordial dwarf or an initial merging event of smaller dwarfs.

The radio quasar NVSS J141922$-$083830 (J1419$-$0838) was initially classified as an uncategorised blazar-type object, following its detection in the $\gamma$-ray band with the Fermi space telescope. Later, using multi-waveband observations and modeling, its was found to be a flat-spectrum radio quasar (FSRQ). However, its radio emission has never been discussed in depth in the literature. Here we present a detailed analysis on the radio properties of J1419$-$0838 using archival interferometric imaging data at pc and kpc scales. We conclude that the flux density variations, the flat radio spectrum, the compact nature of the quasar structure at all scales, and the relativistic Doppler enhancement of the radio emission all support the previous classification as an FSRQ. We also investigated the short- and long-term mid-infrared (MIR) light curve of the quasar based on observations by the Wide-field Infrared Survey Explorer satellite, and found that there is significant variability on time-scales of days as well as years. Comparison of the MIR light curve to the times of previously reported $\gamma$-ray and optical flares shows no clear correlation between the events at different wavebands.

Precise measurements of the stellar orbits around Sagittarius A* have established the existence of a supermassive black hole (SMBH) at the Galactic centre (GC). Due to the interplay between the SMBH and dark matter (DM), the DM density profile in the innermost region of the Galaxy, which is crucial for the DM indirect detection, is still an open question. Among the most popular models in the literature, the theoretical spike profile proposed by Gondolo and Silk (1999; GS hereafter) is well adopted. In this work, we investigate the DM spike profile using updated data from the Keck and VLT telescopes considering that the presence of such an extended mass component may affect the orbits of the S-stars in the Galactic center. We examine the radius and slope of the generalized NFW spike profile, analyze the Einasto spike, and discuss the influence of DM annihilation on the results. Our findings indicate that an initial slope of $\gamma \gtrsim 0.92$ for the generalized NFW spike profile is ruled out at a 95% confidence level. Additionally, the spike radius $R_{\rm sp}$ larger than 21.5 pc is rejected at 95% probability for the Einasto spike with $\alpha=0.17$, which also contradicts the GS spike model. The constraints with the VLT/GRAVITY upper limits are also projected. Although the GS NFW spike is well constrained by the Keck and VLT observation of S2, an NFW spike with a weak annihilation cusp may still be viable, as long as the DM annihilation cross section satisfies $\left< \sigma v \right> \gtrsim 7.7\times 10^{-27}~{\rm cm^3\,s^{-1}} (m_{\rm DM}/100~{\rm GeV})$ at 95% level.

Context. It is common for massive stars to engage in binary interaction. In close binaries, the components can enter a contact phase, where both stars overflow their respective Roche lobes simultaneously. While there exist observational constraints on the stellar properties of such systems, the most detailed stellar evolution models that feature a contact phase are not fully reconcilable with those measurements. Aims. We aim to consistently model contact phases of binary stars in a 1D stellar evolution code. To this end, we develop the methodology to account for energy transfer in the common contact layers. Methods. We implement an approximative model for energy transfer between the components of a contact binary based on the von Zeipel theorem in the stellar evolution code MESA. We compare structure and evolution models with and without this transfer and analyze the implications for the observable properties of the contact phase. Results. Implementing energy transfer helps eliminating baroclinicity in the common envelope between the components of a contact binary, which, if present, would drive strong thermal flows. We find that accounting for energy transfer in massive contact binaries significantly alters the mass ratio evolution and can extend the lifetime of an unequal mass ratio contact system.

The detection of exoplanets with the radial velocity method consists in detecting variations of the stellar velocity caused by an unseen sub-stellar companion. Instrumental errors, irregular time sampling, and different noise sources originating in the intrinsic variability of the star can hinder the interpretation of the data, and even lead to spurious detections. In recent times, work began to emerge in the field of extrasolar planets that use Machine Learning algorithms, some with results that exceed those obtained with the traditional techniques in the field. We seek to explore the scope of the neural networks in the radial velocity method, in particular for exoplanet detection in the presence of correlated noise of stellar origin. In this work, a neural network is proposed to replace the computation of the significance of the signal detected with the radial velocity method and to classify it as of planetary origin or not. The algorithm is trained using synthetic data of systems with and without planetary companions. We injected realistic correlated noise in the simulations, based on previous studies of the behaviour of stellar activity. The performance of the network is compared to the traditional method based on null hypothesis significance testing. The network achieves 28 % fewer false positives. The improvement is observed mainly in the detection of small-amplitude signals associated with low-mass planets. In addition, its execution time is five orders of magnitude faster than the traditional method. The superior performance exhibited by the algorithm has only been tested on simulated radial velocity data so far. Although in principle it should be straightforward to adapt it for use in real time series, its performance has to be tested thoroughly. Future work should permit evaluating its potential for adoption as a valuable tool for exoplanet detection.

We show how magnetic accretion of positronium (electron-positron) plasma by primordial black holes might significantly contribute to the mass of dark matter in the present Universe. Assuming that background gamma radiation is primordial black hole Hawking radiation rules out Bondi accretion, while magnetic accretion known from studies of active galactic nuclei could explain the abundance of dark matter. Various accretion scenarios are discussed.

The transformation and evolution of a galaxy is strongly influenced by interactions with its environment. Neutral hydrogen (\HI) is an excellent way to trace these interactions. Here, we present \HI\ observations of the spiral galaxy NGC~895, which was previously thought to be isolated. High-sensitivity \HI\ observations from the MeerKAT large survey project MIGHTEE reveal possible interaction features, such as extended spiral arms, and the two newly discovered \HI\ companions, that drive us to change the narrative that it is an isolated galaxy. We combine these observations with deep optical images from the Hyper Suprime Camera to show an absence of tidal debris between NGC 895 and its companions. We do find an excess of light in the outer parts of the companion galaxy MGTH$\_$J022138.1-052631 which could be an indication of external perturbation and thus possible sign of interactions. Our analysis shows that NGC~895 is an actively star-forming galaxy with a SFR of $\mathrm{1.75 \pm 0.09 [M_{\odot}/yr]}$, a value typical for high stellar mass galaxies on the star forming main sequence. It is reasonable to state that different mechanisms may have contributed to the observed features in NGC~895 and this emphasizes the need to revisit the target with more detailed observations. Our work shows the high potential and synergy of using state-of-the-art data in both \HI\ and optical to reveal a more complete picture of galaxy environments.

The cluster R136 in the giant star-forming region 30 Doradus in the Large Magellanic Cloud (LMC) offers a unique opportunity to resolve a stellar population in a starburst-like environment. We obtain the near-infrared to ultraviolet extinction towards 50 stars in the core of R136, employing the `extinction without standards' method. To assure good fits over the full wavelength range, we combine and modify existing extinction laws. We detect a strong spatial gradient in the extinction properties across the core of R136, coinciding with a gradient in density of cold gas that is part of a molecular cloud lying northeast of the cluster. In line with previous measurements of R136 and the 30 Doradus region, we obtain a high total-to-relative extinction ($R_V = 4.38 \pm 0.87$). However, the high values of $R_V$ are accompanied by relatively strong extinction in the ultraviolet, contrary to what is observed for Galactic sightlines. The relatively strong ultraviolet extinction suggests that the properties of the dust towards R136 differ from those in the Milky Way. For $R_{V} \sim 4.4$, about three times fewer ultraviolet photons can escape from the ambient dust environment relative to the canonical Galactic value of $R_{V} \sim 3.1$ at the same $A_{V}$. Therefore, if dust in the R136 star-bursting environment is characteristic for cosmologically distant star-bursting regions, the escape fraction of ultraviolet photons from such regions is overestimated by a factor of three relative to the standard Milky Way assumption for the total-to-selective extinction. Furthermore, a comparison with average curves tailored to other regions of the LMC shows that large differences in ultraviolet extinction exist within this galaxy. Further investigation is required in order to decipher whether or not there is a relation between $R_V$ and ultraviolet extinction in the LMC.

Ultra-hot Jupiters, with their high equilibrium temperatures and resolved spectral lines, have emerged as a perfect testbed for new analysis techniques in the study of exoplanet atmospheres. In particular, the resolved sodium doublet as a resonant line has proven a powerful indicator to probe the atmospheric structure over a wide pressure range. We explore an atmospheric origin of the observed blueshifted feature next to the sodium doublet of the ultra-hot Jupiter WASP-121~b, using a partial transit obtained with the 4-UT mode of ESPRESSO. We study its atmospheric dynamics visible across the terminator by splitting the data into mid-transit and egress. We determine that the blueshifted high-velocity absorption component is generated only during the egress part of the transit when a larger fraction of the day side of the planet is visible. For the egress data, MERC retrieves the blueshifted high-velocity absorption component as an equatorial day-to-night side wind across the evening limb, with no zonal winds visible on the morning terminator with weak evidence compared to a model with only vertical winds. For the mid-transit data, the observed line broadening is attributed to a vertical, radial wind. We attribute the equatorial day-to-night side wind over the evening terminator to a localised jet and restrain its existence between the substellar point and up to $10^\circ$ to the terminator in longitude, an opening angle of the jet of at most $60^\circ$ in latitude, and a lower boundary in altitude between [1.08, 1.15] $R_p$. Due to the partial nature of the transit, we cannot make any statements on whether the jet is truly super-rotational and one-sided or part of a symmetric day-to-night side atmospheric wind from the hotspot.

[Abridged] Aim: We aim to derive a new mass-loss rate prescription for RSGs that is not afflicted with some uncertainties inherent in preceding studies. Methods: We have observed CO rotational line emission towards a sample of RSGs in the open cluster RSGC1 that all are of similar initial mass. The ALMA CO(2-1) line detections allow to retrieve the gas mass-loss rates (Mdot_CO). In contrast to mass-loss rates derived from the analysis of dust spectral features (Mdot_SED), the data allow a direct determination of the wind velocity and no uncertain dust-to-gas correction factor is needed. Results: Five RSGs in RSGC1 have been detected in CO(2-1). The retrieved Mdot_CO values are systematically lower than Mdot_SED. Although only five RSGs in RSGC1 have been detected, the data allow to propose a new mass-loss rate relation for M-type red supergiants that is dependent on luminosity and initial mass. The new mass-loss rate relation is based on the new Mdot_CO values for the RSGs in RSGC1 and on prior Mdot_SED values for RSGs in 4 clusters, including RSGC1. The new Mdot-prescription yields a good prediction for the mass-loss rate of some well-known Galactic RSGs that are observed in multiple CO rotational lines, including alpha Ori, mu Cep and VX Sgr. However, there are indications that a stronger, potentially eruptive, mass-loss process - different than captured by our new mass-loss rate prescription - is occurring during some fraction of the RSG lifetime. Implementing a lower mass-loss rate in evolution codes for massive stars has important consequences for the nature of their end-state. A reduction of the RSG mass-loss rate implies that quiescent RSG mass loss is not enough to strip a single star's hydrogen-rich envelope. Upon core-collapse such single stars would explode as RSG.

A group of massive galaxies at redshifts of $z\geq 6.5$ have been recently detected by James Webb Space Telescope (JWST), which were unexpected to form at such early times within the standard Big Bang cosmology. In this work we propose that the formation of some $\sim 50~M_\odot$ primordial black holes (PBHs) formed in the early Universe via super-Eddington accretion within the dark matter halo can explain these observations. These PBHs may act as seeds for early galaxies formation with masses of $\sim 10^{9}-10^{10}~M_\odot$ at $z\sim 8$, hence accounting for the JWST observations. We use a hierarchical Bayesian inference framework to constrain the PBH mass distribution models, and find that the Lognormal model with the $M_{\rm c}\sim 35M_\odot$ is strongly preferred over other hypotheses. These rapidly growing BHs are expected to have strong radiation and may appear as the high-redshift compact objects, similar to those recently discovered by JWST.

The light curve of Gaia23bab (= SPICY 97589) shows two significant ($\Delta G>2$ mag) brightening events, one in 2017 and an ongoing event starting in 2022. The source's quiescent spectral energy distribution indicates an embedded ($A_V>5$ mag) pre-main-sequence star, with optical accretion emission and mid-infrared disk emission. This characterization is supported by the source's membership in an embedded cluster in the star-forming cloud DOBASHI 1604 at a distance of $900\pm45$~pc. Thus, the brightening events are probable accretion outbursts, likely of EX Lup-type.

Access to knowledge of the point spread function (PSF) of adaptive optics(AO)-assisted observations is still a major limitation when processing AO data. This limitation is particularly important when image analysis requires the use of deconvolution methods. As the PSF is a complex and time-varying function, reference PSFs acquired on calibration stars before or after the scientific observation can be too different from the actual PSF of the observation to be used for deconvolution, and lead to artefacts in the final image. We improved the existing PSF-estimation method based on the so-called marginal approach by enhancing the object prior in order to make it more robust and suitable for observations of resolved extended objects. Our process is based on a two-step blind deconvolution approach from the literature. The first step consists of PSF estimation from the science image. For this, we made use of an analytical PSF model, whose parameters are estimated based on a marginal algorithm. This PSF was then used for deconvolution. In this study, we first investigated the requirements in terms of PSF parameter knowledge to obtain an accurate and yet resilient deconvolution process using simulations. We show that current marginal algorithms do not provide the required level of accuracy, especially in the presence of small objects. Therefore, we modified the marginal algorithm by providing a new model for object description, leading to an improved estimation of the required PSF parameters. Our method fulfills the deconvolution requirement with realistic system configurations and different classes of Solar System objects in simulations. Finally, we validate our method by performing blind deconvolution with SPHERE/ZIMPOL observations of the Kleopatra asteroid.

The discovery of a broad, $\sim$1.5$^{\circ}$ long filamentary [OIII] 5007 emission $\sim$1.2$^{\circ}$ south-east of the M31 nucleus has recently been reported. More than 100 hours of exposures of a wide field (3.48$^{\circ} \times 2.32^{\circ}$) have allowed this pioneering detection based on 30 \AA\ narrow-band filters and several small refractors equipped with large cameras. We report a first velocity measurement in this extensive [OIII] emission line region. We used the low-resolution spectrograph MISTRAL (R $\sim$ 750), a facility of the Haute-Provence Observatory 193 cm telescope. The velocity measurement is based on the H$\alpha$, [NII], [SII] and [OIII] lines. The best solution to fit the spectrum indicates that the H$\alpha$ and [OIII] emissions are at the same heliocentric line-of-sight velocity of -96$\pm$4 km s$^{-1}$. This was measured within an area of $\sim$250 arcsec$^2$ selected on a bright knot along the long filament of $\sim$1.5$^{\circ}$, together with a [OIII]5007 surface brightness of 4.2$\pm$2.1 10$^{-17}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$. This agrees moderately well with the previous measurement. We also estimated the H$\alpha$/[NII] line ratio as $\sim$1.1. The radial velocities at which the H$\alpha$ and [OIII] lines were detected seem to show that these hydrogen and oxygen atoms belong to the same layer, but we cannot exclude that another weaker [OIII] line, belonging to another structure, that is, at another velocity, is below our detection threshold. Different scenarios have been considered to explain this filamentary structure...

The innermost ejecta of core-collapse supernovae are considered to be the sources of some iron-group and heavier nuclei. The ejecta are predominantly driven by neutrino heating, principally due to neutrino capture on free neutrons and protons. Such neutrino interaction plays a crucial role for setting neutron richness in the ejecta. Recent hydrodynamics work with sophisticated neutrino transport indicates that the ejecta are only mildly neutron rich or even proton rich. In such conditions a wide variety of trans-iron isotopes are synthesized, while the neutron richness is insufficient for the production of r-process nuclei. In this capter, basic concepts of nucleosynthesis in neutrino-heated ejecta and neutrino-driven winds of core-collapse supernovae are presented along with latest studies. Neutrino-heated ejecta indicate the early ejecta component within the first few seocnds in which anisotropic convective activities of material above the neutrinosphere become important for nucleosynthesis. Then, neutrino-driven winds follow, which are approximately isotropic outflows emerging from the surface of a proto-neutron star . According to such characteristics, studies of nucleosynthesis here are based on recent multi-dimentional hydrodynamics simulations and semi-alalytic wind solutions, respectively. These studies suggest that trans-iron species up to the atomic mass number of 90, as well as some rare isotopes such as 48Ca and 92Mo, are produced in the neutrino-heated ejecta. Neutrino-driven winds are unlikely sources of r-process nuclei, but rather promising sources of proton-rich isotopes up to the atomic number of 110.

The Euclid mission of the European Space Agency will perform a survey of weak lensing cosmic shear and galaxy clustering in order to constrain cosmological models and fundamental physics. We expand and adjust the mock Euclid likelihoods of the MontePython software in order to match the exact recipes used in previous Euclid Fisher matrix forecasts for several probes: weak lensing cosmic shear, photometric galaxy clustering, the cross-correlation between the latter observables, and spectroscopic galaxy clustering. We also establish which precision settings are required when running the Einstein-Boltzmann solvers CLASS and CAMB in the context of Euclid. For the minimal cosmological model, extended to include dynamical dark energy, we perform Fisher matrix forecasts based directly on a numerical evaluation of second derivatives of the likelihood with respect to model parameters. We compare our results with those of other forecasting methods and tools. We show that such MontePython forecasts agree very well with previous Fisher forecasts published by the Euclid Collaboration, and also, with new forecasts produced by the CosmicFish code, now interfaced directly with the two Einstein-Boltzmann solvers CAMB and CLASS. Moreover, to establish the validity of the Gaussian approximation, we show that the Fisher matrix marginal error contours coincide with the credible regions obtained when running Monte Carlo Markov Chains with MontePython while using the exact same mock likelihoods. The new Euclid forecast pipelines presented here are ready for use with additional cosmological parameters, in order to explore extended cosmological models.

At the highest stellar masses (log(\mstar) $\gtrsim$ 11.5 \msun), only a small fraction of galaxies are disk-like and actively star-forming objects. These so-called `super spirals' are ideal objects to better understand how galaxy evolution proceeds and to extend our knowledge about the relation between stars and gas to a higher stellar mass regime. We present new CO(1-0) data for a sample of 46 super spirals and for 18 slightly lower-mass (log(\mstar) $>$ 11.0 \msun ) galaxies with broad HI lines -- HI fast-rotators (HI-FRs). We analyze their molecular gas mass, derived from CO, in relation to their star formation rate (SFR) and stellar mass, and compare the results to values and scaling relations derived from lower-mass galaxies. We confirm that super spirals follow the same star-forming main sequence (SFMS) as lower-mass galaxies. We find that they possess abundant molecular gas, which lies above the extrapolation of the scaling relation with stellar mass derived from lower-mass galaxies, but within the relation between \mmol/\mstar and the distance to the SFMS. The molecular gas depletion time, \taudep = \mmol/SFR, is higher than for lower-mass galaxies on the SFMS (\taudep = 9.30 $\pm$ 0.03, compared to \taudep = 9.00 $\pm$ 0.02 for the comparison sample) and seems to continue an increasing trend with stellar mass. HI-FR galaxies have an atomic-to-molecular gas mass ratio that is in agreement with that of lower-mass galaxies, indicating that the conversion from the atomic to molecular gas proceeds in a similar way. We conclude that the availability of molecular gas is a crucial factor to enable star formation to continue and that, if gas is present, quenching is not a necessary destiny for high-mass galaxies. The difference in gas depletion time suggests that the properties of the molecular gas at high stellar masses are less favorable for star formation.

We aim to narrow down the origin of the non-detection of the metastable HeI triplet at about 10830 A obtained for the hot Jupiter WASP-80b. We measure the X-ray flux of WASP-80 from archival observations and use it as input to scaling relations accounting for the coronal [Fe/O] abundance ratio to infer the extreme-ultraviolet (EUV) flux in the 200-504 A range, which controls the formation of metastable HeI. We run three dimensional (magneto) hydrodynamic simulations of the expanding planetary upper atmosphere interacting with the stellar wind to study the impact on the HeI absorption of the stellar high-energy emission, the He/H abundance ratio, the stellar wind, and the possible presence of a planetary magnetic field up to 1 G. For a low stellar EUV emission, which is favoured by the measured logR'HK value, the HeI non-detection can be explained by a solar He/H abundance ratio in combination with a strong stellar wind, or by a sub-solar He/H abundance ratio, or by a combination of the two. For a high stellar EUV emission, the non-detection implies a sub-solar He/H abundance ratio. A planetary magnetic field is unlikely to be the cause of the non-detection. The low EUV stellar flux, driven by the low [Fe/O] coronal abundance, is the likely primary cause of the HeI non-detection. High-quality EUV spectra of nearby stars are urgently needed to improve the accuracy of high-energy emission estimates, which would then enable one to employ the observations to constrain the planetary He/H abundance ratio and the stellar wind strength. This would greatly enhance the information that can be extracted from HeI atmospheric characterisation observations.

We present an in-depth study of the late-time near-infrared plateau in Type Ia supernovae (SNe Ia), which occurs between 70-500 d. We double the existing sample of SNe Ia observed during the late-time near-infrared plateau with new observations taken with the Hubble Space Telescope, Gemini, New Technology Telescope, the 3.5m Calar Alto Telescope, and the Nordic Optical Telescope. Our sample consists of 24 nearby SNe Ia at redshift < 0.025. We are able to confirm that no plateau exists in the Ks band for most normal SNe Ia. SNe Ia with broader optical light curves at peak tend to have a higher average brightness on the plateau in J and H, most likely due to a shallower decline in the preceding 100 d. SNe Ia that are more luminous at peak also show a steeper decline during the plateau phase in H. We compare our data to state-of-the-art radiative transfer models of nebular SNe Ia in the near-infrared. We find good agreement with the sub-Mch model that has reduced non-thermal ionisation rates, but no physical justification for reducing these rates has yet been proposed. An analysis of the spectral evolution during the plateau demonstrates that the ratio of [Fe II] to [Fe III] contribution in a near-infrared filter determines the light curve evolution in said filter. We find that overluminous SNe decline slower during the plateau than expected from the trend seen for normal SNe Ia

We revisit a method to incorporate the Vainshtein screening mechanism in N-body simulations proposed by R. Scoccimarro in~\cite{Scoccimarro:2009eu}. We further extend this method to cover a subset of Horndeski theories that evade the bound on the speed of gravitational waves set by the binary neutron star merger GW170817. The procedure consists of the computation of an effective gravitational coupling that is time and scale dependent, $G_{\rm eff}\left(k,z\right)$, where the scale dependence will incorporate the screening of the fifth-force. This is a fast procedure that when contrasted to the alternative of solving the full equation of motion for the scalar field inside N-body codes, reduces considerably the computational time and complexity required to run simulations. To test the validity of this approach in the non-linear regime, we have implemented it in a COmoving Lagrangian Approximation (COLA) N-body code, and ran simulations for two gravity models that have full N-body simulation outputs available in the literature, nDGP and Cubic Galileon. We validate the combination of the COLA method with this implementation of the Vainshtein mechanism with full N-body simulations for predicting the boost function: the ratio between the modified gravity non-linear matter power spectrum and its General Relativity counterpart. This quantity is of great importance for building emulators in beyond-$\Lambda$CDM models, and we find that the method described in this work has an agreement of below $2\%$ for scales down to $k \approx 3h/$Mpc with respect to full N-body simulations.

Sub-GeV dark matter particles can annihilate or decay producing e^\pm pairs which upscatter the low-energy photon fields in the Galaxy and generate an X-ray emission (via the Inverse Compton effect). Using X-ray data from Xmm-Newton, Integral, NuStar and Suzaku, we derive new constraints on this class of dark matter (DM). In the annihilation case, our new bounds are the strongest available for DM masses above 180 MeV, reaching <sigma v> < 10^-28 cm^3/s for m_DM ~ 1 GeV. In the decay case, our bounds are the strongest to date essentially in the whole considered mass range, constraining tau > 10^28 s for m_DM ~ 1 GeV and improving by up to 3 orders of magnitude upon existing limits.

Particles with masses much larger than the inflationary Hubble scale, $H_I$, can be pair-produced non-adiabatically during inflation. Due to their large masses, the produced particles modify the curvature perturbation around their locations. These localized perturbations eventually give rise to localized signatures on the Cosmic Microwave Background (CMB), in particular, pairwise hotspots (PHS). In this work, we show that Convolutional Neural Networks (CNN) provide a powerful tool for identifying PHS on the CMB. While for a given hotspot profile a traditional Matched Filter Analysis is known to be optimal, a Neural Network learns to effectively detect the large variety of shapes that can arise in realistic models of particle production. Considering an idealized situation where the dominant background to the PHS signal comes from the standard CMB fluctuations, we show that a CNN can isolate the PHS with $\mathcal{O}(10)\%$ efficiency even if the hotspot temperature is $\mathcal{O}(10)$ times smaller than the average CMB fluctuations. Overall, the CNN search is sensitive to heavy particle masses $M_0/H_I=\mathcal{O}(200)$, and constitutes one of the unique probes of very high energy particle physics.

In the framework of a simple gravitational theory that contains a scalar field minimally coupled to gravity, we investigate the emergence of analytic black-hole solutions with non-trivial scalar hair of secondary type. Although it is possible for one to obtain asymptotically (A)dS solutions using our setup, in the context of the present work, we are solely interested in asymptotically flat solutions. At first, we study the properties of static and spherically symmetric black-hole solutions emanating from both regular and phantom scalar fields. We find that the regular-scalar-field-induced solutions are solutions describing ultra-compact black holes, while the phantom scalar fields generate ultra-sparse black-hole solutions. The latter are black holes that can be potentially of very low density since, contrary to ultra-compact ones, their horizon radius is always greater than the horizon radius of the corresponding Schwarzschild black hole of the same mass. Then, we generalize the above static solutions to slowly rotating ones and compute their angular velocities explicitly. Finally, the study of the axial perturbations of the derived solutions takes place, in which we show that there is always a region in the parameter space of the free parameters of our theory that allows the existence of both ultra-compact and ultra-sparse black holes.

ESA EUCLID mission will be launched in 2020 to understand the nature of the dark energy responsible of the accelerated expansion of the Universe and to map the geometry of the dark matter. The map will investigate the distanceredshift relationship and the evolution of cosmic structures thanks to two instruments: the NISP and the VIS. The NISP (Near Infrared Spectro-Photometer) is operating in the near-IR spectral range (0.9-2$\mu$m) with two observing modes: the photometric mode for the acquisition of images with broad band filters, and the spectroscopic mode for the acquisition of slitless dispersed images on the detectors. The spectroscopic mode uses four low resolution grisms to cover two spectral ranges: three ''red'' grisms for 1250-1850nm range, with three different orientations, and one ''blue'' grism for 920- 1300nm range. The NISP grisms are complex optical components combining four main optical functions: a grism function (dispersion without beam deviation of the first diffracted order) done by the grating on the prism hypotenuse, a spectral filter done by a multilayer filter deposited on the first face of the prism to select the spectral bandpass, a focus function done by the curved filter face of the prism (curvature radius of 10m) and a spectral wavefront correction done by the grating which grooves paths are nor parallel, neither straight. The development of these components have been started since 10 years at the Laboratoire d'Astrophysique de Marseille (LAM) and was linked to the project phases: prototypes have been developed to demonstrate the feasibility, then engineering and qualification models to validate the optical and mechanical performance of the component, finally the flight models have been manufactured and tested and will be installed on NISP instrument. In this paper, we present the optical performance of the four EUCLID NISP grisms flight models characterized at LAM: wavefront error, spectral transmission and grating groove profiles. The test devices and the methods developed for the characterization of these specific optical components are described. The analysis of the test results have shown that the grisms flight models for NISP are within specifications with an efficiency better than 70% on the spectral bandpass and a wavefront error on surfaces better than 30nm RMS. The components have withstood vibration qualification level up to 11.6g RMS in random test and vacuum cryogenics test down to 130K with measurement of optical quality in transmission. The EUCLID grisms flight models have been delivered to NISP project in November 2017 after the test campaign done at LAM that has demonstrated the compliance to the specifications.

We study dynamical dark energy models within Einstein's theory by means of matter perturbations and the growth index $\gamma$. Within four-dimensional General Relativity, we assume that dark energy does not cluster, and we adopt a linear ansatz for the growth index to investigate its impact on the deceleration parameter, $q$, and on the dark energy equation-of-state parameter, $w$. Following this approach, we identify a relationship between $q_0$ (today's value of $q$) and $\gamma$, which to the best of our knowledge is new. For $w(z)$, we find that in most of the cases considered it crosses the -1 line (quintom) ending at a present day value $w_0 > -1$. Furthermore, we show that an analytic expression for $w(z)$ may be obtained in the form of order (4,4) (or higher) Pad{\'e} parameterizations.

Mars, lacking an intrinsic dynamo, is an ideal laboratory to comparatively study induced magnetospheres, which can be found in other terrestrial bodies as well as comets. Additionally, Mars is of particular interest to further exploration due to its loss of habitability by atmospheric escape and possible future human exploration. In this context, we propose the Mars Magnetospheric Multipoint Measurement Mission (M$^5$), a multi-spacecraft mission to study the dynamics and energy transport of the Martian induced magnetosphere comprehensively. Particular focus is dedicated to the largely unexplored magnetotail region, where signatures of magnetic reconnection have been found. Furthermore, a reliable knowledge of the upstream solar wind conditions is needed to study the dynamics of the Martian magnetosphere, especially the different dayside boundary regions but also for energy transport phenomena like the current system and plasma waves. This will aid the study of atmospheric escape processes of planets with induced magnetospheres. In order to resolve the three-dimensional structures varying both in time and space, multi-point measurements are required. Thus, M$^5$ is a five spacecraft mission, with one solar wind monitor orbiting Mars in a circular orbit at 5 Martian radii, and four smaller spacecraft in a tetrahedral configuration orbiting Mars in an elliptical orbit, spanning the far magnetotail up to 6 Mars radii with a periapsis within the Martian magnetosphere of 1.8 Mars radii. We not only present a detailed assessment of the scientific need for such a mission but also show the resulting mission and spacecraft design taking into account all aspects of the mission requirements and constraints such as mass, power, and link budgets. This mission concept was developed during the Alpbach Summer School 2022.