The large constituent quark mass of bottomonium, the bottom quark/anti-quark bound state $(bbbar)$, affords a rich spectroscopy in which the perturbative (non-relativistic) limit of Quantum Chromodynamics may be theoretically described and experimentally investigated. The radial excitations of bottomonia---with radial quantum number $n$, one unit of total angular momentum $(J=1)$, and orbital angular momentum $L=0$, labeled $Upsilon(nS)$---are copiously produced in electron--positron $(epem)$ collisions. The Belle Collaboration is a high energy physics experiment located at the KEKB B-Factory epem collider, based at KEK in Tsukuba, Japan. Belle has accumulated a large dataset near the FourS and ThreeS resonances, collectively containing more than 28 million ThreeS and 556 million FourS. Some of these decay to other bbbar states---with one unit of orbital angular momentum and total angular momentum $J=0,1,2$, labeled cbj{n} ---via the emission of a photon, with subsequent transition to the OneS with the emission of one or more gluons, which hadronize to form an om meson. This dissertation presents an analysis of the hadronic transitions $chi_{bJ}(nP) rightarrow omega Upsilon(1S)$, where $Upsilon(1S) rightarrow ell^{+}ell^{-}$ with $ell=e,mu$, at Belle. The transitions of the $n=2$ triplet states provide a unique laboratory in which to study nonrelativistic quantum chromodynamics (NRQCD), as the kinematic threshold for production of an $omega$ and $Upsilon(1S)$ lies between the $J=0$ and $J=1$ states. The results presented herein constitute the first confirmation measurement of the $omega$ transitions of the $chi_{bJ}(2P)$ states since their discovery in 2004, with evidence---in excess of three standard deviations---for the sub-threshold transition of the $J=0$ state. The branching fraction $mathcal{B}big( chi_{b0}(2P) rightarrow omega Upsilon(1S) big)$ is found to be as large as the corresponding rate for the $J=2$ transition. The ratio of the $J=2$ to $J=1$ transitions is also measured and compared with the expectation from NRQCD, which we compute, revealing a $3.3sigma$ tension between experiment and theory. This work is leveraged to perform a search for radiative transitions of the $Upsilon(4S)$ to the $chi_{bJ}(2P)$ and $chi_{bJ}(3P)$ states, which are reconstructed in an inclusive $omega Upsilon(1S)$ final state. With no significant signal seen, limits are set on the corresponding branching fractions. / Doctor of Philosophy / Atoms, the stuff of everyday matter, consist of a number of electrons bound to a compact nucleus. This nucleus, in turn, contains one or more protons and neutrons, which are themselves made up of constituent particles called quarks that interact with one another by exchanging particles called gluons. Although great strides were made during the last century to further our understanding of the fundamental structure of matter, a comprehensive description of nuclear structure, at the quark level, eludes us. What we do know is that the force responsible for binding the large number of positively charged protons within the narrowly confined nucleus of, say, a gold atom is incredibly strong---in reality, more than 137 times as strong as the electromagnetic (EM) interaction, which is responsible for binding electrons around the nucleus in atoms. Unlike the EM force, which has one charge that can be either positive or negative, the strong interaction has three. This leads to a manifestly more complicated phenomena whose mathematical descriptions are computationally intractable. To study the strong interaction, we seek out the simplest of strongly bound states---called the meson---which consist of a quark and its anti-particle counterpart. The meson made up of a bottom quark/anti-quark pair, called bottomonium, provides an ideal laboratory for our investigations. In bottomonium, the quarks are very heavy (about 4.5 times the mass of a proton) and move relatively slowly compared to the quarks within a proton. This allows for some simplifications in the mathematical description of the bottomonium system, making it possible to compute predictions that can be tested in the lab. In this low energy regime, the strong interaction gives rise to a family of excited bottomonium states that have a structure similar to the excited states of an atom. Just as scientists learned about the EM interaction by studying the decays of excited atomic states, so too do we study the strong force by measuring the decays of bottomonium states. We call this study heavy quarkonium spectroscopy. When excited bottomonium states transition to lower-energy states, they may emit photons (as excited atoms do) or gluons. These emitted gluons, in turn, produce other particles. Measurements of the decay rates of bottomonium states may be predicted from the mathematical description of the strong interaction, providing direct experimental tests of the theoretical models. This dissertation presents a study of the decays of several bottomonium states, which are produced at the Belle experiment at the KEKB electron--positron collider. The decay rates, called the branching fractions, of these transitions are measured and used to test the prediction from theory, which we calculate. This work is leveraged to search for several previously unobserved decays, which are expected to be exceptionally rare.

In der Z-Boson-Resonanzphysik sind mehrere Präzisionsobservablen in einem perfekten Zustand, bei dem die theoretische Unsicherheit niedriger ist als die gegenwärtige experimentelle Unsicherheit. Das Konzept für den zukünftigen Teilchenbeschleuniger Future Circular Collider (FCC), will eine Verbesserung der Messungen für die Präzisionsobservablen um ein bis zwei signifikante Stellen erreichen. Damit werden die Vorhersagen des elektroschwachen Standardmodells in eine Situation versetzt, in der vollständige Zweischleifenkorrekturen zusammen mit den führenden Dreischleifenkorrekturen obligatorisch werden. 2016 wurden die vollständigen Zweischleifenkorrekturen für den effektiven schwachen Mischungswinkel für die bottom Quarks sin^2/theta/^b_eff berechnet, indem die fehlenden bosonischen Zweischleifenkorrekturen bereitgestellt wurden. Dabei liegt die Schwierigkeit in der Berechnung der entsprechenden Zwei-Schleifen Vertex-Feynman-Integrale, die mehrere massive Teilchen einschließen. Gegenwärtig ist die analytische Rechnung der meisten dieser Integrale schwierig und deswegen werden rein numerische Techniken, mittels Sektorzerlegungsansatz und der Integralansatz nach Mellin-Barnes, angewandt. Es war bis vor kurzem nicht bekannt, wie Mellin-Barnes-Integraldarstellungen in den minkowskischen Integrationsgebieten numerisch behandelt werden können. Um dieses Problem anzugehen, stellen wir eine Vielzahl von ein- und mehrdimensionaler Techniken vor, die ein Teil des neuen Programms MBnumerics.m sind, welches in dieser Dissertation entwickelt wurde. Der Sektorzerlegungsansatz und der Integralansatz nach Mellin-Barnes sind zusammen ausreichend, um elektroschwache Zweischleifenkorrekturen für die Präzisionsobservablen der Annihilation von e^+e^- in zwei Fermionen in der Z-Bosonresonanz auszurechnen. Aktuell führt dies zu der genauesten Vorhersage für den effektiven elektroschwachen Mischungswinkel für bottom Quarks sin^2/theta/^b_eff = 0.232312. / In the Z-boson resonance physics several precision observables are in a perfect state, where the theory uncertainty is lower than the present experimental uncertainty. The ambitious concepts for the future collider, Future Circular Collider (FCC), aim for an improvement of measurements for the precision observables by one to two significant digits. This will put the Electroweak Standard Model predictions in a situation where complete two-loop corrections together with the leading three-loop corrections will become mandatory. The complete two-loop corrections for effective weak mixing angle for bottom quarks sin^2/theta/^b_eff were reported recently, by providing the missing bosonic two-loop corrections. The difficult task in this computation is the calculation of the corresponding two-loop vertex Feynman integrals which include several massive particles. At present the analytic evaluation for most of these integrals is out of reach and purely numerical techniques were applied. Only two methods, sector decomposition approach and the Mellin-Barnes integral approach, are known to extract infrared and ultraviolet singularities in a systematic way for a general Feynman integral with fully automatized algorithms. It was not known until recently how to treat Mellin-Barnes integral representations in Minkowskian regions numerically. To address this problem we introduce and discuss in detail a variety of one- and multi-dimensional techniques, which are part of a new program MBnumerics.m developed in this thesis work. Two techniques, sector decomposition and Mellin-Barnes integral approach, are together sufficient to treat electroweak two-loop corrections to the precision observables for the e^+e^- annihilation into two fermions at the Z-boson resonance. This leads to the most precise prediction at present for the effective weak mixing angle for bottom quarks: sin^2/theta/^b_eff=0.232312.