We establish the existence of solutions to the following semilinear Neumann problem for fractional Laplacian and critical exponent:{(−Δ)su+λu=|u|p−1uinΩ,Nsu(x)=0inRn∖Ω‾,u≥0inΩ, where λ>0 is a constant and Ω⊂Rn is a bounded domain with smooth boundary. Here, p=n+2sn−2s is a critical exponent, n>max{4s,8s+23}, s∈(0,1). Due to the critical exponent in the problem, the corresponding functional Jλ does not satisfy the Palais-Smale (PS)-condition and therefore one cannot use standard variational methods to find the critical points of Jλ. We overcome such difficulties by establishing a bound for Rayleigh quotient and with the aid of nonlocal version of the Cherrier's optimal Sobolev inequality in bounded domains. We also show the uniqueness of these solutions in small domains.