Paris Observatory

Abstract Gammapy 1is a community-driven 2open-source Python package to analyze very high-energy gamma-ray astronomical data. Created in 2014, it has since expanded to support analysis methods in multiwavelength and multimessenger astrophysics. Currently, in version 1.2, Gammapy is utilized by various instruments like HESS, VERITAS, HAWC, and Fermi-LAT, and tested with X-ray and neutrino data, and is the base for the science analysis tools of the Cherenkov Telescope Array Observatory 3. Operating on the common open FITS-based data format GADF, it ensures interoperability between different experiments. The package’s analysis pipeline handles data, instrument response functions, and dependencies on time, energy, and sky location. It offers diverse methods for background estimation and data reduction, producing standardized datasets. Gammapy includes spatial, spectral, and temporal models, supporting custom model integration and permitting maximum-likelihood fitting on datasets. Additionally, it provides tools for accessing flux points, likelihood profiles, and light curves, facilitating comprehensive time-domain analyses. Gammapy has de facto become the reference software for VHE gamma-ray analysis, strives to be integrated into the open science framework with the FAIR4RS guidelines [20], and will adopt the VODF format in future versions.

Understanding the transport of energy within space plasmas, particularly in the solar wind, remains a complex challenge. Accurate measurement of electron temperatures and their nonthermal characteristics is crucial for comprehending energy transport properties in plasmas. Quasi-thermal-noise (QTN) spectroscopy has emerged as a dependable tool for precise electron parameters assessment as it is less susceptible to spacecraft effects than particle detectors. In this study, we apply a QTN spectroscopy fitting method to analyze data from the Parker Solar Probe FIELDS radio instrument obtained during Encounters 2 through 13, under unbiased antenna conditions. We use the kappa function to characterize the electron velocity distribution and employ a fitting technique to derive the changes in each parameter across heliocentric distances ranging from 12 Rs to 76 Rs. Specifically, we find that the electron density scales as n e ∝ r −2.09±0.04 and the T e ∝ r −0.65±0.02. The distribution of the kappa index has three distinct regions as a function of radial distance from the Sun. Furthermore, we conduct a statistical analysis of solar wind energy flux which we finds follows a power-law relationship w total ∝ r −1.92±0.04.

We present the results of a search for continuous gravitational wave signals (CGWs) in the second data release (DR2) of the European Pulsar Timing Array (EPTA) Collaboration. The most significant candidate event from this search has a gravitational wave frequency of 4–5 nHz. Such a signal could be generated by a supermassive black hole binary (SMBHB) in the local Universe. We present the results of a follow-up analysis of this candidate using both Bayesian and frequentist methods. The Bayesian analysis gives a Bayes factor of 4 in favour of the presence of the CGW over a common uncorrelated noise process. In contrast, the frequentist analysis estimates the p-value of the candidate to be < 1%, also assuming the presence of common uncorrelated red noise. However, comparing a model that includes both a CGW and a gravitational wave background (GWB) to a GWB only, the Bayes factor in favour of the CGW model is only 0.7. Therefore, we cannot conclusively determine the origin of the observed feature, nor can we rule it out as a CGW source. We present results of simulations that demonstrate that data containing a weak gravitational wave background can be misinterpreted as data including a CGW and vice versa, providing two plausible explanations for the EPTA DR2 data. Further investigations combining data from all PTA collaborations will be needed to reveal the true origin of this feature.

Abstract Measurements of ground and near surface atmospheric temperatures at Jezero obtained during 700 sols by the Mars Environmental Dynamics Analyzer (MEDA) characterize the thermal behavior of the near surface Martian atmosphere during a full Martian Year. The seasonal evolution of MEDA measurements is compared with predictions from the Mars Climate Database and the solar irradiance at the surface. Thermal tides observed in the daily cycle of temperatures follow a seasonal cycle with additional variations greater than 2 K on time‐scales of tens of sols. We also observe sol‐to‐sol variations of about 1 K in mean daily air temperatures in autumn and winter with periodicities of 4–7 sols that might be related to baroclinic disturbances that are frequent in those seasons at high latitudes. We examine the evolution of the vertical thermal gradient and temperature fluctuations without finding a seasonal response to irradiance and dust load. We find that the convective boundary layer becomes isothermal and collapses 1 hr before sunset except during northern hemisphere winter, when the collapse occurs closer to sunset, implying a longer duration of the daytime convective instability. Around this period, the rover was located in the delta front in a location of complex topography where we observed stronger thermal gradients and intense daytime air temperature fluctuations. We also find in this place a nighttime event of gravity waves on near‐surface air temperatures, with amplitudes of 2 K and periods of 10 min. These waves possibly propagate downward through a near isothermal stable layer.

Abstract Mercury's planetary magnetic field models (PMFMs) agree on a majorly dipolar field structure with a northward shift of the magnetic equator. However, due to the northerly biased orbit coverage of past spacecraft missions and different data analyzing methods, the available PMFMs differ in the determined multipole magnitudes for the dipole, quadrupole and octupole moments. While the PMFMs agree well with northern observations, we find that the predicted magnetic field values differ in the unexplored equatorial and southern regions. In the forward modeling approach of this study, we apply three different PMFM representatives, differing in their values of the dipole, quadrupole and octupole moments, and model the resulting solar wind interaction via a global hybrid model with sets of the 4 most common interplanetary magnetic field (IMF) directions under otherwise average solar wind conditions. Extracting our modeled fields along the flybys of MESSENGER and BepiColombo, as well as along four representative orbits of the Mercury Planetary Orbiter (MPO), allows us to estimate local field variations due to IMF and PMFM influence. Our modeled magnetic field components separate significantly by up to 60 nT between the PMFMs in low‐altitude southern regions, leading to displacements of magnetopause, cusps and current sheet locations. We find that MESSENGER flyby observations agree best with modeled field ranges resulting from a quadrupolar PMFM and that certain nightside regions are useful to estimate upstream IMF polarity. We also demonstrate how comparing BepiColombo swingby and MPO orbit phase observations with modeled results can help determine the correct PMFM for Mercury.

Abstract We propose high-Q factors sub-stoichiometric TiN-based Lumped Element MKIDs on sapphire with vacuum gap suspended meander. The aim is to improve detector response by trapping phonons within the meander to increase the number of Cooper pairs breaking into quasiparticles. Furthermore, either Nb/Au or Nb/Al reflective thin bilayer can be placed beneath \textcolor{black}{the meander} to respectively improve the photon peak absorption at the near-infrared and visible wavelengths.

We modeled the geometry and the three-dimensional orientation of the ejecta cone triggered by the impact of the DART spacecraft on the asteroid Dimorphos. We used eight LUKE images of the impact acquired by the CubeSat LICIACube that flew by the Didymos system shortly after the impact. These images, which show the ejecta cone in both face-on and side-on profiles, enabled us to reconstruct the ejecta cone in inertial space. We started our model as a simple cone with a circular base and developed it to a rotated cone with an elliptical base that best fit the data. The cone axis points to R.A., decl. (in J2000): , +. The cone is characterized by two perpendicular half-angles of and a rotation of ω = 12° around its axis. The apex of the cone is located near the center of Dimorphos within 15 m. The intersection of the cone and the surface of Dimorphos (surface enclosed by the cone) would correspond to a crater with a maximum radius of about 65 m. The characterization of the cone axis is directly related to the computation of the momentum enhancement factor (β) of the impact, and it hence proves the crucial need of studying impacts in the context of planetary defence scenarios. The results of this work could potentially be used to constrain whether the impact took place in a strength-dominated or a gravity-dominated regime. This work shows the important scientific return of the LICIACube CubeSat in the context of planetary defence.

High-resolution near-infrared ground-based spectroscopic observations of comet 67P/Churyumov–Gerasimenko near its maximum activity in 2021 were conducted from the W. M. Keck Observatory, using the facility spectrograph NIRSPEC. 67P is the best-studied comet to date because of the unprecedented detail and insights provided by the Rosetta mission during 2014–2016. Because 67P is the only comet where the detailed abundances of many coma volatiles were measured in situ, determining its composition from the ground provides a unique opportunity to interpret Rosetta results within the context of the large database of ground-based compositional measurements of comets. However, previous apparitions, including in 2015, have been unfavorable for in-depth ground-based studies of parent volatiles in 67P. The 2021 apparition of 67P was thus the first-ever opportunity for such observations. We report gas spatial distributions, rotational temperatures, production rates, and relative abundances (or stringent upper limits) among seven volatile species: C2H2, C2H6, HCN, NH3, CH3OH, H2CO, and H2O. The measured abundances of trace species relative to water reveal near average or below average values compared to previous comets studied at infrared wavelengths. Both gas rotational temperatures and the spatial distributions of H2O, C2H6, and HCN measured with Keck-NIRSPEC in 2021 are consistent with the outgassing patterns revealed by Rosetta in 2015 at very similar heliocentric distance (post-perihelion). These results can be integrated with both Rosetta mission findings and ground-based cometary studies of the overall comet population, for which we encourage a wide-scale collaboration across measurement techniques.

The Event Horizon Telescope (EHT) is a millimeter very long baseline interferometry (VLBI) array that has imaged the apparent shadows of the supermassive black holes M87* and Sagittarius A*. Polarimetric data from these observations contain a wealth of information on the black hole and accretion flow properties. In this work, we develop polarimetric geometric modeling methods for mm-VLBI data, focusing on approaches that fit data products with differing degrees of invariance to broad classes of calibration errors. We establish a fitting procedure using a polarimetric “m-ring” model to approximate the image structure near a black hole. By fitting this model to synthetic EHT data from general relativistic magnetohydrodynamic models, we show that the linear and circular polarization structure can be successfully approximated with relatively few model parameters. We then fit this model to EHT observations of M87* taken in 2017. In total intensity and linear polarization, the m-ring fits are consistent with previous results from imaging methods. In circular polarization, the m-ring fits indicate the presence of event-horizon-scale circular polarization structure, with a persistent dipolar asymmetry and orientation across several days. The same structure was recovered independently of observing band, used data products, and model assumptions. Despite this broad agreement, imaging methods do not produce similarly consistent results. Our circular polarization results, which imposed additional assumptions on the source structure, should thus be interpreted with some caution. Polarimetric geometric modeling provides a useful and powerful method to constrain the properties of horizon-scale polarized emission, particularly for sparse arrays like the EHT.