The discovery of the Higgs boson at the Large Hadron Collider (LHC) in 2012 marked a significant milestone in particle physics. Since then, extensive investigations have been carried out to unravel the properties of the Higgs boson and understand its interactions with other fundamental particles. Monte Carlo simulations play an important role in studies of particle physics. They generate theoretical predictions and can be compared with experimental data to test our understanding of physics. In this dissertation, the anomalous Higgs boson coupling to vector bosons was implemented in VBFNLO for Higgs boson plus three jets via VBF at the next-to-leading order (NLO) accuracy. I present results for a number of collider observables that are sensitive to the CP structure, whether in different CP scenarios and in the Standard Model. I also investigate the effects of NLO multijet merging and NLO QCD corrections matched to a parton shower on selected observables through the framework provided by Herwig 7.