LAUSR

We theoretically study a stabilization mechanism of the skyrmion crystal in centrosymmetric magnets with a bilayer structure. We show that the interplay between a layer-dependent staggered Dzyaloshinskii-Moriya interaction that arises from the absence of local inversion symmetry and the interlayer exchange interaction gives rise to a plethora of multiple-Q states including the skyrmion crystal with a quantized topological number. By performing the simulated annealing for the bilayer triangular-lattice model under an external magnetic field, we demonstrate that the skyrmion forms the triangular-shaped crystals with different helicities in each layer owing to the staggered Dzyaloshinskii-Moriya interaction. Although the relative positions of the skyrmion core in each layer are different depending on the sign of the interlayer exchange interactions, the skyrmion crystal phases robustly appear under both ferromagnetic and antiferromagnetic interlayer interactions. We also find another two triple-Q states with a uniform scalar chirality but without a quantized topological number in the lowand high-field regions. Especially, the low-field triple-Q state exhibits the opposite sign of the scalar chirality to the skyrmion crystal, which are not found in the single-layer system. Our results indicate that the lack of local inversion symmetry in the lattice structure is another key ingredient to induce topological spin textures in centrosymmetric magnets.