Papers for Wednesday, Mar 20 2024

We present an analysis of high angular resolution images of the microlensing target MOA-2007-BLG-192 using Keck adaptive optics and the Hubble Space Telescope. The planetary host star is robustly detected as it separates from the background source star in nearly all of the Keck and Hubble data. The amplitude and direction of the lens-source separation allows us to break a degeneracy related to the microlensing parallax and source radius crossing time. Thus, we are able to reduce the number of possible solutions by a factor of ${\sim}2$, demonstrating the power of high angular resolution follow-up imaging for events with sparse light curve coverage. Following Bennett et al. 2023, we apply constraints from the high resolution imaging on the light curve modeling to find host star and planet masses of $M_{\textrm{host}} = 0.28 \pm 0.04M_{\odot}$ and $m_p = 12.49^{+65.47}_{-8.03}M_{\oplus}$ at a distance from Earth of $D_L = 2.16 \pm 0.30\,$kpc. This work illustrates the necessity for the Nancy Grace Roman Galactic Exoplanet Survey (RGES) to use its own high resolution imaging to inform light curve modeling for microlensing planets that the mission discovers.

The rise of synoptic sky surveys has ushered in an era of big data in time-domain astronomy, making data science and machine learning essential tools for studying celestial objects. Tree-based (e.g. Random Forests) and deep learning models represent the current standard in the field. We explore the use of different distance metrics to aid in the classification of objects. For this, we developed a new distance metric based classifier called DistClassiPy. The direct use of distance metrics is an approach that has not been explored in time-domain astronomy, but distance-based methods can aid in increasing the interpretability of the classification result and decrease the computational costs. In particular, we classify light curves of variable stars by comparing the distances between objects of different classes. Using 18 distance metrics applied to a catalog of 6,000 variable stars in 10 classes, we demonstrate classification and dimensionality reduction. We show that this classifier meets state-of-the-art performance but has lower computational requirements and improved interpretability. We have made DistClassiPy open-source and accessible at https://pypi.org/project/distclassipy/ with the goal of broadening its applications to other classification scenarios within and beyond astronomy.

We perform simulations here that include the gravitational effects of the primordial planetesimal belt consisting of ~10^5 massive bodies. In our simulations, Neptune unlocks from resonance with the other giant planets and begins to migrate outward due to interactions with planetesimals before a planetary orbital instability is triggered, and afterward, residual Neptunian migration completes the formation of the modern Kuiper belt. Our present work exhibits a number of notable differences from prior work. First, Neptune's planetary resonance unlocking requires the Neptunian 3:2 mean motion resonance to sweep much of the primordial disk interior to 30 au prior to the giant planet instability. The pre-instability population of planetesimals is consequently lower in semimajor axis, eccentricity, and inclination, and this effect persists after the instability. Second, direct scattering between Pluto-mass bodies and other small bodies removes material from Neptunian resonances more efficiently than resonant dropout resulting from small changes in Neptune's semimajor axis during scattering between Pluto-mass bodies and Neptune. Thus, the primordial population of Pluto-mass bodies may be as few as ~200 objects. Finally, our simulation end states display a wide variety of orbital distributions, and clear relationships between final bulk Kuiper belt properties and Neptune's migration or initial planetesimal properties largely elude us. In particular, we find that the rapid, stochastic planetary orbital evolution occurring during the giant planet instability can significantly alter final Kuiper belt properties such as its inclination dispersion and the prominence of resonant populations. This complicates using modern Kuiper belt properties to confidently constrain early solar system events and conditions, including planetary orbital migration and the primordial Kuiper belt's characteristics.

Eruptive stars are a class of young stellar objects that show an abrupt increase in their luminosity. These burst-like episodes are thought to dominate the stellar accretion process during the class 0/class I stage. We present an overview of a survey of seven episodically accreting protostars aimed at studying their potentially complex circumstellar surroundings. The observations were performed with the instrument SPHERE, mounted at the VLT. We observed the eruptive stars in $H$-band with the near-infrared imager IRDIS and used the polarimeter to extract the polarized light scattered from the stars' surroundings. We produced polarized light images for three FUor objects, Z CMa, V960 Mon, and FU Ori, and four EXor objects, XZ Tau, UZ Tau, NY Ori, and EX Lup. We calculated the intrinsic polarization fraction for all the observed stars. In all systems we registered scattered light from around the primary star. FU Ori and V960 Mon are surrounded by complex structures including spiral-like features. In Z CMa, we detected a point source 0.7 arcsec to the northeast of the primary. Based on the astrometric measurements from archival Keck/NIRC2 data, we find this source to be a third member of the system. Further, Z CMa displays an outflow extending thousands of au. Unlike the other EXor objects in our sample, XZ Tau shows bright, extended scattered light structures, also associated with an outflow on a scale of hundreds of au. The other EXors show relatively faint disk-like structures in the immediate vicinity of the coronagraph. Asymmetric arms were only found around FUor objects, while faint disks seem to predominantly occur around EXors. Importantly, for Z CMa the detection of the faint extended structure questions previous interpretations of the system's dynamic state. The streamer which was associated with a fly-by object turned out to be part of a huge outflow extending 6000 au.

We present a machine-learning-based model for the total density profiles of subhaloes with masses $M \gtrsim 7\times 10^8\,h^{-1}{\rm M}_\odot$ in the IllustrisTNG100 simulation. The model is based on an interpretable variational encoder (IVE) which returns the independent factors of variation in the density profiles within a low-dimensional representation, as well as the predictions for the density profiles themselves. The IVE returns accurate and unbiased predictions on all radial ranges, including the outer region profile where the subhaloes experience tidal stripping; here its fit accuracy exceeds that of the commonly used Einasto profile. The IVE discovers three independent degrees of freedom in the profiles, which can be interpreted in terms of the formation history of the subhaloes. In addition to the two parameters controlling the normalization and inner shape of the profile, the IVE discovers a third parameter that accounts for the impact of tidal stripping onto the subhalo outer profile; this parameter is sensitive to the mass loss experienced by the subhalo after its infall onto its parent halo. Baryonic physics in the IllustrisTNG galaxy formation model does not impact the number of degrees of freedom identified in the profile compared to the pure dark matter expectations, nor their physical interpretation. Our newly proposed profile fit can be used in strong lensing analyses or other observational studies which aim to constrain cosmology from small-scale structures.

Gyrochronology, a valuable tool for determining ages of low-mass stars where other techniques fail, relies on accurate calibration. We present a sample of 327 wide ($>$$100$\,au) white dwarf + main sequence (WD + MS) binary systems. Total ages of WDs are computed using all-sky survey photometry, Gaia parallaxes, and current hydrogen atmosphere WD models. Using a magnetic braking law calibrated against open clusters, along with assumptions about initial conditions and angular momentum transport, we construct gyrochrones to predict the rotation periods of the MS stars. Both data and models show that, near the fully convective boundary, MS stars with WD ages up to 7.5\,Gyr experience a rotation period increase by up to a factor of $\approx$$3$ within a $<50\,\mathrm{K}$ effective temperature range. We suggest that rapid braking at this boundary is driven by a sharp rise in the convective overturn timescale ($\tau_{\mathrm{cz}}$) caused by structural changes between partially and fully convective stars and the $^3 \textrm{He}$ instability occurring at this boundary. While the specific location in mass (or temperature) of this feature varies with model physics, we argue that its existence remains consistent. Stars along this feature exhibit rotation periods that can be mapped, within 1$\sigma$, to a range of gyrochrones spanning $\approx 6$\, Gyr. Due to current temperature errors ($\simeq$$50\,\mathrm{K}$), this implies that a measured rotation period cannot be uniquely associated to a single gyrochrone, implying that gyrochronology may not be feasible for M dwarfs very close to the fully convective boundary.

We report the discovery of eight optical counterparts of ALFALFA extragalactic objects from DECaLS, five of which are discovered for the first time. These objects were flagged as HI emission sources with no optical counterparts in SDSS before. Multi-band data reveal their unusual physical properties. They are faint and blue ($g-r=-0.35\sim0.55$), with quite low surface brightness ($\mu_{\rm g,peak}=24.88\sim26.41\,{\rm mag}/{\rm arcsec}^2$), irregular morphologies, low stellar masses ($log_{10}(M_{*}/M_\odot)=5.27\sim7.15$), low star formation rates ($SFR=0.21\sim9.24\times10^{-3}\,{M_\odot}\,{\rm yr}^{-1}$), and remarkably high HI-to-stellar mass ratios ($log_{10}(M_{\rm HI}/M_{*}) = 1.72\sim3.22$, except AGC\,215415). They deviate from the scaling relations between HI and optical properties defined by the ALFALFA sample and the baryonic Tully-Fisher relation. They agree well with the main sequence of star-forming galaxies but exhibit low star-forming efficiency. Based on their physical properties and environments, we speculate that six of these objects may have originated from tidal processes, while the remaining two appear to have isolated origins. They may have had a relatively calm evolutionary history and only begun to form stars recently.

We study the reionization of the Universe due to haloes that host galaxies undergoing bursts of star formation. By comparing the recent results from the James Webb Space Telescope (JWST) with the cosmological hydrodynamical simulation EAGLE at $z\ge 6$, we find that bursty galaxies have specific star formation rate, sSFR $>10^{-2}$ Myr$^{-1}$, and magnitude, $M_{\rm UV}\leq -17$. Most of them reside in haloes of mass $\sim 10^9$ M$_\odot$ and some in more massive haloes. We then construct the models of escape fraction and find that a skewed Gaussian function with a flat tail towards the high mass end best describes the mean dependence of escape fraction on halo mass, considering the haloes hosting bursty galaxies as the primary drivers of reionization. We implement the models of escape fraction in the code 21cmFAST to study the progress of reionization and derive the evolution of the mean ionized fraction that agrees well with observations. We also calculate the brightness temperature, spin temperature, and kinetic temperature and further study the spatial fluctuations in these quantities to gain insights into the progress of reionization. We compute the 21 cm power spectrum and predict a peak in power at $180$ MHz corresponding to redshift, $z\approx6.8$, that is testable by the upcoming Square Kilometre Array (SKA). Our findings suggest that the Universe was reionized by the haloes of $\sim 10^{9}$ M$_\odot$, with haloes more massive than that also provided an essential contribution.

The gravitational capture of a stellar-mass object by a supermassive black hole represents a unique probe of warped spacetime. The small object, typically a stellar-mass black hole, describes a very large number of cycles before crossing the event horizon. Because of the mass difference, we call these captures extreme-mass ratio inspirals (EMRIs). Their merger event rate at the Galactic Centre is negligible, but the amount of time spent in the early inspiral is not. Early EMRIs (E-EMRIs) spend hundreds of thousands of years in band during this phase. At very early stages, the peak of the frequency will not change during an observational time. At later stages, it will change a bit and finally the EMRI explores a wide range of them when it is close to merger. We distinguish between ``monocromatic'' E-EMRIs, which do not change their (peak) frequency, oligochromatic E-EMRIs, which explore a short range and polychromatic ones, the EMRIs which have been discussed so far in the literature. We derive the number of E-EMRIs at the Galactic Centre, and their signal-to-noise ratios (SNR) and perform a study of parameter extraction. We show that parameters such as the spin and the mass can be extracted with an error which can be as small as $10^{-11}$ and $10^{-5}\,M_{\odot}$. There are between hundreds and thousands of E-EMRIs in their monochromatic stage at the GC, and tens in their oligochromatic phase. The SNR ranges from a minimum of $10$ (larger likelihood) to a maximum of $10^6$ (smaller likelihood). Moreover, we derive the contribution signal corresponding to the incoherent sum of continuous with two representatives masses; $10\,M_{\odot}$ and $40\,M_{\odot}$ and show that their curves will cover a significant part of LISA's sensitivity curve. Depending on their level of circularisation, they might be detected as individual sources or form a foreground population.

Stellar feedback plays a crucial role in regulating baryon cycles of a galactic ecosystem, and may manifest itself in the formation of superbubbles in the interstellar medium. In this work, we used a set of high-resolution simulations to systematically study the properties and evolution of superbubbles in galactic environments. The simulations were based on the SMUGGLE galaxy formation framework using the hydrodynamical moving-mesh code Arepo, reaching a spatial resolution of $\sim 4 \, \rm pc$ and mass resolution of $\sim 10^3 \, \rm M_{\odot}$. We identified superbubbles and tracked their time evolution using the parent stellar associations within the bubbles. The X-ray luminosity-size distribution of superbubbles in the fiducial run is largely consistent with the observations of nearby galaxies. The size of superbubbles shows a double-peaked distribution, with the peaks attributed to early feedback (radiative and stellar wind feedback) and supernova feedback. The early feedback tends to suppress the subsequent supernova feedback, and it is strongly influenced by star formation efficiency, which regulates the environmental density. Our results show that the volume filling factor of hot gas ($T > 10^{5.5} ~\mathrm{K}$) is about $12 \%$ averaged over a region of 4 kpc in height and 20 kpc in radius centered on the disk of the galaxy. Overall, the properties of superbubbles are sensitive to the choice of subgrid galaxy formation models and can, therefore, be used to constrain these models.

We investigate the coupling between the temporal variation from galaxy-formation feedback and the bar instability. We show that fluctuations from mass outflow on star-formation time scales affect the radial motion of disk orbits. The resulting incoherence in orbital phase leads to the disruption of the bar-forming dynamics. Bar formation is suppressed in starburst galaxies that have fluctuation time scales within the range 10 Myr to 200 Myr with repeated events with wind mass 15% of the disk within 0.5 scale lengths or 1.4% of the total disk mass. The work done by feedback is capable of reducing the amplitude or, with enough amplitude, destroying an existing bar. AGN feedback with similar amplitude and timescales would have similar behavior. To model the dynamics of the coupling and interpret the results of the full N-body simulations, we introduce a generalization of the Hamiltonian mean-field (HMF) model, drawing inspiration from the Lynden-Bell (1979) mechanism for bar growth. Our non-linear 'BarHMF' model is designed to reproduce linear perturbation theory in the low-amplitude limit. Notably, without star-formation feedback, this model exhibits exponential growth whose rate depends on disk mass and reproduces the expected saturation of bar growth observed in N-body simulations. We describe several promising applications of the BarHMF model beyond this study.

Recent observations from the EIGER JWST program have measured for the first time the quasar-galaxy cross-correlation function at $z\approx6$. The auto-correlation function of faint $z\approx6$ quasars was also recently estimated. These measurements provide key insights into the properties of quasars and galaxies at high redshift and their relation with the host dark matter halos. In this work, we interpret these data building upon an empirical quasar population model that has been applied successfully to quasar clustering and demographic measurements at $z\approx2-4$. We make use of a new, large-volume N-body simulation with more than a trillion particles, FLAMINGO-10k, to model quasars and galaxies simultaneously. We successfully reproduce observations of $z\approx6$ quasars and galaxies (i.e., their clustering properties and luminosity functions), and infer key quantities such as their luminosity-halo mass relation, the mass function of their host halos, and their duty cycle/occupation fraction. Our key findings are: (i) quasars reside on average in $\approx10^{12.4}\,{\rm M}_\odot$ halos (corresponding to $\approx5\sigma$ fluctuations in the initial conditions of the linear density field), but the distribution of host halo masses is quite broad; (ii) the duty cycle of (UV-bright) quasar activity is relatively low ($\lesssim1\%$); (iii) galaxies (that are bright in [OIII]) live in much smaller halos ($\approx10^{10.9}\,{\rm M}_\odot$) and have a larger duty cycle (occupation fraction) of $\approx20\%$. Finally, we focus on the inferred properties of quasars and present a homogeneous analysis of their evolution with redshift. The picture that emerges reveals a strong evolution of the host halo mass and duty cycle of quasars at $z\approx2-6$, and calls for new investigations of the role of quasar activity across cosmic time.

We report the discovery of a $(1.0 \pm 0.28) \times 10^{10}$ M$_\odot$ Supermassive Black Hole (BH) at the centre of NGC 708, the Brightest Cluster Galaxy of Abell 262. Such high BH masses are very rare and allow to investigate BH - host galaxy scaling relations at the high mass end, which in turn provide hints about the (co)evolution of such systems. NGC~708 is found to be an outlier in all the canonical scaling relations except for those linking the BH mass to the core properties. The galaxy mass-to-light ratio points to a Kroupa IMF rather than Salpeter, with this finding confirmed using photometry in two different bands. We perform this analysis using our novel triaxial Schwarzschild code to integrate orbits in a 5-dimensional space, using a semi-parametric deprojected light density to build the potential and non-parametric Line-of-Sight Velocity Distributions (LOSVDs) derived from long-slit spectra recently acquired at Large Binocular Telescope (LBT) to exploit the full information in the kinematic. We find that the galaxy geometry changes as a function of the radius going from prolate, nearly spherical in the central regions to triaxial at large radii, highlighting the need to go beyond constant shape profiles. Our analysis is only the second of its kind and will systematically be used in the future to hunt Supermassive Black Holes in giant ellipticals.

Stellar multiplicity is among the oldest and richest problems in astrophysics. Binary stars are a cornerstone of stellar mass and radius measurements that underpin modern stellar evolutionary models. Binaries are the progenitors of many of the most interesting and exotic astrophysical phenomena, ranging from type Ia supernovae to gamma ray bursts, hypervelocity stars, and most detectable stellar black holes. They are also ubiquitous, accounting for about half of all stars in the Universe. In the era of gravitational waves, wide-field surveys, and open-source stellar models, binaries are coming back stronger than a nineties trend. Much of the progress in the last decade has been enabled by the Gaia mission, which provides high-precision astrometry for more than a billion stars in the Milky Way. The Gaia data probe a wider range of binary separations and mass ratios than most previous surveys, enabling both an improved binary population census and discovery of rare objects. I summarize recent results in the study of binary stars brought about by Gaia, focusing in particular on developments related to wide ($a \gtrsim 100$ au) binaries, evidence of binarity from astrometric noise and proper motion anomaly, astrometric and radial velocity orbits from Gaia DR3, and binaries containing non-accreting compact objects. Limitations of the Gaia data, the importance of ground-based follow-up, and anticipated improvements with Gaia DR4 are also discussed.

This study aims to analyze Phobos' photometric properties using Mars Express mission observations to support the Martian Moons eXploration mission. We analyzed resolved images of Phobos acquired between 2004 and 2022 by the HRSC and the SRC cameras on board the Mars Express spacecraft. We performed photometric analysis using the Hapke model for both integrated and disk-resolved data. The Phobos phase function has a strong opposition effect due to shadow hiding, with an amplitude and a half-width of the opposition surge of 2.28$\pm$0.03 and 0.0573$\pm$0.0001, respectively. Overall, the surface of Phobos is dark, with a geometric albedo of 6.8 % in the green filter and backscattering. We also found a surface porosity of 87\%, indicating the presence of a thick dust mantle or of fractal aggregates on the top surface. The SSA maps revealed high reflectance variability, with the blue unit area in the northeast Stickney rim being up to 65\% brighter than average, while the Stickney floor is among the darkest regions, with reflectance 10 to 20% lower than average. Photometric modeling of the regions of interest selected in the red and blue units indicates that red unit terrains have a stronger opposition effect and a smaller SSA value than the blue ones, but they have similar porosity and backscattering properties. The HRSC data provide a unique investigation of the Phobos phase function and opposition surge, which is valuable information for the MMX observational planning. The Phobos opposition surge, surface porosity, phase integral, and spectral slope are very similar to the values observed for the comet 67P and for Jupiter family comets in general. Based on these similarities, we formulate a hypothesis that the Mars satellites might be the results of a binary or bilobated comet captured by Mars.

Recent surveys have demonstrated the widespread presence of UV emission in early-type (elliptical/S0) galaxies, suggesting the presence of star formation in many of these systems. However, potential UV contributions from old and young stars, together with model uncertainties, makes it challenging to confirm the presence of young stars using integrated photometry alone. This is particularly true in ETGs that are fainter in the UV and have red UV-optical colours. An unambiguous way of disentangling the source of the UV is to look for structure in UV images. Optical images of ETGs, which are dominated by old stars, are smooth and devoid of structure. If the UV is also produced by these old stars, then the UV images will share this smoothness, while, if driven by young stars, they will exhibit significant structure. We compare the UV and optical morphologies of 32 ETGs (93 percent of which are at $z<0.03$) using quantitative parameters (concentration, asymmetry, clumpiness and the S\'ersic index), calculated via deep UV and optical images with similar resolution. Regardless of stellar mass, UV-optical colour or the presence of interactions, the asymmetry and clumpiness of ETGs is significantly (often several orders of magnitudes) larger in the UV than in the optical, while the UV S\'ersic indices are typically lower than their optical counterparts. The ubiquitous presence of structure demonstrates that the UV flux across our entire ETG sample is dominated by young stars and indicates that star formation exists in all ETGs in the nearby Universe.

Recently, nearby active galactic nuclei (AGN) have been subject to long X-rays/UV/optical monitoring campaigns. These campaigns reveal a strong correlation between the various UV and optical bands, with time lags increasing with wavelength. In a series of papers, we demonstrated that a scenario in which a central X-ray source illuminates the accretion disc explains the observed correlations. However, some of the monitored AGN show low/moderate X-rays-UV correlations, which could challenge this scenario. In this paper, we study the broadband X-ray/UV/optical spectral energy distributions (SEDs) of NGC 5548, one of the most intensively monitored AGN. We aim to test if the X-ray illumination model explains the broadband spectral behaviour of the source, despite the moderate X-ray-UV/optical correlation. We model the broadband time-averaged SED, from the STORM monitoring campaign of the source, using the KYNSED model which assumes an X-ray illuminated disc. We assume that the accretion process powers the X-ray corona. We also model 15 time-resolved SEDs from the same campaign to check whether this scenario can account for the observed spectral variability. The proposed model describes well the time-averaged and the time-resolved SEDs of NGC 5548. In this scenario, the corona height, the X-ray photon index, and the power transferred to the corona all vary. This explains the variability behaviour at different wavelengths. The best-fit model is obtained for a non-rotating black hole accreting at a constant rate of 5% its Eddington limit. Since each of the variable parameters affects the observed flux in a particular way, the combined variability of all parameters explains the moderate X-ray-UV/optical correlation. The X-ray illuminated disc model provides a complete description of the behaviour of NGC 5548, explaining its broadband SEDs, time-lag spectrum, and its power spectral distribution.

Aims. We use near-infrared, ground-based data from the VISTA Variables in the Via Lactea (VVV) survey to indirectly extend the astrometry provided by the Gaia catalog to objects in heavily-extincted regions towards the Galactic bulge and plane that are beyond Gaia's reach. Methods. We make use of the state-of-the-art techniques developed for high-precision astrometry and photometry with the Hubble Space Telescope to process the VVV data. We employ empirical, spatially-variable, effective point-spread functions and local transformations to mitigate the effects of systematic errors, like residual geometric distortion and image motion, and to improve measurements in crowded fields and for faint stars. We also anchor our astrometry to the absolute reference frame of the Gaia Data Release 3. Results. We measure between 20 and 60 times more sources than Gaia in the region surrounding the Galactic center, obtaining an single-exposure precision of about 12 mas and a proper-motion precision of better than 1 mas yr$^{-1}$ for bright, unsaturated sources. Our astrometry provides an extension of Gaia into the Galactic center. We publicly release the astro-photometric catalogs of the two VVV fields considered in this work, which contain a total of $\sim$ 3.5 million sources. Our catalogs cover $\sim$ 3 sq. degrees, about 0.5% of the entire VVV survey area.

It is well-known that Einstein's equations constrain only the total energy-momentum tensor of the cosmic substratum, without specifying the characteristics of its individual constituents. Consequently, cosmological models featuring distinct decompositions within the dark sector, while sharing identical values for the sum of dark components' energy-momentum tensor, remain indistinguishable when assessed through observables based on distance measurements. Notably, it has been already demonstrated that cosmological models with dynamical descriptions of dark energy, characterized by a time-dependent equation of state (EoS), can always be mapped into a model featuring a decaying vacuum ($w=-1$) coupled with dark matter. We explore the possibility of breaking this degeneracy by using measurements of the gas mass fraction observed in massive and relaxed galaxy clusters. This data is particularly interesting for this purpose because it isolates the matter contribution, possibly allowing the degeneracy breaking. We study the particular case of the $w$CDM model with its interactive counterpart. We compare the results obtained from both descriptions with a non-parametric analysis obtained through Gaussian Process. Even though the degeneracy may be broken from the theoretical point of view, we find that current gas mass fraction data seems to be insufficient for a final conclusion about which approach is favored, even when combined with SNIa, BAO and CMB.

It is well established that stellar discs are destabilized by sharp features in their phase space, driving recurrent spiral modes. We explore the extent to which surface-density breaks in disc galaxies - which represent sharp changes in the gradient of the disc density - drive new spiral modes. We employ linear perturbation theory to investigate how disc breaks alter the eigenmode spectrum of an otherwise pure exponential disc. We find that the presence of a density break gives rise to a set of new, vigorously growing, modes. For a given multiplicity, these edge modes occur in pairs, with closely separated resonances between each pair. The growth rate of edge modes decreases when the break is weakened or moved outward to lower-density regions of the disc. Both down- and up-bending profiles excite edge modes, whose origin can be best understood via the gravitational torques they exert on the underlying disc. When the profile is down-bending (Type II) the faster growing mode is the inner one while in the up-bending (Type III) case the outer mode is faster growing. In both cases the faster growing mode has a corotation almost coincident with the break. We show that the torques of the edge modes tend to smoothen the break.

Until the advent of the SOLARNET recommendations, metadata sharing of simulated data within the Solar Physics community has been mostly on a "private communication" basis, with the description of the data format and content conveyed in an ad hoc manner. This document aims to amend this situation by establishing recommendations for representing such data and the associated metadata, based on the SOLARNET Metadata Recommendations for Solar Observations (arXiv:2011.12139)

Supernova remnants (SNRs) exhibit varying degrees of anisotropy, which have been extensively modeled using numerical methods. We implement a technique to measure anisotropies in SNRs by calculating power spectra from their high-resolution images. To test this technique, we develop 3D hydrodynamical models of supernova remnants and generate synthetic x-ray images from them. Power spectra extracted from both the 3D models and the synthetic images exhibit the same dominant angular scale, which separates large scale features from small scale features due to hydrodynamic instabilities. The angular power spectrum at small length scales during relatively early times is too steep to be consistent with Kolmogorov turbulence, but it transitions to Kolmogorov turbulence at late times. As an example of how this technique can be applied to observations, we extract a power spectrum from a \textit{Chandra} observation of Tycho's SNR and compare with our models. Our predicted power spectrum picks out the angular scale of Tycho's fleece-like structures and also agrees with the small-scale power seen in Tycho. We use this to extract an estimate for the density of the circumstellar gas ($n \sim 0.28/\mathrm{cm^3}$), consistent with previous measurements of this density by other means. The power spectrum also provides an estimate of the density profile of the outermost ejecta. Moreover, we observe additional power at large scales which may provide important clues about the explosion mechanism itself.

We present a sample of 305 QSO candidates having $|b| < 30^{\circ}$, the majority with GALEX magnitudes NUV < 18.75. To generate this sample, we apply UV-IR color selection criteria to photometric data from the Ultraviolet GAlactic Plane Survey (UVGAPS) as part of GALEX-CAUSE, the Million Quasars Catalog, Gaia DR2, and Pan-STARRS DR1. 165 of these 305 candidate UV-bright AGN (54%) have published spectroscopic redshifts from 45 different surveys, confirming them as AGN. We further obtained low-dispersion, optical, longslit spectra with the APO 3.5-m, MDM 2.4-m, and MDM 1.3-m telescopes for 84 of the candidates, and confirm 86% (N = 72) as AGN, generally with z < 0.6. These sources fill a gap in the Galactic latitude coverage of the available samples of known UV-bright QSO background probes. Along with a description of the confirmed QSO properties, we provide the fully-reduced, flux and wavelength-calibrated spectra of 84 low-latitude QSOs through the Mikulski Archive for Space Telescopes. Future HST/COS spectroscopy of these low-Galactic-latitude QSOs has the potential to transform our view of the Milky Way and Local Group circumgalactic medium.

The combined strength of asteroseismology and empirical stellar basic parameter determinations for in-depth asteroseismic analysis of massive pulsators in eclipsing binaries shows great potential for treating the challenging and mysterious discrepancies between observations and models of stellar structure and evolution of massive stars. This paper compiles a comprehensive list of massive pulsators in eclipsing binary systems observed with TESS. The TESS light curves and Discrete Fourier Transforms (DFT) of a sample of 8055 stars of spectral type B0--B3 were examined for eclipses and stellar pulsations and the ephemerides of the resulting sub-sample of massive pulsators in eclipsing binaries were computed. This sub-sample was also cross-matched with existing catalogues of massive pulsators. Until now, fewer than 30 $\beta$ Cep pulsators in eclipsing binaries have been reported in the literature. Here we announce a total of 78 pulsators of the $\beta$ Cephei type in eclipsing binaries, 59 of which are new discoveries. Forty-three are recognized as definite and 35 are candidate pulsators. Our sample of pulsating massive stars in eclipsing binaries allows for future asteroseismic modelling to better understand the internal mixing profile and to resolve the mass discrepancy in massive stars. We have already started follow-up of some of the most interesting candidates.

We present an X-ray analysis of three different XMM-Newton observations together with simultaneous NICER and NuSTAR observations of the ultraluminous X-ray source NGC 4190 ULX-1. Our goal is to constrain the structure of the accretion disk and the geometrical properties of the source. We performed a temporal and spectral analyses in the 0.4--30 keV energy range where the source is significantly detected in dedicated XMM-Newton, NICER and NuSTAR observations. The temporal analysis shows no flaring activity in the light curves. No pulsation is detected throughout. The source exhibits a typical ULX spectrum, which can be fitted with two thermal blackbody components plus a Comptonization tail at high energies. The luminosity-temperature relation of each thermal spectral component is consistent with the $L \propto T^{2}$ relation expected from an advection-dominated supercritical disk. We interpret these results as a super-Eddington accreting black hole seen almost face-on. A dense wind ejected from the disk obscures the central source, and a hot electron plasma is evacuated through the funnel formed above the hole. Geometric beaming is responsible for the ULX soft emission, whereas the hard tail is the result of Comptonization of soft photons by the electrons ejected through the funnel.

The next generation of cosmological surveys is expected to generate unprecedented high-quality data, consequently increasing the already substantial computational costs of Bayesian statistical methods. This will pose a significant challenge to analyzing theoretical models of cosmology. Additionally, new mitigation techniques of baryonic effects, intrinsic alignment, and other systematic effects will inevitably introduce more parameters, slowing down the convergence of Bayesian analyses. In this scenario, machine-learning-based accelerators are a promising solution, capable of reducing the computational costs and execution time of such tools by order of thousands. Yet, they have not been able to provide accurate predictions over the wide prior ranges in parameter space adopted by Stage III/IV collaborations in studies employing real-space two-point correlation functions. This paper offers a leap in this direction by carefully investigating the modern transformer-based neural network (NN) architectures in realistic simulated Rubin Observatory year one cosmic shear $\Lambda$CDM inferences. Building on the framework introduced in Part I, we generalize the transformer block and incorporate additional layer types to develop a more versatile architecture. We present a scalable method to efficiently generate an extensive training dataset that significantly exceeds the scope of prior volumes considered in Part I, while still meeting strict accuracy standards. Through our meticulous architecture comparison and comprehensive hyperparameter optimization, we establish that the attention-based architecture performs an order of magnitude better in accuracy than widely adopted NN designs. Finally, we test and apply our emulators to calibrate tension metrics.

Investigations into the propagation characteristics, specifically loss and wave velocity, of superconducting coplanar waveguides and microstrip lines were conducted at a 2 mm wavelength. This was achieved through the measurement of on-chip half-wavelength resonators, employing superconductor-insulator-superconductor tunnel junctions as detectors. A continuous wave millimeter wave probe signal was introduced to the chip via a silicon membrane-based orthomode transducer. This setup not only facilitated the injection of the probe signal but also provided a reference path essential for differential measurements. The observed resonance frequencies aligned closely with theoretical predictions, exhibiting a discrepancy of only several percent. However, the measured losses significantly exceeded those anticipated from quasi-particle loss mechanisms, suggesting the presence of additional loss factors. Notably, the measurement results revealed that the tangential loss attributable to the dielectric layer, specifically silicon dioxide, was approximately $\rm{7\pm 2 \times 10^{-3}}$. This factor emerged as the dominant contributor to overall loss at temperatures around 4 K.

We detected a giant X-ray flare from the RS-CVn type binary star UX Ari using MAXI on 2020 August 17 and started a series of NICER observations 89 minutes later. For a week, the entire duration of the flare was covered with 32 snapshot observations including the rising phase. The X-ray luminosity reached 2$\times$10$^{33}$ erg s$^{-1}$ and the entire energy release was $\sim 10^{38}$ erg in the 0.5--8.0~keV band. X-ray spectra characterized by continuum emission with lines of Fe XXV He$\alpha$ and Fe XXVI Ly$\alpha$ were obtained. We found that the temperature peaks before that of the flux, which suggests that the period of plasma formation in the magnetic flare loop was captured. Using the continuum information (temperature, flux, and their delay time), we estimated the flare loop size to be $\sim 3 \times 10^{11}$ cm and the peak electron density to be $\sim 4\times10^{10}$ cm$^{-3}$. Furthermore, using the line ratio of Fe XXV and Fe XXVI, we investigated any potential indications of deviation from collisional ionization equilibrium (CIE). The X-ray spectra were consistent with CIE plasma throughout the flare, but the possibility of an ionizing plasma away from CIE was not rejected in the flux rising phase.

An expansion of cross-sectional area directly impacts the mass flow along a coronal loop, and significantly alters the radiative and hydrodynamic evolution of that loop as a result. Previous studies have found that an area expansion from chromosphere to corona significantly lengthens the cooling time of the corona, and appears to suppress draining from the corona. In this work, we examine the fluid dynamics to understand how the mass flow rate, the energy balance, and the cooling and draining timescales are affected by a non-uniform area. We find that in loops with moderate or large expansion (cross-sectional area expansion factors of 2, 3, 10, 30, 100 from photosphere to apex), impulsive heating, for either direct thermal heating or electron beam heating, induces a steady flow into the corona, so that the coronal density continues to rise during the cooling phase, whereas a uniform loop drains during the cooling phase. The induced upflow carries energy into the corona, balancing the losses from thermal conduction, and continues until thermal conduction weakens enough so that it can no longer support the radiative losses of the transition region (TR). As a result, the plasma cools primarily radiatively until the onset of catastrophic collapse. The speed and duration of the induced upflow both increase in proportion to the rate of area expansion. We argue that observations of blue-shifted spectral lines, therefore, could place a constraint on a loop's area expansion.

Binary Cepheids with giant companions are crucial for studying the physical properties of Cepheid variables, providing the best means to measure their masses. Systems composed of two Cepheids are even more important but to date, only one such system in the Large Magellanic Cloud (LMC) was known. Our current aim is to increase the number of these systems tenfold and provide their basic characteristics. The final goal is to obtain the physical properties of the component Cepheids, including their masses and radii, and to learn about their evolution in the multiple systems, also revealing their origin. We started a spectroscopic monitoring of nine unresolved pairs of Cepheids from the OGLE catalog, to check if they are gravitationally bound. Two of these so-called double Cepheids are located in the LMC, five in the Small Magellanic Cloud (SMC), and two in the Milky Way (MW). We report the spectroscopic detection of binarity of all 9 of these double Cepheids with orbital periods from 2 to 18 years. This increases the number of known binary double (BIND) Cepheids from 1 to 10 and triples the number of all confirmed double-lined binary (SB2) Cepheids. For five BIND Cepheids disentangled pulsational light curves of the components show anti-correlated phase shifts due to orbital motion. We show the first empirical evidence that typical period-luminosity relations (PLRs) are rather binary Cepheid PLRs that include the companion's light. The statistics of pulsation period ratios of BIND Cepheids do not agree with those expected for pairs of the same-age Cepheids. These ratios together with the mass ratios far from unity suggest merger-origin of at least one component for about half of the systems. The SMC and MW objects are the first found in SB2 systems composed of giants in their host galaxies. The Milky Way BIND Cepheids are also the closest such systems, being located at about 11 and 26 kpc.

This work shows details of an evaluation of an observational system comprising a CMOS detector, 60-cm telescope, and filter complement. The system's photometric precision and differential photometric precision, and extinction coefficients were assessed through observations of Supersky flat fields, open clusters, standard stars, and exoplanets. Photometry was precision achieved at the 0.02 mag level, while differential photometry of 0.004 mag precision. Extinction was found to be agreed with previous studies conducted at Xinglong Observatory. Ultimately, the results demonstrate this observing system is capable of precision scientific observations with CCD across the optical wavelengths.

Observations of five Y dwarfs with three optical and near-infrared instruments at the 10.4 m Gran Telescopio Canarias are reported. Deep images of the five targets and a low-resolution far-red optical spectrum for one of the targets were obtained. One of the Y dwarfs, WISE J173835+273258 (Y0), was clearly detected in the optical (z- and i-bands) and another, WISE J182831+265037 (Y2), was detected only in the z-band. We measured the colours of our targets and found that the z-J and i-z colours of the Y dwarfs are bluer than those of mid- and late-T dwarfs. This optical blueing has been predicted by models, but our data indicates that it is sharper and happens at temperatures about 150 K warmer than expected. Likely, the culprit is the K I resonance doublet, which weakens more abruptly in the T- to Y-type transition than expected. We show that the alkali resonance lines (Cs I and K I) are weaker in Y dwarfs than in T dwarfs; the far-red optical spectrum of WISE J173835+273258 is similar to that of late-T dwarfs, but with stronger methane and water features; and we noted the appearance of new absorption features that we propose could be due to hydrogen sulphide. The optical properties of Y dwarfs presented here pose new challenges to understanding grain sedimentation in extremely cool objects. The weakening of the very broad K I resonance doublet due to condensation in dust grains is more abrupt than theoretically anticipated. Consequently, the observed blueing of the z-J and i-z colours of Y dwarfs with respect to T dwarfs is more pronounced than predicted by models and could boost the potential of upcoming deep large-area optical surveys regarding their ability to detect extremely cool objects

To understand the magnetic fields of the polar crown filaments (PCFs) at high latitudes near polar regions of the Sun, we perform magnetofrictional numerical simulations on the long-term magnetic evolution of bipolar fields with roughly east-west polarity inversion lines (PILs) in a three-dimensional (3D) spherical wedge domain near polar regions. The Coriolis effect induced vortical motions at the boundaries of several supergranular cells inject magnetic helicity from the photospheric boundary into the solar atmosphere. Supergranular-scale helicity injection, transfer, and condensation produce strongly sheared magnetic fields. Magnetic reconnections at footpoints of the sheared fields produce magnetic flux ropes (MFRs) with helicity signs consistent with the observed Hemispheric Helicity Rule (HHR). The cross-sectional area of MFRs exhibits an uneven distribution, resembling a "foot-node-foot" periodic configuration. Experiments with different tilt directions of PILs indicate that the PCFs preferably form along PILs with the western end close to the polar region. The bending of PILs caused by supergranular flows, forming S-shape (Z-shape) PIL segments, promotes the formation of dextral (sinistral) MFRs. The realistic magnetic models we got can serve as starting points for the study of the plasma formation and eruption of PCFs.

Uranus and Neptune are commonly considered ice giants, and it is often assumed that, in addition to a solar mix of hydrogen and helium, they contain roughly twice as much water as rock. This classical picture has led to successful models of their internal structure and has been understood to be compatible with the composition of the solar nebula during their formation (Reynolds and Summers 1965; Podolak and Cameron 1974; Podolak and Reynolds 1984; Podolak et al. 1995; Nettelmann et al. 2013). However, the dominance of water has been recently questioned (Teanby et al. 2020; Helled and Fortney 2020; Podolak et al. 2022). Planetesimals in the outer solar system are composed mainly of refractory materials, leading to an inconsistency between the icy composition of Uranus and Neptune and the ice-poor planetesimals they accreted during formation (Podolak et al. 2022). Here we elaborate on this problem, and propose a new potential solution. We show that chemical reactions between planetesimals dominated by organic-rich refractory materials and the hydrogen in gaseous atmospheres of protoplanets can form large amounts of methane 'ice'. Uranus and Neptune could thus be compatible with having accreted refractory-dominated planetesimals, while still remaining icy. Using random statistical computer models for a wide parameter space, we show that the resulting methane-rich internal composition could be a natural solution, giving a good match to the size, mass and moment of inertia of Uranus and Neptune, whereas rock-rich models appear to only work if a rocky interior is heavily mixed with hydrogen. Our model predicts a lower than solar hydrogen to helium ratio, which can be tested. We conclude that Uranus, Neptune and similar exoplanets could be methane-rich, and discuss why Jupiter and Saturn cannot.

We implement a sub-resolution prescription for the unresolved dynamical friction onto black holes (BHs) in the OpenGadget3 code. We carry out cosmological simulations of a volume of 16 cMpc3 and zoom-ins of a galaxy group and of a galaxy cluster. The advantages of our new technique are assessed in comparison to commonly adopted methods to hamper spurious BH displacements, i.e. repositioning onto a local minimum of the gravitational potential and ad-hoc boosting of the BH particle dynamical mass. The newly-introduced dynamical friction correction provides centering of BHs on host halos which is at least comparable with the other techniques. It predicts half as many merger events with respect to the repositioning prescription, with the advantage of being less prone to leave sub-structures without any central BH. Simulations featuring our dynamical friction prescription produce a smaller (by up to 50% with respect to repositioning) population of wandering BHs and final BH masses in good agreement with observations. As for individual BH-BH interactions, our dynamical friction model captures the gradual inspiraling of orbits before the merger occurs. By contrast, the repositioning scheme, in its most classical renditions considered, describes extremely fast mergers, while the dynamical mass misrepresents the BHs' dynamics, introducing numerical scattering between the orbiting BHs. Given its performances in describing the centering of BHs within host galaxies and the orbiting of BH pair before their merging, our dynamical friction correction opens interesting applications for an accurate description of the evolution of BH demography within cosmological simulations of galaxy formation at different cosmic epochs and within different environments.

In the Orion Nebula Cluster (ONC), protoplanetary disks exhibit ionized gas clouds in the form of a striking teardrop shape as massive stars irradiate the disk material. We present the first spatially and spectrally resolved observations of 12 proplyds, using Integral Field Spectroscopy observations performed with the MUSE instrument in Narrow Field Mode (NFM) on the VLT. We present the morphology of the proplyds in seven emission lines and measure the radius of the ionization front of the targets in four tracers, covering transitions of different ionization states for the same element. We also derive stellar masses for the targets. The measurements follow a consistent trend of increasing I-front radius for a decreasing strength of the far-UV radiation as expected from photoevaporation models. By analyzing the ratios of the I-front radii as measured in the emission lines of Ha, [OI] 6300, [OII] 7330, and [OIII] 5007, we observe the ionization stratification, that is, the most ionized part of the flow being the furthest from the disk (and closest to the UV source). The ratios of I-front radii scale in the same way for all proplyds in our sample regardless of the incident radiation. We show that the stratification can help constrain the densities near the I-front by using a 1D photoionization model. We derive the upper limits of photoevaporative mass-loss rates by assuming ionization equilibrium, and estimate values decreasing towards lower impinging radiation. We do not find a correlation between Mloss and stellar mass. The highest mass-loss rate is for the proplyd 244-440. These values of Mloss, combined with estimates of the disk mass with ALMA, confirm previous estimates of the short lifetime of these proplyds. This work demonstrates the potential of this dataset and offers a new set of observables to be used to test current and future models of external photoevaporation.

The Crab Nebula is a unique laboratory for studying the acceleration of electrons and positrons through their non-thermal radiation. Observations of very-high-energy $\gamma$ rays from the Crab Nebula have provided important constraints for modelling its broadband emission. We present the first fully self-consistent analysis of the Crab Nebula's $\gamma$-ray emission between 1 GeV and $\sim$100 TeV, that is over five orders of magnitude in energy. Using the open-source software package Gammapy, we combine 11.4 yr of data from the Fermi Large Area Telescope and 80 h of High Energy Stereoscopic System (H.E.S.S.) data at the event level and provide a measurement of the spatial extension of the nebula and its energy spectrum. We find evidence for a shrinking of the nebula with increasing $\gamma$-ray energy. Furthermore, we fit several phenomenological models to the measured data, finding that none of them can fully describe the spatial extension and the spectral energy distribution at the same time. Especially the extension measured at TeV energies appears too large when compared to the X-ray emission. Our measurements probe the structure of the magnetic field between the pulsar wind termination shock and the dust torus, and we conclude that the magnetic field strength decreases with increasing distance from the pulsar. We complement our study with a careful assessment of systematic uncertainties.

Detecting an atmosphere on nearby temperate planets is one of the most important scientific objectives of the Webb mission, an endeavour in practice limited to a handful of well-characterized planets: Trappist-1d, e, f, g, LHS1140b, and the mini-Neptune K2- 18b. The first 18 months of atmospheric characterization with JWST have confirmed both its power and versatility to probe exoplanet atmospheres, and have highlighted the challenge of stellar activity in studying those atmospheres through transmission spectroscopy. Assessing the prevalence of atmosphere in temperate planets with a minimal degree of confidence will require a multi-cycle program of order of a few 1000 hours involving both eclipse photometry and transmission spectroscopy, a program that would have to be executed over a significant fraction of JWST's lifetime. The forthcoming 500 hours of Cycle 3 Director Discretionary Time dedicated to exoplanet programs represent a unique opportunity to initiate a deep reconnaissance of habitability of the best keystone temperate planets.

Solar filament eruptions are often characterized by stepwise evolution due to the involvement of multiple mechanisms, such as magnetohydrodynamic instabilities and magnetic reconnection. In this article, we investigated a confined filament eruption with a distinct two-stage evolution by using the imaging and spectroscopic observations from the Interface Region Imaging Spectrograph (IRIS) and the Solar Dynamics Observatory (SDO). The eruption originated from a kinked filament thread that separated from an active region filament. In the first stage, the filament thread rose slowly and was obstructed due to flux pile-up in its front. This obstruction brought the filament thread into reconnection with a nearby loop-like structure, which enlarged the flux rope and changed its connectivity through the foot-point migration. The newly formed flux rope became more kink unstable and drove the rapid eruption in the second stage. It ascended into the upper atmosphere and initiated the reconnection with the overlying field. Finally, the flux rope was totally disintegrated, producing several solar jets along the overlying field. These observations demonstrate that the external reconnection between the flux rope and overlying field can destroy the flux rope, thus playing a crucial role in confining the solar eruptions.

Radio signals propagating via the solar corona and solar wind are significantly affected by compressive waves, impacting solar burst properties as well as sources viewed through the turbulent atmosphere. While static fluctuations scatter radio waves elastically, moving, turbulent or oscillating density irregularities act to broaden the frequency of the scattered waves. Using a new anisotropic density fluctuation model based on solar radio bursts, we deduce the plasma velocities required to explain observations of spacecraft signal frequency broadening. The frequency broadening is consistent with motions that are dominated by the solar wind at distances $\gtrsim 10$ $R_\odot$, but the levels of frequency broadening for $\lesssim 10$ $R_\odot$ require additional radial speeds $\sim (100-300)$ km s$^{-1}$ and/or transverse speeds $\sim (20-70)$ km s$^{-1}$. The inferred radial velocities appear consistent with the sound or proton thermal speeds, while the speeds perpendicular to the radial direction are consistent with non-thermal motions measured via coronal Doppler-line broadening, interpreted as Alfv\'enic fluctuations. Landau damping of parallel propagating ion-sound (slow MHD) waves allow an estimate of the proton heating rate. The energy deposition rates due to ion-sound wave damping peak at a heliocentric distance of $\sim(1-3)$ $R_\odot$ are comparable to the rates available from a turbulent cascade of Alfv\'enic waves at large scales, suggesting a coherent picture of energy transfer, via the cascade or/and parametric decay of Alfv\'en waves to the small scales where heating takes place.

GHZ2/GLASS-z12 has been recently observed by JWST with both NIRSpec and MIRI spectrographs, making it the most distant galaxy ($z_{spec}=12.34$) with complete spectroscopic coverage from rest-frame UV to optical. It is identified as a strong CIV$_{1549}$ emitter with many other detected emission lines (NIV], HeII, OIII], NIII], CIII], [OII], [NeIII], [OIII], and H$\alpha$), including a remarkable OIII$_{1333}$ Bowen fluorescence line. We analyze in this paper the joint NIRSpec+MIRI spectral data set. Combining six optical diagnostics (namely R2, R3, R23, O32, Ne3O2, and Ne3O2Hd), we find extreme ionization conditions, with O32 $=1.39 \pm 0.19$ and Ne3O2 $=0.37 \pm 0.18$ in stark excess compared to typical values in the ISM at lower redshifts. These line properties are compatible either with an AGN or with a compact, very dense star-forming environment ($\Sigma_{\rm SFR}$ $\sim 10^2$-$10^3$ Msun/yr/kpc$^2$), with a high ionization parameter ($\log_{10}$(U) $=-1.75 \pm 0.16$), a high ionizing photon production efficiency $\log(\xi_{\rm ion}) = 25.7_{-0.1}^{+0.2}$, and a low, although not pristine, metal content ranging between $5\%$ and $11\%$ Z$_\odot$, indicating a rapid enrichment of the ISM in the last few Myrs. These properties also suggest that a substantial amount of ionizing photons ($\sim 10\%$) are leaking outside. The general lessons learned from GHZ2 are the following: (i) the UV to optical combined nebular indicators are broadly in agreement with UV-only or optical-only indicators. (ii) UV+optical diagnostics fail to discriminate between an AGN and star-formation in a low metallicity, high density, and extreme ionization environment. (iii) comparing the nebular line ratios with local analogs may be approaching its limits at $z \gtrsim 10$, as this approach is potentially challenged by the unique conditions of star formation experienced by galaxies at these extreme redshifts.

We have conducted an analysis of nebulae around Wolf-Rayet (WR) stars in M33 using data collected by the imaging Fourier transform spectrometer SITELLE at the Canada-France-Hawaii telescope as part of the SIGNALS Large Program. Of the 211 known Wolf-Rayet stars in M33, 178 are located in the fields observed in this study. We present the results of this analysis in the form of a comprehensive summary of all nebulae found around the observed WR stars. Based on three criteria we find to be the most effective for their detection, we detect a clear association with a circumstellar bubble around 33 of them (19\%). Our results show that the presence of bubbles does not correlate with the spectral type of the central star. The mean diameter of the WR nebulae we have found is 21 parsecs.

Hierarchical massive quadruple systems are ideal laboratories for examining the theories of star formation, dynamical evolution, and stellar evolution. The successive mergers of hierarchical quadruple systems might explain the mass gap between neutron stars and black holes. Looking for light curves of O-type binaries identified by LAMOST, we find a (2+2) quadruple system: TYC 3340-2437-1, located in the stellar bow-shock nebula (SBN). It has a probability of over 99.99\% being a quadruple system derived from the surface density of the vicinity stars. Its inner orbital periods are 3.390602(89) days and 2.4378(16) days, respectively, and the total mass is about (11.47 + 5.79) + (5.2 + 2.02) = 24.48 $M_{\odot}$. The line-of-sight inclinations of the inner binaries, B$_1$ and B$_2$, are 55.94 and 78.2 degrees, respectively, indicating that they are not co-planar. Based on observations spanning 34 months and the significance of the astrometric excess noise ($D>2$) in Gaia DR3 data, we guess that its outer orbital period might be a few years. If it were true, the quadruple system might form through the disk fragmentation mechanism with outer eccentric greater than zero. This eccentricity could be the cause of both the arc-like feature of the SBN and the noncoplanarity of the inner orbit. The outer orbital period and outer eccentric could be determined with the release of future epoch astrometric data of Gaia.

The dwarf planet Haumea is one of the most compelling transneptunian objects (TNOs) to study, hosting two small, dynamically interacting satellites, a family of nearby spectrally unique objects, and a ring system. Haumea itself is extremely oblate due to its 3.9 hour rotation period. Understanding the orbits of Haumea's satellites, named Hi'iaka and Namaka, requires detailed modeling of both satellite-satellite gravitational interactions and satellite interactions with Haumea's nonspherical gravitational field (parameterized here as $J_2$). Understanding both of these effects allows for a detailed probe of the satellites' masses and Haumea's $J_2$ and spin pole. Measuring Haumea's $J_2$ provides information about Haumea's interior, possibly determining the extent of past differentation. In an effort to understand the Haumea system, we have performed detailed non-Keplerian orbit fitting of Haumea's satellites using a decade of new ultra-precise observations. Our fits detect Haumea's $J_2$ and spin pole at $\gtrsim2.5\sigma$ confidence. Degeneracies present in the dynamics prevent us from precisely measuring Haumea's $J_2$ with the current data, but future observations should enable a precise measurement. Our dynamically determined spin pole shows excellent agreement with past results, illustrating the strength of non-Keplerian orbit fitting. We also explore the spin-orbit dynamics of Haumea and its satellites, showing that axial precession of Hi'iaka may be detectable over decadal timescales. Finally, we present an ephemeris of the Haumea system over the coming decade, enabling high-quality observations of Haumea and its satellites for years to come.

About 40 transneptunian binaries (TNBs) have fully determined orbits with about 10 others being solved except for breaking the mirror ambiguity. Despite decades of study almost all TNBs have only ever been analyzed with a model that assumes perfect Keplerian motion (e.g., two point masses). In reality, all TNB systems are non-Keplerian due to non-spherical shapes, possible presence of undetected system components, and/or solar perturbations. In this work, we focus on identifying candidates for detectable non-Keplerian motion based on sample of 45 well-characterized binaries. We use MultiMoon, a non-Keplerian Bayesian inference tool, to analyze published relative astrometry allowing for non-spherical shapes of each TNB system's primary. We first reproduce the results of previous Keplerian fitting efforts with MultiMoon, which serves as a comparison for the non-Keplerian fits and confirms that these fits are not biased by the assumption of a Keplerian orbit. We unambiguously detect non-Keplerian motion in 8 TNB systems across a range of primary radii, mutual orbit separations, and system masses. As a proof of concept for non-Keplerian fitting, we perform detailed fits for (66652) Borasisi-Pabu, possibly revealing a $J_2 \approx 0.44$, implying Borasisi (and/or Pabu) may be a contact binary or an unresolved compact binary. However, full confirmation of this result will require new observations. This work begins the next generation of TNB analyses that go beyond the point mass assumption to provide unique and valuable information on the physical properties of TNBs with implications for their formation and evolution.

Many details of the formation and evolution of the solar system are best inferred by understanding the orbital and physical properties of small bodies in the solar system. For example, small body binaries are particularly valuable for measuring masses. By extending the models of small body binaries beyond point masses, new information about the shape and spin orientation becomes available. This is particularly informative for Trans-Neptunian multiples (two or more components) where shapes and spin orientations are poorly understood. Going beyond point masses requires modeling tools that no longer assume fixed Keplerian orbits. To this end, we have developed a new n-quadrupole integrator SPINNY (SPIN+N-bodY) and pair it with a Bayesian parameter inference tool MultiMoon, both of which are publicly available. We describe these tools and how they can be used to learn more about solar system small body multiple systems. We then apply them to the unique Trans-Neptunian hierarchical triple system (47171) Lempo, finding a three-point-mass solution for the first time. This solution has two surprises: unequal densities of the inner components and a dynamical configuration apparently unstable on the age of the solar system.

Dynamically studying Trans-Neptunian Object (TNO) binaries allows us to measure masses and orbits. Most of the known objects appear to have only two components, except (47171) Lempo which is the single known hierarchical triple system with three similar-mass components. Though hundreds of TNOs have been imaged with high-resolution telescopes, no other hierarchical triples (or trinaries) have been found among solar system small bodies, even though they are predicted in planetesimal formation models such as gravitational collapse after the streaming instability. By going beyond the point-mass assumption and modeling TNO orbits as non-Keplerian, we open a new window into the shapes and spins of the components, including the possible presence of unresolved ``inner'' binaries. Here we present evidence for a new hierarchical triple, (148780) Altjira (2001 UQ$_{18}$), based on non-Keplerian dynamical modeling of the two observed components. We incorporate two recent Hubble Space Telescope (HST) observations, leading to a 17 year observational baseline. We present a new open-source Bayesian Point Spread Function (PSF) fitting code called nPSF that provides precise relative astrometry and uncertainties for single images. Our non-Keplerian analysis measures a statistically-significant ($\sim$2.5-$\sigma$) non-spherical shape for Altjira. The measured $J_2$ is best explained as an unresolved inner binary and an example hierarchical triple model gives the best fit to the observed astrometry. Using an updated non-Keplerian ephemeris (which is significantly different from the Keplerian predictions), we show that the predicted mutual event season for Altjira has already begun with several excellent opportunities for observations through $\sim$2030.

For the first time, we investigate the non-perturbative dynamics of single field inflation with a departure from slow-roll. Using simulations, we find that oscillatory features in the potential can drastically alter the course of inflation, with major phenomenological implications. In certain cases, the entire Universe gets trapped in a forever inflating de Sitter state. In others, only some regions get stuck in a false vacuum, offering an alternative channel for primordial black hole formation. Analogous to the flap of a butterfly, these results show that small-scale phenomena can have profound consequences on the evolution of the entire Universe. This demonstrates the necessity of a non-perturbative approach in the exploration of the small-scale physics of inflation, particularly in the regime relevant for gravitational-wave astronomy.

The Milky Way is widely considered to exhibit features of a rotational bar or quadrupole bar. In either case, the feature of the resonance of the Galactic bar should be present in the properties of the chemistry and kinematics, over a large area of the disk. With a sample of over 170,000 red clump (RC) stars from LAMOST-APOGEE data, we attempt to detect the chemical and kinematic signatures of the resonances of the Galactic bar, within 4.0 $\leq$ $R$ $\leq$ 15.0 kpc and $|Z|$ $\leq$ 3.0 kpc. The measurement of the $\Delta$[Fe/H]/$\Delta|Z|$ $-$ $R$ with subtracted the global profiles trends, shows that the thin and thick disks values Cor_$\Delta$[Fe/H]/$\Delta|Z|$ = 0.010 $\mathrm{sin}$ (1.598 $R$ + 2.551) and Cor_$\Delta$[Fe/H]/$\Delta|Z|$ = 0.006 $\mathrm{sin}$ (1.258 $R$ $-$ 0.019), respectively. The analysis of the tilt angle of the velocity ellipsoid indicates that the thin and thick disks are accurately described as $\alpha$ = $\alpha_{0}$ arctan (Z/R), with $\alpha_{0}$ = 0.198 $\mathrm{sin}$ (0.853 $R$ + 1.982) + 0.630 and $\alpha_{0}$ = 0.220 $\mathrm{sin}$ (0.884 $R$ + 2.012) + 0.679 for thin and thick disks, respectively. These periodic oscillations in Cor_$\Delta$[Fe/H]/$\Delta|Z|$ and $\alpha_{0}$ with $R$ appear in both thin and thick disks, are the most likely chemical and kinematic signatures of the resonance of the Galactic bar. The difference in the phase of the functions of the fitted periodic oscillations for the thin and thick disks may be related to the presence of a second Galactic bar.

The protostar IRAS 15398-3359 is associated with a bipolar molecular outflow ejected in an nearly northeast-southwest (NE-SW) direction which has been extensively studied. It has been suggested previous episodic accretion events by this source. Furthermore, the analysis of the morphology and kinematics of the molecular outflow revealed the presence of four $^{12}$CO(2-1) bipolar elliptical shock-like structures identified in both lobes. These structures seem to trace different ejections inclined $\sim$10\deg on the plane of the sky from each other. This led to the hypothesis that the outflow axis likely precesses and launches material episodically. We analyze ALMA archive observations in Band 6, revealing the presence of low-velocity ($<3.5$km s$^{-1}$) emission from the line $^{12}$CO(2-1) to the south and north of the protostar. We study the morphology and kinematics of the gas, which seems to support the hypothesis of a precessing episodic outflow. The ALMA observations reveal a north-south (N-S) outflow most likely associated with the IRAS 15398-3359 protostellar system. This outflow could be older than the well-studied NE-SW outflow. The orientation of the N-S outflow is 50\deg - 60\deg on the plane of the sky away from that of the NE-SW outflow. We also analyze the Spectral Energy Distribution of a far away young star and preliminary discard it as the driver of the SE outflow remnants. The new observations support the hypothesis of strong episodic accretion-ejection events in IRAS 15398-3359, accompanied by dramatic changes in the orientation of its ejection axis, implying that all the outflows in the region may have been driven by the same protostar.

As the brightest gamma-ray burst ever observed, GRB 221009A provided a precious opportunity to explore spectral line features. In this paper, we performed a comprehensive spectroscopy analysis of GRB 221009A jointly with GECAM-C and Fermi/GBM data to search for emission and absorption lines. For the first time we investigated the line feature throughout this GRB including the most bright part where many instruments suffered problems, and identified prominent emission lines in multiple time intervals. The central energy of the Gaussian emission line evolves from about 37 MeV to 6 MeV, with a nearly constant ratio (about 10\%) between the line width and central energy. Particularly, we find that both the central energy and the energy flux of the emission line evolve with time as a power law decay with power law index of -1 and -2 respectively. We suggest that the observed emission lines most likely originate from the blue-shifted electron positron pair annihilation 511 keV line. We find that a standard high latitude emission scenario cannot fully interpret the observation, thus we propose that the emission line comes from some dense clumps with electron positron pairs traveling together with the jet. In this scenario, we can use the emission line to directly, for the first time, measure the bulk Lorentz factor of the jet ($\Gamma$) and reveal its time evolution (i.e. $\Gamma \sim t^{-1}$) during the prompt emission. Interestingly, we find that the flux of the annihilation line in the co-moving frame keeps constant. These discoveries of the spectral line features shed new and important lights on the physics of GRB and relativistic jet.

Context. Investigating the diagnostic power of H$\alpha$/H$\beta$ Charge-Exchange (CE) emission as Low-Energy Galactic Cosmic Rays(LECRs) tracer in diffuse regions. Aims. In this work, we define and test a spectroscopic indicator of CE reactions between LECRs protons and neutral hydrogen atoms of the diffuse medium. This indicator can be used for mapping LECRs density in diffuse clouds and can lead to the identification of new LECRs sources as we expect density variations caused by the distance between an observed cloud and the nearest site of particle acceleration. We also lay the foundations for the definition of a photometric indicator to be used in the next full-sky photometric surveys such as the Vera Rubin 10-year Legacy Survey of Space and Time (LSST). Methods. Based on literature cross-sections, we calculate H$\alpha$/H$\beta$ line profile ratio in the case of CE and compare it with the recombination ratio. We then test our results on the Balmer-dominated filaments of the SNR RCW 86 and we explore how the spectroscopic constraints can turn into a photometric indicator based on colour indices. Results. We find that, in shocked environments, CE between LECRS and neutral hydrogen become the dominant process for Balmer lines emission. The hydrogen spectroscopic emission is expected to be modified, with respect to the recombination Balmer decrement,to result in double the H$\alpha$/H$\beta$ with respect to a similar but quiescent region. The test on the known Balmer-dominated filaments of the SNR RCW 86 confirm the efficiency of our spectroscopic indicator. Therefore we explore possible conversions of the spectroscopic indicator into colour indices combinations. This is the first step toward the definition and test of a photometric indicator for tracing LECRs to be applied in the LSST pipelines to photometrically identify new LECRs accelerators in the whole Galaxy.

We report a new hot-PAGB star identified in the Galactic GC E3, which is one of the first to show a binary signature among the identified PAGB stars of GCs. We present a detailed photometric and spectroscopic analysis to study the evolutionary status of the newly identified hot-PAGB star. We use Gaia EDR3 proper motion (PM) and parallax measurements to confirm the cluster membership. We supplement the photometric observations with radial velocities (RVs) from multi-epoch high resolution (R$\sim$28000) spectroscopic observations. We fitted the spectral energy distribution of 30 photometric fluxes of the PAGB star from ultraviolet to infra-red passbands with the TLUSTY BSTAR2006 stellar atmosphere model to derive its physical parameters, \Teff, $\log g$, \Lbol, and R. We derive the chemical abundances of nine elements (He, C, N, O, Ne, Al, Si, S, and Fe) using spectroscopic measurements. The derived chemical abundances, kinematic information, and stellar parameters are used to study the evolution history of the star. The PM and parallax of the identified PAGB star match well with cluster parameters, which confirms its cluster membership. We find that the RVs vary over $\sim$6 \kms\ between the two epochs. This is an indication of the star being in a binary orbit. A simulation of possible binary systems with the observed RVs suggests a binary period of either 39.12 days or 17.83 days with mass ratio, q$\geq$1.0. The [Fe/H] derived using the high-resolution spectra is $\sim -$0.7 dex matches well with the cluster metallicity. Various PAGB evolutionary tracks on the H-R diagram suggest a current mass of the identified PAGB star in the range 0.51$-$0.55 \Msun. The star is enriched with C and O abundances, showing similar CNO abundances compared to the other PAGB stars of GCs with the evidence of 3rd dredge-up on the AGB phase.

The correlation between interstellar turbulent speed and local star formation rate surface density, Sigma_SFR, is studied using CO observations in the PHANGS survey. The local velocity dispersion of molecular gas, sigma, increases with Sigma_SFR, but the virial parameter, alpha_vir, is about constant, suggesting the molecular gas remains self-gravitating. The correlation arises because sigma depends on the molecular surface density, Sigma_mol, and object cloud mass, M_mol, with the usual molecular cloud correlations, while Sigma_SFR increases with both of these quantities because of a nearly constant star formation efficiency for CO. Pressure fluctuations with Delta Sigma_SFR are also examined. Azimuthal variations of molecular pressure, Delta P_mol, have a weaker correlation with Delta Sigma_SFR than expected from the power-law correlation between the total quantities, suggesting slightly enhanced SFR efficiency per molecule in spiral arms. Dynamical equilibrium pressure and star formation rate correlate well for the whole sample, as P_DE~Sigma_SFR^1.3, which is steeper than in other studies. The azimuthal fluctuations, Delta P_DE(Delta Sigma_SFR), follow the total correlation P_DE(Sigma_SFR) closely, hinting that some of this correlation may be a precursor to star formation, rather than a reaction. Galactic dynamical processes correlate linearly such that Sigma_SFR~(Sigma_gas R)^(1.0\pm0.3) for total gas surface density Sigma_gas and galactic dynamical rates, R, equal to kappa, A, or Omega, representing epicyclic frequency, shear rate A, and orbit rate Omega. These results suggest important roles for both feedback and galactic dynamics.

The interstellar medium (ISM) in starburst galaxies contains plenty of chemical elements synthesised by core-collapse supernova explosions. By measuring the abundances of these metals, we can study the chemical enrichment within galaxies and the transportation of metals into circumgalactic environments through powerful outflows. We perform the spectral analysis of the X-ray emissions from the M82 core using the Reflection Grating Spectrometer (RGS) onboard XMM-Newton to accurately estimate the metal abundances in the ISM. We analyse over 300 ks of RGS data observed with fourteen position angles, covering an 80 arcsec cross-dispersion width. We employ multi-temperature thermal plasma components in collisional ionisation equilibrium (CIE) to reproduce the observed spectra, each exhibiting different spatial broadenings. The O vii band CCD image shows a broader distribution compared to those for O viii and Fe-L bands. The O viii line profiles have a prominent double-peaked structure, corresponding to the northward and southward outflows. The O vii triplet feature exhibits marginal peaks, and a single CIE component, convolved with the O vii band image, approximately reproduces the spectral shape. Combining a CIE model with a charge-exchange emission model also successfully reproduces the O vii line profiles. However, the ratio of these two components varies significantly with the observed position angles, which is physically implausible. Spectral fitting of the broadband spectra suggests a multi-temperature phase in the ISM, approximated by three components at 0.1, 0.4, and 0.7 keV. Notably, the 0.1 keV component exhibits a broader distribution than the 0.4 and 0.7 keV plasmas. The derived abundance pattern shows super-solar N/O, solar Ne/O and Mg/O, and half-solar Fe/O ratios. These results indicate the chemical enrichments by core-collapse supernovae in starburst galaxies.

We present results of a spectral survey towards a dense molecular condensation and young stellar objects (YSOs) projected on the border of the HII region RCW 120 and discuss emission of 20 molecules which produce the brightest lines. The survey was performed with the APEX telescope in the frequency range 200 -- 260 GHz. We provide evidences for two outflows in the dense gas. The first one is powered by the RCW 120 S2 YSO and oriented along the line of sight. The second outflow around RCW 120 S1 is aligned almost perpendicular to the line of sight. We show that area with bright emission of CH$_3$OH, CH$_3$CCH and CH$_3$CN are organised into an onion-like structure where CH$_3$CN traces warmer regions around the YSOs than the other molecules. Methanol seems to be released to the gas phase by shock waves in the vicinity of the outflows while thermal evaporation still does not work towards the YSOs. We find only a single manifestation of the UV radiation to the molecules, namely, enhanced abundances of small hydrocarbons CCH and c-C$_3$H$_2$ in the photo-dissociation region.

Recent developments in multi-dimensional simulations of core-collapse supernovae have considerably improved our understanding of this complex phenomenon. In addition to that, one-dimensional (1D) studies have been employed to study the explosion mechanism and its causal connection to the pre-collapse structure of the star, as well as to explore the vast parameter space of supernovae. Nonetheless, many uncertainties still affect the late stages of the evolution of massive stars, their collapse, and the subsequent shock propagation. In this review, we will briefly summarize the state-of-the-art of both 1D and 3D simulations and how they can be employed to study the evolution of massive stars, supernova explosions, and shock propagation, focusing on the uncertainties that affect each of these phases. Finally, we will illustrate the typical nucleosynthesis products that emerge from the explosion.

Large-scale structure surveys have reported measurements of the density of matter, $\Omega_\mathrm{m}$, and the amplitude of clustering, $\sigma_8$, that are in tension with the values inferred from observations of the cosmic microwave background. While this may be a sign of new physics that slows the growth of structure at late times, strong astrophysical feedback processes could also be responsible. In this work, we argue that astrophysical processes are not independent of cosmology and that their coupling naturally leads to stronger baryonic feedback in cosmological models with suppressed structure formation or when combined with a mechanism that removes dark matter from halos. We illustrate this with two well-motivated extensions of the Standard Model known to suppress structure formation: massive neutrinos and decaying dark matter. Our results, based on the FLAMINGO suite of hydrodynamical simulations, show that the combined effect of baryonic and non-baryonic suppression mechanisms is greater than the sum of its parts, particularly for decaying dark matter. We also show that the dependence of baryonic feedback on cosmology can be modelled as a function of the ratio $f_\mathrm{b}/c^2_\mathrm{v}\sim f_\mathrm{b}/(\Omega_\mathrm{m}\sigma_8)^{1/4}$ of the universal baryon fraction, $f_\mathrm{b}$, to a velocity-based definition of halo concentration, $c^2_\mathrm{v}$, giving an accurate fitting formula for the baryonic suppression of the matter power spectrum. Although the combination of baryonic and non-baryonic suppression mechanisms can resolve the tension, the models with neutrinos and decaying dark matter are challenged by constraints on the expansion history.

We study the physics of photon rings in a wide range of axisymmetric black holes admitting a separable Hamilton-Jacobi equation for the geodesics. Utilizing the Killing-Yano tensor, we derive the Penrose limit of the black holes, which describes the physics near the photon ring. The obtained plane wave geometry is directly linked to the frequency matrix of the massless wave equation, as well as the instabilities and Lyapunov exponents of the null geodesics. Consequently, the Lyapunov exponents and frequencies of the photon geodesics, along with the quasinormal modes, can be all extracted from a Hamiltonian in the Penrose limit plane wave metric. Additionally, we explore potential bounds on the Lyapunov exponent, the orbital and precession frequencies, in connection with the corresponding inverted harmonic oscillators and we discuss the possibility of photon rings serving as holographic horizons in a holographic duality framework for astrophysical black holes. Our formalism is applicable to spacetimes encompassing various types of black holes, including stationary ones like Kerr, Kerr-Newman, as well as static black holes such as Schwarzschild, Reissner-Nordstr\"om, among others.

Compact stellar objects like supernovae and neutron stars are believed to cool by emitting axions predominantly via axion bremsstrahlung ($NN \to NNa$), pion conversion ($\pi^- p^+ \to N a$) and photo-production ($\gamma N \to N a$). In this $\textit{letter}$, we study a previously overlooked contribution to the photo-production channel, originating from the $\textit{unavoidable}$ anomaly-induced Wess-Zumino-Witten term $\propto \epsilon^{\mu \nu \alpha \beta}\, F_{\mu \nu}\, \partial_\alpha a \, \omega_\beta$. As a result, supernovae can cool up to three times faster than what was previously thought, translating into stronger bounds on the axion decay constant. Furthermore, the spectrum of axions emitted in the process is $\textit unique$ as it is significantly harder than those originating from bremsstrahlung. As a consequence, we find that these axions are more likely to show up in near future water Cherenkov detectors as compared to traditional ones.

A dense neutrino gas exhibiting angular crossings in the electron lepton number is unstable and develops fast flavor conversions. Instead of assuming an unstable configuration from the onset, we imagine that the system is externally driven toward instability. We use the simplest model of two neutrino beams initially of different flavor that either suddenly appear or one or both slowly build up. Flavor conversions commence well before the putative unstable state is fully attained, and the final outcome depends on how the system is driven. Our results suggest that in an astrophysical setting, one should focus less on flavor instabilities in the neutrino radiation field and more on the external dynamics that leads to the formation of the unstable state.

We compute the rate with which unobserved fields decohere other fields to which they couple, both in flat space and in de Sitter space, for scalar fields prepared in their standard adiabatic vacuum. The process is very efficient in de Sitter space once the modes in question pass outside the Hubble scale, displaying the tell-tale phenomenon of secular growth that indicates the breakdown of perturbative methods on a time scale parameterically long compared with the Hubble time. We show how to match the perturbative evolution valid at early times onto a late-time Lindblad evolution whose domain of validity extends to much later times, thereby allowing a reliable resummation of the perturbative result beyond the perturbative regime. Super-Hubble modes turn out to be dominantly decohered by unobserved modes that are themselves also super-Hubble. When applied to curvature perturbations during inflation this observation closes a potential loophole in recent calculations of the late-time purity of the observable primordial fluctuations.

The Amaterasu cosmic ray particle appears to have come from the direction of the local cosmic void. We take this as evidence that it is a magnetic monopole rather than a proton or nucleus. This in turn strongly suggests physics at high energy is described by a quiver gauge theory.

Electromagnetic observations reveal that almost all galaxies have supermassive black holes (SMBHs) at their centers, but their properties, especially their spins, are not fully understood. Some of the authors have recently shown [Oshita and Tsuna (2023)] that rapid spins of $>0.9$, inferred for masses around $10^7\ M_\odot$ from observations of local SMBHs and cosmological simulations, source {\it long-lived} ringdowns that enhance the precision of black hole spectroscopy to test gravity in the near-extreme Kerr spacetime. In this work, we estimate the statistical errors in the SMBH mass-spin inference in anticipation of the LISA's detection of extreme mass-ratio mergers. We show that for rapidly spinning SMBHs, more precise mass and spin measurements are expected due to the excitations of higher angular modes. For a near-extremal SMBH of mass $10^7M_\odot$ merging with a smaller BH with mass ratio $10^{-3}$ at a luminosity distance of $\lesssim 10\:\mathrm{Gpc}$ (redshift $z \lesssim 1.37$), the measurement errors in the mass and spin of the SMBH would be $\sim 1\:\mathrm{\%}$ and $\sim 10^{-1}\:\mathrm{\%}$ respectively.

In gravitational wave astronomy, non-Gaussian noise, such as scattered light noise disturbs stable interferometer operation, limiting the interferometer's sensitivity, and reducing the reliability of the analyses. In scattered light noise, the non-Gaussian noise dominates the sensitivity in a low frequency range of less than a few hundred Hz, which is sensitive to gravitational waves from compact binary coalescence. This non-Gaussian noise prevents reliable parameter estimation, since several analysis methods are optimized only for Gaussian noise. Therefore, identifying data contaminated by non-Gaussian noise is important. In this work, we extended the conventional method to evaluate non-Gaussian noise, Rayleigh statistic, by using a statistical hypothesis test to determine a threshold for non-Gaussian noise. First, we estimated the distribution of the Rayleigh statistic against Gaussian noise, called the background distribution, and validated that our extension serves as the hypothetical test. Moreover, we investigated the detection efficiency by assuming two non-Gaussian noise models. For example, for the model with strong scattered light noise, the true positive rate was always above 0.7 when the significance level was 0.05. The results showed that our extension can contribute to an initial detection of non-Gaussian noise and lead to further investigation of the origin of the non-Gaussian noise.

A critical discussion on the $H_0$ Hubble constant tension is presented by considering both early and late-type observations. From recent precise measurements, discrepancies emerge when comparing results for some cosmological quantities obtained at different redshifts. We highlight the most relevant measurements of $H_0$ and propose potential ideas to solve its tension. These solutions concern the exploration of new physics beyond the $\Lambda$CDM model or the evaluation of $H_0$ by other methods. In particular, we focus on the role of the look-back time.

We consider the Penrose process with the charged particles in the Reissner-Nordstr\"{o}m background. Let parent particle 0 decay to particles 1 and 2. With the assumption that all three particles move in the same plane, the exact formulas for characteristics of particles 1 and 2 in terms of those of particle 0 are derived. We concentrate on scenarios in which particle 1 and 2 are ejected along the trajectory of particle 0. It is shown that such scenarios correspond to the extrema of energies $E_{1}$ or $E_{2}$ of daughter particles with respect to the angular momentum $L_{1}$ or $L_{2}$ of one of them. We derive bounds on the values of angular momenta $L_{1}$ and $L_{2}$. We give classification of these scenarios and discuss their properties including decay in the near-horizon region. The results are reformulated in terms of velocities of daughter particles in the center of mass frame. The approach is applicable also to collisional Penrose process, in which a combination of particles 1 and 2 is considered as one effective particle. If the mass of particle 0 $m_{0}\rightarrow \infty $, the situation corresponds to the Ba\~{n}ados-Silk-West effect, the results agree with the ones known in literature before. In addition, we consider special cases when decay occurs in the turning point for one or all three particles.

Motivated by string-theoretic swampland conjectures, the existence of a dark fifth dimension, whose size is roughly 1 -- 10 microns, has been proposed. A great deal of supporting evidence has been presented, and definitive experimental tests are likely to be carried out. The basic idea is that the four-dimensional spacetime that we observe lives on a brane that is localized in the dark dimension. This short note points out that there are two distinct ways to realize such a scenario in string theory/M-theory. In the one considered previously the dark dimension is topologically a circle and our observable 4d spacetime is confined to a brane that is localized in a GUT-scale region of the circle. An alternative possibility is that the dark dimension is a line interval with branes attached at each end. This option would imply the existence of a parallel 4d spacetime microns away from us!