Linking Multiwavelength Pulsar Emission to Super-Dense Physics inside neutron star

Neutron-star observations are now precise enough that the next major advances will depend not only on better data, but also on more physical forward models and more powerful inference tools. In this talk,  I will present a research program aimed at connecting pulsar magnetospheric physics,  X-ray pulse-profile modeling, and dense-matter inference within a unified framework. I will first highlight new results that derive,  from first principles and in analytic form, the surface return-current distribution in multipolar pulsar magnetospheres, providing a physically consistent route from global magnetic geometry to polar-cap heating and X-ray pulse profile. I will then briefly introduce two computational frameworks that make this broader program practical: a new public GPU-accelerated X-ray pulse-profile modeling framework, which reduces high-fidelity waveform evaluations from minutes to milliseconds, and CompactObject, an open-source full-scope platform for neutron-star equation-of-state inference that combines X-ray, gravitational-wave, radio, and nuclear constraints. Together, these developments provide a pathway from first-principles pulsar physics to robust neutron star dense matter constraints, with ongoing applications including physics-motivated hotspot inference for PSR J0740+6620 using NICER data and future NewAthena mission.