Strongly correlated low-dimensional quantum spin models provide a well-established frame- work to study magnetic properties of insulators, and are of great theoretical interest and experimental relevance in condensed-matter physics. In this thesis, I use quantum Monte Carlo methods to numerically study quantum critical behavior in low-dimensional quantum spin models and wavefunctions. First, I study spinons &ndash emergent spin-1/2 bosonic excitations &ndash at certain one- and two-dimensional quantum phase transitions (QPTs) in spin models, by characterizing their size and confinement length quantitatively. In particular, I focus on the QPT from an antiferromagnetic (AFM) phase into a valence-bond solid (VBS) phase, which is an example of a violation of the standard Landau-Ginzburg-Wilson paradigm for phase transitions. This transition in two dimensions (2D) is instead likely described by a novel theory called "deconfined quantum criticality" (DQC). According to the theory, spinons should be deconfined. The degree of deconfinement is quantified in my calculations. Second, I present a comprehensive study of so-called short-bond resonating-valence-bond (RVB) spin liquids in 2D, which have been suggested as a good starting point for understanding the spin physics of high-temperature cuprates. I find that these RVB states can also be classified as quantum-critical VBS states, which indicates that RVB is less disordered than expected. This work suggests a possible mapping from the quantum RVB states to classical dimer models via a classical continuum field theory--the height model. This map explicitly bridges well-established classical results to future quantum studies. Third, I consider 1D amplitude product (AP) states, which are generalized versions of RVB states, with different wavefunction weightings of bonds according to their lengths. AP states constitute a good ansatz for certain Hamiltonians and are of broad interest in quantum magnetism. I study phase transitions from AFM-VBS phases in AP states by tuning their amplitudes, and obtain continuously varying critical exponents. In addition, I classify the 1D AP states through entanglement entropy calculations of the central charge in (1+1)D conformal field theory. This new classification could serve as guide for AP states as trial wavefunctions to search for ground states of corresponding quantum spin models.

In the study of strongly interacting condensed-matter systems controlled microscopic theories hold a key position. Spin-wave theory, large-N expansion, and $epsilon$-expansion are some of the few successful cornerstones. In this doctoral thesis work, we have developed a novel large-$d$ expansion method, $d$ being the spatial dimension, to study model Hamiltonians hosting a quantum phase transition between a paramagnet and a magnetically ordered phase. A highlight of this technique is that it can consistently describe the entire phase diagram of the above mentioned models, including the quantum critical point. Note that most analytical techniques either efficiently describe only one of the phases or suffer from divergences near the critical point. The idea of large-$d$ formalism is that in this limit, non-local fluctuations become unimportant and that a suitable product state delivers exact expectation values for local observables, with corrections being suppressed in powers of $1/d$. It turns out that, due to momentum summation properties of the interaction structure factor, all diagrams are suppressed in powers of $1/d$ leading to an analytic expansion. We have demonstrated this method in two important systems namely, the coupled-dimer magnets and the transverse-field Ising model. Coupled-dimer magnets are Heisenberg spin systems with two spins, coupled by intra-dimer antiferromagnetic interaction, per crystallographic unit cell (dimer). In turn, spins from neighboring dimers interact via some inter-dimer interaction. A quantum paramagnet is realized for a dominant intra-dimer interaction, while a magnetically ordered phase exists for a dominant (or of the same order as intra-dimer interaction) inter-dimer interaction. These two phases are connected by a quantum phase transition, which is in the Heisenberg O(3) universality class. Microscopic analytical theories to study such systems have been restricted to either only one of the phases or involve uncontrolled approximations. Using a non-linear bond-operator theory for spins with S=$1/2$, we have calculated the $1/d$ expansion of static and dynamic observables for coupled dimers on a hypercubic lattice at zero temperature. Analyticity of the $1/d$ expansion, even at the critical point, is ensured by correctly identifying suitable observables using the mean-field critical exponents. This method yields gapless excitation modes in the continuous symmetry broken phase, as required by Goldstone\'s theorem. In appropriate limits, our results match with perturbation expansion in small ratio of inter-dimer and intra-dimer coupling, performed using continuous unitary transformations, as well as the spin-wave theory for spin-$1/2$ in arbitrary dimensions. We also discuss the Brueckner approach, which relies on small quasiparticle density, and derive the same $1/d$ expansion for the dispersion relation in the disordered phase. Another success of our work is in describing the amplitude (Higgs) mode in coupled-dimer magnets. Our novel method establishes the popular bond-operator theory as a controlled approach. In $d=2$, the results from our calculations are in qualitative agreement with the quantum Monte Carlo study of the square-lattice bilayer Heisenberg AF spin-$1/2$ model. In particular, our results are useful to identify the amplitude (Higgs) mode in the QMC data. The ideas of large-$d$ are also successfully applied to the transverse-field Ising model on a hypercubic lattice. Similar to bond operators, we have introduced auxiliary Bosonsic operators to set up our method in this case. We have also discussed briefly the bilayer Kitaev model, constructed by antiferromagnetically coupling two layers of the Kitaev model on a honeycomb lattice. In this case, we investigate the dimer quantum paramagnetic phase, realized in the strong inter-layer coupling limit. Using bond-operator theory, we calculate the mode dispersion in this phase, within the harmonic approximation. We also conjecture a zero-temperature phase diagram for this model.

In the study of strongly interacting condensed-matter systems controlled microscopic theories hold a key position. Spin-wave theory, large-N expansion, and $epsilon$-expansion are some of the few successful cornerstones. In this doctoral thesis work, we have developed a novel large-$d$ expansion method, $d$ being the spatial dimension, to study model Hamiltonians hosting a quantum phase transition between a paramagnet and a magnetically ordered phase. A highlight of this technique is that it can consistently describe the entire phase diagram of the above mentioned models, including the quantum critical point. Note that most analytical techniques either efficiently describe only one of the phases or suffer from divergences near the critical point. The idea of large-$d$ formalism is that in this limit, non-local fluctuations become unimportant and that a suitable product state delivers exact expectation values for local observables, with corrections being suppressed in powers of $1/d$. It turns out that, due to momentum summation properties of the interaction structure factor, all diagrams are suppressed in powers of $1/d$ leading to an analytic expansion. We have demonstrated this method in two important systems namely, the coupled-dimer magnets and the transverse-field Ising model. Coupled-dimer magnets are Heisenberg spin systems with two spins, coupled by intra-dimer antiferromagnetic interaction, per crystallographic unit cell (dimer). In turn, spins from neighboring dimers interact via some inter-dimer interaction. A quantum paramagnet is realized for a dominant intra-dimer interaction, while a magnetically ordered phase exists for a dominant (or of the same order as intra-dimer interaction) inter-dimer interaction. These two phases are connected by a quantum phase transition, which is in the Heisenberg O(3) universality class. Microscopic analytical theories to study such systems have been restricted to either only one of the phases or involve uncontrolled approximations. Using a non-linear bond-operator theory for spins with S=$1/2$, we have calculated the $1/d$ expansion of static and dynamic observables for coupled dimers on a hypercubic lattice at zero temperature. Analyticity of the $1/d$ expansion, even at the critical point, is ensured by correctly identifying suitable observables using the mean-field critical exponents. This method yields gapless excitation modes in the continuous symmetry broken phase, as required by Goldstone\'s theorem. In appropriate limits, our results match with perturbation expansion in small ratio of inter-dimer and intra-dimer coupling, performed using continuous unitary transformations, as well as the spin-wave theory for spin-$1/2$ in arbitrary dimensions. We also discuss the Brueckner approach, which relies on small quasiparticle density, and derive the same $1/d$ expansion for the dispersion relation in the disordered phase. Another success of our work is in describing the amplitude (Higgs) mode in coupled-dimer magnets. Our novel method establishes the popular bond-operator theory as a controlled approach. In $d=2$, the results from our calculations are in qualitative agreement with the quantum Monte Carlo study of the square-lattice bilayer Heisenberg AF spin-$1/2$ model. In particular, our results are useful to identify the amplitude (Higgs) mode in the QMC data. The ideas of large-$d$ are also successfully applied to the transverse-field Ising model on a hypercubic lattice. Similar to bond operators, we have introduced auxiliary Bosonsic operators to set up our method in this case. We have also discussed briefly the bilayer Kitaev model, constructed by antiferromagnetically coupling two layers of the Kitaev model on a honeycomb lattice. In this case, we investigate the dimer quantum paramagnetic phase, realized in the strong inter-layer coupling limit. Using bond-operator theory, we calculate the mode dispersion in this phase, within the harmonic approximation. We also conjecture a zero-temperature phase diagram for this model.