Papers for Monday, Nov 20 2023

Self-supervised learning has been known for learning good representations from data without the need for annotated labels. We explore the simple siamese (SimSiam) architecture for representation learning on strong gravitational lens images. Commonly used image augmentations tend to change lens properties; for example, zoom-in would affect the Einstein radius. To create image pairs representing the same underlying lens model, we introduce a lens augmentation method to preserve lens properties by fixing the lens model while varying the source galaxies. Our research demonstrates this lens augmentation works well with SimSiam for learning the lens image representation without labels, so we name it LenSiam. We also show that a pre-trained LenSiam model can benefit downstream tasks. We open-source our code and datasets at https://github.com/kuanweih/LenSiam .

We present MEGARA (Multi-Espectr\'ografo en GTC de Alta Resoluci\'on para Astronom\'ia) Integral Field Unit (IFU) observations of 5 local type-2 quasars (QSO2s, z $\sim 0.1$) from the Quasar Feedback (QSOFEED) sample. These active galactic nuclei (AGN) have bolometric luminosities of 10$^{45.5-46}$ erg/s and stellar masses of $\sim$10$^{11}$ M$_{\odot}$. We explore the kinematics of the ionized gas through the [O~III]$\lambda$5007 $\r{A}$ emission line. The nuclear spectra of the 5 QSO2s, extracted in a circular aperture of $\sim$ 1.2" ($\sim$ 2.2 kpc) in diameter, show signatures of high velocity winds in the form of broad (full width at half maximum; 1300$\leq$FWHM$\leq$2240 km/s and blueshifted components. We find that 4 out of the 5 QSO2s present outflows that we can resolve with our seeing-limited data, and they have radii ranging from 3.1 to 12.6 kpc. In the case of the two QSO2s with extended radio emission, we find that it is well-aligned with the outflows, suggesting that low-power jets might be compressing and accelerating the ionized gas in these radio-quiet QSO2s. In the four QSO2s with spatially resolved outflows, we measure ionized mass outflow rates of 3.3-6.5 Msun/yr when we use [S~II]-based densities, and of 0.7-1.6 Msun/yr when trans-auroral line-based densities are considered instead. We compare them with the corresponding molecular mass outflow rates (8 - 16 Msun/yr), derived from CO(2-1) ALMA observations at 0.2" resolution. Both phases show lower outflow mass rates than those expected from observational scaling relations where uniform assumptions on the outflow properties were adopted. This might be indicating that the AGN luminosity is not the only driver of massive outflows and/or that these relations need to be re-scaled using accurate outflow properties. We do not find a significant impact of the outflows on the global star formation rates.

The orbit distribution of young stars in the Galactic disk is highly structured, from well-defined clusters to streams of stars that may be widely dispersed across the sky, but are compact in orbital action-angle space. The age distribution of such groups can constrain the timescales over which co-natal groups of stars disperse into the `field'. Gaia data have proven powerful to identify such groups in action-angle space, but the resulting member samples are often too small and have too narrow a CMD coverage to allow robust age determinations. Here, we develop and illustrate a new approach that can estimate robust stellar population ages for such groups of stars. This first entails projecting the predetermined action-angle distribution into the 5D space of positions, parallaxes and proper motions, where much larger samples of likely members can be identified over a much wider range of the CMD. It then entails isochrone fitting that accounts for a) widely varying distances and reddenings; b) outliers and binaries; c) sparsely populated main sequence turn-offs, by incorporating the age information of the low-mass main sequence; and d) the possible presence of an intrinsic age spread in the stellar population. When we apply this approach to 92 nearby stellar groups identified in 6D orbit space, we find that they are predominately young ($\lesssim 1$ Gyr), mono-age populations. Many groups are established (known) localized clusters with possible tidal tails, others tend to be widely dispersed and manifestly unbound. This new age-dating tool offers a stringent approach to understanding on which orbits stars form in the solar neighborhood and how quickly they disperse into the field.

The recently-discovered stellar system Ursa Major III/UNIONS 1 (UMa3/U1) is the faintest known Milky Way satellite to date. With a stellar mass of $16^{+6}_{-5}\,\rm M_\odot$ and a half-light radius of $3\pm1$ pc, it is either the darkest galaxy ever discovered or the faintest self-gravitating star cluster known to orbit the Galaxy. Its line-of-sight velocity dispersion suggests the presence of dark matter, although current measurements are inconclusive because of the unknown contribution to the dispersion of potential binary stars. We use N-body simulations to show that, if self-gravitating, the system could not survive in the Milky Way tidal field for more than a single orbit (roughly 0.4 Gyr), which strongly suggests that the system is stabilized by the presence of large amounts of dark matter. If UMa3/U1 formed at the centre of a ~$10^9\rm M_\odot$ cuspy LCDM halo, its velocity dispersion would be predicted to be of order ~1 km/s. This is roughly consistent with the current estimate, which, neglecting binaries, places $\sigma_{\rm los}$ in the range 1 to 4 km/s. Because of its dense cusp, such a halo should be able to survive the Milky Way tidal field, keeping UMa3/U1 relatively unscathed until the present time. This implies that UMa3/U1 is in all likelihood the faintest and densest dwarf galaxy satellite of the Milky Way, with important implications for alternative dark matter models, and for the minimum halo mass threshold for luminous galaxy formation in the LCDM cosmology.

We present analysis of the mass-metallicity relation (MZR) for a sample of 67 [OIII]-selected star-forming galaxies at a redshift range of $z=1.99 - 2.32$ ($z_{\text{med}} = 2.16$) using \emph{Hubble Space Telescope} Wide Field Camera 3 grism spectroscopy from the Quasar Sightline and Galaxy Evolution (QSAGE) survey. Metallicities were determined using empirical gas-phase metallicity calibrations based on the strong emission lines [OII]3727,3729, [OIII]4959,5007 and H$\beta$. Star-forming galaxies were identified, and distinguished from active-galactic nuclei, via Mass-Excitation diagrams. Using $z\sim0$ metallicity calibrations, we observe a negative offset in the $z=2.2$ MZR of $\approx -0.51$ dex in metallicity when compared to locally derived relationships, in agreement with previous literature analysis. A similar offset of $\approx -0.46$ dex in metallicity is found when using empirical metallicity calibrations that are suitable out to $z\sim5$, though our $z=2.2$ MZR, in this case, has a shallower slope. We find agreement between our MZR and those predicted from various galaxy evolution models and simulations. Additionally, we explore the extended fundamental metallicity relation (FMR) which includes an additional dependence on star formation rate (SFR). Our results consistently support the existence of the FMR, as well as revealing an offset of $0.28\pm0.04$ dex in metallicity compared to locally-derived relationships, consistent with previous studies at similar redshifts. We interpret the negative correlation with SFR at fixed mass, inferred from an FMR existing for our sample, as being caused by the efficient accretion of metal-poor gas fuelling SFR at cosmic noon.

The Event Horizon Telescope (EHT) Collaboration has recently published the first horizon-scale images of the supermassive black holes M87* and Sgr A* and provided some first information on the physical conditions in their vicinity. The comparison between the observations and the three-dimensional general-relativistic magnetohydrodynamic (GRMHD) simulations has enabled the EHT to set initial constraints on the properties of these black-hole spacetimes. However, accurately distinguishing the properties of the accretion flow from those of the spacetime, most notably, the black-hole mass and spin, remains challenging because of the degeneracies the emitted radiation suffers when varying the properties of the plasma and those of the spacetime. The next-generation EHT (ngEHT) observations are expected to remove some of these degeneracies by exploring the complex interplay between the disk-jet dynamics, which represents one of the most promising tools for extracting information on the black-hole spin. By using GRMHD simulations of magnetically arrested disks (MADs) and general-relativistic radiative-transfer (GRRT) calculations of the emitted radiation, we have studied the properties of the jet and the accretion-disk dynamics on spatial scales that are comparable with the horizon. In this way, we are able to highlight that the radial and azimuthal dynamics of the jet are well correlated with the black-hole spin. Based on the resolution and image reconstruction capabilities of the ngEHT observations of M87*, we can assess the detectability and associated uncertainty of this correlation. Overall, our results serve to assess what are the prospects for constraining the black-hole spin with future EHT observations.

Physically motivated measurements are crucial for understanding galaxy growth and the role of the environment on their evolution. In particular, the growth of galaxies as measured by their size or radial extent provides an empirical approach for addressing this issue. However, the established definitions of galaxy size used for nearly a century are ill-suited for these studies because of a previously ignored bias. The conventionally-measured radii consistently miss the diffuse, outer extensions of stellar emission which harbour key signatures of galaxy growth, including star formation and gas accretion or removal. This issue is addressed by examining low surface brightness truncations or galaxy "edges" as a physically motivated tracer of size based on star formation thresholds. Our total sample consists of $\sim900$ galaxies with stellar masses ranging from $10^5 M_{\odot} < M_{\star} < 10^{11} M_{\odot}$. This sample of nearby cluster, group satellite and nearly isolated field galaxies was compiled using multi-band imaging from the Fornax Deep Survey, deep IAC Stripe 82 and Dark Energy Camera Legacy Surveys. Across the full mass range studied, we find that compared to the field, the edges of galaxies in the Fornax Cluster are located at 50% smaller radii and the average stellar surface density at the edges are two times higher. These results are consistent with the rapid removal of loosely bound neutral hydrogen (HI) in hot, crowded environments which truncates galaxies outside-in earlier, preventing the formation of more extended sizes and lower density edges. In fact, we find that galaxies with lower HI fractions have edges with higher stellar surface density. Our results highlight the importance of deep imaging surveys to study the low surface brightness imprints of the large scale structure and environment on galaxy evolution.

We present the discovery of Ursa Major III/UNIONS 1, the least luminous known satellite of the Milky Way, which is estimated to have an absolute V-band magnitude of $+2.2^{+0.4}_{-0.3}$ mag, equivalent to a total stellar mass of 16$^{+6}_{-5}$ M$_{\odot}$. Ursa Major III/UNIONS 1 was uncovered in the deep, wide-field Ultraviolet Near Infrared Optical Northern Survey (UNIONS) and is consistent with an old ($\tau > 11$ Gyr), metal-poor ([Fe/H] $\sim -2.2$) stellar population at a heliocentric distance of $\sim$ 10 kpc. Despite being compact ($r_{\text{h}} = 3\pm1$ pc) and composed of so few stars, we confirm the reality of Ursa Major III/UNIONS 1 with Keck II/DEIMOS follow-up spectroscopy and identify 11 radial velocity members, 8 of which have full astrometric data from $Gaia$ and are co-moving based on their proper motions. Based on these 11 radial velocity members, we derive an intrinsic velocity dispersion of $3.7^{+1.4}_{-1.0}$ km s$^{-1}$ but some caveats preclude this value from being interpreted as a direct indicator of the underlying gravitational potential at this time. Primarily, the exclusion of the largest velocity outlier from the member list drops the velocity dispersion to $1.9^{+1.4}_{-1.1}$ km s$^{-1}$, and the subsequent removal of an additional outlier star produces an unresolved velocity dispersion. While the presence of binary stars may be inflating the measurement, the possibility of a significant velocity dispersion makes Ursa Major III/UNIONS 1 a high priority candidate for multi-epoch spectroscopic follow-ups to deduce to true nature of this incredibly faint satellite.

Understanding surface temperature is important for habitability. Recent work on Mars has found that the dependence of surface temperature on elevation (surface lapse rate) converges to zero in the limit of a thin CO2 atmosphere. However, the mechanisms that control the surface lapse rate are still not fully understood. It remains unclear how the surface lapse rate depends on both greenhouse effect and surface pressure. Here, we use climate models to study when and why "mountaintops are cold". We find the tropical surface lapse rate increases with the greenhouse effect and with surface pressure. The greenhouse effect dominates the surface lapse rate transition and is robust across latitudes. The pressure effect is important at low latitudes in moderately opaque atmospheres. A simple model provides insights into the mechanisms of the transition. Our results suggest that topographic cold-trapping may be important for the climate of arid planets.

RCB stars are $L\approx10^4\,L_{\odot}$ solar-mass objects that can exhibit large periods of extinction from dust ejection episodes. Many exhibit semiregular pulsations in the range of $30-50$ days with semi-amplitudes of $0.05-0.3$ magnitude. Space-based photometry has discovered that solar-like oscillations are ubiquitous in hydrogen-dominated stars that have substantial outer convective envelopes, so we explore the hypothesis that the pulsations in RCB stars and the closely related dustless hydrogen-deficient carbon (dLHdC) stars, which have large convective outer envelopes of nearly pure helium, have a similar origin. Through stellar modeling and pulsation calculations, we find that the observed periods and amplitudes of these pulsations follows the well-measured phenomenology of their H-rich brethren. In particular, we show that the observed modes are likely of angular orders $l=0,1$ and $2$ and predominantly of an acoustic nature (i.e. $p$-modes with low radial order). The modes with largest amplitude are near the acoustic cut-off frequency appropriately rescaled to the helium-dominated envelope, and the observed amplitudes are consistent with that seen in high luminosity ($L>10^3\,L_{\odot}$) H-rich giants. We also find that for $T_{\mathrm{eff}}\gtrsim5400\,\mathrm{K}$, an HdC stellar model exhibits a radiative layer between two outer convective zones, creating a $g$-mode cavity that supports much longer period ($\approx 100$ days) oscillations. Our initial work was focused primarily on the adiabatic modes, but we expect that subsequent space-based observations of these targets (e.g. with TESS or Plato) are likely to lead to a larger set of detected frequencies that would allow for a deeper study of the interiors of these rare stars.

With hydrodynamical simulations we examine the evolution of a highly misaligned circumbinary disc around a black hole binary including the effects of general relativity. We show that a disc mass of just a few percent of the binary mass can significantly increase the binary eccentricity through von-Zeipel--Kozai-Lidov (ZKL) like oscillations provided that the disc lifetime is longer than the ZKL oscillation timescale. The disc begins as a relatively narrow ring of material far from the binary and spreads radially. When the binary becomes highly eccentric, disc breaking forms an inner disc ring that quickly aligns to polar. The polar ring drives fast retrograde apsidal precession of the binary that weakens the ZKL effect. This allows the binary eccentricity to remain at a high level and may significantly shorten the black hole merger time. The mechanism requires the initial disc inclination relative to the binary to be closer to retrograde than to prograde.

We announce the C23.01 update of Cloudy. This corrects a simple coding error, present since $\sim$ 1990, in one routine that required a conversion from the line-center to the mean normalization of the Ly$\alpha$ optical depth. This affects the destruction of H I Ly$\alpha$ by background opacities. Its largest effect is upon the Ly$\alpha$ intensity in high-ionization dusty clouds, where the predicted intensity is now up to three times stronger. Other properties that depend on Ly$\alpha$ destruction, such as grain infrared emission, change in response.

Gamma-ray burst GRB221009A was of unprecedented brightness in gamma-rays and X-rays, and through to the far ultraviolet, allowing for identification within a host galaxy at redshift z=0.151 by multiple space and ground-based optical/near-infrared telescopes and enabling a first association - via cosmic-ray air-shower events - with a photon of 251 TeV. That is in direct tension with a potentially observable phenomenon of quantum gravity (QG), where spacetime "foaminess" accumulates in wavefronts propagating cosmological distances, and at high-enough energy could render distant yet bright pointlike objects invisible, by effectively spreading their photons out over the whole sky. But this effect would not result in photon loss, so it remains distinct from any absorption by extragalactic background light. A simple multiwavelength average of foam-induced blurring is described, analogous to atmospheric seeing from the ground. When scaled within the fields of view for the Fermi and Swift instruments, it fits all z<5 GRB angular-resolution data of 10 MeV or any lesser peak energy and can still be consistent with the highest-energy localization of GRB221009A: a limiting bound of about 1 degree is in agreement with a holographic QG-favored formulation.

Context: The measurement of the Yarkovsky effect on near-Earth asteroids (NEAs) is common practice in orbit determination today, and the number of detections will increase with the developments of new and more accurate telescopic surveys. However, the process of finding new detections and identifying spurious ones is not yet automated, and it often relies on personal judgment. Aims: We aim to introduce a more automated procedure that can search for NEA candidates to measure the Yarkovsky effect, and that can identify spurious detections. Methods: The expected semi-major axis drift on an NEA caused by the Yarkovsky effect was computed with a Monte Carlo method on a statistical model of the physical parameters of the asteroid that relies on the most recent NEA population models and data. The expected drift was used to select candidates in which the Yarkovsky effect might be detected, according to the current knowledge of their orbit and the length of their observational arc. Then, a nongravitational acceleration along the transverse direction was estimated through orbit determination for each candidate. If the detected acceleration was statistically significant, we performed a statistical test to determine whether it was compatible with the Yarkovsky effect model. Finally, we determined the dependence on an isolated tracklet. Results: Among the known NEAs, our procedure automatically found 348 detections of the Yarkovsky effect that were accepted. The results are overall compatible with the predicted trend with the the inverse of the diameter, and the procedure appears to be efficient in identifying and rejecting spurious detections. This algorithm is now adopted by the ESA NEO Coordination Centre to periodically update the catalogue of NEAs with a measurable Yarkovsky effect, and the results are automatically posted on the web portal.

In the context of this paper, microlenses present as oblate clusters of dark matter structures called massive astrophysical compact halo object(MACHO) in the galactic halo are considered. The NFW density profile [1] is derived from the observational data and works best in the halo region of the exterior part of the galaxy. Hence this profile is used to plot the potential, deflection angle, and critical and caustic curves for the aforementioned microlenses using numerical methods. Moreover, this model is compared with an older density profile model [2], and the differences in their caustic and critical curves are pointed out. This leads to the conclusion that the NFW model produces caustic and critical curves that occupy a smaller region than portrayed by the caustic and critical curves produced by the older density profile model. However, the differences are not that significant to these structures act as microlenses as they are so small and beyond the scope of modern telescopes, and thus only the light curves can be detected from these structures.

In recent years, certain luminous extragalactic optical transients have been observed to last only a few days. Their short observed duration implies a different powering mechanism from the most common luminous extragalactic transients (supernovae) whose timescale is weeks. Some short-duration transients, most notably AT2018cow, display blue optical colours and bright radio and X-ray emission. Several AT2018cow-like transients have shown hints of a long-lived embedded energy source, such as X-ray variability, prolonged ultraviolet emission, a tentative X-ray quasiperiodic oscillation, and large energies coupled to fast (but subrelativistic) radio-emitting ejecta. Here we report observations of minutes-duration optical flares in the aftermath of an AT2018cow-like transient, AT2022tsd (the "Tasmanian Devil"). The flares occur over a period of months, are highly energetic, and are likely nonthermal, implying that they arise from a near-relativistic outflow or jet. Our observations confirm that in some AT2018cow-like transients the embedded energy source is a compact object, either a magnetar or an accreting black hole.

We recently presented the first 3D numerical simulation of the solar interior for which tachocline confinement was achieved by a dynamo-generated magnetic field. In this followup study, we analyze the degree of confinement as the magnetic field strength changes (controlled by varying the magnetic Prandtl number) in a coupled radiative zone (RZ) and convection zone (CZ) system. We broadly find three solution regimes, corresponding to weak, medium, and strong dynamo magnetic field strengths. In the weak-field regime, the large-scale magnetic field is mostly axisymmetric with regular, periodic polarity reversals (reminiscent of the observed solar cycle), but fails to create a confined tachocline. In the strong-field regime, the large-scale field is mostly non-axisymmetric with irregular, quasi-periodic polarity reversals, and creates a confined tachocline. In the medium-field regime, the large-scale field resembles a strong-field dynamo for extended intervals, but intermittently weakens to allow temporary epochs of strong differential rotation. In all regimes, the amplitude of poloidal field strength in the RZ is very well explained by skin-depth arguments, wherein the oscillating field that gives rise to the skin depth (in the medium- and strong-field cases) is a non-axisymmetric field structure rotating with respect to the RZ. These new simulations reaffirm that tachocline confinement by the solar dynamo (the so-called fast magnetic confinement scenario) is possible, but suggest a new picture in which non-axisymmetric field components rotating with respect to the RZ play the primary role, instead of the regularly reversing axisymmetic field associated with the 22-year cycle.

The observation of GW170817, the first binary neutron star merger observed in both gravitational waves (GW) and electromagnetic (EM) waves, kickstarted the age of multi-messenger GW astronomy. This new technique presents an observationally rich way to probe extreme astrophysical processes. With the onset of the LIGO-Virgo-KAGRA Collaboration's O4 observing run and wide-field EM instruments well-suited for transient searches, multi-messenger astrophysics has never been so promising. We review recent searches and results for multi-messenger counterparts to GW events, and describe existing and upcoming EM follow-up facilities, with a particular focus on WINTER, a new near-infrared survey telescope, and TESS, an exoplanet survey space telescope.

The latest generation of high-resolution spectrograph instruments on 10m-class telescopes continue to pursue challenging science cases. Consequently, ever more precise calibration methods are necessary to enable trail-blazing science methodology. We present the High-resolution Infrared SPectrograph for Exoplanet Characterization (HISPEC) Calibration Unit (CAL), designed to facilitate challenging science cases such as Doppler imaging of exoplanet atmospheres, precision radial velocity, and high-contrast high-resolution spectroscopy of nearby exoplanets. CAL builds upon the heritage from the pathfinder instrument Keck Planet Imager and Characterizer (KPIC) and utilizes four near-infrared (NIR) light sources encoded with wavelength information that are coupled into single-mode fibers. They can be used synchronously during science observations or asynchronously during daytime calibrations. A hollow cathode lamp (HCL) and a series of gas absorption cells provide absolute calibration from 0.98 {\mu}m to 2.5 {\mu}m. A laser frequency comb (LFC) provides stable, time-independent wavelength information during observation and CAL implements a lower finesse astro-etalon as a backup for the LFC. Design lessons from instrumentation like HISPEC will serve to inform the requirements for similar instruments for the ELTs in the future.

Galaxy gas kinematics are sensitive to the physical processes that contribute to a galaxy's evolution. It is expected that external processes will cause more significant kinematic disturbances in the outer regions, while internal processes will cause more disturbances for the inner regions. Using a subsample of 47 galaxies ($0.27<z<0.36$) from the Middle Ages Galaxy Properties with Integral Field Spectroscopy (MAGPI) survey, we conduct a study into the source of kinematic disturbances by measuring the asymmetry present in the ionised gas line-of-sight velocity maps at the $0.5R_e$ (inner regions) and $1.5R_e$ (outer regions) elliptical annuli. By comparing the inner and outer kinematic asymmetries, we aim to better understand what physical processes are driving the asymmetries in galaxies. We find the local environment plays a role in kinematic disturbance, in agreement with other integral field spectroscopy studies of the local universe, with most asymmetric systems being in close proximity to a more massive neighbour. We do not find evidence suggesting that hosting an Active Galactic Nucleus (AGN) contributes to asymmetry within the inner regions, with some caveats due to emission line modelling. In contrast to previous studies, we do not find evidence that processes leading to asymmetry also enhance star formation in MAGPI galaxies. Finally, we find a weak anti-correlation between stellar mass and asymmetry (ie. high stellar mass galaxies are less asymmetric). We conclude by discussing possible sources driving the asymmetry in the ionised gas, such as disturbances being present in the colder gas phase (either molecular or atomic) prior to the gas being ionised, and non-axisymmetric features (e.g., a bar) being present in the galactic disk. Our results highlight the complex interplay between ionised gas kinematic disturbances and physical processes involved in galaxy evolution.

Ionized nebulae are key to understanding the chemical composition and evolution of the Universe. Among these nebulae, H~{\sc ii} regions and planetary nebulae are particularly important as they provide insights into the present and past chemical composition of the interstellar medium, along with the nucleosynthetic processes involved in the chemical evolution of the gas. However, the heavy-element abundances derived from collisional excited lines (CELs) and recombination lines (RLs) do not align. This longstanding abundance-discrepancy problem calls into question our absolute abundance determinations. Which of the lines (if any) provides the correct heavy-element abundances? Recently, it has been shown that there are temperature inhomogeneities concentrated within the highly ionized gas of the H~{\sc ii} regions, causing the reported discrepancy. However, planetary nebulae do not exhibit the same trends as the H~{\sc ii} regions, suggesting a different origin for the abundance discrepancy. In this proceedings, we briefly discuss the state-of-the-art of the abundance discrepancy problem in both H~{\sc ii} regions and planetary nebulae.

Millimter (mm) frequencies are primarily sensitive to thermal emission from layers across the stellar chromosphere up to the transition region, while metrewave (radio) frequencies probe the coronal heights. Together the mm and radio band spectroscopic snapshot imaging enables the tomographic exploration of the active atmospheric layers of the cool main-sequence stars (spectral type: FGKM), including our Sun. Sensitive modern mm and radio interferometers let us explore solar/stellar activity covering a range of energy scales at sub-second and sub-MHz resolution over wide operational bandwidths. The superior uv-coverage of these instruments facilitate high dynamic range imaging, letting us explore the morphological evolution of even energetically weak events on the Sun at fine spectro-temporal cadence. This article will introduce the current advancements, the data analysis challenges and available tools. The impact of these tools and novel data in field of solar/stellar research will be summarised with future prospects.

We propose a pebble-driven core accretion scenario to explain the formation of giant planets around the late-M dwarfs of $M_{\star}{=}0.1{-}0.2 \ M_{\odot}$. In order to explore the optimal disk conditions for giant planet, we perform N-body simulations to investigate the growth and dynamical evolution of both single and multiple protoplanets in the disks with both inner viscously heated and outer stellar irradiated regions. The initial masses of the protoplanets are either assumed to be equal to $0.01 \ M_{\oplus}$ or calculated based on the formula derived from streaming instability simulations. Our findings indicate that massive planets are more likely to form in disks with longer lifetimes, higher solid masses, moderate to high levels of disk turbulence, and larger initial masses of protoplanets. In the single protoplanet growth cases, the highest planet core mass that can be reached is generally lower than the threshold necessary to trigger rapid gas accretion, which impedes the formation of giant planets. Nonetheless, in multi-protoplanet cases, the cores can exceed the pebble isolation mass barrier aided by frequent planet-planet collisions. This consequently speeds up their gas accretion and promotes giant planet formation, making the optimal parameter space to grow giant planets substantially wider. Taken together, our results suggest that even around very low-mass stellar hosts, the giant planets with orbital periods of ${\lesssim}100$ days are still likely to form when lunar-mass protoplanets first emerge from planetesimal accretion and then grow rapidly by a combination of pebble accretion and planet-planet collisions in disks with a high supply of pebble reservoir ${>}50 \ M_{\oplus}$ and turbulent level of $\alpha_{\rm t} {\sim} 10^{-3}{-}10^{-2}$.

Heavy (Z>26) solar energetic particles (SEPs) with energies ~1 MeV/nucleon are known to leave visible damage tracks in meteoritic materials. The density of such solar flare tracks in lunar and asteroidal samples has been used as a measure of a sample's exposure time to space, yielding critical information on planetary space weathering rates, the dynamics and lifetimes of interplanetary dust grains, and the long-term history of solar particle fluxes. Knowledge of the SEP track accumulation rate in planetary materials at 1 au is critical for properly interpreting observed track densities. Here, we use in-situ particle observations of the 0.50-3.0 MeV/nuc Fe-group SEP flux taken by NASA's Advanced Composition Explorer (ACE) to calculate a flux of track-inducing particles at 1 au of 6.0x10^5 /cm2/yr/str. Using the observed energy spectrum of Fe-group SEPs, we find that the depth distribution of SEP-induced damage tracks inferred from ACE measurements matches closely to that recently measured in lunar sample 64455; however, the magnitude of the ACE-inferred rate is approximately 25x higher than that observed in the lunar sample. We discuss several hypotheses for the nature of this discrepancy, including inefficiencies in track formation, thermal annealing of lunar samples, erosion via space weathering processing, and variations in the SEP flux at the Moon, yet find no satisfactory explanation. We encourage further research on both the nature of SEP track formation in meteoritic materials and the flux of Fe-group SEPs at the lunar surface in recent and geologic times to resolve this discrepancy.

To ensure that magnetohydrodynamical (MHD) turbulence simulations accurately reflect the physics, it is critical to understand numerical dissipation. Here we determine the hydrodynamic and magnetic Reynolds number (Re and Rm) as a function of linear grid resolution N, in MHD simulations with purely numerical viscosity and resistivity (implicit large eddy simulations; ILES). We quantify the numerical viscosity in the subsonic and supersonic regime, via simulations with sonic Mach numbers of Mach=0.1 and Mach=10, respectively. We find Re=(N/N_Re)^p_Re, with p_Re=[1.2,1.4] and N_Re=[0.8,1.7] for Mach=0.1, and p_Re=[1.5,2.0] and N_Re=[0.8,4.4] for Mach=10, and we find Rm=(N/N_Rm)^p_Rm, with p_Rm=[1.3,1.5] and N_Rm=[1.1,2.3] for Mach=0.1, and p_Rm=[1.2,1.6] and N_Rm=[0.1,0.7] for Mach=10. The resulting magnetic Prandtl number (Pm=Rm/Re) is consistent with a constant value of Pm=1.3+/-1.1 for Mach=0.1, and Pm=2.0+/-1.4 for Mach=10. We compare our results with an independent study in the subsonic regime and find excellent agreement in p_Re and p_Rm, and agreement within a factor of ~2 for N_Re and N_Rm (due to differences in the codes and solvers). We compare these results to the target Re and Rm set in direct numerical simulations (DNS, i.e., using explicit viscosity and resistivity) from the literature. This comparison and our ILES relations can be used to determine whether a target Re and Rm can be achieved in a DNS for a given N. We conclude that for the explicit (physical) dissipation to dominate over the numerical dissipation, the target Reynolds numbers must be set lower than the corresponding numerical values.

Conventional galaxy mass estimation methods suffer from model assumptions and degeneracies. Machine learning, which reduces the reliance on such assumptions, can be used to determine how well present-day observations can yield predictions for the distributions of stellar and dark matter. In this work, we use a general sample of galaxies from the TNG100 simulation to investigate the ability of multi-branch convolutional neural network (CNN) based machine learning methods to predict the central (i.e., within $1-2$ effective radii) stellar and total masses, and the stellar mass-to-light ratio $M_*/L$. These models take galaxy images and spatially-resolved mean velocity and velocity dispersion maps as inputs. Such CNN-based models can in general break the degeneracy between baryonic and dark matter in the sense that the model can make reliable predictions on the individual contributions of each component. For example, with $r$-band images and two galaxy kinematic maps as inputs, our model predicting $M_*/L$ has a prediction uncertainty of 0.04 dex. Moreover, to investigate which (global) features significantly contribute to the correct predictions of the properties above, we utilize a gradient boosting machine. We find that galaxy luminosity dominates the prediction of all masses in the central regions, with stellar velocity dispersion coming next. We also investigate the main contributing features when predicting stellar and dark matter mass fractions ($f_*$, $f_{\rm DM}$) and the dark matter mass $M_{DM}$, and discuss the underlying astrophysics.

We present an extensive recent multi-band optical photometric observations of the blazar S5 0716+714 carried out over 53 nights with two telescopes in India, two in Bulgaria, one in Serbia, and one in Egypt during 2019 November -- 2022 December. We collected 1401, 689, 14726, and 165 photometric image frames in B, V, R, and I bands, respectively. We montiored the blazar quasi-simultaneously during 3 nights in B, V, R, and I bands; 4 nights in B, V, and R; 2 nights in V, R, and I; 5 nights in B and R; and 2 nights in V and R bands. We also took 37 nights of data only in R band. Single band data are used to study intraday flux variability and two or more bands quasi-simultaneous observations allow us to search for colour variation in the source. We employ the power-enhanced F-test and the nested ANOVA test to search for genuine flux and color variations in the light curves of the blazar on intraday timescales. Out of 12, 11, 53, and 5 nights observations, intraday variations with amplitudes between ~3% and ~20% are detected in 9, 8, 31 and 3 nights in B, V, R, and I bands, respectively, corresponding to duty cycles of 75%, 73%, 58% and 60%. These duty cycles are lower than those typically measured at earlier times. On these timescales color variations with both bluer-when-brighter and redder-when-brighter are seen, though nights with no measurable colour variation are also present. We briefly discuss possible explanations for this observed intraday variability.

Lyman-Werner (LW) radiation photodissociating molecular hydrogen (H$_2$) influences the thermal and dynamical evolution of the Population III (Pop III) star-forming gas cloud. The effect of powerful LW radiation has been well investigated in the context of supermassive black hole formation in the early universe. However, the average intensity in the early universe is several orders of magnitude lower. For a comprehensive study, we investigate the effects of LW radiation at $18$ different intensities, ranging from $J_{\rm LW}/J_{21}=0$ (no radiation) to $30$ (H-cooling cloud), on the primordial star-forming gas cloud obtained from a three-dimensional cosmological simulation. The overall trend with increasing radiation intensity is a gradual increase in the gas cloud temperature, consistent with previous works. Due to the HD-cooling, on the other hand, the dependence of gas cloud temperature on $J_{\rm LW}$ deviates from the aforementioned increasing trend for a specific range of intensities ($J_{\rm LW}/J_{21}=0.025-0.09$). In HD-cooling clouds, the temperature remained below $200$ K during $10^5$ yr after the first formation of the high-density region, maintaining a low accretion rate. Finally, the HD-cooling clouds have only a low-mass dense core (above $10^8\,{\rm cm^{-3}}$) with about $1-16\, M_{\odot}$, inside which a low-mass Pop III star with $\leq\!0.8\,M_{\odot}$ (so-called "surviving star") could form. The upper limit of star formation efficiency $M_{\rm core}/M_{\rm vir, gas}$ significantly decreases from $10^{-3}$ to $10^{-5}$ as HD-cooling becomes effective. This tendency indicates that, whereas the total gas mass in the host halo increases with the LW radiation intensity, the total Pop III stellar mass does not increase similarly.

We present XMM-Newton observations of a representative X-ray selected sample of 31 galaxy clusters at moderate redshift $(0.4<z<0.6)$, spanning the mass range $10^{14} < M_{\textrm 500} < 10^{15}$~M$_\odot$. This sample, EXCPRES (Evolution of X-ray galaxy Cluster Properties in a Representative Sample), is used to test and validate a new method to produce optimally-binned cluster X-ray temperature profiles. The method uses a dynamic programming algorithm, based on partitioning of the soft-band X-ray surface brightness profile, to obtain a binning scheme that optimally fulfills a given signal-to-noise threshold criterion out to large radius. From the resulting optimally-binned EXCPRES temperature profiles, and combining with those from the local REXCESS sample, we provide a generic scaling relation between the relative error on the temperature and the [0.3-2] keV surface brightness signal-to-noise ratio, and its dependence on temperature and redshift. We derive an average scaled 3D temperature profile for the sample. Comparing to the average scaled 3D temperature profiles from REXCESS, we find no evidence for evolution of the average profile shape within the redshift range that we probe.

The IceCube Neutrino Observatory is uniquely sensitive to MeV neutrinos emitted during a core-collapse supernova. The Supernova Data Acquisition System (SNDAQ) monitors in real-time the detector rate deviation searching for bursts of MeV neutrinos. We present a new analysis stream that uses SNDAQ to respond to external alerts from gravitational waves detected in LIGO-Virgo-KAGRA.

Supernova (SN) 2023ixf was discovered on May 19th, 2023. The host galaxy, M101, was observed by the Hobby Eberly Telescope Dark Energy Experiment (HETDEX) collaboration over the period April 30, 2020 -- July 10, 2020, using the Visible Integral-field Replicable Unit Spectrograph (VIRUS; $3470\lesssim\lambda\lesssim5540$ \r{A}) on the 10-m Hobby-Eberly Telescope (HET). The fiber filling factor within $\pm$ 30 arcsec of SN 2023ixf is 80% with a spatial resolution of 1 arcsec. The r<5.5 arcsec surroundings are 100% covered. This allows us to analyze the spatially resolved pre-explosion local environments of SN 2023ixf with nebular emission lines. The 2-dimensional (2D) maps of the extinction and the star-formation rate (SFR) surface density ($\Sigma_{\rm SFR}$) show weak increasing trends in the radial distributions within the r<5.5 arcsec regions, suggesting lower values of extinction and SFR in the vicinity of the progenitor of SN 2023ixf. The median extinction and that of the surface density of SFR within r<3 arcsec are $E(B-V)=0.06\pm0.14$, and $\Sigma_{\rm SFR}=10^{-5.44\pm0.66}~\rm M_{\odot}\cdot yr^{-1}\cdot arcsec^{-2}$. There is no significant change in extinction before and after the explosion. The gas metallicity does not change significantly with the separation from SN 2023ixf. The metal-rich branch of the $R_{23}$ calculations indicates that the gas metallicity around SN 2023ixf is similar to the solar metallicity ($\sim Z_{\odot}$). The archival deep images from the Canada-France-Hawaii Telescope Legacy Survey (CFHTLS) show a clear detection of the progenitor of SN 2023ixf in the $z$-band at $22.778\pm0.063$ mag, but non-detections in the remaining four bands of CFHTLS ($u,g,r,i$). The results suggest a massive progenitor of $\approx$ 22 $M_\odot$.

Various types of inertial modes have been observed and identified on the Sun, including the equatorial Rossby modes, critical-latitude modes, and high-latitude modes. Recent observations further report a detection of equatorially-antisymmetric radial vorticity modes which propagate in a retrograde direction about three times faster than those of the equatorial Rossby modes when seen in the corotating frame with the Sun. Here, we study the properties of these equatorially-antisymmetric vorticity modes using a realistic linear model of the Sun's convection zone. We find that they are essentially non-toroidal, involving a substantial radial flow at the equator. Thus, the background density stratification plays a critical role in determining their dispersion relation. The solar differential rotation is also found to have a significant impact by introducing the viscous critical layers and confining the modes near the base of the convection zone. Furthermore, we find that their propagation frequencies are strikingly sensitive to the background superadiabaticity $\delta$ because the buoyancy force acts as an additional restoring force for these non-toroidal modes. The observed frequencies are compatible with the linear model only when the bulk of the convection zone is weakly subadiabatic ($-5\times 10^{-7} \lesssim \delta \lesssim -2.5\times 10^{-7}$). Our result is consistent with but tighter than the constraint independently derived in a previous study ($\delta<2 \times 10^{-7}$) employing the high-latitude inertial mode. It is implied that, below the strongly superadiabatic near-surface layer, the bulk of the Sun's convection zone might be much closer to adiabatic than typically assumed and may even be weakly subadiabatic.

Simulation based inference has seen increasing interest in the past few years as a promising approach to model the non linear scales of galaxy clustering. The common approach using Gaussian process is to train an emulator over the cosmological and galaxy-halo connection parameters independently for every scales. We present a new Gaussian process model allowing to extend the input parameter space dimensions and to use non-diagonal noise covariance matrix. We use our new framework to emulate simultaneously every scales of the non-linear clustering of galaxies in redshift space from the AbacusSummit N-body simulations at redshift $z=0.2$. The model includes nine cosmological parameters, five halo occupation distribution (HOD) parameters and one scale dimension. Accounting for the limited resolution of the simulations, we train our emulator on scales from $0.3~h^{-1}\mathrm{Mpc}$ to $60~h^{-1}\mathrm{Mpc}$ and compare its performance with the standard approach of building one independent emulator for each scales. The new model yields more accurate and precise constraints on cosmological parameters compared to the standard approach. As the new model is able to interpolate over the scales space, we are also able to account for the Alcock-Paczynski distortion effect leading more accurate constraints on the cosmological parameters.

Extrasolar debris disks are detected by observing dust, which is thought to be released during planetesimal collisions. This implies that planetesimals are dynamically excited ("stirred"), such that collisions are sufficiently common and violent. The most frequently considered stirring mechanisms are self-stirring by disk self-gravity, and planet-stirring via secular interactions. However, these models face problems when considering disk mass, self-gravity, and planet eccentricity, leading to the possibility that other, unexplored mechanisms instead stir debris. We hypothesize that planet-stirring could be more efficient than the traditional secular model implies, due to two additional mechanisms. First, a planet at the inner edge of a debris disk can scatter massive bodies onto eccentric, disk-crossing orbits, which then excite debris ("projectile stirring"). Second, a planet can stir debris over a wide region via broad mean-motion resonances, both at and between nominal resonance locations ("resonant stirring"). Both mechanisms can be effective even for low-eccentricity planets, unlike secular-planet-stirring. We run N-body simulations across a broad parameter space, to determine the viability of these new stirring mechanisms. We quantify stirring levels using a bespoke program for assessing Rebound debris simulations, which we make publicly available. We find that even low-mass projectiles can stir disks, and verify this with a simple analytic criterion. We also show that resonant stirring is effective for planets above ~0.5 MJup. By proving that these mechanisms can increase planet-stirring efficiency, we demonstrate that planets could still be stirring debris disks even in cases where conventional (secular) planet-stirring is insufficient.

We use the GALFORM semi-analytical galaxy formation model implemented in the Planck Millennium N-body simulation to build a mock galaxy catalogue on an observer's past lightcone. The mass resolution of this N-body simulation is almost an order of magnitude better than in previous simulations used for this purpose, allowing us to probe fainter galaxies and hence build a more complete mock catalogue at low redshifts. The high time cadence of the simulation outputs allows us to make improved calculations of galaxy properties and positions in the mock. We test the predictions of the mock against the Physics of the Accelerating Universe Survey, a narrow band imaging survey with highly accurate and precise photometric redshifts, which probes the galaxy population over a lookback time of 8 billion years. We compare the model against the observed number counts, redshift distribution and evolution of the observed colours and find good agreement; these statistics avoid the need for model-dependent processing of the observations. The model produces red and blue populations that have similar median colours to the observations. However, the bimodality of galaxy colours in the model is stronger than in the observations. This bimodality is reduced on including a simple model for errors in the GALFORM photometry. We examine how the model predictions for the observed galaxy colours change when perturbing key model parameters. This exercise shows that the median colours and relative abundance of red and blue galaxies provide a constraint on the strength of the feedback driven by supernovae used in the model.

Comets are pristine remnants of the Solar system, composed of dust and ice. They remain inactive and undetectable for most of their orbit due to low temperatures. However, as they approach the Sun, volatile materials sublimate, expelling dust and creating a visible coma. Spectroscopic observations of comets help the simultaneous study of both the gas emissions and reflected sunlight from dust particles. By implementing a long slit, the spatial variations in molecular emissions can be analysed to be further used for other computations. Additionally, spatial information aids in extracting the characteristic profile of the Af(rho) parameter, revealing insights into the behaviour of dust emissions. A sufficiently long slit would prove advantageous in extracting information about the emissions occurring at different parts of the coma or even the tail. We can gain an overall comprehensive understanding of a comet's chemical composition and dust emission by constructively utilising low-resolution spectroscopy with the help of a long slit.

At present, the field of astronomical machine learning lacks widely-used benchmarking datasets; most research employs custom-made datasets which are often not publicly released, making comparisons between models difficult. In this paper we present CRUMB, a publicly-available image dataset of Fanaroff-Riley galaxies constructed from four "parent" datasets extant in the literature. In addition to providing the largest image dataset of these galaxies, CRUMB uses a two-tier labelling system: a "basic" label for classification and a "complete" label which provides the original class labels used in the four parent datasets, allowing for disagreements in an image's class between different datasets to be preserved and selective access to sources from any desired combination of the parent datasets.

The disk of an active galactic nucleus (AGN) is widely regarded as a prominent formation channel of binary black hole (BBH) mergers that can be detected through gravitational waves (GWs). Besides, the presence of dense environmental gas offers the potential for an embedded BBH merger to produce electromagnetic (EM) counterparts. In this paper, we investigate EM emission powered by the kicked remnant of a BBH merger occurring within the AGN disk. The remnant BH will launch a jet via accreting magnetized medium as it traverses the disk. The resulting jet will decelerate and dissipate energy into a lateral cocoon during its propagation. We explore three radiation mechanisms of the jet cocoon system: jet breakout emission, disk cocoon cooling emission, and jet cocoon cooling emission, and find that the jet cocoon cooling emission is more likely to be detected in its own frequency bands. We predict a soft X-ray transient, lasting for O($10^3$) s, to serve as an EM counterpart, of which the time delay O(10) days after the GW trigger contributes to follow-up observations. Consequently, BBH mergers in the AGN disk represent a novel multimessenger source. In the future, enhanced precision in measuring and localizing GWs, coupled with diligent searches for such associated EM signal, will effectively validate or restrict the origin of BBH mergers in the AGN disk.

Current observations do not rule out the possibility that the Universe might end up in an abrupt event. Different such scenarios may be explored through suitable parameterizations of the dark energy and then confronted to cosmological background data. Here we parameterize a pseudorip scenario using a particular sigmoid function and carry an in-depth multifaceted examination of its evolutionary features and statistical performance. This depiction of a non violent final fate of our cosmos seems to be arguably statistically favoured over the consensus {\Lambda}CDM model according to some Bayesian discriminators.

In this work, we derive for the first time observational constraints on the extended Minimal Theory of Massive Gravity (eMTMG) framework in light of Planck-CMB data, geometrical measurements from Baryon Acoustic Oscillation (BAO), Type Ia supernovae from the recent Pantheon+ samples, and also using the auto and cross-correlations cosmic shear measurements from KIDS-1000 survey. Given the great freedom of dynamics choice for the theory, we consider an observationally motivated subclass in which the background evolution of the Universe goes through a transition from a (positive or negative) value of the effective cosmological constant to another value. From the statistical point of view, we did not find evidence of such a transition, i.e. deviation from the standard $\Lambda$CDM behavior, and from the joint analysis using Planck + BAO + Pantheon+ data, we constrain the graviton mass to $< 6.6 \times 10^{-34}$ eV at 95% CL. We use KIDS-1000 survey data to constrain the evolution of the scalar perturbations of the model and its limits for the growth of structure predicted by the eMTMG scenario. In this case, we find small evidence at 95% CL for a non-zero graviton mass. We interpret and discuss these results in light of the current tension on the $S_8$ parameter. We conclude that, within the subclass considered, the current data are only able to impose upper bounds on the eMTMG dynamics. Given its potentialities beyond the subclass, eMTMG can be classified as a good candidate for modified gravity, serving as a framework in which observational data can effectively constrain (or confirm) the graviton mass and deviations from the standard $\Lambda$CDM behavior.

The usual approach to classify the ionizing source using optical spectroscopy is based on the use of diagnostic diagrams that compares the relative strength of pairs of collisitional metallic lines (e.g., [O iii] and [N ii]) with respect to recombination hydrogen lines (e.g., H{\beta} and H{\alpha}). Despite of being accepted as the standard procedure, it present known problems, including confusion regimes and/or limitations related to the required signal-to-noise of the involved emission lines. These problems affect not only our intrinsic understanding of inter-stellar medium and its poroperties, but also fundamental galaxy properties, such as the star-formation rate and the oxygen abundance, and key questions just as the fraction of active galactic nuclei, among several others. We explore the existing alternatives in the literature to minimize the confusion among different ionizing sources and proposed a new simple diagram that uses the equivalent width and the velocity dispersion from one single emission line, H{\alpha}, to classify the ionizing sources. We use aperture limited and spatial resolved spectroscopic data in the nearby Universe (z{\sim}0.01) to demonstrate that the new diagram, that we called WHaD, segregates the different ionizing sources in a more efficient way that previously adopted procedures. A new set of regions are defined in this diagram to select betweeen different ionizing sources. The new proposed diagram is well placed to determine the ionizing source when only H{\alpha} is available, or when the signal-to-noise of the emission lines involved in the classical diagnostic diagrams (e.g., H{\beta}).

Electronically Assisted Astronomy consists in capturing deep sky images with a digital camera coupled to a telescope to display views of celestial objects that would have been invisible through direct observation. This practice generates a large quantity of data, which may then be enhanced with dedicated image editing software after observation sessions. In this study, we show how Image Quality Assessment can be useful for automatically rating astronomical images, and we also develop a dedicated model by using Automated Machine Learning.

Our understanding of the surface porosity of icy moons and its evolution with depth remains limited, including the precise scale at which ice compaction occurs under self-weight pressure. This parameter is of crucial interest for the correct interpretation of current remote sensing data (spectroscopy in the visible, infrared to passive microwave) but also for planetary exploration when designing a lander, a rover or a cryobot. In situ exploration of the ice crust would require knowledge about subsurface porosity. This study employs a compaction model solely driven by overburden pressure based on prior research. The formulation for density as a function of depth, incorporates an essential parameter: the ice compaction coefficient. To determine this coefficient, we fit our depth-dependent density model to existing data obtained from Earth-based measurements of ice cores in Antarctica and North Greenland. Our results yield a typical lengthscale for ice compaction on Earth of approximately 20.1 $\pm$ 0.6 m , consistent with the existing literature. We apply the model to Europa, which due to its lower gravity, has a typical ice compaction scale of 150 $\pm$ 4 m. We compare it with the depths scanned by current spaceborne data and find that porosity can be considered constant when accounting only for gravity-induced compaction.

We present LinearResponse.jl, an efficient, versatile public library written in julia to compute the linear response of self-gravitating (3D spherically symmetric) stellar spheres and (2D axisymmetric razor-thin) discs. LinearResponse.jl can scan the whole complex frequency plane, probing unstable, neutral and (weakly) damped modes. Given a potential model and a distribution function, this numerical toolbox estimates the modal frequencies as well as the shapes of individual modes. The libraries are validated against a combination of previous results for the spherical isochrone model and Mestel discs, and new simulations for the spherical Plummer model. Beyond linear response theory, the realm of applications of LinearResponse.jl also extends to the kinetic theory of self-gravitating systems through a modular interface.

Since the discovery of dark matter-deficient galaxies, numerous studies have shown that these exotic galaxies naturally occur in the $\Lambda$CDM model due to stronger tidal interactions. They are typically satellites, with stellar masses in the $10^8-10^9\;\mathrm{M}_\odot$ range, of more massive galaxies. The recent discovery of a massive galaxy lacking dark matter and also lacking a more massive neighbor is puzzling. Two possible scenarios have been suggested in the literature: either the galaxy lost its dark matter early or it had been lacking ab initio. As a proof of concept for the former assumption, we present an example from IllustrisTNG300. At present, the galaxy has a stellar mass of $M_\star \simeq 6.8 \cdot 10^9\; \mathrm{M}_\odot$, with no gas, $M_\mathrm{DM}/M_\mathrm{B} \simeq 1.31$, and a stellar half-mass radius of $R_{0.5,\star} = 2.45\;\mathrm{kpc}$. It lost the majority of its dark matter early, between $z = 2.32$ and $z = 1.53$. Since then, it has continued to dwell in the cluster environment, interacting with the cluster members without merging, while accelerating on its orbit. Eventually, it left the cluster and it has spent the last $\sim 2\;\mathrm{Gyr}$ in isolation, residing just outside the most massive cluster in the simulation. Thus, the galaxy represents the first example found in simulations of both an isolated dark matter-poor galaxy that lost its extended envelope early and a fairly compact stellar system that has managed to escape.

An ESO internal ALMA development study, BRAIN, is addressing the ill-posed inverse problem of synthesis image analysis employing astrostatistics and astroinformatics. These emerging fields of research offer interdisciplinary approaches at the intersection of observational astronomy, statistics, algorithm development, and data science. In this study, we provide evidence of the benefits of employing these approaches to ALMA imaging for operational and scientific purposes. We show the potential of two techniques, RESOLVE and DeepFocus, applied to ALMA calibrated science data. Significant advantages are provided with the prospect to improve the quality and completeness of the data products stored in the science archive and overall processing time for operations. Both approaches evidence the logical pathway to address the incoming revolution in data rates dictated by the planned electronic upgrades. Moreover, we bring to the community additional products through a new package, ALMASim, to promote advancements in these fields, providing a refined ALMA simulator usable by a large community for training and/or testing new algorithms.

The terrestrial detection of a neutrino burst from the next galactic core-collapse supernova (CCSN) will provide profound insight into stellar astrophysics, as well as fundamental neutrino physics. Using Time-Of-Flight (ToF) effects, a CCSN signal can be used to constrain the absolute neutrino mass. In this work, we study the case where a black hole forms during core-collapse, abruptly truncating the neutrino signal. This sharp cutoff is a feature that can be leveraged in a ToF study, enabling strict limits to be set on the neutrino mass which are largely model-independent. If supernova neutrinos are detected on Earth in liquid scintillator detectors, the exceptional energy resolution would allow an energy-dependent sampling of the ToF effects at low neutrino energies. One promising experimental program is the Jiangmen Underground Neutrino Observatory (JUNO), a next-generation liquid scintillator detector currently under construction in China. Using three-dimensional black hole-forming core-collapse supernova simulations, the sensitivity of a JUNO-like detector to the absolute neutrino mass is conservatively estimated to be $m_\nu < 0.39^{+0.06}_{-0.01}$ eV for a 95% CL bound. A future-generation liquid scintillator observatory like THEIA-100 could even achieve sub-0.2 eV sensitivity.

In the present work, we apply consistency relation tests to several cosmological models, including the flat and non-flat $\Lambda$CDM models, as well as the flat XCDM model. The analysis uses a non-parametric Gaussian Processes method to reconstruct various cosmological quantities of interest, such as the Hubble parameter $H(z)$ and its derivatives from $H(z)$ data, as well as the comoving distance and its derivatives from SNe Ia data. We construct consistency relations from these quantities which should be valid only in the context of each model and test them with the current data. We were able to find a general method of constructing such consistency relations in the context of $H(z)$ reconstruction. In the case of comoving distance reconstruction, there were not a general method of constructing such relations and this work had to write an specific consistency relation for each model. From $H(z)$ data, we have analyzed consistency relations for all the three above mentioned models, while for SNe Ia data we have analyzed consistency relations only for flat and non-flat $\Lambda$CDM models. Concerning the flat $\Lambda$CDM model, some inconsistency was found, at more than $2\sigma$ c.l., with the $H(z)$ data in the interval $1.8\lesssim z\lesssim2.4$, while the other models were all consistent at this c.l. Concerning the SNe Ia data, the flat $\Lambda$CDM model was consistent in the $0<z<2.5$ interval, at $1\sigma$ c.l., while the nonflat $\Lambda$CDM model was consistent in the same interval, at 2$\sigma$ c.l.

The key to detecting neutral hydrogen during the epoch of reionization (EoR) is to separate the cosmological signal from the dominating foreground radiation. We developed direct optimal mapping (Xu et al. 2022) to map interferometric visibilities; it contains only linear operations, with full knowledge of point spread functions from visibilities to images. Here we present an FFT-based image power spectrum and its window functions based on direct optimal mapping. We use noiseless simulation, based on the Hydrogen Epoch of Reionization Array (HERA) Phase I configuration, to study the image power spectrum properties. The window functions show $<10^{-11}$ power leakage from the foreground-dominated region into the EoR window; the 2D and 1D power spectra also verify the separation between the foregrounds and the EoR. Furthermore, we simulated visibilities from a $uv$-complete array and calculated its image power spectrum. The result shows that the foreground--EoR leakage is further suppressed below $10^{-12}$, dominated by the tapering function sidelobes; the 2D power spectrum does not show signs of the horizon wedge. The $uv$-complete result provides a reference case for future 21cm cosmology array designs.

Predicting the impacts of coronal mass ejections (CMEs) is a major focus of current space weather forecasting efforts. Typically, CME properties are reconstructed from stereoscopic coronal images and then used to forward model a CME's interplanetary evolution. Knowing the uncertainty in the coronal reconstructions is then a critical factor in determining the uncertainty of any predictions. A growing number of catalogs of coronal CME reconstructions exist, but no extensive comparison between these catalogs has yet been performed. Here we develop a Living List of Attributes Measured in Any Coronal Reconstruction (LLAMACoRe), an online collection of individual catalogs, which we intend to continually update. In this first version, we use results from 24 different catalogs with 3D reconstructions using STEREO observations between 2007--2014. We have collated the individual catalogs, determining which reconstructions correspond to the same events. LLAMACoRe contains 2954 reconstructions for 1863 CMEs. Of these, 510 CMEs contain multiple reconstructions from different catalogs. Using the best-constrained values for each CME, we find that the combined catalog reproduces the generally known solar cycle trends. We determine the typical difference we would expect between two independent reconstructions of the same event and find values of 4.0 deg in the latitude, 8.0 deg in the longitude, 24.0 deg in the tilt, 9.5 deg in the angular width, 0.1 in the shape parameter kappa, 115 km/s in the velocity, and 2.5e15 g in the mass. These remain the most probable values over the solar cycle, though we find more extreme outliers in the deviation toward solar maximum.

The angular correlation of pulsar residuals observed by NANOGrav and other pulsar timing array (PTA) collaborations show evidence in support of the Hellings-Downs correlation expected from stochastic gravitational waves (SGW). In this paper, we offer a non-gravitational wave explanation of the observed pulsar timing correlations as caused by an ultra-light $L_{\mu} - L_{\tau}$ gauge boson dark matter (ULDM). ULDM can affect the pulsar correlations in two ways. The gravitational potential of vector ULDM gives rise to a Shapiro time-delay of the pulsar signals and a non-trivial angular correlation (as compared to the scalar ULDM case). In addition, if the pulsars have a non-zero charge of the dark matter gauge group then the electric field of the local dark matter causes an oscillation of the pulsar and a corresponding Doppler shift of the pulsar signal. We point out that pulsars carry a significant charge of muons and thus the $L_{\mu} - L_{\tau}$ vector dark matter contributes to both the Doppler oscillations and the time-delay of the pulsar signals. Our analysis shows that the NANOGrav data has a better fit to the $L_{\mu} - L_{\tau}$ ULDM scenario compared to the SGW or the SGW with Shapiro time-delay hypotheses.

We systematically study a family of loop quantizations for the classical Kruskal spacetimes using the effective description motivated from loop quantum gravity for four generic parameters, $c_o, m, \delta_b$ and $\delta_c$, where the latter two denote the polymerization parameters which capture the underlying quantum geometry. We focus on the family where polymerization parameters are constant on dynamical trajectories, and of which the Ashtekar-Olmedo-Singh (AOS) and Corichi-Singh (CS) models appear as special cases. We study general features of singularity resolution in all these models due to quantum gravity effects and analytically extend the solutions across the white hole (WH) and black hole (BH) horizons to the exterior. We find that the leading term in the asymptotic expansion of the Kretschmann scalar is $r^{-4}$. However, for AOS and CS models black holes with masses greater than solar mass the dominant term behaves as $r^{-6}$ for the size of the observable universe and {our analysis can be used to phenomenologically constrain the choice of parameters for other models.} In addition, one can uniquely fix the parameter $c_o$ by requiring that the Hawking temperature at the BH horizon to the leading order be consistent with its classical value for a macroscopic BH. Assuming that the BH and WH masses are of the same order, we are able to identify a family of choices of $\delta_b$ and $\delta_c$ which share all the desired properties of the AOS model.

In this brief contribution I will highlight some directions where the developments in the frontier of (quantum) metrology may be key for fundamental high energy physics (HEP). I will focus on the detection of dark matter and gravitational waves, and introduce ideas from atomic clocks and magnetometers, large atomic interferometers and detection of small fields in electromagnetic cavities. Far from being comprehensive, this contribution is an invitation to everyone in the HEP and quantum technologies communities to explore this fascinating topic.

We show that squeezing is a crucial resource for interferometers based on the spatial separation of ultra-cold interacting matter. Atomic interactions lead to a general limitation for the precision of these atom interferometers, which can neither be surpassed by larger atom numbers nor by conventional phase or number squeezing. However, tailored squeezed states allow to overcome this sensitivity bound by anticipating the major detrimental effect that arises from the interactions. We envisage applications in future high-precision differential matter-wave interferometers, in particular gradiometers, e.g., for gravitational-wave detection.

We propose a novel method to significantly enhance the signal rate in the qubit-based dark matter detection experiments with the help of quantum interference. Various quantum sensors possess ideal properties for detecting wave-like dark matter, and qubits, commonly employed in quantum computers, are excellent candidates for dark matter detectors. We demonstrate that, by designing an appropriate quantum circuit to manipulate the qubits, the signal rate scales proportionally to $n_{\rm q}^2$, with $n_{\rm q}$ being the number of sensor qubits, rather than linearly with $n_{\rm q}$. Consequently, in the dark matter detection with a substantial number of sensor qubits, a significant increase in the signal rate can be expected. We provide a specific example of a quantum circuit that achieves this enhancement by coherently combining the phase evolution in each individual qubit due to its interaction with dark matter. We also demonstrate that the circuit is fault tolerant to de-phasing noises, a critical quantum noise source in quantum computers. The enhancement mechanism proposed here is applicable to various modalities for quantum computers, provided that the quantum operations relevant to enhancing the dark matter signal can be applied to these devices.

Superradiance is the effect of field waves being amplified during reflection from a charged or rotating black hole. In this paper, we study the low-energy dynamics of super-radiant scattering of massive scalar and massless higher spin field perturbations in a generic axisymmetric stationary Kerr-Taub-NUT (Newman-Unti-Tamburino) spacetime, which represents sources with both gravitomagnetic monopole moment (magnetic mass) and gravitomagnetic dipole moment (angular momentum). We obtain a generalized Teukolsky master equation for all spin perturbation fields. The equations are separated into their angular and radial parts. The angular equations lead to spin-weighted spheroidal harmonic functions that generalize those in Kerr spacetime. We identify an effective spin as a coupling between frequency (or energy) and the NUT parameter. The behaviors of the radial wave function near the horizon and at the infinite boundary are studied. We provide analytical expressions for low-energy observables such as emission rates and cross sections of all massless fields with spin, including scalar, neutrino, electromagnetic, Rarita-Schwinger, and gravitational waves.

We explore the potential for detecting rotational instabilities in the post-merger phase of binary neutron star mergers using different network configurations of upgraded and next-generation gravitational wave detectors. Our study employs numerically generated post-merger waveforms, which reveal the re-excitation of the $l=m=2$ $f$-mode at a time of $O(10{\rm})$ms after merger. We evaluate the detectability of these signals by injecting them into colored Gaussian noise and performing a reconstruction as a sum of wavelets using Bayesian inference. Computing the overlap between the reconstructed and injected signal, restricted to the instability part of the post-merger phase, we find that one could infer the presence of rotational instabilities with a network of planned 3rd-generation detectors, depending on the total mass and distance to the source. For a recently suggested high-frequency detector design, we find that the instability part would be detectable even at 200 Mpc, significantly increasing the anticipated detection rate. For a network consisting of the existing HLV detectors, but upgraded to twice the A+ sensitivity, we confirm that the peak frequency of the whole post-merger gravitational-wave emission could be detectable with a network signal-to-noise ratio of 8 at a distance of 40Mpc.