We investigated a wealth of X-ray and gamma-ray spectral energy distribution (SED) and multiband light-curve (LC) data of the gamma-ray binary HESS J0632+057 using a phenomenological intrabinary shock (IBS) model. Our baseline model assumes that the IBS is formed by colliding winds from a putative pulsar and its Be companion and that particles accelerated in the IBS emit broadband radiation via synchrotron (SY) and inverse Compton upscattering (ICS) processes. Adopting the latest orbital solution and system geometry, we reproduced the global X-ray and TeV LC features, two broad bumps at ϕ ∼ 0.3 and ∼0.7, with the SY and ICS model components. We found that these TeV LC peaks originate from ICS emission caused by the enhanced seed photon density near periastron and superior conjunction or Doppler-beamed emission of bulk-accelerated particles in the IBS at inferior conjunction. While our IBS model successfully explained most of the observed SED and LC data, we found that phase-resolved SED data in the TeV band require an additional component associated with ICS emission from preshock particles (produced by the pulsar wind). This finding indicates a possibility of delineating the IBS emission components and determining the bulk Lorentz factors of the pulsar wind at certain orbital phases.

HESS J0632+057 belongs to a rare subclass of binary systems that emit gamma rays above 100 GeV. It stands out for its distinctive high-energy light curve, which features a sharp “primary” peak and broader “secondary” peak. We present the results of contemporaneous observations by NuSTAR and VERITAS during the secondary peak between 2019 December and 2020 February, when the orbital phase (ϕ) is between 0.55 and 0.75. NuSTAR detected X-ray spectral evolution, while VERITAS detected TeV emission. We fit a leptonic wind-collision model to the multiwavelength spectra data obtained over the four NuSTAR and VERITAS observations, constraining the pulsar spin-down luminosity and the magnetization parameter at the shock. Despite long-term monitoring of the source from 2019 October to 2020 March, the MDM observatory did not detect significant variation in Hα and Hβ line equivalent widths, an expected signature of Be-disk interaction with the pulsar. Furthermore, fitting folded Swift-XRT light-curve data with an intrabinary shock model constrained the orbital parameters, suggesting two orbital phases (at ϕ D = 0.13 and 0.37), where the pulsar crosses the Be-disk, as well as phases for the periastron (ϕ 0 = 0.30) and inferior conjunction (ϕ IFC = 0.75). The broadband X-ray spectra with Swift-XRT and NuSTAR allowed us to measure a higher neutral hydrogen column density at one of the predicted disk-passing phases.