This thesis presents the search of the Higgs boson produced by Vector Boson Fusion (VBF) and decaying to bottom quarks. A search using the ATLAS detector was performed with 2016 proton-proton collision data. The multi-variate analysis measured the signal strengths of both the inclusive Higgs production and the vector-boson fusion production relative to the Standard Model prediction. This analysis led to the observation of Higgs coupling to b-quarks in the summer of 2018. Furthermore, potential improvements of the analysis techniques using complex neural networks are investigated. In order to understand better the Quantum Chromodynamics (QCD) backgrounds of the Higgs search, the characteristic variables of the gluon splitting vertex are measured.