A bstract The doubly Cabibbo-suppressed decay Ξ c + → pϕ with ϕ → K + K − is observed for the first time, with a statistical significance of more than fifteen standard deviations. The data sample used in this analysis corresponds to an integrated luminosity of 2 fb −1 recorded with the LHCb detector in pp collisions at a centre-of-mass energy of 8 TeV. The ratio of branching fractions between the decay Ξ c + → pϕ and the singly Cabibbo-suppressed decay Ξ c + → pK − π + is measured to be $$ \frac{\mathrm{\mathcal{B}}\left({\varXi}_c^{+}\to p\phi \right)}{\mathrm{\mathcal{B}}\left({\varXi}_c^{+}\to p{K}^{-}{\pi}^{+}\right)}=\left(19.8\pm 0.7\pm 0.9\pm 0.2\right)\times {10}^{-3}, $$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mfrac> <mml:mrow> <mml:mi>ℬ</mml:mi> <mml:mfenced> <mml:mrow> <mml:msubsup> <mml:mi>Ξ</mml:mi> <mml:mi>c</mml:mi> <mml:mo>+</mml:mo> </mml:msubsup> <mml:mo>→</mml:mo> <mml:mi>p</mml:mi> <mml:mi>ϕ</mml:mi> </mml:mrow> </mml:mfenced> </mml:mrow> <mml:mrow> <mml:mi>ℬ</mml:mi> <mml:mfenced> <mml:mrow> <mml:msubsup> <mml:mi>Ξ</mml:mi> <mml:mi>c</mml:mi> <mml:mo>+</mml:mo> </mml:msubsup> <mml:mo>→</mml:mo> <mml:mi>p</mml:mi> <mml:msup> <mml:mi>K</mml:mi> <mml:mo>−</mml:mo> </mml:msup> <mml:msup> <mml:mi>π</mml:mi> <mml:mo>+</mml:mo> </mml:msup> </mml:mrow> </mml:mfenced> </mml:mrow> </mml:mfrac> <mml:mo>=</mml:mo> <mml:mfenced> <mml:mrow> <mml:mn>19.8</mml:mn> <mml:mo>±</mml:mo> <mml:mn>0.7</mml:mn> <mml:mo>±</mml:mo> <mml:mn>0.9</mml:mn> <mml:mo>±</mml:mo> <mml:mn>0.2</mml:mn> </mml:mrow> </mml:mfenced> <mml:mo>×</mml:mo> <mml:msup> <mml:mn>10</mml:mn> <mml:mrow> <mml:mo>−</mml:mo> <mml:mn>3</mml:mn> </mml:mrow> </mml:msup> <mml:mtext>,</mml:mtext> </mml:math> where the first uncertainty is statistical, the second systematic and the third due to the knowledge of the ϕ → K + K − branching fraction.