Models of magnetically-driven accretion and outflows reproduce many observational properties of T Tauri stars. This concept is not well established for the more massive Herbig Ae/Be stars. We intend to examine the magnetospheric accretion in Herbig Ae/Be stars and search for rotational modulation using spectroscopic signatures, in this first paper concentrating on the well-studied Herbig Ae star HD101412. We used near-infrared spectroscopic observations of the magnetic Herbig Ae star HD101412 to test the magnetospheric character of its accretion disk/star interaction. We reduced and analyzed 30 spectra of HD101412, acquired with the CRIRES and X-shooter spectrographs installed at the VLT (ESO, Chile). The spectroscopic analysis was based on the He I lambda 10,830 and Pa gamma lines, formed in the accretion region. We found that the temporal behavior of these diagnostic lines in the near-infrared spectra of HD101412 can be explained by rotational modulation of line profiles generated by accreting gas with a period P = 20.53+-1.68 d. The discovery of this period, about half of the magnetic rotation period P_m = 42.076 d previously determined from measurements of the mean longitudinal magnetic field, indicates that the accreted matter falls onto the star in regions close to the magnetic poles intersecting the line-of-sight two times during the rotation cycle. We intend to apply this method to a larger sample of Herbig Ae/Be stars.

We present the results of a study of the temporal behaviour of several diagnostic lines formed in the region of the accretion-disk/star interaction in the three magnetic Herbig Ae stars HD101412, HD104237, and HD190073. More than 100 spectra acquired with the ISAAC, X-shooter, and CRIRES spectrographs installed at the VLT-8m telescope (ESO, Chile), as well as at other observatories (OHP, Crimean AO) were analyzed. The spectroscopic data were obtained in the He I lambda10830, Pa gamma and He I lambda5876 lines. We found that the temporal behaviour of the diagnostic lines in the spectra of all program stars can be widely explained by a rotational modulation of the line profiles generated by a local accretion flow. This result is in good agreement with the predictions of the magnetospheric accretion model. For the first time, the rotation period of HD104237 (P_rot = 5.37+-0.03 days), as well as the inclination angle (i = 21+-4deg) were determined. Additional analysis of the HARPSpol spectra of HD104237 and HD190073, taken from the ESO archive, with the use of the SVD method shows that the magnetic field structure of HD190073 is likely more complex than a simple dipole and contains a circumstellar component. For the first time, the magnetic field of the secondary component of the binary system HD104237 was also detected (<B_z> = 128+-10G).