We study the effect of trigonal warping on the focussing of electrons by n-p junctions in graphene. We find that perfect focussing, which was predicted for massless Dirac fermions, is only preserved for one specific sample orientation. In the general case, trigonal warping leads to the formation of cusp caustics, with a different position of the focus for graphene's two valleys. We develop a semiclassical theory to compute these positions and find very good agreement with tight-binding simulations. Considering the transmission as a function of potential strength, we find that trigonal warping splits the single Dirac peak into two distinct peaks, leading to valley polarization. We obtain the transmission curves from tight-binding simulations and find that they are in very good agreement with the results of a billiard model that incorporates trigonal warping. Furthermore, the positions of the transmission maxima and the scaling of the peak width are accurately predicted by our semiclassical theory. Our semiclassical analysis can easily be carried over to other Dirac materials, which generally have different Fermi surface distortions.Veselago lenses [1] are special types of lenses, which are made of materials with a negative refractive index. Such lenses can overcome the diffraction limit [2] and can nowadays be realized in metamaterials [3][4][5], chiral metamaterials [6-9] and photonic crystals [10,11]. An electronic analog of a Veselago lens can be created using n-p junctions, with the classical trajectories of the charge carriers playing the role of the rays in geometrical optics. Such junctions focus electrons, because the group velocity for holes is in the direction opposite to their phase velocity, whereas the two velocities are in the same direction for electrons. However, in conventional semiconductors, such interfaces are unsuitable because of their high reflectivity, owing to the presence of a depletion region.Cheianov et al. [12] realized that graphene does not have this drawback. It has zero bandgap, as the valence and conduction bands touch at two non-equivalent corners of the Brillouin zone, known as K and K . The lowenergy charge carriers are ballistic over large distances and follow the Dirac equation [13][14][15][16][17][18]. This gives rise to Klein tunneling: normally incident electrons are transmitted with unit probability [19][20][21][22][23][24][25][26][27], which makes the interface exceptionally transparent. Recently, Veselago lensing in graphene was experimentally observed [27,28].Theoretical studies have used the Dirac equation to investigate focussing by flat [12,29] and circular [30][31][32] junctions and in zigzag nanoribbons [33]. It was shown that when the electron and hole charge carrier densities are equal, a flat interface is able to focus all trajectories into a single point [12]. According to catastrophe theory [34][35][36][37], such a situation is exceptional: any perturbation of the Hamiltonian will ruin this ideal focus. Indeed, in tight-binding simulations of an n-p junction...