The objective of this thesis is to study the fracture behavior of bulk metallic glasses. For this purpose, detailed finite element investigation of the mode I and mixed mode (I and II) stationary crack tip fields under plane strain, small scale yielding conditions is carried out. An implicit backward Euler finite element implementation of the Anand and Su constitutive model [Anand, L. and Su, C., 2005, J. Mech. Phys. Solids 53, 1362] is used in the simulations. The effects of internal friction (μ), strain softening, Poisson's ratio (ν) and elastic mode mixity (Me) on the near-tip stress and deformation fields are examined. The results show that under mode I loading, a higher μ leads to a larger normalized plastic zone size and higher plastic strain level near the notch tip, but causes a substantial decrease in the opening stress. The brittle crack trajectories and shear band patterns around the notch are also simulated. An increase in ν reduces the extent of plastic zone and plastic strain levels in front of the notch tip. The results from mixed mode simulations show that increase in the mode II component of loading dramatically increases the maximum plastic zone extent, lowers the stresses and significantly enhances the plastic strain levels near the notch tip. Higher μ causes the peak magnitudes of tensile tangential stress to decrease. The implications of the above results on the fracture response of bulk metallic glasses are discussed. The possible variations of fracture toughness with mode mixity predicted by employing two simple fracture criteria are examined. Finally, mixed mode (I and II) fracture experiments on a Zr-based bulk metallic glass are performed. It is found that the fracture toughness increases with Me and Jc under mode I is higher than that under mode II loading by a factor of 4. The operative failure mechanism and fracture process zone size are discerned based on observations of incipient crack growth and fractographs. Lastly, a fracture criterion is proposed which predicts the experimentally observed variation of fracture toughness with mode mixity.