We used two XMM-Newton and six Neutron Star Interior Composition Explorer (NICER) observations to investigate the fractional rms amplitude of the millihertz quasi-periodic oscillations (mHz QPOs) in the neutron-star low-mass X-ray binary 4U 1636-53. We studied, for the first time, the fractional rms amplitude of the mHz QPOs vs. energy in 4U 1636-53 down to 0.2 keV. We find that, as the energy increases from 0.2 keV to 3 keV, the rms amplitude of the mHz QPOs increases, different from the decreasing trend that has been previously observed above 3 keV. This finding has not yet been predicted by any current theoretical model, however, it provides an important observational feature to speculate whether a newly discovered mHz oscillation originates from the marginally stable nuclear burning process on the neutron star surface.

There is a persistent $H_0$-tension, now at more than $\gtrsim 4\sigma$ level, between the local distance ladder value and the \emph{Planck} cosmic microwave background measurement, in the context of flat $\Lambda$CDM model. We reconstruct $H(z)$ in a cosmological-model-independent way using three low-redshift distance probes including the latest data from baryon acoustic oscillation, Type Ia supernova and four gravitational lensing Time-Delay observations. We adopt general parametric models of $H(z)$ and assume a Gaussian prior on the sound horizon at drag epoch, $r_{\mathrm s}$, from \emph{Planck} measurement. The reconstructed $H_0$ using Pantheon SN Ia and BAO data are consistent with the \emph{Planck} flat $\Lambda$CDM value. When including the GLTD data, $H_0$ increases mildly, yet remaining discrepant with the local measurement at $\sim 2.5\sigma$ level. Our reconstructions being blind to the dark sectors at low redshift, we reaffirm the earlier claims that the Hubble tension is not likely to be solved by modifying the energy budget of the low-redshift universe. We further forecast the constraining ability of future realistic mock BAO data from DESI and GLTD data from LSST, combining which, we anticipate that the uncertainty of the inferred $H_0$ would be improved by $\sim 38\%$, reaching $\sigma_{H_0} \approx 0.56$ uncertainty level.

We investigate the frequency and amplitude of the millihertz quasi-periodic oscillations (mHz QPOs) in the neutron-star low-mass X-ray binary (NS LMXB) 4U 1636-53 using Rossi X-ray Timing Explorer observations. We find that no mHz QPOs appear when the source is in the hard spectral state. We also find that there is no significant correlation between the frequency and the fractional rms amplitude of the mHz QPOs. Notwithstanding, for the first time, we find that the absolute RMS amplitude of the mHz QPOs is insensitive to the parameter Sa, which measures the position of the source in the colour-colour diagram and is usually assumed to be an increasing function of mass accretion rate. This finding indicates that the transition from marginally stable burning to stable burning or unstable burning could happen very rapidly since, before the transition, the mHz QPOs do not gradually decay as the rate further changes.

We report the discovery of millihertz quasi-periodic oscillations (mHz QPOs) from the bursting, high-inclination atoll neutron-star low-mass X-ray binary (NS LMXB) EXO 0748--676 with the Rossi X-ray Time Explorer (RXTE). This class of QPO, originally discovered in three NS LMXBs, has been interpreted as a consequence of a special mode of nuclear burning on the NS surface. Using all the RXTE archival observations of the source, we detected significant ($>3\sigma$) mHz QPOs in 11 observations. The frequency of the oscillations was between $\sim$ 5 and $\sim$ 13 mHz. We also found a decrease of the QPO frequency with time in two occasions; in one of these the oscillations disappeared with the onset of an X-ray burst, similar to what was reported in other sources. Our analysis of the X-ray colours revealed that EXO 0748--676 was in a soft spectral state when it exhibited the QPOs. This makes EXO 0748--676 the sixth source with mHz oscillations associated to marginally stable burning, and the second one that shows a systematic frequency drift. Our results suggest that the mechanism that produces the drift might always be present if the mHz QPOs are observed in the so-called intermediate state.

We analysed an XMM-Newton plus a simultaneous Rossi X-ray Timing Explorer observation and a separate Suzaku observation of the neutron-star low-mass X-ray binary 4U 1728-34. We fitted the X-ray spectra with the self-consistent reflection model relxill. We found that the inclination angle of 4U 1728-34 is 49 degrees, consistent with the upper limit of 60 degrees deduced from the absence of eclipses or dips in this source. The inclination angle in the fit of the XMM-Newton/RXTE observation is larger than 85 degrees, which may be due to the possible calibration issues of the PN instrument in timing mode. We also found that the thermal emission from the accretion disc is not significant. This could be explained either by the relatively high column density of the interstellar medium along the line of sight to the source, which decreases the number of soft disc photons, or if most of the soft thermal photons from the disc are reprocessed in the corona. The ionisation parameter derived from the fits is larger than the value predicted in the framework of the standard reflection model, wherein the disc is irradiated by an X-ray source above the compact object. This inconsistency suggests that irradiation from the neutron star and its boundary layer may play an important role in the ionisation of the accretion disc, and hence in the reflection component in this source.

We investigated the convexity of all type I X-ray bursts with millihertz quasi-periodic oscillations (mHz QPOs) in 4U 1636-53 using archival observations with the Rossi X-ray Timing Explorer. We found that, at a 3.5 $\sigma$ confidence level, in all 39 cases in which the mHz QPOs disappeared at the time of an X-ray burst, the convexity of the burst is positive. The convexity measures the shape of the rising part of the burst light curve and, according to recent models, it is related to the ignition site of bursts on the neutron star surface. This finding suggests that in 4U 1636-53 these 39 bursts and the marginally-stable nuclear burning process responsible for the mHz QPOs take place at the neutron-star equator. This scenario could explain the inconsistency between the high accretion rate required for triggering mHz QPOs in theoretical models and the relatively low accretion rate derived from observations.

We investigate the spectral and timing properties of the millihertz quasi-periodic oscillations (mHz QPOs) in neutron-star low-mass X-ray binary 4U 1636-53 using XMM-Newton and Rossi X-ray Timing Explorer (RXTE) observations. The mHz QPOs in the XMM-Newton/RXTE observations show significant frequency variation and disappear right before type I X-ray bursts. We find no significant correlation between the mHz QPO frequency and the temperature of the neutron-star surface, which is different from theoretical predictions. For the first time we observed the full lifetime of a mHz QPO lasting 19 ks. Besides, we also measure a frequency drift timescale around 15 ks, we speculate that this is the cooling timescale of a layer deeper than the burning depth, possibly heated by the previous burst. Moreover, the analysis of all X-ray bursts in this source shows that all type I X-ray bursts associated with the mHz QPOs are short, bright and energetic, suggesting a potential connection between mHz QPOs and He-rich X-ray bursts.

We detected millihertz quasi-periodic oscillations (QPOs) in an XMM-Newton observation of the neutron-star low-mass X-ray binary 4U 1636-53. These QPOs have been interpreted as marginally-stable burning on the neutron-star surface. At the beginning of the observation the QPO was at around 8 mHz, together with a possible second harmonic. About 12 ks into the observation a type I X-ray burst occurred and the QPO disappeared; the QPO reappeared ~25 ks after the burst and it was present until the end of the observation. We divided the observation into four segments to study the evolution of the spectral properties of the source during intervals with and without mHz QPO. We find that the temperature of the neutron-star surface increases from the QPO segment to the non-QPO segment, and vice versa. We also find a strong correlation between the frequency of the mHz QPO and the temperature of a black-body component in the energy spectrum representing the temperature of neutron-star surface. Our results are consistent with previous results that the frequency of the mHz QPO depends on the variation of the heat flux from the neutron star crust, and therefore supports the suggestion that the observed QPO frequency drifts could be caused by the cooling of deeper layers.

We used six simultaneous XMM-Newton and Rossi X-ray Timing Explorer plus five Suzaku observations to study the continuum spectrum and the iron emission line in the neutron-star low-mass X-ray binary 4U 1636-53. We modelled the spectra with two thermal components (representing the accretion disc and boundary layer), a Comptonised component (representing a hot corona), and either a Gaussian or a relativistic line component to model an iron emission line at about 6.5 keV. For the relativistic line component we used either the diskline, laor or kyrline model, the latter for three different values of the spin parameter. The fitting results for the continuum are consistent with the standard truncated disc scenario. We also find that the flux and equivalent width of the iron line first increase and then decrease as the flux of the Comptonised component increases. This could be explained either by changes in the ionisation state of the accretion disc where the line is produced by reflection, or by light bending of the emission from the Comptonised component if the height at which this component is produced changes with mass accretion rate.