2 D. Toyouchi et al. optically emittting magnetars and pulsars. In this paper, we quan- 107 titatively investigate the feasibility of hunting HVNSs in wide-area Low mass optical surveys and discuss how the future observations can be used High mass to constrain the formation and evolution of NSs. DM heating The rest of the paper is organized as follows. We first model the optical luminosity of NSs in Section 2.1 and describe the method 106 to quantify the detectability of HVNSs in Section 2.2. Based on the model calculations, we present the prediction of the number of [K] detectable HVNSs and clarify the properties of the obtained solution eff in Section 3.1, and also discuss the dependence of our results on T several model parameters in Section 3.2. Finally, the summary and 105 conclusion are given in Section4. 2 METHOD 104 103 104 105 106 107 108 109 2.1 Modeling Optical Luminosity of Neutron Stars age [yr] In this paper we consider the possibility of finding neutron stars (NS) with an optical observation. An NS is extremely faint due to its small Figure 1. Three NS cooling curves considered in this paper. The cyan, orange, size ('NS ∼ 10 km for the radius), but it is not impossible to find and green lines represent Low-mass, High-mass, and DM heating models, re- if an NS exists in the solar neighborhood, say within a few 100 pc spectively. For comparison, the observational data are overplotted (see Yanagi from the Earth. To make a quantitative estimate of detectability of et al. 2020, for references). NSs with an optical observation, we first need to model a luminosity of NS in optical wavelengths as a function of age. A radiation spectrum of NSs would generally consist of two contri- presented by Beniamini et al.(2019), which successfully reproduces butions: thermal radiation emitted from the neutron surface and non- the observed distributions of the spin period % and the derivative %¤ thermal radiation from magnetosphere. Each contribution changes for known magnetars. In this model, satisfying the general properties with time. We assume that the thermal emission from the surface of spin-down, Ω¤ / 퐵2Ω=, and the magnetic field decay, 퐵¤ / 퐵1¸U consists of a single black body spectrum of an effective temperature, (Colpi et al. 2000), the magnetic field and spin are assumed to evolve )NS, which declines with age of NSs, gNS, assuming the energy loss with time as via neutrino emission (see, e.g., Tsuruta et al. 2002; Page et al. 2004; UC −1/U Potekhin et al. 2020, for details). Since the cooling history of NSs 퐵 = 퐵0 1 ¸ , (1) is still uncertain, we here examine the uncertainty in our calculation g퐵 by considering three NS cooling curves as displayed in Figure1. and Here, “Low-” and “High-mass” models are taken from Figure 11 of 1 " ( 1 ) # 1−= Yakovlev & Pethick(2004), and dark matter (DM) heating model is ¹= − 1ºg퐵 UC 1−= Ω = Ω0 1 ¸ − 1 ¸ 1 , (2) taken from Yanagi et al.(2020). Low-mass model assumes an NS 2¹U − 2ºgΩ g퐵 mass of "NS = 1.35 " , in which the core density is too small to cause neutrino emission efficiently. This cooling curve is gen- where we set U = −1 and = = 3 as adopted in Beniamini et al. erally in good agreement with the observed temperature of young (2019), quantities with subscript 0 denote the quantities at the NS 6 birth, and g and g are the decay timescales of magnetic field and NSs with gNS . 10 yr. Note here that the absence of NSs with 퐵 Ω 5 spin, respectively. Assuming a dipole magnetic field, we describe the )NS . 5 × 10 K would reflect the fact that the observed thermal flux is quickly dimmed in X-ray band, so that optical observations spin-down timescale as, are advantageous in exploring the low-temperature regime. To argue Ω 322퐼 g = 0 = , (3) the detectability of such low-temperature NSs, we adopt High-mass Ω ¤ 2 6 2 2Ω0 4퐵 ' Ω model corresponding to greater NS mass case with "NS = 1.98 " , 0 NS 0 in which the thermal energy can be efficiently released via direct where 2 is the speed of light, and 퐼 is the inertia moment of NS. Urca process (Nadyozhin 1995). The spin-down evolution highly depends on the initial strength Note that both Low and High-mass models predict rapid decre- of magnetic field and its decay timescale. In our calculation, we 6 15 ments of )NS at ¡ 10 yr, in contrast to the thermal emission reported consider two NS models; a magnetar-like NS with 퐵0 = 10 G and 4 12 7 for very old millisecond pulsars (Kargaltsev et al. 2004; Webb et al. g퐵 = 10 yr, and a pulsar-like NS with 퐵0 = 10 G and g퐵 = 10 yr, −1 2019). Several possibilities for the internal heating mechanism as well and we set Ω0 = 100 s for both cases. The spin-down evolution as the heating due to non-thermal particles in the magnetosphere are of two NS models is shown in Figure2. In general, the spin and suggested (e.g. Gonzalez & Reisenegger 2010; Yanagi et al. 2020, magnetic field of NSs slowly declines when gNS g퐵, and then ¤ and references therein). Motivated by these, here we consider a dark suddenly decays at gNS ∼ g퐵. The magnetar-like NSs have higher % matter (DM) and rotochemical heating model, in which the surface at any % and more rapidly lose their rotational and magnetic energies temperature of NSs stays & 105 K for ∼ 1 Gyr (Hamaguchi et al. than the pulsar-like NSs do. ¤ ¤ 2019). Including such heating sources will enhance the detectability With the spin-down energy loss, 퐸sd = 퐼ΩΩ, we assess the NS op- of thermally emitting isolated NSs in optical wavelengths. tical luminosity adopting the conversion factor of the spin-down en- ¤ −2 −7 We also calculate the non-thermal radiation that arises from the de- ergy to the optical radiation, !opt/퐸sd = 10 and 10 for magnetars cay of the spin and the magnetic field of NSs using an analytic model and pulsars, respectively, based on the observational implication for MNRAS 000,1–8 (2021) Hunting neutron stars with wide-area optical surveys 3 0 6 1 15 4 10− Magnetar (Ω0 = 100 s− ,B0 = 10 G, τB = 10 yr) Thermal 1 12 7 Pulsar (Ω0 = 100 s− ,B0 = 10 G, τB = 10 yr) Magnetar 8 10− 1 yr 5 Pulsar 2 ] 10 10 yr 1 10− − 4 10 12 10 yr [s s 10− [mag] ˙ P 6 10 yr g 14 10− M 15 108 yr 16 10− 1010 yr 20 18 10− 2 1 0 1 2 10− 10− 10 10 10 25 P [s] 102 103 104 105 106 107 age [yr] Figure 2. Blue and orange lines show the spin evolution of magnetar- and pulsar-like NSs, respectively, calculated with the theoretical model by Be- Figure 3. Blue and orange lines represents the absolute magnitude in the niamini et al.(2019). Dotted lines represent the corresponding age of NSs 6-band of magnetar- and pulsar-like NSs, respectively, calculated from the calculated with %/%¤. composite spectra of the spin-down and the black body radiation. The black dashed line denotes the optical luminosity only from the thermal radiation isolated NSs (Mereghetti 2008; Mignani 2009; Enoto et al. 2019)2. calculated with Low-mass cooling model (see Figure1). Blue circles, orange With this model we consider the composite spectra of the spin-down squares, and green triangles represent the observed magnetars, pulsars, and X-ray dim NSs, respectively. The plotted magnitudes are dust extinction cor- and thermal radiation: Figure3 shows the absolute magnitude for rected although the uncorrected magnitude is also indicated for the magnetar magnetar- and pulsar-like NSs in optical wavelengths, compared to suffering from significant dust extinction. the observations of known NSs (Mignani 2009). While young mag- 4 netars with gNS < 10 yr are brighter than "6 ∼ 15 mag due to the greater spin-down radiation, the magnetars quickly become fainter detect with a high-quality ground-based data (Qiu et al. 2021) or after the magnetic field decays and then afterwards have a black body would be relatively easy to detect with a follow-up observation by an 6 radiation until gNS ∼ 10 yr. The spin-down radiation of pulsars is extremely large telescope such as the Thirty Meter Telescope (TMT) much weaker than the young magnetars, but the pulsars can be as or space-based observations such as the Roman Space Telescope. 5 bright as "6 ∼ 18 mag for ∼ 10 yr. The optical luminosity of mag- Hence, a task of finding an NS with an optical observation is to 6 netars and pulsars at gNS ¡ 10 yr is almost the same and arises from search for a very faint point-source object(s) that has a large proper the black body radiation, which fades with the NS cooling. Note that motion, together with a possible detection of the parallax. our model is generally in good agreement with the observed opti- In order to estimate the expected number of detectable NSs in cal luminosity of NSs except for young pulsars with ∼ 103 yr that the solar neighborhood, we consider the polar coordinate around the are thought to be extremely bright with the interaction between a Sun as illustrated in Figure4. We assume that each NS formed in central NS and the dense circum-stellar medium. Since the radiation the Galactic disk plane (i.e. ignore a possible offset of the newly mechanism powering such young pulsars is poorly understood, in this formed NS in the vertical direction from the disk plane) because NS paper we focus on the spin-down and thermal radiation and provide progenitors are short-lived enough to be regarded as a member of a conservative estimate in the number of detectable NSs.