We simulate the three-dimensional quantum Heisenberg model with a spatially anisotropic ladder pattern using the first principles Monte Carlo method. Our motivation is to investigate quantitatively the newly established universal relation TN / √ c 3 ∝ Ms near the quantum critical point (QCP) associated with dimerization. Here TN , c, and Ms are the Néel temperature, the spinwave velocity, and the staggered magnetization density, respectively. For all the physical quantities considered here, such as TN and Ms, our Monte Carlo results agree nicely with the corresponding results determined by the series expansion method. In addition, we find it is likely that the effect of a logarithmic correction, which should be present in (3+1)-dimensions, to the relation TN / √ c 3 ∝ Ms near the investigated QCP only sets in significantly in the region with strong spatial anisotropy.Introduction.-While being the simplest models, Heisenberg-type models provide qualitatively, or even quantitatively useful information regarding the properties of cuprate materials. For example, the spatially anisotropic quantum Heisenberg model with different antiferromagnetic couplings in the 1 and 2 directions is demonstrated to be relevant for the underdoped cuprate superconductor YBa 2 Cu 3 O 6.45 [1,2]. Specifically, it is argued that this model provides a possible mechanism for the newly discovered pinning effects of the electronic liquid crystal in YBa 2 Cu 3 O 6.45 [3]. Because of their phenomenological importance, these models continue to attract a lot of attention analytically and numerically. In addition to being relevant to real materials, Heisenbergtype models on geometrically nonfrustrated lattices are important from a theoretical point of view as well. This is because these models can be simulated very efficiently using first principles Monte Carlo methods. Hence they are very useful in exploring ideas and examining theoretical predictions [4][5][6][7][8][9][10][11][12].