Physics and Astronomy

Permanent URI for this collectionhttps://uwspace.uwaterloo.ca/handle/10012/9949

This is the collection for the University of Waterloo's Department of Physics and Astronomy.

Research outputs are organized by type (eg. Master Thesis, Article, Conference Paper).

Waterloo faculty, students, and staff can contact us or visit the UWSpace guide to learn more about depositing their research.

Browse

Now showing 1 - 20 of 818

The Dynamical States and Mass Accretion Histories of Galaxy Clusters in IllustrisTNG
(University of Waterloo, 2025-09-19) Reid, Rashaad
The concordance cosmological model describes the history and large-scale structure of the universe using a few key parameters. Two of these parameters, σ8 and Ωm, determine the clustering of matter due to the growth of density fluctuations in the early universe. Current constraints on these parameters measured from nearby large structures and from the early universe are in statistical tension. Our ability to resolve this tension is limited by the degeneracy between the parameters when measured from observations of nearby structure. Since galaxy clusters are the most massive gravitationally bound objects in the universe, their formation is sensitive to σ8 and Ωm. An improved understanding of the formation histories of galaxy clusters can break the measurement degeneracy, thus providing new insights into this tension in our cosmological model. Since the formation time scales of galaxy clusters are unobservable, we must use the structure of clusters to probe their formation histories. In this thesis, we relate the observable structural properties of galaxy clusters to the mass accretion histories of their surrounding dark matter halos in the IllustrisTNG cosmo- logical simulations. Structure formation in the universe is hierarchical, so recently formed galaxy clusters will have experienced recent mergers with other systems. We examine a set of structural properties that are related to the dynamical states of clusters as indications of recent mergers to relate the structures of clusters to their formation histories. Using the cluster formation history information that is available in IllustrisTNG, we classify clusters as dynamically relaxed or unrelaxed based on their structural properties and compare the mean mass accretion histories of the resulting groups. We establish in this work that the stellar mass asymmetry and magnitude gaps of galaxy clusters are readily observable structural parameters that most effectively predict the mass accretion histories of halos. By comparing the gravitational lensing profiles of dynamically relaxed and unrelaxed clusters classified using different structural parameters, we demonstrate that the stellar mass asymmetry most reliably distinguishes between halos in different dynamical states with different density profiles. We also show that line of sight galaxy projection does not significantly affect IllustrisTNG cluster samples and that differences between 3D-identified clusters and optically selected clusters can be accounted for with accurate cluster mass estimates. However, we find that the density profiles traced by the weak gravitational lensing around relaxed and unrelaxed clusters in IllustrisTNG simulations and DLIS x UNIONS observations are discrepant. The structural differences between the simulated and observed galaxy clusters will be further explored in future work to better relate this work’s findings to real astrophysical systems. Overall, we find through cosmological simulations that the structural properties of galaxy clusters can be used to effectively trace their mass accretion histories. The findings of this thesis establish which observable properties of clusters can be targeted in both observations and simulations to grant us insight into the formation histories of the largest structures in the universe.
Fluctuation-induced order and thermal transport in frustrated quantum magnets
(University of Waterloo, 2025-09-18) Hickey, Alexander
Order-by-disorder is a mechanism of "fluctuation-induced" ordering that occurs in many frustrated magnetic systems where magnetic moments, or spins, are subject to competing interactions. So far, this phenomenon has been discussed in systems where the quantum ground state is not a "classical" product state. In such a case, both thermal and quantum fluctuations act to lift the accidental classical degeneracy, raising the question of whether one mechanism of order-by-disorder is possible without the other. In this thesis, we present results exposing a novel route to order-by-disorder, one without quantum zero-point fluctuations, in the ferromagnetic pyrochlore Heisenberg model with the Dzyaloshinskii-Moriya (DM) interaction as the leading perturbation. We show that any collinear ferromagnetic state is an exact eigenstate even in the presence of the anisotropic DM interaction, while thermal fluctuations give rise to a preferred magnetization direction. Using linear spin wave theory, we find that the anisotropy appears at lowest order as a sub-leading term in the low-temperature expansion of the free energy. Our results thus show that the phenomenon of thermal order-by-disorder can, in principle, occur even in the absence of quantum zero-point fluctuations driving quantum order-by-disorder. By extending our calculations to non-linear spin wave theory, we find that the ferromagnetic ground state becomes unstable for a spin-1/2 system when the DM interaction is large. Next, we ask the question of how to adequately characterize order-by-disorder in real materials, and how to distinguish it from conventional energetic ordering. Currently, the only clear and universal signature that has been proposed is a characteristic temperature dependence of the fluctuation-induced pseudo-Goldstone gap. Thus far, this temperature dependence of the pseudo-Goldstone gap has only been characterized in the classical limit. Here, we use non-linear spin wave theory to characterize the pseudo-Goldstone gap in quantum magnets at low temperature, to leading order in 1/S. Using exact sum-rules for the magnon spectral functions, we find that the gap exhibits a distinct power-law temperature dependence. We examine the implications of our results for several candidate materials. The final part of this thesis examines the thermodynamic and transport properties of the ferromagnetic pyrochlore Lu₂V₂O₇. Over the last decade, there has been immense interest in magnetic materials that host topologically non-trivial excitations. In ordered magnetic insulators, features analogous to those of topological insulators and semimetals can arise in the magnon band structure, and the associated Berry phases can manifest in observable heat and spin transport phenomena. This was unambiguously observed in Lu₂V₂O₇ in the form of a magnon thermal Hall signal, and proposed to arise from the DM interaction. A precise value of the DM interaction is not known, as the values obtained from fitting both thermal transport and inelastic neutron scattering data, as well as from density functional theory, are all mutually inconsistent. Motivated by this, we investigate the effect of additional symmetry allowed perturbations to the spin Hamiltonian of Lu₂V₂O₇ in an attempt to reconcile the different experimental probes of this material. We find that the thermal transport and neutron scattering measurements are consistent with the addition of a small second-nearest-neighbour DM interaction to the model. Conversely, we argue that existing specific heat measurements are inconsistent with neutron scattering experiments and cannot be reconciled with any additional exchange couplings to the bilinear spin-model in the perturbative regime. Our results motivate future thermal transport and specific heat measurements of this material.
Asymptotic Higher Spin Symmetries: Noether Realization & Algebraic Structure in Einstein-Yang-Mills Theory
(University of Waterloo, 2025-09-17) Cresto, Nicolas
This thesis deals with the phase space realization of asymptotic higher spin symmetries, in 4-dimensional asymptotically flat spacetimes. These symmetries live on the conformal null boundary, namely null infinity, and were first revealed few years ago via the study of conformally soft gluons and gravitons operator product expansions, in the context of celestial holography. Connections with twistor theory and phase space realization followed soon after. In the gravitational case for instance, these symmetries generalize the BMS algebra to include an infinite tower of symmetry generators constructed from tensors on the sphere of arbitrary high rank s. The bulk interpretation of these transformation parameters is still largely under investigation. Building on a first series of results on their canonical representation, we develop the necessary framework to define Noether charges for all degrees s. Importantly, we construct these charges out of a `holomorphic' asymptotic symplectic potential, such that we obtain an infinite collection of charges conserved in the absence of radiation. The classical symmetry is then realized non-perturbatively and non-linearly in the so-called holomorphic coupling constant, generalizing the perturbative linear and quadratic approach known so far. The infinitesimal action defines a symmetry algebroid which reduces to a symmetry algebra at non-radiative cuts of null infinity. The key ingredient for our construction is to consider field and time dependent symmetry parameters constrained to evolve according to equations of motion dual to (a truncation of) the asymptotic equations of motion in vacuum. We expose our results for Yang-Mills, General Relativity, and Einstein-Yang-Mills theories. This canonical analysis comes hand in hand with an in-depth study of the algebraic structure underlying the symmetry. We reveal several Lie algebroid and Lie algebra brackets, which connect the Carrollian, celestial and twistorial realizations. We show that these brackets are a deformation of the soft celestial algebra, where the deformation parameters are the radiative asymptotic data. On the one hand they allow us to define the symmetry algebra at non-radiative cuts of null infinity, for arbitrary values of the asymptotic shear and gauge potential. On the other hand, we can accommodate for radiation using the algebroid framework. For the specific case of non-abelian gauge theory, we also investigate how the asymptotic expansion around null infinity of the full Yang-Mills equations of motion in vacuum can be recast in terms of the higher spin charge aspects. Since the analysis of higher spin symmetries is for now inherent to a truncation of the latter equations of motion, this paves the way towards the understanding of the relevance of these symmetries in the full theory.
Statistical Mechanics of Finite Length Semiflexible Wormlike Polymers
(University of Waterloo, 2025-09-03) Andersen, Nigel
The wormlike chain is a fundamental model in polymer physics used for describing the statistics of semiflexible polymers. It is a general model applicable to many polymers, with double-stranded DNA being a notable example. Previous theoretical work has focused on long polymers, where subtle differences in the choice of statistical ensemble are not relevant. This worked well in explaining early data on polymer stretching, but as more modern techniques have allowed for the stretching of shorter, more rodlike polymers, finite length ensemble effects have become visible in experimental measurements and computer simulations. Here, these effects are computed, and a full picture of wormlike chain statistical mechanics is presented. The new results are compared against several decades of polymer stretching literature, and are found to describe previously unexplained behaviour. Finally, the broad use of the model is demonstrated by applying it to DNA wrapped carbon nanotubes, a system that requires considering polymer flexibility, and predicting the pitch of the helical structure formed by these complexes.
Characterization of prostate tumor spheroid growth and response to treatment using dynamic optical coherence tomography
(University of Waterloo, 2025-08-18) Swanson, Steph
Prostate cancer, the most prevalent cancer in North American men, is often treated with radiation and in the case of advanced disease, with the first-line chemotherapeutic docetaxel (DTX). The use of in vitro cell culture as a model for in vivo patient tumors has been instrumental in understanding the cancer biology underpinning tumor development, progression, and treatment. Conventionally, in vitro cells were cultured as 2D monolayers on hard, flat surfaces. However, it has become well understood in the last half-century that the behavior of in vivo tumors is better replicated by small 3D aggregates of in vitro cancer cells called tumor spheroids. While cells proliferate along the spheroid periphery, the center of the spheroid succumbs to starvation and waste build up. In between forms an intermediate layer of hypoxic and non-proliferative quiescent cells: two characteristics associated with treatment resistance. Despite the physiological pertinence of spheroid culture, its 3D nature challenges conventional biological methods. Moreover, cell culture geometry influences cell behavior not only during treatment, but in post-treatment recovery and throughout measurement as well. The clonogenic assay is a gold standard method for quantifying cell survival following treatment; however, it requires spheroids to be disaggregated and cultured in 2D for colony formation, which alters the cellular response. Proliferation assays quantify cellular activity inside intact spheroids, but similarly lack spatial resolution. Although fluorescence microscopy (FM) enables spatially resolved spheroid evaluation, it remains largely qualitative rather than quantitative. Proliferation assays and FM also require the addition of exogenous agents that are invasive to the sample and struggle to penetrate large spheroids. Alternatively, optical coherence tomography (OCT) enables non-invasive, high-speed, volumetric imaging of biological tissues with cellular resolution. Analysis of temporal OCT intensity fluctuations generates dynamic OCT (dOCT) images that can provide a quantitative and spatially resolved measurement of cellular activity throughout the spheroid. Nevertheless, investigating treatment response in spheroids is known to be tedious and time consuming, and some questions of interest cannot be accessed experimentally with adequate accuracy. As such, mathematical in silico models of spheroids have become increasingly prevalent. The heterogeneity inherent to in vivo tumors and in vitro spheroids is recapitulated by discrete in silico methods like agent-based modeling (ABM). These methods simulate each cancer cell as an individual “agent” that responds to neighboring cells, nutrients, drugs, and other components of its local environment. In particular, the cellular Potts model (CPM) is a form of ABM with phenomenological utility, given adequate validation to experimental observations. As such, it can describe, visualize, interpolate, and potentially even extrapolate experimental data to guide future experiments. This dissertation investigated and validated emerging and under-used in vitro tools for characterizing prostate tumor spheroid growth and response to treatment with radiation and DTX. First, I created a semi-automated masking process to isolate spheroids from volumetric OCT images for high resolution morphological analysis. For analysis of dynamic motion, I generated dOCT images using a frequency banding method that had only been previously subject to qualitative evaluation. Then, I developed a technique to quantify cellular activity in volumetric dOCT images of masked spheroids. Visual and quantitative comparison of live and formaldehyde-fixed spheroids imaged with dOCT and FM confirmed that dynamic motion measured with the dOCT method was associated with cellular activity. To mitigate influences unrelated to cellular activity, the average measurement of formaldehyde-fixed spheroids was subtracted from quantitative dOCT measurement as background. Post-treatment recovery of prostate tumor spheroids exposed to DTX was investigated with repeated measurement via Alamar Blue (AB) proliferation assay and longitudinal observation of clonogenic assay colony formation. Excellent agreement was observed between the quantitative dOCT method and AB proliferation assay over two weeks of longitudinal spheroid growth. Volumetric morphological analysis supported the measured cellular activity trends. However, fixation-subtraction could not be performed in a spatially sensitive manner and the dOCT images failed to resolve the longitudinal formation of a necrotic spheroid core that was observed via FM. Nonetheless, dOCT images and quantification of spheroids post-radiation demonstrated good agreement with FM and AB, respectively. Spheroids treated with radiation and DTX demonstrated better survival compared to monolayer culture. Monolayer culture treated with low dose DTX demonstrated higher post-treatment cellular activity and faster colony formation, but clonogenic survival remained lower than the untreated control. This effect was also observed in spheroids, albeit to a lesser extent. Experimental observations were probed with novel in silico CPMs of prostate tumor spheroid treatment. This thesis serves as a step towards validating emerging and under-used in vitro tools for spheroid evaluation in well-studied conventional cancer treatments. Once validated, these in vitro tools are particularly well suited for discovery and testing of novel targeted cancer treatments since spheroid protein and gene expressions are more physiologically representative than monolayer culture.
Quantum Error Correction and Quantum Metrology with Non-Markovian Noise
(University of Waterloo, 2025-08-12) Mann, Zachary
Quantum technologies have the potential to solve many important problems across science and industry. An important example is quantum computation. Quantum simulators promise to better model chemistry. Further, Shor’s factoring algorithm solves a problem in exponentially less time than what it would take now on our classical computers. This has led many to believe that quantum computers could bring exponential speedups to other difficult, real-world problems, such as optimization. Another example is quantum sensing, where quantum mechanical effects can be leveraged to increase measurement precision beyond the classical state of the art. This has many applications in both fundamental science, such as the LIGO experiment, and in industry, such as Nitrogen vacancy magnetometers. For these quantum technologies to reach their full potential, however, the barrier of noise must be overcome. Quantum effects usually live at very small system sizes or very cold temperatures, making them extra sensitive to thermal noise or small perturbations of the environment. A proposed solution to this problem, for both computation and sensing, is to use quantum error correction. Quantum error correction encodes a few quantum degrees of freedom into many physical degrees of freedom, building in redundancy. This redundancy allows for the detection and correction of unwanted errors in our protocol. Most of the literature on quantum error correction focuses on Markovian noise models, i.e., models where the noise is not temporally correlated. The temporally correlated, or non-Markovian, regime remains relatively unexplored. In this thesis, we explore quantum error correction for non-Markovian noise models. We first present a few of the many definitions and models for quantum non-Markovian phenomena present in the literature. We then generalize the Knill-Laflamme quantum error conditions to the hidden Markov model, an experimentally motivated model of non-Markovian noise. These conditions allow one to guarantee that a quantum error-correcting code will still do its job for more realistic noise models. Finally, we apply our notion of non-Markovian error correction to quantum sensing. We generalize previous Markovian results and derive conditions for guaranteeing Heisenberg limited precision scaling in the presence of temporally correlated noise using quantum error correction. The Heisenberg limit is the fundamental precision limit allowed by quantum mechanics for parameter estimation in a physical system. We also study the next-best achievable precision scaling when the Heisenberg limit is unattainable.
Quantum Monte Carlo Simulations of Rydberg Atom Arrays
(University of Waterloo, 2025-07-07) Merali, Ejaaz
Rydberg atom arrays form a promising platform for quantum computation. Through their strong, long-range interaction, they are able to encode various difficult combinatorial problems, as well as hosting a plethora of intriguing physical phenomena. In this thesis, we develop and apply a Stochastic Series Expansion Quantum Monte Carlo method to simulate Rydberg systems at zero-temperature and above. We then apply this simulation method alongside variational models to verify correctness of both methods. The data produced from the simulations is also used to train Neural Network wavefunctions, which we find are effectively able to grasp some of the physics of the Rydberg atom array on a square lattice.
A Machine Learning Model for Trapped-ion State Classification
(University of Waterloo, 2025-05-14) Balaniuk, Severyn
Academia and industry have been working to build a quantum computer that is able to perform certain tasks significantly better than classical computers. This thesis focuses on improving a trapped-ion-based approach to quantum computing. This platform has advanced significantly over the last 10 years, but there are numerous issues we need to resolve to make this architecture scalable. We consider experiments that use a high-sensitivity photon detection module as the readout tool. This setup lets us see qubit measurement outcomes as a fluorescent signal on a digital camera. This thesis addresses the problem of classifying fluorescence states for experimental systems of up to four ions, representing them as binary sequences. Our model hopefully will be able to classify systems with 8, 16, and 32 ions, allowing for further application of this methodology as quantum computers grow in scale. Our datasets are mainly unlabelled, which is the biggest challenge for training an accurate machine learning (ML) model. Nevertheless, we showed a significant improvement in classification accuracy with reduced bias over legacy models on unlabelled data containing mixed states of a four-ion system. We tested different machine learning architectures like feed-forward neural network (FFNN), convolutional neural networks (CNN), and semi-supervised learning to evaluate their efficiency for our specific dataset and tested their performance. In addition, we also developed and improved our own FFNN architecture with custom loss functions.
Towards Optical Simulation of Topological Phenomena with Ring Resonators
(University of Waterloo, 2025-04-30) Kuchhal, Bharat
In recent years, there has been a growing interest in synthetic dimensions. Unlike physical dimensions where one is restricted with three possible dimensions, one can tailor a lattice system with several higher order dimensions even with a lower-order physical system. In this thesis, we will explore, both theoretically and experimentally, how equidistant frequency modes in a resonator—specifically a ring resonator—can be used as a synthetic dimension. In this respect, we will first explore basics of resonators by reviewing a Fabry Perot resonator and developing a coupled-mode theory to understand its transmission and reflection spectrum. We will then expand and develop a coupled-mode theory for the case of a ring resonator, and explore experimental results for a basic fiber-beam splitter based ring resonator. Next, we will modify the ring resonator with a few off the shelf fiber optic components like a Dense Division Wavelength Multiplexing filter, a fiber electro-optic modulator, and an optical amplifier, and see how they transform the transmission characteristics. Lastly, we will modulate the ring resonator at the same frequency as the free-spectral range of the resonator to couple the modes to observe a band structure in the reciprocal space, much like atoms in a periodic solid-state crystal coupled via Coulomb interactions lead to energy bands in the conjugate momentum space. For this purpose, we will develop a coupled-mode theory for a ring resonator modulated by a phase modulator, look at some theoretical results for symmetric and asymmetric band structures, and observe some experimental results for the symmetric case. Thus, we highlight the potential of fiber ring resonators as versatile platforms for optical simulation of topological systems. Further, the ability to replicate higher-dimensional phenomena using frequency modes opens new pathways for investigating exotic physical behaviors, such as non-Hermitian systems and synthetic magnetic fields, in compact and experimentally accessible setups.
A Multi-Phase Analysis of Gas Dynamics and Perturbations in the Galaxy Cluster Cores
(University of Waterloo, 2025-02-14) Li, Muzi; McNamara, Brian
This thesis provides a detailed analysis of gas kinematics and their interactions across various phases within galaxy cluster cores. It examines the processes that generate gas perturbations and the factors that contribute to the thermal stability of the intracluster medium (ICM). A focus is placed on exploring the origins of multi-phase gas and the mechanisms—particularly AGN feedback—that either couple or decouple their motions. Radio-mechanical AGN feedback is identified as one of the most promising heating mechanisms that prevent the cooling of gas. However, the debate on the details of the heating transport processes has remained open. The atmospheres of 5 cool-core clusters, Abell 2029, Abell 2107, Abell 2151, RBS0533 and RBS0540, have short central cooling times but little evidence of cold gas, and jet-inflated bubbles. The amplitudes of gas density fluctuations were measured using a new statistical analysis of X-ray surface brightness fluctuations within the cool cores of these ‘spoil’ clusters in Chapter 2. The derived velocities of gas motions, typically around 100 - 200 km/s, are comparable to those in atmospheres around central galaxies experiencing energetic feedback, such as in the Perseus Cluster, and align well with the turbulent velocities expected in the ICM. Regardless of the mechanisms driving these perturbations, turbulent heating appears sufficient to counteract radiative losses in four of the five spoiler cluster cores. We thus suggest that other mechanisms, such as gas sloshing, may be responsible for generating turbulence, offering a plausible solution to suppress cooling in these structureless atmospheres. Multiphase filaments, key byproducts of AGN feedback, are frequently observed near central galaxies, with their morphologies and kinematics closely linked to bubbles. In Chapter 3, we analyzed the velocity structure functions (VSFs) of warm ionized gas and cold molecular gas, identified through [OII] emission and CO emissions observed by the Keck Cosmic Web Imager (KCWI) and the Atacama Large Millimeter/submillimeter Array (ALMA), respectively, in four clusters: Abell 1835, PKS 0745-191, Abell 262, and RXJ0820.9+0752. Excluding Abell 262, where gas forms a circumnuclear disk, the remaining clusters exhibit VSFs steeper than the Kolmogorov slope. The VSFs of CO and [OII] in RXJ0820 and Abell 262 show close alignment, whereas in PKS 0745 and Abell 1835, were differentiated across most scales, likely due to the churning caused by the radio-AGN. The large-scale consistency in Abell 1835 and RXJ0820, together with scale-dependent velocity amplitudes of the hot atmospheres obtained from Chandra X-ray data, may support the idea of cold gas condensation from the hot atmospheres. X-ray observations have previously been constrained by low energy resolution, which has impeded direct measurements of velocity fields in galaxy clusters. However, the recent release of initial data from the X-ray Imaging and Spectroscopy Mission (XRISM) provides a non-dispersive energy resolution of about 5 eV, facilitating the measurement of line broadening and shifts. In Chapter 4 of this thesis, I detail my contributions to calibrating the optical blocking filters for XRISM using synchrotron beamlines at the Canadian Light Source (CLS) and Advanced Light Source (ALS) prior to its launch, and I discuss the model-based estimation of the parameters of the calibrated filters. This capability for direct measurement of plasma velocities is expected to greatly improve our understanding of the ICM dynamics with high accuracy.
Investigating Abundances in Galaxy Clusters and Gas Motions in M87 using XRISM
(University of Waterloo, 2025-01-24) Dizdar, Neo; McNamara, Brian
Galaxy clusters are the forefront of extragalactic diffuse X-ray astrophysics, yet there are still many questions about their formation and evolution. The creation of XRISM, a new X-ray imaging and spectroscopy mission, will study the metal abundance history of clusters and the conversion of jet energy into atmospheric kinetic energy. XRISM’s payload contains an instrument with the highest spectral resolution (5 eV) in the field of X-ray astronomy so far. With this resolution, we observed metal abundances and the broadening of metal lines through turbulent motions in the intracluster medium. In this thesis I present the conversion of data from the Chandra X-ray Observatory to XRISM’s high-resolution format. This includes the preparation and selection of clusters in Chandra, simulating selected clusters for XRISM and applying for proposals. Finally, we extracted abundance and velocity information from the Virgo cluster’s early XRISM data.
On the Initial Boundary Value Problem in Numerical Relativity
(University of Waterloo, 2025-01-23) Dailey, Conner; Afshordi, Niayesh; Schnetter, Erik
The principal goal of this thesis is to properly understand, characterize, and numerically implement initial boundary value problems in numerical relativity. Throughout the history of solving Einstein's field equations on computers, boundaries have been mostly dealt with in an approximate way. For example, boundaries might be placed far away from strongly gravitating sources, where approximations like linearized gravity are valid. It has become necessary however to place boundaries in the strong gravity regime of a dynamical spacetime to model complicated and interesting physics, which necessitates a complete understanding of the initial boundary value problem of Einstein's field equations. One motivation for this comes from a need to simulate black hole echoes. In classical general relativity, black holes are perfectly absorbing objects, where the mass of radially incoming wavepackets of matter or gravitational waves is absorbed by the black hole. Thus conclusive evidence of modifications to general relativity, such as quantum gravity, could include partial reflections of radially incoming wavepackets, called black hole echoes. To properly understand the modifications this would bring to detectable gravitational wave signals, we require simulations where reflecting boundary conditions are imposed close to the horizon of a black hole. Another motivation comes from recent advances in Cauchy characteristic matching, which combines state of the art numerical techniques to obtain physically accurate gravitational waveforms from simulations. This can allow numerical relativists to dramatically save on the computational cost of black hole merger simulations, but only if boundaries can be placed in the strong gravity regime. This thesis presents advances in simulating initial boundary value problems in numerical relativity. Starting with spherical symmetry, a framework for reflecting a scalar field in a fully dynamical spacetime is developed and implemented numerically using the Einstein-Christoffel formulation. The evolution of a wave packet and its numerical convergence, including when the location of a reflecting boundary is very close to the horizon of a black hole, is studied. Next, this approach is generalized to spacetimes with no symmetries and implemented numerically using the generalized harmonic formulation. The evolution equations are cast into a summation by parts scheme, which seats the numerical method closer to a class of provably numerically stable systems. State of the art numerical methods are demonstrated, including an embedded boundary numerical method that allows for arbitrarily shaped domains on a rectangular grid and even boundaries that evolve and move across the grid. As a demonstration of these frameworks, the evolution of gravitational wave scattering off of a boundary either inside or just outside the horizon of a black hole, is studied. Finally, a boundary condition framework designed to control quasi-local energy flux is proposed motivated by examples from electromagnetism.
Studies on the 1/f noise in Shunted Josephson Junctions and Effective Masses in 2DHG GaAs Hall Bars
(University of Waterloo, 2025-01-22) Coschizza, Andree; Kycia, Jan
Two projects are presented with the motivation of improving materials and devices for superconducting and semiconducting qubit applications. The first project focuses on quantifying the 1/f noise in SIS Josephson junctions using a dc SQUID. The current-voltage characteristics of two junctions are measured and the degree of thermal rounding is analyzed as a function of temperature. The current-voltage characteristics are then used to compute the voltage bias dependence of both the 1/f and white noise of the junction. Finally, the temperature dependence of the 1/f noise is investigated and found to exhibit a linear relationship within the region of 100mK-4.5K. This is then compared to similar systems in literature. The second project explores a GaAs/AlGaAs 2DHG system using quantum Hall effect measurements. The mobility and charge carrier density dependence is calculated, and an experiment is conducted to investigate their dependence on highly negative top gates. It is discovered that the mobility-density curve increases as the top gate is pushed to more negative values. Furthermore, the spin-orbit interaction is parameterized and compared to previously measured samples. Finally, the band structure of GaAs/AlGaAs holes is discussed in reference to the observed Shubnikov-de Haas oscillations. The effective mass corresponding to the light heavy hole subband is measured using the Ando, Fowler, and Stern equation at low magnetic field oscillations.
Time-resolved and equilibrium resonant X-ray studies of nematicity in cuprate superconductors
(University of Waterloo, 2025-01-16) Gupta, Naman Kumar; Hawthorn, David
Attracting four decades of experimental and theoretical investigations since the discovery of high-temperature superconductivity, cuprates have become an essential benchmark for the exploration of intertwined electronic symmetry-breaking phases such as superconductivity, magnetism, nematicity, charge or spin density wave order in strongly correlated systems. Understanding how these phases are intertwined is of general and fundamental interest to a wide range of quantum materials. This doctoral work employs the element- and orbital-specific sensitivity of time-resolved and equilibrium resonant X-ray scattering to investigate electronic nematicity and charge density wave (CDW) order in cuprates. By probing their responses to photoexcitation using ultrafast laser pump-probe setups, varying hole doping levels, and applying compressive uniaxial stress, this thesis attempts to disentangle, understand, and manipulate intertwined phases in stripe-ordered cuprates.
Neurovascular Coupling in healthy human retina evaluated with Optical Coherence Tomography
(University of Waterloo, 2024-12-21) Dhaliwal, Khushmeet; Bizheva, Kostadinka
Retinal neurodegenerative diseases such as glaucoma, age-related macular degeneration (AMD), diabetic retinopathy, and retinitis pigmentosa affect millions of people worldwide and pose a significant burden on public health and the economy. Glaucoma, impacting approximately 80 million people globally, is a leading cause of irreversible blindness. In 2019, an estimated 19.8 million Americans (12.6%) were living with AMD, of which about 1.49 million people faced vision-threatening conditions. An estimated 9.6 million people were living with diabetic retinopathy in 2021, with about 1.84 million of them experiencing vision-threatening stages. In the United States alone, vision impairments- including those resulting from these retinal diseases- cost an estimated $139 billion annually. Retinal neurodegenerative diseases not only cause progressive damage to the retinal morphology and vascular network, but also cause acute and transient metabolic, physiological, and blood flow changes at the early stages of the disease development which become permanent and chronic at the advanced stages of the disease. Neurovascular Coupling (NVC) refers to the transient vasodilation and increased retinal blood flow resulting from the increased metabolic activity of retinal neurons in response to visual stimulation. Over the past few decades, a range of imaging techniques from clinical ophthalmoscope and confocal microscopy to adaptive optics scanning laser ophthalmoscope and optical coherence tomography (OCT) have been used ex vivo and in vivo to study components of the neurovascular coupling and its underlying mechanisms. Techniques such as Laser Doppler Velocimetry, Optical Coherence Tomography Angiography (OCTA), and Doppler Optical Coherence Tomography (D-OCT) have been used to observe the vascular responses of the retina caused by visual stimulation. Additionally, Electroretinography (ERG) has been widely used in clinical settings to evaluate the electrical activity of the neuronal retina. More recently, an optical equivalent to ERG, Optoretinography (ORG) was developed and OCT technology, imaging protocols, and image processing algorithms were designed to conduct OCT-based ORG studies in the human and animal retina. However, most of the Doppler OCT, OCTA, and ORG studies have examined components of the neurovascular coupling separately, potentially overlooking the dynamic interactions and comprehensive responses inherent in neurovascular coupling. OCT, which acquires simultaneously both intensity and phase information, is particularly well-suited for investigating neurovascular coupling in the retina, as it enables a completely non-invasive approach for simultaneous monitoring of retinal blood flow dynamics and neuronal responses. The integration of a commercial ERG system with a research-grade OCT modality adds further value by offering easy control of the visual stimulus, use of clinically established ERG protocols designed to elicit responses from specific types of retinal neuronal cells, and using the ERG recordings to validate the visually-evoked neuronal responses. The main objectives of this PhD thesis were: 1. To develop a combined OCT+ERG imaging system to conduct in vivo and simultaneously morphological and functional imaging that can be utilized for investigating neurovascular coupling in the human retina. 2. To evaluate the performance and capabilities of the OCT+ERG system, imaging protocols, and image processing algorithms by conducting a pilot study on healthy human subjects. 3. To utilize the OCT+ERG technology to explore the neurovascular coupling mechanisms in the healthy human retina by extracting vascular and neuronal responses from different retinal layers simultaneously. 4. To examine the effects of different wavelengths and flicker frequencies on the dynamic retinal blood flow changes evoked by visual stimulation, providing deeper insights into the mechanisms of neurovascular coupling. Results from this PhD research have been summarized in three manuscripts that are either under review or under preparation for submission. Therefore, this PhD thesis was prepared in such a way that individual manuscripts represent separate thesis chapters.
The Evolution of Halo Properties Through Binary Major Mergers
(University of Waterloo, 2024-10-17) Marchioni, Justin; Taylor, James
Galaxy clusters are massive, gravitationally bound objects composed of a large population of galaxies. Each of these galaxies occupies a dark matter halo and collectively the cluster has its own extended halo. Cluster halos can be described by many structural properties including their mass, concentration, shape, spin, and asymmetry. These properties, among others, can be used as proxies to constrain cosmology. The issue with galaxy clusters is that they are still assembling at the present-day. These clusters primarily grow through mergers, where smaller systems coalesce to form larger ones. If the mass ratio between the two merging components is sufficiently large (i.e. 3:1 or below), this is known as a major merger. The effect of major mergers is to significantly redistribute the matter distribution in the host system. This leads to pronounced fluctuations in the cluster's structural properties during the merger, making measurements of these properties hard to interpret. Therefore, accurately predicting how cluster properties vary during mergers is important in order to use them as a cosmological tool. In this thesis, we use simulations to study how the structure of remnant systems evolves during mergers. These simulations consider the merger of two isolated components, each represented by truncated Navarro-Frenk-White (NFW) profiles. We find that mergers produce oscillations in structural parameters for both the overall remnant and the host system. For example, the host halo's concentration experiences one of two types of responses to the satellite's motion depending primarily on the pericentric passage distance of the orbit. Given the simulation results, we present a semi-analytic model for the evolution of structure in remnant systems due to isolated, binary mergers. The model consists of two components, a treatment for the orbital evolution of the satellite and a prescription for changes in the host halo's potential. This second component is often neglected when modeling satellite orbits in minor mergers. Interestingly, we find that adding a host halo response model has little impact on the orbital evolution of the satellite and its mass loss. In contrast, this model must be incorporated in order to accurately predict how the remnant's structure changes after the satellite first passes pericentre. While our model generally works well at replicating the median concentration for the first two orbits, it is unable to recreate any of the remnant's anisotropy properties (i.e. shape, spin, and asymmetry). Overall, our results provide a framework for analyzing the response of cluster halo properties to mergers in more realistic scenarios.
Characterization and First Results of an Inverse Photoemission Spectrometer
(University of Waterloo, 2024-09-26) Bouliane, Michael; Hawthorn, David
This thesis discusses our efforts to characterize our home built inverse photoemission spectrometer. We review the relevant theoretical and practical considerations for the technique of inverse photoemssion spectroscopy. We then detail our efforts to characterize our low energy electron gun, presenting a method for determining the total current delivered by the beam using just a Faraday cup. Measurements of the beam’s profile are presented and are used to calculate the parallel momentum resolution of the beam. Equally important is the characterization of our photon detectors, which we show are operating in the proportional region with a minimal dark count rate. We ascertained a spectrometer energy resolution of 415(55) meV by performing inverse photoemission measurements on polycrystalline gold foil, single crystal Cu (111), and pyrolytic graphite. As a final demonstration of our spectrometer’s capabilities we provide a full unoccupied band mapping for pyrolytic graphite showing its ability to resolve dispersive electronic features in reciprocal space.
Towards Large Scale Quantum Simulations with Trapped Ions: Programmable XY model, Precise Light Sensing, and Extreme High Vacuum
(University of Waterloo, 2024-09-24) Kotibhaskar, Nikhil; Islam, Rajibul
We are currently witnessing a revolution in quantum technologies. Today's controllable quantum devices have reached a complexity that makes it practically intractable to fully simulate their dynamics using current classical supercomputers. Decades of fundamental research and development have led us to this point. In the coming years, billions of dollars in investments from governments and private entities are expected worldwide. Although general-purpose fault-tolerant quantum computers are expected to impact computing profoundly, today's quantum devices are best suited for their analog quantum operation, where a well-controlled quantum simulator mimics the dynamics of the other quantum system being studied. This affords an advantage over classical simulators at the cost of a restricted set of physical phenomena that can be studied. Today's quantum devices are already providing insights into large-scale entanglement, the underlying physics of high-temperature superconductivity, disordered quantum systems, and much more. Enhancing the capabilities of today's analog quantum simulators requires adding more classes of interactions, reducing errors due to calibration and noise, and increasing the system size to allow larger-scale simulations. The work described in this thesis directly addresses these core points for a system of trapped ions, which are ideal quantum simulators of the coupled dynamics of a large number of magnetic spins. First, the theory and experiment pertaining to the simulation of the anisotropic XY model on trapped ions has been presented. The theoretical proposal does not require added technical improvements over what has existed in the field for over a decade. The experimental validation is performed on a system with two 171Yb+ ions. This directly enhances the repertoire of trapped ions simulators and opens avenues to the exploration of high-temperature superconductivity, superfluidity, and spin liquids. The second result is the demonstration of the highest resolution readout of optical intensity and polarization using a single 171Yb+ ion as the field probe. The technique utilized the intensity- and polarization-dependent optical pumping of the ions as a signature to detect light parameters. This will be useful for the characterization of the optical addressing fields in trapped ion quantum simulators and hence for the calibration of large-scale quantum devices. Finally, the design and construction of a large-scale ion trapping apparatus for quantum simulation are described. The ion trap allows for the trapping of more than 50 ions, and the vacuum chamber used to house the trap with pressure below 1.5E-12 mbar (measurement limited by pressure gauge saturation) likely sets a record for the lowest pressure achieved on a room-temperature trapped ions system. This increases the useful simulation time of large-scale trapped-ion devices and paves the way for further enhancement of the scale of the simulations performed. Together, these results are another step in advancing the capabilities of today's quantum devices to explore physical phenomena far beyond the capability of classical supercomputers.
Development of the CHORD Galaxy Search Strategy
(University of Waterloo, 2024-09-24) Hopkins, Hans; Lang, Dustin; Taylor, James
This thesis presents my contribution to the CHORD galaxies science case. I helped develop an algorithm that can automatically pick out galaxies from CHORD driftscan data. The method used is a matched filter, and it acts on spatial data and spectral data. On the spatial side, it searches for point sources. Because CHORD is a highly redundant interferometer, it suffers from spatial aliasing. I present a program that is able to predict the severity of spatial aliasing. It predicts that integrating over periods of time is required to break the alias degeneracy, and that "dithering" CHORD (rotating it by a couple degrees) further helps in reducing the alias issue. I offer a framework for estimating the probability of spatial alias confusion. On the frequency side, I present a method of running the matched filter quickly. CHORD frequency data undergoes a process called upchannelization, which would distort the shape of a galaxy profile. I show how this can be accounted for without incurring a large time-cost penalty. Lastly, I discuss how a full matched filter program would be put together, and implications that my research has on selecting search parameters for a future CHORD galaxy survey.
Studying Unmodeled Physics from Gravitational Wave Data
(University of Waterloo, 2024-09-24) Dideron, Guillaume; Lehner, Luis
This thesis explores the detection and analysis of unmodeled physics in Gravitational Wave (GW) data. To this end, we develop the SCoRe framework, which uses the Correlated Residual Power Spectrum (CRPS) between pairs of detectors to identify deviations from our Standard Model (SM) of GW. This model includes General Relativity (GR) as the theory describing gravity, binary Black Holes (BHs) and Neutron Star (NS) merging as the sources of GWs, our model of the noise in the detectors, and the template waveform models used for data analysis. The thesis starts with a theoretical overview of GW physics, including an overview of GR, and how it describes the way GWs are generated and how they propagate and interact with matter. We then discuss the practical aspects of GW detection: the modelling of the noise in the detectors and the data analysis techniques used to extract and interpret GW signals. Next, we describe the SCoRe framework in Chapter 2, which is designed to distinguish between noise and deviations from the SM, while also shedding light on the underlying physics of the deviation. We detail its three main components: cross-correlating residual power between detectors, projecting onto physically motivated or agnostic bases, and combining information from multiple events by assuming a dependence of the unmodeled physics on the source parameters. To illustrate the method, we apply the SCoRe framework to toy models in Chapter 3. We demonstrate how the method can recover unmodeled signals without prior assumptions about their form, how to choose the timescale of cross-correlation, and how the method can be used to perform a null test of the SM. In Chapter 4, we then forecast the precision with which the SCoRe method can recover a deviation from the SM from a population of Binary Black Hole (BBH) mergers observed by a network of third-generation GW detectors. As the method leverages the dependence of the deviation from the SM on the source parameters, we investigate the effect the distribution of these parameters has on the method. For a model where the deviation scales with the chirp mass as a decaying power law, we show that the precision of the constraints on the deviation decreases as the power law becomes steeper. This has implications for constraining higher-dimensional operators in Effective Field Theories (EFTs) of gravity: higher dimensional operators correspond to steeper power laws and are, therefore, harder to constrain with the method. Finally, in Chapter 5, we illustrate another approach to testing the SM of GWs, where the GW signal in an alternative theory of gravity is numerically computed. We give an overview of the mathematical challenges and describe a method, the “fixing the equations” method, which aims to reduce pathologies in evolving EFTs of gravity by controlling energy flow to high frequencies.

Browse

Recent Submissions