skip to main content

Title: Measuring the Nonaxially Symmetric Surface Temperature Distribution of the Central Compact Object in Puppis A
Abstract The surface temperature distributions of central compact objects (CCOs) are powerful probes of their crustal magnetic field strengths and geometries. Here we model the surface temperature distribution of RX J0822−4300, the CCO in the Puppis A supernova remnant, using 471 ks of XMM-Newton data. We compute the energy-dependent pulse profiles in 16 energy bands, fully including the general relativistic effects of gravitational redshift and light bending, to accurately model the two heated surface regions of different temperatures and areas, in addition to constraining the viewing geometry. This results in precise measurements of the two temperatures: kT warm = ( 1 + z ) × 0.222 − 0.019 + 0.018 keV and kT hot = (1 + z ) × 0.411 ± 0.011 keV. The two heated surface regions are likely located very close to the rotational poles, with the most probable position of the hotter component ≈ 6° from the rotational pole. For the first time, we are able to measure a deviation from a pure antipodal hot-spot geometry, with a longitudinal offset δ γ = 11 .° 7 − 2 .° 5 + 2 .° 6 . The discovery of this asymmetry, along with the factor of ≈2 more » temperature difference between the two emitting regions, may indicate that RX J0822−4300 was born with a strong, tangled crustal magnetic field. « less
; ; ;
Award ID(s):
Publication Date:
Journal Name:
The Astrophysical Journal
Page Range or eLocation-ID:
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Magnetars, isolated neutron stars with magnetic-field strengths typically ≳10 14 G, exhibit distinctive months-long outburst epochs during which strong evolution of soft X-ray pulse profiles, along with nonthermal magnetospheric emission components, is often observed. Using near-daily NICER observations of the magnetar SGR 1830-0645 during the first 37 days of a recent outburst decay, a pulse peak migration in phase is clearly observed, transforming the pulse shape from an initially triple-peaked to a single-peaked profile. Such peak merging has not been seen before for a magnetar. Our high-resolution phase-resolved spectroscopic analysis reveals no significant evolution of temperature despite the complex initial pulse shape, yet the inferred surface hot spots shrink during peak migration and outburst decay. We suggest two possible origins for this evolution. For internal heating of the surface, tectonic motion of the crust may be its underlying cause. The inferred speed of this crustal motion is ≲100 m day −1 , constraining the density of the driving region to ρ ∼ 10 10 g cm −3 , at a depth of ∼200 m. Alternatively, the hot spots could be heated by particle bombardment from a twisted magnetosphere possessing flux tubes or ropes, somewhat resembling solar coronal loops, thatmore »untwist and dissipate on the 30–40 day timescale. The peak migration may then be due to a combination of field-line footpoint motion (necessarily driven by crustal motion) and evolving surface radiation beaming. This novel data set paints a vivid picture of the dynamics associated with magnetar outbursts, yet it also highlights the need for a more generic theoretical picture where magnetosphere and crust are considered in tandem.« less
  2. Abstract We report on NICER X-ray monitoring of the magnetar SGR 1830−0645 covering 223 days following its 2020 October outburst, as well as Chandra and radio observations. We present the most accurate spin ephemerides of the source so far: ν = 0.096008680(2) Hz, ν ̇ = − 6.2 ( 1 ) × 10 − 14 Hz s −1 , and significant second and third frequency derivative terms indicative of nonnegligible timing noise. The phase-averaged 0.8–7 keV spectrum is well fit with a double-blackbody (BB) model throughout the campaign. The BB temperatures remain constant at 0.46 and 1.2 keV. The areas and flux of each component decreased by a factor of 6, initially through a steep decay trend lasting about 46 days, followed by a shallow long-term one. The pulse shape in the same energy range is initially complex, exhibiting three distinct peaks, yet with clear continuous evolution throughout the outburst toward a simpler, single-pulse shape. The rms pulsed fraction is high and increases from about 40% to 50%. We find no dependence of pulse shape or fraction on energy. These results suggest that multiple hot spots, possibly possessing temperature gradients, emerged at outburst onset and shrank as the outburst decayed.more »We detect 84 faint bursts with NICER, having a strong preference for occurring close to the surface emission pulse maximum—the first time this phenomenon is detected in such a large burst sample. This likely implies a very low altitude for the burst emission region and a triggering mechanism connected to the surface active zone. Finally, our radio observations at several epochs and multiple frequencies reveal no evidence of pulsed or burst-like radio emission.« less
  3. ABSTRACT We present results from a deep (174 ks) Chandra observation of the FR-II radio galaxy 3C 220.1, the central brightest cluster galaxy (BCG) of a kT ∼ 4 keV cluster at z = 0.61. The temperature of the hot cluster medium drops from ∼5.9 to ∼3.9 keV at ∼35 kpc radius, while the temperature at smaller radii may be substantially lower. The central active galactic nucleus (AGN) outshines the whole cluster in X-rays, with a bolometric luminosity of 2.0 × 1046 erg s−1 (∼10 per cent of the Eddington rate). The system shows a pair of potential X-ray cavities ∼35 kpc east and west of the nucleus. The cavity power is estimated within the range of 1.0 × 1044 and 1.7 × 1045 erg s−1, from different methods. The X-ray enhancements in the radio lobes could be due to inverse Compton emission, with a total 2–10 keV luminosity of ∼8.0 × 1042 erg s−1. We compare 3C 220.1 with other cluster BCGs, including Cygnus A, as there are few BCGs in rich clusters hosting an FR-II galaxy. We also summarize the jet power of FR-II galaxies from different methods. The comparison suggests that the cavity power of FR-II galaxies likely underestimates the jet power. The properties of 3C 220.1 suggest that it is at the transition stage from quasar-mode feedback to radio-mode feedback.

    We present investigations of rapidly rotating convection in a thick spherical shell geometry relevant to planetary cores, comparing results from quasi-geostrophic (QG), 3-D and hybrid QG-3D models. The 170 reported calculations span Ekman numbers, Ek, between 10−4 and 10−10, Rayleigh numbers, Ra, between 2 and 150 times supercritical and Prandtl numbers, Pr, between 10 and 10−2. The default boundary conditions are no-slip at both the ICB and the CMB for the velocity field, with fixed temperatures at the ICB and the CMB. Cases driven by both homogeneous and inhomogeneous CMB heat flux patterns are also explored, the latter including lateral variations, as measured by Q*, the peak-to-peak amplitude of the pattern divided by its mean, taking values up to 5. The QG model is based on the open-source pizza code. We extend this in a hybrid approach to include the temperature field on a 3-D grid. In general, we find convection is dominated by zonal jets at mid-depths in the shell, with thermal Rossby waves prominent close to the outer boundary when the driving is weaker. For the thick spherical shell geometry studied here the hybrid method is best suited for studying convection at modest forcing, $Ra \le 10 \,more »Ra_c$ when Pr = 1, and departs from the 3-D model results at higher Ra, displaying systematically lower heat transport characterized by lower Nusselt and Reynolds numbers. We find that the lack of equatorially-antisymmetric motions and z-correlations between temperature and velocity in the buoyancy force contributes to the weaker flows in the hybrid formulation. On the other hand, the QG models yield broadly similar results to the 3-D models, for the specific aspect ratio and range of Rayleigh numbers explored here. We cannot point to major disagreements between these two data sets at Pr ≥ 0.1, with the QG model effectively more strongly driven than the hybrid case due to its cylindrically averaged thermal boundary conditions. When Pr is decreased, the range of agreement between the hybrid and 3-D models expands, for example up to $Ra \le 15 \, Ra_c$ at Pr = 0.1, indicating the hybrid method may be better suited to study convection in the low Pr regime. We thus observe a transition between two regimes: (i) at Pr ≥ 0.1 the QG and 3-D models agree in the studied range of Ra/Rac while the hybrid model fails when $Ra\gt 15\, Ra_c$ and (ii) at Pr = 0.01 the QG and 3-D models disagree for $Ra\gt 10\, Ra_c$ while the hybrid and 3-D models agree fairly well up to $Ra \sim 20\, Ra_c$. Models that include laterally varying heat flux at the outer boundary reproduce regional convection patterns that compare well with those found in similarly forced 3-D models. Previously proposed scaling laws for rapidly rotating convection are tested; our simulations are overall well described by a triple balance between Coriolis, inertia and Archimedean forces with the length-scale of the convection following the diffusion-free Rhines-scaling. The magnitude of Pr affects the number and the size of the jets with larger structures obtained at lower Pr. Higher velocities and lower heat transport are seen on decreasing Pr with the scaling behaviour of the convective velocity displaying a strong dependence on Pr. This study is an intermediate step towards a hybrid model of core convection also including 3-D magnetic effects.

    « less
  5. ABSTRACT FU Ori is the prototype of FU Orionis systems that are outbursting protoplanetary discs. Magnetic fields in FU Ori’s accretion discs have previously been detected using spectropolarimetry observations for Zeeman effects. We carry out global radiation ideal MHD simulations to study FU Ori’s inner accretion disc. We find that (1) when the disc is threaded by vertical magnetic fields, most accretion occurs in the magnetically dominated atmosphere at z ∼ R, similar to the ‘surface accretion’ mechanism in previous locally isothermal MHD simulations. (2) A moderate disc wind is launched in the vertical field simulations with a terminal speed of ∼300–500 km s−1 and a mass-loss rate of 1–10 per cent the disc accretion rate, which is consistent with observations. Disc wind fails to be launched in simulations with net toroidal magnetic fields. (3) The disc photosphere at the unit optical depth can be either in the wind launching region or the accreting surface region. Magnetic fields have drastically different directions and magnitudes between these two regions. Our fiducial model agrees with previous optical Zeeman observations regarding both the field directions and magnitudes. On the other hand, simulations indicate that future Zeeman observations at near-IR wavelengths or towards other FU Orionis systems maymore »reveal very different magnetic field structures. (4) Due to energy loss by the disc wind, the disc photosphere temperature is lower than that predicted by the thin disc theory, and the previously inferred disc accretion rate may be lower than the real accretion rate by a factor of ∼2–3.« less