Le projet LIGO a pour but la détection et l'étude d'ondes gravitationnelles via un réseau de détecteurs. LIGO possède deux détecteurs d'architecture et de fonctionnement identiques, situés aux États-Unis. Chaque détecteur est une version améliorée d'un interféromètre de Michelson avec des bras optiques de 4 km de long. Ces interféromètres ont observé une onde gravitationnelle pour la première fois en septembre 2015, suivi par cinq autres détections à ce jour. Ces détections marquent le début d™une nouvelle ère pour l'astrophysique, en liaison étroite avec la physique des trous noirs et des étoiles à neutrons. Depuis, un grand nombre d'activités sont en développement pour perfectionner les interféromètres. Cette thèse a pour objectif d'améliorer le temps de service des détecteurs, en répondant en particulier à deux problématiques majeures : le problème des impacts environnementaux, et notamment celui des tremblements de terre, ainsi que le problème lié à des couplages opto-mécaniques instables dans les cavités optiques, appelés instabilités paramétriques. Les stratégies de contrôle et les outils développés pour résoudre ces problématiques sont présentés. Les résultats prémilinaires montrent une réduction du temps d'arrêt généré par les tremblements de terre d'environ 40%. De plus, le dispositif ‚Acoustic Mode Damper™ développé pendant la thèse devrait complètement résoudre le problème des instabilités paramétriques pour LIGO. En conclusion, il sera démontré en quoi les problématiques résolues ont permis d'améliorer le cycle de service des détecteurs de LIGO de 4,6%, ce qui correspond à une augmentation du nombre d'ondes gravitationnelles détectées par an de 14%. / The LIGO project is a large-scale physics experiment the goal of which is to detect and study gravitational waves of astrophysical origin. It is composed of two instruments identical in design, located in the United States. The two instruments are specialized versions of a Michelson interferometer with 4km-long arms. They observed a gravitational-wave signal for the first time in September 2015 from the merger of two stellar-mass black holes. This is the first direct detection of a gravitational wave and the first direct observation of a binary black hole merger. Five more detections from binary black hole mergers and neutron stars merger have been reported to date, marking the beginning of a new era in astrophysics. As a result of these detections, many activities are in progress to improve the duty cycle and sensitivity of the detectors. This thesis addresses two major issues limiting the duty cycle of the LIGO detectors: environmental impacts, especially earthquakes, and the issue of unstable opto-mechanical couplings in the cavities, referred to as parametric instabilities. The control strategies and tools developed to tackle these issues are presented. Early results have shown a downtime reduction during earthquakes of ~40% at one of the LIGO sites. Moreover, the electro-mechanical device called ‚Acoustic Mode Damper™ designed and tested during the thesis should completely solve the issue of parametric instabilities for LIGO. In conclusion, we will show that the problems tackled in this thesis improved the overall duty cycle of LIGO by 4.6%, which corresponds to an increase of the gravitational-wave detection rate by 14%.

Virgo est un interféromètre de Michelson dont les bras contiennent des cavités Fabry-Perot. Il a été construit pour détecter directement les ondes gravitationnelles. Le projet Advanced Virgo est une amélioration majeure de Virgo pour atteindre une sensibilité encore plus élevée avec laquelle la détection des ondes gravitationnelles deviendra probable. On prévoit un système laser mono-fréquence de 175 Watts de puissance optique présentant des stabilités accrues pour le bruit relatif de puissance et pour le bruit de fréquence. Ce travail de thèse a pour objet la réalisation de ce système laser de haute-puissance et de haute-stabilité basée sur l'utilisation d'amplificateurs à fibre combinés de façon cohérente. Des amplificateurs à fibre disponibles dans le commerce sont caractérisés en termes de qualité de faisceau, de bruit de puissance, de bruit de fréquence, de stabilité de pointé du faisceau, et également en terme de stabilité à long terme sur quelques milliers d'heures. On implémente l'interférométrie de Mach-Zehnder pour la combinaison cohérente de faisceaux. Les techniques de caractérisation de faisceaux laser sont aussi développées en considérant leurs limites ultimes. Hormis un déficit de puissance optique, le système laser développé dans cette étude sur la base de la combinaison cohérente de Master Oscillator Fiber Power Amplifiers, remplit les conditions posées par Advanced Virgo. / Virgo is a cavity-enhanced Michelson interferometer built for the direct detection of gravitational waves. The Advanced Virgo project consists of major upgrades to the Virgo gravitational wave detector for an order of magnitude improvement in differential strain sensitivity, one of which is the tenfold increase in injected laser power to 175 Watts. The use of fiber laser amplifiers and their coherent combination are foreseen to deliver the required high-power low-noise beam. In this thesis work, we review the laser requirements for gravitational wave detectors, introduce the design of the laser system for Advanced Virgo, and develop the means for laser characterization in accordance with the stringent noise specifications. We then present the results to date, notably the quasi-continuous long-term operation of two 40-Watt fiber laser amplifiers over thousands of hours and their coherent combination with Mach-Zehnder interferometry. Although the targeted power for Advanced Virgo is not yet attained, the developed system shows decent noise performance and is promising for further power-scaling efforts.

En décembre 2015, les détecteurs LIGO ont pour la première fois détecté une onde gravitationnelle émise lors de la coalescence d'une paire de trous noirs il y a de cela 1.3 milliards d'années. Une telle première dans la toute nouvelle astronomie gravitationnelle a été suivie par plusieurs autres observations. La dernière en date est la fusion de deux étoiles à neutron dont la contrepartie électromagnétique a pu être observée par plusieurs observatoires à travers le monde. A cette occasion, les ondes gravitationnelles se sont inscrites dans l'astronomie multi-messager. Ces observations ont été rendues possibles par des techniques avancées d'analyse de données. Grâce à elles, la faible empreinte laissée par une onde gravitationnelle dans les données de détecteurs peut être isolée. Le travail de cette thèse est dédié au développement d'une technique de détection d'ondes gravitationnelles ne reposant que sur une connaissance minimale du signal à isoler. Le développement de cette méthode consiste plus précisément à introduire une information sur la phase du signal d'onde gravitationnelle selon un contexte astrophysique déterminé. La première partie de cette thèse est consacrée à la présentation de la méthode. Dans une seconde partie cette méthode est appliquée à la recherche de signaux d'ondes gravitationnelles en provenance de systèmes binaires de trous noirs de masse stellaire dans du bruit Gaussien. Puis l'étude est répétée dans du bruit de détecteurs collecté pendant la première période de prise de données. Enfin la troisième partie est dédiée à la recherche de binaires de trous noirs dont l'orbite montre un écart à la géométrie circulaire, ce qui complexifie la morphologie du signal. De telles orbites sont qualifiées d'excentriques. Cette troisième analyse permet d'établir de premiers résultats quant à la méthode proposée lorsque le signal d'intérêt est peu connu / In december 2015 the LIGO detectors have first detected a gravitational wave emitted by a pair of coalescing black holes 1.3 billion years ago. Many more observations have been realised since then and heralded gravitational waves as a new messenger in astronomy. The latest detection is the merge of two neutron stars whose electromagnetic counterpart has been followed up by many observatories around the globe. These direct observations have been made possible by the developpement of advanced data analysis techniques. With them the weak gravitational wave inprint in detectors may be recovered. The realised work during this thesis aims at developping an existing gravitational wave detection method which relies on minimal assumptions of the targeted signal. It more precisely consists in introducing an information on the signal phase depending on the astrophysical context. The first part is dedicated to a presentation of the method. The second one presents the results obtained when applying the method to the search of stellar mass binary black holes in simulated Gaussian noise data. The study is repeated in real instrumental data collected during the first run of LIGO. Finally, the third part presents the method applied in the search for eccentric binary black holes. Their orbit exhibits a deviation from the quasi-circular orbit case considered so far and thus complicates the signal morphology. This third analysis establishes first results with the proposed method in the case of a poorly modeled signal

Binary systems comprising of compact objects like neutron stars (NS) and/or black holes (BH) lose their energy and angular momentum via gravitational waves (GW). Radiation reaction due to the emission of GW results in a gradual shrinking of the binary orbit and an accompanying gradual increase in the orbital frequency. The preliminary phase of the binary evolution when the radiation-reaction time-scale is much larger than the orbital time-scale is called the inspiral phase. GW emitted during the final stages of the inspiral phase constitute one of the most important sources for the ground-based laser interferometric GW detectors like LIGO, VIRGO and the proposed space-based detector LISA. For the ground-based detectors, NS and/or stellar mass BH binaries are primary sources, while for LISA super-massive BH (SMBH) binaries are potential targets. Inspiralling compact binaries (ICB) are among the prime targets for interferometric detectors because using approximation schemes in general relativity (GR) like the post-Minkowskian (PM) and the post-Newtonian (PN) approximations one can compute the GW emitted by them with sufficient accuracy both for their detection and parameter estimation leading to GW astronomy. The extreme weakness of gravitational interactions implies that if a GW signal from an ICB is incident on a detector, it will be buried in the noisy detector output. Therefore, sophisticated data analysis techniques are required for detecting the signal in presence of the dominant noise and also estimating the parameters of the signal. From the pre-calculated theoretical waveforms called templates, one already knows the structure of the waveform from an ICB. The technique for detecting signals which are of known form in a noisy detector is matched filtering. This technique consists of cross-correlating the output of a noisy detector assumed to contain the signal of known form with a set of templates. It then finds an ‘optimal’ template that would produce, on average, the highest signal-to-noise ratio (SNR). The efficient performance of matched filtering as a data-analysis strategy for GW signals from ICB presupposes very accurate theoretical templates. Slight mismatches between the signal and the template will result in a loss of signal to noise ratio. Computing very accurate theoretical templates and including effects such as eccentricity are challenging tasks for the theoreticians. This thesis addresses some of the issues related to the waveform modelling of the ICB and their implications for GW data analysis. It is known theoretically that compact binaries reduce their eccentricity through the emission of GW. When GW signals from prototype ICB reach the GW detector bandwidth, their orbits are almost circular. Hence one usually models the binary orbit to be circular for computation of the search templates. The waveform from an ICB in a circular orbit is, at any given PN order of approximation, a linear combination of a finite number of harmonics of the orbital frequency. At the lowest order of approximation, called the Newtonian order, the waveform comprises a single harmonic at twice the orbital frequency. Inclusion of higher order PN corrections lead to the appearance of higher harmonics of the orbital frequency. Since the amplitudes of the higher harmonics contain higher powers of the PN expansion parameter, relative to the Newtonian order, they are referred to as amplitude corrections. The phase of each harmonic, determined by the orbital phase, is known upto 3.5PN order (nPN is the order of approximation equivalent to terms ~(v/c)2n beyond the Newtonian order, where v denotes the binary’s orbital velocity and c is the speed of light). Matched filtering is more sensitive to the phase of the signal rather than its amplitude, since the correlation builds up as long as the signal and the template remain in phase. Motivated by this fact, search templates so far have been a waveform model involving only the dominant harmonic (at twice the orbital frequency), although the phase evolution itself is included upto the maximum available PN order. Such waveforms, in which all amplitude corrections are neglected, but the phase is treated to the maximum available order, are called restricted waveforms (RWF) and these are generally used in the data-analysis of ground-based detectors and also simulated searches for the planned LISA. However, recent studies, in the case of ground-based interferometers, showed that going beyond the RWF approximation could improve the efficiency of detection as well as parameter estimation of the inspiral signal. After a brief overview of the properties of GW and their detection strategies in chapter 1, in chapters 2 and 3, we investigate the implications of going beyond the RWF, in the context of the planned space-based Laser Interferometric Space Antenna (LISA). The sensitivity of ground-based detectors is limited by seismic noise below 20Hz. On the other hand, the space-based LISA will be designed to be sensitive to GWs of frequency (10−4 _1)Hz. The most important source in this frequency band are supermassive BH (SMBH) binaries. There is strong observational evidence for the existence of SMBH with masses in the range of in most galactic nuclei. Mergers of such galaxies result in SMBH binaries whose evolution is governed by the emission of GW. Observation of the GW from SMBH binaries at high redshifts is one of the major science goals of LISA. These observations will allow us to probe the evolution of SMBHs and structure formation and provide an unique opportunity to test General Relativity (and its alternatives) in the strong field regime of the theory. Observing SMBH coalescences with high (100-1000) SNR is crucial for performing all the aforementioned tests. The LISA bandwidth (10−4_ 1)Hz determines the range of masses accessible to LISA because the inspiral signal would end when the system’s orbital frequency reaches the mass-dependent last stable orbit (LSO). In the test-mass approximation, the angular velocity ι at LSO is given by where M is the total mass of the binary. Search templates using the RWF, which contains only the dominant harmonic at twice the orbital frequency, cannot extract power in the signal beyond This further implies that the frequency range [0.1, 100] mHz corresponds to the range for the total mass of BH binaries that would be accessible to LISA. In chapter 2, we show that inclusion of higher harmonics will enhance the mass-range of LISA (for the same frequency range) and allow for the detection of SMBH binaries with total masses higher than The template employed in chapter 2 includes amplitude corrections upto 2.5PN order, while keeping the phase upto 3.5PN order. We call this template the full waveform (FWF). The FWF defined above contains higher harmonics of the orbital frequency, the highest of them being 7 times the orbital frequency. For a SMBH binary with total mass the dominant harmonic at LSO is less than the lower cut-off of the LISA bandwidth. Therefore, if one uses the RWF as a search template, this system is ‘invisible’ to LISA. However, the seventh harmonic can still enter the LISA bandwidth and produce a significant SNR and thus allow its detection. With the FWF, LISA can observe sources which are favoured by astronomical observations, but not observable with the RWF. More specifically, with the inclusion of all known harmonics LISA will be able to observe SMBH coalescences with total mass (and mass-ratio 0.1) for a low frequency cut-off of 10−4Hz (10−5Hz) with an SNR up to ~ 60 (~30) at a distance of 3 Gpc. The orbital motion of LISA around the Sun induces frequency, phase and amplitude modulations in the observed GW signal. These modulations carry information about both the source’s location and orientation. Determination of the angular coordinates of the source also allows determination of the luminosity distance of SMBH binaries. Therefore, SMBH binaries are often referred to as GW “standard sirens” (analogous to the electromagnetic “standard candles”). LISA would also be able to measure the “redshifted” masses of the component black holes with good accuracy for sources up to redshifts of a few. However, GW observations alone cannot provide any information about the redshift of the source. If the host galaxy or galaxy cluster is known one can disentangle the redshift from the masses by optical measurement of the redshift. This would not only allow one to extract the “physical” masses, but also provide an exciting possibility to study the luminosity distance-redshift relation providing a totally independent confirmation of the cosmological parameters. Further, this combined observation can be used to map the distribution of black hole masses as a function of redshift. Another outstanding issue in present day cosmology in which LISA can play a role is the dark energy and its physical origin. Probing the equation-of-state-ratio (w(z)) provides an important clue to the question of whether dark energy is truly a cosmological constant (i.e., w = -1). Assuming the Universe to be spatially flat, a combination of WMAP and Supernova Legacy Survey (SNLS) data yields significant constraints on Without including the spatial flatness as a prior, WMAP, large-scale structure and supernova data place a stringent constraint on the dark energy equation of state, For this to be possible, LISA should (a) measure the luminosity distance to the source with a good accuracy and (b) localize the coalescence event on the sky with good angular resolution so that the host galaxy/galaxy cluster can be uniquely identified. Based on analysis with the RWF, it is found that LISA’s angular resolution is not good enough to identify the source galaxy or galaxy cluster, and that other forms of identification would be needed. Secondly, weak lensing effects would corrupt the distance estimation to the same level as LISA’s systematic error. In chapter 3, we study the problem of parameter estimation in the context of LISA, but using the FWF. We investigate systematically the variation in parameter estimation with PN orders by critically examining the role of higher harmonics in the fast GW phasing and their interplay with the slow modulations induced due to LISA’s motion. More importantly, we explore the improvement in the estimation of the luminosity distance and the angular parameters due to the inclusion of higher harmonics in the waveform. We translate the error in the angular resolution to obtain the number of galaxies (or galaxy clusters) within the error box on the sky. We find that independent of the angular position of the source on the sky, higher harmonics improve LISA’s performance on both counts raised in earlier works based on the RWF. We show that the angular resolution enhances typically by a factor of ~2-500 (greater at higher masses) and the error on the estimation of the luminosity distance goes down by a factor of ~ 2-100 (again, larger at higher masses). For many possible sky positions and orientations of the source, the inaccuracy in our measurement of the dark energy would be at the level of a few percent, so that it would only be limited by weak lensing. We conclude that LISA could provide interesting constraints on cosmological parameters, especially the dark energy equation-of-state, and yet circumvent all the lower rungs of the cosmic distance ladder. Having emphasized the need to consider the FWF as a more powerful template, in chapter 4 we calculate a higher order term in the amplitude corrections of the waveform. In chapters 2 & 3, the FWF incorporated amplitude corrections upto 2.5PN order. In chapter 4 the waveform is calculated upto 3PN order. Recent progress in Numerical Relativity (NR) has resulted in computation of the late inspiral and subsequent merger and ringdown phases of the binary evolution (where PN theory does not hold good) by a full-fledged numerical integration of the Einstein field equations. A new field has emerged recently consisting of high-accuracy comparisons between the PN predictions and the numerically-generated waveforms. Such comparisons and matching to the PN results have proved currently to be very successful. They clearly show the need to include high PN corrections not only for the evolution of the binary’s orbital phase but also for the modulation of the gravitational amplitude. This leads to one more motivation for the work in this chapter: providing the associated spin-weighted spherical harmonic decomposition to facilitate comparison and match of the high PN prediction for the inspiral waveform to the numerically-generated waveforms for the merger and ringdown. For the computation of waveforms from the inspiralling compact binaries one needs to solve the two-body problem in general relativity. The nonlinear structure of general relativity prevents one from obtaining a general solution to this problem. The two-body problem is tackled using the multipolar post-Minkowskian (MPM) wave generation formalism. The MPM formalism describes the radiation field of any isolated post-Newtonian source. The radiation field is first of all parametrized by means of two sets of radiative multipole moments. These moments are then related (by means of an algorithm for solving the non-linearities of the field equations) to the so-called canonical moments which constitute some useful intermediaries for describing the external field of the source. The canonical moments are then expressed in terms of the operational source moments obtained by matching to a PN source and are given by explicit integrals extending over the matter source and gravitational field. The extension of the waveform by half a PN order requires as inputs the relations between the radiative, canonical and source multipole moments for general sources at 3PN order. We also require the 3PN extension of the source multipole moments in the case of compact binaries. The waveform in the far-zone consists of two types of terms, instantaneous and hereditary. The instantaneous terms are determined by the dynamical state of the binary at the retarded time. The hereditary terms, on the other hand, depend on the entire past history of the source. These terms originate from the nonlinear interactions between the various multipole moments and also from backscattering off the curved spacetime generated by the waves themselves. In this chapter, we compute the contributions of all the instantaneous and hereditary terms (which include tails, tails-of-tails and memory integrals) up to 3PN order. The end results of this chapter are given in terms of both the 3PN plus and cross polarizations and the separate spin-weighted spherical harmonic modes. Though most of the sources will be in circular orbits by the time the GWs emitted by the system enter the sensitivity band of the laser interferometers, astrophysical scenarios such as Kozai mechanism could produce binaries which have nonzero eccentricity. Studies have shown that filtering the signal from an eccentric binary with circular orbit templates could significantly degrade the SNR. For constructing a phasing formula for eccentric binaries one has to compute the energy and angular momentum fluxes carried away by the GWs and then compute how the orbital elements evolve with time under gravitational radiation reaction. The far-zone energy and angular momentum fluxes, like the waveform, contain both instantaneous and hereditary contributions. The complete 3PN energy flux and instantaneous terms in the 3PN angular momentum flux are already known. In chapter 5, the hereditary terms in the 3PN angular momentum flux from an ICB moving in quasi-elliptical orbits are computed. A semi-analytic method in the frequency domain is used to compute the hereditary contributions. At 3PN order, the quasi-Keplerian representation of elliptical orbits at 1PN order is required. To calculate the tail contributions we exploit the doubly periodic nature of the motion to average the 3PN fluxes over the binary’s orbit. The hereditary part of the angular momentum flux provided here has to be supplemented with the instantaneous part to obtain the final input needed for the construction of templates for binaries moving in elliptical orbits, a class of sources for both the space based detectors and the ground based ones. Using the hereditary contributions in the 3PN energy flux, we also compute the 3PN accurate hereditary contributions to the secular evolution of the orbital elements of the quasi-Keplerian orbit description.

Ce mémoire présente les efforts effectués pour réduire l'angle de perte de couches de pentoxyde de tantale amorphes telles qu'utilisées pour les miroirs de LIGO. Afin d'améliorer le niveau de relaxation des couches, celles-ci ont été déposées par pulvérisation cathodique magnétron à des températures allant de 50 °C à 480 °C, elles ont subi un recuit thermique rapide, elles ont été implantées par des ions d'oxygène, elles ont été déposées par pulvérisation cathodique magnétron en appliquant une tension de polarisation sur le substrat lors du dépôt allant de 0 V à -450 V et elles ont été déposées par pulvérisation cathodique magnétron pulsée à haute puissance dans le cadre de différentes expériences. L'angle de perte, l'épaisseur, la rugosité, l'indice de réfraction, la composition atomique, la contrainte, l'état de relaxation et le module de Young des couches ont par la suite été trouvés à l'aide de l'ellipsométrie spectralement résolue, la spectrométrie de rétrodiffusion de Rutherford, la détection des reculs élastiques, la spectroscopie Raman, la diffraction de rayons X et la nano-indentation. Il a été trouvé que la température de dépôt améliorait légèrement le degré de relaxation des couches jusqu'à 250 °C, mais qu'elle avait peu d'impact après recuit. Aussi, lors de dépôt à température de la pièce, une forte tension de polarisation réduit l'angle de perte, mais cet effet est encore une fois perdu suite au recuit. Les autres méthodes mentionnées ci-dessus n'ont pas influencé le degré de relaxation des couches selon l'angle de perte, la spectroscopie Raman et la diffraction de rayons X. Cette recherche a été réalisée avec le support financier du CRSNG et du FRQNT (numéro de dossier 206976). / This master's thesis presents the experiments made to reduce the loss angle of tantala coatings similar to those used in LIGO. To improve the relxation level of the coatings they were deposited by magnetron sputtering at temperatures varying from 50 °C to 480 °C. They were also subjected to rapid thermal annealing, and oxygen implantation. In another experiment, the coatings were deposited by magnetron sputtering with substrate biasing varying from 0 V to -450 V at room temperature and at 250 °C. Finally, the coatings of tantala were deposited by high power impulse magnetron sputtering. The loss angle, thickness, roughness, refractive index, atomic composition, stress, the relaxation state and Young's modulus of the coatings were characterized using spectroscopic ellipsometry, Rutherford backscattering, elastic recoil detection, Raman spectroscopy, X-ray diffraction and nanoindentation. It was found that the deposition temperature improved the loss angle until it reached 250 °C. However, annealing the coatings had a superior impact and the influence of the deposition temperature was not visible after annealing. When was applied a high bias to the susbtrate at room temperature, the obtained coating was slightly more relaxed than when a low bias was applied but this effect is, once again, insignificant after annealing. The other methods of deposition mentioned did not improve the loss angle or modify the relaxation state found by Raman spectroscopy and X-ray diffraction of the tantala coatings. This research was made with the financial support of the NSERC and of the FRQNT (file number 206976).