Abstract visibility
Abstract text size

Papers for Friday, Mar 29 2024

Papers with local authors

Valentina La Torre, Anna Sajina, Andy D. Goulding, Danilo Marchesini, Rachel Bezanson, Alan N. Pearl, Laerte Sodré Jr

22 pages, 12 figures, 1 table, accepted for publication in The Astronomical Journal

0 votes
Paper 1 — arXiv:2403.18888
0 votes
Paper 1 — arXiv:2403.18888

The current and upcoming large data volume galaxy surveys require the use of machine learning techniques to maximize their scientific return. This study explores the use of Self-Organizing Maps (SOMs) to estimate galaxy parameters with a focus on handling cases of missing data and providing realistic probability distribution functions for the parameters. We train a SOM with a simulated mass-limited lightcone assuming a ugrizYJHKs+IRAC dataset, mimicking the Hyper Suprime-Cam (HSC) Deep joint dataset. For parameter estimation, we derive SOM likelihood surfaces considering photometric errors to derive total (statistical and systematic) uncertainties. We explore the effects of missing data including which bands are particular critical to the accuracy of the derived parameters. We demonstrate that the parameter recovery is significantly better when the missing bands are "filled-in" rather than if they are completely omitted. We propose a practical method for such recovery of missing data.

M. Tanious, R. Le Gal, R. Neri, A. Faure, A. Gupta, C.J. Law, J. Huang, N. Cuello, J.P. Williams, F. Ménard

16 pages, 12 figures, accepted for publication in Astronomy & Astrophysics

0 votes
Paper 7 — arXiv:2403.18905
0 votes
Paper 7 — arXiv:2403.18905

Over the past few years, chemical studies have revealed multiple structures in the vicinity of young stellar objects (YSOs). It has become evident that specific physical conditions are associated with the emission of particular molecular lines, allowing us to use molecular probes of the YSO physics. Consequently, chemical surveys are now necessary to fully constrain the origin of the observed structures. Several surveys have been conducted to explore the chemistry of YSOs, focusing on Class 0 and Class II objects. However, our knowledge of intermediate objects, that are Class I objects, remains limited. To bridge the gap and establish the relationship between observed structures and molecular line emission at the Class I evolutionary stage, we investigate the spatial distribution of key molecular gas species in the low-mass Class I protostar L1489 IRS (IRAS 04016+2610), a source part of the ChemYSO survey. We performed a 3mm line survey at high spatial and high spectral resolution using the NOEMA interferometer and the IRAM-30m telescope. We present here the ten brightest lines of our survey, in which we identified a new ~ 3 000 au long streamer in HC3N, C2H, and c-C3H2 emission, likely associated with more localized accretion shocks probed in SO. In addition, two ~ 10 000 au bubbles are seen with the dense molecular tracers HCO+, CS, and HCN around the YSO. Additionally, potential indicators of a second outflow appear in CS and HCN emission, but its nature remains to be confirmed. The late infall identified at large scales may originate from the nearby prestellar core L1489 and is likely responsible for the formation of an external warped disk in this system. The detection of a potential second outflow could be the direct evidence of a binary system. Finally, we hypothesize that the bubbles may result from the magnetic pressure as observed in numerical simulations.

J. Richard Bond, George M. Fuller, Evan Grohs, Joel Meyers, Matthew James Wilson

50 pages, 14 figures

0 votes
Paper 21 — arXiv:2403.19038
0 votes
Paper 21 — arXiv:2403.19038

Very different processes characterize the decoupling of neutrinos to form the cosmic neutrino background (C$\nu$B) and the much later decoupling of photons from thermal equilibrium to form the cosmic microwave background (CMB). The C$\nu$B emerges from the fuzzy, energy-dependent neutrinosphere and encodes the physics operating in the early universe in the temperature range $T\sim 10\,{\rm MeV}$ to $T\sim10\,{\rm keV}$. This is the epoch where beyond Standard Model (BSM) physics may be influential in setting the light element abundances and the necessarily distorted fossil neutrino energy spectra. Here we use techniques honed in extensive CMB studies to analyze the C$\nu$B as calculated in detailed neutrino energy transport and nuclear reaction simulations. Our moment method, relative entropy, and differential visibility approach can leverage future high precision CMB and primordial abundance measurements to provide new insights into the C$\nu$B and any BSM physics it encodes. We demonstrate that the evolution of the energy spectrum of the C$\nu$B throughout the weak decoupling epoch is accurately captured in the Standard Model by only three parameters per species, a non-trivial conclusion given the deviation from thermal equilibrium. Furthermore, we can interpret each of the three parameters as physical characteristics of a non-equilibrium system. The success of our compact description within the Standard Model motivates its use also in BSM scenarios. We demonstrate how observations of primordial light element abundances can be used to place constraints on the C$\nu$B energy spectrum, deriving response functions that can be applied for general C$\nu$B spectral distortions. Combined with the description of those deviations that we develop here, our methods provide a convenient and powerful framework to constrain the impact of BSM physics on the C$\nu$B.

Ting-Ting Ge, Xiao-Na Sun, Rui-Zhi Yang, Pak-Hin Thomas Tam, Ming-Xuan Lu, En-Wei Liang

8 pages, 5 figures, accepted for publication in MNRAS. arXiv admin note: text overlap with arXiv:2210.01352

0 votes
Paper 38 — arXiv:2403.19362
0 votes
Paper 38 — arXiv:2403.19362

We report the detection of gamma-ray emission by the Fermi Large Area Telescope (Fermi-LAT) towards the young massive star cluster RCW 38 in the 1-500 GeV photon energy range. We found spatially extended GeV emission towards the direction of RCW 38, which is best modelled by a Gaussian disc of 0.23$\deg$ radius with a significance of the extension is $\sim 11.4 \sigma$. Furthermore, the spatial correlation with the ionized and molecular gas content favors the hadronic origin of the gamma-ray emission. The gamma-ray spectrum of RCW 38 has a relatively hard photon index of $2.44 \pm 0.03$, which is similar to other young massive star clusters. We argue that the diffuse GeV gamma-ray emission in this region likely originates from the interaction of accelerated protons in the stellar cluster with the ambient gas.

All other papers

Samuel H. C. Cabot, Nikku Madhusudhan, Savvas Constantinou, Diana Valencia, Johanna M. Vos, Thomas Masseron, Connor J. Cheverall

Accepted for publication in ApJ Letters

The study of temperate sub-Neptunes is the new frontier in exoplanetary science. A major development in the past year has been the first detection of carbon-bearing molecules in the atmosphere of a temperate sub-Neptune, K2-18 b, a possible Hycean world, with the James Webb Space Telescope (JWST). The JWST is poised to characterise the atmospheres of several other such planets with important implications for planetary processes in the temperate regime. Meanwhile, ground-based high-resolution spectroscopy has been highly successful in detecting chemical signatures of giant exoplanets, though low-mass planets have remained elusive. In the present work, we report the atmospheric reconnaissance of a temperate sub-Neptune using ground-based high-resolution transmission spectroscopy. The long orbital period and the low systemic velocity results in a low planetary radial velocity during transit, making this system a valuable testbed for high-resolution spectroscopy of temperate sub-Neptunes. We observe high-resolution time-series spectroscopy in the H- and K-bands during the planetary transit with the IGRINS instrument (R$\sim$45,000) on Gemini-South. Using observations from a single transit we find marginal evidence (2.2$\sigma$) for the presence of methane (CH$_4$) in the atmosphere and no evidence for ammonia (NH$_3$) despite its strong detectability for a cloud-free H$_2$-rich atmosphere. We assess our findings using injection tests with different atmospheric scenarios, and find them to be consistent with a high CH$_4$/NH$_3$ ratio and/or the presence of high-altitude clouds. Our results demonstrate the capability of Gemini-S/IGRINS for atmospheric characterization of temperate sub-Neptunes, and the complementarity between space- and ground-based facilities in this planetary regime.

Himanish Ganjoo, M. Sten Delos

32 pages, 17 figures. Comments welcome

In cosmologies with hidden sector dark matter, the lightest hidden sector species can come to dominate the energy budget of the universe and cause an early matter-dominated era (EMDE). EMDEs amplify the matter power spectrum on small scales, leading to dense, early-forming microhalos which massively boost the dark matter annihilation signal. We use the Fermi-LAT measurement of the isotropic gamma-ray background to place limits on the parameter space of hidden sector models with EMDEs. We calculate the amplified annihilation signal by sampling the properties of prompt cusps, which reside at the centers of these microhalos and dominate the signal on account of their steep $\rho\propto r^{-3/2}$ density profiles. We also include the portions of the parameter space affected by the gravitational heating that arises from the formation and subsequent destruction of nonlinear structure during the EMDE. We are able to rule out significant portions of the parameter space, particularly at high reheat temperatures. Long EMDEs remain poorly constrained despite large structure-induced boosts to the annihilation signal.

Connor Cheverall, Nikku Madhusudhan

Accepted for publication in AJ

In recent years, high-resolution transmission spectroscopy in the near-infrared has led to detections of prominent molecules in several giant exoplanets on close-in orbits. This approach has traditionally relied on the large Doppler shifts of the planetary spectral lines induced by the high velocities of the close-in planets, which were considered necessary for separating them from the quasi-static stellar and telluric lines. In this work we demonstrate the feasibility of high-resolution transmission spectroscopy for chemical detections in atmospheres of temperate low-mass exoplanets around M dwarfs with low radial velocity variation during transit. We pursue this goal using model injection and recovery tests with H- and K- band high-resolution spectroscopy of the temperate sub-Neptune TOI-732 c, observed using the IGRINS spectrograph on Gemini-S. We show that planetary signals in transit may be recovered when the change in the planet's radial velocity is very small, down to sub-pixel velocities. This is possible due to the presence of the planetary signal in only a subset of the observed spectra. A sufficient number of out-of-transit spectra can create enough contrast between the planet signal and telluric/stellar contaminants that the planet signal does not constitute a principal component of the time-series spectra and can therefore be isolated using PCA-based detrending without relying on a significant Doppler shift. We additionally explore novel metrics for finding such signals, and investigate trends in their detectability. Our work extends the scope of high-resolution transmission spectroscopy and creates a pathway towards the characterisation of habitable sub-Neptune worlds with ground-based facilities.

Massimo Pascale, Brenda L. Frye, Justin D.R. Pierel, Wenlei Chen, Patrick L. Kelly, Seth H. Cohen, Rogier A. Windhorst, Adam G. Riess, Patrick S. Kamieneski, Jose M. Diego, Ashish K. Meena, Sangjun Cha, Masamune Oguri, Adi Zitrin, M. James Jee, Nicholas Foo, Reagen Leimbach, Anton M. Koekemoer, C. J. Conselice, Liang Dai, Ariel Goobar, Matthew R. Siebert, Lou Strolger, S. P. Willner

Submitted to ApJ. 22 pages, 7 Figures

The first James Webb Space Telescope ({\it JWST}) Near InfraRed Camera (NIRCam) imaging in the field of the galaxy cluster PLCK G165.7+67.0 ($z=0.35$) uncovered a Type Ia supernova (SN~Ia) at $z=1.78$, called ``SN H0pe." Three different images of this one SN were detected as a result of strong gravitational lensing, each one traversing a different path in spacetime, thereby inducing a relative delay in the arrival of each image. Follow-up {\it JWST} observations of all three SN images enabled photometric and rare spectroscopic measurements of the two relative time delays. Following strict blinding protocols which oversaw a live unblinding and regulated post-unblinding changes, these two measured time delays were compared to the predictions of seven independently constructed cluster lens models to measure a value for the Hubble constant, $H_0=71.8^{+9.8}_{-7.6}$~km~s$^{-1}$~Mpc$^{-1}$. The range of admissible $H_0$ values predicted across the lens models limits further precision, reflecting the well-known degeneracies between lens model constraints and time delays. It has long been theorized that a way forward is to leverage a standard candle, however this has not been realized until now. For the first time, the lens models are evaluated by their agreement with the SN absolute magnification, breaking these degeneracies and producing our best estimate, $H_0=75.4^{+8.1}_{-5.5}$~km~s$^{-1}$~Mpc$^{-1}$. This is the first precision measurement of $H_0$ from a multiply-imaged SN~Ia, and provides a measurement in a rarely utilized redshift regime. This result agrees with other local universe measurements, yet exceeds the value of $H_0$ derived from the early Universe with $\gtrsim90\%$ confidence, increasing evidence of the Hubble tension. With the precision provided by only four more events, this approach could solidify this disagreement to $>3\sigma$.

The sample of host stars with multiple transiting planets has illuminated the orbital architectures of exoplanetary systems. These architectures may be shaped mostly by formation conditions, be continually sculpted by ongoing dynamical processes, or both. As more studies place planet occurrence within a galactic context, evidence has emerged for variable planet multiplicity over time. In this manuscript, we investigate the use of transit multiplicity as a tool to constrain longer-timescale (>1 Gyr) dynamical sculpting. First, with a suite of injection-and-recovery tests, we quantify sensitivity to sculpting laws across different regimes. We employ a forward modeling framework in which we generate synthetic planetary systems, according to a prescribed sculpting speed and timescale, around the FGK dwarfs studied by the Kepler Mission. Some sculpting scenarios are hypothetically detectable in the Kepler sample, while others can be disfavored from Kepler transit statistics alone. Secondly, we apply our analysis to reverse-engineer the sculpting laws consistent with the true yield from Kepler. We confirm the present-day fraction of host stars containing dynamically cool "systems with tightly-packed inner planets" (STIPs) is 4-13%. A variety of Gyr-timescale sculpting laws successfully predict the transit multiplicity of the Kepler sample, but none of these laws succeeds in also producing a detectable trend with transit multiplicity and stellar age. An improvement to measured stellar age precision may help uncover such a sculpting law, but nevertheless reflects limitations in transit multiplicity as an observable. Therefore other phenomena, apart from Gyr-timescale dynamical sculpting, may be required to explain the Kepler yield.

Davide Lazzati, Rosalba Perna, Taeho Ryu

Submitted to ApJL, comments most welcome

The majority of massive stars are members of binary systems, a fraction of which is expected to remain bound after the first star goes off as a supernova. When the second star also explodes, the SN ejecta are bound to interact with the compact object companion. We explore the consequences of this interaction in the case of a black hole (BH) companion. We show that accretion of the SN ejecta by the BH generally occurs with highly super-Eddington rates, which can be strongly modulated if the ejecta are clumpy, as typically observed in supernova remnants. This late accretion produces transient flaring activity with a time delay from the SN explosion of days to months, depending on the orbital separation and the velocity of the ejecta. The flares are expected to have a non-thermal, broad band spectrum, but their high-frequency emission (UV and X-rays) would be absorbed within the remnant. Observed flares should therefore be in the optical and near-infrared. We propose this model as an explanation to the late-time flaring activity observed in the fast blue optical transient AT2022tsd.

Raphael Galicher, Axel Potier, Johan Mazoyer, Zahed Wahhaj, Pierre Baudoz, Gaël Chauvin

6 pages, 4 figures, 1 table

Imaging exoplanetary systems is essential to characterizing exoplanetary systems and to studying planet-disk interactions to understand planet formation. Such imaging in the visible and near-infrared is challenging because these objects are very faint relative to their star and only fractions of an arcsecond away. Coronagraphic instruments have already allowed the imaging of a few exoplanets, but their performance is limited by wavefront aberrations. Adaptive optics systems partly compensate for the Earth's atmosphere turbulence, but they cannot fully control the wavefront. Some of the starlight leaks through the coronagraph and forms speckles in the image. Focal plane wavefront control, used as a second stage after the adaptive optics system, can minimize the speckle intensity within an area called the dark hole. We demonstrated the on-sky performance of dark hole techniques, pairwise probing coupled with electric field conjugation, using the apodized pupil Lyot coronagraph of the VLT/SPHERE instrument. In this paper, we probe their performance using the SPHERE four-quadrant phase mask coronagraph. We demonstrate the interest of combining dark hole techniques and reference differential imaging (RDI). We create a dark hole on-sky in the narrow band around~$1.7\,\mu$m observing HR\,4796. We then record broadband images of HR\,4796 and a reference star at the H band. The dark hole techniques improve the H-band detection limit by a factor of three. The dark hole is stable from one star to a nearby star enabling RDI. This stability offers two new strategies of observation. First, one can quickly create a dark hole observing a bright star before pointing to a faint target star. Furthermore, one can couple dark hole techniques and RDI. A very interesting point is that the performance of these methods does not depend on the astrophysical signal.

Gabriella Contardo, David W. Hogg

23 pages, 17 figures, 2 tables. Candidates' IDs and (some) code available at this https URL . Comments welcomed!

Stellar infrared excesses can indicate various phenomena of interest, from protoplanetary disks to debris disks, or (more speculatively) techno-signatures along the lines of Dyson spheres. In this paper, we conduct a large search for such excesses, designed as a data-driven contextual anomaly detection pipeline. We focus our search on FGK stars close to the main sequence to favour non-young host stars. We look for excess in the mid-infrared, unlocking a large sample to search in while favouring extreme IR excess akin to the ones produced by Extreme Debris Disks (EDD). We combine observations from ESA Gaia DR3, 2MASS, and the unWISE of NASA WISE, and create a catalogue of 4,898,812 stars with $G < 16$ mag. We consider a star to have an excess if it is substantially brighter in $W1$ and $W2$ bands than what is predicted from an ensemble of machine-learning models trained on the data, taking optical and near-infrared information as input features. We apply a set of additional cuts (derived from the ML models and the objects' astronomical features) to avoid false-positive and identify a set of 53 objects (a rate of $1.1\times 10^{-5}$), including one previously identified EDD candidate. Typical infrared-excess fractional luminosities we find are in the range 0.005 to 0.1, consistent with known EDDs.

Lorenzo Roberti, Marco Limongi, Alessandro Chieffi

39 pages, 6 figures, 4 tables. Accepted for publication in ApJS

We present an extension of the set of models published in Limongi & Chieffi, 2018, ApJS, 237, 13, at metallicity two times solar, i.e. [Fe/H]=0.3. The key physical properties of these models at the onset of the core collapse are mainly due to the higher mass loss triggered by the higher metallicity: the super solar metallicity (SSM) models reach the core collapse with smaller He- and CO-core masses, while the amount of 12C left by the central He burning is higher. These results are valid for all the rotation velocities. The yields of the neutron capture nuclei expressed per unit mass of Oxygen (i.e. the X/O) are higher in the SSM models than in the SM ones in the non rotating case while the opposite occurs in the rotating models. The trend shown by the non rotating models is the expected one, given the secondary nature of the n-capture nucleosynthesis. Vice versa, the counter intuitive trend obtained in the rotating models is the consequence of the higher mass loss present in the SSM models that removes the H rich envelope faster than in the SM ones while the stars are still in central He burning, dumping out the entanglement (activated by the rotation instabilities) and therefore a conspicuous primary n-capture nucleosynthesis.

Bernhard Müller (School of Physics and Astronomy, Monash University)

37 pages, 4 figures. Invited chapter for the edited book "New Frontiers in GRMHD Simulations" (Eds. C. Bambi, Y. Mizuno, S. Shashank and F. Yuan, Springer Singapore, expected in 2024)

Magnetohydrodynamic simulations of core-collapse supernovae have become increasingly mature and important in recent years. Magnetic fields take center stage in scenarios for explaining hypernova explosions, but are now also considered in supernova theory more broadly as an important factor even in neutrino-driven explosions, especially in the context of neutron star birth properties. Here we present an overview of simulation approaches currently used for magnetohydrodynamic supernova simulations and sketch essential physical concepts for understanding the role of magnetic fields in supernovae of slowly or rapidly rotating massive stars. We review progress on simulations of neutrino-driven supernovae, magnetorotational supernovae, and the relevant field amplification processes. Recent results on the nucleosynthesis and gravitational wave emission from magnetorotational supernovae are also discussed. We highlight efforts to provide better initial conditions for magnetohydrodynamic supernova models by simulating short phases of the progenitor evolution in 3D to address uncertainties in the treatment of rotation and magnetic fields in current stellar evolution models.

J.D.R. Pierel, B. L. Frye, M. Pascale, G. B. Caminha, W. Chen, S. Dhawan, D. Gilman, M. Grayling, S. Huber, P. Kelly, S. Thorp, N. Arendse, S. Birrer, M. Bronikowski, R. Canameras, D. Coe, S. H. Cohen, C. J. Conselice, S. P. Driver, J. C. J. Dsilva, M. Engesser, N. Foo, C. Gall, N. Garuda, C. Grillo, N. A. Grogin, J. Henderson, J. Hjorth, R. A. Jansen, J. Johansson, P. S. Kamieneski, A. M. Koekemoer, C. Larison, M. A.Marshall, L. A. Moustakas, M. Nonino, R. Ortiz Iii, T. Petrushevska, N. Pirzkal, A. Robotham, R. E. Ryan, Jr., S. Schuldt, L. G. Strolger, J. Summers, S. H. Suyu, T. Treu, C. N. A. Willmer, R. A. Windhorst, H. Yan, A. Zitrin, A. Acebron, S. Chakrabarti, D. A. Coulter, O. D. Fox, X. Huang, S. W. Jha, G. Li, P. A. Mazzali, A. K. Meena, I. Perez-Fournon, F. Poidevin, A. Rest, A. G. Riess

Accepted in ApJ

Supernova (SN) H0pe is a gravitationally lensed, triply-imaged, Type Ia SN (SN Ia) discovered in James Webb Space Telescope imaging of the PLCK G165.7+67.0 cluster of galaxies. Well-observed multiply-imaged SNe provide a rare opportunity to constrain the Hubble constant ($H_0$), by measuring the relative time delay between the images and modeling the foreground mass distribution. SN H0pe is located at $z=1.783$, and is the first SN Ia with sufficient light curve sampling and long enough time delays for an $H_0$ inference. Here we present photometric time-delay measurements and SN properties of SN H0pe. Using JWST/NIRCam photometry we measure time delays of $\Delta t_{ab}=-116.6^{+10.8}_{-9.3}$ and $\Delta t_{cb}=-48.6^{+3.6}_{-4.0}$ observer-frame days relative to the last image to arrive (image 2b; all uncertainties are $1\sigma$), which corresponds to a $\sim5.6\%$ uncertainty contribution for $H_0$ assuming $70 \rm{km s^{-1} Mpc^{-1}}$. We also constrain the absolute magnification of each image to $\mu_{a}=4.3^{+1.6}_{-1.8}$, $\mu_{b}=7.6^{+3.6}_{-2.6}$, $\mu_{c}=6.4^{+1.6}_{-1.5}$ by comparing the observed peak near-IR magnitude of SN H0pe to the non-lensed population of SNe Ia.

In this paper we have performed a comparative study of different types of oscillating dark energy models using the Metropolis algorithm of MCMC. Eight different oscillating parameterization being examined herein which have demonstrated considerable ability to fit the overall cosmological observations including Pantheon sample of SnIa, Baryon Acoustic Oscillations, Cosmic Chronometer Hubble data and distance priors of Planck CMB. In order to compare the consistency of these models with observations, we have used both Akaike and Deviation information criteria. Although, the values of Akaike information criteria for different models indicate that there is no support for oscillating DE models, Deviation information criteria showed that there is significant support for some of these models. Our results showed that these models are capable of solving cosmic coincidence problem and alleviating the Hubble tension. Comparing the $H_0$ values obtained for different oscillating scenarios, with that of $\Lambda$CDM, we observe that our oscillating models led to $\bar{H_0}=69.78$ which is $0.29$ greater than $H_{0,\Lambda}$ and thus reduce the Hubble tension. Among all of models, Model(1) with $H_0=70.00 \pm 0.71$ is the most capable of alleviating the $H_0$ tension. Furthermore, we examined our models assuming $H_0=73.0 \pm 1.4$ from SHoES measurements. We find that adding this data point to our data combination, led to a $\Delta \bar{H_0}=0.95$ increase in $H_0$ value for different models.

Hanlin Yang, Chunlan Jin, Zifan Wang, Jingxiu Wang

12 pages, 9 figures, 1 table and accepted for publication in the Astrophysical Journal

As a relatively active region, ephemeral region (ER) exhibits highly complex pattern of magnetic flux emergence. We aim to study detailed secondary flux emergences (SFEs) which we define as bipoles that they appear close to ERs and finally coalesce with ERs after a period. We study the SFEs during the whole process from emergence to decay of 5 ERs observed by the Helioseismic and Magnetic Imager (HMI) aboard Solar Dynamics Observatory (SDO) . The maximum unsigned magnetic flux for each ER is around $10^{20}$ Mx. Each ER has tens of SFEs with an average emerging magnetic flux of approximately 5$\times10^{18}$ Mx. The frequency of normalized magnetic flux for all the SFEs follows a power law distribution with an index of -2.08 . The majority of SFEs occur between the positive and negative polarities of ER , and their growth time is concentrated within one hour. The magnetic axis of SFE is found to exhibit a random distribution in the 5 ERs. We suggest that the relationship between SFEs and ERs can be understood by regarding the photospheric magnetic field observations as cross-sections of an emerging magnetic structure. Tracking the ERs' evolution, we propose that these SFEs in ERs may be sequent emergences from the bundle of flux tube of ERs, and that SFEs are partially emerged $\Omega$-loops.

Trey V. Wenger, Daniel R. Rybarczyk, Snežana Stanimirović

13 pages, 3 figures, 3 tables. Accepted for publication in the Astrophysical Journal

The vertical distribution of cold neutral hydrogen (HI) clouds is a constraint on models of the structure, dynamics, and hydrostatic balance of the interstellar medium. In 1978, Crovisier pioneered a method to infer the vertical distribution of HI absorbing clouds in the solar neighborhood. Using data from the Nan\c{c}ay 21-cm absorption survey, they determine the mean vertical displacement of cold HI clouds, $\langle|z|\rangle$. We revisit Crovisier's analysis and explore the consequences of truncating the HI absorption sample in Galactic latitude. For any non-zero latitude limit, we find that the quantity inferred by Crovisier is not the mean vertical displacement but rather a ratio involving higher moments of the vertical distribution. The resultant distribution scale heights are thus ${\sim}1.5$ to ${\sim}3$ times smaller than previously determined. In light of this discovery, we develop a Bayesian Monte Carlo Markov Chain method to infer the vertical distribution of HI absorbing clouds. We fit our model to the original Nan\c{c}ay data and find a vertical distribution moment ratio $\langle|z|^3\rangle/\langle|z|^2\rangle = 97 \pm 15\,\text{pc}$, which corresponds to a Gaussian scale height $\sigma_z = 61 \pm 9\,\text{pc}$, an exponential scale height $\lambda_z = 32 \pm 5\,\text{pc}$, and a rectangular half-width $W_{z, 1/2} = 129 \pm 20\,\text{pc}$. Consistent with recent simulations, the vertical scale height of cold HI clouds appears to remain constant between the inner-Galaxy and the Galactocentric distance of the solar neighborhood. Local fluctuations might explain the large scale height observed at the same Galactocentric distance on the far side of the Galaxy.

Yuto Bekki, Robert H. Cameron, Laurent Gizon

52 pages, 21 figures, 2 tables. Published in Science Advances on 27 Mar 2024. The authors' version

Rapidly rotating fluids have a rotation profile which depends only on the distance from the rotation axis, in accordance with the Taylor-Proudman theorem. Although the Sun was expected to be such a body, helioseismology showed that the rotation rate in the convection zone is closer to constant on radii. It has been postulated that this deviation is due to the poles being warmer than the equator by a few degrees. Using numerical simulations, we show that the pole-to-equator temperature difference cannot exceed 7 Kelvin as a result of the back-reaction of the high-latitude baroclinically unstable inertial modes. The observed amplitudes of the modes further indicate that this maximum temperature difference is reached in the Sun. We conclude that the Sun's latitudinal differential rotation reaches its maximum allowed value.

Varun Chaturmutha, Bernhard Fleck, Stuart Jefferies

12 pages, 4 figures, accepted for publication in ApJ

We present evidence supporting wave reflection in the lower solar chromosphere based on helioseismic analysis of multi-height Doppler data. This evidence is derived through a wave propagation model that incorporates both upward- and downward-traveling (reflected) waves. Moreover, we find that the height of the reflecting region varies with magnetic field strengths in a way that suggests a connection with the the plasma $\beta\sim1$ region. We measure an effective reflection coefficient of $13\,\%$ in a magnetically quiet region of the Sun.

The tidal disruption of old, compact stellar structures provides strong constraints on macroscopic dark matter candidates such as primordial black holes. In view of recent, new observational data on the Eridanus II dwarf galaxy and on its central stellar cluster, we employ, for the first time, $N$-body simulations to assess the impact of compact massive dark matter candidates on the gravitational stability of the cluster. We find evidence that such candidates must be lighter than about one solar mass if they constitute the totality of the dark matter. We additionally derive robust constraints on the fraction of the dark matter in macroscopic objects as a function of mass, by suitably modeling the remainder of the dark matter as standard fluid-like cold dark matter.

Wenlei Chen, Patrick L. Kelly, Brenda L. Frye, Justin Pierel, S. P. Willner, Massimo Pascale, Seth H. Cohen, Christopher J. Conselice, Michael Engesser, Lukas J. Furtak, Daniel Gilman, Norman A. Grogin, Simon Huber, Saurabh W. Jha, Joel Johansson, Anton M. Koekemoer, Conor Larison, Ashish K. Meena, Matthew R. Siebert, Rogier A. Windhorst, Haojing Yan, Adi Zitrin

27 pages (including appendices), 11 figures, 13 supplemental figures

SN H0pe is a triply imaged supernova (SN) at redshift $z=1.78$ discovered using the James Webb Space Telescope (JWST). In order to classify the SN spectroscopically and measure the relative time delays of its three images (designated A, B, and C), we acquired NIRSpec follow-up spectroscopy spanning 0.6 to 5 microns. From the high signal-to-noise spectra of the two bright images B and C, we first classify the SN, whose spectra most closely match those of SN 1994D and SN 2013dy, as a Type Ia SN. We identify prominent blueshifted absorption features corresponding to Si II $\lambda6355$ and Ca II H $\lambda3970$ and K $\lambda3935$. We next measure the absolute phases of the three images from our spectra, which allows us to constrain their relative time delays. The absolute phases of the three images, determined by fitting the three spectra to Hsiao07 SN templates, are $6.5_{-1.8}^{+2.4}$d, $24.3_{-3.9}^{+3.9}$d, and $50.6_{-15.3}^{+16.1}$d for the brightest to faintest images. These correspond to relative time delays between Image A and Image B and between Image B and Image C of $-122.3_{-43.8}^{+43.7}$d and $49.3_{-14.7}^{+12.2}$d, respectively. The SALT3-NIR model yields phases and time delays consistent with these values. After unblinding, we additionally explored the effect of using Hsiao07 template spectra for simulations through eighty instead of sixty days past maximum, and found a small (11.5 and 1.0 days, respectively) yet statistically insignificant ($\sim$0.25$\sigma$ and $\sim$0.1$\sigma$) effect on the inferred image delays.

L. J. A. Scott, C. S. Jeffery, C. M. Byrne, M. Dorsch

11 pages, 11 figures, accepted for publication in MNRAS

Heavy metal subdwarfs are a class of hot subdwarfs with very high abundances of heavy elements, typically around 10 000 times solar. They include stars which are strongly enhanced in either lead or zirconium, as well as other elements. Vertical stratification of the enhanced elements, where the element is concentrated in a thin layer of the atmosphere, has been proposed as a mechanism to explain the apparent high abundances. This paper explores the effects of the vertical stratification of lead on theoretical spectra of hot subdwarfs. The concentration of lead in different regions of the model atmosphere is found to affect individual lines in a broadly wavelength-dependent manner, with the potential for lines to display modified profiles depending on the location of lead enhancement in the atmosphere. This wavelength dependence highlights the importance of observations in both the optical and the UV for determining whether stratification is present in real stars.

The age of the Universe at a given redshift is a fundamental relationship in cosmology. For many years, the uncertainties in it were dauntingly large, close to a factor of 2. In this age of precision cosmology, they are now at the percent level and dominated solely by the Hubble constant. The uncertainties due to the parameters that describe the cosmological model are must less important. In decreasing order they are: uncertainty due to the dark energy equation-of-state parameter $w$, at most 0.9%; uncertainty due to the matter density $\Omega_M$, at most 0.5% and uncertainty due to the curvature parameter $\Omega_k$, at most 0.07%.

Julian Gerber, Jeremy Hare, John A. Tomsick, Benjamin M. Coughenour, Aarran W. Shaw, Maïca Clavel, Francesca Fornasini, Jules Halpern, Alyson Joens, Roman Krivonos, Koji Mukai

Accepted for publication in MNRAS

IGR J18434-0508 is a Galactic Intermediate Polar (IP) type Cataclysmic Variable (CV) previously classified through optical spectroscopy. The source is already known to have a hard Chandra spectrum. In this paper, we have used follow-up XMM-Newton and NuSTAR observations to measure the white dwarf (WD) mass and spin period. We measure a spin period of P = 304.4 +/- 0.3 s based on the combined MOS1, MOS2, and pn light curve. Although this is twice the optical period found previously, we interpret this value to be the true spin period of the WD. The source has an 8 +/- 2% pulsed fraction in the 0.5-10 keV XMM-Newton data and shows strong dips in the soft energy band (0.5-2 keV). The XMM-Newton and NuSTAR joint spectrum is consistent with a thermal bremsstrahlung continuum model with an additional partial covering factor, reflection, and Fe line Gaussian components. Furthermore, we fit the joint spectrum with the post-shock region "ipolar" model which indicates a high WD mass $>$ $\sim$ 1.36 Msun, approaching the Chandrasekhar limit.

Wenzhi Ruan, Rony Keppens, Limei Yan, Patrick Antolin

22 pages, 10 figures

The hour-long, gradual phase of solar flares is well-observed across the electromagnetic spectrum, demonstrating many multi-phase aspects, where cold condensations form within the heated post-flare system, but a complete three-dimensional (3D) model is lacking. Using a state-of-the-art 3D magnetohydrodynamic simulation, we identify the key role played by the Lorentz force through the entire flare lifespan, and show that slow variations in the post-flare magnetic field achieve the bulk of the energy release. Synthetic images in multiple passbands closely match flare observations, and we quantify the role of conductive, radiative and Lorentz force work contributions from flare onset to decay. This highlights how the non-force-free nature of the magnetic topology is crucial to trigger Rayleigh-Taylor dynamics, observed as waving coronal rays in extreme ultraviolet observations. Our C-class solar flare reproduces multi-phase aspects such as post-flare coronal rain. In agreement with observations, we find strands of cooler plasma forming spontaneously by catastrophic cooling, leading to cool plasma draining down the post-flare loops. As there is force balance between magnetic pressure and tension and the plasma pressure in gradual-phase flare loops, this has potential for coronal seismology to decipher the magnetic field strength variation from observations.

Chuan He, Cong Kevin Xu, Ute Lisenfeld, Y Sophia Dai, Taotao Fang, Jia-Sheng Huang, Wei Wang, Qingzheng Yu

Accepted for publication in RAA, ?? pages, 10 figures, 4 tables

We present a study of star-forming galaxies (SFGs) with pseudobulges (bulges with S\'ersic index $\rm n < 2$) in a local close major-merger galaxy pair sample (H-KPAIR). With data from new aperture photometries in the optical and near-infrared bands (aperture size of 7\;kpc) and from the literature, we find that the mean Age of central stellar populations in Spirals with pseudobulges is consistent with that of disky galaxies and is nearly constant against the bulge-to-total ratio (B/T). Paired Spirals have a slightly lower fraction of pure disk galaxies ($\rm B/T \leq 0.1$) than their counterparts in the control sample. Compared to SFGs with classical bulges, those with pseudobulges have a higher ($>2\;\sigma$) mean of specific star formation rate (sSFR) enhancement ($\rm sSFR_{enh} = 0.33\pm0.07$ vs $\rm sSFR_{enh} = 0.12\pm0.06$) and broader scatter (by $\sim 1$\;dex). The eight SFGs that have the highest $\rm sSFR_{enh}$ in the sample all have pseudobulges. A majority (69\%) of paired SFGs with strong enhancement (having sSFR more than 5 times the median of the control galaxies) have pseudobulges. The Spitzer data show that the pseudobulges in these galaxies are tightly linked to nuclear/circum-nuclear starbursts. Pseudobulge SFGs in S+S and in S+E pairs have significantly ($>3\;\sigma$) different sSFR enhancement, with the means of $\rm sSFR_{enh} = 0.45\pm0.08$ and $-0.04\pm0.11$, respectively. We find a decrease in the sSFR enhancements with the density of the environment for SFGs with pseudobulges. Since a high fraction (5/11) of pseudobulge SFGs in S+E pairs are in rich groups/clusters (local density $\rm N_{1Mpc} \geq 7$), the dense environment might be the cause for their low $\rm sSFR_{enh}$.

Jiazheng Dou, Shamik Ghosh, Larissa Santos, Wen Zhao

19 pages, 6 figures, 2 tables. arXiv admin note: text overlap with arXiv:2310.19627

The correlations between T, E modes and B modes in cosmic microwave background (CMB) radiation, which are expected to vanish under parity symmetry, have become a sensitive probe of the new physics beyond the standard model. In this paper, we forecast the estimation of TB and EB cross power spectra using NILC and cILC on AliCPT-1 simulations. We find that, NILC performs better than cILC on TB and EB correlations in light of its lower uncertainties. In terms of the birefringence angle estimation without assuming systematic errors, the combination of CMB TB and EB spectrum from NILC cleaned simulations could reach a sensitivity of $-0.049^\circ<\beta<0.056^\circ$ (95% CL).

Alex Krolewski, Will J. Percival, Alex Woodfinden

8 pages, 2 figures, comments welcome!

We introduce a new method for measuring the Hubble parameter from low-redshift large-scale observations that is independent of the comoving sound horizon. The method uses the baryon-to-photon ratio determined by the primordial deuterium abundance, together with Big Bang Nucleosynthesis (BBN) calculations and the present-day CMB temperature to determine the physical baryon density $\Omega_b h^2$. The baryon fraction $\Omega_b/\Omega_m$ is measured using the relative amplitude of the baryonic signature in galaxy clustering measured by the Baryon Oscillation Spectroscopic Survey, scaling the physical baryon density to the physical matter density. The physical density $\Omega_mh^2$ is then compared with the geometrical density $\Omega_m$ from Alcock-Paczynski measurements from Baryon Acoustic Oscillations (BAO) and voids, to give $H_0$. Including type Ia supernovae and uncalibrated BAO, we measure $H_0 = 67.1^{+6.3}_{-5.3}$ km s$^{-1}$ Mpc$^{-1}$. We find similar results when varying analysis choices, such as measuring the baryon signature from the reconstructed correlation function, or excluding supernovae or voids. This measurement is currently consistent with both the distance-ladder and CMB $H_0$ determinations, but near-future large-scale structure surveys will obtain 3--4$\times$ tighter constraints.

Alex Krolewski, Will J. Percival

39 pages, 20 figures, comments welcome!

The amplitude of the baryon signature in galaxy clustering depends on the cosmological baryon fraction. We consider two ways to isolate this signal in galaxy redshift surveys. First, we extend standard template-based Baryon Acoustic Oscillation (BAO) models to include the amplitude of the baryonic signature by splitting the transfer function into baryon and cold dark matter components with freely varying proportions. Second, we include the amplitude of the split as an extra parameter in Effective Field Theory (EFT) models of the full galaxy clustering signal. We find similar results from both approaches. For the Baryon Oscillation Spectroscopic Survey (BOSS) data we find $f_b\equiv\Omega_b/\Omega_m=0.173\pm0.027$ for template fits post-reconstruction, $f_b=0.153\pm0.029$ for template fits pre-reconstruction, and $f_b=0.154\pm0.022$ for EFT fits, with an estimated systematic error of 0.013 for all three methods. Using reconstruction only produces a marginal improvement for these measurements. Although significantly weaker than constraints on $f_b$ from the Cosmic Microwave Background, these measurements rely on very simple physics and, in particular, are independent of the sound horizon. In a companion paper we show how they can be used, together with Big Bang Nucleosynthesis measurements of the physical baryon density and geometrical measurements of the matter density from the Alcock-Paczynski effect, to constrain the Hubble parameter. While the constraints on $H_0$ based on density measurements from BOSS are relatively weak, measurements from DESI and Euclid will lead to errors on $H_0$ that are competitive with those from local distance ladder measurements.

Sonu Tabitha Paulson, K. K. Mallick, D. K. Ojha

18 pages,16 figures. Accepted for publication in MNRAS

In this paper, we have conducted an investigation focused on a segment of the $Spitzer$ mid-infrared bubble N59, specifically referred to as R1 within our study. Situated in the inner Galactic plane, this region stands out for its hosting of five 6.7 GHz methanol masers, as well as numerous compact H II regions, massive clumps, filaments, and prominent bright rims. As 6.7 GHz masers are closely linked to the initial phases of high-mass star formation, exploring regions that exhibit a high abundance of these maser detections provides an opportunity to investigate relatively young massive star-forming sites. To characterize the R1 region comprehensively, we utilize multi-wavelength (archival) data from optical to radio wavelengths, together with $^{13}$CO and C$^{18}$O data. Utilizing the $Gaia$ DR3 data, we estimate the distance towards the bubble to be $4.66 \pm 0.70$ kpc. By combining near-infrared (NIR) and mid-infrared (MIR) data, we identify 12 Class I and 8 Class II sources within R1. Furthermore, spectral energy distribution (SED) analysis of selected sources reveals the presence of four embedded high-mass sources with masses ranging from 8.70-14.20 M$_\odot$. We also identified several O and B-type stars from radio continuum analysis. Our molecular study uncovers two distinct molecular clouds in the region, which, although spatially close, occupy different regions in velocity space. We also find indications of a potential hub-filament system fostering star formation within the confines of R1. Finally, we propose that the feedback from the H II regions has led to the formation of prominent Bright Rimmed Clouds (BRC) within our region of interest.

O. Feher, S. E. Ragan, F. D. Priestley, P. C. Clark, T. J. T. Moore

18 pages, 9 figures, accepted by MNRAS

Recent advances in identifying giant molecular filaments in galactic surveys allow us to study the interstellar material and its dense, potentially star forming phase on scales comparable to resolved extragalactic clouds. Two large filaments detected in the CHIMPS $^{13}$CO(3-2) survey, one in the Sagittarius-arm and one in an inter-arm region, were mapped with dense gas tracers inside a 0.06 deg$^2$ area and with a spatial resolution of around 0.4 and 0.65 pc at the distance of the targets using the IRAM 30m telescope, to investigate the environmental dependence of the dense gas fraction. The N$_2$H$^+$(1-0) transition, an excellent tracer of the dense gas, was detected in parsec-scale, elliptical clumps and with a filling factor of around 8.5% in our maps. The N$_2$H$^+$-emitting areas appear to have higher dense gas fraction (e.g. the ratio of N$_2$H$^+$ and $^{13}$CO emission) in the inter-arm than in the arm which is opposite to the behaviour found by previous studies, using dust emission rather than N$_2$H$^+$ as a tracer of dense gas. However, the arm filament is brighter in $^{13}$CO and the infrared emission of dust, and the dense gas fraction determined as above is governed by the $^{13}$CO brightness. We caution that measurements regarding the distribution and fraction of dense gas on these scales may be influenced by many scale- and environment-dependent factors, as well as the chemistry and excitation of the particular tracers, then consider several scenarios that can reproduce the observed effect.

Rapidly growing datasets from stellar spectroscopic surveys are providing unprecedented opportunities to analyse the chemical evolution history of our Galaxy. However, spectral analysis requires accurate modelling of synthetic stellar spectra for late-type stars, for which the assumption of local thermodynamic equilibrium (LTE) has been shown to be insufficient in many cases. Errors associated with LTE can be particularly large for Ti I, which is susceptible to over-ionisation, particularly in metal-poor stars. The aims of this work are to study and quantify the 1D non-LTE effects on titanium abundances across the Hertzsprung-Russell diagram for a large sample of stars. A large grid of departure coefficients, $\beta_\nu$, were computed on standard MARCS model atmospheres. The grid extends from 3000K to 8000K in T$_{\mathrm{eff}}$, -0.5 to +5.5 dex in log(g), and -5.0 to +1.0 in [Fe/H], with non-LTE effects in this grid reaching up to 0.4 dex. This was used to compute abundance corrections that were subsequently applied to the LTE abundances of over 70,000 stars selected from the GALAH survey and additional metal-poor dwarfs. The non-LTE effects grow towards lower [Fe/H], lower log(g), and higher T$_{\mathrm{eff}}$, with a minimum and maximum $\Delta$A(Ti I) of 0.02 and 0.19 in the GALAH sample. For metal-poor giants, the non-LTE modelling reduces the average ionisation imbalance from -0.11 dex to -0.01 dex at [Fe/H] = -1.7, and the enhancement in titanium abundances from Ti I lines results in a [Ti/Fe] versus [Fe/H] trend that more closely resembles the behaviour of Ti II at low metallicities. Non-LTE effects on titanium abundances are significant. Neglecting them may alter our understanding of Galactic chemical evolution. We have made our grid of departure coefficients publicly available, with the caveat that the Ti abundances of metal-poor dwarfs need further study in 3D non-LTE.

Wooseok Kang, Ho Seong Hwang, Hyunmi Song, Changbom Park, Narae Hwang, Byeong-Gon Park

21 pages, 17 figures, accepted for publication in ApJS

We study the global kinematics of the Perseus galaxy cluster (Abell 426) at redshift z = 0.017 using a large sample of galaxies from our new MMT/Hectospec spectroscopic observation for this cluster. The sample includes 1447 galaxies with measured redshifts within 60' from the cluster center (1148 from this MMT/Hectospec program and 299 from the literature). The resulting spectroscopic completeness is 67% at r-band apparent magnitude $r_{\rm{Petro, 0}}\leq 18.0$ within 60' from the cluster center. To identify cluster member galaxies in this sample, we develop a new open-source Python package, CausticSNUpy. This code implements the algorithm of the caustic technique and yields 418 member galaxies within 60' of the cluster. We study the cluster using this sample of member galaxies. The cluster shows no significant signal of global rotation. A statistical test shows that the cluster does not have a noticeable substructure within 30'. We find two central regions where the X-ray emitting intracluster medium and galaxies show significant velocity differences ($>7\sigma$). On a large scale, however, the overall morphology and kinematics between the intracluster medium and galaxies agree well. Our results suggest that the Perseus cluster is a relaxed system and has not experienced a recent merger.

Gabriele Mainetti, Fabio Hernandez, Fabrice Jammes, Quentin Le Boulc'h

8 pages, 4 figures, 26th International Conference on Computing in High Energy & Nuclear Physics

The Vera C. Rubin Observatory is preparing for the execution of the most ambitious astronomical survey ever attempted, the Legacy Survey of Space and Time (LSST). Currently in its final phase of construction in the Andes mountains in Chile and due to start operations in 2025 for 10 years, its 8.4-meter telescope will nightly scan the southern sky and collect images of the entire visible sky every 4 nights using a 3.2 Gigapixel camera, the largest imaging device ever built for astronomy. Automated detection and classification of celestial objects will be performed by sophisticated algorithms on high-resolution images to progressively produce an astronomical catalog eventually composed of 20 billion galaxies and 17 billion stars and their associated physical properties. In this paper, we briefly present the infrastructure deployed at the French Rubin data facility (operated by IN2P3 computing center, CC-IN2P3) to deploy the Rubin Science Platform, a set of web-based services to provide effective and convenient access to LSST data for scientific analysis. We describe the main services of the platform, the components that provide those services and our deployment model. We also present the Kubernetes-based infrastructure we are experimenting with for hosting the LSST astronomical catalog, a petabyte-scale relational database developed for the specific needs of the project.

J. M. Coloma-Nadal, F.-S. Kitaura, J. E. García-Farieta, F. Sinigaglia, G. Favole, D. Forero Sánchez

23 pages, 5 figures

Accurate modeling of galaxy distributions is paramount for cosmological analysis using galaxy redshift surveys. However, this endeavor is often hindered by the computational complexity of resolving the dark matter halos that host these galaxies. To address this challenge, we propose the development of effective assembly bias models down to small scales, i.e., going beyond the local density dependence capturing non-local cosmic evolution. We introduce a hierarchical cosmic web classification that indirectly captures up to third-order long- and short-range non-local bias terms. This classification system also enables us to maintain positive definite parametric bias expansions. Specifically, we subdivide the traditional cosmic web classification, which is based on the eigenvalues of the tidal field tensor derived, with an additional classification based on the Hessian matrix of the negative density contrast. We obtain the large-scale dark matter field on a mesh with $\sim3.9\,h^{-1}$ Mpc cell side resolution through Augmented Lagrangian Perturbation Theory. To assess the effectiveness of our model, we conduct tests using a reference halo catalogue extracted from the UNIT project simulation, which was run within a cubical volume of 1 $h^{-1}$ Gpc side. The resulting mock halo catalogs, generated trough our approach, exhibit a high level of accuracy in terms of the one-, two- and three-point statistics. They reproduce the reference power spectrum within better than 2 percent accuracy up to wavenumbers $k\sim0.8\,h$ Mpc$^{-1}$ and provide accurate bispectra within the scales that are crucial for cosmological analysis. This effective bias approach provides a forward model appropriate for field-level cosmological inference and holds significant potential for facilitating cosmological analysis of galaxy redshift surveys, particularly in the context of projects such as DESI, EUCLID, and LSST.

M. Kuźniak, S. Pawłowski, A. Abramowicz, A. F. V. Cortez, M. Kumosiński, T. Łęcki, G. Nieradka

Proceedings of LIght Detection In Noble Elements - LIDINE 2023

Polyethylene naphthalate (PEN) foils have been demonstrated as a wavelength shifter suitable for operation in liquid argon. At the same time, wavelength shifting efficiency of technical grades of PEN, commercially available on the market, is lower than that of tetraphenyl butadiene (TPB). This paper reports on an R&D program focused on exploring the intrinsic limitations of PEN and optimizing it for the highest achievable wavelength shifting efficiency.

Neal Jackson (1), Shruti Badole (1), Thomas Dugdale (1), Hannah R. Stacey (2), Philippa Hartley (3), J.P. McKean (4, 5, 6). ((1) Jodrell Bank Centre for Astrophysics, Department of Physics and Astronomy, University of Manchester, Manchester, UK, (2) European Southern Observatory, Garching bei Muenchen, Germany, (3) SKA Observatory, Jodrell Bank, Macclesfield, UK, (4) Kapteyn Astronomical Institute, University of Groningen, AV Groningen, The Netherlands, (5) South African Radio Astronomy Observatory (SARAO), Krugersdorp, South Africa, (6) Department of Physics, University of Pretoria, Hatfield, Pretoria, South Africa)

Accepted by MNRAS

We present 6-GHz Very Large Array radio images of 70 gravitational lens systems at 300-mas resolution, in which the source is an optically-selected quasar, and nearly all of which have two lensed images. We find that about in half of the systems (40/70, with 33/70 secure), one or more lensed images are detected down to our detection limit of 20microJy/beam, similar to previous investigations and reinforcing the conclusion that typical optically-selected quasars have intrinsic GHz radio flux densities of a few microJy ($\sim10^{23}$WHz$^{-1}$) at redshifts of 1--2. In addition, for ten cases it is likely that the lensing galaxies are detected in the radio. Available detections of, and limits on the far-infrared luminosities from the literature, suggest that nearly all of the sample lie on the radio-FIR correlation typical of star-forming galaxies, and that their radio luminosities are at least compatible with the radio emission being produced by star formation processes. One object, WISE2329$-$1258, has an extra radio component that is not present in optical images, and is difficult to explain using simple lens models. In-band spectral indices, where these can be determined, are generally moderately steep and consistent with synchrotron processes either from star-formation/supernovae or AGN. Comparison of the A/B image flux ratios at radio and optical wavelengths suggests a 10 per cent level contribution from finite source effects or optical extinction to the optical flux ratios, together with sporadic larger discrepancies that are likely to be due to optical microlensing.

Valentin Delabrosse, Catherine Dougados, Sylvie Cabrit, Benoit Tabone, Lukasz Tychoniec, Tom Ray, Linda Podio, Melissa McClure

23 pages, 21 figures, 2 tables, submitted and reviewed to A&A

The origin and impact of outflows on proto-planetary disks and planet formation are key open questions. DG Tau B, a Class I protostar with a structured disk and a striking rotating conical CO outflow, recently identified with ALMA as one of the best MHD disk wind candidate, is an ideal target for studying these phenomena. Our aim is to analyse the outflow components intermediate between the fast axial jet and the wider molecular CO outflow to discriminate between the different scenarios at their origin (irradiated/shocked disk wind or swept-up material). Using observations from JWST NIRSpec-IFU, NIRCam and SINFONI/VLT, we investigate the morphology, kinematics and excitation conditions of H$_2$ emission lines of the red-shifted outflow lobe. We find an onion-like structure of the outflows with increasing temperature, velocity and collimation towards the flow axis. The red-shifted H$_2$ emission reveals a narrow conical cavity nested inside the CO outflow and originating from the inner disk regions (< 6 au). The H$_2$ shell exhibits a constant vertical velocity ($\simeq$22 km/s), twice faster that of the CO flow and an average mass flux of $\dot{M}$(H$_2$) = 3e-11 M$_\odot$/yr significantly lower than the jet and CO values, suggesting low H$_2$ abundance. The global layered structure of the H$_2$/CO outflows is consistent with an MHD disk wind scenario, with the hot H$_2$ possibly tracing an inner dense photodissociation layer of the wind coming from a launching radius in the disk of 0.2-0.4 au. Further analysis, including MIRI observations will provide additional insights into the H$_2$ excitation mechanisms and the origin of the layered outflows observed in DG Tau B.

Thea Hood, Florian Debras, Claire Moutou, Baptiste Klein, Pascal Tremblin, Vivien Parmentier, Andres Carmona, Annabella Meech, Olivia Vénot, Adrien Masson, Pascal Petit, Sandrine Vinatier, Eder Martioli, Flavien Kiefer, Martin Turbet, the ATMOSPHERIX consortium

20 pages, 19 figures, accepted for publication in A&A

Measuring the abundances of C- and O-bearing species in exoplanet atmospheres enables us to constrain the C/O ratio, that contains indications about the planet formation history. With a wavelength coverage going from 0.95 to 2.5 microns, the high-resolution (R$\sim$70 000) spectropolarimeter SPIRou can detect spectral lines of major bearers of C and O in exoplanets. Here we present our study of SPIRou transmission spectra of WASP-76 b acquired for the ATMOSPHERIX program. We applied the publicly available data analysis pipeline developed within the ATMOSPHERIX consortium, analysing the data using 1-D models created with the petitRADTRANS code, with and without a grey cloud deck. We report the detection of H$_2$O and CO at a Doppler shift of around -6 km.s$^{-1}$, consistent with previous observations of the planet. Finding a deep cloud deck to be favoured, we measured in mass mixing ratio (MMR) log(H$_2$O)$_{MMR}$ = -4.52 $\pm$ 0.77 and log(CO)$_{MMR}$ = -3.09 $\pm$ 1.05 consistent with a sub-solar metallicity to more than 1$\sigma$. We report 3$\sigma$ upper limits for the abundances of C$_2$H$_2$, HCN and OH. We estimated a C/O ratio of 0.94 $\pm$ 0.39 ($\sim$ 1.7 $\pm$ 0.7 x solar, with errors indicated corresponding to the 2$\sigma$ values) for the limbs of WASP-76 b at the pressures probed by SPIRou. We used 1-D ATMO forward models to verify the validity of our estimation. Comparing them to our abundance estimations of H$_2$O and CO, as well as our upper limits for C$_2$H$_2$, HCN and OH, we found that our results were consistent with a C/O ratio between 1 and 2 x solar, and hence with our C/O estimation. Finally, we found indications of asymmetry for both H$_2$O and CO when investigating the dynamics of their signatures, pointing to a complex scenario involving possibly both a temperature difference between limbs and clouds being behind the asymmetry this planet is best known for.

Niankun Yu, Zheng Zheng, Chao-Wei Tsai, Pei Zuo, Sara L. Ellison, David V. Stark, Di Li, Jingwen Wu, Karen L. Masters, Ting Xiao, Yinghui Zheng, Zongnan Li, Kai Zhang, Hongying Chen, Shu Liu, Sihan Jiao, Fanyi Meng

Accepted by SCPMA. We sincerely appreciate the constructive suggestions and comments from Lihwai Lin

The atomic-to-molecular gas conversion is a critical step in the baryon cycle of galaxies, which sets the initial conditions for subsequent star formation and influences the multi-phase interstellar medium. We compiled a sample of 94 nearby galaxies with observations of multi-phase gas contents by utilizing public H I, CO, and optical IFU data from the MaNGA survey together with new FAST H I observations. In agreement with previous results, our sample shows that the global molecular-to-atomic gas ratio ($R_{\rm mol} \equiv$ log $M_{\rm H_2}/M_{\rm H\ I}$) is correlated with the global stellar mass surface density $\mu_*$ with a Kendall's $\tau$ coefficient of 0.25 and $p < 10^{-3}$, less tightly but still correlated with stellar mass and NUV$-$ r color, and not related to the specific star formation rate (sSFR). The cold gas distribution and kinematics inferred from the H I and CO global profile asymmetry and shape do not significantly rely on $R_{\rm mol}$. Thanks to the availability of kpc-scale observations of MaNGA, we decompose galaxies into H II, composite, and AGN-dominated regions by using the BPT diagrams. With increasing $R_{\rm mol}$, the fraction of H II regions within 1.5 effective radius decreases slightly; the density distribution in the spatially resolved BPT diagram also changes significantly, suggesting changes in metallicity and ionization states. Galaxies with high $R_{\rm mol}$ tend to have high oxygen abundance, both at one effective radius with a Kendall's $\tau$ coefficient of 0.37 ($p < 10^{-3}$) and their central regions. Among all parameters investigated here, the oxygen abundance at one effective radius has the strongest relation with global $R_{\rm mol}$, but the dependence of gas conversion on gas distribution and galaxy ionization states is weak.

Peijin Zhang, Diana E. Morosan, Pietro Zucca, Sanna Normo, Bartosz Dabrowski, Andrzej Krankowski, Christian Vocks

Context. Observations of solar type II radio bursts provide a unique opportunity to analyze the non-thermal electrons accelerated by coronal shocks and also to diagnose the plasma density distribution in the corona. However, there are very rare high-frequency resolution interferometric observations for type II radio bursts that are capable of tracking these electrons. Aims. Recently, more spatially resolved high-resolution observations of type II radio bursts have been recorded with the Low-Frequency Array (LOFAR). Using these observations, we aim to track the location of a type II radio burst that experiences a sudden spectral bump. Methods. Here, we present the first radio imaging observations for a type II burst with a spectral bump. We measure the variation in source location and frequency drift of the type II burst, and deduct the density distribution along its propagation direction. Results. We identified a type II burst that experiences a sudden spectral bump in its frequency-time profile. The overall frequency drift rate is 0.06 MHz/s and it corresponds to an estimated speed of 295 km/s. The projected speed of the radio source obtained from imaging is 380 km/s towards the east direction. At the spectral bump, a deviation in the source locations of the type II split bands is observed. The band separation increases significantly in the north-south direction. Conclusions. The spectral bump shows an 8 MHz deviation at 60 MHz which corresponds to a 25% decrease in the plasma density. The estimated crossing distance during the spectrum bump in type II is 29 Mm suggesting that this density variation occurs in a confined area. This indicates that the shock most likely encounters the upper extent of a coronal hole.

Zehao Jin, Andrea V. Macciò, Nicholas Faucher, Mario Pasquato, Tobias Buck, Keri L. Dixon, Nikhil Arora, Marvin Blank, Pavle Vulanović

13 pages, 16 figures, published in MNRAS

Cosmological galaxy formation simulations are powerful tools to understand the complex processes that govern the formation and evolution of galaxies. However, evaluating the realism of these simulations remains a challenge. The two common approaches for evaluating galaxy simulations is either through scaling relations based on a few key physical galaxy properties, or through a set of pre-defined morphological parameters based on galaxy images. This paper proposes a novel image-based method for evaluating the quality of galaxy simulations using unsupervised deep learning anomaly detection techniques. By comparing full galaxy images, our approach can identify and quantify discrepancies between simulated and observed galaxies. As a demonstration, we apply this method to SDSS imaging and NIHAO simulations with different physics models, parameters, and resolution. We further compare the metric of our method to scaling relations as well as morphological parameters. We show that anomaly detection is able to capture similarities and differences between real and simulated objects that scaling relations and morphological parameters are unable to cover, thus indeed providing a new point of view to validate and calibrate cosmological simulations against observed data.

Sz. Kálmán, A. Derekas, Sz. Csizmadia, A. Pál, R. Szabó, A. M. S. Smith, K. Nagy, V. Hegedűs, T. Mitnyan, L. Szigeti, Gy. M. Szabó

19 pages, 15 figures, resubmitted to Astronomy and Astrophysics after minor revision

Ultra-hot Jupiters (UHJs) orbiting pulsating A/F stars represent an important subset of the exoplanetary demographic, as they are excellent candidates for the study of exoplanetary atmospheres, as well as being astrophysical laboratories for the investigation of planet-to-star interactions. We analyse the \texttt{TESS} (Transiting Exoplanet Survey Satellite) light curve of the WASP-167 system, consisting of an F1V star and a substellar companion on a $\sim 2.02$ day orbit. We model the combination of the ellipsoidal variability and the Doppler beaming to measure the mass of WASP-167b, and the reflection effect to obtain constraints on the geometric albedo, while placing a special emphasis on noise separation. We implement a basic model to determine the dayside ($T_{\rm Day}$), nightside ($T_{\rm Night}$) and intrinsic ($T_{\rm Internal}$) temperatures of WASP-167b and put a constraint on its Bond albedo. We confirm the transit parameters of the planet seen in the literature. We find that a resonant $\sim 2P^{-1}$ stellar signal (which may originate from planet-to-star interactions) interferes with the phase curve analysis. After considerate treatment of this signal, we find $M_p = 0.34 \pm 0.22$~$M_J$. We measure a dayside temperature of $2790 \pm 100$ K, classifying WASP-167b as an UHJ. We find a $2\sigma$ upper limit of $0.51$ on its Bond albedo, and determine the geometric albedo at $0.34 \pm 0.11$ ($1 \sigma$ uncertainty). With an occultation depth of $106.8 \pm 27.3$ ppm in the \texttt{TESS} passband, the UHJ WASP-167b will be an excellent target for atmospheric studies, especially those at thermal wavelength ranges, where the stellar pulsations are expected to be be less influential.

Prasanta K. Nayak, Mayank Narang, P. Manoj, D. K. Ojha, Blesson Mathew, T. Baug, S. Chandra, S. Vig, G. Maheswar, U. S. Kamath

9 pages, 3 figures, 2 tables

The paper demonstrates the spectroscopic and photometric capabilities of the Ultra-Violet Imaging Telescope (UVIT) to study T-Tauri stars (TTSs). We present the first UVIT/Far-UV spectrum of a TTS, TW Hya. Based on C IV line luminosity, we estimated accretion luminosity (0.1 $L_\odot$) and mass accretion rate (2.2 $\times$ $10^{-8} M_\odot /yr$) of TW Hya, and compared these values with the accretion luminosity (0.03 $L_\odot$) and mass accretion rate (0.6 $\times$ $10^{-8} M_\odot /yr$) derived from spectral energy distribution (SED). From the SED, we derive best-fitted parameters for TW Hya: $T_{eff}$ = 3900$\pm$50 K, radius = 1.2$\pm$0.03 $R_\odot$, $\mathrm{log}\, g = 4.0$ and equivalent black-body temperatures corresponding to accretion luminosity as 14100$\pm$25 K. The parameters of TW Hya derived from UVIT observations were found to be matched well with the literature. Comparison with IUE spectra also suggests that UVIT can be used to study the spectroscopic variability of young stars. This study proposes leveraging the FUV spectroscopic capabilities of UVIT to contribute to the advancement of upcoming UV spectroscopic missions, including the Indian Spectroscopic Imaging Space Telescope (INSIST).

E. Shablovinskaya, C. Ricci, C-S. Chang, A. Tortosa, S. del Palacio, T. Kawamuro, S. Aalto, Z. Arzoumanian, M. Balokovic, F. E. Bauer, K. C. Gendreau, L. C. Ho, D. Kakkad, E. Kara, M. J. Koss, T. Liu, M. Loewenstein, R. Mushotzky, S. Paltani, G. C. Privon, K. Smith, F. Tombesi, B. Trakhtenbrot

Submitted to A&A

The origin of a compact millimeter (mm, 100-250 GHz) emission in radio-quiet active galactic nuclei (RQ AGN) remains debated. Recent studies propose a connection with self-absorbed synchrotron emission from the accretion disk X-ray corona. We present the first joint ALMA ($\sim$100 GHz) and X-ray (NICER/XMM-Newton/Swift; 2-10 keV) observations of the unobscured RQ AGN, IC 4329A ($z = 0.016$). The time-averaged mm-to-X-ray flux ratio aligns with recently established trends for larger samples (Kawamuro et al. 2022, Ricci et al. 2023), but with a tighter scatter ($\sim$0.1 dex) compared to previous studies. However, there is no significant correlation on timescales of less than 20 days. The compact mm emission exhibits a spectral index of $-0.23 \pm 0.18$, remains unresolved with a 13 pc upper limit, and shows no jet signatures. Notably, the mm flux density varies significantly (factor of 3) within 4 days, exceeding the contemporaneous X-ray variability (37% vs. 18%) and showing the largest mm variations ever detected in RQ AGN over daily timescales. The high amplitude variability rules out scenarios of heated dust and thermal free-free emission, pointing toward a synchrotron origin for the mm radiation in a source of $\sim$1 light day size. While the exact source is not yet certain, an X-ray corona scenario emerges as the most plausible compared to a scaled-down jet or outflow-driven shocks.}

Eunseong Lee, Richard Battye, Boris Bolliet

17 pages, 14 figures

We have adapted the Planck cluster likelihood in such a way that it can be applied to the sample of clusters detected by the Atacama Cosmology Telescope (ACT). Applying it to the 2016 sample from Planck and the 2018 sample from ACT we find, by fixing the cosmology using CMB observations and the cluster model adopted by Planck, that the mass bias required by the two are $1-b_{\rm Planck}=0.61\pm 0.03$ and $1-b_{\rm ACT}=0.75\pm 0.06$. These are broadly in agreement but hint that the model could be adapted to reach a better agreement. By normalizing the cluster model using weak lensing observations, we find evidence for either evolution in the cluster model, quantified by the cluster modeling parameter describing redshift dependence $\beta=0.86 \pm 0.07$ using an updated CCCP-based normalization, or evolution in the cosmological model quantified by the dark energy equation of state parameter $w=-0.82 \pm 0.07$.

Gravitational waves (GW) are expected to interact with dark energy and dark matter, affecting their propagation on cosmological scales. In order to model this interaction, we derive a gauge invariant effective equation and action valid for all GWs polarizations, based on encoding the effects of the interaction of GWs at different order in perturbations, in a polarization, frequency and time dependent effective speed. The invariance of perturbations under time dependent conformal transformations and the gauge invariance of the GWs allow to obtain the unitary gauge effective action in any conformally related frame, making transparent the relation between Einstein and Jordan frame. The propagation time and luminosity distance of different GWs polarizations allow to probe at different frequencies and redshift the dark Universe, which act as an effective medium, whose physical properties can be modeled by the GWs effective speed.

Akhil Antony, Fabio Finelli, Dhiraj Kumar Hazra, Daniela Paoletti, Arman Shafieloo

10 pages, 9 figures, 2 tables

We search for super-imposed oscillations linearly or logarithmically spaced in Fourier wavenumbers in Planck and South Pole Telescope (SPT-3G) 2018 temperature and polarization data. The SPT-3G temperature and polarization data provide a new window to test these oscillations at high multipoles beyond the Planck angular resolution and sensitivity. We consider both models with a constant and a Gaussian modulated amplitude, which correspond to three and four additional parameters beyond power-law primordial power spectrum for the templates considered, respectively. We find that each of the four models considered can provide an improved fit to Planck data, consistently with previous findings, and to SPT-3G data, always compared to $\Lambda$CDM. For a constant amplitude of the superimposed oscillations, we find tighter constraints on the amplitude of the oscillations from the combined Planck/SPT-3G data set than in each individual data sets. When the ranges of parameters which provide a better fit to Planck and SPT-3G data overlap, as in the case of Gaussian modulated oscillations, we find a larger $\Delta \chi^2$ - $- 17 \, (-16)$ for logarithmic (linear) oscillations - in a combined Planck/SPT-3G data set than in each individual data sets. These findings will be further tested with upcoming CMB temperature and polarization measurements at high multipoles provided by ongoing ground experiments.

Emanuelly Silva, Ubaldo Zúñiga-Bolaño, Rafael C. Nunes, Eleonora Di Valentino

14 pages, 6 figures and 1 table

Understanding the behavior of the matter power spectrum on non-linear scales beyond the $\Lambda$CDM model is crucial for accurately predicting the large-scale structure (LSS) of the Universe in non-standard cosmologies. In this work, we present an analysis of the non-linear matter power spectrum within the framework of interacting dark energy-dark matter cosmologies (IDE). We employ N-body simulations and theoretical models to investigate the impact of IDE on these non-linear scales. Beginning with N-body simulations characterized by a fixed parameter space delineated by prior observational research, we adeptly fit the simulated spectra with a simple parametric function, achieving accuracy within 5\%. Subsequently, we refine a modified halo model tailored to the IDE cosmology, exhibiting exceptional precision in fitting the simulations down to scales of approximately 1 h/Mpc. To assess the model's robustness, we conduct a forecast analysis for the Euclid survey, employing our refined model. We find that the coupling parameter $\xi$ will be constrained to $\sigma(\xi) = 0.0046$. This marks a significant improvement by an order of magnitude compared to any other current observational tests documented in the literature. These primary findings pave the way for a novel preliminary approach, enabling the utilization of IDE models for observational constraints concerning LSS data on non-linear scales.

Alexios Liakos, Alceste Z. Bonanos, Emmanouil M. Xilouris, Detlef Koschny, Ioannis Bellas-Velidis, Panagiotis Boumis, Athanassios Maroussis, Richard Moissl

32 pages, accepted for publication in A&A

We present results of the NELIOTA campaign for lunar impact flashes observed with the 1.2 m Kryoneri telescope. From August 2019 to August 2023, we report 113 validated and 70 suspected flashes. For the validated flashes, we calculate the physical parameters of the corresponding projectiles, the temperatures developed during the impacts, and the expected crater sizes. For the multiframe flashes we present light curves and thermal evolution plots. Using the whole sample of NELIOTA that encompasses 192 validated flashes in total from 2017, the statistics of the physical parameters of the meteoroids, the peak temperatures of the impacts and the expected crater sizes has been updated. Using this large sample, empirical relations correlating the luminous energies per photometric band were derived and used to roughly estimate the parameters of 92 suspected flashes of the NELIOTA archive. For a typical value of the luminous efficiency, we found that more than the 75% of the impacting meteoroids have masses between 1-200 g, radii between 0.5-3 cm and produced craters up to 3.5 m. 85% of the peak temperatures of the impacts range between 2000 and 4500 K. Statistics regarding the magnitude decline and the cooling rates of the multiframe flashes are also presented. The recalculation of the appearance frequency of meteoroids, lying within the aforementioned ranges of physical parameters, on the Moon yields that the total lunar surface is bombarded with 7.4 sporadic meteoroids/ hour and up to 12.6 meteoroids/hour when the Earth-Moon system passes through a strong meteoroid stream. By extrapolating these rates on Earth, the respective rates for various distances from its surface are calculated and used to estimate the probability of an impact of a meteoroid with a hypothetical infrastructure on the Moon, or with a satellite orbiting Earth for various impact surfaces and duration times of the missions.

Abhijeet Borkar, Romana Grossová, Jiří Svoboda, Emily Moravec, Konstantinos Kouroumpatzakis, Peter G. Boorman, Barbora Adamcová, Beatriz Mingo, Matthias Ehle

9 pages, 4 figures, accepted for publication to A&A

Green Peas (GPs) are young, compact, star-forming dwarf galaxies, and local (z~0.3) analogues of the early galaxies (z>6) considered to be mainly responsible for the reionisation of the Universe. Recent X-ray observations of GPs have detected high excess emission which cannot be accounted for by star formation alone, and implies presence of an active galactic nucleus (AGN). We employ radio observations to study the radio properties of GPs, build their radio spectral energy distributions, and verify the presence of AGNs. We performed new radio observations of three GPs with the Karl G. Jansky Very Large Array (JVLA) in the L, C, and X bands (1.4, 6 and 10 GHz resp.), supplemented by the data from archival observations and large radio surveys. We also analysed the archival radio data for a larger sample of GPs and blueberry (BBs) galaxies, which are lower-mass and lower-redshift analogues of the GPs. To understand the significance of the radio observations, we assess the detectability of these sources, and compare the detected radio luminosities with the expectations from theoretical and empirical relations. Two of the three targeted GPs were strongly detected ($>10\sigma$) in the JVLA observations and their fluxes are consistent with star formation, while the third source was undetected. Although a large fraction (~75%) of the sources from the larger archival sample of GPs and BBs could be detected with archival surveys, only a small number (<40%) are detected and their radio luminosity is significantly lower than the expectation from empirical relations. Our results show that the majority of the dwarf galaxy sample is highly underluminous. Especially towards the lower end of galaxy mass and star formation rate (SFR), the radio luminosity-SFR relation deviates from the empirical relations, suggesting that the relations established for larger galaxies may not hold towards the low-mass end.

John Ellis, Malcolm Fairbairn, Gert Hütsi, Juan Urrutia, Ville Vaskonen, Hardi Veermäe

JWST observations have opened a new chapter in studies of supermassive black holes (SMBHs), stimulating discussion of two puzzles: the abundance of SMBHs in the early Universe and the fraction of dual AGNs. In this paper we argue that the answers to these puzzles may be linked to an interpretation of the data on the nHz gravitational wave (GWs) discovered by NANOGrav and other Pulsar Timing Arrays (PTAs) in terms of SMBH binaries losing energy by interactions with their environments as well as by GW emission. According to this interpretation, the SMBHs in low-$z$ AGNs are the tip of the iceberg of the local SMBH population, which are mainly in inactive galaxies. This interpretation would favour the observability of GW signals from BH binaries in LISA and deciHz GW detectors.

Dominic Agius, Rouven Essig, Daniele Gaggero, Francesca Scarcella, Gregory Suczewski, Mauro Valli

34 pages, 10 figures

If present in the early universe, primordial black holes (PBHs) will accrete matter and emit high-energy photons, altering the statistical properties of the {Cosmic Microwave Background (CMB)}. This mechanism has been used to constrain the fraction of dark matter that is in the form of PBHs to be much smaller than unity for PBH masses well above one solar mass. Moreover, the presence of dense dark matter mini-halos around the PBHs has been used to set even more stringent constraints, as these would boost the accretion rates. In this work, we critically revisit CMB constraints on PBHs taking into account the role of the local ionization of the gas around them. We discuss how the local increase in temperature around PBHs can prevent the dark matter mini-halos from strongly enhancing the accretion process, in some cases significantly weakening previously derived CMB constraints. We explore in detail the key ingredients of the CMB bound and derive a conservative limit on the cosmological abundance of massive PBHs.

The parametrized post-Einsteinian (ppE) framework and its variants are widely used to probe gravity through gravitational-wave tests that apply to a large class of theories beyond general relativity. However, the ppE framework is not truly theory-agnostic as it only captures certain types of deviations from general relativity: those that admit a post-Newtonian series representation in the inspiral of coalescencing compact objects. Moreover, each type of deviation in the ppE framework has to be tested separately, making the whole process computationally inefficient and expensive, possibly obscuring the theoretical interpretation of potential deviations that could be detected in the future. We here present the neural post-Einsteinian (npE) framework, an extension of the ppE formalism that overcomes the above weaknesses using deep-learning neural networks. The core of the npE framework is a variantional autoencoder that maps the discrete ppE theories into a continuous latent space in a well-organized manner. This design enables the npE framework to test many theories simultaneously and to select the theory that best describes the observation in a single parameter estimation run. The smooth extension of the ppE parametrization also allows for more general types of deviations to be searched for with the npE model. We showcase the application of the new npE framework to future tests of general relativity with the fifth observing run of the LIGO-Virgo-KAGRA collaboration. In particular, the npE framework is demonstrated to efficiently explore modifications to general relativity beyond what can be mapped by the ppE framework, including modifications coming from higher-order curvature corrections to the Einstein-Hilbert action at high post-Newtonian order, and dark-photon interactions in possibly hidden sectors of matter that do not admit a post-Newtonian representation.

Several theoretical waveform models have been developed over the years to capture the gravitational wave emission from the dynamical evolution of compact binary systems of neutron stars and black holes. As ground-based detectors improve their sensitivity at low frequencies, the real-time computation of these waveforms can become computationally expensive, exacerbating the steep cost of rapidly reconstructing source parameters using Bayesian methods. This paper describes an efficient numerical algorithm for generating high-fidelity interpolated compact binary waveforms at an arbitrary point in the signal manifold by leveraging computational linear algebra techniques such as singular value decomposition and meshfree approximation. The results are presented for the time-domain \texttt{NRHybSur3dq8} inspiral-merger-ringdown (IMR) waveform model that is fine tuned to numerical relativity simulations and parameterized by the two component-masses and two aligned spins. For demonstration, we target a specific region of the intrinsic parameter space inspired by the previously inferred parameters of the \texttt{GW200311\_115853} event -- a binary black hole system whose merger was recorded by the network of advanced-LIGO and Virgo detectors during the third observation run. We show that the meshfree interpolated waveforms can be evaluated in $\sim 2.3$ ms, which is about $\times 38$ faster than its brute-force (frequency-domain tapered) implementation in the \textsc{PyCBC} software package at a median accuracy of $\sim \mathcal{O}(10^{-5})$. The algorithm is computationally efficient and scales favourably with an increasing number of dimensions of the parameter space. This technique may find use in rapid parameter estimation and source reconstruction studies.

TianQin is a proposed space-based gravitational-wave observatory mission that critically relies on the stability of an equilateral-triangle constellation. Comprising three satellites in high Earth orbits of a $ 10^5 $ km radius, this constellation's geometric configuration is significantly affected by gravitational perturbations, primarily originating from the Moon and the Sun. In this paper, we present an analytical model to quantify the effects of lunisolar perturbations on the TianQin constellation, derived using Lagrange's planetary equations. The model provides expressions for three kinematic indicators of the constellation: arm-lengths, relative line-of-sight velocities, and breathing angles. Analysis of these indicators reveals that lunisolar perturbations can distort the constellation triangle, resulting in three distinct variations: linear drift, bias, and fluctuation. Furthermore, it is shown that these distortions can be optimized to display solely fluctuating behavior, under certain predefined conditions. These results can serve as the theoretical foundation for numerical simulations and offer insights for engineering a stable constellation in the future.

Luis A. Anchordoqui, Ignatios Antoniadis, Dieter Lust

32 pages revtex with 3 figures

In the last two years the dark dimension scenario has emerged as focal point of many research interests. In particular, it functions as a stepping stone to address the cosmological hierarchy problem and provides a colosseum for dark matter contenders. We reexamine the possibility that primordial black holes (PBHs) perceiving the dark dimension could constitute all of the dark matter in the universe. We re-assess limits on the abundance of PBHs as dark matter candidates from $\gamma$-ray emission resulting from Hawking evaporation. We re-evaluate constraints from the diffuse $\gamma$-ray emission in the direction of the Galactic center which offer the best and most solid upper limits on the dark matter fraction composed of PBHs. The revised mass range which allows PBHs to assemble all cosmological dark matter is estimated to be $10^{15} \alt M_{\rm BH}/{\rm g} \alt 10^{21}$. We demonstrate that due to the constraints from $\gamma$-ray emission, quantum corrections due to the speculative memory burden effect do not modify this mass range. We also investigate the main characteristics of PBHs which are localized in the bulk. We show that PBHs localized in the bulk can make all cosmological dark matter if $10^{11} \alt M_{\rm BH}/{\rm g} \alt 10^{21}$. Finally, we comment on the black holes that could be produced if one advocates a space with two boundaries for the dark dimension.