Papers for Monday, May 01 2023

We present our study on the reconstruction of photoelectron tracks in gas pixel detectors used for astrophysical X-ray polarimetry. Our work aims to maximize the performance of convolutional neural networks (CNNs) to predict the impact point of incoming X-rays from the image of the photoelectron track. A very high precision in the reconstruction of the impact point position is achieved thanks to the introduction of an artificial sharpening process of the images. We find that providing the CNN-predicted impact point as input to the state-of-the-art analytic analysis improves the modulation factor ($\sim 1 \%$ at 3 keV and $\sim 6 \%$ at 6 keV) and naturally mitigates a subtle effect appearing in polarization measurements of bright extended sources known as "polarization leakage".

We have studied gas-grain chemical models of interstellar clouds to search for nonlinear dynamical evolution. A prescription is given for producing oscillatory solutions when a bistable solution exists in the gas-phase chemistry and we demonstrate the existence of limit cycle and relaxation oscillation solutions. As the autocatalytic chemical processes underlying these solutions are common to all models of interstellar chemistry, the occurrence of these solutions should be widespread. We briefly discuss the implications for interpreting molecular cloud composition with time-dependent models and some future directions for this approach.

Tidal disruption events (TDEs), in which a star is destroyed by the gravitational field of a supermassive black hole (SMBH), are being observed at a high rate owing to the advanced state of survey science. One of the properties of TDEs that is measured with increasing statistical reliability is the TDE luminosity function, $d\dot{N}_{\rm TDE}/dL$, which is the TDE rate per luminosity (i.e., how many TDEs are within a given luminosity range). Here we show that if the luminous emission from a TDE is directly coupled to the rate of return of tidally destroyed debris to the SMBH, then the TDE luminosity function is in good agreement with observations and scales as $\propto L^{-2.5}$ for high luminosities, provided that the SMBH mass function $dN_{\bullet}/dM_{\bullet}$ -- the number of SMBHs ($N_{\bullet}$) per SMBH mass ($M_{\bullet}$) -- is approximately flat in the mass range over which we observe TDEs. We also show that there is a cutoff in the luminosity function at low luminosities that is a result of direct captures, and this cutoff has been tentatively observed. If $dN_{\bullet}/dM_{\bullet}$ is flat, which is in agreement with some observational campaigns, these results suggest that the fallback rate feeds the accretion rate in TDEs. Contrarily, if $dN_{\bullet}/d\log M_{\bullet}$ is flat, which has been found theoretically and is suggested by other observational investigations, then the emission from TDEs is likely powered by another mechanism. Future observations and more TDE statistics, provided by the Rubin Observatory/LSST, will provide additional evidence as to the reality of this tension.

We present a new suite of over 1,500 cosmological N-body simulations with varied Warm Dark Matter (WDM) models ranging from 2.5 to 30 keV. We use these simulations to train Convolutional Neural Networks (CNNs) to infer WDM particle masses from images of DM field data. Our fiducial setup can make accurate predictions of the WDM particle mass up to 7.5 keV at a 95% confidence level from small maps that cover an area of (25 h$^{-1}$ Mpc)$^2$. We vary the image resolution, simulation resolution, redshift, and cosmology of our fiducial setup to better understand how our model is making predictions. Using these variations, we find that our models are most dependent on simulation resolution, minimally dependent on image resolution, not systematically dependent on redshift, and robust to varied cosmologies. We also find that an important feature to distinguish between WDM models is present with a linear size between 100 and 200 h$^{-1}$ kpc. We compare our fiducial model to one trained on the power spectrum alone and find that our field-level model can make 2x more precise predictions and can make accurate predictions to 2x as massive WDM particle masses when used on the same data. Overall, we find that the field-level data can be used to accurately differentiate between WDM models and contain more information than is captured by the power spectrum. This technique can be extended to more complex DM models and opens up new opportunities to explore alternative DM models in a cosmological environment.

The mean matter density within the turnaround radius, which is the boundary that separates a nonexpanding structure from the Hubble flow, was recently proposed as a novel cosmological probe. According to the spherical collapse model, the evolution with cosmic time of this turnaround density, $\rm \rho_{ta}(z)$, can be used to determine both $\rm \Omega_m$ and $\Omega_\Lambda$, independently of any other currently used probe. The properties of $\rm \rho_{ta}$ predicted by the spherical collapse model were also shown to persist in the presence of full three-dimensional effects in $\rm \Lambda$CDM N-body cosmological simulations when considering galaxy clusters at the present time, $z=0$. However, a small offset was discovered between the spherical-collapse prediction of the value of $\rho_{ta}$ at $z=0$ and its value measured in simulations. In this letter, we explore whether this offset evolves with cosmic time; whether it differs in different cosmologies; whether its origin can be confidently identified; and whether it can be corrected. We found that the offset does evolve slightly with redshift, and that it correlates strongly with the deviation from spherical symmetry of the dark matter halo distribution inside and outside of the turnaround radius. We used an appropriate metric to quantify deviations in the environment of a structure from spherical symmetry. We found that using this metric, we can construct a sphericity-selected sample of halos for which the offset of $\rho_{ta}$ from the spherical collapse prediction is zero, independently of redshift and cosmology. We found that a sphericity-selected halo sample allows us to recover the simulated cosmology, and we conclude that the turnaround density evolution indeed encodes the cosmology in N-body simulations.

The demographic of circumstellar disks, the birthplaces of planets, is diverse and rich in disks featuring rings, gaps, spirals, filaments, and arcs. Many studies revealing these disk structures have focused on objects around single stars and disks in isolation. The scenario is more complex if binarity or multiplicity is involved; most stars are part of multiple systems in crowded star-forming regions. How does the presence of one or more stellar companions affect the shape and size of the circumstellar disks? Here we review the landscape of results from optical, infrared, and (sub-) millimeter observations of the effects of multiplicity on protoplanetary disks, emphasizing the demographic studies of nearby molecular clouds and the high-resolution studies of multiple disk systems.

Since identifying the gap in the H-R Diagram (HRD) marking the transition between partially and fully convective interiors, a unique type of slowly pulsating M dwarf has been proposed. These unstable M dwarfs provide new laboratories in which to understand how changing interior structures result in potentially observable activity at the surface. In this work, we report the results of the largest high-resolution spectroscopic H$\alpha$ emission survey to date spanning this transition region, including 480 M dwarfs observed using the CHIRON spectrograph at CTIO/SMARTS 1.5-m. We find that M dwarfs with H$\alpha$ in emission are almost entirely found 0 to 0.5 magnitude above the top edge of the gap in the HRD, whereas effectively no stars in and below the gap show emission. Thus, the top edge of the gap marks a relatively sharp activity transition and there is no anomalous H$\alpha$ activity for stars in the gap. We also identify a new region at 10.3 $<M_{G}<$ 10.8 on the main sequence where fewer M dwarfs exhibit H$\alpha$ emission compared to M dwarfs above and below this magnitude range. Careful evaluation of literature results indicates that 1) rotation and H$\alpha$ activity distributions on the main sequence are closely related, and 2) fewer stars in this absolute magnitude range rotate in less than $\sim$13 days than populations surrounding this region. This result suggests that the most massive fully convective stars lose their angular momentum faster than both partially convective stars and less massive fully convective stars.

We present a new determination of the evolving galaxy UV luminosity function (LF) over the redshift range $9.5<z<12.5$ based on a wide-area ($>250$ arcmin$^2$) data set of JWST NIRCam near-infrared imaging assembled from thirteen public JWST surveys. Our relatively large-area search allows us to uncover a sample of 61 robust $z>9.5$ candidates detected at $\geq 8\sigma$, and hence place new constraints on the intermediate-to-bright end of the UV LF. When combined with our previous JWST+UltraVISTA results, this allows us to measure the form of the LF over a luminosity range corresponding to four magnitudes ($M_{1500}$). At these early times we find that the galaxy UV LF is best described by a double power-law function, consistent with results obtained from recent ground-based and early JWST studies at similar redshifts. Our measurements provide further evidence for a relative lack of evolution at the bright-end of the UV LF at $z=9-11$, but do favour a steep faint-end slope ($\alpha\leq-2$). The luminosity-weighted integral of our evolving UV LF provides further evidence for a gradual, smooth (exponential) decline in co-moving star-formation rate density ($\rho_{\mathrm{SFR}}$) at least out to $z\simeq12$, with our determination of $\rho_{\mathrm{SFR}}(z=11)$ lying significantly above the predictions of many theoretical models of galaxy evolution.

We investigate the effects of an oscillating gas filament on the dynamics of its embedded stellar clusters. Motivated by recent observational constraints, we model the host gas filament as a cylindrically symmetrical potential, and the star cluster as a Plummer sphere. In the model, the motion of the filament will produce star ejections from the cluster, leaving star cluster remnants that can be classified into four categories: a) Filament Associated clusters, which retain most of their particles (stars) inside the cluster and inside the filament; b) destroyed clusters, where almost no stars are left inside the filament, and there is no surviving bound cluster; c) ejected clusters, that leave almost no particles in the filament, since the cluster leaves the gas filament; and d) transition clusters, corresponding to those clusters that remain in the filament, but that lose a significant fraction of particles due to ejections induced by filament oscillation. Our numerical investigation predicts that the Orion Nebula Cluster is in the process of being ejected, after which it will most likely disperse into the field. This scenario is consistent with observations which indicate that the Orion Nebula Cluster is expanding, and somewhat displaced from the Integral Shaped Filament ridgeline.

We have investigated the chemistry of dense interstellar clouds and found new bistable solutions in the nitrogen and carbon chemistries. We identify the autocatalytic processes that are present in the pure, reduced, chemical networks and, as previously found for oxygen chemistry, that He$^+$ plays an important role. The applicability of these results to astronomical environments is briefly discussed. The bistable solutions found for carbon chemistry occur for low densities and high ionization fractions that are not compatible with that found cold, dense clouds. Bistability in the pure nitrogen chemistry occurs for conditions that are relevant for prestellar cores in which significant CO depletion has taken place. We conclude that several autocatalyses are embedded in gas-phase interstellar chemistry and that many more are potentially present.

Small solar system bodies have widely dispersed orbital poles, posing challenges to dynamical models of solar system origin and evolution. To characterize the orbit pole distribution of dynamical groups of small bodies it helps to have a functional form for a model of the distribution function. Previous studies have used the small-inclination approximation and adopted variations of the normal distribution to model orbital inclination dispersions. Because the orbital pole is a directional variable, its distribution can be more appropriately modeled with directional statistics. We describe the von Mises-Fisher (vMF) distribution on the surface of the unit sphere for application to small bodies' orbital poles. We apply it to the orbit pole distribution of the observed Plutinos. We find a mean pole located at inclination of 3.57 degrees and a longitude of ascending node of 124.38 degrees (in the J2000 reference frame), with a 99.7 per cent confidence cone of half-angle 1.68 degrees. We also estimate a debiased mean pole located 4.6 degrees away, at an inclination of 2.26 degrees and a longitude of ascending node of 292.69 degrees, of similar-size confidence cone. The vMF concentration parameter of Plutino inclinations (relative to either mean pole estimate) is 31.6. This resembles a Rayleigh distribution function with a width parameter of 10.2 degrees. Unlike previous models, the vMF model naturally accommodates all physical inclinations (and no others), whereas Rayleigh or Gaussian models must be truncated to the physical inclination range 0-180 degrees. Further work is needed to produce a theory for the mean pole of the Plutinos against which to compare the observational results.

The Ultraviolet Transient Astronomy Satellite (ULTRASAT) is scheduled to be launched to geostationary orbit in 2026. It will carry a telescope with an unprecedentedly large field of view (204 deg$^2$) and NUV (230-290nm) sensitivity (22.5 mag, 5$\sigma$, at 900s). ULTRASAT will conduct the first wide-field survey of transient and variable NUV sources and will revolutionize our ability to study the hot transient universe: It will explore a new parameter space in energy and time-scale (months long light-curves with minutes cadence), with an extra-Galactic volume accessible for the discovery of transient sources that is $>$300 times larger than that of GALEX and comparable to that of LSST. ULTRASAT data will be transmitted to the ground in real-time, and transient alerts will be distributed to the community in $<$15 min, enabling a vigorous ground-based follow-up of ULTRASAT sources. ULTRASAT will also provide an all-sky NUV image to $>$23.5 AB mag, over 10 times deeper than the GALEX map. Two key science goals of ULTRASAT are the study of mergers of binaries involving neutron stars, and supernovae: With a large fraction ($>$50%) of the sky instantaneously accessible, fast (minutes) slewing capability and a field-of-view that covers the error ellipses expected from GW detectors beyond 2025, ULTRASAT will rapidly detect the electromagnetic emission following BNS/NS-BH mergers identified by GW detectors, and will provide continuous NUV light-curves of the events; ULTRASAT will provide early (hour) detection and continuous high (minutes) cadence NUV light curves for hundreds of core-collapse supernovae, including for rarer supernova progenitor types.

The BL Lac Markarian 501 exhibited two flaring activities in the very-high-energy (VHE) band in May 2009. The lack of correlation between X-rays and TeV gamma-rays without increasing in other bands suggested that more than one emission zone could be involved. Moreover, fast variability in the flaring state was observed, indicating that the emission zones responsible must have small sizes. We use a lepto-hadronic model with two-zone emission to explain the spectral energy distribution during quiescent and these flaring states. In the proposed scenario, the photopion processes explain the VHE flaring activities successfully, and variability constraints place the activity in a zone located near the jet's base or named inner blob, while synchrotron self-Compton emission describing the X-ray signature during that flaring state occurs in the zone situated far the central engine or named outer blob.

We study interacting galaxy pairs in the TNG100-1 and TNG300-1 cosmological simulations using previously generated closest companion samples. We study the specific star formation rates (sSFR) of massive ($10^{10} M_{\odot} < M_* < 10^{12} M_{\odot}$) galaxies at $z \leq 0.2$ as a function of separation from the closest companion galaxy. We split our sample based on whether the companion galaxy is star-forming or passive. We find that galaxies with close star-forming companions have sSFRs that are enhanced (on average) by a factor of $2.9 \pm 0.3$ in TNG100-1 and $2.27 \pm 0.06$ in TNG300-1 compared to controls, with enhancements present out to separations of $\sim 300$ kpc. Galaxies with passive companions in TNG300-1 exhibit mild sSFR suppression ($\sim12$ percent) at 100-300 kpc and small sSFR enhancements at separations below 50 kpc. sSFR suppression is strongest in pairs where the galaxy's stellar mass is more than 2 times that of its passive companion. By generating a stellar mass-matched ("twinned") sample in TNG300-1, we show that differences in sSFR trends between companion types are not a result of intrinsic stellar mass differences in star-forming vs. passive galaxies. We compare with an analogous sample of galaxy pairs from SDSS, finding consistent results between observations and simulations. Overall, we find that star-forming galaxies show enhanced sSFRs regardless of companion type, but that galaxies with close passive companions are more likely to be passive themselves.

In September 2021, a site scouting mission known as the TRIDENT pathfinder experiment (TRIDENT EXplorer, T-REX for short) was conducted in the South China Sea with the goal of envisaging a next-generation multi-cubic-kilometer neutrino telescope. One of the main tasks is to measure the \textit{in-situ} optical properties of seawater at depths between $2800~\mathrm{m}$ and $3500~\mathrm{m}$, where the neutrino telescope will be instrumented. To achieve this, we have developed a light emitter module equipped with a clock synchronization system to serve as the light source, which could be operated in pulsing and steady modes. Two light receiver modules housing both photomultiplier tubes (PMTs) and cameras are employed to detect the photons emitted by the light source. This paper presents the instrumentation of the light source in T-REX, including its design, calibration, and performance.

The Kepler and TESS space missions significantly expanded our knowledge of what types of stars display flaring activity by recording a vast amount of super-flares from solar-like stars, as well as detecting flares from hotter stars of A-F spectral types. Currently, we know that flaring occurs in the stars as hot as B-type ones. However, the structures of atmospheres of hot B-A stars crucially differ from the ones of late types, and thus the occurrence of flaring in B-A type stars requires some extension of our theoretical views of flare formation and therefore a detailed study of individual objects. Here we present the results of our spectral and photometric study of HD 36030, which is a B9 V star with flares detected by the TESS satellite. The spectra we acquired suggest that the star is in a binary system with a low-mass secondary component, but the light curve lacks any signs of periodic variability related to orbital motion or surface magnetic fields. Because of that, we argue that the flares originate due to magnetic interaction between the components of the system.

We report the observations of FRB 20220912A using the Five-hundred-meter Aperture Spherical radio Telescope (FAST). We conducted 17 observations totaling 8.67 hours and detected a total of 1076 bursts with an event rate up to 390 hr$^{-1}$. The cumulative energy distribution can be well described using a broken power-law function with the lower and higher-energy slopes of $-0.38\pm0.02$ and $-2.07\pm0.07$, respectively. We also report the L band ($1-1.5$ GHz) spectral index of the synthetic spectrum of FRB~20220912A bursts, which is $-2.6\pm0.21$. The average rotation measure (RM) value of the bursts from FRB~20220912A is $-0.08\pm5.39\ \rm rad\,m^{-2}$, close to 0 $\rm rad\,m^{-2}$ and maintain relatively stable over two months. Most bursts have nearly 100\% linear polarization. About 45\% of the bursts have circular polarization with SNR $>$ 3, and the highest circular polarization degree can reach 70\%. Our observations suggest that FRB~20220912A is located in a relatively clean local environment with complex circular polarization characteristics. These various behaviors imply that the mechanism of circular polarization of FRBs likely originates from an intrinsic radiation mechanism, such as coherent curvature radiation or inverse Compton scattering inside the magnetosphere of the FRB engine source (e.g. a magnetar).

IRAS 00500+6713 is a hypothesized remnant of a type Iax supernova SN 1181. Multi-wavelength observations have revealed its complicated morphology; a dusty infrared ring is sandwiched by the inner and outer X-ray nebulae. We analyze the archival X-ray data taken by XMM-Newton and Chandra to constrain the angular size, mass, and metal abundance of the X-ray nebulae, and construct a theoretical model describing the dynamical evolution of IRAS 00500+6713, including the effects of the interaction between the SN ejecta and the intense wind enriched with carbon burning ashes from the central white dwarf (WD) J005311. We show that the inner X-ray nebula corresponds to the wind termination shock while the outer X-ray nebula to the shocked interface between the SN ejecta and the interstellar matter. The observed X-ray properties can be explained by our model with an SN explosion energy of $E_\mathrm{ej} = (0.77 \mbox{--} 1.1)\times 10^{48}$~erg, an SN ejecta mass of $M_\mathrm{ej} = 0.18\mbox{--}0.53~M_\odot$, if the currently observed wind from WD J005311 started to blow $t_\mathrm{w} \gtrsim 810$ yr after the explosion, i.e., approximately after A.D. 1990. The inferred SN properties are compatible with those of Type Iax SNe and the timing of the wind launch may correspond to the Kelvin-Helmholtz contraction of the oxygen-neon core of WD J005311 that triggered a surface carbon burning. Our analysis supports that IRAS 00500+6713 is the remnant of SN Iax 1181 produced by a double degenerate merger of oxygen-neon and carbon-oxygen WDs, and WD J005311 is the surviving merger product.

Fast radio bursts (FRBs) are bright millisecond radio bursts at cosmological distances. Only three FRBs have exhibited extreme activities, such as achieving a peak event rate $\gtrsim 100$ hr$^{-1}$ or being persistently active. Only these three among $\sim 50$ known repeating FRBs have circular polarization. We observed the FRB 20220912A with the Robert C. Byrd Green Bank Telescope (GBT) at L-band on 24 October 2022 and detected 128 bursts in 1.4 hours, corresponding to a burst rate of about 90 hr$^{-1}$, which is the highest yet for FRBs observed by the GBT and makes it the fourth extremely active FRB. The median energy of the bursts is $4.0\times10^{37}$ erg, close to the characteristic energy of FRB 20121102A. The average rotation measure (RM) was $-$0.4 rad m$^{-2}$ with unnoticeable intraday RM change, indicating a likely clean environment, in contrast to the other three extremely active repeating FRBs. Most bursts have nearly 100% linear polarization. Approximately 56% of the bright bursts have circular polarization, the highest such fraction among all FRBs. A downward drift in frequency and polarization angle swings were found in our sample. The discovery and characterization of FRB 20220912A support the view that the downward drift in frequency, polarization angle swings, and circular polarization are intrinsic to radiation physics, which may be shared by active repeaters regardless of the environments.

This study aims to investigate the systematics in planetary ephemeris reference frames through pulsar timing observations. We used the published data sets from several pulsar timing arrays and performed timing analyses for each pulsar using different planetary ephemerides retrieved from the Jet Propulsion Laboratory's Development Ephemeris (DE), Ephemeris of Planets and the Moon (EPM), and INPOP (Int\'egration Num\'erique Plan\'etaire de l'Observatoire de Paris). Then, we compared the timing solutions and modeled the differences in position and proper motion by vector spherical harmonics of the first degree. The timing solutions were also compared with those determined by very long baseline interferometry (VLBI) astrometry. The orientation offsets between the latest editions of the DE, EPM, and INPOP series do not exceed 0.4 milliarcseconds (mas), while the relative spins between these ephemerides are less than 5 microarcseconds per year ($\mathrm{\mu as\,yr^{-1}}$). We do not detect significant glides in either position or proper motion between these ephemerides. The orientation of the pulsar timing frames deviates from that of the VLBI frame from zero by approximately $\mathrm{0.4\,mas}$ when considering the formal uncertainty and possible systematics. The orientation of current planetary ephemeris frames is as accurate as at least 0.4 mas, and the nonrotating is better than $\mathrm{5\,\mu as\,yr^{-1}}$.

We present MeerKAT L-band (886-1682 MHz) observations of the extended radio structure of the peculiar galaxy pair PKS 2130-538 known as the "Dancing Ghosts". The complex of bending and possibly interacting jets and lobes originate from two Active Galactic Nuclei hosts in the Abell 3785 galaxy cluster, one of which is the brightest cluster galaxy. The radio properties of the PKS 2130-538 flux density, spectral index and polarization - are typical for large, bent-tail galaxies. We also investigate a number of thin extended low surface brightness filaments originating from the lobes. Southeast from the Dancing Ghosts, we detect a region of low surface brightness emission that has no clear origin. While it could originate from the Abell 3785 radio halo, we investigate the possibility that it is a radio relic or emission associated with the two PKS 2130-538 hosts. We find no evidence of interaction between the two PKS 2130-538 hosts.

In this second communication we continue our analysis of the turbulence in the Huygens Region of the Orion Nebula (M 42). We calculate the associated transverse structure functions up to order 8-th and find that the higher-order transverse structure functions are almost proportional to the second-order transverse structure function: we find that after proper normalisation, the higher-order transverse structure functions only differ by very small deviations from the second-order transverse structure function in a sub-interval of the inertial range. We demonstrate that this implies that the turbulence in the Huygens Region is quasi-log-homogeneous, or to a better degree of approximation, binomially weighted log-homogeneous in the statistical sense, this implies that there is some type of invariant statistical structure in the velocity field of the Huygens Region. We also obtain and analyse the power-spectrum of the turbulent field and find that it displays a large tail that follows very approximately two power-laws, one of the form $E(k)\propto k^{-2.7}$ for the initial side of the tail, and one of the form $E(k)\propto k^{-1}$ for the end of the tail. We find that the power-law with exponent $\beta\sim -2.7$ corresponds to spatial scales of 0.0301--0.6450 pc. We find that the exponent of the first power-law $\beta\sim -2.7$ is related to the exponent $\alpha_2$ of the second-order structure function in the inertial range. We interpret the second power-law with exponent $\beta \sim -1$ as an indicator of viscous-dissipative processes occurring at scales of $\delta r=1$--5 pixels which correspond to spatial scales of 0.00043--0.00215 pc.

We present the interstellar scintillation analysis of fast radio burst (FRB) 20220912A during its extremely active episode in 2022 using data from the Five-hundred-meter Aperture Spherical Radio Telescope (FAST). We detect a scintillation arc in the FRB's secondary spectrum, which describes the power in terms of the scattered FRB signals' time delay and Doppler shift. The arc indicates that the scintillation is caused by a highly localized region of the ionized interstellar medium (IISM). Our analysis favors a Milky Way origin for the localized scattering medium but cannot rule out a host galaxy origin. We present our method for detecting the scintillation arc, which can be applied generally to sources with irregularly spaced bursts or pulses. These methods could help shed light on the complex interstellar environment surrounding the FRBs and in our Galaxy.

We present the results of a direct measurement of the cosmic-ray helium spectrum with the CALET instrument in operation on the International Space Station since 2015. The observation period covered by this analysis spans from October 13, 2015 to April 30, 2022 (2392 days). The very wide dynamic range of CALET allowed to collect helium data over a large energy interval, from ~40 GeV to ~250 TeV, for the first time with a single instrument in Low Earth Orbit. The measured spectrum shows evidence of a deviation of the flux from a single power-law by more than 8$\sigma$ with a progressive spectral hardening from a few hundred GeV to a few tens of TeV. This result is consistent with the data reported by space instruments including PAMELA, AMS-02, DAMPE and balloon instruments including CREAM. At higher energy we report the onset of a softening of the helium spectrum around 30 TeV (total kinetic energy). Though affected by large uncertainties in the highest energy bins, the observation of a flux reduction turns out to be consistent with the most recent results of DAMPE. A Double Broken Power Law (DBPL) is found to fit simultaneously both spectral features: the hardening (at lower energy) and the softening (at higher energy). A measurement of the proton to helium flux ratio in the energy range from 60 GeV/n to about 60 TeV/n is also presented, using the CALET proton flux recently updated with higher statistics.

We report the mass determination of TOI-519 b, a transiting substellar object around a mid-M dwarf. We carried out radial velocity measurements using Subaru / InfraRed Doppler (IRD), revealing that TOI-519 b is a planet with a mass of $0.463^{+0.082}_{-0.088}~M_{\rm Jup}$. We also find that the host star is metal rich ($\rm [Fe/H] = 0.27 \pm 0.09$ dex) and has the lowest effective temperature ($T_{\rm eff}=3322 \pm 49$ K) among all stars hosting known close-in giant planets based on the IRD spectra and mid-resolution infrared spectra obtained with NASA Infrared Telescope Facility / SpeX. The core mass of TOI-519 b inferred from a thermal evolution model ranges from $0$ to $\sim30~M_\oplus$, which can be explained by both the core accretion and disk instability models as the formation origins of this planet. However, TOI-519 is in line with the emerging trend that M dwarfs with close-in giant planets tend to have high metallicity, which may indicate that they formed in the core accretion model. The system is also consistent with the potential trend that close-in giant planets around M dwarfs tend to be less massive than those around FGK dwarfs.

Solar activity and solar wind parameters decreased significantly in solar cycles (SCs) 23-24. In this paper, we analyze solar wind measurements at the rising phase of SC 25 and compare them with similar data from the previous cycles. For this purpose, we simultaneously selected the OMNI database data for 1976-2022, both by phases of the 11-year solar cycle and by large-scale solar wind types (in accordance with IKI's catalog, see this http URL ), and calculated the mean values of the parameters for the selected datasets. The obtained results testify in favor of the hypothesis that the continuation of this cycle will be similar to the previous cycle 24, i.e. SC 25 will be weaker than SCs 21 and 22.

ASCC 123 is a little-studied young and dispersed open cluster. Recently, we conducted the first research devoted to it. In this paper, we complement our previous work with TESS photometry for the 55 likely members of the cluster. We pay special attention to seven of these high-probability members, all with FGK spectral types, for which we have high-resolution spectra from our preceding work. By studying the TESS light curves of the cluster members we determine the rotational period and the amplitude of the rotational modulation for 29 objects. The analysis of the distribution of the periods allows us to estimate a gyrochronogical age for ASCC 123 similar to that of the Pleiades, confirming the value obtained in our previous investigation. A young cluster age is also suggested by the distribution of variation amplitudes. In addition, for those stars with spectroscopic data we calculate the inclination of their rotation axis. These values appear to follow a random distribution, as already observed in young clusters, with no indication of spin alignment. However, our sample is too small to confirm this on more solid statistical grounds. Finally, for these seven stars we study the level of magnetic activity from the H$\alpha$ and CaII H&K lines. Despite the small number of data points, we find a correlation of the H$\alpha$ and CaII flux with Rossby number. The position of these stars in flux--flux diagrams follows the general trends observed in other active late-type stars.

Age-metallicity relation for the Galactic disc is a crucial tool and to constrain the Galactic chemical evolution models. We investigate the age-metallicity relation of the Galactic disc using the red giant branch stars in the Solar neighbourhood. The data cover the Galactocentric radius of $7\leq R_{\rm gc} (\rm kpc) \leq9.5$, but extends up to 4 kpc in height from the Galactic plane. We use kinematic age derived from highly precise astrometric data of Gaia Data Release 2 and element abundance ratios from high-resolution spectroscopic data of APOGEE-2 catalogues. We apply a two-component Gaussian mixture model to chemically separate the programme stars into thin and thick disc populations. The stars in each population are grouped into different distance intervals from the Galactic plane. The mean metal abundances and velocity dispersions of the stars in the groups were calculated and the kinematic ages were determined from their kinematic parameters. We found a steep relation for the thin disc with -0.057$\pm$0.007 dex Gyr$^{-1}$, and even a steeper value of -0.103$\pm$0.009 dex Gyr$^{-1}$ for the thick disc. These age-metallicity relations along with the prominent differences in age, metallicity, and kinematic behaviours seen from the data, clearly show it is important to consider the distinct formation scenarios of the Galactic disc components in modelling the Milky Way.

Aims. We address the problem of young stellar objects that are found too far away from possible star formation sites. Different mechanisms have been proposed before to explain this unexpected circumstance. The idea of high-velocity protostars is one of these mechanisms, although observational support is not always easy to obtain. We aim to shed light on this issue after the serendipitous discovery of a related stellar system. Methods. Following the inspection of archival infrared data, a peculiar anonymous star was found that apparently heads a long tail that resembles a wake-like feature. We conducted a multiwavelength analysis including photometry, astrometry, and spectroscopy. Together with theoretical physical considerations, this approach provided a reasonable knowledge of the stellar age and kinematic properties, together with compelling indications that the extended feature is indeed the signature of a high-velocity, or runaway, newborn star. Results. Our main result is the discovery of a low-mass young stellar object that fits the concept of a runaway T-Tauri star that was hypothesized several decades ago. In this peculiar star, nicknamed UJT-1, the interaction of the stellar wind with the surrounding medium becomes extreme. Under reasonable assumptions, this unusual degree of interaction has the potential to encode the mass-loss history of the star on timescales of several $10^5$ years

Fast radio bursts (FRBs) have been suggested as an excellent celestial laboratory for testing the zero-mass hypothesis of the photon. In this work, we use the dispersion measure (DM)--redshift measurements of 23 localized FRBs to revisit the photon rest mass $m_{\gamma}$. As an improvement over previous studies, here we take into account the more realistic probability distributions of DMs contributed by the FRB host galaxy and intergalactic medium (IGM) from the IllustrisTNG simulation. To better account for the systematic uncertainty induced by the choices of priors of cosmological parameters, we also combine the FRB data with the cosmic microwave background data, the baryon acoustic oscillation data, and type Ia supernova data to constrain the cosmological parameters and $m_{\gamma}$ simultaneously. We derive a new upper limit of $m_{\gamma}\le3.8\times 10^{-51}\;\rm{kg}$, or equivalently $m_{\gamma}\le2.1 \times 10^{-15} \, \rm{eV/c^2}$ ($m_{\gamma} \le 7.2 \times 10^{-51} \, \rm{kg}$, or equivalently $m_{\gamma}\le4.0 \times 10^{-15} \, \rm{eV/c^2}$) at $1\sigma$ ($2\sigma$) confidence level. Meanwhile, our analysis can also lead to a reasonable estimation for the IGM baryon fraction $f_{\rm IGM}=0.873^{+0.061}_{-0.050}$. With the number increment of localized FRBs, the constraints on both $m_{\gamma}$ and $f_{\rm IGM}$ will be further improved. A caveat of constraining $m_{\gamma}$ within the context of the standard $\Lambda$CDM cosmological model is also discussed.

BE Cet is a young solar analog with an age of 0.6 Gyr and a rotation period of 7.655 days. According to chromospheric and photospheric indices, its activity is higher than the solar one. An analysis of photometric data on the time interval between 1977 and 2019 shows the presence of only 6.76 yr cyclic variations in the mean brightness with an amplitude of 0.02 mag. The obtained cycle is 1-2 yr shorter in comparison with the chromospheric cycle determined earlier, whose length was estimated to be 9 or 7.6 yr. Parameters of the cycle, its amplitude and duration change slightly in different epochs. The short-term light variations due to rotational modulation occur with an increase in amplitude up to 0.05 mag near the activity cycle minimum and a decrease in its maximum. Some events of a rapid increase in brightness of 0.2-0.6 mag may be considered as flares.

Aims. With the recent preliminary release of the LOFAR LBA Sky Survey (LoLSS), the first wide-area, ultra-low frequency observations from LOFAR were published. Our aim is to combine this data set with other surveys at higher frequencies to study the spectral properties of a large sample of radio sources. Methods. We present a new cross-matching algorithm taking into account the sizes of the radio sources and apply it to the LoLSS-PR, LoTSS-DR1, LoTSS-DR2 (all LOFAR), TGSS-ADR1 (GMRT), WENSS (WSRT) and NVSS (VLA) catalogues. We then study the number of matched counterparts for LoLSS radio sources and their spectral properties. Results. We find counterparts for 22 607 (89.5%) LoLSS sources. The remaining 2 640 sources (10.5%) are identified either as an artefact in the LoLSS survey (3.6%) or flagged due to their closeness to bright sources (6.9%). We find an average spectral index of $\alpha = -0.77 \pm 0.18$ between LoLSS and NVSS. Between LoLSS and LoTSS-DR2 we find $\alpha = -0.71 \pm 0.31$. The average spectral index is flux density independent above $S_{54} = 181$ mJy. Comparison of the spectral slopes from LoLSS--LoTSS-DR2 with LoTSS-DR2--NVSS indicates that the probed population of radio sources exhibits evidence for a negative spectral curvature.

We perform 3D hydrodynamical simulations to study recombination and ionization during the common envelope (CE) phase of binary evolution, and develop techniques to track the ionic transitions in time and space. We simulate the interaction of a $2\,M_\odot$ red giant branch primary and a $1\,M_\odot$ companion modeled as a particle. We compare a run employing a tabulated equation of state (EOS) that accounts for ionization and recombination, with a run employing an ideal gas EOS. During the first half of the simulations, $\sim15$ per cent more mass is unbound in the tabulated EOS run due to the release of recombination energy, but by simulation end the difference has become negligible. We explain this as being a consequence of (i) the tabulated EOS run experiences a shallower inspiral and hence smaller orbital energy release at late times because recombination energy release expands the envelope and reduces drag, and (ii) collision and mixing between expanding envelope gas, ejecta and circumstellar ambient gas assists in unbinding the envelope, but does so less efficiently in the tabulated EOS run where some of the energy transferred to bound envelope gas is used for ionization. The rate of mass unbinding is approximately constant in the last half of the simulations and the orbital separation steadily decreases at late times. A simple linear extrapolation predicts a CE phase duration of $\sim2\,\mathrm{yr}$, after which the envelope would be unbound.

Recently, a new class of white dwarfs (dubbed ``slowly cooling WDs'') has been identified in two globular clusters (namely M13 and NGC 6752) showing a horizontal branch (HB) morphology with an extended blue tail. The cooling rate of these WDs is reduced by stable thermonuclear hydrogen burning in their residual envelope, and they are thought to be originated by stars that populate the blue tail of the HB and then skip the asymptotic giant branch phase. Consistently, no evidence of such kind of WDs has been found in M3, a similar cluster with no blue extension of the HB. To further explore this phenomenon, we took advantage of deep photometric data acquired with the Hubble Space Telescope in the near-ultraviolet and investigate the bright portion of the WD cooling sequence in M5, another Galactic globular cluster with HB morphology similar to M3. The normalized WD luminosity function derived in M5 turns out to be impressively similar to that observed in M3, in agreement with the fact that the stellar mass distribution along the HB of these two systems is almost identical. The comparison with theoretical predictions is consistent with the fact that the cooling sequence in this cluster is populated by canonical (fast cooling) WDs. Thus, the results presented in this paper provide further support to the scenario proposing a direct causal connection between the slow cooling WD phenomenon and the horizontal branch morphology of the host stellar cluster.

Solar coronal waves frequently appear as bright disturbances that propagate globally from the eruption center in the solar atmosphere, just like the tsunamis in the ocean on Earth. Theoretically, coronal waves can sweep over the underlying chromosphere and leave an imprint in the form of Moreton wave, due to the enhanced pressure beneath their coronal wavefront. Despite the frequent observations of coronal waves, their counterparts in the chromosphere are rarely detected. Why the chromosphere rarely bears the imprints of solar tsunamis remained a mystery since their discovery three decades ago. To resolve this question, all coronal waves and associated Moreton waves in the last decade have been initially surveyed, though the detection of Moreton waves could be hampered by utilising the low-quality H$\alpha$ data from Global Oscillations Network Group. Here, we present 8 cases (including 5 in Appendix) of the coexistence of coronal and Moreton waves in inclined eruptions where it is argued that the extreme inclination is key to providing an answer to address the question. For all these events, the lowest part of the coronal wavefront near the solar surface appears very bright, and the simultaneous disturbances in the solar transition region and the chromosphere predominantly occur beneath the bright segment. Therefore, evidenced by observations, we propose a scenario for the excitation mechanism of the coronal-Moreton waves in highly inclined eruptions, in which the lowest part of a coronal wave can effectively disturb the chromosphere even for a weak (e.g., B-class) solar flare.

Context. Solar extreme ultraviolet (EUV) waves are propagating disturbances in the corona, and they usually accompany with various solar eruptions, from large-scale coronal mass ejections to small-scale coronal jets. Aims. Generally, it is believed that EUV waves are driven by the rapid expansion of coronal loops overlying the erupting cores. In this Letter, we present an exception of EUV wave that was not triggered by the expansion of coronal loops overlying the erupting core. Methods. Combining the multiwavelength observations from multiple instruments, we studied the event in detail. Results. The eruption was restricted in the active region (AR) and disturbed the nearby sheared arcades (SAs) connecting the source AR to a remote AR. Interestingly, following the disturbance, an EUV wave formed close to the SAs, but far away from the eruption source. Conclusions. All the results showed that the EUV wave had a closer temporal and spatial relationship with the disappearing part of SAs than the confined eruption. Hence, we suggest that the EUV wave was likely triggered by the expansion of some strands of SAs, rather than the expansion of erupting loops. It can be a possible complement for the driving mechanisms of EUV waves.

In the beginning of 2023 the Be transient X-ray pulsar RX J0440.9+4431 underwent a fist-ever giant outburst observed from the source peaking in the beginning of February and reaching peak luminosity of $\sim 4.3\times10^{37}$ erg s$^{-1}$. Here we present the results of a detailed spectral and temporal study of the source based on NuSTAR, INTEGRAL, Swift, and NICER observations performed during this period and covering wide range of energies and luminosities. We find that both the pulse profile shape and spectral hardness change abruptly around $\sim2.8\times10^{37}$ erg s$^{-1}$, which we associate with a transition to super-critical accretion regime and erection of the accretion column. The observed pulsed fraction decreases gradually with energy up to 20 keV (with a local minimum around fluorescence iron line), which is unusual for an X-ray pulsar, and then rises rapidly at higher energies with the pulsations significantly detected up to $\sim120$ keV. The broadband energy spectra of RX J0440.9+4431 at different luminosity states can be approximated with a two-hump model with peaks at energies of about 10-20 and 50-70 keV previously suggested for other pulsars without additional features. In particular an absorption feature around 30 keV previously reported and interpreted as a cyclotron line in the literature appears to be absent when using this model, so the question regarding the magnetic field strength of the neutron star remains open. Instead, we attempted to estimate field using several indirect methods and conclude that all of them point to a relatively strong field of around $B\sim 10^{13}$ G.

Rotating massive stars with initial progenitor masses $M_{\rm prog} \sim$ 25 $M_{\odot}$ -- $\sim$140 $M_{\odot}$ can leave rapidly rotating black holes to become collapsars. The black holes and the surrounding accretion disks may develop powerful jets by magneto-hydrodynamics instabilities. The propagation of the jet in the stellar envelope provides the necessary shock heating for triggering nucleosynthesis unseen in canonical core-collapse supernovae. Yet, the energy budget of the jet and its effects on the final chemical abundance pattern are unclear. In this exploratory work, we present a survey on the parameter dependence of collapsar nucleosynthesis on jet energetics. We use the zero-metallicity star with $M_{\rm prog} \sim$ 40 $M_{\odot}$ as the progenitor. The parameters include the jet duration, its energy deposition rate, deposited energy, and the opening angle. We examine the correlations of following observables: (1) the ejecta and remnant masses, (2) the energy deposition efficiency, (3) the $^{56}$Ni production and its correlation with the ejecta velocity, deposited energy, and the ejected mass, (4) the Sc-Ti-V correlation as observed in metal-poor stars, and (5) the [Zn/Fe] ratio as observed in some metal-poor stars. We also provide the chemical abundance table of these explosion models for the use of the galactic chemical evolution and stellar archaeology.

The aim of this work is to verify the stability of the proposed orbital solutions for the third moonlet (Delta) taking into account a realistic gravitational potential for the central body of the quadruple system (Alpha). We also aim to estimate the location and size of a stability region inside the orbit of Gamma. First, we created a set of test particles with intervals of semi-major axis, eccentricities, and inclinations that covers the region interior to the orbit of Gamma, including the proposed orbit of Delta and a wide region around it. We considered three different models for the gravitational potential of Alpha: irregular polyhedron, ellipsoidal body and oblate body. For a second scenario, Delta was considered a massive spherical body and Alpha an irregular polyhedron. Beta and Gamma were assumed as spherical massive bodies in both scenarios. The simulations showed that a large region of space is almost fully stable only when Alpha was modeled as simply as an oblate body. For the scenario with Delta as a massive body, the results did not change from those as massless particles. Beta and Gamma do not play any relevant role in the dynamics of particles interior to the orbit of Gamma. Delta's predicted orbital elements are fully unstable and far from the nearest stable region. The primary instability source is Alpha's elongated shape. Therefore, in the determination of the orbital elements of Delta, it must be taken into account the gravitational potential of Alpha assuming, at least, an ellipsoidal shape.

The $^{12}$CO rotational lines in protoplanetary discs are good tracers of the total spatial extension of the gas component, and potentially planet-disc interactions. We present ALMA long baseline observations of the $^{12}$CO(2-1) line of ten protoplanetary discs from the Ophiuchus DIsc Survey Employing ALMA (ODISEA) project, aiming to set constraints on the gas distribution of these sources. The position angle of the gaseous disc can be inferred for five sources using high-velocity channels, which trace the gas in the inner part of the disc. We compare the high-velocity PAs to the orientations inferred from the continuum, representative of the orientation over $\sim$ 53 to 256 au in these resolved discs. We find a significant difference in orientation for DoAr 44, which is evidence of a tilted inner disc. Eight discs show evidence of gas inside inner dust cavities or gaps, and the disc of ISO-Oph 196 is not detected in $^{12}$CO(2-1), except for the compact signal located inside its dust cavity. Our observations also point out a possible outflow in WLY 2-63.

For decades, the application of muon tomography to spent nuclear fuel (SNF) cask imaging has been theoretically evaluated and experimentally verified by many research groups around the world, including Los Alamos National Laboratory in the United States, Canadian Nuclear Laboratory in Canada, the National Institute for Nuclear Physics in Italy, and Toshiba in Japan. Although monitoring of SNF using cosmic ray muons has attracted significant attention as a promising nontraditional nondestructive radiographic technique, the wide application of muon tomography is often limited because of the natural low cosmic ray muon flux at sea level: 100 m-2min-1sr-1. Recent studies suggest measuring muon momentum in muon scattering tomography (MST) applications to address this challenge. Some techniques have been discussed; however, an imaging algorithm for momentum-coupled MST had not been developed. This paper presents a new imaging algorithm for MST which integrates muon scattering angle and momentum in a single M-value. To develop a relationship between muon momentum and scattering angle distribution, various material samples (Al, Fe, Pb, and U) were thoroughly investigated using a Monte Carlo particle transport code GEANT4 simulation. Reconstructed images of an SNF cask using the new algorithm are presented herein to demonstrate the benefit of measuring muon momentum in MST. In this analysis a missing fuel assembly (FA) was located in the dry storage cask.

Scalar polarization modes of gravitational waves, which are often introduced in the context of the viable extension of gravity, have been actively searched. However, couplings of the scalar modes to the matter are strongly constrained by the fifth-force experiments. Thus, the amplitude of scalar polarization in the observed gravitational-wave signal must be significantly suppressed compared to that of the tensor modes. Here, we discuss the implications of the experiments in the solar system on the detectability of scalar modes in gravitational waves from compact binary coalescences, taking into account the whole processes from the generation to the observation of gravitational waves. We first claim that the energy carried by the scalar modes at the generation is, at most, that of the tensor modes from the observed phase evolution of the inspiral gravitational waves. Next, we formulate general gravitational-wave propagation and point out that the energy flux hardly changes through propagation as long as the background changes slowly compared to the wavelength of the propagating waves. Finally, we show that the possible magnitude of scalar polarization modes detected by the ground-based gravitational-wave telescopes is already severely constrained by the existing gravity tests in the solar system.

Recently, a new field of study called fractional cosmology has emerged. It uses fractional calculus to modify the standard derivative equations and change the Friedmann equations. The evolution of cosmic species densities is also affected by the $\mu$ fractional parameter and the age of the Universe $t_0$. This new approach to cosmology modifies the Friedmann equations and allows for a late cosmic acceleration without the need for a dark energy component. This could be a breakthrough in solving longstanding problems in cosmology. By analyzing observational Hubble data and Type Ia supernovae, we have been able to place strict constraints on the fractional and cosmological parameters. Our results suggest that the Universe may be older than previously estimated. We also explore whether fractional cosmology can help resolve the $H_0$ tension.

Neutrino telescopes are gigaton-scale neutrino detectors comprised of individual light-detection units. Though constructed from simple building blocks, they have opened a new window to the Universe and are able to probe center-of-mass energies that are comparable to those of collider experiments. \prometheus{} is a new, open-source simulation tailored for this kind of detector. Our package, which is written in a combination of \texttt{C++} and \texttt{Python} provides a balance of ease of use and performance and allows the user to simulate a neutrino telescope with arbitrary geometry deployed in ice or water. \prometheus{} simulates the neutrino interactions in the volume surrounding the detector, computes the light yield of the hadronic shower and the out-going lepton, propagates the photons in the medium, and records their arrival times and position in user-defined regions. Finally, \prometheus{} events are serialized into a \texttt{parquet} file, which is a compact and interoperational file format that allows prompt access to the events for further analysis.

The Rossby mode (r-mode) perturbations in pulsars as a steady gravitational wave (GW) sources have been explored. The time evolution and the intensity of the emitted GWs in terms of the strain tensor amplitude have been estimated with the approximation of slow rotation adopting the equation of state derived using the Skyrme effective interaction with NRAPR parameter set. The core of the neutron star has been considered to be $\beta$-equilibrated nuclear matter composed of neutrons, protons, electrons and muons, which is surrounded by a solid crust. Calculations have been made for the critical frequencies, the evolution of frequencies and frequency change rates with time as well as the fiducial viscous and gravitational timescales, across a broad spectrum of pulsar masses. Our findings reveal that the r-mode instability region is associated with rotating young and hot pulsars. Furthermore, it is noteworthy that pulsars with low $L$ value emit gravitational radiation and fall within the r-mode instability region if the primary dissipative mechanism is shear viscosity along the crust-core interface boundary layer. The r-mode perturbation amplitude increases because of GW emissions, in contrast to other non-radial perturbations which transport to infinity the star's angular momentum. Thus the presence of these stellar perturbations implies a non-negative rate of change in transfer of rotational angular momentum. This observation suggests that for a glitching pulsar, the GW emission intensity evolves increasingly with time till the angular frequency diminishes to a value that is below a crucial threshold, after which the compact star ceases to emit radiation.

The new geodetic satellite LARES 2, cousin of LAGEOS and sharing with it almost the same orbital parameters apart from the inclination, displaced by 180 deg, was launched last year. Its proponents suggest using the sum of the nodes of LAGEOS and of LARES 2 to measure the sum of the Lense-Thirring node precessions independently of the systematic bias caused by the even zonal harmonics of the geopotential, claiming a final $\simeq 0.2$ percent total accuracy. In fact, the actual orbital configurations of the two satellites do not allow one to attain the sought for mutual cancellation of their classical node precessions due to the Earth's quadrupole mass moment, as their sum is still $\simeq 5000$ times larger than the added general relativistic rates. This has important consequences. One is that the current uncertainties in the eccentricities and the inclinations of both satellites do not presently allow the stated accuracy goal to be met, needing improvements of 3-4 orders of magnitude. Furthermore, the imperfect knowledge of the Earth's angular momentum $S$ impacts the uncancelled sum of the node precessions, from 150 to 4900 percent of the relativistic signal depending on the uncertainty assumed in $S$. It is finally remarked that the real breakthrough in reliably testing the gravitomagnetic field of the Earth would consist in modeling it and simultaneously estimating one or more dedicated parameter(s) along with other ones characterising the geopotential, as is customarily performed for any other dynamical feature.

We investigate Fermi gases at finite temperature for which the in-medium effective mass may not be constant as a function of the density, the temperature, or the chemical potential. We suggest a formalism that separates the terms for which the mass is constant from the terms which explicitly treat the correction due to the in-medium effective mass. We employ the ensemble equivalence in infinite matter in order to treat these different terms. Our formalism is applied in nuclear matter and we show its goodness by comparing it to an exact treatment based on the numerical calculation of the Fermi integrals.