Weak gravitational lensing (WL) peak statistics capture cosmic non-linear structures and can provide additional cosmological information complementary to cosmic shear two-point correlation analyses. They have been applied to different WL surveys successfully. To further facilitate their high precision applications, it is very timely to investigate the impacts of different systematics on WL peak statistics and how to mitigate them. Concerning the influence from galaxy intrinsic alignments (IAs), in this paper, we develop a theoretical model for WL high peaks taking into account the IA effects. It is an extension of our previous halo-based model. The IA corrections mainly include the modification of the lensing profile of clusters of galaxies due to the alignments of satellite galaxies and the additional shape noise correlations. We validate our model using simulations with the semi-analytical galaxy formation. We consider the cases where the satellite galaxies are averagely radially aligned toward the centers of their host clusters but with different dispersions σ_θ. We show that our model works well for σ_θ>45^. If the IA corrections are not included in the model, for the Euclid/CSST-like source galaxy distribution and the survey area of 1000 deg², the IA induced bias on S₈ can reach 8σ even for σ_θ=75^. With our model, not only the bias can be well mitigated, but also we can constrain the satellite IA to the level of σ (σ_θ) 24^{} simultaneously from WL high peak analyses alone using data from such a survey.
Zhang et al. (Sat,) studied this question.