The molecular dipole moment and magnetic hyperfine structure constant demand an accurate wavefunction far from the nucleus and in near nuclear region, respectively. We, therefore, employ the so-called Z-vector method in the domain of relativistic coupled cluster theory to calculate the first order property of molecular systems in their open-shell ground state configuration. The implemented method is applied to calculate molecular dipole moment and parallel component of the magnetic hyperfine structure constant of SrF molecule. The results of our calculation are compared with the experimental and other available theoretically calculated values. We are successful in achieving good accordance with the experimental results. The result of our calculation of molecular dipole moment is in the accuracy of ∼ 0.5 %, which is clearly an improvement over the previous calculation based on the expectation value method in the four component coupled cluster framework [V. S. Prasannaa et al, Phys. Rev. A 90, 052507 (2014)] and it is the best calculated value till date. Thus, it can be inferred that the Z-vector method can provide an accurate wavefunction in both near and far nuclear region, which is evident from our calculated results.PACS numbers: 31.15.bw, 31.15.vn, 32.10.Fn Theoretical physicists find it very challenging to calculate the spectroscopic properties of atoms and molecules. The precise description of the spectroscopic properties demands the wavefunction to be accurate both in the nuclear region and the region far from the nucleus. The calculation of an accurate wavefunction involving heavy atoms and molecules needs to include the relativistic and electron correlation effects simultaneously, as these two effects are non-additive in nature [1,2]. The best possible way to include the effects of relativity in a single determinantal theory is to solve the Dirac-Hartree-Fock (DHF) Hamiltonian in its four component formalism. The DHF Hamiltonian converts the complicated many electron problem into a sum of many one-electron problems by assuming an average electron-electron interaction. Therefore, the DHF Hamiltonian lacks the correlation of opposite spin electrons. The missing electron correlation can be included by adding orthogonal space to the DHF wavefunction. On the other hand, normal coupled cluster (NCC) [3][4][5] method is known to be the most elegant many-body theory to effectuate the dynamic part of the electron correlation.The calculations of one electron response properties in the NCC framework, can either be done by taking expectation value of the desired property operator or as a derivative of energy. These two approaches are not same as the NCC is by nature non-variational. In fact, the first order derivative of energy is the corresponding expectation value plus some additional terms, which makes the derivative approach closer to the full configuration * sk.sasmal@ncl.res.in interaction (FCI) property value. It is worth to mention that the expectation value approach in the NCC leads to a nonterminating ...