Abstract visibility
Abstract text size

Papers for Monday, Apr 22 2024

Itai Linial, Brian D. Metzger

24 pages, 6 figures. Submitted to ApJ. Comments welcome!

A modest fraction of the stars in galactic nuclei fed towards the central supermassive black hole (SMBH) approach on low-eccentricity orbits driven by gravitational-wave radiation (extreme mass ratio inspiral, EMRI). In the likely event that a gaseous accretion disk is created in the nucleus during this slow inspiral (e.g., via an independent tidal-disruption event; TDE), star-disk collisions generate regular short-lived flares consistent with the observed quasi-periodic eruption (QPE) sources. We present a model for the coupled star-disk evolution which self-consistently accounts for mass and thermal energy injected into the disk from stellar collisions and associated mass ablation. For weak collision/ablation heating, the disk is thermally-unstable and undergoes limit-cycle oscillations which modulate its properties and lead to accretion-powered outbursts on timescales of years to decades, with a time-averaged accretion rate $\sim 0.1 \dot{M}_{\rm Edd}$. Stronger collision/ablation heating acts to stabilize the disk, enabling roughly steady accretion at the EMRI-stripping rate. In either case, the stellar destruction time through ablation, and hence the maximum QPE lifetime, is $\sim 10^{2}-10^{3}$ yr, far longer than fall-back accretion after a TDE. The quiescent accretion disks in QPE sources may at the present epoch be self-sustaining and fed primarily by EMRI ablation. Indeed, the observed range of secular variability broadly match those predicted for collision-fed disks. Changes in the QPE recurrence pattern following such outbursts, similar to that observed in GSN 069, could arise from temporary misalignment between the EMRI-fed disk and the SMBH equatorial plane as the former regrows its mass after a state transition.

Aleksandra Olejak, Jakub Klencki, Xiao-Tian Xu, Chen Wang, Krzysztof Belczynski, Jean-Pierre Lasota

14 pages + appendix, comments welcome, to be submitted to A&A

The growing database of gravitational-wave (GW) detections with the binary black holes (BHs) merging in the distant Universe contains subtle insights into their formation scenarios. One of the puzzling properties of detected GW sources is the possible (anti)correlation between mass ratio q of BH-BH binaries and their effective spin. We use rapid binary evolution models to demonstrate that the isolated binary evolution followed by efficient tidal spin-up of stripped helium core produces a similar pattern in Xeff vs q distributions of BH-BH mergers. In our models, the progenitors of unequal BH-BH systems in the stable mass transfer formation scenario are more likely to efficiently shrink their orbits during the second Roche-lobe overflow than the binaries that evolve into nearly equal-mass component systems. This makes it easier for unequal-mass progenitors to enter the tidal spin-up regime and later merge due to GW emission. Our results are, however, sensitive to some input assumptions, especially, the stability of mass transfer and the angular momentum loss during non-conservative mass transfer. We note that mass transfer prescriptions widely adopted in rapid codes favor the formation of BH-BH merger progenitors with unequal masses and moderate separations. We compare our results with detailed stellar model grids and find reasonable agreement after appropriate calibration of the physics models. We anticipate that future detections of unequal-mass BH-BH mergers could provide valuable constraints on the role of the stable mass transfer formation channel. A significant fraction of BH-BH detections with mass ratio q in range (0.4 - 0.7) would be consistent with the mass ratio reversal scenario during the first, relatively conservative mass transfer, and a non-enhanced angular momentum loss during the second, highly non-conservative mass transfer phase.

Collin T. Christy, Kate D. Alexander, Yvette Cendes, Ryan Chornock, Tanmoy Laskar, Raffaella Margutti, Edo Berger, Michael Bietenholz, Deanne Coppejans, Fabio De Colle, Tarraneh Eftekhari, Thomas W.-S. Holoien, Tatsuya Matsumoto, James C. A. Miller-Jones, Enrico Ramirez-Ruiz, Richard Saxton, Sjoert van Velzen, Mark Wieringa

25 pages. Submitted to ApJ

We present detailed radio observations of the tidal disruption event (TDE) ASASSN-19bt/AT2019ahk, obtained with the Australia Telescope Compact Array (ATCA), the Atacama Large Millimeter/submillimeter Array (ALMA), and the MeerKAT radio telescopes, spanning 40 to 1464 days after the onset of the optical flare. We find that ASASSN-19bt displays unusual radio evolution compared to other TDEs, as the peak brightness of its radio emission increases rapidly until 457 days post-optical discovery and then plateaus. Using a generalized approach to standard equipartition techniques, we estimate the energy and corresponding physical parameters for two possible emission geometries: a non-relativistic spherical outflow and a relativistic outflow observed from an arbitrary viewing angle. We find that the non-relativistic solution implies a continuous energy rise in the outflow from $E\sim10^{46}$ erg to $E\sim10^{49}$ erg with $\beta \approx 0.05$, while the off-axis relativistic jet solution instead suggests $E\approx10^{52}$ erg with $\Gamma\sim10$ erg at late times in the maximally off-axis case. We find that neither model provides a holistic explanation for the origin and evolution of the radio emission, emphasizing the need for more complex models. ASASSN-19bt joins the population of TDEs that display unusual radio emission at late times. Conducting long-term radio observations of these TDEs, especially during the later phases, will be crucial for understanding how these types of radio emission in TDEs are produced.

Martje Slob, Mariska Kriek, Aliza G. Beverage, Katherine A. Suess, Guillermo Barro, Rachel Bezanson, Chloe M. Cheng, Charlie Conroy, Anna de Graaff, Natascha M. Förster Schreiber, Marijn Franx, Brian Lorenz, Pavel E. Mancera Piña, Danilo Marchesini, Adam Muzzin, Andrew B. Newman, Sedona H. Price, Alice E. Shapley, Mauro Stefanon, Pieter van Dokkum, Daniel R. Weisz

Submitted to ApJ; 24 pages, 13 figures, 2 tables (excluding appendices)

We present an overview and first results from the Spectroscopic Ultradeep Survey Probing Extragalactic Near-infrared Stellar Emission (SUSPENSE), executed with NIRSpec on JWST. The primary goal of the SUSPENSE program is to characterize the stellar, chemical, and kinematic properties of massive quiescent galaxies at cosmic noon. In a single deep NIRSpec/MSA configuration, we target 20 distant quiescent galaxy candidates ($z=1-3$, $H_{AB}<23$), as well as 53 star-forming galaxies at $z=1-4$. With 16hr of integration and the G140M-F100LP dispersion-filter combination, we observe numerous Balmer and metal absorption lines for all quiescent candidates. We derive stellar masses (log$M_*/M_{\odot}\sim10.3-11.5$) and detailed star-formation histories (SFHs) and show that all 20 candidate quiescent galaxies indeed have quenched stellar populations. These galaxies show a variety of mass-weighted ages ($0.8-3.0$Gyr) and star formation timescales ($\sim0.5-4$Gyr), and four out of 20 galaxies were already quiescent by $z=3$. On average, the $z>1.75$ $[z<1.75]$ galaxies formed 50% of their stellar mass before $z=4$ $[z=3]$. Furthermore, the typical SFHs of galaxies in these two redshift bins ($z_{\text{mean}}=2.2$ and $z_{\text{mean}}=1.3$) indicate that galaxies at higher redshift formed earlier and over shorter star-formation timescales compared to lower redshifts. Although this evolution is naturally explained by the growth of the quiescent galaxy population over cosmic time, we cannot rule out that mergers and late-time star formation also contribute to the evolution. In future work, we will further unravel the early formation, quenching, and late-time evolution of these galaxies by extending this work with studies on their chemical abundances,

The comparison of the properties of galaxy cluster samples selected using observations in different wavebands may shed light on potential biases of the way in which the samples are assembled. For this comparison, we introduce a new observable that does not require previous knowledge of the cluster mass: the X-ray mean surface brightness within the central 300 kpc. We found that clusters with low surface brightness, defined as those with a mean surface brightness below 43.35 \subr , are about one quarter of the whole cluster population in a sample of 32 clusters in the nearby Universe, selected independently of the intracluster medium properties. Almost no example of a low central surface brightness cluster exists instead in two X-ray selected samples, one sample based on XMM-Newton XXL-100 survey data and the other on full-depth eROSITA eFEDS data, although these clusters are known to exist in the range of redshift and mass as probed by these two surveys. Furthermore, the Sunayev-Zeldovich Atacama Cosmology Telescope cluster survey is even more selective than the previous two samples because it does not even include clusters with intermediate surface brightness, which are instead present in X-ray selected samples that explore the same volume of the Universe. Finally, a measure of the mean surface brightness, which is obtained without knowledge of the mass, proves to be effective in narrowing the number of clusters to be followed-up because it recognizes those with a low gas fraction or with a low X-ray luminosity for their mass. Identifying these would otherwise require knowledge of the mass for all clusters.

Aman Khalid, Sarah Brough, Garreth Martin, Lucas C. Kimmig, Claudia Del P. Lagos, Rhea-Silvia Remus, Cristina Martinez-Lombilla

25 pages, 14 figures, accepted in MNRAS

Tidal features provide signatures of recent mergers and offer a unique insight into the assembly history of galaxies. The Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) will enable an unprecedentedly large survey of tidal features around millions of galaxies. To decipher the contributions of mergers to galaxy evolution it will be necessary to compare the observed tidal features with theoretical predictions. Therefore, we use cosmological hydrodynamical simulations NewHorizon, EAGLE, IllustrisTNG, and Magneticum to produce LSST-like mock images of $z\sim0$ galaxies ($z\sim0.2$ for NewHorizon) with $M_{\scriptstyle\star,\text{ 30 pkpc}}\geq10^{9.5}$ M$_{\scriptstyle\odot}$. We perform a visual classification to identify tidal features and classify their morphology. We find broadly good agreement between the simulations regarding their overall tidal feature fractions: $f_{\text{NewHorizon}}=0.40\pm0.06$, $f_{\text{EAGLE}}=0.37\pm0.01$, $f_{\text{TNG}}=0.32\pm0.01$ and $f_{\text{Magneticum}}=0.32\pm0.01$, and their specific tidal feature fractions. Furthermore, we find excellent agreement regarding the trends of tidal feature fraction with stellar and halo mass. All simulations agree in predicting that the majority of central galaxies of groups and clusters exhibit at least one tidal feature, while the satellite members rarely show such features. This agreement suggests that gravity is the primary driver of the occurrence of visually-identifiable tidal features in cosmological simulations, rather than subgrid physics or hydrodynamics. All predictions can be verified directly with LSST observations.

Manuela Bischetti, Hyunseop Choi, Fabrizio Fiore, Chiara Feruglio, Stefano Carniani, Valentina D'Odorico, Eduardo Bañados, Huanqing Chen, Roberto Decarli, Simona Gallerani, Julie Hlavacek-Larrondo, Samuel Lai, Karen M. Leighly, Chiara Mazzucchelli, Laurence Perreault-Levasseur, Roberta Tripodi, Fabian Walter, Feige Wang, Jinyi Yang, Maria Vittoria Zanchettin, Yongda Zhu

Submitted to ApJ

Although the mass growth of supermassive black holes during the Epoch of Reionisation is expected to play a role in shaping the concurrent growth of their host-galaxies, observational evidence of feedback at z$\gtrsim$6 is still sparse. We perform the first multi-scale and multi-phase characterisation of black-hole driven outflows in the $z\sim6.6$ quasar J0923+0402 and assess how these winds impact the cold gas reservoir. We employ the SimBAL spectral synthesis to fit broad absorption line (BAL) features and find a powerful ionized outflow on $\lesssim210$ pc scale, with a kinetic power $\sim2-100$\% of the quasar luminosity. ALMA observations of [CII] emission allow us to study the morphology and kinematics of the cold gas. We detect high-velocity [CII] emission, likely associated with a cold neutral outflow at $\sim0.5-2$ kpc scale in the host-galaxy, and a bright extended [CII] halo with a size of $\sim15$ kpc. For the first time at such an early epoch, we accurately constrain the outflow energetics in both the ionized and the atomic neutral gas phases. We find such energetics to be consistent with expectations for an efficient feedback mechanism, and both ejective and preventative feedback modes are likely at play. The scales and energetics of the ionized and atomic outflows suggest that they might be associated with different quasar accretion episodes. The results of this work indicate that strong black hole feedback is occurring in quasars at $z\gtrsim6$ and is likely responsible for shaping the properties of the cold gas reservoir up to circum-galactic scales.

Haocheng Zhang (1, 2), Markus Böttcher (3), Ioannis Liodakis (4) ((1) University of Maryland Baltimore County Baltimore, (2) NASA Goddard Space Flight Center, (3) North-West University, (4) NASA Marshall Space Flight Center)

11 pages, 3 figures, accepted by ApJ

X-ray and MeV polarization can be powerful diagnostics for leptonic and hadronic blazar models. Previous predictions are mostly based on a one-zone framework. However, recent IXPE observations of Mrk~421 and 501 strongly favor a multi-zone framework. Thus, the leptonic and hadronic polarization predictions need to be revisited. Here we identify two generic radiation transfer effects, namely, double depolarization and energy stratification, that can have an impact on the leptonic and hadronic polarization. We show how they are generalized from previously known multi-zone effects of the primary electron synchrotron radiation. Under our generic multi-zone model, the leptonic polarization degree is expected to be much lower than the one-zone prediction, unlikely detectable in most cases. The hadronic polarization degree can reach a value as high as the primary electron synchrotron polarization during simultaneous multi-wavelength flares, consistent with the one-zone prediction. Therefore, IXPE and future X-ray and MeV polarimeters such as eXTP, COSI, and AMEGO-X, have good chances to detect hadronic polarization during flares. However, the hadronic polarization cannot be well constrained during the quiescent state. Nonetheless, if some blazar jets possess relatively stable large-scale magnetic structures, as suggested by radio observations, a non-trivial polarization degree may show up for the hadronic model after a very long exposure time ($\gtrsim 1$ year).

Haylee N. Archer, Deidre A. Hunter, Bruce G. Elmegreen, Monica Rubio, Phil Cigan, Rogier A. Windhorst, Juan R. Cortés, Rolf A. Jansen

In press in the Astronomical Journal

Wolf-Lundmark-Melotte (WLM) is a Local Group dwarf irregular (dIrr) galaxy with a metallicity 13% of solar. At 1 Mpc, the relative isolation of WLM provides a unique opportunity to investigate the internal mechanisms of star formation at low metallicities. The earliest stages of star formation in larger spirals occur in embedded clusters within molecular clouds, but dIrrs lack the dust, heavy metals, and organized structure of spirals believed necessary to collapse the molecular clouds into stars. Despite actively forming stars, the early stages of star formation in dIrrs is not well understood. We examine the relationship between early star formation and molecular clouds at low metallicities. We utilize ALMA-detected CO cores, $\textit{JWST}$ near-infrared (NIR) images (F090W, F150W, F250M, and F430M), and $\textit{GALEX}$ far-ultraviolet (FUV) images of WLM to trace molecular clouds, early star formation, and longer star formation timescales respectively. We compare clumps of NIR-bright sources (referred to as objects) categorized into three types based on their proximity to FUV sources and CO cores. We find objects, independent of their location, have similar colors and magnitudes and no discernible difference in temperature. However, we find that objects near CO have higher masses than objects away from CO, independent of proximity to FUV. Additionally, objects near CO are coincident with Spitzer 8 $\mu$m sources at a higher frequency than objects elsewhere in WLM. This suggests objects near CO may be embedded star clusters at an earlier stage of star formation, but accurate age estimates for all objects are required for confirmation.

Jets are endemic to both Galactic solar mass and extragalactic supermassive black holes. A recent 86 GHz image of M\,87 shows a jet emerging from the accretion ring around a black hole, providing the first direct observational constraint on the kinematics of the jet-launching region in any black hole jetted system. The very wide ($\sim280\mu\rm{as}$), highly collimated, limb-brightened cylindrical jet base is not predicted in current numerical simulations. The emission was shown to be consistent with that of a thick-walled cylindrical source that apparently feeds the flow that produces the bright limbs of the outer jet at an axial distance downstream of $0.4 \,\rm{mas}<z<0.65\, \rm{mas}$. The analysis here applies the conservation laws of energy, angular momentum, and magnetic flux to the combined system of the outer jet, the cylindrical jet, and the launch region. It also uses the brightness asymmetries of the jet and counterjet to constrain the Doppler factor. The only global solutions have a source that is located $<34 \mu\rm{as}$ from the event horizon. This includes the Event Horizon Telescope annulus of emission and the regions interior to this annulus. The axial jet begins as a magnetically dominated flow that spreads laterally from the launch radius ($<34 \mu\rm{as}$). It becomes super-magnetosonic before it reaches the base of the cylindrical jet. The flow is ostensibly redirected and collimated by a cylindrical nozzle formed in a thick accretion disk. The flow emerges from the nozzle as a mildly relativistic ($0.3c<v<0.4c$) jet with a significant protonic kinetic energy flux.

Hidekazu Tanaka, Kohei Okada

16 pages, 6 figures, accepted for publication in The Astrophysical Journal

We performed linear calculations to determine the Type I planetary migration rate for three-dimensional locally isothermal disks with radial temperature gradients. For 3D disks with radial temperature gradients, the linear wave equation has a divergent term of the third pole, which makes corotation a non-removal singularity. We suppressed the divergence with the Landau prescription to obtain the wave solutions. Despite the singularity at corotation, we derived a definite torque on the planet because the divergent term amplifies the waves only in the neighborhood of corotation and has little effect on the planetary torque. Consequently, we derived the formulas for the total, Lindblad, and corotation torques for locally isothermal disks. The resulting torque term due to the disk temperature gradient agrees well with the results of previous 3D hydrodynamical simulations for locally isothermal disks. Our linear calculation also provides the 3D horseshoe torque, which is close to the results of previous 3D hydrodynamical simulations.

N. Alvarez-Baena (1), R. Carrera (2, 3), H. Thompson (4), L. Balaguer-Nuñez (5, 6), A. Bragaglia (2), C. Jordi (5, 6), E. Silva-Villa (1), A. Vallenari (3) ((1) Instituto de Física, Universidad de Antioquia, Medellín, Colombia, (2) INAF-Osservatorio di Astrofisica e Scienza dello Spazio di Bologna, Bologna, Italy, (3) INAF-Osservatorio Astronomico di Padova, Padova, Italy, (4) Department of Physics, University of Surrey, Guildford, UK, (5) Institut d'Estudis Espacials de Catalunya (IEEC), Barcelona, Spain, (6) Institut de Ciències del Cosmos (ICCUB), Universitat de Barcelona (UB), Barcelona, Spain, (7) Departament de Física Quàntica i Astrofísica (FQA), Universitat de Barcelona (UB), Barcelona, Spain)

16 pages, 11 figures, accepted for publication in A&A

Context: Open clusters' dynamical evolution is driven by stellar evolution, internal dynamics and external forces, which according to dynamical simulations, will evaporate them in a timescale of about 1 Ga. However, about 10\% of the known open clusters are older. They are special systems whose detailed properties are related to their dynamical evolution and the balance between mechanisms of cluster formation and dissolution. Aims:We investigate the spatial distribution and structural parameters of six open clusters older than 1 Ga in order to constrain their dynamical evolution, and longevity.} Methods: We identify members using Gaia EDR3 data up to a distance of 150 pc from each cluster's centre. We investigate the spatial distribution of stars inside each cluster to understand their degree of mass segregation. Finally, in order to interpret the obtained radial density profiles we reproduced them using the lowered isothermal model explorer with PYTHON Limepy and spherical potential escapers stitched SPES. Results: All the studied clusters seem more extended than previously reported in the literature. The spatial distributions of three of them show some structures aligned with their orbits. They may be related to the existence of extra tidal stars. In fact, we find that about 20 % of their members have enough energy to leave the systems or are already unbound. Together with their initial masses, their distances to the Galactic plane may play significant roles in their survival. We found clear evidences that the most dynamically evolved clusters do not fill their Roche volumes, appearing more concentrated than the others. Finally, we find a cusp-core dichotomy in the central regions of the studied clusters, which shows some similarities to the one observed among globular clusters.

The evolution of Cataclysmic Variables (CVs) is driven by period-changes ($\dot{P}$), for which the long-venerable consensus is the Magnetic Braking Model (MBM). The MBM has its only distinctive assumption being a power-law `recipe' describing the angular momentum loss (AML) in the binary, producing a single unique evolutionary track with $\dot{P}$ as a function of the orbital period. This required prediction can be used to test the most-fundamental assumption of MBM, but it has never been tested previously. In this paper, I collect $\dot{P}$ measures for 52 CVs of all types. First, 44 per cent of the CVs have positive-$\dot{P}$, with such being impossible in MBM. Second, even amongst the CVs with negative-$\dot{P}$, their $\dot{P}$ measures are always more-negative than required by MBM, with an average deviation of 110$\times$. Third, three CVs have large chaotic variations in $\dot{P}$ that are impossible for MBM, proving that some unknown mechanism exists and is operating that dominates for these systems. Fourth, the MBM does not account for the long-term effects on evolution arising from the large sudden period decreases seen across many nova events, with this unaccounted effect dominating for the majority of nova systems and changing the sign of the overall evolutionary $\dot{P}$. Fifth, three recurrent novae are observed to suddenly change $\dot{P}$ by an order-of-magnitude across a nova event, with this being impossible in the MBM. In all, the required MBM $\dot{P}$ predictions all fail for my 52 CVs, usually by orders-of-magnitude, so the MBM AML-recipe is wrong by orders-of-magnitude.

Abhinna Sundar Samantaray, H.K. Jassal, Kulinder Pal Singh, G.C. Dewangan

Paper accepted for publication in the Astronomy & Astrophysics (A&A) journal

We have used AstroSat's Ultra Violet Imaging Telescope (UVIT) to obtain NUV and FUV images of NGC~7469 in various filters. We have carried out photometry of star-forming regions in the two galaxies and obtained their distribution. We have obtained the distributions of star formation rates (SFR) in NGC 7469 and IC 5283 using the estimates obtained from the FUV and NUV bands and carried out Kolmogorov-Smirnov tests to look for differences in the SFRs in the two galaxies. We have derived the spectral energy distribution (SED), leading to the determination of physical parameters, including overall SFR, stellar mass ($\text{M}_{*}$), dust mass ($\text{M}_\text{Dust}$), and specific star formation rates in both the galaxies. Our NUV and FUV images show the presence of an outer spiral arm which is better resolved. We have identified 33 new star-forming regions out of 51 total identified in the UVIT composite image. Enhanced star formation activity is observed coinciding with the interaction, and KS tests show that there are no significant differences in the SFR distributions of NGC 7469 and IC 5283 indicating that the interaction between the galaxies has not influenced their star formation processes differently. The SED plots and the photometric results show that most of the star formation activity is confined inside the central starburst (SB) ring.

Innocent Okwudili Eya, Evaristus Uzochukwu Iyida

7 figures, 32 pages

Glitch activity refers to the mean increase in pulsar spin frequency per year due to rotational glitches. It is an important tool for studying super-nuclear matter using neutron star interiors as templates. Glitch events are typically observed in the spin frequency ($\nu$) and frequency derivative ($\dot{\nu}$) of pulsars. The rate of glitch recurrence decreases as the pulsar ages, and the activity parameter is usually measured by linear regression of cumulative glitches over a given period. This method is effective for pulsars with multiple regular glitch events. However, due to the scarcity of glitch events and the difficulty of monitoring all known pulsars, only a few have multiple records of glitch events. This limits the use of the activity parameter in studying neutron star interiors with multiple pulsars. In this study, we examined the relationship between the activity parameters and pulsar spin parameters (spin frequency, frequency derivative, and pulsar characteristic age). We found that a quadratic function provides a better fit for the relationship between activity parameters and spin parameters than the commonly used linear functions. Using this information, we were able to estimate the activity parameters of other pulsars that do not have records of glitches. Our analysis shows that the relationship between the estimated activity parameters and pulsar spin parameters is consistent with that of the observed activity parameters in the ensemble of pulsars.

Dante S. Lauretta, Harold C. Connolly, Jr., Joseph E. Aebersold, Conel M. O. D. Alexander, Ronald-L. Ballouz, Jessica J. Barnes, Helena C. Bates, Carina A. Bennett, Laurinne Blanche, Erika H. Blumenfeld, Simon J. Clemett, George D. Cody, Daniella N. DellaGiustina, Jason P. Dworkin, Scott A. Eckley, Dionysis I. Foustoukos, Ian A. Franchi, Daniel P. Glavin, Richard C. Greenwood, Pierre Haenecour, Victoria E. Hamilton, Dolores H. Hill, Takahiro Hiroi, Kana Ishimaru, Fred Jourdan, Hannah H. Kaplan, Lindsay P. Keller, Ashley J. King, Piers Koefoed, Melissa K. Kontogiannis, Loan Le, Robert J. Macke, Timothy J. McCoy, Ralph E. Milliken, Jens Najorka, Ann N. Nguyen, Maurizio Pajola, Anjani T. Polit, Heather L. Roper, Sara S. Russell, Andrew J. Ryan, Scott A. Sandford, Paul F. Schofield, Cody D. Schultz, et al. (9 additional authors not shown)

73 pages, 22 figures

On 24 September 2023, the NASA OSIRIS-REx mission dropped a capsule to Earth containing approximately 120 g of pristine carbonaceous regolith from Bennu. We describe the delivery and initial allocation of this asteroid sample and introduce its bulk physical, chemical, and mineralogical properties from early analyses. The regolith is very dark overall, with higher-reflectance inclusions and particles interspersed. Particle sizes range from sub-micron dust to a stone about 3.5 cm long. Millimeter-scale and larger stones typically have hummocky or angular morphologies. A subset of the stones appears mottled by brighter material that occurs as veins and crusts. Hummocky stones have the lowest densities and mottled stones have the highest. Remote sensing of the surface of Bennu detected hydrated phyllosilicates, magnetite, organic compounds, carbonates, and scarce anhydrous silicates, all of which the sample confirms. We also find sulfides, presolar grains, and, less expectedly, Na-rich phosphates, as well as other trace phases. The sample composition and mineralogy indicate substantial aqueous alteration and resemble those of Ryugu and the most chemically primitive, low-petrologic-type carbonaceous chondrites. Nevertheless, we find distinct hydrogen, nitrogen, and oxygen isotopic compositions, and some of the material we analyzed is enriched in fluid-mobile elements. Our findings underscore the value of sample return, especially for low-density material that may not readily survive atmospheric entry, and lay the groundwork for more comprehensive analyses.

Hyosun Kim, Mark R. Morris, Jongsoo Kim, Jinhua He

16 pages, 7 figures, 1 table, ApJ in press

We develop a physical framework for interpreting complex circumstellar patterns whorled around asymptotic giant branch (AGB) stars by investigating stable, coplanar triple systems using hydrodynamic and particle simulations. The introduction of a close tertiary body causes an additional periodic variation in the orbital velocity and trajectory of the AGB star. As a result, the circumstellar outflow builds a fine non-Archimedean spiral pattern superimposed upon the Archimedean spiral produced by the outer binary alone. This fine spiral can be approximated by off-centered circular rings that become tangent to each other at the location of the Archimedean spiral. The superimposed fine pattern fades out relatively quickly as a function of distance from the center of the system, in contrast to the dominant Archimedean spiral pattern, which presents a much slower fractional density decrease with radius. The different rates of radial decrease of the density contrast in the two superimposed patterns, coupled with their different time and spatial scales, lead to an apparent, but illusory radial change in the observed pattern interval, as has been reported, for example, in CW Leo. The function describing the detailed radial dependence of the expansion velocity is different in the two patterns, which may be used to distinguish them. The shape of the circumstellar whorled pattern is further explored as a function of the orbital eccentricity and the inner companion's mass. Although this study is confined to stable, coplanar triple systems, the results are likely applicable to moderately noncoplanar systems and open interesting avenues for studying noncoplanar systems.

Alan M. Knee, Jess McIver, Smadar Naoz, Isobel M. Romero-Shaw, Bao-Minh Hoang

10 pages, 5 figures. To be submitted to ApJL

Stellar-mass black hole binaries (BHBs) in galactic nuclei are gravitationally perturbed by the central supermassive black hole (SMBH) of the host galaxy, potentially inducing strong eccentricity oscillations through the eccentric Kozai-Lidov (EKL) mechanism. These highly eccentric binaries emit a train of gravitational-wave (GW) bursts detectable by the Laser Interferometer Space Antenna (LISA) -- a planned space-based GW detector -- with signal-to-noise ratios (SNRs) up to ${\sim}100$ per burst. In this work, we study the GW signature of BHBs orbiting our galaxy's SMBH, Sgr A$^*$, which are consequently driven to very high eccentricities. We demonstrate that an unmodeled approach using a wavelet decomposition of the data effectively yields the time-frequency properties of each burst, provided that the GW frequency peaks between $10^{-3}\,\,\mathrm{Hz}$--$10^{-1}\,\,\mathrm{Hz}$. The wavelet parameters may be used to infer the eccentricity of the binary, measuring $\log_{10}(1-e)$ within an error of $20\%$. Our proposed search method can thus constrain the parameter space to be sampled by complementary Bayesian inference methods, which use waveform templates or orthogonal wavelets to reconstruct and subtract the signal from LISA data.

Bradley Greig, Andrei Mesinger, Eduardo Bañados, George D. Becker, Sarah E. I. Bosman, Huanqing Chen, Frederick B. Davies, Valentina D'Odorico, Anna-Christina Eilers, Simona Gallerani, Martin G. Haehnelt, Laura Keating, Samuel Lai, Yuxiang Qin, Emma Ryan-Weber, Sindhu Satyavolu, Feige Wang, Jinyi Yang, Yongda Zhu

21 pages, 9 figures and 2 tables. Accepted for publication in MNRAS

The attenuation of Ly$\alpha$ photons by neutral hydrogen in the intergalactic medium (IGM) at $z\gtrsim5$ continues to be a powerful probe for studying the epoch of reionisation. Given a framework to estimate the intrinsic (true) Ly$\alpha$ emission of high-$z$ sources, one can infer the ionisation state of the IGM during reionisation. In this work, we use the enlarged XQR-30 sample of 42 high-resolution and high-SNR QSO spectra between $5.8\lesssim\,z\lesssim\,6.6$ obtained with VLT/X-Shooter to place constraints on the IGM neutral fraction. This is achieved using our existing Bayesian QSO reconstruction framework which accounts for uncertainties such as the: (i) posterior distribution of predicted intrinsic Ly$\alpha$ emission profiles (obtained via covariance matrix reconstruction of the Ly$\alpha$ and N V emission lines from unattenuated high-ionisation emission line profiles; C IV, Si IV + O IV] and C III]) and (ii) distribution of ionised regions within the IGM using synthetic damping wing profiles drawn from a $1.6^3$ Gpc$^3$ reionisation simulation. Following careful quality control, we used 23 of the 42 available QSOs to obtain constraints/limits on the IGM neutral fraction during the tail-end of reionisation. Our median and 68th percentile constraints on the IGM neutral fraction are: $0.20\substack{+0.14\\-0.12}$ and $0.29\substack{+0.14\\-0.13}$ at $z = 6.15$~and 6.35. Further, we also report 68th percentile upper-limits of $\bar{x}_{\mathrm{H\,{\scriptscriptstyle I}}} < 0.21$, 0.20, 0.21 and 0.18 at $z = 5.8, 5.95, 6.05$~and 6.55. These results imply reionisation is still ongoing at $5.8\lesssim\,z\lesssim\,6.55$, consistent with previous results from XQR-30 (dark fraction and Ly$\alpha$ forest) along with other observational probes considered in the literature.

Supposing the Laser Interferometer Space Antenna (LISA) gravitational wave (GW) detector, we exhibit the detectability of a hypothetical smooth crossover in the early universe beyond the Standard Model of particle physics through the scalar-induced gravitational wave (SIGW) in terms of the Fisher forecast. A crossover at $\sim100\,\mathrm{TeV}$ can leave a signal on the GW spectrum in the $\sim\mathrm{mHz}$ frequency range, the sweet spot of the LISA sensitivity. These possibilities are also interesting in the primordial black hole (PBH) context as the associated PBH mass $\sim10^{22}\,\mathrm{g}$ lies at the window to explain the whole dark matter. We found that the properties of the crossover can be well determined if the power spectrum of primordial scalar perturbations are as large as $\sim5\times10^{-4}$ on the corresponding scale $\sim10^{12}\,\mathrm{Mpc^{-1}}$.

Adam B. Watts, Luca Cortese, Barbara Catinella, Amelia Fraser-McKelvie, Eric Emsellem, Lodovico Coccato, Jesse van de Sande, Toby H. Brown, Yago Ascasibar, Andrew Battisti, Alessandro Boselli, Timothy A. Davis, Brent Groves, Sabine Thater

16 pages, 7 figures, 1 appendix. Accepted to MNRAS

Stellar feedback-driven outflows are important regulators of the gas-star formation cycle. However, resolving outflow physics requires high resolution observations that can only be achieved in very nearby galaxies, making suitable targets rare. We present the first results from the new VLT/MUSE large program MAUVE (MUSE and ALMA Unveiling the Virgo Environment), which aims to understand the gas-star formation cycle within the context of the Virgo cluster environment. Outflows are a key part of this cycle, and we focus on the peculiar galaxy NGC 4383, which hosts a $\sim6\,$kpc bipolar outflow fuelled by one of Virgo's most HI-rich discs. The spectacular MUSE data reveal the clumpy structure and complex kinematics of the ionised gas in this M82-like outflow at 100 pc resolution. Using the ionised gas geometry and kinematics we constrain the opening half-angle to $\theta=25-35^\circ$, while the average outflow velocity is $\sim210$ kms$^{-1}$. The emission line ratios reveal an ionisation structure where photoionisation is the dominant excitation process. The outflowing gas shows a marginally elevated gas-phase oxygen abundance compared to the disc but is lower than the central starburst, highlighting the contribution of mixing between the ejected and entrained gas. Making some assumptions about the outflow geometry, we estimate an integrated mass outflow-rate of $\sim1.8~$M$_\odot$yr$^{-1}$ and a corresponding mass-loading factor in the range 1.7-2.3. NGC 4383 is a useful addition to the few nearby examples of well-resolved outflows, and will provide a useful baseline for quantifying the role of outflows within the Virgo cluster.

Javier Silva-Farfán, Francisco Förster, Takashi J. Moriya, L. Hernández-García, A. M. Muñoz Arancibia, P. Sanchéz-Sáez, Joseph P. Anderson, John L. Tonry, Alejandro Clocchiatti

We report an analysis of a sample of 186 spectroscopically confirmed Type II supernova (SN) light curves (LCs) obtained from a combination of Zwicky Transient Facility (ZTF) and Asteroid Terrestrial-impact Last Alert System (ATLAS) observations. We implement a method to infer physical parameters from these LCs using hydrodynamic models that take into account the progenitor mass, the explosion energy, and the presence of circumstellar matter (CSM). The CSM is modelled via the mass loss rate, wind acceleration at the surface of the progenitor star with a $\beta$ velocity law, and the CSM radius. We also infer the time of explosion, attenuation (A$_V$), and the redshift for each SN. Our results favor low-mass progenitor stars (M$_{ZAMS}$\,$<$14\,$M_\odot$) with a dense CSM ($\dot{M}$ $>$ 10$^{-3}$ [M$_\odot$ yr$^{-1}$], a CSM radius of $\sim$ 10$^{15}$ cm, and $\beta$ $>$ 2). Additionally, we find that the redshift inferred from the supernova LCs is significantly more accurate than that inferred using the host galaxy photometric redshift, suggesting that this method could be used to infer more accurate host galaxy redshifts from large samples of SNe II in the LSST era. Lastly, we compare our results with similar works from the literature.

Ci Xue, Anthony Remijan, Alexandre Faure, Emmanuel Momjian, Todd R. Hunter, Ryan A. Loomis, Eric Herbst, Brett McGuire

17 figures

At centimeter wavelengths, single-dish observations have suggested that the Sagittarius (Sgr) B2 molecular cloud at the Galactic Center hosts weak maser emission from several organic molecules, including CH$_2$NH, HNCNH, and HCOOCH$_3$. However, the lack of spatial distribution information of these new maser species has prevented us from assessing the excitation conditions of the maser emission as well as their pumping mechanisms. Here, we present a mapping study toward Sgr B2 North (N) to locate the region where the complex maser emission originates. We report the first detection of the Class I methanol (CH$_3$OH) maser at 84 GHz and the first interferometric map of the methanimine (CH$_2$NH) maser at 5.29 GHz toward this region. In addition, we present a tool for modeling and fitting the unsaturated molecular maser signals with non-LTE radiative transfer models and Bayesian analysis using the Markov-Chain Monte Carlo approach. These enable us to quantitatively assess the observed spectral profiles. The results suggest a two-chain-clump model for explaining the intense CH$_3$OH Class I maser emission toward a region with low continuum background radiation. By comparing the spatial origin and extent of maser emission from several molecular species, we find that the 5.29 GHz CH$_2$NH maser has a close spatial relationship with the 84 GHz CH$_3$OH Class I masers. This relationship serves as observational evidence to suggest a similar collisional pumping mechanism for these maser transitions.

The violation of the null energy condition (NEC) may play a crucial role in enabling a scalar field to climb over high potential barriers, potentially significant in the very early universe. We propose a single-field model where the universe sequentially undergoes a first stage of slow-roll inflation, NEC violation, and a second stage of slow-roll inflation. Through the NEC violation, the scalar field climbs over high potential barriers, leaving unique characteristics on the primordial gravitational wave power spectrum, including a blue-tilted nature in the middle-frequency range and diminishing oscillation amplitudes at higher frequencies. Additionally, the power spectrum exhibits nearly scale-invariant behavior on both large and small scales.

M. Dafa Wardana, Masashi Chiba, Kohei Hayashi

19 pages, 14 figures, submitted to ApJ

Distinguishing a core and a cusp within dark matter halos is complexified by the existence of mass-anisotropy degeneracy, where various combinations of velocity anisotropy ($\beta$) and inner density slope ($\gamma$) yield similar observational signatures. We construct a dynamical model that incorporates the 4th-order velocity moments to alleviate this challenge. The inclusion of the 4th-order velocity moments enables stars' line-of-sight velocity distribution (LOSVD) to be flexible. This flexible LOSVD can cover from a thin-tailed to a heavy-tailed distribution that is inaccessible if only the 2nd-order moments are considered. We test our dynamical model using mock galaxies and find that a ratio of the global line-of-sight velocity dispersion of the mock galaxy to the velocity error measurement $\sigma_{\mathrm{los,global}} / \Delta v_{\mathrm{los}} \gtrsim 4$ is required to avoid obtaining systematically biased results. This bias arises from the strong dependency of the 4th-order moments on the LOSVD tails, and not even increasing the sample size to $10^4$ stars can mitigate this effect. In that velocity ratio, $\beta$ is recovered within $\sim 1 \sigma$ even when the sample size is only 500 stars, regardless of the recovery of the other parameters. However, the estimation of $\gamma$ varies, depending on the degree to which LOSVDs deviate from Gaussianity. Because of the more significant change in its LOSVDs, a cored dark halo is more likely to be identified than a cusp.

I. Marini, P. Popesso, G. Lamer, K. Dolag, V. Biffi, S. Vladutescu-Zopp, A. Dev, V. Toptun, E. Bulbul, J. Comparat, N. Malavasi, A. Merloni, T. Mroczkowski, G. Ponti, R. Seppi, S. Shreeram, Y. Zhang

18 pages, 14 figures, submitted to A&A

The eROSITA will deliver an unprecedented volume of X-ray survey observations, 20-30 times more sensitive than ROSAT in the soft band (0.5-2 keV) and for the first time imaging in the hard band (2-10 keV) including galaxy clusters and groups along with obscured and unobscured AGNs. This calls for a powerful theoretical effort to control the systematics and biases that may affect the data analysis. We investigate the detection technique and selection effects in the galaxy group and AGN populations of a mock eROSITA survey at the depth of eRASS:4. We create a $30\times 30$ deg$^{2}$ mock observation based on the cosmological hydrodynamical simulation Magneticum Pathfinder within z=0-0.2. We combine a physical background extracted from the real eFEDS background analysis with realistic simulations of X-ray emission for the hot gas, AGNs and X-ray binaries. We apply a detection procedure equivalent to the reduction done on eRASS data and evaluate the completeness and contamination to reconstruct the luminosity functions of the extended and point sources in the catalogue. We assess the completeness of extended detections as a function of the input X-ray flux and halo. We achieve full recovery of the brightest (most massive) clusters and AGNs. However, a significant fraction of galaxy groups remains undetected. Examining the gas properties between the detected and undetected galaxy groups at fixed halo mass, we observe that the detected population exhibits, on average, higher X-ray brightness compared to the undetected ones. Moreover, we find that X-ray luminosity primarily correlates with the hot gas fraction, rather than temperature or metallicity. Our simulation suggests the presence of a systematic selection effect in current surveys, resulting in X-ray survey catalogues predominantly composed of the lowest-entropy, gas-richest, and highest surface brightness halos on galaxy group scales.

M. Timmermans, G. Dransfield, M. Gillon, A.H.M.J. Triaud, B.V. Rackham, C. Aganze, K. Barkaoui, C. Briceño, A.J. Burgasser, K.A. Collins, M. Cointepas, M. Dévora-Pajares, E. Ducrot, S. Zúñiga-Fernández, S.B. Howell, L. Kaltenegger, C.A. Murray, E.K. Pass, S.N. Quinn, S.N. Raymond, D. Sebastian, K.G. Stassun, C. Ziegler, J.M. Almenara, Z. Benkhaldoun, X. Bonfils, J.L. Christiansen, F. Davoudi, J. de Wit, L. Delrez, B.-O. Demory, W. Fong, G. Fűrész, M. Ghachoui, L.J. Garcia, Y. Gómez Maqueo Chew, M.J. Hooton, K. Horne, M.N. Günther, E. Jehin, J.M. Jenkins, N. Law, A.W. Mann, F. Murgas, F.J. Pozuelos, P.P. Pedersen, D. Queloz, G. Ricker, P. Rowden, R.P. Schwarz, S. Seager, R.L. Smart, G. Srdoc, S. Striegel, S. Thompson, R. Vanderspek, J.N. Winn

27 pages, 19 figures, 7 tables. Accepted for publication in A&A

Small planets transiting bright nearby stars are essential to our understanding of the formation and evolution of exoplanetary systems. However, few constitute prime targets for atmospheric characterization, and even fewer are part of multiple star systems. This work aims to validate TOI-4336 A b, a sub-Neptune-sized exoplanet candidate identified by the TESS space-based transit survey around a nearby M-dwarf. We validate the planetary nature of TOI-4336 A b through the global analysis of TESS and follow-up multi-band high-precision photometric data from ground-based telescopes, medium- and high-resolution spectroscopy of the host star, high-resolution speckle imaging, and archival images. The newly discovered exoplanet TOI-4336 A b has a radius of 2.1$\pm$0.1R$_{\oplus}$. Its host star is an M3.5-dwarf star of mass 0.33$\pm$0.01M$_{\odot}$ and radius 0.33$\pm$0.02R$_{\odot}$ member of a hierarchical triple M-dwarf system 22 pc away from the Sun. The planet's orbital period of 16.3 days places it at the inner edge of the Habitable Zone of its host star, the brightest of the inner binary pair. The parameters of the system make TOI-4336 A b an extremely promising target for the detailed atmospheric characterization of a temperate sub-Neptune by transit transmission spectroscopy with JWST.

Dion Linssen, Jim Shih, Morgan MacLeod, Antonija Oklopčić

Submitted to A&A

Atmospheric escape is thought to significantly influence the evolution of exoplanets, especially for sub-Jupiter planets on short orbital periods. Theoretical models predict that hydrodynamic escape could erode the atmospheres of such gaseous planets, leaving only a rocky core. Deriving atmospheric mass-loss rates from observations is necessary to check these predictions. One of the ways to obtain mass-loss rate estimates is to fit transit spectra of the 10830 {\AA} helium or UV metal lines with Parker wind models. We aim to provide the community with a tool that enables performing this type of analysis, and present sunbather, an open-source Python code to model escaping exoplanet atmospheres and their transit spectra. sunbather incorporates the Parker wind code p-winds and the photoionization code Cloudy, with the ability to calculate any currently known spectral tracer at an arbitrary atmospheric composition. With sunbather, we investigate how the atmospheric structure of a generic hot Neptune planet depends on the metallicity. We find that the mass-loss rate drops by roughly one order of magnitude as we increase the metallicity from solar to 50 times solar. Line cooling by metal species is important already for a solar composition, and more so at higher metallicity. We then demonstrate how sunbather can be used to interpret observations of spectral lines that form in the upper atmosphere. We fit the observed helium spectrum of the mini-Neptune TOI-2134 b and show how even for helium data, the inferred mass-loss rate depends on the metallicity by up to a factor of three.

I review and discuss the possible implications for inflation resulting from considering new physics in light of the Hubble tension. My study is motivated by a simple argument that the constraints on inflationary parameters, most typically the spectral index $n_s$, depend to some extent on the cosmological framework. To avoid broadening the uncertainties resulting from marginalizing over additional parameters (typical in many alternative models), I first adopt the same alternative viewpoint of previous studies and analyze what happens if a physical theory can fix extra parameters to non-standard values. Focusing on the dark energy equation of state $w$ and the effective number of relativistic species $N_{\rm{eff}}$, I confirm that physical theories able to fix $w \approx -1.2$ or $N_{\rm{eff}} \approx 3.9$ produce values of $H_0$ from CMB and BAO in line with the local distance ladder estimate. While in the former case I do not find any relevant implications for inflation, in the latter scenarios, I observe a shift towards $n_s \approx 1$. From a model-selection perspective, both cases are strongly disfavored compared to $\Lambda$CDM. However, models with $N_{\rm{eff}} \approx 3.3 - 3.4$ could bring the $H_0$ tension down to $\sim 3\sigma$ while being moderately disfavored. Yet, this is enough to change the constraints on inflation so that the most accredited models (e.g., Starobinsky inflation) would no longer be favored by data. I then focus on Early Dark Energy (EDE), arguing that an EDE fraction $f_{\rm{EDE}}\sim 0.04 - 0.06$ (only able to mildly reduce the $H_0$-tension down to $\sim 3\sigma$) could already require a similar change in perspective on inflation. In fact, performing a full joint analysis of EDE and Starobinsky inflation, I find that the two models can hardly coexist for $f_{\rm{EDE}}\gtrsim 0.06$.

Large scale, strong magnetic fields are often evoked in black hole accretion flows, for jet launching in the low/hard state and to circumvent the thermal instability in the high/soft state. Here we show how these ideas are strongly challenged by X-ray polarization measurements from IXPE. Quite general arguments show that equipartition fields in the accretion flow should be of order $10^{6-8}$ G. These produce substantial Faraday rotation and/or depolarization for photons escaping the flow in the 2-8 keV IXPE bandpass, which is not consistent with the observed data. While we stress that Faraday rotation should be calculated for each individual simulation (density, field geometry and emissivity), it seems most likely that there are no equipartition strength large scale ordered fields inside the photosphere of the X-ray emitting gas. Strong poloidal fields can still be present if they thread the black hole horizon rather than the X-ray emitting flow, so Blandford-Znajek jets are still possible, but an alternative solution is that the low/hard state jet is dominated by pairs so can be accelerated by lower fields. Fundamentally, polarization data from IXPE means that magnetic fields in black hole accretion flows are no longer invisible and unconstrained.

Callum McCall, Helen Jermak, Iain A. Steele, Iván Agudo, Ulisses Barres de Almeida, Talvikki Hovatta, Gavin P. Lamb, Elina Lindfors, Carole Mundell

32 pages, 35 figures

We present the results of seven years of multicolour photometric monitoring of a sample of 31 $\gamma$-ray bright blazars using the RINGO3 polarimeter on the Liverpool Telescope from 2013--2020. We explore the relationships between simultaneous observations of flux in three optical wavebands along with Fermi $\gamma$-ray data in order to explore the radiation mechanisms and particle populations in blazar jets. We find significant correlations between optical and $\gamma$-ray flux with no detectable time lag, suggesting leptonic emission processes in the jets of these sources. Furthermore, we find the spectral behaviour against optical and $\gamma$-ray flux for many sources is best fit logarithmically. This is suggestive of a transition between bluer-/redder-when-brighter into stable-when-brighter behaviour during high activity states; a behaviour that might be missed in poorly sampled data, resulting in apparent linear relationships.

Mikel F. Barbat, Jürgen Schaffner-Bielich, Laura Tolos

12 pages, 6 figures

We present a comprehensive study of compact stars admixed with non-self annihilating self-interacting fermionic dark matter, delineating the dependence on the nuclear equation of state by considering the two limiting parametrized equations of state for neutron star matter obtained by smoothly matching the low-density chiral effective theory and the high-density perturbative QCD. These two parametrizations are the limiting cases of a wide variety of smooth equations of state, i.e. the softest and stiffest possible one without a phase transition, that generate masses and radii compatible with 2M$_\odot$ observations and the tidal constraint from GW170817. With an exhaustive analysis of the possible stable mass-radius configurations, we determine the quantity of dark matter contained in stars with masses and radii compatible with the aforementioned astrophysical constraints. We find that for dark particle masses of a few tenths of GeV, the dark core collapses and no stable solutions are found irrespective on the chosen nuclear equation of state. For lower masses, the dark matter fraction is limited to 10%, being at most 1% for masses ranging from 0.1 to 10 GeV for the limiting soft nuclear equation of state. For the limiting stiff nuclear equation of state, the dark matter fraction can reach values of more than 10%, but the dark particle mass is being constrained to 0.3 GeV and 10 GeV for the weak self-interacting case and has to be at least 5 GeV for the strong self-interacting one. For dark particle masses of less than 0.1 GeV, stable neutron star configurations should have less than 1% of self-interacting dark matter to be compatible with the constraint of the tidal deformability from GW170817 irrespective on the chosen nuclear equation of state.

Qi Zheng, XueGuang Zhang, QiRong Yuan, Paola Severgnini, Cristian Vignali

9 pages, 5 figures, accepted to be published in MNRAS Letter

In this paper, we presented a detailed analysis of the Sloan Digital Sky Survey optical spectrum of a new sub-kpc scale dual AGN candidate SDSS J222428.53+261423.2 (=SDSS J2224). The target is one of the few AGNs with all the optical narrow emission lines characterized by double-peaked profiles and with peak separations in velocity units of about 930 km/s. If the double-peaked narrow emission lines (DPNELs) are due to a dual AGN in \obj, the estimated physical separation between the two cores is about 500 pc. Meanwhile, three alternative explanations are also discussed in this paper, however, we can not find solid evidence to completely rule them out. Our results support the presence of a sub-kpc dual AGN with DPNELs in all lines, indicating a key episode of galaxy merging evolution at sub-kpc scale.

Ruijing Tang, Dahai Yan, Haiyun Zhang, Qingchang Zhao, Lian Tao, Chengkui Li, Mingyu Ge, Xiaobo Li, Qianqing Yin, Ce Cai

13 pages,17 figures,1 table

Gaussian processes method is employed to analyze the light curves of bursts detected by Insight-HXMT, NICER, and GECAM from SGR 1935+2154 between 2020 to 2022. It is found that a stochastically driven damped simple harmonic oscillator (SHO) is necessary to capture the characteristics of the X-ray bursts. Variability timescale of the X-ray bursts, corresponding to the broken frequencies in the SHO power spectral densities (PSDs), are extracted. In particular, a high broken frequency of 35 Hz where the index of the SHO PSD changes from -4 to -2 is constrained by the HXMT-HE burst associated with FRB 200428. It is suggested that the corresponding timescale of 0.03 s could be the retarding timescale of the system driven by some kind of energy release, and the production of the HE photon should be quasi-simultaneous with the response. The other special event is a NICER burst with a retarding timescale of 1/39 Hz (0.02 s). In the normal X-ray bursts, no retarding timescale is constrained; a long relax/equilibrium timescale (corresponding to a broken frequency of 1-10 Hz where the index of the SHO PSD changing from -4/-2 to 0 in the SHO PSD) is obtained. The results indicate that the FRB-associated HXMT-HE X-ray burst could be produced immediately when the system is responding to the energy disturbance, far before the equilibrium state.

Margarida S. Cunha, Yuri C. Damasceno, Juliana Amaral, Anselmo Falorca, Jørgen Christensen-Dalsgaard, Pedro P. Avelino

16 pages, accepted for publication in Astronomy & Astrophysics

Sharp structural variations induce specific signatures on stellar pulsations that can be studied to infer localised information on the stratification of the star. This information is key to improve our understanding of the physical processes that lead to the structural variations and how to model them. Here we revisit and extend the analysis of the signature of different types of buoyancy glitches in gravity-mode and mixed-mode pulsators presented in earlier works, including glitches with step-like, Gaussian-like, and Dirac-$\delta$-like shapes. In particular, we provide analytical expressions for the perturbations to the periods and show that these can be reliably used in place of the expressions provided for the period spacings, with the advantage that the use of the new expressions does not require modes with consecutive radial orders to be observed. Based on a comparison with two limit cases and on simulated data, we further tested the accuracy of the expression for the Gaussian-like glitch signature whose derivation in an earlier work involved a significant approximation. We find that the least reliable glitch parameter inferred from fitting that expression is the amplitude, which can be up to a factor of two larger than the true amplitude, reaching this limit when the glitch is small. We further discuss the impact on the glitch signature of considering a glitch in the inner and outer half of the g-mode cavity, emphasising the break of symmetry that takes place in the case of mixed-mode pulsators.

Context. Interstellar surface chemistry is a complex process that occurs in icy layers accumulated onto grains of different sizes. Efficiency of surface processes often depends on the immediate environment of adsorbed molecules. Aims. We investigate how gas-grain chemistry changes when surface molecule desorption is made explicitly dependent to the molecular binding energy, which is modified, depending on the properties of the surface. Methods. Molecular binding energy changes gradually for three different environments - bare grain, where polar, water-dominated ices and non-polar, carbon monoxide-dominated ices. In addition to diffusion, evaporation and chemical desorption, photodesorption was also made binding energy-dependent, in line with experimental results. These phenomena occur in a collapsing prestellar core model that considers five grain sizes with ices arranged into four layers. Results. Efficient chemical desorption from bare grains significantly delays ice accumulation. Easier surface diffusion of molecules on non-polar ices promotes the production of carbon dioxide and other species. Conclusions. The composition of interstellar ices is regulated by several binding-energy dependent desorption mechanisms. Their actions overlap in time and space, which explains the ubiquitous proportions of major ice components (water and carbon oxides), observed to be similar in all directions.

Björn J. R. Davidsson

16 pages, 12 figures with 21 panels. This is a pre-copyedited, author-produced PDF of an article accepted for publication in MNRAS following peer review

Cliff collapses on Comet 67P/Churyumov-Gerasimenko expose relatively pristine nucleus matter and offer rare opportunities to characterise ice-rich comet material. Here, Microwave Instrument for \emph{Rosetta} Orbiter (MIRO) observations of two collapsed or crumbling cliffs in the Imhotep and Hathor regions have been assembled. The empirical diurnal antenna temperature curves are analysed with thermophysical and radiative transfer models in order to place constraints on the physical properties and degrees of stratification in the near-surface material. The Imhotep site consists of an exposed dust/water-ice mixture with thermal inertia 100-$160\,\mathrm{J\,m^{-2}\,K^{-1}\,s^{-1/2}}$, having sublimating $\mathrm{CO_2}$ ice located $11\pm 4\,\mathrm{cm}$ below the surface. Its estimated age is consistent with an outburst observed on 2014 April 27-30. The Hathor site has a $0.8\pm 0.2\,\mathrm{cm}$ dust mantle, a thermal inertia of $40\pm 20\,\mathrm{J\,m^{-2}\,K^{-1}\,s^{-1/2}}$, no $\mathrm{CO_2}$ ice to within $1.0\,\mathrm{m}$ depth, and a mantle bulk density of $340\pm 80\,\mathrm{kg\,m^{-3}}$ that is higher than the theoretically expected $180\pm 10\,\mathrm{kg\,m^{-3}}$, suggesting that compression has taken place.

Zhecheng Hu, Wei Zhu, Andrew Gould, Andrzej Udalski, Takahiro Sumi, Ping Chen, Sebastiano Calchi Novati, Jennifer C. Yee, Charles A. Beichman, Geoffery Bryden, Sean Carey, Michael Fausnaugh, B. Scott Gaudi, Calen B. Henderson, Yossi Shvartzvald, Benjamin Wibking, Przemek Mróz, Jan Skowron, Radoslaw Poleski, Michaeł K. Szymański, Igor Soszynśki, Paweł Pietrukowicz, Szymon Kozłowski, Krzysztof Ulaczyk, Krzysztof A. Rybicki, Patryk Iwanek, Marcin Wrona, Mariusz Gromadzki, Fumio Abe, Richard Barry, David P. Bennett, Aparna Bhattacharya, Ian A. Bond, Hirosane Fujii, Akihiko Fukui, Ryusei Hamada, Yuki Hirao, Stela Ishitani Silva, Yoshitaka Itow, Rintaro Kirikawa, Naoki Koshimoto, Yutaka Matsubara, Shota Miyazaki, Yasushi Muraki, Greg Olmschenk, Clément Ranc, et al. (8 additional authors not shown)

13 pages, 10 figures. Submitted to Monthly Notices of the Royal Astronomical Society

We present the analysis of the microlensing event OGLE-2015-BLG-0845, which was affected by both the microlensing parallax and xallarap effects. The former was detected via the simultaneous observations from the ground and Spitzer, and the latter was caused by the orbital motion of the source star in a relatively close binary. The combination of these two effects led to a direct mass measurement of the lens object, revealing a low-mass ($0.14 \pm 0.05 M_{\odot}$) M-dwarf at the bulge distance ($7.6 \pm 1.0$ kpc). The source binary consists of a late F-type subgiant and a K-type dwarf of $\sim1.2 M_{\odot}$ and $\sim 0.9 M_{\odot}$, respectively, and the orbital period is $70 \pm 10$ days. OGLE-2015-BLG-0845 is the first single-lens event in which the lens mass is measured via the binarity of the source. Given the abundance of binary systems as potential microlensing sources, the xallarap effect may not be a rare phenomenon. Our work thus highlights the application of the xallarap effect in the mass determination of microlenses, and the same method can be used to identify isolated dark lenses.

Dori Blakely, Doug Johnstone, Gabriele Cugno, Anand Sivaramakrishnan, Peter Tuthill, Ruobing Dong, Benjamin J. S. Pope, Loïc Albert, Max Charles, Rachel A. Cooper, Matthew De Furio, Louis Desdoigts, René Doyon, Logan Francis, Alexandra Z. Greenbaum, David Lafrenière, James P. Lloyd, Michael R. Meyer, Laurent Pueyo, Shrishmoy Ray, Joel Sánchez-Bermúdez, Anthony Soulain, Deepashri Thatte, Thomas Vandal

Submitted to ApJ

We observed the planet-hosting system PDS 70 with the James Webb Interferometer, JWST's Aperture Masking Interferometric (AMI) mode within NIRISS. Observing with the F480M filter centered at 4.8 $\mu$m, we simultaneously fit a geometric model to the outer disk and the two known planetary companions. We re-detect the protoplanets PDS 70 b and c at an SNR of 21 and 11, respectively. Our photometry of both PDS 70 b and c provide evidence for circumplanetary disk emission through fitting SED models to these new measurements and those found in the literature. We also newly detect emission within the disk gap at an SNR of $\sim$4, at a position angle of $207^{+11}_{-10}$ degrees, and an unconstrained separation within $\sim$200 mas. Follow-up observations will be needed to determine the nature of this emission. We place a 5$\sigma$ upper limit of $\Delta$mag = 7.56 on the contrast of the candidate PDS 70 d at 4.8 $\mu$m, which indicates that if the previously observed emission at shorter wavelengths is due to a planet, this putative planet has a different atmospheric composition than PDS 70 b or c. Finally, we place upper limits on emission from any additional planets in the disk gap. We find an azimuthally averaged 5$\sigma$ upper limit of $\Delta$mag $\approx$ 7.5 at separations greater than 125 mas. These are the deepest limits to date within $\sim$250 mas at 4.8 $\mu$m and the first space-based interferometric observations of this system.

Daniel Marín Pina, Sara Rastello, Mark Gieles, Kyle Kremer, Laura Fitzgerald, Bruno Rando

5 pages, 3 figures, submitted to A&A

Context. The star-black hole (S-BH) binary discovered by the Gaia Collaboration - Gaia BH3 - is chemically and kinematically associated with the metal-poor ED-2 stream in the Milky Way halo. Aims. We explore the possibility that Gaia BH3 was assembled dynamically in the progenitor globular cluster (GC) of the ED-2 stream. Methods. We use a public suite of star-by-star dynamical Monte Carlo models by Kremer et al. (2020) to identify S-BH binaries in GCs with different initial masses and (half-mass) radii. Results. We show that a likely progenitor of the ED-2 stream was a relatively low-mass ($\lesssim10^5M_\odot$) GC with an initial half-mass radius of 2 - 4 pc. Such a GC can dynamically retain a large fraction of its BH population and dissolve on the orbit of ED-2. From the suite of models we find that GCs produce ~ 3 - 30 S-BH binaries, approximately independent of initial GC mass and inversely correlated with initial cluster radius. Scaling the results to the Milky Way GC population, we find that ~$75\%$ of the S-BH binaries formed in GCs are ejected from their host GC in the early phases of evolution ($\lesssim1$ Gyr); these are expected to no longer be close to the stream. The ~25\% of S-BH binaries retained until dissolution are expected to form part of streams, such that for an initial mass of the progenitor of ED-2 of a few $10^4M_\odot$, we expect ~2-3 S-BH to end up in the stream. GC models with metallicities similar to Gaia BH3 ($\lesssim1\%$ solar) include S-BH binaries with similar BH masses ($\gtrsim30M_\odot$), orbital periods, and eccentricities. Conclusions. We predict the Galactic halo contains of order $10^5$ S-BH binaries that formed dynamically in GCs, a fraction of which may readily be detected in Gaia DR4. The detection of these sources provides valuable tests of BH dynamics in clusters and the possible role in formation of gravitational wave sources.

Theoretical models and observational evidence suggest that high-redshift galaxies grow under the bursty mode of star formation, with large temporal star formation rate (SFR) fluctuations around some mean value. From an observational perspective, it has not been clear at which redshift and stellar population characteristics the transition from bursty to smooth star formation occurs. Here, we investigate these using a uniformly reduced sample of NIRSpec prism spectra of 631 galaxies at $3 < z_{\rm spec} < 14$, stacked in 8 redshift and 8 UV slope bins. We evaluate the burstiness of star formation histories using the Balmer break strengths as well as the ratios of SFRs as measured from the emission lines to those measured from the UV continua. The break strength increases monotonically from $z = 10$ to $z = 3$, and from $\beta_{\rm UV} = -3.0$ to $\beta_{\rm UV} = 0.0$. The break strength is tightly anti-correlated with specific SFR (sSFR), and in dusty galaxies, strongly correlated with dust attenuation. Based on the SFR ratios, we find that bursty star formation thrives in the highest redshift, bluest, and lowest stellar mass galaxies, which exhibit the highest sSFRs. The burstiness appears to plateau at $z > 6$, suggesting that we might be observing the peak of star formation burstiness at these redshifts. The $z < 4$ galaxies do not appear particularly bursty, suggesting that the smooth mode of star formation starts taking over right before cosmic noon. As galaxies mature and develop redder UV colors and more pronounced Balmer breaks, their ability to sustain star formation over longer timescales increases, signalling their transition from bursty to smooth star formation.

Kareem El-Badry

7 pages, 4 figures, submitted to OJAp. Comments welcome

A $33\,M_\odot$ black hole (BH) was recently discovered in an 11.6-year binary only 590 pc from the Sun. The system, Gaia BH3, contains a $0.8\,M_\odot$ low-metallicity giant ($\rm [M/H]=-2.2$) and is kinematically part of the Galactic halo, suggesting that the BH formed from a low-metallicity massive star. I show that orbits similar to that of Gaia BH3 are naturally produced through isolated binary evolution. The system's period and eccentricity can result from a broad range of initial orbits with a modest natal kick ($v_{\rm kick}\lesssim 75\,\rm km\,s^{-1}$) to the BH. I construct MESA models for metal-poor massive stars with initial masses ranging from $35-55\,M_{\odot}$, which reach maximum radii of $1150-1800\,R_{\odot}$ as red supergiants. Stars of this size would fit inside most plausible pre-supernova orbits for the system without overflowing their Roche lobes. In addition, models with moderately rapid initial rotation ($\Omega/\Omega_{\rm crit} \gtrsim 0.45$) undergo chemically homogeneous evolution and never expand to radii larger than $10\,R_{\odot}$. There are thus multiple channels through which a low-metallicity, extreme-mass ratio binary could produce a system like Gaia BH3. Dynamical formation scenarios are also viable, and there is little doubt that both isolated and dynamically-formed BH binaries with orbits similar to Gaia BH3 will be discovered in Gaia DR4. Only about 1 in 10,000 stars in the solar neighborhood have metallicities as low as Gaia BH3. This suggests that BH companions are dramatically over-represented at low-metallicity, though caveats related to small number statistics apply. The fact that the luminous star in Gaia BH3 has been a giant - greatly boosting its detectability - only for $\sim$1% of the time since the system's formation implies that additional massive BHs remain to be discovered with only moderately fainter companions.

Francesco Di Clemente, Marco Casolino, Alessandro Drago, Massimiliano Lattanzi, Claudia Ratti

12 pages, 12 figures

Forty years ago, Witten suggested that dark matter could be composed of macroscopic clusters of strange quark matter. This idea was very popular for several years, but it dropped out of fashion once lattice QCD calculations indicated that the confinement/deconfinement transition, at small baryonic chemical potential, is not first order, which seemed to be a crucial requirement in order to produce large clusters of quarks. Here we revisit both the conditions under which strangelets can be produced in the Early Universe and the many phenomenological implications of their existence. Most of the paper discusses the limits on their mass distribution and a possible and simple scheme for their production. Finally, we discuss the most promising techniques to detect this type of objects.

The phenomena of dark matter and the baryon asymmetry pose two of the most pressing questions in today's fundamental physics. Conversion-driven freeze-out has emerged as a successful mechanism to generate the observed dark matter relic density. It supports thermalization of dark matter despite its very weak couplings aligning with the null results from direct and indirect detection experiments. In this letter, we demonstrate that the departure from equilibrium of dark matter, induced by semi-efficient conversions, satisfies Sakharov's conditions, providing a novel explanation for both dark matter and the baryon asymmetry. Specifically, for a leptophilic model, we establish the mechanism of conversion-driven leptogenesis. The scenario predicts dark matter masses close to the electroweak scale offering testable predictions, such as soft displaced leptons at the LHC and future colliders.

In this work, we study the luminosity that results from the conversion of QCD axion particles into photons in the magnetic field of the plasma accreting onto black holes (BHs). For the luminosities to be large two conditions need to be met: i) there are large numbers of axions in the PBH surroundings as a result of the so-called superradiant instability; ii) there exists a point inside the accreting region where the plasma and axion masses are similar and there is resonant axion-photon conversion. For BHs accreting from the interstellar medium in our galaxy, the above conditions require the black hole to have subsolar masses and we are therefore led to consider a population of primordial black holes (PBHs). In the conservative window, where we stay within the non-relativistic behavior of the plasma and neglect the possibility of non-linear enhancement via magnetic stimulation, the typical frequencies of the emitted photons lie on the low-radio band. We thus study the prospects for detection using the LOFAR telescope, assuming the PBH abundance to be close to the maximal allowed by observations. We find that for PBH and QCD axion with masses in the range $10^{-5}-10^{-4}\, M_\odot$ and $4 \times 10^{-8}$ and $4 \times 10^{-7}$ eV, respectively, the flux density emitted by the closest PBH, assuming it accretes from the warm ionized medium, can be detected at the LOFAR telescope. Coincidently, the PBH mass range coincides with the range that would explain the microlensing events found in OGLE. This might further motivate a dedicated search of these signals in the LOFAR data and other radio telescopes.

Axions and axion-like particles are strongly motivated dark matter candidates that are the subject of many current ground based dark matter searches. We present first results from the Axion Dark-Matter Birefringent Cavity (ADBC) experiment, which is an optical bow-tie cavity probing the axion-induced birefringence of electromagnetic waves. Our experiment is the first optical axion detector that is tunable and quantum noise limited, making it sensitive to a wide range of axion masses. We have iteratively probed the axion mass range 40.9-43.3$\text{ neV/c}^2$, 49.3-50.6$\text{ neV/c}^2$, and 54.4-56.7$\text{ neV/c}^2$, and found no dark matter signal. On average, we constrain the ALP-photon coupling at the level $g_{a\gamma\gamma} \leq 1.9\times 10^{-8} \text{ GeV}^{-1}$. We also present prospects for future axion dark matter detection experiments using optical cavities.

F. Agurto-Sepúlveda, J. Oliva, M. Oyarzo, D.R.G Schleicher

14 pages, 8 figures

In this paper we study the qualitative features induces by corrections to GR coming from String Theory, on the shadows of rotating black holes. We deal with the slowly rotating black hole solutions up to order $\mathcal{O}(a^3)$, to first order in $\alpha'$, including also the dilaton. We provide a detailed characterization of the geometry, as well as the ISCO and photon ring, and then we proceed to obtain the black hole images within the relativistic thin-disk model. We characterize the images by computing the diameter, displacement and asymmetry. A comparison with the Kerr case, indicates that all these quantities grow due to the $\alpha'$ correction, and that the departure from GR for different observable is enhanced depending on the angle of view, namely for the diameter the maximum departure is obtained when the system is face-on, while for the displacement and asymmetry the departure from GR is maximized for edge-on point of view.

We study the Barrow cosmological model, which proposes that quantum gravity effects create a complex, fractal structure for the universe's apparent horizon. We leverage the thermodynamics - gravity conjecture. By applying the Clausius relation to the apparent horizon of the Friedmann - Lema\^itre - Robertson - Walker universe within this framework, we derive modified field equations where the Barrow entropy is linked to the horizon. We assess the Barrow cosmology against current observations - cosmic microwave background , supernovae , and baryon acoustic oscillations data - and include projections for future Laser Interferometer Space Antenna (LISA) standard sirens (SS). Our numerical results suggest a modest improvement in the Hubble tension for Barrow cosmology with phantom dark energy behavior, compared to the standard cosmological model. Furthermore, incorporating simulated LISA SS data alongside existing observational constraints tightens the limitations on cosmological parameters, particularly the deformation exponent.

Embedded walls are domain wall solutions which are unstable in the vacuum but stabilized in a plasma of the early Universe. We show how embedded walls in which the electroweak symmetry is restored can lead to an efficient scenario of electroweak baryogenesis. We construct an extension of the Standard Model of particle physics in which embedded walls exist and are stabilized in an electromagnetic plasma.

We analyze the model put forth by Ref. \cite{Oppenheim:2024rcp}, in which it is claimed that "post-quantum classical gravity" (PQCG), a stochastic theory of gravity, leads to modified Newtonian dynamics (MOND) behavior on galactic scales that reproduces galactic rotation curves. We show that this analysis has two basic problems: (i) the equations of PQCG do not lead to a new large scale force of the form claimed in the paper, and (ii) the form claimed is not of the MONDian form anyhow and so does not correspond to observed galactic dynamics. We also mention other potential problems that might arise from stochastic gravitational behavior.