Papers for Wednesday, Nov 29 2023

Neutron star emission originates typically from its magnetosphere due to radiating electrons. Trajectories of relativistic charged particles under uniform electromagnetic fields can be calculated analytically. However, under more complex and realistic fields, numerical solutions are required. Two common schemes are the Boris method, which solves the full equations of motion, and the guiding center approximation (GCA), which only evolves the orbital center. We compare both methods in a series of tests, discuss their characteristics and quantify their accuracy. We apply the methods to dipolar, quadrupolar, and quadrudipolar magnetic fields, as applicable for magnetospheres. It is essential to consider such realistic neutron star magnetic field geometries to model the emission from magnetars and pulsars. Our work can assist the Neutron Star Interior Composition ExploreR (NICER) to understand emission from the surface and magnetosphere of neutron stars and to study their composition.

Gravitational waves have been instrumental in providing deep insights into the nature of gravity. Next-generation detectors, such as the Einstein Telescope, are predicted to have a higher detection rate given the increased sensitivity and lower cut-off frequency. However, this increased sensitivity raises challenges concerning parameter estimation due to the foreseeable overlap of signals from multiple sources. Overlapping signals (OSs), if not properly identified, may introduce biases in estimating post-Newtonian (PN) coefficients in parameterized tests of general relativity (GR). We investigate how OSs affect $-1$PN to 2PN terms in parameterized GR tests, examining their potential to falsely suggest GR deviations. We estimate the prevalence of such misleading signals in next-generation detectors, and their collective influence on GR tests. We compare the effects of OSs on coefficients at different PN orders, concluding that overall the 1PN coefficient suffers the most. Our findings also reveal that while a non-negligible portion of OSs exhibit biases in PN coefficients that might individually prefer to conclude deviations from GR, collectively, the direction to deviate is random and a statistical combination will still be in favor of GR.

Monte-Carlo techniques are standard numerical tools for exploring non-Gaussian and multivariate likelihoods. Many variants of the original Metropolis-Hastings algorithm have been proposed to increase the sampling efficiency. Motivated by Ensemble Monte Carlo we allow the number of Markov chains to vary by exchanging particles with a reservoir, controlled by a parameter analogous to a chemical potential $\mu$, which effectively establishes a random process that samples microstates from a macrocanonical instead of a canonical ensemble. In this paper, we develop the theory of macrocanonical sampling for statistical inference on the basis of Bayesian macrocanonical partition functions, thereby bringing to light the relations between information-theoretical quantities and thermodynamic properties. Furthermore, we propose an algorithm for macrocanonical sampling, $\texttt{Avalanche Sampling}$, and apply it to various toy problems as well as the likelihood on the cosmological parameters $\Omega_m$ and $w$ on the basis of data from the supernova distance redshift relation.

We characterize the Magellanic Corona model of the formation of the Magellanic Stream, which we introduced in Lucchini et al. (2020, 2021). Using high-resolution hydrodynamic simulations, we constrain the properties of the primordial Magellanic Clouds, including the Magellanic Corona -- the gaseous halo around the Large Magellanic Cloud (LMC). With an LMC mass of $1.75\times10^{11}$ M$_\odot$, a Magellanic Corona of $>5\times10^9$ M$_\odot$ at $3\times10^5$ K, a total Small Magellanic Cloud mass $<10^{10}$ M$_\odot$, and a Milky Way corona of $2\times10^{10}$ M$_\odot$, we can reproduce the observed total mass of the neutral and ionized components of the Trailing Stream, ionization fractions along the Stream, morphology of the neutral gas, and on-sky extent of the ionized gas. The inclusion of advanced physical routines in the simulations allow the first direct comparison of a hydrodynamical model with UV absorption-line spectroscopic data. Our model reproduces O I, O VI, and C IV observations from HST/COS and FUSE. The stripped material is also nearby ($<50$ kpc from the Sun), as found in our prior models including a Magellanic Corona.

We extend our recent work on black hole spin in X-ray binary systems to include an analysis of 189 archival NuSTAR observations from 24 sources. Using self-consistent data reduction pipelines, spectral models, and statistical techniques, we report an unprecedented and uniform sample of 36 stellar-mass black hole spin measurements based on relativistic reflection. This treatment suggests that prior reports of low spins in a small number of sources were generally erroneous; our comprehensive treatment finds that those sources also tend to harbor black holes with high spin values. Overall, within $1\sigma$ uncertainty, $\sim86\%$ of the sample is consistent with $a \geq 0.95$, $\sim94\%$ of the sample is consistent with $a\geq 0.9$, and $100\%$ is consistent with $a\geq 0.7$ (the theoretical maximum for neutron stars; $a = cJ/GM^{2}$). We also find that the high-mass X-ray binaries (HMXBs; those with A, B, or O-type companions) are consistent with $a\geq 0.9$ within $1\sigma$ errors; this is in agreement with the low-mass X-ray binary population and may be especially important for comparisons to black holes discovered in gravitational wave events. In some cases, different spectra from the same source yield similar spin measurements but conflicting values for the inclination of the inner disk; we suggest that this is due to variable disk winds obscuring the blue wing of the relativistic Fe K emission line. We discuss the implications of our measurements, the unique view of systematic uncertainties enabled by our treatment, and future efforts to characterize black hole spins with new missions.

We test the possibility of using a convolutional neural network to infer the inclination angle of a black hole directly from the incomplete image of the black hole's shadow in the $uv$-plane. To this end, we develop a proof-of-concept network and use it to explicitly find how the error depends on the degree of coverage, type of input and coverage pattern. We arrive at a typical error of $10^\circ$ at a level of absolute coverage $1\%$ (for a pattern covering a central part of the $uv$-plane), $0.3\%$ (pattern covering the central part and the periphery, the $0.3\%$ referring to the central part only), and $14\%$ (uniform pattern). These numbers refer to a network that takes both amplitude and phase of the visibility function as inputs. We find that this type of network works best in terms of the error itself and its distribution for different angles. In addition, the same type of network demonstrates similarly good performance on highly blurred images mimicking sources nearing being unresolved. In terms of coverage, the magnitude of the error does not change much as one goes from the central pattern to the uniform one. We argue that this may be due to the presence of a typical scale which can be mostly learned by the network from the central part alone.

High resolution infrared data has revealed several young stars in close proximity to Sgr A*. These stars may encounter extremely high dark matter densities. We examine scenarios where dark matter scatters on stellar gas, accumulates in stellar cores, and then annihilates. We study the stars S2, S62, S4711 and S4714 and find three observable effects. First, dark matter interactions can inhibit in situ star-formation close to Sgr A*, favoring scenarios where these stars migrate into the Galactic Center. Second, dark matter interactions can delay main sequence evolution, making stars older than they appear. Third, very high dark matter densities can inject enough energy to disrupt main sequence stars, allowing S-star observations to constrain the dark matter density near Sgr A*.

``Quasi-periodic eruptions'' (QPE) are recurrent nuclear transients with periods of several hours to almost a day, which thus far have been detected exclusively in the X-ray band. We have shown that many of the key properties of QPE flares (period, luminosity, duration, emission temperature, alternating long-short recurrence time behavior, source rates) are naturally reproduced by a scenario involving twice-per-orbit collisions between a solar-type star on a mildly eccentric orbit, likely brought into the nucleus as an extreme mass-ratio inspiral (EMRI), and the gaseous accretion disk of a supermassive black hole (SMBH). The flare is generated by the hot shocked debris expanding outwards from either side of the disk midplane, akin to dual miniature supernovae. Here, we consider the conditions necessary for disk-star collisions to generate lower-temperature flares which peak in the ultraviolet (UV) instead of the X-ray band. We identify a region of parameter space at low SMBH mass $M_{\bullet} \sim 10^{5.5}M_{\odot}$ and QPE periods $P \gtrsim 10$ hr for which the predicted flares are sufficiently luminous $L_{\rm UV} \sim 10^{41}$ erg s$^{-1}$ to outshine the quiescent disk emission at these wavelengths. The prospects to discover such ``UV QPEs'' with future satellite missions such as ULTRASAT and UVEX depends on the prevalence of very low-mass SMBH and the occurrence rate of stellar EMRIs onto them. For gaseous disks produced by the tidal disruption of stars, we predict that X-ray QPEs will eventually shut off, only to later reappear as UV-QPEs as the accretion rate continues to drop.

The central star and its energetic radiation fields play a vital role in setting the vertical and radial chemical structure of planet-forming disks. We present observations that, for the first time, clearly reveal the UV-irradiated surface of a protoplanetary disk. Specifically, we spatially resolve the atomic-to-molecular (C I-to-CO) transition in the IM Lup disk with ALMA archival observations of [C I] $^3$P$_1$-$^3$P$_0$. We derive a C I emitting height of z/r $\gtrsim$ 0.5 with emission detected out to a radius of ${\approx}$600 au. Compared to other systems with C I heights inferred from unresolved observations or models, the C I layer in the IM Lup disk is at scale heights almost double that of other disks, confirming its highly flared nature. C I arises from a narrow, optically-thin layer that is substantially more elevated than that of $^{12}$CO (z/r $\approx$ 0.3-0.4), which allows us to directly constrain the physical gas conditions across the C I-to-CO transition zone. We also compute a radially-resolved C I column density profile and find a disk-averaged C I column density of 2$\times10^{16}$ cm$^{-2}$, which is ${\approx}$3-20$\times$ lower than that of other disks with spatially-resolved C I detections. We do not find evidence for vertical substructures or spatially-localized deviations in C I due, e.g., to either an embedded giant planet or a photoevaporative wind that have been proposed in the IM Lup disk, but emphasize that deeper observations are required for robust constraints.

Understanding the entire history of the ionization state of the intergalactic medium (IGM) is at the frontier of astrophysics and cosmology. A promising method to achieve this is by extracting the damping wing signal from the neutral IGM. As hundreds of redshift $z>6$ quasars are observed, we anticipate determining the detailed time evolution of the ionization fraction with unprecedented fidelity. However, traditional approaches to parameter inference are not sufficiently accurate. We assess the performance of a simulation-based inference (SBI) method to infer the neutral fraction of the universe from quasar spectra. The SBI method adeptly exploits the shape information of the damping wing, enabling precise estimations of the neutral fraction $\left<x_{\rm HI}\right>_{\rm v}$ and the wing position $w_p$. Importantly, the SBI framework successfully breaks the degeneracy between these two parameters, offering unbiased estimates of both. This makes the SBI superior to the traditional method using a pseudo-likelihood function. We anticipate that SBI will be essential to determine robustly the ionization history of the Universe through joint inference from the hundreds of high-$z$ spectra we will observe.

We derive galaxy colour selections from Euclid and ground-based photometry, aiming to accurately define background galaxy samples in cluster weak-lensing analyses. These selections are implemented in the Euclid data analysis pipelines for galaxy clusters. Given any set of photometric bands, we develop a method for the calibration of optimal galaxy colour selections that maximises the selection completeness, given a threshold on purity. Such colour selections are expressed as a function of the lens redshift. We calibrate galaxy selections using ground-based $griz$ and Euclid $Y_{\rm E}J_{\rm E}H_{\rm E}$ bands. Both selections produce a purity higher than 97%. The $griz$ selection completeness ranges from 30% to 84% in the lens redshift range $z_{\rm l}\in[0.2,0.8]$. With the full $grizY_{\rm E}J_{\rm E}H_{\rm E}$ selection, the completeness improves by up to $25$ percentage points, and the $z_{\rm l}$ range extends up to $z_{\rm l}=1.5$. The calibrated colour selections are stable to changes in the sample limiting magnitudes and redshift, and the selection based on $griz$ bands provides excellent results on real and simulated external data sets. Furthermore, the calibrated selections provide stable results using alternative photometric aperture definitions obtained from different ground-based telescopes. The $griz$ selection is also purer at high redshift and more complete at low redshift compared to colour selections found in the literature. We show that the calibrated colour selections provide robust results even when observations from a single band are missing from the ground-based data. Finally, we show that colour selections imply variations within the 1$\sigma$ uncertainty in the mean multiplicative shear bias, $m$, for stage III surveys. The first Euclid data releases will provide further insights into the impact of background selections on $m$.

We present a comprehensive spectral analysis of the ultraluminous X-ray source Holmberg II X-1 using broadband and high-resolution X-ray spectra taken with the XMM-Newton satellite over a period of 19 years benefiting from a recent campaign. We tested several models for the broadband spectra among which a double thermal component provided a reasonable description for the continuum between 0.3-10 keV and enabled us to constrain the properties of the accretion disc. The Luminosity-Temperature trends of the inner and outer disc components broadly agree with the expectations for a thin disc, although the exact values of the slopes are slightly sensitive to the adopted model. However, all tested models show L-T trends which deviate from a power law above a bolometric luminosity of about 5 $\times \ 10^{39} $erg/s, particularly for the hot thermal component associated to the inner accretion flow. Assuming that such deviations are due to the accretion rate exceeding its Eddington limit or, most likely, the super-critical rate, a compact object with a mass 16-36 Msun, i.e. a stellar-mass black hole, is inferred. The time-averaged (2021) high resolution spectra present narrow emission lines at 1 keV primarily from Ne IX-X and a very strong at 0.5 keV from N VII, which indicate Ne-N-rich gas with non-Solar abundances. This favours a nitrogen-rich donor star, such as a blue/red supergiant, which has escaped from its native stellar cluster characterised by a low-metallicity environment.

It is widely accepted that the width-luminosity relation used to standardize normal Type Ia supernovae (SNe Ia) breaks down in underluminous, 1991bg-like SNe Ia. This breakdown may be due to the choice of parameter used as a stand-in for the width of the SN Ia light curve. Using the colour stretch parameter $s_\mathrm{BV}$ instead of older parameters resolves this issue. Here, I assemble a sample of 13 nearby 1991bg-like SNe Ia from the literature, all of which have independent host-galaxy distance moduli and little to no reddening. I use Gaussian process regression to fit the light curves of these SNe in $U/u$, $g$, $B$, $V$, $R/r$, and $I/i$, and measure their peak absolute magnitudes. I find statistically significant ($>5\sigma$ confidence level) correlations between the peak absolute magnitudes of the underluminous SNe and their $s_\mathrm{BV}$ values in the range $0.2<s_\mathrm{BV}<0.6$. These correlations are broadly consistent with fits to $s_\mathrm{BV}<0.7$ SNe Ia with preliminary $B$- and $V$-band peak absolute magnitudes from the Carnegie Supernova Project and significantly inconsistent with similar fits to normal and transitional SNe Ia (with $0.7<s_\mathrm{BV}<1.1$). The underluminous width-luminosity relation shown here needs to be properly calibrated with a homogeneous sample of 1991bg-like SNe Ia, after which it could be used as a rung on a new cosmological distance ladder. With surface-brightness fluctuations (or another non-Cepheid method) used to calibrate distances to nearby 1991bg-like SNe, such a ladder could produce an independent measurement of the Hubble-Lema\^{i}tre Constant, $H_\mathrm{0}$.

We present new high-spectral resolution observations (R = $\lambda/\Delta\lambda$ = 7000) of the filamentary nebula surrounding NGC 1275, the central galaxy of the Perseus cluster. These observations have been obtained with SITELLE, an imaging Fourier transform spectrometer installed on the Canada-France-Hawai Telescope (CFHT) with a field of view of $11\text{ arcmin }\times 11 \text{ arcmin}$ encapsulating the entire filamentary structure of ionised gas despite its large size of $80 \text{ kpc}\times50 \text{ kpc}$. Here, we present renewed flux, velocity and velocity dispersion maps that show in great detail the kinematics of the optical nebula at \sii$\lambda6716$, \sii$\lambda6731$, \nii$\lambda6584$, H$\alpha$(6563\AA), and \nii$\lambda6548$. These maps reveal the existence of a bright flattened disk-shaped structure in the core extending to r $\sim 10$ kpc and dominated by a chaotic velocity field. This structure is located in the wake of X-ray cavities and characterised by a high mean velocity dispersion of $134$ km/s. The disk-shaped structure is surrounded by an extended array of filaments spread out to $r\sim 50$ kpc that are 10 times fainter in flux, remarkably quiescent and has a uniform mean velocity dispersion of $44$ km/s. This stability is puzzling given that the cluster core exhibits several energetic phenomena. Based on these results, we argue that there are two mechanisms to form multiphase gas in clusters of galaxies: a first triggered in the wake of X-ray cavities leading to more turbulent multiphase gas and a second, distinct mechanism, that is gentle and leads to large-scale multiphase gas spread throughout the core.

We have discovered a $\delta$ Scuti pulsator in a tight binary (P = 1.053 d) with nine pulsation modes whose frequencies are between 38 and 56 d$^{-1}$. Each of these modes exhibits amplitude modulations and $\pi$-rad phase shifts twice per orbital cycle. Five of these modes exhibit amplitude and phase shifts that are readily explained by dipole pulsations along an axis that is aligned with the binary's tidal axis. The novelty of the system lies in the remaining four pulsation modes, which we show are dipole pulsations along an axis that is perpendicular to both the tidal axis and the binary's orbital angular momentum axis. There are additionally two pulsation modes whose amplitudes and phases do not change significantly with orbital phase; they are explained as dipole modes along an axis aligned with the orbital/rotation axis. Hence, we propose that TIC 184743498 is a tri-axial pulsator, the first of its kind.

Context. Regularly-spaced, star-forming regions along the spiral arms of nearby galaxies provide insight into the early stages and initial conditions of star formation. The regular separation of these star-forming regions suggests spiral arm instability as their origin. Aims. We explore the effects of magnetic fields on the spiral arm instability. Methods. We use three-dimensional global magnetohydrodynamical simulations of isolated spiral galaxies, comparing three different initial plasma $\beta$ values (ratios of thermal to magnetic pressure) of $\beta=\infty$, $50$, and $10$. We perform Fourier analysis to calculate the separation of the over-dense regions formed from the spiral instability. We then compare the separations with observations. Results. We find that the spiral arms in the hydro case ($\beta = \infty$) are unstable, with the fragments initially connected by gas streams, reminiscent of Kelvin-Helmholtz instability. In the $\beta = 50$ case, the spiral arms fragment, but the fragments separate earlier and tend to be elongated in the direction perpendicular to the spiral arms. However, in the $\beta = 10$ run the arms are stabilised against fragmentation by magnetic pressure. The spiral arms in the unstable cases fragment into regularly-spaced, over-dense regions. We determine their separation to be $\sim 0.5$ kpc in the hydro and $\sim 0.65$ kpc in the $\beta = 50$ case, both in agreement with the observations of nearby galaxies. We find a smaller median characteristic wavelength of the over-densities to be $0.73^{+0.31}_{-0.36}$ kpc in the hydro case, compared to $0.98^{+0.49}_{-0.46}$ kpc in the $\beta = 50$ case. Moreover, we find a higher growth rate of the over-densities in the $\beta = 50$ run compared to the hydro run. We observe magnetic hills and valleys along the fragmented arms in the $\beta = 50$ run, which is characteristic of the Parker instability.

Future x-ray observatories will require imaging detectors with fast readout speeds that simultaneously achieve or exceed the other high performance parameters of x-ray charge-coupled devices (CCDs) used in many missions over the past three decades. Fast readout will reduce the impact of pile-up in missions with large collecting areas while also improving performance in other respects like timing resolution. Event-driven readout, in which only pixels with charge from x-ray events are read out, can be used to achieve these faster operating speeds. Speedster-EXD550 detectors are hybrid CMOS detectors (HCDs) capable of event-driven readout, developed by Teledyne Imaging Sensors and Penn State University. We present initial results from measurements of the first of these detectors, demonstrating their capabilities and performance in both full-frame and event-driven readout modes. These include dark current, read noise, gain variation, and energy resolution measurements from the first two engineering-grade devices.

State-of-the-art 19th century spectroscopy led to the discovery of quantum mechanics, and 20th century spectroscopy led to the confirmation of quantum electrodynamics. State-of-the-art 21st century astrophysical spectrographs, especially ANDES at ESO's ELT, have another opportunity to play a key role in the search for, and characterization of, the new physics which is known to be out there, waiting to be discovered. We rely on detailed simulations and forecast techniques to discuss four important examples of this point: big bang nucleosynthesis, the evolution of the cosmic microwave background temperature, tests of the universality of physical laws, and a real-time model-independent mapping of the expansion history of the universe (also known as the redshift drift). The last two are among the flagship science drivers for the ELT. We also highlight what is required for the ESO community to be able to play a meaningful role in 2030s fundamental cosmology and show that, even if ANDES only provides null results, such `minimum guaranteed science' will be in the form of constraints on key cosmological paradigms: these are independent from, and can be competitive with, those obtained from traditional cosmological probes.

Understanding the gas content in galaxies, its consumption and replenishment, remains pivotal in our comprehension of the evolution of the Universe. Numerous studies have addressed this, utilizing various observational tools and analytical methods. These include examining low-transition $^{12}$CO millimeter rotational lines and exploring the far-infrared and the (sub-)millimeter emission of galaxies. With the capabilities of present-day facilities, much of this research has been centered on relatively bright galaxies. We aim at exploring the gas reservoirs of a more general type of galaxy population at $1.0\leq z\leq 3.0$. We stack ALMA 1.1 mm data to measure the gas content of a mass-complete sample down to $\sim10^{8.6}$ M$_\odot$ at $z=1$ ($\sim10^{9.2}$ M$_\odot$ at $z=3$), extracted from the HST/CANDELS sample in GOODS-S. The sample is composed of 5,530 on average blue ($<b-i>\sim0.12$ mag, $<i-H>\sim0.81$ mag), star-forming main sequence objects ($\Delta$MS$\sim-0.03$). We report measurements at $10^{10-11}$ M$_\odot$ and upper limits for the gas fractions at $10^{8-10}$ M$_\odot$. At $10^{10-11}$ M$_\odot$, our f$_{\mathrm{gas}}$, ranging from 0.32 to 0.48, agree well with other studies based on mass-complete samples down to $10^{10}$ M$_\odot$, and are lower than expected according to other works more biased to individual detections. At $10^{9-10}$ M$_\odot$, we obtain 3$\sigma$ upper limits for f$_{\mathrm{gas}}$ ranging from 0.69 to 0.77. These upper limits are on the level of the extrapolations of scaling relations based on mass-complete samples down to $10^{10}$ M$_\odot$. As such, it suggests that the gas content of low-mass galaxies is at most what is extrapolated from literature scaling relations. The comparison of our results with previous works reflects how the inclusion of bluer, less obscured, and more MS-like objects progressively pushes the gas level to lower values.

The eFEDS is a wide $\approx$ 140 deg$^2$ field that has extensive multiwavelength coverage. To improve the utility of the existing data, we use CIGALE to fit source Spectral Energy Distributions (SEDs) from X-rays to far-infrared (FIR) mainly to derive stellar masses (M*) and star-formation rates (SFRs) for normal galaxies and X-ray Active Galactic Nuclei (AGNs). The catalog consists of 2,057,027 galaxies and 10,373 X-ray AGNs located in the $\approx$ 60 deg$^2$ GAMA09 sub-field. Comparing our M* with other available catalogs and our SFRs with FIR-derived SFRs, we demonstrate the general reliability of our SED-fitting measurements. Our catalog is publicly available at 10.5281/zenodo.10127224.

We present an algorithm that mitigates the effects of charge migration due to the "brighter-fatter effect'' (BFE) that occurs for highly illuminated stars in the Teledyne HAWAII-2RG detectors used in the NIRCam, NIRISS, and NIRSpec science instruments aboard the James Webb Space Telescope (JWST). The impact of this effect is most significant for photometry and spectrophotometry of bright stars in data for which the point spread function (PSF) is undersampled, which is the case for several observing modes of the NIRISS instrument. The main impact of the BFE to NIRISS data is incorrect count rate determinations for pixels in the central regions of PSFs of bright stars due to jump detections that are caused by charge migration from peak pixels to surrounding pixels. The effect is especially significant for bright compact sources in resampled, distortion-free images produced by the drizzle algorithm: quantitatively, apparent flux losses of $>$ 50% can occur in such images due to the BFE. We describe the algorithm of the "charge_migration'' mitigation step that has been implemented in version 10.0 of the operational JWST calibration pipeline as of Dec 5, 2023. We illustrate the impact of this step in terms of the resulting improvements of the precision of imaging photometry of point sources. The algorithm renders the effects of the BFE on photometry and surface brightness measurements to stay within 1%.

N-body systems characterized by inverse square attractive forces may display a self similar collapse known as the gravo-thermal catastrophe. In star clusters, collapse is halted by binary stars, and a large fraction of Milky Way clusters may have already reached this phase. It has been speculated -- with guidance from simulations -- that macroscopic variables such as central density and velocity dispersion are governed post-collapse by an effective, low-dimensional system of ODEs. It is still hard to distinguish chaotic, low dimensional motion, from high dimensional stochastic noise. Here we apply three machine learning tools to state-of-the-art dynamical simulations to constrain the post collapse dynamics: topological data analysis (TDA) on a lag embedding of the relevant time series, Sparse Identification of Nonlinear Dynamics (SINDY), and Tests of Accuracy with Random Points (TARP).

How astrophysical systems translate the kinetic energy of bulk motion into the acceleration of particles to very high energies is a pressing question. SS 433 is a microquasar that emits TeV gamma-rays indicating the presence of high-energy particles. A region of hard X-ray emission in the eastern lobe of SS 433 was recently identified as an acceleration site. We observed this region with the Imaging X-ray Polarimetry Explorer and measured a polarization degree in the range 38% to 77%. The high polarization degree indicates the magnetic field has a well ordered component if the X-rays are due to synchrotron emission. The polarization angle is in the range -12 to +10 degrees (east of north) which indicates that the magnetic field is parallel to the jet. Magnetic fields parallel to the bulk flow have also been found in supernova remnants and the jets of powerful radio galaxies. This may be caused by interaction of the flow with the ambient medium.

Low dust opacity spectral indices ($\beta < 1$) measured in the inner envelopes of class 0/I young stellar objects (age $\sim 10^{4-5}$ yr) have been interpreted as the presence of (sub-)millimetre dust grains in these environments. The density conditions and the lifetimes of collapsing envelopes have proven unfavorable for the growth of solids up to millimetre sizes. As an alternative, magneto-hydrodynamical simulations suggest that protostellar jets and outflows might lift grains from circumstellar discs and diffuse them in the envelope. We reframe available data for the CALYPSO sample of Class 0/I sources and show tentative evidence for an anti-correlation between the value of $\beta_{1-3mm}$ measured in the inner envelope and the mass loss rate of their jets and outflows, supporting a connection between the two. We discuss the implications that dust transport from the disc to the inner envelope might have for several aspects of planet formation. Finally, we urge for more accurate measurements of both correlated quantities and extension of this work to larger samples, necessary to further test the transport scenario.

The ArmazoNes high Dispersion Echelle Spectrograph (ANDES) is the optical and near-infrared high-resolution echelle spectrograph envisioned for the European Extremely Large Telescope (ELT). We present a selection of science cases, supported by new calculations and simulations, where ANDES could enable major advances in the fields of stars and stellar populations. We focus on three key areas, including the physics of stellar atmospheres, structure, and evolution; stars of the Milky Way, Local Group, and beyond; and the star-planet connection. The key features of ANDES are its wide wavelength coverage at high spectral resolution and its access to the large collecting area of the ELT. These features position ANDES to address the most compelling and potentially transformative science questions in stellar astrophysics of the decades ahead, including questions which cannot be anticipated today.

We studied the temporal evolution of quasi-biennial oscillations (QBOs) using acoustic mode oscillation frequencies from the Global Oscillation Network Group. The data used here span over more than 25 yr, covering solar cycles 23 and 24 and the ascending phase of cycle 25. The analysis reveals that the QBO-like signals are present in both the cycles, but with different periods. The dominant QBO period in cycle 23 is found to be about 2 yr while it is about 3 yr in cycle 24. Furthermore, the quasi-biennial oscillatory signals are present only during the ascending and high-activity phases of cycle 23 and quickly weaken around 2005 during the declining phase. In comparison, the QBO signals are present throughout the cycle 24, starting from 2009 to 2017. We also explored the depth dependence in QBO signals and obtained a close agreement at all depths, except in the near-surface shear layer. A detailed analysis of the near-surface shear layer suggests that the source region of QBOs is probably within a few thousand kilometers just below the surface.

Milky Way-type galaxies are surrounded by a warm-hot gaseous halo containing a considerable amount of baryons and metals. The kinematics and spatial distribution of highly-ionized ion species such as \ion{O}{6} can be significantly affected by supernova (SN) explosions and early (pre-SN) stellar feedback (e.g., stellar winds, radiation pressure). Here, we investigate effects of stellar feedback on \ion{O}{6} absorptions in Milky Way-like galaxies by analyzing the suites of high-resolution hydrodynamical simulations under the framework of {\it SMUGGLE}, a physically motivated subgrid interstellar medium and stellar feedback model for the moving-mesh code {\sc Arepo}. We find that the fiducial run with the full suite of stellar feedback and moderate star formation activities can reasonably reproduce Galactic \ion{O}{6} absorptions observed by space telescopes such as {\it FUSE}, including the scale height of low-velocity ($|v_{\rm LSR}|< 100\, \rm km~s^{-1}$) \ion{O}{6}, the column density $-$ line width relation for high-velocity ($100 \leq |v_{\rm LSR}|< 400\, \rm km~s^{-1}$) \ion{O}{6}, and the cumulative \ion{O}{6} column densities. In contrast, model variations with more intense star formation activities deviate from observations further. Additionally, we find that the run considering only SN feedback is in broad agreement with the observations, whereas in runs without SN feedback this agreement is absent, which indicates a dominant role of SN feedback in heating and accelerating interstellar \ion{O}{6}. This is consistent with the current picture that interstellar \ion{O}{6} is predominantly produced by collisional ionization where mechanical feedback can play a central role. In contrast, photoionization is negligible for \ion{O}{6} production due to the lack of high-energy ($\gtrsim114\ {\rm eV}$) photons required.

The linear polarization of the Cosmic Microwave Background (CMB) is highly sensitive to parity-violating physics at the surface of last scattering, which might cause mixing of E and B modes, an effect known as {\it cosmic birefringence}. This has until recently been problematic to detect due to its degeneracy with the instrument polarization miscalibration angle. Recently, a possible detection of a non-zero cosmic-birefringence angle was reported at $\beta={0.35^\circ}\pm 0.14^\circ$, where the miscalibration angle was simultaneously determined and subtracted from the analysis. Starting from this claim, we exploit a simple map of $\beta$ to the coupling constant of a parity-violating term in a generic effective-field theory for Lorentz and CPT violation. We show that the reported constraint on $\beta$ is consistent with current one-sided upper bounds from CMB studies of spacetime-symmetry breaking, and we discuss the implications and interpretation of this detection.

We find that combined Planck cosmic microwave background, baryon acoustic oscillations and supernovae data analyzed under $\Lambda$CDM are in 4.9$\sigma$ tension with eBOSS Ly$\alpha$ forest in inference of the linear matter power spectrum at wavenumber $\sim 1 h\,\mathrm{Mpc}^{-1}$ and redshift = 3. Model extensions can alleviate this tension: running in the tilt of the primordial power spectrum ($\alpha_\mathrm{s} \sim -0.01$); a fraction $\sim (1 - 5)\%$ of ultra-light axion dark matter (DM) with particle mass $\sim 10^{-25}$ eV or warm DM with mass $\sim 90$ eV. The new DESI survey, coupled with high-accuracy modeling, will help distinguish the source of this discrepancy.

The neighborhood of the Galactic black hole boasts a plethora of extended interstellar gas and dust features as well as populations of compact (unresolved, or marginally resolved) features such as the G objects. Most are well manifested in the infrared. To disentangle and characterize the infrared structure of extended features and identify compact sources, we used 3.8~$\mu$m (L' filter) data from the NIRC2 imager at the Keck Observatory and 8.6~$\mu$m (PAH1 filter) data from the VISIR imager at the Very Large Telescope (VLT) to produce the highest-resolution mid-IR color-temperature map of the inner half-parsec of the Galactic Center to date. From this map, we compile a catalog of features that stand out from their background. In particular, we identify 33 compact sources that stand out against the local background temperature, 11 of which are newly identified and are candidates for being members of the G objects population. Additionally, we resolve and newly characterize the morphology of several known extended features. These results prepare the way for ongoing and future JWST studies that have access to a greater range of mid-infrared wavelengths, and thus will allow for refined estimation of the trends of dust temperatures.

We present a new method for measuring the mean age of old/intermediate stellar populations in resolved, metal-poor ($\rm \langle[Fe/H]\rangle \lesssim -1.5$) galaxies using only the morphology of the horizontal branch (HB) and an estimate of the average metallicity. We calculate the ratio of blue-to-red HB stars and the mass-weighted mean ages of 27 M31 satellite galaxies that have star formation histories (SFHs) measured from Hubble Space Telescope-based color-magnitude diagrams (CMDs) that include the oldest Main Sequence Turn-off (MSTO) ages. We find a strong correlation between mean age, metallicity, and HB morphology, for stellar populations older than $\sim6$~Gyr. The correlation allows us to predict a galaxy's mean age from its HB morphology to a precision of $\lesssim 1$~Gyr. We validate our method by recovering the correct ages of Local Group galaxies that have robust MSTO-based ages and are not in our calibration sample. We also use our technique to measure the mean ages of isolated field galaxies KKR25 ($11.21^{+0.70}_{-0.65}$~Gyr) and VV124 ($11.03^{+0.73}_{-0.68}$~Gyr), which indicate that their main star formation episodes may have lasted several Gyr and support the picture that they achieved their early-type characteristics (e.g., low gas content, low star formation activity) in isolation and not through environment. Because the HB is $\sim80\times$ brighter than the oldest MSTO, our method can provide precise characteristic ages of predominantly old galaxies at distances $\sim 9$ times farther. We provide our calibrations in commonly used HST/ACS filters.

Optical bandpass-filter observations can be simply processed to determine similar horizontal ammonia distributions above the Jovian cloud tops as mid-infrared and microwave observations. Current understanding of this distribution and its relationship to aerosol opacity, cloud-top pressure, and circulation is provided by atmospheric retrieval models using observations from major ground-based facilities and spacecraft. These techniques recover high fidelity information on the ammonia distribution but are limited in spatial and temporal coverage. Part of this coverage gap - upper tropospheric abundance - can be bridged by using continuum-divided ammonia and methane absorption images as suggested by Combes and Encrenaz [1979]. In 2020-21, Jupiter was imaged in the 645 nm ammonia absorption band and adjacent continuum bands, demonstrating that the spatially-resolved optical depth in that band could be determined with a 0.28-m Schmidt-Cassegrain telescope (SCT). In 2022, a 620 nm filter was added to include methane absorption images in the same wavelength range. Methane abundance provides a constant reference against which to determine the ammonia abundance, specifically the column-averaged mole fraction above the clouds. VLT/MUSE results are compared to these SCT results and those from the TEXES mid-infrared spectrometer used on the IRTF and the Gemini telescopes. Meridional and longitudinal features are examined, including the Equatorial Zone (EZ) ammonia enhancement, the North Equatorial Belt (NEB) depletion, depletion above the Great Red Spot (GRS), and suggested enhancements over bright plumes in the northern EZ. This work demonstrates meaningful ammonia monitoring that can provide synoptic coverage and continuity between spacecraft or major ground-based facility campaigns.

We construct a new catalog of the blue horizontal-branch (BHB) stars from the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) DR5 dataset, which contains 5355+81 BHB stars at high Galactic latitude (($|Glat|>20^{\circ}$). We combine the spectral line indices with a set of Balmer line profile selection criteria to identify the BHB stars. During the selection process, we use the line index of \ion{Ca}{2}\,K to exclude the metal-rich A-type dwarfs. We obtain their atmospheric parameters by cross-matching our BHB stars with the catalog provided by \citet{Xiang2022}. The results show that our sample is consistent with the theoretical $T_{\rm eff}$-log\,$g$ evolutionary tracks of the BHB stars, indicating that our method is robust for identifying BHB stars from the LAMOST spectra. Their spatial distribution indicates that most of our BHB stars are located in the inner halo or the disk of the Milky Way. Combined with other BHB samples from the literature, the BHB stars can cover a large Galactic volume, which makes it a better probe for studying the kinematics, dynamics, and structural characteristics of the Milky Way.

The measurements of masses and luminosities of massive stars play an important role in understanding the formation and evolution of their host galaxies. In this work, we present the measurement of masses and luminosities of 2,946 OB-type stars, including 78 O-type stars and 2,868 B-type stars, based on their stellar parameters (effective temperature, surface gravity, and metallicity) and PARSEC isochrones model. Our results show that the median mass and luminosity of the 2,946 OB-type stars are 5.4 M$_{\odot}$ and log(L/L$_{\odot}$)=3.2 with the median relative error of 21.4$\%$ and 71.1$\%$, respectively. A good agreement between our results estimated by using our method and those derived by using the orbital motions of binary stars from the literature is found for some B-type stars. In addition, we also fit the mass-luminosity relation of B-type stars by using our derived mass and the luminosity from $Gaia$ DR3.

Nonlinear gravitational evolution induces strong nonlinearities in the observed cosmological density fields, leading to positive off-diagonal correlations in the power spectrum covariance. This has caused the information saturation in the power spectrum, e.g., the neutrino mass constraints from the nonlinear power spectra are lower than their linear counterparts by a factor of $\sim2$ at $z=0$. In this paper, we explore how nonlinear reconstruction methods improve the cosmological information from nonlinear cosmic fields. By applying nonlinear reconstruction to cold dark matter fields from the Quijote simulations, we find that nonlinear reconstruction can improve the constraints on cosmological parameters significantly, nearly reaching the linear theory limit. For neutrino mass, the result is only $12\%$ lower than the linear power spectrum, i.e., the theoretical best result. This makes nonlinear reconstruction an efficient and useful method to extract neutrino information from current and upcoming galaxy surveys.

Lateral motions of spicules serve as vital indicators of transverse waves in the solar atmosphere, and their study is crucial for understanding the wave heating process of the corona. Recent observations have focused on "high-frequency" transverse waves (periods < 100 s), which have the potential to transport sufficient energy for coronal heating. These high-frequency spicule oscillations are distinct from granular motions, which have much longer time scales of 5-10 min. Instead, it is proposed that they are generated through the mode conversion from high-frequency longitudinal waves that arise from a shock steepening process. Therefore, these oscillations may not solely be produced by the horizontal buffeting motions of granulation but also by the leakage of p-mode oscillations. To investigate the contribution of p-modes, our study employs a two-dimensional magneto-convection simulation spanning from the upper convection zone to the corona. During the course of the simulation, we introduce a p-mode-like driver at the bottom boundary. We reveal a notable increase in the mean velocity amplitude of the transverse oscillations in spicules, ranging from 10% to 30%, and attribute this to the energy transfer from longitudinal to transverse waves. This effect results in an enhancement of the estimated energy flux by 30-80%.

In recent years, a Gaussian Process Regression (GPR) based framework has been developed for foreground mitigation from data collected by the LOw-Frequency ARray (LOFAR), to measure the 21-cm signal power spectrum from the Epoch of Reionization (EoR) and Cosmic Dawn. However, it has been noted that through this method there can be a significant amount of signal loss if the EoR signal covariance is misestimated. To obtain better covariance models, we propose to use a kernel trained on the {\tt GRIZZLY} simulations using a Variational Auto-Encoder (VAE) based algorithm. In this work, we explore the abilities of this Machine Learning based kernel (VAE kernel) used with GPR, by testing it on mock signals from a variety of simulations, exploring noise levels corresponding to $\approx$10 nights ($\approx$141 hours) and $\approx$100 nights ($\approx$1410 hours) of observations with LOFAR. Our work suggests the possibility of successful extraction of the 21-cm signal within 2$\sigma$ uncertainty in most cases using the VAE kernel, with better recovery of both shape and power than with previously used covariance models. We also explore the role of the excess noise component identified in past applications of GPR and additionally analyse the possibility of redshift dependence on the performance of the VAE kernel. The latter allows us to prepare for future LOFAR observations at a range of redshifts, as well as compare with results from other telescopes.

The mechanism of pair creation in electrosphere of compact astrophysical objects such as quark stars or neutron stars is revisited, paying attention to evaporation of electrons and acceleration of electrons and positrons, previously not addressed in the literature. We perform a series of numerical simulations using the Vlasov-Maxwell equations. The rate of pair creation strongly depends on electric field strength in the electrosphere. Despite Pauli blocking is explicitly taken into account, we find no exponential suppression of the pair creation rate at low temperatures. The luminosity in pairs increases with temperature and it may reach up to $L_\pm\sim 10^{52}$ erg/s, much larger than previously assumed.

The EBLM project aims to characterise very low-mass stars that are companions to solar-type stars in eclipsing binaries. We describe the history and motivation for this project, the methodology we use to obtain precise mass, radius and effective temperature estimates for very low-mass M-dwarfs, and review results of the EBLM study and those from related projects. We show that radius inflation in fully-convective stars is a more subtle effect than was previously thought based on less precise measurements, i.e. the mass-radius-effective temperature relations we observe for fully-convective stars in single-line eclipsing binaries show reasonable agreement with theoretical models, particularly if we account for the M-dwarf metallicity, as inferred from the analysis of the primary star spectrum.

We studied the light curve of the star CD-36 3202, observed by TESS for the presence of stellar spots and to analyze the rotationally modulated flare. We mainly wanted to model the light curve of this flare and estimate its location regarding stellar spots. The flare lasted approximately 27$\,$h. Using our tool new \texttt{findinc\_mc} we managed to estimate the inclination angle of the star to $70^\circ\pm8^\circ$. With \texttt{BASSMAN} we modeled the light curve of the CD-36 3202 and we estimated that three spots are present on the surface of this star. The mean temperature of the spots was about $4000\pm 765\,$K, and the total spottedness was on average $11.61\%\pm0.13\,$\%. We created a new tool named \texttt{MFUEA} to model rotationally modulated flares. Using this software we estimated the latitude of the flare long-duration event equal to $69^{+2}_{-1}\,$deg in latitude. Our estimation of the flare's location was the first recreation of the exact position of a flare compared with the spots. The flare is placed 12$^\circ$ from the center of the coolest spot. This makes the flare related to the magnetic processes above the active region represented by the spot. Removing the effects of rotational modulation from the flare light curve allowed us to correct the estimation of bolometric energy released during the event from $(1.15\pm 0.35)\times 10^{35}\,$erg to $(3.99\pm 1.22)\times 10^{35}\,$erg.

Ongoing improvements of sub-mm- and mm-range interferometers and single-dish radiotelescopes are progressively allowing the detailed study of planetary nebulae (PNe) in molecular species other than 12CO and 13CO. We are implementing a new set of tables for extending the capabilities of the morpho-kinematical modelling tool SHAPE+shapemol, so radiative transfer in molecular species beyond 12CO and 13CO, namely C17O, C18O, HCN, HNC, CS, SiO, HCO+, and N2H+, are enabled under the Large Velocity Gradient approximation with the ease of use of SHAPE. We present preliminary results on the simultaneous analysis of a plethora of IRAM-30m and HERSCHEL/HIFI spectra, and NOEMA maps of different species in the pre-PN nebula M~1-92, which show interesting features such as a previously undetected pair of polar, turbulent, high-temperature blobs, or a 17O/18O isotopic ratio of 1.7, which indicates the AGB should have turned C-rich, as opposed to the apparent nature of its O-rich nebula.

Close to 100 per cent of massive stars are thought to be in binary systems. The multiplicity of massive stars seems to be intrinsically linked to their formation and evolution, and Massive Young Stellar Objects are key in observing this early stage of star formation. We have surveyed three samples totalling hundreds of MYSOs ($>8M_\odot$) across the Galaxy from the RMS catalogue, using UKIDSS and VVV point source data, and UKIRT $K-$band imaging to probe separations between 0.8-9 arcsec (approx 1000-100,000 au). We have used statistical methods to determine the binary statistics of the samples, and we find binary fractions of $64\pm 4$ per cent for the UKIDSS sample, $53\pm 4$ per cent for the VVV sample, and $49\pm 8$ per cent for the RMS imaging sample. Also we use the $J-$ and $K-$band magnitudes as a proxy for the companion mass, and a significant fraction of the detected systems have estimated mass ratios greater than 0.5, suggesting a deviation from the capture formation scenario which would be aligned with random IMF sampling. Finally, we find that YSOs located in the outer Galaxy have a higher binary fraction than those in the inner Galaxy. This is likely due to a lower stellar background density than observed towards the inner Galaxy, resulting in higher probabilities for visual binaries to be physical companions. It does indicate a binary fraction in the probed separation range of close to 100 per cent without the need to consider selection biases.

With JWST slated to gain high fidelity time dependent data on brown dwarf atmospheres, it is highly anticipated to do the same for directly imaged, sub-Jupiter exoplanets. With this new capability, the need for a full 3D understanding to explain spectral features and their time dependence is becoming an vital aspect for consideration. To examine the atmospheric properties of directly imaged sub-Jupiter exoplanets, we use the three dimensional Exo-FMS general circulation model (GCM) to simulate a metal enhanced generic young sub-Jupiter object. We couple Exo-FMS to a kinetic chemistry scheme, a tracer based cloud formation scheme and a spectral radiative-transfer model to take into account the chemical and cloud feedback on the atmospheric thermochemical and dynamical properties. Our results show a highly complex feedback between clouds and chemistry onto the 3D temperature structure of the atmosphere, bringing about latitudinal differences and inducing time-dependent stormy features at photospheric pressures. This suggests a strong connection and feedback between the spatial cloud coverage and chemical composition of the atmosphere, with the temperature changes and dynamical motions induced by cloud opacity feedback driving chemical species behaviour. In addition, we also produce latitude dependent and time dependent spectra of our model to investigate atmospheric variability and periodicity in commonly used photometric bands. Overall, our efforts put the included physics in 3D simulations of exoplanets on par with contemporary 1D radiative-convective equilibrium modelling.

We present a noteworthy transient event in the optical light curves of ztf18aanlzzf (SDSS J161259.83+421940.3), identified as a Narrow Line Seyfert 1 (NLS1) exhibiting merging patterns in the optical image. The 16-year long-term archived light curve revealed that this target stays in a steady state, while three flares occurred within the past 5 years with time separations ranging from 1 year to 3.5 years. The flare patterns of rapid brightening and slow decline following the peak, coupled with distinctive spectral features with strong He {\sc ii} and rare appearance of Bowen fluorescence line emissions, indicate at least two Tidal Eruption Event (TDE) flares in ztf18aanlzzf with a time separation of 3.5 years. We also apply TiDE light curve modeling and yield a Black Hole (BH) mass of $\sim 10^{6}\ M_{\odot}$, which is consistent with the BH mass measured from single-epoch spectra. Besides, the observed time lags $3.90_{-2.00}^{+2.06}$ days between the g and r bands strongly disagree with the prediction of the standard accretion disk model, highlighting the intricate interaction in the inner region related to the TDE. The reoccurrence gap of these TDEs, surpassing the previously reported repeated TDEs, can be attributed to binary star tidal disruption by a binary SMBH. Notably, the frequent TDE flares observed in this ULIRG-like target align with findings in a previous report for another ULIRG, suggesting a potentially elevated TDE rate in ULIRGs. Systematic variability studies of ULIRGs may help verify whether ULIRGs indeed have higher TDE rates.

The Red Rectangle is a nebula surrounding the post-AGB star HD 44179. It is the prototype of a particular class of nebulae associated with post-AGB binaries characterised by the presence of stable circumbinary disks in (quasi-)Keplerian rotation. Here we present the results of new high-resolution (20-50 mas) ALMA observations of continuum and line emissions at 0.9 mm. The continuum maps are analysed through a simple model of dust emission, which can reproduce the observational data. We find that most dust emission in the Red Rectangle is concentrated in the central regions of the rotating disk and that the settlement of dust grains onto the equatorial plane is very significant, particularly in comparison with the much larger scale height displayed by the gas distribution. The diameter of the dust-emitting region is about 250 au, with a total width of about 50 au. This region coincides with the warm PDR where certain molecules (like HCN), CI, and CII are presumably formed, as well as probably PAHs. From the spectral index, we confirm the presence in the disk of large grains, with a typical radius of about 0.150 mm, which supports the long-lived hypothesis for this structure. We also confirm the existence of a compact ionised wind at the centre of the nebula, probably emerging from the accretion disk around the companion, for which we derive an extent of about 10 au and a total flux of 8 mJy. We also briefly present the results on molecular lines of 12CO, 13CO, and other less abundant species.

High-resolution absorption spectroscopy toward bright background sources has had a paramount role in understanding early galaxy formation, the evolution of the intergalactic medium and the reionisation of the Universe. However, these studies are now approaching the boundaries of what can be achieved at ground-based 8-10m class telescopes. The identification of primeval systems at the highest redshifts, within the reionisation epoch and even into the dark ages, and of the products of the first generation of stars and the chemical enrichment of the early Universe, requires observing very faint targets with a signal-to-noise ratio high enough to detect very faint spectral signatures. In this paper, we describe the giant leap forward that will be enabled by ANDES, the high-resolution spectrograph for the ELT, in these key science fields, together with a brief, non-exhaustive overview of other extragalactic research topics that will be pursued by this instrument, and its synergistic use with other facilities that will become available in the early 2030s.

Identification and follow up observations of the host galaxies of fast radio bursts (FRBs) not only help us understand the environments in which the FRB progenitors reside, but also provide a unique way of probing the cosmological parameters using the dispersion measures of FRBs and distances to their origin. A fundamental requirement is an accurate distance measurement to the FRB host galaxy, but for some sources viewed through the Galactic plane, optical/NIR spectroscopic redshifts are extremely difficult to obtain due to dust extinction. Here we report the first radio-based spectroscopic redshift measurement for an FRB host galaxy, through detection of its neutral hydrogen (HI) 21-cm emission using MeerKAT observations. We obtain an HI-based redshift of z = 0.0357 for the host galaxy of FRB 20230718A, an apparently non-repeating FRB detected in the CRAFT survey and localized at a Galactic latitude of -0.367 deg. Our observations also reveal that the FRB host galaxy is interacting with a nearby companion, which is evident from the detection of an HI bridge connecting the two galaxies. This result demonstrates the value of HI to obtain redshifts of FRBs at low Galactic latitudes and redshifts. Such nearby FRBs whose dispersion measures are dominated by the Milky Way can be used to characterise these components and thus better calibrate the remaining cosmological contribution to dispersion for more distant FRBs that provide a strong lever arm to examine the Macquart relation between cosmological DM and redshift.

In the context of cosmological weak lensing studies, intrinsic alignments (IAs) are one the most complicated astrophysical systematic to model, given the poor understanding of the physical processes that cause them. A number of modelling frameworks for IAs have been proposed in the literature, both purely phenomenological or grounded on a perturbative treatment of symmetry-based arguments. However, the accuracy with which any of these approaches is able to describe the impact of IAs on cosmic shear data, particularly on the comparatively small scales ($k\simeq 1\,{\rm Mpc}^{-1}$) to which this observable is sensitive, is not clear. Here we quantify the level of disagreement between the true underlying intrinsic alignments and the theoretical model used to describe them that can be allowed in the context of cosmic shear analyses with future Stage-IV surveys. We consider various models describing this ``IA residual'', covering both physics-based approaches, as well as completely agnostic prescriptions. The same qualitative results are recovered in all cases explored: for a Stage-IV cosmic shear survey, a mis-modelling of the IA contribution at the $\sim10\%$ level produces shifts of $\lesssim0.5\sigma$ on the final cosmological parameter constraints. Current and future IA models should therefore aim to achieve this level of accuracy, a prospect that is not unfeasible for models with sufficient flexibility.

Context. Bright heads of penumbral filaments, penumbral grains (PGs), show apparent horizontal motions inwards, towards the umbra, or outwards, away from the umbra. Aims. We aim to prove statistically whether the direction of PGs' apparent motion is related to the inclination of the surrounding magnetic field. Methods. We use spectropolarimetric observations of five sunspot penumbrae to compare magnetic inclinations inside PGs with those in their surroundings. The data are taken by three observatories: Hinode satellite, Swedish Solar Telescope, and GREGOR solar telescope. The direction of PGs' motion is determined by feature tracking. The atmospheric conditions in PGs and their surroundings, including the magnetic field information, are retrieved by means of height-stratified spectropolarimetric inversions. Results. On a sample of 444 inward- and 269 outward-moving PGs we show that 43% of the inward-moving PGs have magnetic inclination larger by $8^\circ \pm 4^\circ$ than the inclination in their surroundings and 51% of the outward-moving PGs have the inclination smaller by $13^\circ \pm 7^\circ$ than the surrounding one. The opposite relation of inclinations is observed at only one-fifth of the inward- and outward-moving PGs. Conclusions. Rising hot plasma in PGs surrounded by a less inclined magnetic field may adapt its trajectory to be more vertical, causing an inward apparent motion of PGs. Oppositely, it may be dragged by a more horizontal surrounding magnetic field such that an outward apparent motion is observed.

In the conventional / most studied local distance ladder measurements, Type Ia supernovae (SNe Ia) are used in two of the three rungs. In the second rung, their luminosities are calibrated by standard candles like Cepheids or Tip of the Red Giant Branch (TRGB). In the third rung, the high luminosities and standardizability allow SNe to be used to calibrate the `Hubble' relation between distances and redshifts. Locally, the majority of distance ladder analyses find a high value of the Hubble Constant $H_0$ of $>70$ km/s/Mpc. Given the discrepancy with the inferred value using CMB observations, great scrutiny must be given to the role supernovae play in measuring $H_0$. Here, we review the main methodology, the many crosschecks for the supernova component of the distance ladder, and the various systematics studied. We review the important role supernovae play to explain the small disagreements seen from various local analyses. We also discuss analyses that employ an inverse distance ladder, which use similar sets of supernovae, but in the reverse direction, and yield a low value of $H_0$. We conclude given all available evidence, it is difficult to find a way that a systematic in supernovae measurements, or a non-$\Lambda$CDM component of the universe which could be measured with supernovae, can help explain the Hubble tension.

In this work we present a unified model for the cosmological dark sector. The theory is based on a simple minimally coupled scalar field whose action only contains a canonical kinetic term and is invariant under transverse diffeomorphisms (TDiff). The model has the same number of free parameters as $\Lambda$CDM. We confront the predictions of the model at the background level with data from Planck 2018 CMB distance priors, Pantheon+ and SH0ES SNIa distance moduli, BAO data points from 6dFGS, BOSS, eBOSS and DES and measurements of the Hubble parameter from cosmic chronometers. The model provides excellent results in the joint fit analysis, showing very strong evidence compared to $\Lambda$CDM in the deviance information criterion (DIC). We also show that the Hubble tension between Planck 2018 and SH0ES measurements can be alleviated in the unified TDiff model although further analysis is still needed.

We present an analysis of ALMA multi-band dust-continuum observations for 28 spectroscopically-confirmed bright Lyman-break galaxies at $5<z<8$. Our sample consists of 11 galaxies at $z\sim6$ newly observed in our ALMA program, which substantially increases the number of $5<z<8$ galaxies with both rest-frame 88 and 158 $\mu{\rm m}$ continuum observations, allowing us to simultaneously measure the IR luminosity and dust temperature for a statistical sample of $z\gtrsim5$ galaxies for the first time. We derive the relationship between the UV slope ($\beta_{\rm UV}$) and infrared excess (IRX) for the $z\sim6$ galaxies, and find a shallower IRX-$\beta_{\rm UV}$ relation compared to the previous results at $z\sim2$--4. Based on the IRX-$\beta_{\rm UV}$ relation consistent with our results and the $\beta_{\rm UV}$-$M_{\rm UV}$ relation including fainter galaxies in the literature, we find a limited contribution of the dust-obscured star formation to the total SFR density, $\sim30\%$ at $z\sim6$. Our measurements of the dust temperature at $z\sim6-7$, $T_{\rm dust}=40.9_{-9.1}^{+10.0}\,{\rm K}$ on average, supports a gentle increase of $T_{\rm dust}$ from $z=0$ to $z\sim6$--7. Using an analytic model with parameters consistent with recent {\it{JWST}} results, we discuss that the observed redshift evolution of the dust temperature can be reproduced by an $\sim0.6\,{\rm dex}$ increase in the gas depletion timescale and $\sim0.4\,{\rm dex}$ decrease of the metallicity. The variety of $T_{\rm dust}$ observed at high redshifts can also be naturally explained by scatters around the star-formation main sequence and average mass-metallicity relation, including an extremely high dust temperature of $T_{\rm dust}>80\,{\rm K}$ observed in a galaxy at $z=8.3$.

The current expansion rate of the Universe, the Hubble constant $H_0$, is an important cosmological quantity. However, two different ways to measure its value do not agree -- building a low-redshift distance ladder leads to a higher value of $H_0$ than inferring it from high-redshift observations in a $\Lambda$CDM cosmology. Most approaches to solve this tension either act at very low redshift by modifying the local distance ladder, or at high redshift by introducing new physics that changes the normalization of the inverse distance ladder. Here we discuss a way to address the Hubble tension at intermediate redshifts instead. By keeping the low- and high-redshift normalizations unchanged, we find a violation of the distance duality in the redshift range where luminosity and angular diameter distances overlap. We 'solve' this problem by introducing a redshift-dependent systematic effect that brings the luminosity distance into agreement with the angular diameter distance. The resulting expansion history is no longer compatible with $\Lambda$CDM, but this can be fixed with a dynamical dark energy component. In this way, we are able to solve the Hubble tension at intermediate redshifts.

Transmission spectroscopy is a proven technique to study a transiting exoplanet's atmosphere. However, stellar surface inhomogeneities, spots and faculae, alter the observed transmission spectra: the stellar contamination effect. The variable nature of the stellar activity also makes it difficult to stitch together multi-epoch observations and evaluate any potential variability in the exoplanet's atmosphere. This paper introduces SAGE, a tool to correct for the time-dependent impact of stellar activity on transmission spectra. It uses a pixelation approach to model the stellar surface with spots and faculae, while fully accounting for limb-darkening and rotational line-broadening. The current version is designed for low to medium-resolution spectra. We used SAGE to evaluate stellar contamination for F to M-type hosts, testing various spot sizes and locations, and quantify the impact of limb-darkening. We find that limb-darkening enhances the importance of the spot location on the stellar disk, with spots close to the disk center impacting the transmission spectra more strongly than spots near the limb. Moreover, due to the chromaticity of limb darkening, the shape of the contamination spectrum is also altered. Additionally, SAGE can be used to retrieve the properties and distribution of active regions on the stellar surface from photometric monitoring. We demonstrate this for WASP-69 using TESS data, finding that two spots at mid-latitudes and a combined coverage fraction of $\sim$1% are favoured. SAGE allows us to connect the photometric variability to the stellar contamination of transmission spectra, enhancing our ability to jointly interpret transmission spectra obtained at different epochs.

Using deep, narrowband imaging of the nearby spiral galaxy M101, we present stellar age information across the full extent of the disk of M101. Our narrowband filters measure age-sensitive absorption features such as the Balmer lines and the slope of the continuum between the Balmer break and 4000 \r{A} break. We interpret these features in the context of inside-out galaxy formation theories and dynamical models of spiral structure. We confirm the galaxy's radial age gradient, with the mean stellar age decreasing with radius. In the relatively undisturbed main disk, we find that stellar ages get progressively older with distance across a spiral arm, consistent with the large-scale shock scenario in a quasi-steady spiral wave pattern. Unexpectedly, we find the same pattern across spiral arms in the outer disk as well, beyond the corotation radius of the main spiral pattern. We suggest that M101 has a dynamic, or transient, spiral pattern with multiple pattern speeds joined together via mode coupling to form coherent spiral structure. This scenario connects together the radial age gradient inherent to inside-out galaxy formation with the across-arm age gradients predicted by dynamic spiral arm theories across the full radial extent of the galaxy.

The evolved stages of massive stars are poorly understood, but invaluable constraints can be derived from spatially resolved observations of nearby red supergiants, such as Betelgeuse. ALMA observations of Betelgeuse showing a dipolar velocity field have been interpreted as evidence for a rotation rate of $v\sin i \sim 5\, \mathrm{km\, s^{-1}}$. This is two orders of magnitude larger than predicted by single-star evolution, leading to the suggestion that Betelgeuse is a binary merger product. We propose instead that the velocity field could be due to large-scale convective motions. The resulting surface velocity maps can sometimes be mistaken for rotation, especially when the turbulent motions are only partially resolved, as is the case for the current ALMA beam. We support this claim with 3D CO5BOLD simulations of non-rotating red supergiants post-processed to predict synthetic ALMA images and SiO spectra to compare with observed radial velocity maps. Our simulations show a $\sim 50\%$ chance to be interpreted as evidence for a rotation rate as high as claimed for Betelgeuse. We conclude that we need at least another ALMA observation to firmly establish whether Betelgeuse is indeed rapidly rotating. Such observations would also provide insight into the role of angular momentum and binary interaction in the late evolutionary stages. The data will further probe the structure and complex physical processes in the atmospheres of red supergiants, which are immediate progenitors of supernovae and are believed to be essential in the formation of gravitational wave sources.

We combine photometric data from GALEX GR6+7 AIS and Gaia EDR3 with stellar parameters from the SAGA and PASTEL catalogs to construct high-quality training samples for dwarfs ($\rm 0.4< BP-RP<1.6$) and giants ($\rm 0.6< BP-RP <1.6$). We apply careful reddening corrections using empirical temperature- and extinction-dependent extinction coefficients. Using the two samples, we establish a relationship between stellar loci (NUV$-$BP vs. BP$-$RP colors), metallicity, and $\rm M_G$. For a given BP$-$RP color, a 1 dex change in [Fe/H] corresponds to an approximately 1 magnitude change in NUV$-$BP color for solar-type stars. These relationships are employed to estimate metallicities based on NUV$-$BP, BP$-$RP, and $\rm M_G$. Thanks to the strong metallicity dependence in the GALEX NUV-band, our models enable a typical photometric-metallicity precision of approximately $\sigma_{\rm [Fe/H]}$ = 0.11 dex for dwarfs and $\sigma_{\rm [Fe/H]}$ = 0.17 dex for giants, with an effective metallicity range extending down to [Fe/H] $= -3.0$ for dwarfs and [Fe/H] $= -4.0$ for giants. We also find that the NUV-band based photometric-metallicity estimate is not as strongly affected by carbon enhancement as previous photometric techniques. With the Gaia and GALEX data, we have estimated metallicities for about 5 million stars across almost the entire sky, including approximately 4.5 million dwarfs and 0.5 million giants. This work demonstrates the potential of the NUV-band for estimating photometric metallicities, and sets the groundwork for utilizing the NUV data from space telescopes such as the upcoming Chinese Space Station Telescope.

The spectroscopic survey of the China Space Station Telescope (CSST) is expected to obtain a huge number of slitless spectra, including more than one hundred million galaxy spectra and millions of active galactic nuclei (AGN) spectra. By making use of these spectra, we can measure the Baryon Acoustic Oscillation (BAO) signals over large redshift ranges with excellent precisions. In this work, we predict the CSST measurements of the post-reconstruction galaxy power spectra at 0<z<1.2 and pre-reconstruction AGN power spectra at 0<z<4, and derive the BAO signals at different redshift bins by constraining the BAO scaling parameters using the Markov Chain Monte Carlo method. Our result shows that the CSST spectroscopic survey can provide accurate BAO measurements with precisions higher than 1% and 3% for the galaxy and AGN surveys, respectively. By comparing with current measurements in the same range at low redshifts, this can improve the precisions by a factor of $2\sim3$, and similar precisions can be obtained in the pessimistic case. We also investigate the constraints on the cosmological parameters using the measured BAO data by the CSST, and obtain stringent constraint results for the energy density of dark matter, Hubble constant, and equation of state of dark energy.

The origin and maintenance of spiral structure in galaxies, the correlation between different types of spiral structure and several proposed mechanisms for their generation, and the evolution of spiral arms of galaxies with time are questions that are still controversial. In this note we study the spiral structure in a sample of distant galaxies in order to infer the evolution of spiral arm characteristics with redshift. We considered a sample of 171 face-on spiral galaxies in the Hubble Space Telescope COSMOS (The Cosmic Evolution Survey) field. The galaxies are distributed up to $z \approx 1$ with a mean value of 0.44. For all galaxies, we determined the pitch angles of the spiral arms and analysed their dependence on redshift; a total of 359 arms were measured. Analyses of our measurements combined with the literature data suggest a possible evolution of the pitch angles of spiral galaxies: by the modern epoch the spiral pattern, on average, becomes more tightly wound. This may be a consequence of the general evolution of the structure of galaxies as galaxies become more massive over time and their bulges grow. In addition, the distribution of the cotangent of pitch angle of galaxies indicates the possibility that the dominant mechanism of spiral pattern generation changes over time.

The ability to robustly determine galaxy properties such as masses, ages and star-formation rates is critically limited by the ability to accurately measure dust attenuation. Dust reddening is often characterized by comparing observations to models of either nebular recombination-lines or the ultra violet (UV) continuum. Here, we use a new technique to measure dust reddening by exploiting the HeII $\lambda$1640 to $\lambda$4686 emission lines originating from the stellar winds of Wolf-Rayet stars. The intrinsic line ratio is determined by atomic physics, enabling an estimate of the stellar reddening similar to how the Balmer lines probe reddening of gas emission. The HeII line ratio is measured from UV and optical spectroscopy using the Space Telescope Imaging Spectrograph (STIS) on board the Hubble Space Telescope (HST) for eight nearby galaxies hosting young massive star clusters. We compare our results to dust reddening values estimated from UV spectral slopes and from Balmer line ratios and find tentative evidence for systematic differences. The reddening derived from the He II lines tends to be higher, whereas that from the UV continuum tends to be lower. A larger sample size is needed to confirm this trend. If confirmed, this may indicate an age sequence probing different stages of dust clearing. Broad HeII lines have also been detected in galaxies more distant than our sample, providing the opportunity to estimate the dust reddening of the youngest stellar populations out to distances of $\sim$100 Mpc.

Asteroseismic modelling of isolated star presents significant challenges due to the difficulty in accurately determining stellar parameters, particularly the stellar age. These challenges can be overcomed by observing stars in open clusters, whose coeval members share an initial chemical composition. The light curves by TESS allow us to investigate and analyse stellar variations in clusters with an unprecedented level. We aim to detect gravity-mode oscillations in the early-type main-sequence members of the young open cluster NGC 2516. We selected the 301 member stars as our sample and analysed the TESS FFI light curves. We also collected high-resolution spectra using the FEROS for the g-mode pulsators. By fitting the theoretical isochrones to the colour-magnitude diagram (CMD) of a cluster, we determined an age of 102 $\pm$ 15 Myr and inferred the extinction at 550 nm ($A_0$) is 0.53 $\pm$ 0.04 mag. We identified 147 stars with surface brightness modulations, 24 with g-mode pulsations ($\gamma$ Doradus or Slowly Pulsating B stars), and 35 with p-mode pulsations ($\delta$ Sct stars). When sorted by colour index, the amplitude spectra of the $\delta$ Sct stars show a distinct ordering and reveal a discernible frequency-temperature relationship. The near-core rotation rates, measured from period spacing patterns in two SPB and nine $\gamma$ Dor stars, reach up to 3/d . This is at the high end of the values found from Kepler data of field stars of similar variability type. The $\gamma$ Dor stars have internal rotation rates as high as 50% of their critical value, whereas the SPB stars exhibit rotation rates close to their critical rate. We did not find long-term brightness and colour variations in the mid-infrared, which suggests that there are no disk or shell formation events in our sample. We also discussed the results of our spectroscopic observations for the g-mode pulsators.

Super-metal-rich (SMR) stars, currently in the solar neighbourhood, are expected to originate only in the inner Galaxy and have definitely migrated. We aim at studying a large sample of SMR stars to provide constraints on the epoch of the bar formation and its impact on the MW disc stellar populations. We investigate a sample of 169,701 MSTO and SGB stars with 6D phase space information and high-quality stellar parameters coming from the hybrid-CNN analysis of the Gaia-DR3 RVS stars. We compute distances and ages using the StarHorse code with a mean precision of 1% and 11%, respectively. From these, 11,848 stars have metallicity ([Fe/H]) above 0.15 dex. We report a metallicity dependence of spatial distribution of stellar orbits shown by the bimodal distribution in the guiding radius at 6.9 and 7.9 kpc, first appearing at [Fe/H]~0.1 dex, becoming very pronounced at larger [Fe/H]. In addition, we've observed a trend where the most metal-rich stars, with [Fe/H]~0.4 dex, are predominantly old (9-12 Gyrs) but show a gradual decline in [Fe/H] with age, reaching around 0.25 dex at about 4 Gyrs ago, followed by a sharp drop around 3 Gyrs ago. Furthermore, our full dataset reveals a clear peak in the age-metallicity relationship during the same period, indicating a SF burst around 3-4 Gyrs ago with slightly sub-solar [Fe/H] and enhanced [$\alpha$/Fe]. We show the SMR stars are good tracers of the bar activity. We interpret the steep decrease in number of SMR stars at around 3 Gyr as the end of the bar formation epoch. In this scenario, the peak of bar activity also coincides with a peak in the SF activity in the disc. Although the SF burst around 3 Gyr ago has been reported previously, its origin was unclear. Here, we suggest the SF burst to have been triggered by the high bar activity, 3-4 Gyr ago. According to these results and interpretation, the MW bar could be young.

The JWST Disk Infrared Spectral Chemistry Survey (JDISCS) aims to understand the evolution of the chemistry of inner protoplanetary disks using the Mid-InfraRed Instrument (MIRI) on the James Webb Space Telescope (JWST). With a growing sample of >30 disks, the survey implements a custom method to calibrate the MIRI Medium Resolution Spectrometer (MRS) to contrasts of better than 1:300 across its 4.9-28 micron spectral range. This is achieved using observations of Themis-family asteroids as precise empirical reference sources. High spectral contrast enables precise retrievals of physical parameters, searches for rare molecular species and isotopologues, and constraints on the inventories of carbon- and nitrogen-bearing species. JDISCS also offers significant improvements to the MRS wavelength and resolving power calibration. We describe the JDISCS calibrated data and demonstrate its quality using observations of the disk around the solar-mass young star FZ Tau. The FZ Tau MIRI spectrum is dominated by strong emission from warm water vapor. We show that the water and CO line emission originates from the disk surface and traces a range of gas temperatures of ~500-1500 K. We retrieve parameters for the observed CO and H2O lines, and show that they are consistent with a radial distribution represented by two temperature components. A high water abundance of n(H2O)~10^-4 fills the disk surface at least out to the 350 K isotherm at 1.5 au. We search the FZ Tau environs for extended emission detecting a large (radius of ~300 au) ring of emission from H2 gas surrounding FZ Tau, and discuss its origin.

It is believed that the first star-forming galaxies are the main drivers of cosmic reionization. It is usually assumed that there is a one-to-one relationship between the star formation rate (SFR) inside a galaxy and the host halo mass in semi-analytical/numerical modeling of large-scale ionization maps during the epoch of reionization. However, more accurate simulations and observations suggest that the SFR and ionizing luminosity in galaxies may vary considerably even if the host halo mass is the same. This astrophysical scatter can introduce an additional non-Gaussianity in the cosmological [H I]$_{\text{21cm}}$ signal, which is already expected to have an inherent time-evolving non-Gaussianity due to reionization. Here, we have studied the impact of the scatter on the [H I]$_{\text{21cm}}$ bispectrum using semi-numerical simulations. The scatter primarily affects small ionized regions, whereas the large ionized bubbles remain largely unaffected. Although, the fractional change in the [H I]$_{\text{21cm}}$ bispectra due to the scatter is found to be more than a factor of $10$ at large scales ($\lesssim 1\, {\rm Mpc}^{-1}$), it is found to be statistically not significant. However, we have found the impact due to the scatter to be significant ($|\langle \Delta B \rangle/B_{\text{no-scatter}}| \sim 1$) at small scales ($k\sim2.55\, {\rm Mpc}^{-1}$) and at $\overline{x}_{\rm HI}\sim 0.8$. We have also found that 1000 hours of SKA1-Low like observation is unlikely to detect the signatures of the scatter at these small scales.

Type II Supernova 1987A (SN 1987A), observed in 1987, released an energy of \mbox{$Q \approx 3 \times 10^{53}$ erg}}. This huge energy is essentially the magnitude of gravitational potential or self-gravitational energy (PE) of a new born cold neutron star having a gravitational compactness or redshift $z_b \approx 0.15$. One may wonder what could be the upper limit on the amount of energy that might be released with the formation of a cold Ultra Compact Object (UCO) with an arbitrary high $z_b$. Accordingly, here, for the first time, we obtain an analytical expression for the PE of a homogeneous general relativistic UCO assuming it to be cold and static. It is found that the PE of a homogeneous UCO of mass $M$ may exceed Mc$^2$ and be as large as 1.34 Mc$^2$. This result, though surprising, follows from an \textit{exact and correct} analytical calculation based on the standard General Theory of Relativity (GTR). Further, UCOs supported by tangential stresses may be inhomogeneous and much more massive than neutron stars with PE $\sim$ 2.1 Mc$^2$ Thus, in principle, formation of an UCO of a few solar masses ($M_\odot$) might release an energy $Q\sim10^{55}$ erg.

This study presents a recently developed two-way coupled magnetosphere-ionosphere model named MagPIE that enables the investigation of the impact of flux transfer events (FTEs) on the ionosphere. Our findings highlight the prominent role of cusp-FTE reconnection in influencing the ionosphere. The typical morphology of an FTE signal, represented by field-aligned currents (FACs) on the ionosphere, is shown to exhibit a distinct pattern characterized by an I-shaped patch surrounded by a U-shaped patch. Furthermore, we demonstrate that the effects of FACs resulting from FTEs may extend well into the region of closed field lines on the ionosphere. These FACs are seen to exhibit a remarkable resemblance to discrete dayside auroral arcs, providing further evidence that FTEs can be considered as a probable cause of such phenomena. Additionally, FTEs generate vortex-like patterns of ionospheric flow, which can manifest as either twin vortices or a combination of multiple vortices, depending on the characteristics of the FACs producing them. Furthermore, we present compelling evidence of morphological similarity between the simulated ionospheric signatures obtained from the MagPIE model and an observation made by the SWARM satellites. The agreement between our model and observational data further strengthens the credibility of our model and opens up new avenues to theoretically explore the complex ionospheric effects caused by FTEs.

We consider first order cosmological phase transitions (PT) happening at late times, below Standard Model (SM) temperatures $T_{\rm PT} \lesssim$ GeV. The inherently stochastic nature of bubble nucleation and the finite number of bubbles associated with a late-time PT lead to superhorizon fluctuations in the PT completion time. We compute how such fluctuations eventually source curvature fluctuations with universal properties, independent of the microphysics of the PT dynamics. Using Cosmic Microwave Background (CMB) and Large Scale Structure (LSS) measurements, we constrain the energy released in a dark-sector PT. For 0.1 eV $\lesssim T_{\rm PT} \lesssim$ keV this constraint is stronger than both the current bound from additional neutrino species $\Delta N_{\rm eff}$, and in some cases, even CMB-S4 projections. Future measurements of CMB spectral distortions and pulsar timing arrays will also provide competitive sensitivity for keV $\lesssim T_{\rm PT} \lesssim$ GeV.

We determine the solar neutrino fluxes from the global analysis of the most up-to-date terrestrial and solar neutrino data including the final results of the three phases of Borexino. The analysis are performed in the framework of three-neutrino mixing with and without accounting for the solar luminosity constraint. We discuss the independence of the results on the input from the Gallium experiments. The determined fluxes are then compared with the predictions provided by the latest Standard Solar Models. We quantify the dependence of the model comparison with the assumptions about the normalization of the solar neutrino fluxes produced in the CNO-cycle as well as on the particular set of fluxes employed for the model testing.

We consider an inflationary kinetic function with an integrable pole that is traversed during inflation. This scenario leads to enhanced spectra of primordial scalar inhomogeneities with detectable signals: formation of primordial black holes (that could explain Dark Matter) and scalar-induced gravitational waves (that could reproduce the recent Pulsar Timing Array observation, or predict signals in future detectors such as LISA or ET). Spectral signatures depend on whether the inflaton mass dimension at the pole is above or below 2. Values mildly below 2 allow a big power spectrum enhancement with a mild tuning. Finally, we discuss the possibility that a kinetic pole can arise as anomalous dimension of the inflaton due to quantum effects of Planckian particles that become light at some specific inflaton field value.

Axions can be stimulated to decay to photons by ambient photons of the right frequency or by photons from decay of neighboring axions. If the axion density is high enough the photon intensity can be amplified, which is a type of lasing or an axion maser. Here we review the astrophysical situations where axion lasing can appear and possibly be detected.

The abnormally fast orbital decay observed in the black hole (BH) Low-Mass X-ray binaries (BH-LMXB) A0620-00 and XTE J1118+480 can be explained by the dynamical friction between Dark Matter (DM) and the companion star orbiting around the low-mass BH (of a few $M_\odot$) of the system. In this case the value of the index $\gamma_{\rm sp}$ of the DM spike surrounding the BH can be pinned down with an accuracy of a few percent, way better than that for much bigger systems such as the super massive BHs (SMBHs) in the Galactic Center or in M87. We have used the data from XTE J1118+480 to put bounds on the WIMP annihilation cross section times velocity $\langle \sigma v\rangle$, assuming that DM annihilation is driven by the $b\bar{b}$ channel and that it proceeds in $s$-wave. The bounds are driven by the radio synchrotron signal produced by $e^\pm$ final states propagating in the magnetic field near the BH. For DM masses $m_\chi$ up to the TeV scale XTE J1118+480 allows to constrain $\langle \sigma v\rangle$ well below $\langle\sigma v\rangle_{\rm thermal}$, corresponding to the observed DM relic density in the Universe for a thermal WIMP. On the other hand, for $m_\chi \gtrsim$ 15 GeV the bounds from the SMBHs in the GC or in M87 do not reach $\langle\sigma v\rangle_{\rm thermal}$ when the very large uncertainties on the corresponding spike indices are taken into account, in spite of potentially producing much larger DM densities compared to XTE J1118+480. Our bounds for XTE J1118+480 have a mild sensitivity on spatial diffusion, but diffusion enhances the sensitivity of the results upon the intensity of the magnetic field. Taken at face value the bound from XTE J1118+480 on $\langle \sigma v\rangle$ is the most constraining compared to all others for $m_\chi\lesssim$ 1 TeV, unless the intensity of the magnetic field is smaller than its equipartition estimation.

Gravitational Wave (GW) observations from Neutron Stars (NS) in a binary system provide an excellent scenario to constrain the nuclear parameters. The investigation of Pratten et al. (2022) has shown that the ignorance of f-mode dynamical tidal correction in the GW waveform model of the binary neutron star (BNS) system can lead to substantial bias in the measurement of NS properties and NS equations of state (EOS). In this work, we investigate the bias in the nuclear parameters resulting from the ignorance of dynamical tidal correction. In addition, this work demonstrates the sensitivity of the nuclear parameters and the estimated constraints on them from future GW observations. We infer the nuclear parameters from GW observations by describing the NS matter within the relativistic mean field model. For a population of GW events, we notice that the ignorance of dynamical tide predicts a lower median for nucleon effective mass ($m^*$) by $\sim6\%$ compared to the scenario when dynamical tidal correction is considered. Whereas at a 90\% credible interval(CI), $m^*$ gets constrained up to $\sim 5\%$ and $\sim 3\%$ in A+ (the LIGO-Virgo detectors with a sensitivity of 5th observing run) and Cosmic Explorer (CE) respectively. We also discuss the resulting constraints on all other nuclear parameters, including compressibility, symmetry energy, and slope of symmetry energy, considering an ensemble of GW events. We do not notice any significant impact in analyzing nuclear parameters other than $m^*$ due to the ignorance of f-mode dynamical tides.

We compute the stochastic gravitational wave (GW) background generated by black hole-black hole (BH-BH) hyperbolic encounters with eccentricities close to one and compare them with the respective sensitivity curves of planned GW detectors. We use the Keplerian potential to model the orbits of the encounters and the quadrupole formula to compute the emitted GWs. We take into account hyperbolic encounters that take place in clusters up to redshift $5$, with BH's masses spanning from $5 M_{\odot}$ to $55 M_{\odot}$. We assume the clusters to be virialized and we study several cluster models with different mass and virial velocity, and finally, provide an accumulative result to display the background in an average sense. Using the maxima and the minima of our accumulative result for each frequency, we provide analytical expressions for both optimistic and pessimistic scenarios. Our results imply that the background from these encounters are likely to be detected by the third generation detectors Cosmic explorer, and Einstein telescope, while the tail section at lower frequencies intersects with DECIGO -- making it a possible source for both ground and space based future GW detectors.

We propose a novel scenario to explain the matter-antimatter asymmetry by twofold leptogenesis, wherein heavy Majorana neutrinos exhibit temperature-dependent masses and engage in $CP$-violating decays. This scenario envisages two distinct phases of leptogenesis: one occurring above the electroweak scale and the other below it. The sphaleron process converts the first lepton asymmetry to baryon asymmetry, but not the second one due to its decoupling. This mechanism potentially explains the significant discrepancy between baryon and lepton asymmetries, as suggested by recent observations of Helium-4. Furthermore, our model implies that the present masses of Majorana neutrinos are lighter than the electroweak scale, offering a tangible avenue for experimental verification in various terrestrial settings.

Tracking quintessence in a spatially flat and isotropic space-time, with a minimally coupled canonical scalar field and an asymptotically inverse power-law potential $V(\varphi)\propto\varphi^{-p}$, $p>0$, as $\varphi\rightarrow0$, is investigated by introducing a new three-dimensional \emph{regular} dynamical system. This enables a rigorous explanation of the tracking feature: 1) the dynamical system has a tracker fixed point $\mathrm{T}$ with a two-dimensional stable manifold that pushes an open set of nearby solutions toward a single tracker solution originating from $\mathrm{T}$; 2) all solutions, including the tracker solution and the solutions that track/shadow it, end at a common future attractor fixed point that depends on the potential. Thus, the open set of solutions that shadow the tracker solution share its properties during the tracking quintessence epoch. We also discuss similarities and differences of underlying mechanisms for tracking, thawing and scaling freezing quintessence, and, moreover, we illustrate with state space pictures that all of these types of quintessence exist simultaneously for certain potentials.

We propose a new method to determine the physical parameters of Kerr black holes (namely, specific angular momentum $a$, inclination angle $i$, and distance $D$ from an observer) only from the shadow's information such as the size and shape. Key points in our method are (i) to treat the distance as a free parameter, (ii) to expand the shadow's outline as a Fourier series, and (iii) to construct principal components from the Fourier coefficients. These points enable us to obtain a one-to-one mapping between three principal components, being observables characterizing the size and deviation of shadow's shape from a circular disk, and the values of three parameters ($D/M, a/M, i$), where $M$ is the mass of black hole. Our method is applicable to various type of black holes and even to ones with accretion disks.

We introduce STRAUSS (Sonification Tools and Resources for Analysis Using Sound Synthesis) a modular, self-contained and flexible Python sonification package, operating in a free and open source (FOSS) capacity. STRAUSS is intended to be a flexible tool suitable for both scientific data exploration and analysis as well as for producing sonifications that are suitable for public outreach and artistic contexts. We explain the motivations behind STRAUSS, and how these lead to our design choices. We also describe the basic code structure and concepts. We then present output sonification examples, specifically: (1) multiple representations of univariate data (i.e., single data series) for data exploration; (2) how multi-variate data can be mapped onto sound to help interpret how those data variables are related and; (3) a full spatial audio example for immersive Virtual Reality. We summarise, alluding to some of the future functionality as STRAUSS development accelerates.

We present a detailed investigation of the properties of the galactic rotation curves in the Weyl geometric gravity model, in which the gravitational action is constructed from the square of the Weyl curvature scalar, and of the strength of the Weyl vector. The theory admits a scalar-vector tensor representation, obtained by introducing an auxiliary scalar field. By assuming that the Weyl vector has only a radial component, an exact solution of the field equations can be obtained, which depends on three integration constants, and, as compared to the Schwarzschild solution, contains two new terms, linear and quadratic in the radial coordinate. In the framework of this solution we obtain the exact general relativistic expression of the tangential velocity of the massive test particles moving in stable circular orbits in the galactic halo. We test the theoretical predictions of the model by using 175 galaxies from the Spitzer Photometry \& Accurate Rotation Curves (SPARC) database. We fit the theoretical predictions of the rotation curves in conformal gravity with the SPARC data by using the Multi Start and Global Search methods. In the total expression of the tangential velocity we also include the effects of the baryonic matter, and the mass to luminosity ratio. Our results indicate that the simple solution of the Weyl geometric gravity can successfully account for the large variety of the rotation curves of the SPARC sample, and provide a satisfactory description of the particle dynamics in the galactic halos, without the need of introducing the elusive dark matter particle.

A superradiant cloud of ultralight bosons near a rotating black hole provides a smoking gun for particle physics in the infrared. However, tidal perturbations from a nearby binary companion can destabilise the boson cloud and even terminate superradiance. In this work, we consider the backreaction of superradiance termination to the dynamics of general binary orbits parametrised by their semi-latus rectum, eccentricity and inclination angle. Our analysis focuses on Extreme Mass Ratio Inspiral (EMRI) systems and employs the period-average approximation to derive evolution equations of these binary parameters in the Newtonian limit. We find that the binary evolution history can be significantly modulated by the backreaction towards large circular equatorial orbits with reduced termination rate. This process can generically happen even away from the resonance bands. Our work therefore serves as a first step towards probing ultralight bosons through the statistics of EMRI binary parameters in the future.