Solving the Helmholtz equation provides wavefield solutions that are dimensionally compressed, per frequency, compared to the time domain which is useful for many applications, like full waveform inversion (FWI). However, the efficiency in attaining such wavefield solutions depends often on the size of the model, which tends to be large at high frequencies and for 3D problems. Thus, we use here a recently introduced framework based on neural networks to predict such solutions through setting the underlying physical equation as a loss function to optimize the neural network parameters. For an input point in the model space, the network learns to predict the wavefield value at that point, and its partial derivatives using a concept referred to as automatic differentiation, to fit, in our case, a form of Helmholtz equation. We specifically seek the solution of the scattered wavefield considering a simple homogeneous background model that allows for analytical solutions. Feeding a reasonable number of random points from the model space will ultimately train a fully connected 8-layer deep neural network with each layer having a dimension of 20, to predict the scattered wavefield function. Initial tests on a two-box-shaped scatterer model with a source in the middle, as well as, a layered model with a source on the surface demonstrate the successful training of the NN for this application and provide us with a peek into the potential of such an approach.