Lock Acquisition and Sensitivity Analysis of Advanced LIGO Interferometers

Laser interferometer gravitational wave observatory (LIGO) consists of two complex large-scale laser interferometers designed for direct detection of gravitational waves from distant astrophysical sources in the frequency range 10Hz - 5kHz. Direct detection of space-time ripples will support Einstein's general theory of relativity and provide invaluable information and new insight into physics of the Universe.

Initial phase of LIGO started in 2002, and since then data was collected during six science runs. Instrument sensitivity was improving from run to run due to the effort of commissioning team. Initial LIGO has reached designed sensitivity during the last science run, which ended in October 2010.

In parallel with commissioning and data analysis with the initial detector, LIGO group worked on research and development of the next generation detectors. Major instrument upgrade from initial to advanced LIGO started in 2010 and lasted till 2014.

This thesis describes results of commissioning work done at LIGO Livingston site from 2013 until 2015 in parallel with and after the installation of the instrument. This thesis also discusses new techniques and tools developed at the 40m prototype including adaptive filtering, estimation of quantization noise in digital filters and design of isolation kits for ground seismometers.

The first part of this thesis is devoted to the description of methods for bringing interferometer to the linear regime when collection of data becomes possible. States of longitudinal and angular controls of interferometer degrees of freedom during lock acquisition process and in low noise configuration are discussed in details.

Once interferometer is locked and transitioned to low noise regime, instrument produces astrophysics data that should be calibrated to units of meters or strain. The second part of this thesis describes online calibration technique set up in both observatories to monitor the quality of the collected data in real time. Sensitivity analysis was done to understand and eliminate noise sources of the instrument.

Coupling of noise sources to gravitational wave channel can be reduced if robust feedforward and optimal feedback control loops are implemented. The last part of this thesis describes static and adaptive feedforward noise cancellation techniques applied to Advanced LIGO interferometers and tested at the 40m prototype. Applications of optimal time domain feedback control techniques and estimators to aLIGO control loops are also discussed.

Commissioning work is still ongoing at the sites. First science run of advanced LIGO is planned for September 2015 and will last for 3-4 months. This run will be followed by a set of small instrument upgrades that will be installed on a time scale of few months. Second science run will start in spring 2016 and last for about 6 months. Since current sensitivity of advanced LIGO is already more than factor of 3 higher compared to initial detectors and keeps improving on a monthly basis, upcoming science runs have a good chance for the first direct detection of gravitational waves.

FDTD analysis of Talbot effect of a high density grating

Talbot effect of a grating with different flaws is analyzed with the finite-difference time-domain (FDTD) method. The FDTD method can show the exact near-field distribution of different flaws in a high-density grating, which is impossible to obtain with the conventional Fourier transform method. The numerical results indicate that if a grating is perfect, its Talbot imaging should also be perfect; if the grating is distorted, its Talbot imaging would also be distorted. Furthermore, we can evaluate high density gratings by detecting the near-field distribution.

Rigorous transmittance analysis of a sub-wavelength Sb thin film lens

We quantitatively analysed the factors contributing to the optical transmission enhancement of a sub-wavelength Sb thin film lens, using the finite-difference time-domain (FDTD) method. The results show that the transmission enhancement of the dielectric with a Gaussian distributed refractive index loaded in a sub-wavelength circular hole is not only due to the high refractive index dielectric, but also due to the specific distributions of refractive index. It is the first study about the effects of the refractive index distribution on the transmission of a sub-wavelength aperture. This kind of lens has practical applications in the very small aperture lasers and for near-field optical storage and lithography.

A coupled hydro-mechanical analysis for prediction of hydraulic fracture propagation in saturated porous media using EFG mesh-less method

The details of the Element Free Galerkin (EFG) method are presented with the method being applied to a study on hydraulic fracturing initiation and propagation process in a saturated porous medium using coupled hydro-mechanical numerical modelling. In this EFG method, interpolation (approximation) is based on nodes without using elements and hence an arbitrary discrete fracture path can be modelled.The numerical approach is based upon solving two governing partial differential equations of equilibrium and continuity of pore water simultaneously. Displacement increment and pore water pressure increment are discretized using the same EFG shape functions. An incremental constrained Galerkin weak form is used to create the discrete system of equations and a fully implicit scheme is used for discretization in the time domain. Implementation of essential boundary conditions is based on the penalty method. In order to model discrete fractures, the so-called diffraction method is used.Examples are presented and the results are compared to some closed-form solutions and FEM approximations in order to demonstrate the validity of the developed model and its capabilities. The model is able to take the anisotropy and inhomogeneity of the material into account. The applicability of the model is examined by simulating hydraulic fracture initiation and propagation process from a borehole by injection of fluid. The maximum tensile strength criterion and Mohr-Coulomb shear criterion are used for modelling tensile and shear fracture, respectively. The model successfully simulates the leak-off of fluid from the fracture into the surrounding material. The results indicate the importance of pore fluid pressure in the initiation and propagation pattern of fracture in saturated soils. © 2013 Elsevier Ltd.

Analysis of coupled microcircular resonators coupled to a bus waveguide with high output efficiency

Coupled microcircular resonators tangentially coupled to a bus waveguide, which is between the resonators, are numerically investigated by the finite-difference time-domain technique. For symmetrically coupled microcircular resonators with refractive index of 3.2, radius of 2 mu m, and width of the bus waveguide of 0.4 mu m, a mode Q factor of the order of 105 is obtained for a mode at the frequency of 243 THz. An output coupling efficiency of as high as 0.99 is calculated for a mode with a Q factor ranging from 10(3) to 10(4). The mode Q factor is 2 orders larger than that of the modes confined in a single circular resonator tangentially coupled to the same bus waveguide. Furthermore, the high Q traveling modes in the coupled microcircular resonators are suitable for optical single processing.

Analysis of mode characteristics for equilateral-polygonal resonators with a center hole

Mode characteristics for equilateral triangles, squares, and hexagonal resonators with a center hole are numerically simulated by the finite-different time domain (FDTD) technique. The center hole does not break the symmetry behavior of the original resonators and can result in modification of the mode field patterns and mode Q factors. In an equilateral triangle resonator the center hole can suppress the symmetry state of degenerate states with the merit of single mode operation. In a square resonator, the Q factor can be enhanced for some modes with a suitable size of the hole. For a hexagonal resonator with a side length of 1 mu m and a refractive index of 3.2, the mode Q factors first gradually decrease with the increase of the hole diameter for modes at a wavelength of about 1500 nm, then the modes transform to that of a microdisk with a jump of the mode wavelength as the hole diameter approaches 0.7 mu m. Finally, the mode Q factors greatly enhance as the hole diameter reaches about 1 mu m. The results indicate that the center hole can greatly modify mode characteristics, especially that of the mode Q factor. (C) 2009 Optical Society of America

Mode Analysis for Equilateral-Triangle-Resonator Microlasers with Metal Confinement Layers

Mode characteristics are analyzed for electrically injected equilateral-triangle-resonator (ETR) semiconductor microlasers, which are laterally confined by insulating barrier SiO2 and electrode metals Ti-Au. For the ETR without metal layers, the totally confined mode field patterns are derived based on the reflection phase shifts, and the Q-factors are calculated from the far-field emission of the analytical near field distribution, which are agreement very well with the numerical results of the finite-difference time-domain (FDTD) simulation. The polarization dependence reflections for light rays incident on semiconductor-SiO2 -Ti-Au multi-layer structures are accounted in considering the confinement of TE and TM modes in the ETR with the metal layers. The reflectivity will greatly reduce with a Ti layer between SiO2 and Au for light rays with incident angle less than 30 especially for the TE mode, even the thickness of the Ti layer is only 10 nm. If the ETR is laterally confined by SiO2-Au layers without the Ti layer, the Fabry-Perot type modes with an incident angle of zero on one side of the ETR can also have high Q-factor. The FDTD simulation for the ETR confined by metal layers verifies the above analysis based on multi-layer reflections. The output spectra with mode intervals of whispering-gallery modes and Fabry-Perot type modes are observed from different ETR lasers with side length of 10 m, respectively.

Analysis of mode characteristics for equilateral triangle semiconductor microlasers with imperfect boundaries

The influence of imperfect boundaries on the mode quality factor is investigated for equilateral-triangle-resonator (ETR) semiconductor microlasers by the finite difference time domain technique and the Pade approximation with Baker's algorithm. For 2-D ETR with a refractive index of 3.2 and side length of 5 mum, the confined modes can still have a quality factor of about 1000 as small triangles with side length of 1 mum are cut from the vertices of the ETR. For a deformed 5 mum ETR with round vertices and curve sides, the simulated mode quality factors are comparable to the measured results.

Symmetry analysis and numerical simulation of mode characteristics for equilateral-polygonal optical microresonators

Mode characteristics for two-dimensional equilateral-polygonal microresonators are investigated based on symmetry analysis and finite-difference time-domain numerical simulation. The symmetries of the resonators can be described by the point group C-Nv, accordingly, the confined modes in these resonators can be classified into irreducible representations of the point group C-Nv. Compared with circular resonators, the modes in equilateral-polygonal resonators have different characteristics due to the break of symmetries, such as the split of double-degenerate modes, high field intensity in the center region, and anomalous traveling-wave modes, which should be considered in the designs of the polygonal resonator microlasers or optical add-drop filters.

Mode analysis and Q-factor enhancement due to mode coupling in rectangular resonators

Modes in rectangular resonators are analyzed and classified according to symmetry properties, and quality factor (Q-factor) enhancement due to mode coupling is observed. In the analysis, mode numbers p and q are used to denote the number of wave nodes in the direction of two orthogonal sides. The even and odd mode numbers correspond to symmetric and antisymmetric field distribution relative to the midlines of sides, respectively. Thus, the modes in a rectangle resonator can be divided into four classes according to the parity of p and q. Mode coupling between modes of different classes is forbidden; however, anti-crossing mode coupling between the modes in the same class exists and results in new modes due to the combination of the coupled modes. One of the combined modes has very low power loss and high Q-factor based on far-field emission of the analytical field distribution, which agrees well with the numerical results of the finite-difference time-domain (FDTD) simulation. Both the analytical and FDTD results show that the Q-factors of the high Q-factor combined modes are over one order larger than those of the original modes. Furthermore, the general condition required to achieve high-Q modes in the rectangular resonator is given based on the analytical solution.

Analysis of mode characteristics for deformed square resonators by FDTD technique

The mode frequencies and quality factors (Q-factors) in two-dimensional (2-D) deformed square resonators are analyzed by finite-difference time-domain (FDTD) technique. The results show that the deformed square cavities with circular and cut corners have larger Q-factors than the perfect ones at certain conditions. For a square cavity with side length of 2 mu m and refractive index of 3.2, the mode Q-factor can increase 13 times as the perfect corners are replaced by a quarter of circle with radius of 0.3 pm. Furthermore the blue shift with the increasing deformations is found as a result of the reduction in effective resonator area. In square cavities with periodic roughness at sidewalls which maintains the symmetry of the square, the Q-factors of the whisperin gallery (WG)-like modes are still one order of magnitude larger that those of non-WG-like modes. However, the Q-tactors of these two types of modes are of the same order in the square cavity with random roughness. We also find that the rectangular and rhombic deformation largely reduce the Q-factors with the increasing offset and cause the splitting of the doubly degenerate modes due to the breaking of certain symmetry properties.

Analysis of mode quality factors and mode reflectivities for nanowire cavity by FDTD technique

The mode frequency and the quality factor of nanowire cavities are calculated from the intensity spectrum obtained by the finite-difference time-domain (FDTD) technique and the Pade approximation. In a free-standing nanowire cavity with dielectric constant epsilon = 6.0 and a length of 5 mu m, quality factors of 130, 159, and 151 are obtained for the HE11 modes with a wavelength around 375 nm, at cavity radius of 60, 75, and 90 nm, respectively. The corresponding quality factors reduce to 78, 94, and 86 for a nanowire cavity standing on a sapphire substrate with a refractive index of 1.8. The mode quality factors are also calculated for the TE01 and TM01 modes, and the mode reflectivities are calculated from the mode quality factors.

Analysis of mode quality factors for equilateral triangle semiconductor microlasers with rough sidewalls

The eigenmode characteristics for equilateral triangle resonator (ETR) semiconductor microlasers are analysed by the finite-difference time-domain technique and the Pade approximation. The random Gaussian correlation function and sinusoidal function are used to model the side roughness of the ETR. The numerical results show that the roughness can cause the split of the degenerative modes, but the confined modes can still have a high quality factor. For the ETR with a 3 mum side length and the sinusoidal fluctuation, we can have a quality factor of 800 for the fundamental mode in the wavelength of 1500 nm, as the amplitude of roughness is 75 mn.

Analysis of semiconductor microlasers with an equilateral triangle resonator by rate equations

Semiconductor microlasers with an equilateral triangle resonator (ETR) are analyzed by rate equations with the mode lifetimes calculated by the finite-difference time-domain technique and the Pade approximation. A gain spectrum based on the relation of the gain spectrum and the spontaneous emission spectrum is proposed for considering the mode selection in a wide wavelength span. For an ETR microlaser with the side length of about 5 mum, we find that single fundamental mode operation at about 1.55 mum can be obtained as the side length increases from 4.75 to 5.05 mum. The corresponding wavelength tuning range is 93 nm, and the threshold current is about 0.1 to 0.4 mA.

Analysis and numerical simulation of eigenmode characteristics for semiconductor lasers with an equilateral triangle micro-resonator

The eigenmodes confined in the equilateral triangle resonator (ETR) are analyzed by deriving the eigenvalues and the mode field distributions and by the finite difference time domain (FDTD) technique. The analytical results show that the one-period-length for the mode light rays inside the ETR is the perimeter of the ETR, and the number of transverse modes is limited by the condition of total internal reflection. In addition, the sum of the longitudinal mode index and the transverse mode index should be an even number, which limits the number of confined modes again. Based on the FDTD technique and the Pade approximation, we calculate the mode resonant frequencies and the quality factors from the local maximum and the width of the spectral distribution of the intensity The numerical results of mode frequencies agree very well with the analytical results, and the quality factor of the fundamental mode is usually higher than that of the higher order transverse modes. The results show that the ETR is suitable to realize single-made operation as semiconductor microcavity lasers.