Many studies have shown that additional information can be gained on time series by investigating their associated complex networks. In this work, we investigate the multifractal property and Laplace spectrum of the horizontal visibility graphs (HVGs) constructed from fractional Brownian motions. We aim to identify via simulation and curve fitting the form of these properties in terms of the Hurst index H. First, we use the sandbox algorithm to study the multifractality of these HVGs. It is found that multifractality exists in these HVGs. We find that the average fractal dimension $\langle D(0)\rangle $ of HVGs approximately satisfies the prominent linear formula $\langle D(0)\rangle =2-H$ ; while the average information dimension $\langle D(1)\rangle $ and average correlation dimension $\langle D(2)\rangle $ are all approximately bi-linear functions of H when $H\geqslant 0.15$ . Then, we calculate the spectrum and energy for the general Laplacian operator and normalized Laplacian operator of these HVGs. We find that, for the general Laplacian operator, the average logarithm of second-smallest eigenvalue $\langle \ln \left({{u}_{2}}\right)\rangle $ , the average logarithm of third-smallest eigenvalue $\langle \ln \left({{u}_{3}}\right)\rangle $ , and the average logarithm of maximum eigenvalue $\langle \ln \left({{u}_{n}}\right)\rangle $ of these HVGs are approximately linear functions of H; while the average Laplacian energy $\langle {{E}_{\text{nL}}}\rangle $ is approximately a quadratic polynomial function of H. For the normalized Laplacian operator, $\langle \ln \left({{u}_{2}}\right)\rangle $ and $\langle \ln \left({{u}_{3}}\right)\rangle $ of these HVGs approximately satisfy linear functions of H; while $\langle \ln \left({{u}_{n}}\right)\rangle $ and $\langle {{E}_{\text{nL}}}\rangle $ are approximately a 4th and cubic polynomial function of H respectively.