We present a study of ions produced by the perpendicular neutral beam in the high-beta plasma of Large Helical Device (LHD) that uses the Monte-Carlo code based on orbit following in real coordinates with Coulomb collisions. For no neutrals outside the plasma, the results indicate that there are remarkable differences in the distribution functions that depend on whether the loss boundary is set on the vacuum vessel wall or on the last closed flux surface (LCFS). Also, for finite neutral density n H = 10 18 m −3 , we find that the minor radial profile of the distribution function slightly differs from that for the infinite neutral density case, which corresponds to the loss boundary being set at the LCFS. Keywords: re-entering particle, LHD, perpendicular neutral beam, charge exchange loss, Monte-Carlo simulation DOI: 10.1585/pfr.5.014A perpendicular neutral beam (NB) injector and three tangential-NB injectors haves been installed on Large Helical Device (LHD). Using these NBs, the volume-averaged beta β has reached 5% in recent LHD experiments [1]. In the high-beta plasma of LHD, it has been pointed out that the beam pressure and/or the pressure anisotropy significantly affect the properties of magneto-hydrodynamics equilibrium and stability [2]. Therefore, it is important to identify the contribution of the beam-pressure to the total plasma pressure. Since the beam pressure is calculated from the distribution function of high-energy ions produced by the NBs, it is also important to investigate the distribution function.To analyze the distribution function of the highenergy ions in LHD, we developed a Monte-Carlo code (MORH) [3] based on an orbit following in real coordinates with Coulomb collisions [4]. In this code, we set the particle-loss boundary on the vacuum vessel (VV) wall instead of the last closed flux sufrace (LCFS), because we can follow the ion orbit until the VV wall. Thus, reentering particles [5,6], which repeatedly pass into and out of the LCFS, can be treated. Additionally, the effect of charge exchange (CX) loss [7] on the re-entering particles can be investigated. This code has already been applied to the ions produced by the tangential-NBs in LHD, and it was shown that the re-entering particles play an important role in the analyzing the distribution function of ions . On the other hand, it has also been shown that most of the protons produced by the perpendicular-NB become the re-entering particles in high-beta plasmas [8]. Thus, the re-entering particles may play a more important role in the distribution function in the high-beta plasma.In this paper, we apply the MORH code to the protons produced by the perpendicular-NB in the high-beta plasma of LHD in the case of B ax = 3 T for comparison with the results in Ref. [3]. For the calculation, we assume that the NB (40 keV, 1 MW) is injected torus-outside at the horizontally elongated poloidal cross-section in the equatorial plane. The number of the high-energy protons produced by the NB is assumed to be uniform on the beam line, and their ...