Five conical harmonic beams are generated from the interaction of femtosecond mid-infrared (mid-IR) pulses at a nominal input wavelength of 1997 nm with a 2D LiNbO 3 nonlinear photonic crystal with Sierpinski fractal superlattices. The main diffraction orders and the corresponding reciprocal vectors involved in the interaction are ascertained. Second and third harmonics emerging at external angles of 23.82°and 36.75°result from nonlinearČerenkov and Bragg diffractions, respectively. Three pathways of fourth-harmonic generation are observed at external angles of 14.21°, 36.5°, and 53.48°, with the first one resulting from nonlinearČerenkov diffraction, and the other two harmonics are generated via different cascaded processes.OCIS codes: 190.4420, 190.2620, 050.1940. doi: 10.3788/COL201715.051901.When a monochromatic wave passes through a homogeneous medium with its second-order nonlinear susceptibility varying in some regions, e.g., χ ð2Þ in a nonlinear photonic crystal (NPC), nonlinear diffraction can be detected for the harmonic rings or spots of the input fundamental wave [1][2][3][4][5][6] . Nonlinear diffraction, which is based on phase matching, includes nonlinearČerenkov diffraction (NCD), nonlinear Raman-Nath diffraction (NRND), and nonlinear Bragg diffraction (NBD). For NCD, the longitudinal phase matching (LPM) is satisfied when the longitudinal direction is parallel to the incident fundamental laser beam. This phase matching can be achieved by the addition of integer multiples of the incident wave vector or the combination of the incident wave vector and reciprocal vectors. In the case of NRND, the transverse (perpendicular to the incident fundamental laser beam) phase matching (TPM) has to be satisfied. This can be achieved by the combination of the reciprocal vectors. However, for the NBD process to occur, both LPM and TPM need to be satisfied. In the past twenty years, secondharmonic generation (SHG) via various nonlinear diffractions in 1D and 2D NPCs had been the major focus [7][8][9][10][11][12][13] . Although high-order harmonics (HH) are very important to the field of nonlinear optics, until recently, publications on HH generation using NPCs have been almost nonexistent. Recently, there have been several significant efforts to generate third harmonics in periodic, short-range ordered or radially poled NPCs, NPC waveguide and random quadratic media by NCD and NRND [14][15][16][17][18] . Also, high efficiency quasi-phase-matched harmonic generation from the 2nd to 8th order have been observed within a single LiNbO 3 (LN) NPC [19] . However, fourth-harmonic generation (FHG) via nonlinear diffraction is rarely reported. The advantage of using a 2D fractal superlattice is that it may allow the freedom of optimization of many parameters such as input wavelength, harmonic order, angle, and efficiency. In this Letter, second, third, and especially multiple fourth harmonics are achieved by nonlinear diffractions in an LN NPC with Sierpinski fractal superlattices under femtosecond pulses.In our experim...