Abstract visibility
Abstract text size

Papers for Wednesday, Jul 26 2023

This study focuses on attitude and control motion of two bodies (a base-satellite and a sub-satellite) connected by an inextensible and massless tether in a circular orbit under the influence of the Earths gravitational force. The base-satellite is assumed to be far more heavier than the sub-satellite. In such cases, the base-satellite is regarded as the reference spacecraft. Because of the complexity of the problem, no thrusters on the sub-satellite are considered, and the effect of atmospheric drag, Earths oblateness, and electrodynamic force on the spacecraft are neglected.

Lucas Einig (IRAM), Jérôme Pety (IRAM, LERMA (UMR\_8112)), Antoine Roueff (IM2NP), Paul Vandame, Jocelyn Chanussot, Maryvonne Gerin, Jan H. Orkisz, Pierre Palud (CRIStAL, LERMA), Miriam Garcia Santa-Maria, Victor de Souza Magalhaes, Ivana Bešlić, Sébastien Bardeau, Emeric E. Bron (ICMM), Pierre Chainais (CRIStAL), Javier R Goicoechea, Pierre Gratier, Viviana Guzman Veloso, Annie Hughes, Jouni Kainulainen, David Languignon, Rosine Lallement, François Levrier, Dariuscz C. Lis, Harvey Liszt, Jacques Le Bourlot, Franck Le Petit, Karin Danielsson Öberg, Nicolas Peretto, Evelyne Roueff, Albrecht Sievers, Pierre-Antoine Thouvenin (CRIStAL), Pascal Tremblin (MDLS)

Context. The availability of large bandwidth receivers for millimeter radio telescopes allows the acquisition of position-position-frequency data cubes over a wide field of view and a broad frequency coverage. These cubes contain much information on the physical, chemical, and kinematical properties of the emitting gas. However, their large size coupled with inhomogenous signal-to-noise ratio (SNR) are major challenges for consistent analysis and interpretation.Aims. We search for a denoising method of the low SNR regions of the studied data cubes that would allow to recover the low SNR emission without distorting the signals with high SNR.Methods. We perform an in-depth data analysis of the 13 CO and C 17 O (1 -- 0) data cubes obtained as part of the ORION-B large program performed at the IRAM 30m telescope. We analyse the statistical properties of the noise and the evolution of the correlation of the signal in a given frequency channel with that of the adjacent channels. This allows us to propose significant improvements of typical autoassociative neural networks, often used to denoise hyperspectral Earth remote sensing data. Applying this method to the 13 CO (1 -- 0) cube, we compare the denoised data with those derived with the multiple Gaussian fitting algorithm ROHSA, considered as the state of the art procedure for data line cubes.Results. The nature of astronomical spectral data cubes is distinct from that of the hyperspectral data usually studied in the Earth remote sensing literature because the observed intensities become statistically independent beyond a short channel separation. This lack of redundancy in data has led us to adapt the method, notably by taking into account the sparsity of the signal along the spectral axis. The application of the proposed algorithm leads to an increase of the SNR in voxels with weak signal, while preserving the spectral shape of the data in high SNR voxels.Conclusions. The proposed algorithm that combines a detailed analysis of the noise statistics with an innovative autoencoder architecture is a promising path to denoise radio-astronomy line data cubes. In the future, exploring whether a better use of the spatial correlations of the noise may further improve the denoising performances seems a promising avenue. In addition,

A. Pandhi, R. K. Friesen, L. Fissel, J. E. Pineda, P. Caselli, M. C-Y. Chen, J. Di Francesco, A. Ginsburg, H. Kirk, P. C. Myers, S. S. R. Offner, A. Punanova, F. Quan, E. Redaelli, E. Rosolowsky, S. Scibelli, Y. M. Seo, Y. Shirley

33 pages, 28 figures, accepted to MNRAS

Studies of dense core morphologies and their orientations with respect to gas flows and the local magnetic field have been limited to only a small sample of cores with spectroscopic data. Leveraging the Green Bank Ammonia Survey alongside existing sub-millimeter continuum observations and Planck dust polarization, we produce a cross-matched catalogue of 399 dense cores with estimates of core morphology, size, mass, specific angular momentum, and magnetic field orientation. Of the 399 cores, 329 exhibit 2D $\mathrm{v}_\mathrm{LSR}$ maps that are well fit with a linear gradient, consistent with rotation projected on the sky. We find a best-fit specific angular momentum and core size relationship of $J/M \propto R^{1.82 \pm 0.10}$, suggesting that core velocity gradients originate from a combination of solid body rotation and turbulent motions. Most cores have no preferred orientation between the axis of core elongation, velocity gradient direction, and the ambient magnetic field orientation, favouring a triaxial and weakly magnetized origin. We find, however, strong evidence for a preferred anti-alignment between the core elongation axis and magnetic field for protostellar cores, revealing a change in orientation from starless and prestellar populations that may result from gravitational contraction in a magnetically-regulated (but not dominant) environment. We also find marginal evidence for anti-alignment between the core velocity gradient and magnetic field orientation in the L1228 and L1251 regions of Cepheus, suggesting a preferred orientation with respect to magnetic fields may be more prevalent in regions with locally ordered fields.

M. Sten Delos, Michael Korsmeier, Axel Widmark, Carlos Blanco, Tim Linden, Simon D. M. White

11 pages, 10 figures

Recent studies indicate that thermally produced dark matter will form highly concentrated, low-mass cusps in the early universe that often survive until the present. While these cusps contain a small fraction of the dark matter, their high density significantly increases the expected gamma-ray flux from dark matter annihilation, particularly in searches of large angular regions. We utilize 14 years of Fermi-LAT data to set strong constraints on dark matter annihilation through a detailed study of the isotropic gamma-ray background, excluding with 95% confidence dark matter annihilation to $b\bar{b}$ final states for dark matter masses below 120 GeV.

Xavier Rodrigues, Vaidehi S. Paliya, Simone Garrappa, Anastasiia Omeliukh, Anna Franckowiak, Walter Winter

Submitted to A&A. Contains 28 pages, 13 figures. Results available online in machine-readable format (see caption of Tab. B.1.)

The origin of the diffuse astrophysical neutrino flux observed by the IceCube experiment is still under debate. In recent years there have been associations of neutrino events with individual blazars, which are active galaxies with relativistic jets pointing toward Earth, such as the source TXS 0506+056. From a theoretical perspective, the properties of these sources as neutrino emitters are not yet well understood. In this work we model a sample of 324 blazars detected by the Fermi Large Area Telescope (LAT), most of which are flat-spectrum radio quasars (FSRQs). This amounts to 34% of all FSRQs in the latest Fermi catalog. By numerically modelling the interactions of cosmic-ray electrons and protons, we explain the emitted multi-wavelength fluxes from each source and self-consistently predict the emitted neutrino spectrum. We demonstrate that the optical and GeV gamma-ray broadband features are generally well described by electron emission. For 33% of the blazars in our sample, a description of the observed X-ray spectrum benefits from an additional component from proton interactions, in agreement with recent studies of individual IceCube candidate blazars. We conclude that blazars that are brighter in GeV gamma rays tend to have a higher neutrino production efficiency but a lower best-fit baryonic loading. The predicted neutrino luminosity shows a positive correlation with the observed GeV gamma-ray flux and with the predicted MeV gamma-ray flux. By extrapolating the results for this sample, we show that the diffuse neutrino flux from the population of gamma-ray-bright blazars may be at the level of about 20% of the IceCube flux, in agreement with current limits from stacking analyses. We discuss the implications of our results for future neutrino searches and suggest promising sources for potential detections with future experiments.

Janus Kozdon, Sean Brittain, Jeffrey Fung, Josh Kern, Stanley Jensen, John Carr, Joan Najita, Andrea Banzatti

Accepted to AJ

CI Tau is currently the only T Tauri star with an inner protoplanetary disk that hosts a planet, CI Tau b, that has been detected by a radial velocity survey. This provides the unique opportunity to study disk features that were imprinted by that planet. We present multi-epoch spectroscopic data, taken with NASA IRTF in 2022, of the ${}^{12}$CO and hydrogen Pf$\beta$ line emissions spanning 9 consecutive nights, which is the proposed orbital period of CI Tau b. We find that the star's accretion rate varied according to that 9~d period, indicative of companion driven accretion. Analysis of the ${}^{12}$CO emission lines reveals that the disk can be described with an inner and outer component spanning orbital radii 0.05-0.13~au and 0.15-1.5~au, respectively. Both components have eccentricities of about 0.05 and arguments of periapses that are oppositely aligned. We present a proof-of-concept hydrodynamic simulation that shows a massive companion on a similarly eccentric orbit can recreate a similar disk structure. Our results allude to such a companion being located around an orbital distance of 0.14~au. However, this planet's orbital parameters may be inconsistent with those of CI Tau b whose high eccentricity is likely not compatible with the low disk eccentricities inferred by our model.

Romy Rodríguez Martínez, David V. Martin, B. Scott Gaudi, Joseph G. Schulze, Anusha Pai Asnodkar, Kiersten M. Boley, Sarah Ballard

20 pages, 11 figures, 2 tables. Submitted to ApJ and under review. Comments welcome!

We investigate and compare the composition of M-dwarf planets in systems with only one known planet (``singles") to those residing in multi-planet systems (``multis") and the fundamental properties of their host stars. We restrict our analysis to planets with directly measured masses and radii, which comprise a total of 70 planets: 30 singles and 40 multis in 19 systems. We compare the bulk densities for the full sample, which includes planets ranging in size from $0.52 R_{\oplus}$ to $12.8R_\oplus$, and find that single planets have significantly lower densities on average than multis, which we cannot attribute to selection biases. We compare the bulk densities normalized by an Earth model for planets with $R_{p} < 6R_{\oplus}$, and find that multis are also denser with 99\% confidence. We calculate and compare the core/water mass fractions (CMF/WMF) of low-mass planets ($M_p <10 M_{\oplus}$), and find that the likely rocky multis (with $R_p <1.6 R_{\oplus}$) have lower CMFs than singles. We also compare the [Fe/H] metallicity and rotation period of all single versus multi-planet host stars with such measurements in the literature and find that multi-planet hosts are significantly more metal-poor than those hosting a single planet. Moreover, we find that host star metallicity decreases with increasing planet multiplicity. In contrast, we find only a modest difference in the rotation period. The significant differences in planetary composition and metallicity of the host stars point to different physical processes governing the formation of single- and multi-planet systems in M dwarfs.

Mariana P. Júlio, Jarle Brinchmann, Sebastiaan L. Zoutendijk, Justin I. Read, Daniel Vaz, Sebastian Kamann, Davor Krajnović, Leindert A. Boogaard, Matthias Steinmetz, Nicolas Bouché

Accepted for publication in A&A on 21-07-2023. 17 pages, 22 figures, 4 tables

Aims. We use stellar line-of-sight velocities of Antlia B, a faint dwarf galaxy in the NGC 3109 association, to derive constraints on the fundamental properties of scalar field dark matter originally proposed to solve the small-scale problems faced by cold dark matter models. Methods. We use the first spectroscopic observations of Antlia B, a distant (d $\sim$ 1.35 Mpc) faint dwarf ($M_\text{V} = -9.7$, M$_\star \sim 8\times10^5$M$_\odot$), from MUSE-Faint - a survey of ultra-faint dwarfs with the Multi Unit Spectroscopic Explorer. Through measurement of line-of-sight velocities for stars in the $1'\times 1'$ field-of-view, we identify 127 stars as members of Antlia B, allowing us to model its dark matter density profile with the Jeans modelling code GravSphere. We implement a model for scalar field dark matter into GravSphere and use this to place constraints on the self-coupling strength of this model. Results. We find a virial mass of ${M_{200} \approx 1.66^{+2.51}_{-0.92}\times 10^9}$ M$_\odot$ and a concentration parameter of ${c_{200}\approx 17.38^{+6.06}_{-4.20}}$ for Antlia B. These results are consistent with the mass-concentration relations in the literature. We constrain the characteristic length scale of the repulsive self-interaction $R_{\text{TF}}$ of the scalar field dark matter model to $R_{\text{TF}} \lesssim 180$ pc (68% confidence level), which translates to a self-coupling strength of $\frac{g}{m^2c^4}\lesssim 5.2 \times 10^{-20}$ eV$^{-1}$cm$^3$. The constraint on the characteristic length scale of the repulsive self-interaction is inconsistent with the value required to match the observations of cores of dwarf galaxies in the Local Group, suggesting that the cored density profiles of those galaxies are not caused by scalar field dark matter.

Lumen Boco (1), Andrea Lapi (1), Francesco Shankar (2), Hao Fu (2), Francesco Gabrielli (1), Alex Sicilia (1) (1 - SISSA, 2 - University of Southampton)

22 pages, 13 Figures

In recent years, increasing attention has been devoted to semi empirical, data-driven models to tackle some aspects of the complex and still largely debated topic of galaxy formation and evolution. We here present a new semi empirical model whose marking feature is simplicity: it relies on solely two assumptions, one initial condition and two free parameters. Galaxies are connected to evolving dark matter haloes through abundance matching between specific halo accretion rate (sHAR) and specific star formation rate (sSFR). Quenching is treated separately, in a fully empirical way, to marginalize over quiescent galaxies and test our assumption on the sSFR evolution without contaminations from passive objects. Our flexible and transparent model is able to reproduce the observed stellar mass functions up to $z\sim 5$, giving support to our hypothesis of a monotonic relation between sHAR and sSFR. We then exploit the model to test a hypothesis on morphological evolution of galaxies. We attempt to explain the bulge/disk bimodality in terms of the two halo accretion modes: fast and slow accretion. Specifically, we speculate that bulge/spheroidal components might form during the early phase of fast halo growth, while disks form during the later phase of slow accretion. We find excellent agreement with both the observational bulge and elliptical mass functions.

Andrés N. Ruiz, Héctor J. Martínez, Valeria Coenda, Hernán Muriel, Sofía A. Cora, Martín de los Rios, Cristian A. Vega-Martínez

Accepted for publication in MNRAS. Comments are welcome

We study the population of backsplash galaxies at $z=0$ in the outskirts of massive, isolated clusters of galaxies taken from the MDPL2-SAG semi-analytic catalogue. We consider four types of backsplash galaxies according to whether they are forming stars or passive at three stagesin their lifetimes: before entering the cluster, during their first incursion through the cluster, and after they exit the cluster. We analyse several geometric, dynamic, and astrophysical aspects of the four types at the three stages. Galaxies that form stars at all stages account for the majority of the backsplash population ($58\%$) and have stellar masses typically below $M_\star\sim 3\times 10^{10} h^{-1}{\rm M}_\odot$ that avoid the innermost cluster's regions and are only mildly affected by it. In a similar mass range, galaxies that become passive after exiting the cluster ($26\%$) follow orbits characterised by small pericentric distance and a strong deflection by the cluster potential well while suffering a strong loss of both dark matter and gas content. Only a small fraction of our sample ($4\%$) become passive while orbiting inside the cluster. These galaxies have experienced heavy pre-processing and the cluster's tidal stripping and ram pressure provide the final blow to their star formation. Finally, galaxies that are passive before entering the cluster for the first time ($12\%$) are typically massive and are not affected significantly by the cluster. Using the bulge/total mass ratio as a proxy for morphology, we find that a single incursion through a cluster do not result in significant morphological changes in all four types.

Pau Amaro Seoane

17 pages, submitted

We have evidence of X-ray flares in several galaxies consistent with a a star being tidally disrupted by a supermassive black hole (MBH). If the star starts on a nearly parabolic orbit relative to the MBH, one can derive that the fallback rate follows a $t^{-5/3}$ decay in the bolometric luminosity. We have modified the standard version of the smoothed-particle hydrodynamics (SPH) code {\sc Gadget} to include a relativistic treatment of the gravitational forces. We include non-spinning post-Newtonian corrections to incorpore the periapsis shift and the spin-orbit coupling up to next-to-lowest order. We run a set of simulations for different penetration factors in both the Newtonian- and the relativistic regime. We find that tidal disruptions around MBHs in the relativistic cases are underluminous for values starting at $\beta \gtrapprox 2.25$; i.e. the fallback curves produced in the relativistic cases are progressively lower compared to the Newtonian simulations as the penetration parameter increases. This is due to the fact that, contrary to the Newtonian cases, we find that all relativistic counterparts feature a survival core for penetration factors going to values as high as $12.05$. We derive a relativistic calculation which shows that geodesics of the elements in the star converge as compared to the Newtonian case, allowing for a core to survive the tidal disruption. A survival core should consistently emerge from any TDE with $\beta \gtrapprox 2.25$. The higher the value, the lower the colour temperatures than derived from standard accretion models.

R. Abbasi, M. Ackermann, J. Adams, S. K. Agarwalla, J. A. Aguilar, M. Ahlers, J.M. Alameddine, N. M. Amin, K. Andeen, G. Anton, C. Argüelles, Y. Ashida, S. Athanasiadou, S. N. Axani, X. Bai, A. Balagopal V., M. Baricevic, S. W. Barwick, V. Basu, R. Bay, J. J. Beatty, J. Becker Tjus, J. Beise, C. Bellenghi, C. Benning, S. BenZvi, D. Berley, E. Bernardini, D. Z. Besson, E. Blaufuss, S. Blot, F. Bontempo, J. Y. Book, C. Boscolo Meneguolo, S. Böser, O. Botner, J. Böttcher, E. Bourbeau, J. Braun, B. Brinson, J. Brostean-Kaiser, R. T. Burley, R. S. Busse, D. Butterfield, M. A. Campana, K. Carloni, E. G. Carnie-Bronca, S. Chattopadhyay, N. Chau, C. Chen, Z. Chen, D. Chirkin, S. Choi, B. A. Clark, L. Classen, A. Coleman, G. H. Collin, A. Connolly, J. M. Conrad, P. Coppin, et al. (345 additional authors not shown)

To access the list of contributions, please follow the "HTML" link. Links to individual contributions will fill in as authors upload their material

The IceCube Observatory at the South Pole has been operating in its full configuration since May 2011 with a duty cycle of about 99%. Its main component consists of a cubic-kilometer array of optical sensors deployed deep in the Glacial ice designed for the detection of high-energy astrophysical neutrinos. A surface array for cosmic ray air shower detection, IceTop, and a denser inner subdetector, DeepCore, significantly enhance the capabilities of the observatory, making it a multipurpose facility. This list of contributions to the 38th International Cosmic Ray Conference in Nagoya, Japan (July 26 - August 3, 2023) summarizes the latest results from IceCube covering a broad set of key questions in physics and astrophysics. The papers in this index are grouped topically to highlight IceCube contributions related to high-energy neutrino and multi-messenger astrophysics, cosmic-ray physics, low-energy neutrino transients such as Galactic supernovae, fundamental physics, detector calibration and event reconstruction, education and public outreach, and research and development for the IceCube Upgrade, a scheduled dense sensor infill complemented by calibration devices. Contributions related to IceCube-Gen2, the future extension of IceCube, are available in a separate collection.

R. Abbasi, M. Ackermann, J. Adams, S. K. Agarwalla, J. A. Aguilar, M. Ahlers, J.M. Alameddine, N. M. Amin, K. Andeen, G. Anton, C. Argüelles, Y. Ashida, S. Athanasiadou, J. Audehm, S. N. Axani, X. Bai, A. Balagopal V., M. Baricevic, S. W. Barwick, V. Basu, R. Bay, J. Becker Tjus, J. Beise, C. Bellenghi, C. Benning, S. BenZvi, D. Berley, E. Bernardini, D. Z. Besson, A. Bishop, E. Blaufuss, S. Blot, M. Bohmer, F. Bontempo, J. Y. Book, J. Borowka, C. Boscolo Meneguolo, S. Böser, O. Botner, J. Böttcher, S. Bouma, E. Bourbeau, J. Braun, B. Brinson, J. Brostean-Kaiser, R. T. Burley, R. S. Busse, D. Butterfield, M. A. Campana, K. Carloni, E. G. Carnie-Bronca, M. Cataldo, S. Chattopadhyay, N. Chau, C. Chen, Z. Chen, D. Chirkin, S. Choi, B. A. Clark, R. Clark, L. Classen, et al. (394 additional authors not shown)

To access the list of contributions, please follow the "HTML" link. Links to individual contributions will fill in as authors upload their material

IceCube-Gen2 is a planned next-generation neutrino observatory at the South Pole that builds upon the successful design of IceCube. Integrating two complementary detection technologies for neutrinos, optical and radio Cherenkov emission, in combination with a surface array for cosmic ray air shower detection, IceCube-Gen2 will cover a broad neutrino energy range from MeV to EeV. This index of contributions to the 38th International Cosmic Ray Conference in Nagoya, Japan (July 26 - August 3, 2023) describes research and development efforts for IceCube-Gen2. Included are summaries of the design, status, and sensitivity of the IceCube-Gen2 optical, surface, and radio components; performance studies of next-generation optical sensors detecting optical Cherenkov radiation from cosmic ray and neutrino events; reconstruction techniques of radio and optical events in terms of energy, direction, and neutrino flavor; and sensitivity studies of astrophysical neutrino flavors, diffuse neutrino fluxes, and cosmic ray anisotropies. Contributions related to IceCube and the scheduled IceCube Upgrade are available in a separate collection.

F. D. Priestley, P. C. Clark, S. C. O. Glover, S. E. Ragan, O. Fehér, L. R. Prole, R. S. Klessen

14 pages, 13 figures. MNRAS accepted

Much of what we know about molecular clouds, and by extension star formation, comes from molecular line observations. Interpreting these correctly requires knowledge of the underlying molecular abundances. Simulations of molecular clouds typically only model species that are important for the gas thermodynamics, which tend to be poor tracers of the denser material where stars form. We construct a framework for post-processing these simulations with a full time-dependent chemical network, allowing us to model the behaviour of observationally-important species not present in the reduced network used for the thermodynamics. We use this to investigate the chemical evolution of molecular gas under realistic physical conditions. We find that molecules can be divided into those which reach peak abundances at moderate densities ($10^3 \, {\rm cm^{-3}}$) and decline sharply thereafter (such as CO and HCN), and those which peak at higher densities and then remain roughly constant (e.g. NH$_3$, N$_2$H$^+$). Evolving the chemistry with physical properties held constant at their final values results in a significant overestimation of gas-phase abundances for all molecules, and does not capture the drastic variations in abundance caused by different evolutionary histories. The dynamical evolution of molecular gas cannot be neglected when modelling its chemistry.

Ruskin Patel, Kaloyan Penev, Joshua Schussler

Accepted for publication in MNRAS 16 pages, 9 figures

Tidal dissipation in binary systems is the primary source for synchronization and circularization of the objects in the system. The efficiency of the dissipation of tidal energy inside stars or planets results in significant changes in observed properties of the binary system and is often studied empirically using a parameter, commonly known as the modified tidal quality factor ($Q_\star'$). Though often assumed constant, in general that parameter will depend on the particular tidal wave experiencing the dissipation and the properties of the tidally distorted object. In this work we study the frequency dependence of $Q_\star'$ for Sun-like stars. We parameterize $Q_\star'$ as a saturating power-law in tidal frequency and obtain constraints using the stellar rotation period of 70 eclipsing binaries observed by Kepler. We use Bayesian analysis to account for the uncertainties in the observational data required for tidal evolution. Our analysis shows that $Q_\star'$ is well constrained for tidal periods > 15 days, with a value of $Q_\star' \sim 10^8$ for periods > 30 days and a slight suggested decrease at shorter periods. For tidal periods < 15 days, $Q_\star'$ is no longer tightly constrained, allowing for a broad range of possible values that overlaps with the constraints obtained using tidal circularization in binaries, which point to much more efficient dissipation: $Q_\star' \sim 10^6$.

Pablo Lemos, Paul Shah

Invited chapter for the edited book "Hubble Constant Tension" (Eds. E. Di Valentino and D. Brout, Springer Singapore, expected in 2024)

The cosmic microwave background (CMB) offers a unique window into the early universe, providing insights into cosmological parameters like the Hubble constant. Recent precise measurements of the CMB by experiments like Planck seem to point to a lower value for the Hubble constant compared to some other measurements like those from Type Ia supernovae. This discrepancy, known as the Hubble tension, currently lacks a definitive explanation. In this chapter, we provide an overview of how the Hubble constant is determined from detailed measurements of the CMB power spectrum. We explain the physics underlying key features of the CMB spectrum and their connection to cosmological parameters. We then examine the consistency of Planck's Hubble constant determination, both internally within the data itself and externally with other astrophysical probes. While largely consistent, some anomalies like the lensing amplitude parameter $A_L$ remain unresolved. We also explore various theoretical extensions to the standard ${\Lambda}$CDM cosmological model and assess their potential to resolve the Hubble tension. No clear resolution emerges, indicating significant tensions remain between early and late universe probes within simple extensions to ${\Lambda}$CDM. Upcoming CMB experiments promise improved precision and should provide further insights into this cosmic conundrum. A coherent picture bridging measurements across cosmic time remains an open challenge at the forefront of modern cosmology.

Cristian Hernández Cuevas, Roberto E. González, Nelson D. Padilla

17 pages, 16 figures, to be published in MNRAS

Using data from TNG300-2, we train a neural network (NN) to recreate the stellar mass ($M^*$) and star formation rate (SFR) of central galaxies in a dark-matter-only simulation. We consider 12 input properties from the halo and sub-halo hosting the galaxy and the near environment. $M^*$ predictions are robust, but the machine does not fully reproduce its scatter. The same happens for SFR, but the predictions are not as good as for $M^*$. We chained neural networks, improving the predictions on SFR to some extent. For SFR, we time-averaged this value between $z=0$ and $z=0.1$, which improved results for $z=0$. Predictions of both variables have trouble reproducing values at lower and higher ends. We also study the impact of each input variable in the performance of the predictions using a leave-one-covariate-out approach, which led to insights about the physical and statistical relation between input variables. In terms of metrics, our machine outperforms similar studies, but the main discoveries in this work are not linked with the quality of the predictions themselves, but to how the predictions relate to the input variables. We find that previously studied relations between physical variables are meaningful to the machine. We also find that some merger tree properties strongly impact the performance of the machine. %We highlight the value of machine learning (ML) methods in helping understand the information contained in different variables, since with its help we were able to obtain useful insights resulting from studying the impact of input variables on the resulting behaviour of galaxy properties. We conclude that ML models are useful tools to understand the significance of physical different properties and their impact on target characteristics, as well as strong candidates for potential simulation methods.

We present a solution to Liouville's equation for an ensemble of charged particles propagating in magnetic fields. The solution is presented using an expansion in spherical harmonics of the phase space density, allowing a direct interpretation of the distribution of arrival directions of cosmic rays. The results are found for chosen conditions of variability and source distributions. We show there are two conditions for an initially isotropic flux of particles to remain isotropic while traveling through a magnetic field: isotropy and homogeneity of the sources. The formalism is used to analyze the data measured by the Pierre Auger Observatory, contributing to the understanding of the dependence of the dipole amplitude with energy and predicting the energy in which the quadrupole signal should be measured.

Recently, pulsar timing array (PTA) collaborations announced evidence for an isotropic stochastic gravitational wave (GW) background. The origin of the PTA signal can be astrophysical or cosmological. In the latter case, the so-called secondary scalar-induced GW scenario is one of the viable explanations, but it has a potentially serious issue of the overproduction of primordial black holes (PBHs) due to the enhanced curvature perturbation. In this letter, we present a new interpretation of the PTA signal. Namely, it is originated from an extra spectator tensor field that exists on top of the metric tensor perturbation. As the energy density of the extra tensor field is always subdominant, it cannot lead to the formation of PBHs. Thus our primordial-tensor-induced scenario is free from the PBH overproduction issue.

Henric Krawczynski (1), Yajie Yuan (1), Alexander Y. Chen (1), Nicole Rodriguez Cavero (1), Kun Hu (1), Ephraim Gau (1), James F. Steiner (2), Michal Dovčiak (3), ((1) Washington University in St. Louis, (2) Center for Astrophysics, Harvard and Smithsonian, (3) Astronomical Institute, Academy of Sciences of the Czech Republic)

8 pages and 6 figures, submitted to the Astrophysical Journal Letters

The Imaging X-ray Polarimetry Explorer (IXPE) observations of the X-ray binary 4U 1630-47 in the high soft state revealed linear polarization degrees (PDs) rising from 6% at 2 keV to 10% at 8 keV. Explaining the results in the framework of the standard optically thick, geometrically thin accretion disk scenario requires careful fine-tuning of the relevant model parameters. We argue here that the emission of polarized Bremsstrahlung by anisotropic electrons in the accretion disk atmosphere can account for the overall high PDs and the increase of the PDs with energy. We discuss plasma and accretion effects that can generate electron anisotropies at a level required by the 4U 1630-47 results. We conclude by emphasizing that X-ray polarimetry affords us the opportunity to obtain information about the magnetization of the accretion disk atmosphere.

S.Q. Wang, J.B. Wang, D.Z. Li, J.M. Yao, R.N. Manchester, G. Hobbs, N. Wang, S. Dai, H. Xu, R. Luo, Y. Feng, W.Y. Wang, D. Li, Y.W. Yu, Z.X. Du, C.H. Niu, S.B. Zhang, C.M. Zhang

7 pages, 3 figures, accept for publication in ApJ

Black widows are millisecond pulsars ablating their companions. The material blown from the companion blocks the radio emission, resulting in radio eclipses. The properties of the eclipse medium are poorly understood. Here, we present direct evidence of the existence of magnetic fields in the eclipse medium of the black widow PSR J2051$-$0827 using observations made with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). We detect a regular decrease in rotation measure (RM) in the egress of eclipse, changing from $60\,\rm rad\,m^{-2}$ to $-28.7\,\rm rad\,m^{-2}$. The RM gradually changes back to normal when the line-of-sight moves away from the eclipse. The estimated line-of-sight magnetic field strength in the eclipse medium is $\sim 0.1$ G. The RM reversal could be caused by a change of the magnetic field strength along the line of sight due to binary orbital motion. The RM reversal phenomenon has also been observed in some repeating fast radio bursts (FRBs), and the study of spider pulsars may provide additional information about the origin of FRBs.

Satabdwa Majumdar, Ritesh Patel, Vaibhav Pant, Dipankar Banerjee, Aarushi Rawat, Abhas Pradhan, Paritosh Singh

29 Pages, 18 Figures. Accepted in The Astrophysical Journal Supplement Series (APJS)

The primary objective of this study is to connect the coronal mass ejections (CMEs) to their source regions, primarily creating a CME source region (CSR) catalogue, and secondly probing into the influence the source regions have on different statistical properties of CMEs. We create a source region catalogue for 3327 CMEs from 1998 to 2017, thus capturing the different phases of cycle 23 and 24. The identified source regions are segregated into 3 classes, Active Regions (ARs), Prominence Eruptions (PEs) and Active Prominences (APs), while the CMEs are segregated into slow and fast based on their average projected speeds. We find the contribution of these three source region types to the occurrences of slow and fast CMEs to be different in the above period. A study of the distribution of average speeds reveals different power-laws for CMEs originating from different sources, and the power-law is different during the different phases of cycles 23 and 24. A study of statistical latitudinal deflections showed equator-ward deflections, while the magnitude of deflections again bears an imprint of the source regions. An East-West asymmetry is also noted, particularly in the rising phase of cycle 23, with the presence of active longitudes for the CMEs, with a preference towards the Western part of the Sun. Our results show that different aspects of CME kinematics bear a strong imprint of the source regions they originate from, thus indicating the existence of different ejection and/or propagation mechanisms of these CMEs.

Xing Xu, Jianhua Guo, Shen Wang (for the GRAND Collaboration)

Presented at the 38th International Cosmic Ray Conference (ICRC 2023)

GRANDProto300 is the planned 300-antenna pathfinder array of the Giant Radio Array for Neutrino Detection (GRAND), of which the first 100 detection units have been already produced. Its main goal is to demonstrate the viability of the detection of the radio emission from air showers initiated by inclined ultra-high-energy cosmic rays with energies of $10^{16.5}$ to $10^{18.5}$ eV, covering the purported transition region from their Galactic to extragalactic origin. The front-end readout system of each detection unit of GRANDProto300 processes signals from the radio antenna and the particle detector, generates the first-level trigger, and communicates with a central processing station. Based on earlier designs, we have built the first prototype of this system using two development boards and one self-designed front-end board. We present our new design that is improved and more economical than the earlier one, as well as test results and prospects for future work.

Yanan Feng, Xiaohu Li, Tom J. Millar, Ryszard Szczerba, Ke Wang, Donghui Quan, Shengli Qin, Xuan Fang, Juan Tuo, Zhenzhen Miao, Rong Ma, Fengwei Xu, Jingfei Sun, Biwei Jiang, Qiang Chang, Jianchao Yang, Gao-Lei Hou, Fangfang Li, Yong Zhang

Accepted in Frontiers in Astronomy and Space Sciences

Whether SiC$_2$ is a parent species, that is formed in the photosphere or as a by-product of high-temperature dust formation, or a daughter species, formed in a chemistry driven by the photodestruction of parent species in the outer envelope, has been debated for a long time. Here, we analyze the ALMA observations of four SiC$_2$ transitions in the CSEs of three C-rich AGB stars (AI Vol, II Lup, and RAFGL 4211), and found that SiC$_2$ exhibits an annular, shell-like distribution in these targets, suggesting that SiC$_2$ can be a daughter species in the CSEs of carbon-rich AGB stars. The results can provide important references for future chemical models.

Kai-Kai Duan, Peng-Xiong Ma, Ke-Wen Zhang, Xiao-Yuan Huang, Yi Zhang (for the GRAND collaboration)

Presented at the 38th International Cosmic Ray Conference (ICRC 2023), 8 pages, 9 figures

GRANDProto300 is a 300-antenna prototype array of the envisioned GRAND (Giant Radio Array for Neutrino Detection) project. The goal of GRANDProto300 is to detect radio signals emitted by cosmic ray-induced air showers, with energies ranging from $10^{16.5}$~eV to $10^{18.5}$~eV, which covers the transition region between Galactic and extragalactic sources. We use simulations to optimize the layout of GRANDProto300 and develop a shower reconstruction method. Based on them, we present the performance of GRANDProto300 for cosmic-ray detection, by means of its effective area, angular resolution, and energy resolution.

Kai Xiao, Haibo Yuan, Bowen Huang, Shuai Xu, Jie Zheng, Chun Li, Zhou Fan, Wei Wang, Gang Zhao, Guojie Feng, Xuan Zhang, Jinzhong Liu, Ruoyi Zhang, Lin Yang, Yu Zhang, Chunhai Bai, Hubiao Niu, Esamdin Ali, Lu Ma

15 pages in Chinese language, 8 figures, Chinese Science Bulletin accepted and published online ( this https URL ), see main results in Figures 6, 7 and 8

In this paper, a total of approximately 2.6 million dwarfs were constructed as standard stars, with an accuracy of about 0.01-0.02 mag for each band, by combining spectroscopic data from the Large Sky Area Multi-Object Fiber Spectroscopic Telescope Data Release 7, photometric data from the corrected Gaia Early Data Release 3, and photometric metallicities. Using the spectroscopy based stellar color regression method (SCR method) and the photometric-based SCR method (SCR' method), we performed the relative calibration of the Nanshan One-meter Wide-field Telescope imaging data. Based on the corrected Pan-STARRS DR1 photometry, the absolute calibration was also performed. In the photometric calibration process, we analyzed the dependence of the calibration zero points on different images (observation time), different gates of the CCD detector, and different CCD positions. We found that the stellar flat and the relative gain between different gates depend on time. The amplitude of gain variation in three channels is approximately 0.5%-0.7% relative to the other channel, with a maximum value of 4%. In addition, significant spatial variations of the stellar flat fitting residual are found and corrected. Using repeated sources in the adjacent images, we checked and discovered internal consistency of about 1-2 mmag in all the filters. Using the PS1 magnitudes synthesized by Gaia DR3 BP/RP spectra by the synthetic photometry method, we found that the photometric calibration uniformity is about 1-2 mmag for all the bands, at a spatial resolution of 1.3 degree. A detailed comparison between the spectroscopy-based SCR and photometric-based SCR method magnitude offsets was performed, and we achieved an internal consistency precision of about 2 mmag or better with resolutions of 1.3 degree for all the filters. Which is mainly from the position-dependent errors of the E(B-V) used in SCR' method.

Sophia V. Sosa Fiscella, Michael T. Lam, Zaven Arzoumanian, Harsha Blumer, Paul R. Brook, H. Thankful Cromartie, Megan E. DeCesar, Paul B. Demorest, Timothy Dolch, Justin A. Ellis, Robert D. Ferdman, Elizabeth C. Ferrara, Emmanuel Fonseca, Nate Garver-Daniels, Peter A. Gentile, Megan L. Jones, Duncan R. Lorimer, Ryan S. Lynch, Maura A. McLaughlin, Cherry Ng, David J. Nice, Timothy T. Pennucci, Scott M. Ransom, Renée Spiewak, Ingrid H. Stairs, Kevin Stovall, Joseph K. Swiggum, Sarah J. Vigeland, WeiWei Zhu

15 pages, 7 figures

Noise characterization for pulsar-timing applications accounts for interstellar dispersion by assuming a known frequency-dependence of the delay it introduces in the times of arrival (TOAs). However, calculations of this delay suffer from mis-estimations due to other chromatic effects in the observations. The precision in modeling dispersion is dependent on the observed bandwidth. In this work, we calculate the offsets in infinite-frequency TOAs due to mis-estimations in the modeling of dispersion when using varying bandwidths at the Green Bank Telescope. We use a set of broadband observations of PSR J1643-1224, a pulsar with an excess of chromatic noise in its timing residuals. We artificially restricted these observations to a narrowband frequency range, then used both data sets to calculate residuals with a timing model that does not include short-scale dispersion variations. By fitting the resulting residuals to a dispersion model, and comparing the ensuing fitted parameters, we quantify the dispersion mis-estimations. Moreover, by calculating the autocovariance function of the parameters we obtained a characteristic timescale over which the dispersion mis-estimations are correlated. For PSR J1643-1224, which has one of the highest dispersion measures (DM) in the NANOGrav pulsar timing array, we find that the infinite-frequency TOAs suffer from a systematic offset of ~22 microseconds due to DM mis-estimations, with correlations over ~1 month. For lower-DM pulsars, the offset is ~7 microseconds. This error quantification can be used to provide more robust noise modeling in NANOGrav's data, thereby increasing sensitivity and improving parameter estimation in gravitational wave searches.

Payel Nandi, C. S. Stalin, D. J. Saikia, Rogemar A. Riffel, Arijit Manna, Sabyasachi Pal, O. L. Dors, Dominika Wylezalek, Vaidehi S. Paliya, P. Saikia, Pratik Dabhade, Markus-Kissler Patig, Ram Sagar

Submitted. Comments are welcome

Black hole driven outfBlack hole driven outflows in galaxies hosting active galactic nuclei (AGN) may interact with their interstellar medium (ISM) affecting star formation. Such feedback processes, reminiscent of those seen in massive galaxies, have been reported recently in some dwarf galaxies. However, such studies have usually been on kiloparsec and larger scales and our knowledge on the smallest spatial scales to which this feedback processes can operate is unclear. Here we demonstrate radio jet-ISM interaction on the scale of an asymmetric triple radio structure of $\sim$10 parsec size in NGC 4395. This triplet radio structure is seen in the 15 GHz continuum image and the two asymmetric jet like structures are situated on either side of the radio core that coincides with the optical {\it Gaia} position. The high resolution radio image and the extended [OIII]$\lambda$5007 emission, indicative of an outflow, are spatially coincident and are consistent with the interpretation of a low power radio jet interacting with the ISM. Modelling of the spectral lines using CLOUDY and MAPPINGS, and estimation of temperature using Gemini and MaNGA integral field spectroscopic data suggest shock ionization of the gas. The continuum emission at 237 GHz, though weak was found to spatially coincide with the AGN, however, the CO(2-1) line emission was found to be displaced by around 1 arcsec northward of the AGN core. The spatial coincidence of molecular H2$\lambda$2.4085 along the jet direction, the morphology of ionised [OIII]$\lambda$5007 and displacement of CO(2-1) emission argues for conditions less favourable for star formation at $\sim$5 parsec.

Kosuke Nishiwaki, Katsuaki Asano, Kohta Murase

12 pages, 6 figures, accepted for publication in ApJ

A fraction of merging galaxy clusters host diffuse radio emission in their central region, termed as a giant radio halo (GRH). The most promising mechanism of GRHs is the re-acceleration of non-thermal electrons and positrons by merger-induced turbulence. However, the origin of these seed leptons has been under debate, and either protons or electrons can be primarily-accelerated particles. In this work, we demonstrate that neutrinos can be used as a probe of physical processes in galaxy clusters, and discuss possible constraints on the amount of relativistic protons in the intra-cluster medium with the existing upper limits by IceCube. We calculate radio and neutrino emission from massive ($>10^{14}M_\odot$) galaxy clusters, using the cluster population model of Nishiwaki & Asano (2022). This model is compatible with the observed statistics of GRHs, and we find that the contribution of GRHs to the isotropic radio background observed with the ARCADE-2 experiment should be subdominant. Our fiducial model predicts the all-sky neutrino flux that is consistent with IceCube's upper limit from the stacking analysis. We also show that the neutrino upper limit gives meaningful constraints on the parameter space of the re-acceleration model, such as the electron-to-proton ratio of primary cosmic-rays and the magnetic field, and in particular the secondary scenario, where the seed electrons mostly originate from inelastic $pp$ collisions, can be constrained even in the presence of re-acceleration.

Thomas Lang, Samuel T. Spencer, Alison M.W. Mitchell

11 pages, 14 figures, 5 tables. Accepted in A&A. Reproduced with permission from Astronomy & Astrophysics, $\copyright$ ESO

The number of satellites launched into Earth orbit has almost tripled in the last three years (to over 4000) due to the increasing commercialisation of space. Multiple satellite constellations, consisting of over 400,000 individual satellites, have either been partially launched or are proposed for launch in the near future. Many of these satellites are highly reflective, resulting in a high optical brightness that affects ground-based astronomical observations. Despite this, the potential effect of these satellites on gamma-ray-observing Imaging Atmospheric Cherenkov Telescopes (IACTs) has largely been assumed to be negligible due to their nanosecond-scale integration times. This has, however, never been verified. As IACTs are sensitive to optical wavelength light, we aim to identify satellite trails in data taken by the High Energy Stereoscopic System (H.E.S.S.) IACT array. This is to quantify the potential effects on data quality and extensive air shower event classification and reconstruction. Using night sky background measurements from H.E.S.S., we determine which observation times and pointing directions are affected most by these satellite trails, and evaluate the impact on the standard Hillas parameter variables used for event analysis. Due to the brightest trails, false trigger events can occur, however for most modern analyses the effect on astronomical results will be minimal. We observe a mild increase in the rate of trail detections over time (approximately doubling in three years), which is partially correlated with the number of satellite launches. But the fraction of H.E.S.S. data affected ($\sim0.2\%$ of dark time observations) is currently small. Nevertheless, these trails could have a non-negligible effect on future Cherenkov Telescope Array observations if advanced analysis techniques designed to lower the energy threshold of the instrument are used.

Yael Naze (FNRS/ULiege), Jan Robrade (Hamburger Sternwarte)

accepted for publication by MNRAS

Massive stars are known X-ray emitters and those belonging to the Be category are no exception. One type of X-ray emission even appears specific to that category, the gamma Cas phenomenon. Its actual incidence has been particularly difficult to assess. Thanks to four semesters of sky survey data taken by SRG (Spectrum Roentgen Gamma)/eROSITA, we revisit the question of the X-ray properties of Be stars. Amongst a large catalog of Be stars, eROSITA achieved 170 detections (20% of sample), mostly corresponding to the earliest spectral types and/or close objects. While X-ray luminosities show an uninterrupted increasing trend with the X-ray-to-bolometric luminosity ratios, the X-ray hardness was split between a large group of soft (and fainter on average) sources and a smaller group of hard (and brighter on average) sources. The latter category gathers at least 34 sources, nearly all displaying early spectral types. Only a third of them were known before to display such X-ray properties. The actual incidence of hard and bright X-rays amongst early-type Be stars within 100--1000pc appears to be ~12%, which is far from negligible. At the other extreme, no bright supersoft X-ray emission seem to be associated to any of our targets.

Context. Recent spacecraft observations in the inner heliosphere have revealed the presence of local Alfvenic reversals of the magnetic field, while the field magnitude remains almost constant. They are called magnetic switchbacks and are very common in the plasma environment close to the Sun explored by the Parker Solar Probe satellite. Aims. A simple numerical model of a magnetic field reversal with constant magnitude is used in order to explore the influence of switchbacks on the propagation of energetic particles, within a range of energy typical of solar energetic particles. Methods. We model the reversal as a region of space of adjustable size bounded by two rotational discontinuities. By means of test particle simulations, beams of mono-energetic particles can be injected upstream of the switchback with various initial pitch- and gyro-phase angles. In each simulation, the particle energy may also be changed. Results. Particle dynamics is highly affected by the ratio between the particle gyroradius and the size of the switchback, with multiple pitch-angle scatterings when the particle gyroradius is of the order of the switchback size. Further, particle motion is extremely sensitive to the initial conditions implying a transition to chaos; for some parameters of the system, a large share of particles is reflected backwards upstream as they interact with the switchback. These results can have a profound impact on the solar energetic particle transport in the inner heliosphere, thus possible comparisons with in-situ spacecraft data are discussed.

Sunil Malik, Ka Ho Yuen, Huirong Yan

6 pages, 3 figures, 1 table

Magnetic fields are ubiquitous in the interstellar medium, including extended objects such as supernova remnants (SNRs) and Pulsar Wind Nebulae (PWNe). Its turbulent characteristics govern the diffusion of cosmic rays and the multi-wavelength emission from PWNe. However, the geometry and turbulence nature of the magnetic fields in the ambient region of PWN is still unknown. Recent gamma-ray observations from HAWC and synchrotron observations suggest a highly suppressed diffusion coefficient compared to the mean interstellar value. In this letter, we present the first direct observational evidence that the local mean magnetic field is nearly aligned toward the line of sight (LoS) with an inclination angle $\theta_{\lambda} <10^{\circ}$ employing a recently developed statistical recipe known as `Y-parameter'. Furthermore, we report that the magnetic field fluctuations are mostly dominated by compressible modes, with a 2D correlation length of approximately $3 \ {\rm pc}$ in the vicinity of Monogem PWN region. Our study highlights the pivotal role of magnetic field and turbulence in unraveling the physical processes in TeV halos and cosmic ray transport.

T. J. Millar

13 pages, 2 figures, early on-line publication

Fifty years on from the first detailed chemical kinetic modelling of astronomical sources, I provide some introductory comments on the history of astrochemistry, summarise some personal views on the topics covered in this discussion meeting, and conclude with some thoughts on its future development. I have left out the jokes.

W. Mercier, B. Epinat, T. Contini, D. Krajnović, L. Ciesla, B. C. Lemaux, V. Abril-Melgarejo, L. Boogaard, D. Pelliccia

23 pages, including 5 appendices, 14 figures, accepted version before language correction

Aims: At intermediate redshift, galaxy groups/clusters are thought to impact galaxies (e.g. their angular momentum). We investigate whether the environment has an impact on the galaxies' angular momentum and identify underlying driving physical mechanisms. Methods: We derive robust estimates of the stellar angular momentum using Hubble Space Telescope (HST) images combined with spatially resolved ionised gas kinematics from the Multi-Unit Spectroscopic Explorer (MUSE) for a sample of ~200 galaxies in groups and in the field at z~0.7 drawn from the MAGIC survey. Using various environmental tracers, we study the position of the galaxies in the the angular momentum-stellar mass (Fall) relation as a function of environment. Results: We measure a 0.12 dex (2sigma significant) depletion of angular momentum for low-mass galaxies (M* < 10^10 Msun) in groups with respect to the field. Massive galaxies located in dense environments have less angular momentum than expected from the low-mass Fall relation but, without a comparable field sample, we cannot infer whether this effect is mass- or environmentally-driven. Furthermore, massive galaxies are found in the centre of the structures and have low systemic velocities. The observed depletion of angular momentum at low mass does not appear linked with the strength of the over-density around the galaxies but it is strongly correlated with the galaxies' systemic velocity normalised by the dispersion of their host group and with their ionised gas velocity dispersion. Conclusions: Group galaxies seem depleted in angular momentum, especially at low mass. Our results suggest that this depletion might be induced by physical mechanisms that scale with the systemic velocity of the galaxies (e.g. stripping or merging) and that such mechanism might be responsible for enhancing the velocity dispersion of the gas as galaxies lose angular momentum.

Yoon-Hyun Ryu, Andrzej Udalski, Jennifer C. Yee, Weicheng Zang, Yossi Shvartzvald, Cheongho Han, Andrew Gould, Michael D. Albrow, Sun-Ju Chung, Kyu-Ha Hwang, Youn Kil Jung, In-Gu Shin, Hongjing Yang, Sang-Mok Cha, Dong-Jin Kim, Seung-Lee Kim, Chung-Uk Lee, Dong-Joo Lee, Yongseok Lee, Byeong-Gon Park, Richard W. Pogge, Hanyue Wang, Przemek Mróz, Michał K. Szymański, Jan Skowron, Radek Poleski, Igor Soszyński, Paweł Pietrukowicz, Szymon Kozłowski, Krzysztof Ulaczyk, Krzysztof A. Rybicki, Patryk Iwanek, Marcin Wrona, Charles Beichman, Geoffry Bryden, Sean Carey, Calen B. Henderson, Sebastiano Calchi Novati, Wei Zhu, Savannah Jacklin, Matthew T.Penny

67 pages, 13 figures, 16 tables

We complete the analysis of planetary candidates found by the KMT AnomalyFinder for the 2017 prime fields that cover $\sim 13\,{\rm deg}^2$. We report 3 unambiguous planets: OGLE-2017-BLG-0640, OGLE-2017-BLG-1275, and OGLE-2017-BLG-1237. The first two of these were not previously identified, while the last was not previously published due to technical complications induced by a nearby variable. We further report that a fourth anomalous event, the previously recognized OGLE-2017-BLG-1777, is very likely to be planetary, although its light curve requires unusually complex modeling because the lens and source both have orbiting companions. One of the 3 unambiguous planets, OGLE-2017-BLG-1275 is the first AnomalyFinder discovery that has a {\it Spitzer} microlens parallax measurement, $\pi_E \sim 0.045\pm0.015$, implying that this planetary system almost certainly lies in the Galactic bulge. In the order listed, the four planetary events have planet-host mass ratios $q$, and normalized projected separations $s$, of $(\log q,s)$ = $(-2.31,0.61)$, $(-2.06,0.63/1.09)$, $(-2.10,1.04)$, and $(-2.86,0.72)$. Combined with previously published events, the 2017 AnomalyFinder prime fields contain 11 unambiguous planets with well-measured $q$ and one very likely candidate, of which 3 are AnomalyFinder discoveries. In addition to these 12, there are three other unambiguous planets with large uncertainties in $q$.

Mudit Garg, Shubhanshu Tiwari, Andrea Derdzinski, John Baker, Sylvain Marsat, Lucio Mayer

11 pages, 12 figures. Submitted to MNRAS

We explore the eccentricity measurement threshold of LISA for gravitational waves radiated by massive black hole binaries (MBHBs) with redshifted BH masses $M_z$ in the range $10^{4.5}\mbox{-}10^{7.5}~{\rm M}_\odot$ at redshift $z=1$. The eccentricity can be an important tracer of the environment where MBHBs evolve to reach the merger phase. To consider LISA's motion and apply the time delay interferometry, we employ the lisabeta software and produce year-long eccentric waveforms using the inspiral-only post-Newtonian model TaylorF2Ecc. We study the minimum measurable eccentricity ($e_{\rm min}$, defined at one year before the merger) analytically by computing matches and Fisher matrices, and numerically via Bayesian inference by varying both intrinsic and extrinsic parameters. We find that $e_{\rm min}$ has a strong dependence on $M_z$ and a weak dependence on mass ratio and extrinsic parameters. Match-based signal-to-noise ratio criterion suggest that LISA will be able to detect $e_{\rm min}\sim10^{-2.5}$ for lighter systems ($M_z\lesssim10^{5.5}~{\rm M}_\odot$) and $\sim10^{-1.5}$ for heavier MBHBs with a $90\%$ confidence. Bayesian inference with Fisher initialization and a zero noise realization pushes this limit to $e_{\rm min}\sim10^{-2.75}$ for lower-mass binaries assuming a $<50\%$ relative error. Bayesian inference can recover injected eccentricities of $0.1$ and $10^{-2.75}$ for a $10^5~{\rm M}_\odot$ system with a $\sim10^{-2}\%$ and a $\sim10\%$ relative errors, respectively. Both analytical and numerical methodologies provide almost consistent results for our systems of interest. LISA will launch in a decade, making this study valuable and timely to prepare for unlocking the mysteries of the MBHB evolution.

Aneesh P. Naik, Axel Widmark

11 pages, 9 figures (excluding appendices). Submitted to MNRAS. RV catalogues available for download at this http URL . Source code, plotting scripts etc. accessible at this https URL

In an earlier work, we demonstrated the effectiveness of Bayesian neural networks in estimating the missing line-of-sight velocities of Gaia stars, and published an accompanying catalogue of blind predictions for the line-of-sight velocities of stars in Gaia DR3. These were not merely point predictions, but probability distributions reflecting our state of knowledge about each star. Here, we verify that these predictions were highly accurate: the DR3 measurements were statistically consistent with our prediction distributions, with an approximate error rate of 1.5%. We use this same technique to produce a publicly available catalogue of predictive probability distributions for the 185 million stars up to a G-band magnitude of 17.5 still missing line-of-sight velocities in Gaia DR3. Validation tests demonstrate that the predictions are reliable for stars within approximately 7 kpc from the Sun and with distance precisions better than around 20%. For such stars, the typical prediction uncertainty is 25-30 km/s. We invite the community to use these radial velocities in analyses of stellar kinematics and dynamics, and give an example of such an application.

Gravity is the driving force of star formation. Although gravity is caused by the presence of matter, its role in complex regions is still unsettled. One effective way to study the pattern of gravity is to compute the accretion it exerts on the gas by providing gravitational acceleration maps. A practical way to study acceleration is by computing it using 2D surface density maps, yet whether these maps are accurate remains uncertain. Using numerical simulations, we confirm that the accuracy of the acceleration maps $\mathbf a_{\rm 2D}(x,y)$ computed from 2D surface density are good representations for the mean acceleration weighted by mass. Due to the under-estimations of the distances from projected maps, the magnitudes of accelerations will be over-estimated $|\mathbf a_{\rm 2D}(x,y)| \approx 2.3 \pm 1.8 \; |\mathbf a_{\rm 3D}^{\rm proj}(x,y)|$, where $\mathbf a_{\rm 3D}^{\rm proj}(x,y)$ is mass-weighted projected gravitational acceleration, yet $\mathbf a_{\rm 2D}(x,y)$ and $ \mathbf a_{\rm 3D}^{\rm proj}(x,y)$ stay aligned within 20$^{\circ}$. Significant deviations only occur in regions where multiple structures are present along the line of sight. The acceleration maps estimated from surface density provide good descriptions of the projection of 3D acceleration fields. We expect this technique useful in establishing the link between cloud morphology and star formation, and in understanding the link between gravity and other processes such as the magnetic field. A version of the code for calculating surface density gravitational potential is available at \url{https://github.com/zhenzhen-research/phi_2d}.

Leon Ofman, Scott A Boardsen, Lan K Jian, Parisa Mostafavi, Jaye L Verniero, Roberto Livi, Michael McManus, Ali Rahmati, Davin Larson, Michael L Stevens

Accepted for publication in The Astrophysical Journal

Past observations show that solar wind (SW) acceleration occurs inside the sub-Alfvenic region, reaching the local Alfven speed at typical distances ~ 10 - 20 Rs (solar radii). Recently, Parker Solar Probe (PSP) traversed regions of sub-Alfvenic SW near perihelia in encounters E8-E12 for the first time providing data in these regions. It became evident that properties of the magnetically dominated SW are considerably different from the super-Alfvenic wind. For example, there are changes in relative abundances and drift of alpha particles with respect to protons, as well as in the magnitude of magnetic fluctuations. We use data of the magnetic field from the FIELDS instrument, and construct ion velocity distribution functions (VDFs) from the sub-Alfvenic regions using Solar Probe Analyzer Ions (SPAN-I) data, and run 2.5D and 3D hybrid models of proton-alpha sub-Alfvenic SW plasma. We investigate the nonlinear evolution of the ion kinetic instabilities in several case studies, and quantify the transfer of energy between the protons, alpha particles, and the kinetic waves. The models provide the 3D ion VDFs at the various stages of the instability evolution in the SW frame. By combining observational analysis with the modeling results, we gain insights on the evolution of the ion instabilities, the heating and the acceleration processes of the sub-Alfvenic SW plasma and quantify the exchange of energy between the magnetic and kinetic components. The modeling results suggest that the ion kinetic instabilities are produced locally in the SW, resulting in anisotropic heating of the ions, as observed by PSP.

P. Weber, S. Pérez, A. Zurlo, J. Miley, A. Hales, L. Cieza, D. Principe, M. Cárcamo, A. Garufi, Á. Kóspál, M. Takami, J. Kastner, Z. Zhu, J. Williams

Published in ApJL. 11 pages, 4 figures. Images available as fits files at this https URL

The formation of giant planets has traditionally been divided into two pathways: core accretion and gravitational instability. However, in recent years, gravitational instability has become less favored, primarily due to the scarcity of observations of fragmented protoplanetary disks around young stars and low occurrence rate of massive planets on very wide orbits. In this study, we present a SPHERE/IRDIS polarized light observation of the young outbursting object V960 Mon. The image reveals a vast structure of intricately shaped scattered light with several spiral arms. This finding motivated a re-analysis of archival ALMA 1.3 mm data acquired just two years after the onset of the outburst of V960 Mon. In these data, we discover several clumps of continuum emission aligned along a spiral arm that coincides with the scattered light structure. We interpret the localized emission as fragments formed from a spiral arm under gravitational collapse. Estimating the mass of solids within these clumps to be of several Earth masses, we suggest this observation to be the first evidence of gravitational instability occurring on planetary scales. This study discusses the significance of this finding for planet formation and its potential connection with the outbursting state of V960 Mon.

Elena M. Golubeva, Akash Biswas, Anna I. Khlystova, Pawan Kumar, Bidya Binay Karak

Accepted for publication in MNRAS

The polar field reversal is a crucial process in the cyclic evolution of the large-scale magnetic field of the Sun.Various important characteristics of a solar cycle, such as its duration and strength, and also the cycle predictability, are determined by the polar field reversal time. While the regular measurements of solar magnetic field have been accumulated for more than half a century, there is no consensus in the heliophysics community concerning the interpretation of the Sun's polar field measurements and especially the determination of polar field reversal time. There exists a severe problem of non-reproducibility in the reported results even from studies of the same observational dataset, and this causes an obstacle to make more accurate forecasts of solar cycle. Here, we analyze the solar magnetograms from four instruments for the last four cycles, to provide a more correct interpretation of the polar field observations and to find more accurate time of the reversals. We show the absence of triple (multipolar) reversals in Cycles 21 - 24, significant variations in the time interval between reversals in the hemispheres and in the time interval between a reversal and a cycle beginning. In order to understand the origin of the reversal time variation, we perform Surface Flux Transport (SFT) simulations and find out that the presence of the 'anomalous' bipolar magnetic regions (BMRs) in different phases of a cycle can cause cycle-to-cycle variations of the reversal time within the similar range found in observations.

Tiziano Schiavone, Enea Di Dio, Giuseppe Fanizza

30 pages, 4 figures

Starting from a recently proposed framework for the evaluation of the cosmological averages, we evaluate the higher-order moments for the distribution of a given observable. Then, we explicitly discuss the case of the Hubble-Lema\^itre diagram and evaluate its skewness at the leading order in the cosmological perturbative expansion of the gravitational potential. In particular, we focus on perturbations of the luminosity distance due to gravitational lensing. Finally, we discuss our findings in view of recent numerical relativistic simulations, confirming that the skewness in the Hubble-Lema\^itre diagram primarily originates from the late-time matter bispectrum, with other line-of-sight projection effects being sub-dominant.

Yang A. Li, Luis C. Ho, Jinyi Shangguan, Ming-Yang Zhuang, Ruancun Li

36 pages, 24 figures

The broadband spectral energy distribution of a galaxy encodes valuable information on its stellar mass, star formation rate (SFR), dust content, and possible fractional energy contribution from nonstellar sources. We present a comprehensive catalog of panchromatic photometry, covering 17 bands from the far-ultraviolet to 500 $\mu$m, for 2685 low-redshift (z=0.01-0.11), massive ($M_* > 10^{10}\,M_\odot$) galaxies selected from the Stripe 82 region of the Sloan Digital Sky Survey, one of the largest areas with relatively deep, uniform observations over a wide range of wavelengths. Taking advantage of the deep optical coadded images, we develop a hybrid approach for matched-aperture photometry of the multi-band data. We derive robust uncertainties and upper limits for undetected galaxies, deblend interacting/merging galaxies and sources in crowded regions, and treat contamination by foreground stars. We perform spectral energy distribution fitting to derive the stellar mass, SFR, and dust mass, critically assessing the influence of flux upper limits for undetected photometric bands and applying corrections for systematic uncertainties based on extensive mock tests. Comparison of our measurements with those of commonly used published catalogs reveals good agreement for the stellar masses. While the SFRs of galaxies on the star-forming main sequence show reasonable consistency, galaxies in and below the green valley show considerable disagreement between different sets of measurements. Our analysis suggests that one should incorporate the most accurate and inclusive photometry into the spectral energy distribution analysis, and that care should be exercised in interpreting the SFRs of galaxies with moderate to weak star formation activity.

Yang A. Li, Luis C. Ho, Jinyi Shangguan

19 pages, 9 figures

Interactions and mergers play an important role in regulating the physical properties of galaxies, such as their morphology, gas content, and star formation rate (SFR). Controversy exists as to the degree to which these events, even gas-rich major mergers, enhance star formation activity. We study merger pairs selected from a sample of massive ($M_* \ge 10^{10}\,M_\odot$), low-redshift ($z = 0.01-0.11$) galaxies located in the Stripe 82 region of the Sloan Digital Sky Survey, using stellar masses, SFRs, and total dust masses derived from a new set of uniformly measured panchromatic photometry and spectral energy distribution analysis. The dust masses, when converted to equivalent total atomic and molecular hydrogen, probe gas masses as low as $\sim 10^{8.5}\,M_\odot$. Our measurements delineate a bimodal distribution on the $M_{\rm gas}-M_*$ plane: the gas-rich, star-forming galaxies that trace the well-studied gas mass main sequence, and passive galaxies that occupy a distinct, gas-poor regime. These two populations, in turn, map into a bimodal distribution on the relation between SFR and gas mass surface density. Among low-redshift galaxies, galaxy mergers, including those that involve gas-rich and nearly equal-mass galaxies, exert a minimal impact on their SFR, specific SFR, or star formation efficiency. Starbursts are rare. The star formation efficiency of gas-rich, minor mergers even appears suppressed. This study stresses the multiple, complex factors that influence the evolution of the gas and its ability to form stars in mergers.

Rocco D'Agostino, Rafael C. Nunes

12 pages, 6 figures. Matches the version published in PRD

Measurements of the $H_0$ and $\sigma_8$ parameters within the standard cosmological model recently highlighted significant statistical tensions between the cosmic microwave background and low-redshift probes, such as local distance ladder, weak lensing and galaxy clustering surveys. In this work, we frame geometrical distances in a model-independent way by means of cosmographic approximations in the range $z \in (0, 2.3)$ to take into account a robust dataset composed of Baryon Acoustic Oscillations (BAO), type Ia Supernovae (SN), Cosmic Chronometer (CC) data, and measurements from Redshift Space Distortions (RSD). From the joint analysis BAO\,+\,SN\,+\,CC\,+\,RSD, we find an accuracy of $\sim$1.4\% and $\sim$3.7\% on $H_0$ and $\sigma_8$, respectively. Our result for $H_0$ is at 2$\sigma$ tension with local measurements by the SH0ES team, while our $\sigma_8$ estimate is at 2.6$\sigma$ tension with Planck-CMB analysis. This inference shows a tension statistically smaller when compared to those estimated via the $\Lambda$CDM model. We also find that the jerk parameter can deviate more than 3$\sigma$ from the $\Lambda$CDM prediction. Under the same cosmographic setup, we also present results by considering a SH0ES Gaussian prior on $H_0$ that allows for improved accuracy of the parameter space of the models. The present work brings observational constraints on $H_0$ and $\sigma_8$ into a new model-independent perspective, which differs from the predictions obtained within the $\Lambda$CDM paradigm.

The Fred Young Submillimeter Telescope (FYST) and the Simons Observatory Large Aperture Telescope (SO\ LAT) will deliver unprecedented high-resolution measurements of microwave sky emissions. Notably, one of those microwave sky emissions, the thermal Sunyaev-Zeldovich (tSZ) signal, is an essential probe for cluster astrophysics and cosmology. However, an obstacle to its measurement is contamination by the cosmic infrared background (CIB), especially at high frequencies. Our goal is to assess the detection and purity of tSZ power spectrum measurements from these two telescopes. We demonstrate that FYST's high-frequency coverage helps lower CIB contamination and improves signal detection. We simulated the various components of the microwave sky at the frequencies, sensitivities, and beam sizes of the upcoming SO LAT and FYST telescopes using full-sky Hierarchical Equal Area isoLatitude Pixelisation (HEALPix) map templates from the Websky simulations and the Python Sky Model (PySM). We used a map-based internal linear combination (ILC) and a constrained ILC (CILC) to extract the tSZ signal and compute residual noises to assess CIB contamination and signal recovery. We find that the CIB's residual noise power spectrum in the ILC-recovered tSZ is lowered by $\sim 35\%$ on average over the scales $\ell \in [500,5000]$ when SO LAT and FYST are combined compared to when SO LAT is used alone. We find that when using CILC to deproject CIB, the combined abilities of SO LAT and FYST offer a large $\ell \in [1800,3500]$ window in which the recovered tSZ power spectrum is not noise dominated.

Y.X. Jane Yap, Albert Kong, Kwan Lok Li

10 pages, 7 figures, 4 tables (table 4 included as ancillary file), accepted for publication in ApJ

We report optical observations of the millisecond pulsar binary system PSR J1622-0315 with the Lulin 1m telescope in Taiwan and the Lijiang 2.4m telescope in China between 2019 and 2021. The companion of the pulsar, which is of V~19 mag, showed ellipsoidal-distorted orbital variations in its light curves. The best-fit model to the light curves, with the binary code PHOEBE, gives a companion mass of 0.122+/-0.006 M_sun. This places PSR J1622-0315 in the spider-system subclass. We compared the properties of PSR J1622-0315 with other spider pulsar binaries for the scalings between the spin-down luminosity derived for the pulsar, irradiation luminosity of the companion, and X-ray luminosity of the binary. We find that pulsar irradiation in PSR J1622-0315 is insignificant and the irradiation luminosity of the transitional millisecond pulsars PSR J1023+0038 and PSR J1227-4853 are the highest among the redback systems.

The vast majority of Pulsar Wind Nebulae (PWNe) present in the Galaxy is formed by middle-aged systems characterized by a strong interaction of the PWN itself with the supernova remnant (SNR). Unfortunately, modelling these systems can be quite complex and numerically expensive, due to the non-linearity of the PWN-SNR evolution even in the simple 1D / one-zone case when the reverse shock of the SNR reaches the PWN, and the two begin to interact (and reverberation starts). Here we introduce a new numerical technique that couples the numerical efficiency of the one-zone thin shell approach with the reliability of a full ``lagrangian'' evolution, able to correctly reproduce the PWN-SNR interaction during the reverberation and to consistently evolve the particle spectrum beyond. Based on our previous findings, we show that our novel strategy resolves many of the uncertainties present in previous approaches, as the arbitrariness in the SNR structure, and ensure a robust evolution, compatible with results that can be obtained with more complex 1D dynamical approaches. Our approach enable us for the first time to provide reliable spectral models of the later compression phases in the evolution of PWNe. While in general we found that the compression is less extreme than that obtained without such detailed dynamical considerations, leading to the formation of less structured spectral energy distributions, we still find that a non negligible fraction of PWNe might experience a super-efficient phase, with the optical and/or X-ray luminosity exceeding the spin-down one.

Tinne Pauwels, Maddalena Reggiani, Hugues Sana, Alan Rainot, Kaitlin Kratter

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

One of the most remarkable properties of massive stars is that almost all of them are found in binaries or higher-order multiple systems. Observations that cover the full companion mass ratio and separation regime are essential to constrain massive star and binary formation theories. We used VLT/SPHERE to characterise the multiplicity properties of 20 OB stars in the active star-forming region Sco OB1. We simultaneously observed with the IFS and IRDIS instruments, obtaining high-contrast imaging observations that cover a field of view of 1".73 x 1".73 in YJH bands and 11" x 12".5 in $K_1$ and $K_2$ bands, respectively, corresponding to a separation range between $\sim$200 and 9000 AU. The observations reach contrast magnitudes down to $\Delta K_1 \sim 13$, allowing us to detect companions at the stellar-substellar boundary. In total, we detect 789 sources, most of which are likely background or foreground objects. We obtain SPHERE companion fractions of $2.3 \pm 0.4$ and $4.2 \pm 0.8$ for O- and B-type stars, respectively. Including all previously detected companions, we find a total multiplicity fraction of $0.89\pm0.07$ for our sample in the range of $\sim$0-12000 AU. In conclusion, SPHERE explores an as of yet uncharted territory of companions around massive stars, which is crucial to ultimately improve our understanding of massive star and binary formation.

Sandro Mereghetti, Michela Rigoselli, Ruben Salvaterra, Andrea Tiengo, Dominik Pacholski

Submitted to The Astrophysical Journal

230307A is the second brightest gamma ray burst detected in more than 50 years of observations and is located in the direction of the Magellanic Bridge. Despite its long duration, it is most likely the result of the compact merger of a binary ejected from a galaxy in the local universe (redshift z=0.065). Our XMM-Newton observation of its afterglow at 4.5 days shows a power-law spectrum with photon index $\Gamma =1.73 \pm0.10$, unabsorbed flux $F_{0.3-10\,\rm keV}=(8.8\pm0.5)\times 10^{-14}$ erg cm$^{-2}$ s$^{-1}$ and no absorption in excess of that produced in our Galaxy and in the Magellanic Bridge. We derive a limit of $N_{\rm H}^{\rm HOST} < 5\times 10^{20}$ cm$^{-2}$ on the absorption at the GRB redshift, which is a factor $\sim\,$5 below the value measured during the prompt phase. We searched for the presence of dust scattering rings with negative results and set an upper limit of the order of $A_V<0.05$ on the absorption from dust in the Magellanic Bridge.

In support of studies of decadal-timescale evolution of outer solar system atmospheres and ring systems, we present detailed Earth-based stellar occultation predictions for Jupiter, Saturn, Uranus, Neptune, Titan, and Triton for 2023-2050, based on the Gaia DR3 star catalog and near-IR K-band photometry from the 2MASS catalog. We tabulate the number of observable events by year and magnitude interval, reflecting the highly variable frequency of high-SNR events depending on the target's path relative to the star-rich regions of the Milky Way. We identify regions on Earth where each event is potentially observable, and for atmospheric occultations, we determine the latitude of the ingress and egress events. For Saturn, Uranus, and Neptune, we also compute the predicted ring occultation event times. We present representative subsets of the predicted events and highlights particularly promising events. Jupiter occultations with K $\leq$7 occur at a cadence of about one per year, with bright events at higher frequency in 2031 and 2043. Saturn occultations are much rarer, with only two predicted events with K $\leq$5 in 2032 and 2047. Ten Uranus ring occultations are predicted with K$\leq$10 for the period 2023 to 2050. Neptune traverses star-poor regions of the sky until 2068, resulting in only 13 predicted occultations for K$\leq$12 between 2023 and 2050. Titan has several high-SNR events between 2029--2031, whereas Triton is limited to a total of 22 occultations with K$\leq$15 between 2023 and 2050. Details of all predicted events are included in the Supplementary Online Material.

Florent G. Mertens, Jérôme Bobin, Isabella P. Carucci

14 pages, 9 figues, submitted to MNRAS, any comments welcome

Direct detection of the Cosmic Dawn and Epoch of Reionization via the redshifted 21-cm line of neutral Hydrogen will have unprecedented implications for studying structure formation in the early Universe. This exciting goal is challenged by the difficulty of extracting the faint 21-cm signal buried beneath bright astrophysical foregrounds and contaminated by numerous systematics. Here, we focus on improving the Gaussian Process Regression (GPR) signal separation method originally developed for LOFAR observations. We address a key limitation of the current approach by incorporating covariance prior models learnt from 21-cm signal simulations using Variational Autoencoder (VAE) and Interpolatory Autoencoder (IAE). Extensive tests are conducted to evaluate GPR, VAE-GPR, and IAE-GPR in different scenarios. Our findings reveal that the new method outperforms standard GPR in component separation tasks. Moreover, the improved method demonstrates robustness when applied to signals not represented in the training set. It also presents a certain degree of resilience to data systematics, highlighting its ability to effectively mitigate their impact on the signal recovery process. However, our findings also underscore the importance of accurately characterizing and understanding these systematics to achieve successful detection. Our generative approaches provide good results even with limited training data, offering a valuable advantage when a large training set is not feasible. Comparing the two algorithms, IAE-GPR shows slightly higher fidelity in recovering power spectra compared to VAE-GPR. These advancements highlight the strength of generative approaches and optimise the analysis techniques for future 21-cm signal detection at high redshifts.

Strong evidence for a gravitational-wave background (GWB) has been reported in the nano-Hertz band. Interpreting the origin of this background to be scalar-induced gravitational waves (SIGWs), we explore the equation of state (EoS) of the early universe by performing Bayes parameter inferences across the big-bang nucleosynthesis (BBN), cosmic microwave background (CMB), and pulsar timing array (PTA) joint observations for the first time. Assuming a monochromatic power spectrum for primordial curvature perturbations, we obtain the spectral amplitude $A\sim10^{-3}-10^{-1}$ and spectral peak frequency $f_\ast\sim10^{-7}-10^{-6}$ Hz. We find that the radiation domination with EoS $w=1/3$ is compatible with the current observational data, the kination domination with EoS $w=1$ is not forbidden, while the early matter domination with EoS $w=0$ is excluded at more than $2\sigma$ confidence level. These results can be tested with future observations.

M. Razzano, A. Fiori, P. M. Saz Parkinson, R. P. Mignani, A. De Luca, A. K. Harding, M. Kerr, M. Marelli, V. Testa

9 pages, 7 figures, 1 Table. Accepted for publication in Astronomy and Astrophysics (A&A)

Context. We have investigated the multiwavelength emission of PSR J2021+4026, the only isolated gamma-ray pulsar known to be variable, which in October 2011 underwent a simultaneous change in gamma-ray flux and spin-down rate, followed by a second mode change in February 2018. Multiwavelength monitoring is crucial to understand the physics behind these events and how they may have affected the structure of the magnetosphere. Aims.The monitoring of pulse profile alignment is a powerful diagnostic tool for constraining magnetospheric reconfiguration. We aim to investigate timing or flux changes related to the variability of PSR J2021+4026 via multiwavelength observations, including gamma-ray observations from Fermi-LAT, X-ray observations from XMM-Newton, and a deep optical observation with the Gran Telescopio Canarias.Methods. We performed a detailed comparison of the timing features of the pulsar in gamma and X-rays and searched for any change in phase lag between the phaseogram peaks in these two energy bands. Although previous observations did not detect a counterpart in visible light, we also searched for optical emission that might have increased due to the mode change, making this pulsar detectable in the optical. Results.We have found a change in the gamma-to X-ray pulse profile alignment by 0.21$\pm$0.02 in phase, which indicates that the first mode change affected different regions of the pulsar magnetosphere. No optical counterpart was detected down to g'=26.1 and r'=25.3. Conclusions.We suggest that the observed phase shift could be related to a reconfiguration of the connection between the quadrupole magnetic field near the stellar surface and the dipole field that dominates at larger distances. This is consistent with the picture of X-ray emission coming from the heated polar cap and with the simultaneous flux and frequency derivative change observed during the mode changes.

Ali Pourmand, Natalia Ivanova

18 pages. Published in ApJ 2023,952, 126. Data and code are available at this https URL

Evolutionary calculations for stars in close binary systems are in high demand to obtain better constraints on gravitational wave source progenitors, understand transient events from stellar interactions, and more. Modern one-dimensional stellar codes make use of the Roche lobe radius $R_{\rm L}$ concept in order to treat stars in binary systems. If the stellar companion is approaching its $R_{\rm L}$, mass transfer treatment is initiated. However, the effective acceleration also affects the evolution of a star in a close binary system. This is different from the gravity inside a single star, whether that single star is rotating or not. Here, we present numerically obtained tables of properties of stars in a binary system as a function of the effective potential: volume-equivalent radii of the equipotential surfaces, effective accelerations and the inverse effective accelerations averaged over the same equipotential surfaces, and the properties of the L1 plane cross-sections. The tables are obtained for binaries where the ratios of the primary star mass to the companion star mass are from $10^{-6}$ to $10^5$ and include equipotential surfaces up to the star's outer Lagrangian point. We describe the numerical methods used to obtain these quantities and report how we verified the numerical results. We also describe and verify the method to obtain the effective acceleration for non-point mass distributions. We supply a sample code showing how to use our tables to get the average effective accelerations in one-dimensional stellar codes.

Blagoy Rangelov, Hui Yang, Brice Williams, Oleg Kargaltsev, Jeremy Hare, Kean Martinic

16 pages, 4 figures, 5 tables. Submitted to ApJ on 6th July 2023

In the latest data release from the Fermi $\gamma$-Ray Space Telescope (the 4th Fermi LAT 12-year Catalog or 4FGL) more than 50% of the Galactic sources are yet to be identified. We observed thirteen unidentified Fermi LAT sources with Chandra X-Ray Observatory (CXO) to explore their nature. We report the results of the classification of X-ray sources in the fields of these $\gamma$-ray sources and discuss the implications for their nature. We use the multiwavelength (MW) data for machine-learning classification accompanied by a more detailed spectral/variability analysis for brighter sources. Seven 4FGL sources have $\gamma$-ray pulsars within their position error ellipses. Three of these pulsars are either detected in the CXO images or show hints of X-ray emission. Within the positional uncertainties of three 4FGL sources we detect X-ray sources that may be yet unknown pulsars, depending on the MW association. In addition to point sources, we discovered 2 extended sources one of which is likely to be a bowshock pulsar-wind nebula associated with PSR J1358.3-6026. Finally, we classify other X-ray sources detected in these observations and report most interesting classifications.

Angelos Vourlidas, Amir Caspi, Yuan-Kuen Ko, J. Martin Laming, James P. Mason, Mari Paz Miralles, Nour-Eddine Raouafi, John C. Raymond, Daniel B. Seaton, Leonard Strachan, Nicholeen Viall, Juliana Vievering, Matthew J. West

White paper submitted to the Decadal Survey for Solar and Space Physics (Heliophysics) 2024-2033; 6 pages, 1 figure, 2 tables

Our current theoretical and observational understanding suggests that critical properties of the solar wind and Coronal Mass Ejections (CMEs) are imparted within 10 Rs, particularly below 4 Rs. This seemingly narrow spatial region encompasses the transition of coronal plasma processes through the entire range of physical regimes from fluid to kinetic, and from primarily closed to open magnetic field structures. From a physics perspective, therefore, it is more appropriate to refer to this region as the Critical Coronal Transition Region (CCTR) to emphasize its physical, rather than spatial, importance to key Heliophysics science. This white paper argues that the comprehensive exploration of the CCTR will answer two of the most central Heliophysics questions, "How and where does the solar wind form?" and "How do eruptions form?", by unifying hardware/software/modeling development and seemingly disparate research communities and frameworks. We describe the outlines of decadal-scale plan to achieve that by 2050.

Sina Chen, Ari Laor, Ehud Behar, Ranieri D. Baldi, Joseph D. Gelfand

Accepted for publication in MNRAS

The origin of the radio emission in radio-quiet quasars (RQQ) is not established yet. We present new VLBA observations at 1.6 and 4.9 GHz of ten RQQ (nine detected), which together with published earlier observations of eight RQQ (five detected), forms a representative sample of 18 RQQ drawn from the Palomar-Green sample of low z (< 0.5) AGN. The spectral slope of the integrated emission extends from very steep (alpha < -1.98) to strongly inverted (alpha = +2.18), and the slopes of nine of the 14 objects are flat (alpha > -0.5). Most objects have an unresolved flat-spectrum core, which coincides with the optical Gaia position. The extended emission is generally steep-spectrum, has a low brightness temperature (< 10^7 K), and is displaced from the optical core (the Gaia position) by ~ 5-100 pc. The VLBA core flux is tightly correlated with the X-ray flux, and follows a radio to X-ray luminosity relation of log L_R/L_X = -6, for all objects with a black hole mass log M_BH/M_Sun < 8.5. The flatness of the core emission implies a compact source size (< 0.1 pc), which likely originates from the accretion disk corona. The mas-scale extended emission is optically thin and of clumpy structure, and is likely produced by an outflow from the center. Radio observations at higher frequencies can further test the accretion disk coronal emission interpretation for the core emission in RQQ.

Fabo Feng, R. Paul Butler, Steven S. Vogt, Jennifer Burt, Bradford Holden, Yicheng Rui

14 pages, 5 figures, 4 tables, accepted for publication in MNRAS

With its near-to-mid-infrared high contrast imaging capabilities, JWST is ushering us into a golden age of directly imaging Jupiter-like planets. As the two closest cold Jupiters, $\varepsilon$ Ind A b and $\varepsilon$ Eridani b have sufficiently wide orbits and adequate infrared emissions to be detected by JWST. To detect more Jupiter-like planets for direct imaging, we develop a GOST-based method to analyze radial velocity data and multiple Gaia data releases simultaneously. Without approximating instantaneous astrometry by catalog astrometry, this approach enables the use of multiple Gaia data releases for detection of both short-period and long-period planets. We determine a mass of $2.96_{-0.38}^{+0.41}$ $M_{\rm Jup}$ and a period of $42.92_{-4.09}^{+6.38}$ yr for $\varepsilon$ Ind A b. We also find a mass of $0.76_{-0.11}^{+0.14}$ $M_{\rm Jup}$, a period of $7.36_{-0.05}^{+0.04}$ yr, and an eccentricity of 0.26$_{-0.04}^{+0.04}$ for $\varepsilon$ Eridani b. The eccentricity differs from that given by some previous solutions probably due to the sensitivity of orbital eccentricity to noise modeling. Our work refines the constraints on orbits and masses of the two nearest Jupiters and demonstrate the feasibility of using multiple Gaia data releases to constrain Jupiter-like planets.

S. Torres, P. Cruz, R. Murillo-Ojeda, F. M. Jiménez-Esteban, A. Rebassa-Mansergas, E. Solano, M. E. Camisassa, R. Raddi, J. Doliguez Le Lourec

7 pages, 11 pages. Accepted for publication in A&A

The characterization of white dwarf atmospheres is crucial for accurately deriving stellar parameters such as effective temperature, mass, and age. We aim to classify the population of white dwarfs up to 500 pc into hydrogen-rich or hydrogen-deficient atmospheres based on Gaia spectra and to derive an accurate spectral type-temperature distribution of white dwarfs as a function of the effective temperature for the largest observed unbiased sample of these objects. We took advantage of the recent Gaia low-resolution spectra available for 76,657 white dwarfs up to 500 pc. We calculated synthetic J-PAS narrow-band photometry and fitted the spectral energy distribution of each object with up-to-date models for hydrogen-rich and helium-rich white dwarf atmospheres. We estimated the probability for a white dwarf to have a hydrogen-rich atmosphere and validated the results using the Montreal White Dwarf Database. Finally, precise effective temperature values were derived for each object using La Plata evolutionary models. We have successfully classified a total of 65,310 white into DAs and non-DAs with an accuracy of 94%. An unbiased subsample of nearly 34,000 objects was built, from which we computed a precise spectral distribution spanning an effective temperature range from 5,500 to 40,000 K, while accounting for potential selection effects. Some characteristic features of the spectral evolution, such as the deficit of helium-rich stars at T_eff $\approx$35,000-40,000 K and in the range 22,000 < T_eff < 25,000 K, as well as a gradual increase from 18,000K to T_eff $\approx$7,000K, where the non-DA stars percentage reaches its maximum of 41%, followed by a decrease for cooler temperatures, are statistically significant. These findings will provide precise constraints for the proposed models of spectral evolution.

Bea Zenteno-Quinteros, Pablo S Moya, Marian Lazar, Adolfo F Vinas, Stefaan Poedts

18 pages, 8 figures, accepted for publication in The Astrophysical Journal

Temperature anisotropy and field-aligned skewness are commonly observed non-thermal features in electron velocity distributions in the solar wind. These characteristics can act as a source of free energy to destabilize different electromagnetic wave modes, which may alter the plasma state through wave-particle interactions. Previous theoretical studies have mainly focused on analyzing these non-thermal features and self-generated instabilities individually. However, to obtain a more accurate and realistic understanding of kinetic processes in the solar wind, it is necessary to examine the interplay between these two energy sources. By means of linear kinetic theory, in this paper we investigate the excitation of the parallel-propagating whistler mode, when it is destabilized by electron populations exhibiting both temperature anisotropy and field-aligned strahl or skewness. To describe the solar wind electrons, we adopt the Core-Strahlo model as an alternative approach. This model offers the advantage of representing the suprathermal features of halo and strahl electrons, using a single skew-Kappa distribution already known as the strahlo population. Our findings show that when the electron strahlo exhibits an intrinsic temperature anisotropy, this suprathermal population becomes a stronger and more efficient source of free energy for destabilizing the whistler mode. This suggests a greater involvement of the anisotropic strahlo in processes conditioned by wave-particle interactions. Present results also suggest that the contribution of core anisotropy can be safely disregarded when assessing the importance of instabilities driven by the suprathermal population. This allows for a focused study, particularly regarding the regulation of electron heat flux in the solar wind.

We discuss the non-linear corrections entering in the calculation of the primordial black hole abundance from the non-linear radiation transfer function and the determination of the true physical horizon crossing. We show that the current standard techniques to calculate the abundance of primordial black holes suffer from uncertainties and argue that the primordial black hole abundance may be much smaller than what routinely considered. This would imply, among other consequences, that the interpretation of the recent pulsar timing arrays data from scalar-induced gravitational waves may not be ruled out because of an overproduction of primordial black holes.

Suvashis Maity, H. V. Ragavendra, Shiv K. Sethi, L. Sriramkumar

34 pages, 8 figures

[Abridged] In contemporary literature, the calculation of modifications to the inflationary scalar power spectrum due to the loops from the higher order interaction terms in the Hamiltonian have led to a discussion regarding the validity of perturbation theory. Recently, there have been efforts to examine the contributions to the scalar power spectrum due to the loops arising from the cubic order terms in the action describing the perturbations, specifically in inflationary scenarios that permit an epoch of ultra slow roll (USR). A phase of USR inflation leads to significant observational consequences, such as the copious production of primordial black holes. In this work, we study the loop contributions to the scalar power spectrum in a scenario of USR inflation arising due to the quartic order terms in the action describing the scalar perturbations. We compute the loop contributions to the scalar power spectrum due to the dominant term in the action at the quartic order. We consider a scenario wherein a phase of USR is sandwiched between two stages of slow roll inflation and analyze the behavior of the loop contributions in terms of the parameters involved. We examine the late, intermediate and early epochs of USR during inflation. In the inflationary scenario involving a late phase of USR, for reasonable choices of the parameters, we show that the loop corrections are negligible for the entire range of wave numbers. In the intermediate case, the contributions from the loops prove to be scale invariant over large scales, and we find that these contributions can amount to 30% of the leading order power spectrum. In the case wherein USR sets in early, we find that the loop contributions could be negative and can dominate the power spectrum at the leading order, which indicates a breakdown of the perturbative expansion. We conclude with a brief summary and outlook.

João R. T. de Mello Neto

Presented at the 38th International Cosmic Ray Conference (ICRC 2023). 8 pages, 3 figures

Ultra-high-energy cosmic neutrinos (UHE), with energies above 100 PeV, are unparalleled probes of the most energetic astrophysical sources and weak interactions at energies beyond the reach of accelerators. GRAND is an envisioned observatory of UHE particles - neutrinos, cosmic rays, and gamma rays - consisting of 200,000 radio antennas deployed in sub-arrays at different locations worldwide. GRAND aims to detect the radio emission from air showers induced by UHE particle interactions in the atmosphere and underground. For neutrinos, it aims to reach a flux sensitivity of $\sim 10^{-10}$ GeV cm$^{-2}$ s$^{-1}$ sr$^{-1}$, with a sub-degree angular resolution, which would allow it to test the smallest predicted diffuse fluxes of UHE neutrinos and to discover point sources. The GRAND Collaboration operates three prototype detector arrays simultaneously: GRAND@Nan\c{c}ay in France, GRANDProto300 in China, and GRAND@Auger in Argentina. The primary purpose of GRAND@Nan\c cay is to serve as a testbench for hardware and triggering systems. On the other hand, GRANDProto300 and GRAND@Auger are exploratory projects that pave the way for future stages of GRAND. GRANDProto300 is being built to demonstrate autonomous radio-detection of inclined air showers and study cosmic rays near the proposed transition between galactic and extragalactic sources. All three arrays are in the commissioning stages. It is expected that by 2028, the detector units of the final design could be produced and deployed, marking the establishment of two GRAND10k arrays in the Northern and Southern hemispheres. We will survey preliminary designs, simulation results, construction plans, and the extensive research program made possible by GRAND.

K. Böckmann, M. Brüggen, B. Koribalski, A. Veronica, T.H. Reiprich, E. Bulbul, Y. E. Bahar, F. Balzer, J. Comparat, C. Garrel, V. Ghirardini, G. Gürkan, M. Kluge, D. Leahy, A. Merloni, A. Liu, M. E. Ramos-Ceja, M. Salvato, J. Sanders, S. Shabala, X. Zhang

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

The extended ROentgen Survey with an Imaging Telescope Array (eROSITA) telescope onboard the Spectrum-Roentgen-Gamma (SRG) mission has finished the first eROSITA All-Sky Survey (eRASS:1), and detected 10$^4$ galaxy clusters in the western Galactic hemisphere. In the radio band, the Australian Square Kilometre Array Pathfinder (ASKAP) telescope finished its pilot 1 phase of the project 'Evolutionary Map of the Universe' (EMU) with 220.000 sources in a 270 deg$^2$ field overlapping with eRASS:1. These two surveys are used to study radio-mode Active Galactic Nuclei (AGN) in clusters. In order to understand the efficiency of radio-mode feedback at the centers of galaxy clusters, we relate the radio properties of brightest cluster galaxies (BCG) to the X-ray properties of the host clusters. We identify the central radio sources in eRASS:1 clusters or calculate corresponding upper limits on the radio luminosity. Then, we derive relations between the X-ray properties of the clusters and the radio properties of the corresponding central radio source. We also apply a mid-infrared color criterion using WISE colors to identify AGN. In total we investigate a sample of 75 clusters. We find a statistically significant between the X-ray luminosity of the cluster and the 944 MHz radio luminosity of the corresponding central radio galaxy. There is also a positive trend between the radio power and the largest linear size (LLS) of the radio source. The density and the LLS do not show any correlation. We find that in high luminosity clusters with L_X > $10^{43}$ erg s$^{-1}$ the kinetic luminosity of the radio jets is not longer correlated with the X-ray luminosity and discuss various reasons. We find an anti-correlation between the central cooling time t_cool and the radio luminosity L_R indicating a need for more powerful AGN in clusters with short central cooling times.

P. Short, A. Lawrence, M. Nicholl, M. Ward, T. M. Reynolds, S. Mattila, C. Yin, I. Arcavi, A. Carnall, P. Charalampopoulos, M. Gromadzki, P. G. Jonker, S. Kim, G. Leloudas, I. Mandel, F. Onori, M. Pursiainen, S. Schulze, C. Villforth, T. Wevers

Tidal disruption events (TDEs) occur when a star gets torn apart by a supermassive black hole as it crosses its tidal radius. We present late-time optical and X-ray observations of the nuclear transient AT2019qiz, which showed the typical signs of an optical-UV transient class commonly believed to be TDEs. Optical spectra were obtained 428, 481 and 828 rest-frame days after optical lightcurve peak, and a UV/X-ray observation coincided with the later spectrum. The optical spectra show strong coronal emission lines, including [Fe VII], [Fe X], [Fe XI] and [Fe XIV]. The Fe lines rise and then fall, except [Fe XIV] which appears late and rises. We observe increasing flux of narrow H-alpha and H-beta and a decrease in broad H-alpha flux. The coronal lines have FWHMs ranging from ~150 - 300km/s, suggesting they originate from a region between the broad and narrow line emitting gas. Between the optical flare and late-time observation, the X-ray spectrum softens dramatically. The 0.3-1 keV X-ray flux increases by a factor of ~50 while the hard X-ray flux decreases by a factor of ~6. WISE fluxes also rose over the same period, indicating the presence of an infrared echo. With AT2017gge, AT2019qiz is one of two examples of a spectroscopically-confirmed optical-UV TDE showing delayed coronal line emission, supporting speculations that Extreme Coronal Line Emitters in quiescent galaxies can be echos of unobserved past TDEs. We argue that the coronal lines, narrow lines, and infrared emission arise from the illumination of pre-existing material likely related to either a previous TDE or AGN activity.

Alberto Manuel Martínez-García, Andrés del Pino, Ewa L. Łokas, Roeland P. van der Marel, Antonio Aparicio

12 pages, 9 figures, submitted to MNRAS

We present a kinematic study of a thousand of dwarf satellites of MW/M31-like hosts from the IllustrisTNG50 simulation. Internal kinematics were derived for all the snapshots to obtain a historical record of their rotation velocity in the plane of the sky ($|V_T|$) and the amplitude of their velocity gradients along the line of sight ($A_{\rm grad}^{v_z}$) measured from the host. For the majority of the satellites we initially detected rotation in the plane of the sky (65%) or velocity gradients (80%), and this was progressively reduced to 45% and 68% at $z = 0$ respectively. We find that the evolution of the rotation in the plane of the sky and the velocity gradients differs according to type of dwarfs, which could be explained in terms of their different masses and orbital histories. We observe that interaction with the host has an impact on the evolution of the internal kinematics of the satellites. The rotation signal of the satellites is progressively reduced during pericentric passages, the first pericentre being especially disruptive for the initial kinematics. We observe temporary increases in $A_{\rm grad}^{v_z}$ during pericentric passage caused by tidal interaction with the host, $A_{\rm grad}^{v_z}$ increasing as the satellites approach their pericentre and dropping as they move away. In summary, we conclude that the presence of detectable rotation in dwarf satellites is not uncommon, and that the evolution of their internal kinematics is clearly affected by their interaction with the host.

H. Qu, M. Sako, M. Vincenzi, C. Sanchez, D. Brout, R. Kessler, R. Chen, L. Galbany, L. Kelsey, J. Lee, C. Lidman, B. Popovic, B. Rose, M. Smith, M. Sullivan, P. Wiseman, T. M. C. Abbott, M. Aguena, O. Alves, D. Bacon, E. Bertin, D. Brooks, D. L. Burke, A. Carnero Rosell, J. Carretero, L. N. da Costa, M. E. S. Pereira, H. T. Diehl, P. Doel, S. Everett, I. Ferrero, J. Frieman, J. Garcia-Bellido, G. Giannini, D. Gruen, R. A. Gruendl, G. Gutierrez, S. R. Hinton, D. L. Hollowood, K. Honscheid, D. J. James, K. Kuehn, O. Lahav, J. L. Marshall, J. Mena-Fernandez, F. Menanteau, R. Miquel, R. L. C. Ogando, A. Palmese, A. Pieres, A. A. Plazas Malagon, M. Raveri, E. Sanchez, I. Sevilla-Noarbe, M. Soares-Santos, E. Suchyta, G. Tarle, N. Weaverdyck, DES Collaboration

Submitted to ApJ

Redshift measurements, primarily obtained from host galaxies, are essential for inferring cosmological parameters from type Ia supernovae (SNe Ia). Matching SNe to host galaxies using images is non-trivial, resulting in a subset of SNe with mismatched hosts and thus incorrect redshifts. We evaluate the host galaxy mismatch rate and resulting biases on cosmological parameters from simulations modeled after the Dark Energy Survey 5-Year (DES-SN5YR) photometric sample. For both DES-SN5YR data and simulations, we employ the directional light radius method for host galaxy matching. In our SN Ia simulations, we find that 1.7% of SNe are matched to the wrong host galaxy, with redshift difference between the true and matched host of up to 0.6. Using our analysis pipeline, we determine the shift in the dark energy equation of state parameter (Dw) due to including SNe with incorrect host galaxy matches. For SN Ia-only simulations, we find Dw = 0.0013 +/- 0.0026 with constraints from the cosmic microwave background (CMB). Including core-collapse SNe and peculiar SNe Ia in the simulation, we find that Dw ranges from 0.0009 to 0.0032 depending on the photometric classifier used. This bias is an order of magnitude smaller than the expected total uncertainty on w from the DES-SN5YR sample of around 0.03. We conclude that the bias on w from host galaxy mismatch is much smaller than the uncertainties expected from the DES-SN5YR sample, but we encourage further studies to reduce this bias through better host-matching algorithms or selection cuts.

This study explores the implications of Dark Matter in Neutron Stars (DMANS) by focusing on two specific astronomical objects: HESS J1731-347 and PSR J0952-0607. Varying the Fermi momentum k$f^{\rm DM}$ of DM, the study analyzes the EOS for the INRS model with and without DM. Results show the robustness of the model, with most EOS curves within chiral Effective Field Theory bounds. Our model predicts a maximum mass of $2.343 \ M_\odot$ for PSR J0952-0607, satisfying NICER bounds. The analysis suggests HESS J1731-347 could be a DMANS. Constraints on DM within NSs are established, and tidal deformability lies within GW event limits. Nonradial $f$-mode oscillations increase with DM, concluding low mass stars pulsate at higher frequencies.

G. Angloher, S. Banik, G. Benato, A. Bento, A. Bertolini, R. Breier, C. Bucci, J. Burkhart, L. Canonica, A. D'Addabbo, S. Di Lorenzo, L. Einfalt, A. Erb, F. v. Feilitzsch, S. Fichtinger, D. Fuchs, A. Garai, V.M. Ghete, P. Gorla, P.V. Guillaumon, S. Gupta, D. Hauff, M. Jeskovsky, J. Jochum, M. Kaznacheeva, A. Kinast, H. Kluck, H. Kraus, S. Kuckuk, A. Langenkaemper, M. Mancuso, L. Marini, L. Meyer, V. Mokina, A. Nilima, M. Olmi, T. Ortmann, C. Pagliarone, L. Pattavina, F. Petricca, W. Potzel, P. Povinec, F. Proebst, F. Pucci, F. Reindl, J. Rothe, K. Schaeffner, J. Schieck, D. Schmiedmayer, S. Schoenert, C. Schwertner, M. Stahlberg, L. Stodolsky, C. Strandhagen, R. Strauss, I. Usherov, F. Wagner, M. Willers, V. Zema (CRESST Collaboration), F. Ferella, M. Laubenstein, S. Nisi

24 pages, 14 figures, submitted to EPJC

Using CaWO$_4$ crystals as cryogenic calorimeters, the CRESST experiment searches for nuclear recoils caused by the scattering of potential Dark Matter particles. A reliable identification of a potential signal crucially depends on an accurate background model. In this work we introduce an improved normalisation method for CRESST's model of the electromagnetic backgrounds. Spectral templates, based on Geant4 simulations, are normalised via a Bayesian likelihood fit to experimental background data. Contrary to our previous work, no assumption of partial secular equilibrium is required, which results in a more robust and versatile applicability. Furthermore, considering the correlation between all background components allows us to explain 82.7% of the experimental background within [1 keV, 40 keV], an improvement of 18.6% compared to our previous method.

Geraint Pratten, Patricia Schmidt, Hannah Middleton, Alberto Vecchio

14+6 pages, 8 figures

The Laser Interferometer Space Antenna (LISA) will play a vital role in constraining the origin and evolution of massive black holes throughout the Universe. In this study we use a waveform model (IMRPhenomXPHM) that includes both precession and higher multipoles, and full Bayesian inference to explore the accuracy to which LISA can constrain the binary parameters. We demonstrate that LISA will be able to track the evolution of the spins -- magnitude and orientation -- to percent accuracy, providing crucial information on the dynamics and evolution of massive black hole binaries and the galactic environment in which the merger takes place. Such accurate spin-tracking further allows LISA to measure the recoil velocity of the remnant black hole to better than $100\,\mathrm{km}\,\mathrm{s}^{-1}$ (90\% credibility) and its direction to a few degrees, which provides additional important astrophysical information on the post-merger association. Using a systematic suite of binaries, we showcase that the component masses will be measurable at the sub-percent level, the sky area can be constrained to withing $\Delta \Omega_{90} \approx 0.01 \, \rm{deg}^2$, and the binary redshift to less than $0.01$.

Noah Franz, Mitchell Dennis, Jeremy Sakstein

12 pages, 5 figures. Dataset and code available at this https URL

We use a novel method to constrain the neutrino magnetic dipole moment ($\mu_{\nu}$) using the empirically-calibrated tip of the red giant branch I-band magnitude that fully accounts for uncertainties in stellar physics. Our method uses machine learning to emulate the results of stellar evolution codes. This reduces the I-Band magnitude computation time to milliseconds, which enables a Bayesian statistical analysis where $\mu_{\nu}$ is varied simultaneously with the stellar physics, allowing for a complete exploration of parameter space. We find the region $\mu_{\nu} \leq 6\times10^{-12}\mu_{\textrm{B}}$ (with $\mu_{\textrm{B}}$ the Bohr magneton), previously believed to be excluded, is unconstrained after accounting for degeneracies with stellar physics. It is likely that larger values are similarly unconstrained. We discuss the implications of our results for future neutrino magnetic dipole moment searches and for other astrophysical probes.

Salvatore Capozziello, Vittorio De Falco, Carmen Ferrara

7 pages, 1 figure

In the framework of metric-affine gravity, we consider the role of the boundary term in Symmetric Teleparallel Gravity assuming $f(Q,B)$ models where $f$ is a smooth function of the non-metricity scalar $Q$ and the related boundary term $B$. Starting from a variational approach, we derive the field equations and compare them with respect to those of $f(Q)$ gravity in the limit of $B\to0$. It is possible to show that $f(Q,B)=f(Q-B)$ models are dynamically equivalent to $f(R)$ gravity as in the case of teleparallel $f(\tilde{B}-T)$ gravity (where $B\neq \tilde{B}$). Furtherrmore, conservation laws are derived. In this perspective, considering boundary terms in $ f(Q)$ gravity represents the last ingredient towards the Extended Geometric Trinity of Gravity where $f(R)$, $f(T,\tilde{B})$ and $f(Q,B)$ can be dealt under the same standard.

We demonstrated that a weak magnetic field can increase the permittivity, leading to a reduction in the potential barrier within the Debye sphere consisting of electrons and a nucleus. By solving the Boltzmann equation with the inclusion of the magnetic field, we obtained the magnetized permittivity. The resulting enhanced permittivity field inversely decreases the potential barrier, thereby increasing the reaction rate between two fusing nuclei. We compared this Boltzmann kinetic approach with the Debye potential method. We found that they are qualitatively consistent. Further, we also derived the magnetized Debye potential composed of the conventional term with a new magnetic effect. Both approaches indicate that magnetized plasmas, which have existed since the Big Bang, have ultimately influenced permittivity, potential barrier, and nucleosynthesis.

Soichiro Morisaki, Rory Smith, Leo Tsukada, Surabhi Sachdev, Simon Stevenson, Colm Talbot, Aaron Zimmerman

19 pages, 9 figures

We present a rapid parameter estimation framework for compact binary coalescence (CBC) signals observed by the LIGO-Virgo-KAGRA (LVK) detector network. The goal of our framework is to enable optimal source localization of binary neutron star (BNS) signals in low latency, as well as improve the overall scalability of full CBC parameter estimation analyses. Our framework is based on the reduced order quadrature (ROQ) technique, and resolves its shortcomings by utilizing multiple ROQ bases in a single parameter estimation run. We have also developed sets of compact ROQ bases for various waveform models, IMRPhenomD, IMRPhenomPv2, IMRPhenomPv2$\_$NRTidalv2, and IMRPhenomXPHM. We benchmark our framework with hundreds of simulated observations of BNS signals by the LIGO-Virgo detector network, and demonstrate that it provides accurate and unbiased estimates on BNS source location, with a median analysis time of $6$ minutes. The median searched area is reduced by around 30$\%$ compared to estimates produced by BAYESTAR: from $21.8\,\mathrm{deg^2}$ to $16.6\,\mathrm{deg^2}$. Our framework also enables detailed parameter estimation taking into account gravitational-wave higher multipole moments, the tidal deformation of colliding objects, and detector calibration errors of amplitude and phase with the time scale of hours. Our rapid parameter estimation technique has been implemented in one of the LVK parameter estimation engines, BILBY, and is being employed by the automated parameter estimation analysis of the LVK alert system.

N. V. Nunes, N. Bartel, A. Belonenko, G. D. Manucharyan, S. M. Popov, V. N. Rudenko, L. I. Gurvits, G. Cimò, G. Molera Calvés, M. V. Zakhvatkin, M. F. Bietenholz

Nelson V Nunes et al 2023 Class. Quantum Grav. in press this https URL

The Einstein Equivalence Principle (EEP) is a cornerstone of general relativity and predicts the existence of gravitational redshift. We report on new results of measuring this shift with RadioAstron (RA), a space VLBI spacecraft launched into an evolving high eccentricity orbit around Earth with geocentric distances reaching 353,000 km. The spacecraft and ground tracking stations at Pushchino, Russia, and Green Bank, USA, were each equipped with a hydrogen maser frequency standard allowing a possible violation of the predicted gravitational redshift, in the form of a violation parameter $\varepsilon$, to be measured. By alternating between RadioAstron's frequency referencing modes during dedicated sessions between 2015 and 2017, the recorded downlink frequencies can essentially be corrected for the non-relativistic Doppler shift. We report on an analysis using the Doppler-tracking frequency measurements made during these sessions and find $\varepsilon = (2.1 \pm 3.3)\times10^{-4}$. We also discuss prospects for measuring $\varepsilon$ with a significantly smaller uncertainty using instead the time-domain recordings of the spacecraft signals and envision how $10^{-7}$ might be possible for a future space VLBI mission.

Lech Wiktor Piotrowski

Proceedings of the 38th International Cosmic Ray Conference (ICRC2023)

The vast amounts of data to be collected by the Giant Radio Array for Neutrino Detection (GRAND) and its prototype - GRANDProto300 - require the use of a data format very efficient in terms of i/o speed and compression. At the same time, the data should be easily accessible, without the knowledge of the intricacies of the format, both for bulk processing and for detailed event-by-event analysis and reconstruction. We present the format and the structure prepared for GRAND data, the concept of the data-processing chain, and data-oriented and analysis-oriented interfaces written in Python.

For the classical N-body problem, an approach is proposed based on the introduction of some natural in the physical sense optimization problems of mathematical programming for finding a conditional minimum for the characteristics of the system on the set of its possible states. The solution of these problems then makes it possible to construct families of flat stationary and periodic trajectories of the system and also to find relationships and estimates for the characteristics of the system on these trajectories. It is shown that when the system moves on a plane on trajectories generated by the global minimum in these optimization problems, at any time the minimum possible size of the system is achieved at each current level of its "cohesion" (or potential energy). Similar optimization problems are considered for finding a conditional minimum for the characteristics of a system in three-dimensional space. It is shown that the solution of these problems can be achieved only on flat trajectories of the system and is achieved, in particular, on the constructed flat stationary and periodic trajectories. In addition, it is shown that the trajectory of the system in three-dimensional space, at least at one point of which the minimum possible size of the system is achieved at the current value of its cohesion (or potential energy), can only be flat. And such trajectories are, in particular, flat stationary and periodic trajectories generated by the global minimum in the considered optimization problems.

S. Mahesh Chandran, Karthik Rajeev, S. Shankaranarayanan (IIT Bombay)

41 pages, 7 figures, 1 Table

Understanding the emergence of classical behavior from a quantum theory is vital to establishing the quantum origin for the temperature fluctuations observed in the Cosmic Microwave Background (CMB). We show that a real-space approach can comprehensively address the quantum-to-classical transition problem in the leading order of curvature perturbations. To this end, we test spatial bipartitions of quadratic systems for the interplay between three different signatures of classical behavior : i) decoherence, ii) peaking of the Wigner function about classical trajectories, and iii) relative suppression of non-commutativity in observables. We extract these signatures from the covariance matrix of a multi-mode Gaussian state and address them primarily in terms of entanglement entropy and log-classicality. Through a phase-space stability analysis of spatial sub-regions via their reduced Wigner function, we ascertain that the underlying cause for the dominance of classicality signatures is the occurrence of gapped inverted mode instabilities. While the choice of conjugate variables enhances some of these signatures, decoherence studied via entanglement entropy is the stronger and more reliable condition for classicality to emerge. We demonstrate the absence of decoherence, which preempts a quantum-to-classical transition of scalar fluctuations in an expanding background in $(1+1)$-dimensions using two examples : i) a Tanh-like expansion and ii) a de-Sitter expansion. We then extend the analysis to leading order fluctuations in $(3+1)-$dimensions to show that a quantum-to-classical transition occurs in the de-Sitter expansion and discuss the relevance of our analysis in distinguishing cosmological models.