This study proposed fully nonlinear free surface physics-informed neural networks (FNFS-PINNs), an advancement within the framework of PINNs, to tackle wave propagation in fully nonlinear potential flows with the free surface. Utilizing the nonlinear fitting capabilities of neural networks, FNFS-PINNs offer an approach to addressing the complexities of modeling nonlinear free surface flows, broadening the scope for applying PINNs to various wave propagation scenarios. The improved quasi-σ coordinate transformation and dimensionless formulation of the basic equations are adopted to transform the time-dependent computational domain into the stationary one and align variable scale changes across different dimensions, respectively. These innovations, alongside a specialized network structure and a two-stage optimization process, enhance the mathematical formulation of nonlinear water waves and solvability of the model. FNFS-PINNs are evaluated through three scenarios: solitary wave propagation featuring nonlinearity, regular wave propagation under high dispersion, and an inverse problem of nonlinear free surface flow focusing on the back-calculation of an initial state from its later state. These tests demonstrate the capability of FNFS-PINNs to compute the propagation of solitary and regular waves in the vertical two-dimensional scenarios. While focusing on two-dimensional wave propagation, this study lays the groundwork for extending FNFS-PINNs to other free surface flows and highlights their potential in solving inverse problems.