Phase-resolved Spectra Research Articles

Inferring the origins of the pulsed γ-ray emission from the Crab pulsar with ten-year Fermi-LAT data (Corrigendum)

Context: The Crab pulsar is a bright $\gamma$-ray source detected at photon energies up to $\sim$1 TeV. Its phase-averaged and phase-resolved $\gamma$-ray spectra below 10 GeV exhibit exponential cutoffs while those above 10 GeV apparently follow simple power-laws. Aims: We re-visit the $\gamma$-ray properties of the Crab pulsar with 10-year \emph{Fermi} Large Area Telescope (LAT) data in the range of 60 MeV--500 GeV. With the phase-resolved spectra, we investigate the origins and mechanisms responsible for the emissions. Methods: The phaseograms are reconstructed for different energy bands and further analysed using a wavelet decomposition. The phase-resolved energy spectra are combined with the observations of ground-based instruments like MAGIC and VERITAS to achieve a larger energy converage. We fit power-law models to the overlapping energy spectra from 10 GeV to $\sim$1 TeV. We include in the fit a relative cross-calibration of energy scales between air-shower based gamma-ray telescopes with the orbital pair-production telescope of the Fermi mission. Results: We confirm the energy-dependence of the $\gamma$-ray pulse shape, and equivalently, the phase-dependence of the spectral shape for the Crab pulsar. A relatively sharp cutoff at a relatively high energy of $\sim$8 GeV is observed for the bridge-phase emission. The $E>$10 GeV spectrum observed for the second pulse peak is harder than those for other phases. Conclusions: In view of the diversity of phase-resolved spectral shapes of the Crab pulsar, we tentatively propose a multi-origin scenario where the polar-cap, outer-gap and relativistic-wind regions are involved.

Studies of emission regions of the X-ray pulsar Hercules X-1 with pulse-phase-resolved spectra observed with Suzaku

Abstract Pulse-phase-resolved X-ray spectra of the X-ray pulsar Her X-1 observed with Suzaku are analyzed. The analysis is based on a hypothesis that the observed variations of the continuum spectra associated with the neutron star spin originate in the changes of the projection areas of the emission regions viewed from the observer, and that the spectral shape observed from each of the emission regions does not change with the viewing angle variation. We calculate the ratios of the fluxes in each energy bin between arbitrary pairs of the pulse-phase-resolved spectra over the entire energy range. We identify flat portions in these ratio spectra, which are likely to indicate the intrinsic emission components that have invariant spectral shapes, in the three energy bands of below ∼1 keV, 2–6 keV, and above ∼18 keV. Then, we approximately solve the simultaneous equations for the intrinsic spectral components from a pair of the phase-resolved spectra that show flat portions in the ratio spectrum. As a result, spectra of the three components are obtained as the low-, medium-, and high-energy bands corresponding to those for the flat portions. They are found to be well fitted with a soft blackbody, a cut-off power law, and a hard blackbody model, respectively. Finally, we fit all the phase-resolved spectra simultaneously with a composite model of three continuum components together with a cyclotron absorption feature and an iron line, and obtain an acceptable result. The origins of the cut-off power-law and hard blackbody components are discussed in relation to recently developed theoretical models.

NuSTAR observation of GRO J1744–28 at low mass accretion rate

Context. Neutron stars in low-mass X-ray binaries are important systems for studying the physics of accretion onto compact objects. The system GRO J1744–28 is particularly interesting as it usually shows clear pulsations as well as X-ray bursts. Additionally, there are claims for a magnetic field of 5 × 1011 G through the detection of a cyclotron resonant scattering feature (CRSF). Aims. We present the spectral analysis of GRO J1744–28 using ∼29 ks of NuSTAR data taken in 2017 February at a low luminosity of 3.2 × 1036 erg s−1 (3−50 keV). Our goal is to study the variability of the source spectrum with pulse phase and to search for the claimed CRSF. Methods. The continuum spectrum was modeled with an absorbed power law with exponential cutoff, and an additional iron line component. We found no obvious indications for a CRSF, and therefore performed a detailed cyclotron line search using statistical methods. We performed this search on pulse phase-averaged spectra and on phase-resolved spectra. Results. GRO J1744–28 was observed in a low-luminosity state. The previously detected Type II X-ray bursts are absent. Clear pulsations at a period of 2.141124(9) Hz are detected. The pulse profile shows an indication of a secondary peak that was not seen at higher flux. The upper limit for the strength of a CRSF in the 3−20 keV band is 0.07 keV (90% CL), lower than the strength of the line found at higher luminosity. Conclusions. The detection of pulsations shows that the source did not enter the “propeller” regime, even though the source flux of 4.15 × 10−10 erg cm−2 s−1 was almost one order of magnitude below the threshold for the propeller regime claimed in previous studies on this source. The transition into the propeller regime in GRO J1744–28 must therefore be below a luminosity of 3.2 × 1036 erg s−1 (3−50 keV), which implies a surface magnetic field ≲2.9 × 1011 G and mass accretion rate ≲1.7 × 1016 g s−1. A change of the CRSF depth as function of luminosity is not unexpected and has been observed in other sources. This result possibly implies a change in emission geometry as function of mass accretion rate to reduce the depth of the line below our detection limit.

Evolution of timing and spectral characteristics of 4U 1901+03 during its 2019 outburst using the Swift and NuSTAR observatories

ABSTRACT We report the results from a detailed timing and spectral study of a transient X-ray pulsar, 4U 1901+03 during its 2019 outburst. We performed broadband spectroscopy in the 1–70 $\rm keV$ energy band using four observations made with Swift and NuSTAR at different intensity levels. Our timing results reveal the presence of highly variable pulse profiles dependent on both luminosity and energy. Our spectroscopy results showed the presence of a cyclotron resonance scattering feature (CRSF) at ∼ 30 keV. This feature at 30 keV is highly luminosity and pulse phase dependent. Phase-averaged spectra during the last two observations, made close to the declining phase of the outburst, showed the presence of this feature at around $30~\rm {keV}$. The existence of CRSF at 30 keV during these observations is well supported by an abrupt change in the shape of pulse profiles found close to this energy. We also found that 30 keV feature was significantly detected in the pulse phase-resolved spectra of observations made at relatively high luminosities. Moreover, all spectral fit parameters showed a strong pulse phase dependence. In line with the previous findings, an absorption feature at around $10~\rm {keV}$ is significantly observed in the phase-averaged X-ray spectra of all observations and also showed a strong pulse phase dependence.

Phase-resolved Analyses of Millihertz Quasi-periodic Oscillations in 4U 1636-53 using the Hilbert–Huang Transform

Abstract We present phase-resolved spectroscopy based on the Hilbert–Huang transform (HHT) for millihertz quasi-periodic oscillations (mHz QPOs) in 4U 1636-53. The ∼8 mHz QPO can be detected approximately several thousand seconds before a type-I X-ray burst. It has been interpreted as marginally stable burning on the neutron-star surface. In this study, we use the HHT to analyze the data collected by XMM-Newton between 2007 and 2009. The HHT is a powerful tool that enables us to obtain instantaneous frequency, amplitude, and phase of phenomena with nonstationary periodicity, such as QPOs. With well-defined phases, the oscillation profile of the ∼8 mHz QPO for 4U 1636-53 can be precisely revealed. In addition to the oscillation profile, phase-resolved spectra for the complete cycle are constructed. From the correlation between spectral parameters and fluxes, we find that the oscillation is mainly attributed to variations in the area emitting blackbody radiation in three out of four observations with mHz QPO detections, whereas the other one shows a concurrent variation of temperature and flux with a constant emitting area. Although the cause of the difference is not clear, it might be related to the spectral state of the source that can be observed from a hard color difference in the color–color diagram.

Inferring the origins of the pulsed γ-ray emission from the Crab pulsar with ten-year Fermi-LAT data

Context. The Crab pulsar is a bright γ-ray source, which has been detected at photon energies up to ∼1 TeV. Its phase-averaged and phase-resolved γ-ray spectra below 10 GeV exhibit exponential cutoffs, while those above 10 GeV apparently follow simple power laws. Aims. We re-visit the γ-ray properties of the Crab pulsar with ten-year Fermi Large Area Telescope (LAT) data in the range of 60 MeV–500 GeV. With the phase-resolved spectra, we investigate the origins and mechanisms responsible for the emissions. Methods. The phaseograms were reconstructed for different energy bands and further analysed using a wavelet decomposition. The phase-resolved energy spectra were combined with the observations of ground-based instruments, such as MAGIC and VERITAS, to achieve a larger energy converage. We fitted power-law models to the overlapping energy spectra from 10 GeV to ∼1 TeV. In the fit, we included a relative cross-calibration of energy scales between air-shower-based gamma-ray telescopes with the orbital pair-production telescope from the Fermi mission. Results. We confirm the energy-dependence of the γ-ray pulse shape and, equivalently, the phase-dependence of the spectral shape for the Crab pulsar. A relatively sharp cutoff at a relatively high energy of ∼8 GeV is observed for the bridge-phase emission. The E > 10 GeV spectrum observed for the second pulse peak is harder than those for other phases. Conclusions. In view of the diversity of phase-resolved spectral shapes of the Crab pulsar, we tentatively propose a multi-origin scenario where the polar-cap, outer-gap, and relativistic-wind regions are involved.

The unusual behavior of the young X-ray pulsar SXP 1062 during the 2019 outburst

We present the results of the first dedicated observation of the young X-ray pulsar SXP 1062 in the broad X-ray energy band obtained during its 2019 outburst with the NuSTAR and XMM-Newton observatories. The analysis of the pulse-phase averaged and phase-resolved spectra in the energy band from 0.5 to 70 keV did not reveal any evidence for the presence of a cyclotron line. The spin period of the pulsar was found to have decreased to 979.48 ± 0.06 s implying a ∼10% reduction compared to the last measured period during the monitoring campaign conducted about five years ago, and is puzzling considering that the system apparently has not shown major outbursts ever since. The switch of the pulsar to the spin-up regime supports the common assumption that torques acting on the accreting neutron star are nearly balanced and thus SXP 1062 likely also spins with a period close to the equilibrium value for this system. The current monitoring of the source also revealed a sharp drop in its soft X-ray flux right after the outburst, which is in drastic contrast to the behavior during the previous outburst when the pulsar remained observable for years with only a minor flux decrease after the end of the outburst. This unexpected off state of the source lasted for at most 20 days after which SXP 1062 returned to the level observed during previous campaigns. We discuss this and other findings in context of the modern models of accretion onto strongly magnetized neutron stars.

Phase-resolved joint spectra tomography of a ring resonator photon pair source using a silicon photonic chip.

The exponential growth of photonic quantum technologies is driving the demand for tools to measure the quality of their information carriers. One of the most prominent is stimulated emission tomography (SET), which uses classical coherent fields to measure the joint spectral amplitude (JSA) of photon pairs with high speed and resolution. While the modulus of the JSA can be directly addressed from a single intensity measurement, the retrieval of the joint spectral phase (JSP) is far more challenging and received minor attention. However, a wide class of spontaneous sources of technological relevance, as chip integrated micro-resonators, have a JSP with a rich structure that carries correlations hidden in the intensity domain. Here, using a compact and reconfigurable silicon photonic chip, the complex JSA of a micro-ring resonator photon pair source is measured for the first time. The photonic circuit coherently excites the ring and a reference waveguide, and the interferogram formed by their stimulated fields is used to map the ring JSP through a novel phase reconstruction technique. This tool complements the traditionally bulky and sophisticated methods implemented so far, simultaneously minimizing the set of required resources.

A General Approach To Combine the Advantages of Collinear and Noncollinear Spectrometer Designs in Phase-Resolved Second-Order Nonlinear Spectroscopy.

Recent years have seen a huge progress in the development of phase-sensitive second-order laser spectroscopy which has proven to be a very powerful tool for the investigation of interfaces. In these techniques, the nonlinear interaction between two short laser pulses and the sample yields a signal pulse which subsequently interferes with a third pulse, the so-called local oscillator. To obtain accurate phase information, the relative phases between the signal and local oscillator pulses must be stabilized and their timings precisely controlled. Despite much progress made, fulfilling both requirements remains a formidable experimental challenge. The two common approaches employ different beam geometries which each yields its particular advantages and deficiencies. While noncollinear spectrometers allow for a relatively simple timing control they typically yield poor phase stability and require a challenging alignment. Collinear approaches in contrast come with a simplified alignment and improved phase stability but typically suffer from a highly limited timing control. In this contribution we present a general experimental solution which allows for combining the advantages of both approaches while being compatible with most of the common spectrometer types. On the basis of a collinear geometry, we exploit different selected polarization states of the light pulses in well-defined places in the spectrometer to achieve a precise timing control. The combination of this technique with a balanced detection scheme allows for the acquisition of highly accurate phase-resolved nonlinear spectra without any loss in experimental flexibility. In fact, we show that the implementation of this technique allows us to employ advanced pulse timing schemes inside the spectrometer, which can be used to suppress nonlinear background signals and extend the capabilities of our spectrometer to measure phase-resolved sum frequency spectra of interfaces in a liquid cell.

A diagnostic of the orbital spectrum of LS 5039 with Fermi-LAT

LS 5039 is a well-known γ-ray binary system which consists of an unknown compact object and a massive companion O star. It shows rather stable emissions at high energies over years and hence serves as an ideal laboratory to investigate the emission mechanism for such peculiar systems which emit prominent γ-rays. To this end, we take the orbital phase resolved energy spectrum as observed by Fermi-LAT over 10 years. We divide the orbit into four orbital phases, each with an orbital phase range of 0.25, centered at 0.00, 0.25, 0.50 and 0.75 respectively, where the phase 0.0 is the periastron and phase 0.5 is the apastron. The phases around 0.25 and 0.75 are symmetric and hence are supposed to have identical local acceleration environment. The spectral analysis shows that the Fermi-LAT spectra are largely different from these two symmetric orbital phases: the emission from orbital phase 0.25 turns out to be significantly stronger than that from 0.75. This result does not fit a scenario that γ-rays are Doppler boosted emission from bow shock tails if LS 5039 has a shock configuration similar to PSR B1259-63, and indicates that the inverse Compton scatterings between the shock accelerated plasma and the stellar particle environment is the underline procedure. We also find that the previous report for a disappearance of the orbital modulation at 3–20 GeV is due to the similar spectral turn-over energies of the different orbital phases. The spectral properties of periastron and apastron regions are addressed in the context of the measurements in phase regions around 0.25 and 0.75.

NUSTAR and Swift observations of AMXP Swift J1756.9−2508 during its 2018 outburst

ABSTRACT We present here the timing and spectral analysis of the accreting millisecond pulsar (AMXP) SwiftJ1756.9−2508 during its recent outburst in 2018 using Swift and NUSTAR observations. The simultaneous fitting of the Swift and NUSTAR spectra indicates that the source was in the hard state with a cut-off energy of about 74.58 keV. We also study in detail the pulse profile of the AMXP and its dependence on energy. The colour–colour diagram of the source is different from those previously reported. We performed phase- and time-resolved spectral analysis using NUSTAR data. Pulse phase-resolved spectra were fitted with a power-law model and significant changes in the spectral parameters with pulse phase were observed. The orbital phase and time-resolved spectra were fitted with a cut-off power-law model. The column density and photon index obtained from orbital phase spectral analysis were found to show some anticorrelation with the flux. Through time-resolved spectral analysis, we observed that the spectral parameters show positive correlation with each other and with the flux. We do not observe a softening of the spectrum with time. No emission lines or Compton bump were observed in the spectrum of the AMXP.

Temporal and spectral X-ray properties of magnetar SGR 1900+14 derived from observations with NuSTAR and XMM-Newton

Abstract X-ray observations play a crucial role in understanding the emission mechanism and relevant physical phenomena of magnetars. We report on X-ray observations made in 2016 of a young magnetar, SGR 1900+14, which is famous for a giant flare in 1998 August. Simultaneous observations were conducted with XMM-Newton and NuSTAR on 2016 October 20 with 23 and 123 ks exposures, respectively. The NuSTAR hard X-ray coverage enabled us to detect the source up to 70 keV. The 1–10 keV and 15–60 keV fluxes were $3.11(3)\times 10^{-12} \, {\rm erg \, s^{-1} \, cm^{-2}}$ and $6.8(3)\times 10^{-12} \, {\rm erg \, s^{-1} \, cm^{-2}}$, respectively. The 1–70 keV spectra were fitted well by a blackbody plus power-law model with a surface temperature of $kT=0.52(2) \, {\rm keV}$, a photon index of the hard power-law of Γ = 1.21(6), and a column density of $N_{\,\rm H}=1.96(11)\times 10^{22} \, {\rm cm^{-2}}$. Compared with previous observations with Suzaku in 2006 and 2009, the 1–10 keV flux showed a decrease by 25%–40%, while the spectral shape did not show any significant change with differences of kT and NH being within 10% of each other. Through timing analysis, we found that the rotation period of SGR 1900+14 on 2016 October 20 was $5.22669(3) \, {\rm s}$. The long-term evolution of the rotation period shows a monotonic decrease in the spin-down rate $\dot{P}$ lasting for more than 15 yr. We also found characteristic behavior of the hard-tail power-law component of SGR 1900+14. The energy-dependent pulse profiles vary in morphology with a boundary of 10 keV. The phase-resolved spectra show the differences between photon indices (Γ = 1.02–1.44) as a function of the pulse phase. Furthermore, the photon index is positively correlated with the X-ray flux of the hard power-law component, which could not be resolved by the previous hard X-ray observations.

Effects of Compton scattering on the neutron star radius constraints in rotation-powered millisecond pulsars

The aim of this work is to study the possible effects and biases on the radius constraints for rotation-powered millisecond pulsars when using Thomson approximation to describe electron scattering in the atmosphere models, instead of using exact formulation for Compton scattering. We compare the differences between the two models in the energy spectrum and angular distribution of the emitted radiation. We also analyse a self-generated, synthetic, phase-resolved energy spectrum, based on Compton atmosphere and the most X-ray luminous, rotation-powered millisecond pulsars observed by the Neutron star Interior Composition ExploreR (NICER). We derive constraints for the neutron star parameters using both the Compton and Thomson models. The results show that the method works by reproducing the correct parameters with the Compton model. However, biases are found in both the size and the temperature of the emitting hotspot, when using the Thomson model. The constraints on the radius are still not significantly changed, and therefore the Thomson model seems to be adequate if we are interested only in the radius measurements using NICER.

Temporal evolution of electron density and temperature in low pressure transient Ar/N2 plasmas estimated by optical emission spectroscopy

A recently published method for the analysis of phase-resolved optical emission spectra was extended in order to permit estimation of time-resolved electron density profiles. The previously presented method combined collisional-radiative modelling with a self-absorption method to estimate the evolution of Te with sub-cycle time-resolution. However, it was not capable to give similar profiles for ne as the model was insensitive to its variations. The extensions proposed in this work describe a way to also estimate the electron density with sub-cycle time resolution from the changing rates of the argon Paschen 1s states. The method was applied to a low-pressure DBD-jet operated with argon and several argon–nitrogen mixtures with up to 4% N2. Good agreement among evaluation of ne from changing rates of individual 1s states was observed during the collisional phase and the full-cycle temporal profile could be calculated from relative changes in light emission. Electron densities exhibited a drop for larger admixtures of nitrogen and ranged from 1017 m−3 to 1018 m−3. As assumed in a previous work, the electron temperature model worked without explicit consideration of additional processes even when N2 affected the plasma. However, presumably due to collisional quenching by nitrogen, two argon Paschen 2p levels were found to be inappropriate for Te estimation and had to be removed. Values for electron temperature from the remaining levels remained at a similar value as for pure argon.

Insight-HXMT observations of the Crab pulsar

We report on X-ray emission properties of the Crab pulsar (PSR B0531+21) using observations by insight-HXMT during its first year after launch. We obtained high signal-to-noise profiles in the X-ray energy band 11–250 keV. We have confirmed an increase in the flux ratio of the second peak over the main peak with increasing energy, consistent with other missions. The separation of the two peaks shows no significant trend with increasing energy. The phase-averaged spectrum, fitted by a logpar model, and the phase-resolved spectra of the Crab pulsar, fitted by a powerlaw in the different energy bands of HXMT, are consistent with RXTE and NuSTAR in that photon indices evolve as a function of phase as well as a function of energy, contributing to a broadband modeling.

Climate of an ultra hot Jupiter

We present the analysis of a full-orbit, spectroscopic phase curve of the ultra hot Jupiter (UHJ) WASP-18b, obtained with the Wide Field Camera 3 aboard the Hubble Space Telescope. We measured the normalised day-night contrast of the planet as >0.96 in luminosity: the disc-integrated dayside emission from the planet is at 964 ± 25 ppm, corresponding to 2894 ± 30 K, and we place an upper limit on the nightside emission of <32 ppm or 1430 K at the 3σ level. We also find that the peak of the phase curve exhibits a small, but significant offset in brightness of 4.5 ± 0.5° eastward. We compare the extracted phase curve and phase-resolved spectra to 3D global circulation models and find that broadly the data can be well reproduced by some of these models. We find from this comparison several constraints on the atmospheric properties of the planet. Firstly we find that we need efficient drag to explain the very inefficient day-night recirculation observed. We demonstrate that this drag could be due to Lorentz-force drag by a magnetic field as weak as 10 gauss. Secondly, we show that a high metallicity is not required to match the large day-night temperature contrast. In fact, the effect of metallicity on the phase curve is different from cooler gas-giant counterparts because of the high-temperature chemistry in the atmosphere of WASP-18b. Additionally, we compared the current UHJ spectroscopic phase curves, WASP-18b and WASP-103b, and show that these two planets provide a consistent picture with remarkable similarities in their measured and inferred properties. However, key differences in these properties, such as their brightness offsets and radius anomalies, suggest that UHJ could be used to separate between competing theories for the inflation of gas-giant planets.

Thermal Emission and Magnetic Beaming in the Radio and X-Ray Mode-switching PSR B0943+10

PSR B0943+10 is a mode-switching radio pulsar characterized by two emission modes with different radio and X-ray properties. Previous studies, based on simple combinations of blackbody and power-law models, showed that its X-ray flux can be decomposed in a pulsed thermal plus an unpulsed nonthermal components. However, if PSR B0943+10 is a nearly aligned rotator seen pole-on, as suggested by the radio data, it is difficult to reproduce the high observed pulsed fraction unless magnetic beaming is included. In this work, we reanalyze all of the available X-ray observations of PSR B0943+10 with simultaneous radio coverage, modeling its thermal emission with polar caps covered by a magnetized hydrogen atmosphere or with a condensed iron surface. The condensed surface model provides good fits to the spectra of both pulsar modes, but, similarly to the blackbody, it cannot reproduce the observed pulse profiles, unless an additional power law with an ad hoc modulation is added. Instead, the pulse profiles and phase-resolved spectra are well described using the hydrogen atmosphere model to describe the polar cap emission plus an unpulsed power law. For the X-ray brighter state (Q-mode) we obtain a best fit with a temperature , an emitting radius m, a magnetic field consistent with the value of the dipole field of 4 × 1012 G inferred from the timing parameters, and a small angle between the magnetic and spin axis, ξ = 5 The corresponding parameters for the X-ray fainter state (B-mode) are and R ∼ 170 m.

Pulsar gamma-ray emission in the radiation reaction regime

Since the era of the Fermi/LAT and atmospheric Cerenkov telescopes, pulsars are known to emit high and very high-energy photons, in the MeV-GeV range and sometimes up to TeV. To date, it is still unclear where and how these photons are produced. Nevertheless gamma-ray photons require particle acceleration to ultra-relativistic speeds. In this paper, we compute single particle trajectories for leptons in an arbitrary strong electromagnetic field in the so-called radiation reaction limit. In this picture, particle velocity only depends on the local electromagnetic field which we assume to follow the vacuum dipole rotator. From this velocity field, we compute the curvature radiation spectrum and light-curves. Sky maps and phase-resolved spectra are then deduced accounting for realistic pulsar periods and magnetic field strengths. Emission sites within the pulsar magnetosphere where most of radiation emanates are then localized. For standard parameters of millisecond and normal pulsars, we show that a break in the spectrum occurs at several GeV in agreement with the Fermi/LAT second pulsar catalogue. A sample of representative phase-resolved spectra and sky-maps are shown. A pair multiplicity of several tenths to several thousands is required to account for the total gamma-ray luminosity. Moreover depending on the geometry, single or double-peaked light-curves are found. Our model shows that minimalist assumptions are already able to reproduce salient features of pulsar emission.

NICER Discovers mHz Oscillations in the “Clocked” Burster GS 1826−238

Abstract We report the discovery with the Neutron Star Interior Composition Explorer (NICER) of mHz X-ray brightness oscillations from the “clocked burster” GS 1826−238. NICER observed the source in the periods 2017 June 20–29, July 11–13, and September 9–15, for a total useful exposure of 34 ks. Two consecutive dwells obtained on 2017 September 9 revealed highly significant oscillations at a frequency of 8 mHz. The fractional, sinusoidal modulation amplitude increases from 0.7% at 1 keV to ≈2% at 6 keV. Similar oscillations were also detected at lower significance in three additional dwells. The oscillation frequency and amplitude are consistent with those of mHz QPOs reported in other accreting neutron star systems. A thermonuclear X-ray burst was also observed on 2017 June 22. The burst properties and X-ray colors are both consistent with GS 1826 being in a soft spectral state during these observations, findings that are confirmed by ongoing monitoring with MAXI and SWIFT-BAT. Assuming that the mHz oscillations are associated with blackbody emission from the neutron star surface, modeling of the phase-resolved spectra shows that the oscillation is consistent with being produced by modulation of the temperature component of this emission. In this interpretation, the blackbody normalization, proportional to the emitting surface area, is consistent with being constant through the oscillation cycle. We place the observations in the context of the current theory of marginally stable burning and briefly discuss the potential for constraining neutron star properties using mHz oscillations.

Global Climate and Atmospheric Composition of the Ultra-hot Jupiter WASP-103b from HST and Spitzer Phase Curve Observations

Abstract We present thermal phase curve measurements for the hot Jupiter WASP-103b observed with Hubble/WFC3 and Spitzer/IRAC. The phase curves have large amplitudes and negligible hotspot offsets, indicative of poor heat redistribution to the nightside. We fit the phase variation with a range of climate maps and find that a spherical harmonics model generally provides the best fit. The phase-resolved spectra are consistent with blackbodies in the WFC3 bandpass, with brightness temperatures ranging from 1880 ± 40 K on the nightside to 2930 ± 40 K on the dayside. The dayside spectrum has a significantly higher brightness temperature in the Spitzer bands, likely due to CO emission and a thermal inversion. The inversion is not present on the nightside. We retrieved the atmospheric composition and found that it is moderately metal-enriched ( ) and the carbon-to-oxygen ratio is below 0.9 at 3σ confidence. In contrast to cooler hot Jupiters, we do not detect spectral features from water, which we attribute to partial H2O dissociation. We compare the phase curves to 3D general circulation models and find that magnetic drag effects are needed to match the data. We also compare the WASP-103b spectra to brown dwarfs and young, directly imaged companions. We find that these objects have significantly larger water features, indicating that surface gravity and irradiation environment play an important role in shaping the spectra of hot Jupiters. These results highlight the 3D structure of exoplanet atmospheres and illustrate the importance of phase curve observations for understanding their complex chemistry and physics.