Atmospheric extinction coefficients in the Ic\mathrm{I_c} band for several major international observatories: Results from the BiSON telescopes, 1984 to 2016

Over 30 years of solar data have been acquired by the Birmingham Solar Oscillations Network (BiSON), an international network of telescopes used to study oscillations of the Sun. Five of the six BiSON telescopes are located at major observatories. The observational sites are, in order of increasing longitude: Mount Wilson (Hale) Observatory (MWO), California, USA; Las Campanas Observatory (LCO), Chile; Observatorio del Teide, Iza\~{n}a, Tenerife, Canary Islands; the South African Astronomical Observatory (SAAO), Sutherland, South Africa; Carnarvon, Western Australia; and the Paul Wild Observatory, Narrabri, New South Wales, Australia. The BiSON data may be used to measure atmospheric extinction coefficients in the $\mathrm{I_c}$ band (approximately 700-900 nm), and presented here are the derived atmospheric extinction coefficients from each site over the years 1984 to 2016.Comment: 15 pages, 10 figures, 4 tables. Accepted by Astronomical Journal: 2017 July 2

Fixing the BMS Frame of Numerical Relativity Waveforms with BMS Charges

The Bondi-van der Burg-Metzner-Sachs (BMS) group, which uniquely describes the symmetries of asymptotic infinity and therefore of the gravitational waves that propagate there, has become increasingly important for accurate modeling of waveforms. In particular, waveform models, such as post-Newtonian (PN) expressions, numerical relativity (NR), and black hole perturbation theory, produce results that are in different BMS frames. Consequently, to build a model for the waveforms produced during the merging of compact objects, which ideally would be a hybridization of PN, NR, and black hole perturbation theory, one needs a fast and robust method for fixing the BMS freedoms. In this work, we present the first means of fixing the entire BMS freedom of NR waveforms to match the frame of either PN waveforms or black hole perturbation theory. We achieve this by finding the BMS transformations that change certain charges in a prescribed way -- e.g., finding the center-of-mass transformation that maps the center-of-mass charge to a mean of zero. We find that this new method is 20 times faster, and more correct when mapping to the superrest frame, than previous methods that relied on optimization algorithms. Furthermore, in the course of developing this charge-based frame fixing method, we compute the PN expression for the Moreschi supermomentum to 3PN order without spins and 2PN order with spins. This Moreschi supermomentum is effectively equivalent to the energy flux or the null memory contribution at future null infinity $\mathscr{I}^{+}$. From this PN calculation, we also compute oscillatory ($m\not=0$ modes) and spin-dependent memory terms that have not been identified previously or have been missing from strain expressions in the post-Newtonian literature. <br

Improving Gravitational Wave Extraction Using Weyl Characteristic Fields

We present a detailed methodology for extracting the full set of Newman-Penrose Weyl scalars from numerically generated spacetimes without requiring a tetrad that is completely orthonormal or perfectly aligned to the principal null directions. We also describe an extrapolation technique for computing the Weyl scalars' contribution at asymptotic null infinity in post-processing. These methods have continued to be used to produce $\Psi_4$ and $h$ waveforms for the Simulating eXtreme Spacetimes Waveform Catalog and now have been expanded to produce the entire set of Weyl scalars. These new waveform quantities are critical for the future of gravitational wave astronomy in order to understand the finite-amplitude gauge differences that can occur in numerical waveforms. We also present a new analysis of the accuracy of waveforms produced by the Spectral Einstein Code. While ultimately we expect Cauchy Characteristic Extraction to yield more accurate waveforms, the extraction techniques described here are far easier to implement and have already proven to be a viable way to produce production-level waveforms that can meet the demands of current gravitational-wave detectors

Surfaces roughness effects on the transmission of Gaussian beams by anisotropic parallel plates

Influence of the plate surfaces roughness in precise ellipsometry experiments is studied. The realistic case of a Gaussian laser beam crossing a uniaxial platelet is considered. Expression for the transmittance is determined using the first order perturbation theory. In this frame, it is shown that interference takes place between the specular transmitted beam and the scattered field. This effect is due to the angular distribution of the Gaussian beam and is of first order in the roughness over wavelength ratio. As an application, a numerical simulation of the effects of quartz roughness surfaces at normal incidence is provided. The interference term is found to be strongly connected to the random nature of the surface roughness.Comment: 18 pages, Journal of Physics D: Applied Physics, volume 36, issue 21, pages 2697 - 270

Fixing the BMS Frame of Numerical Relativity Waveforms

Understanding the Bondi-Metzner-Sachs (BMS) frame of the gravitational waves produced by numerical relativity is crucial for ensuring that analyses on such waveforms are performed properly. It is also important that models are built from waveforms in the same BMS frame. Up until now, however, the BMS frame of numerical waveforms has not been thoroughly examined, largely because the necessary tools have not existed. In this paper, we show how to analyze and map to a suitable BMS frame for numerical waveforms calculated with the Spectral Einstein Code (SpEC). However, the methods and tools that we present are general and can be applied to any numerical waveforms. We present an extensive study of 13 binary black hole systems that broadly span parameter space. From these simulations, we extract the strain and also the Weyl scalars using both SpECTRE's Cauchy-characteristic extraction module and also the standard extrapolation procedure with a displacement memory correction applied during post-processing. First, we show that the current center-of-mass correction used to map these waveforms to the center-of-mass frame is not as effective as previously thought. Consequently, we also develop an improved correction that utilizes asymptotic Poincar\'e charges instead of a Newtonian center-of-mass trajectory. Next, we map our waveforms to the post-Newtonian (PN) BMS frame using a PN strain waveform. This helps us find the unique BMS transformation that minimizes the $L^{2}$ norm of the difference between the numerical and PN strain waveforms during the early inspiral phase. We find that once the waveforms are mapped to the PN BMS frame, they can be hybridized with a PN strain waveform much more effectively than if one used any of the previous alignment schemes, which only utilize the Poincar\'e transformations

High Precision Ringdown Modeling: Multimode Fits and BMS Frames

Quasi-normal mode (QNM) modeling is an invaluable tool for studying strong gravity, characterizing remnant black holes, and testing general relativity. To date, most studies have focused on the dominant $(2, 2)$ mode, and have fit to standard strain waveforms from numerical relativity. But, as gravitational wave observatories become more sensitive, they can resolve higher-order modes. Multimode fits will be critically important, and in turn require a more thorough treatment of the asymptotic frame at $\mathscr{I}^+$. The first main result of this work is a method for systematically fitting a QNM model containing many modes to a numerical waveform produced using Cauchy-characteristic extraction, which is known to exhibit memory effects. We choose the modes to model based on their power contribution to the residual between numerical and model waveforms. We show that the all-angles mismatch improves by a factor of $\sim 10^5$ when using multimode fitting as opposed to only fitting $(2, \pm2, n)$ modes. Our second main result addresses a critical point that has been overlooked in the literature: the importance of matching the Bondi-van der Burg-Metzner-Sachs (BMS) frame of the simulated gravitational wave to that of the QNM model. We show that by mapping the numerical relativity waveforms to the super rest frame, there is an improvement of $\sim 10^5$ in the all-angles strain mismatch, compared to using the strain whose BMS frame is not fixed. We illustrate the effectiveness of these modeling enhancements by applying them to families of waveforms produced by numerical relativity, and comparing our results to previous studies

High-accuracy numerical models of Brownian thermal noise in thin mirror coatings

Brownian coating thermal noise in detector test masses is limiting the sensitivity of current gravitational-wave detectors on Earth. Therefore, accurate numerical models can inform the ongoing effort to minimize Brownian coating thermal noise in current and future gravitational-wave detectors. Such numerical models typically require significant computational resources and time, and often involve closed-source commercial codes. In contrast, open-source codes give complete visibility and control of the simulated physics and enable direct assessment of the numerical accuracy. In this article, we use the open-source SpECTRE numerical-relativity code and adopt a novel discontinuous Galerkin numerical method to model Brownian coating thermal noise. We demonstrate that SpECTRE achieves significantly higher accuracy than a previous approach at a fraction of the computational cost. Furthermore, we numerically model Brownian coating thermal noise in multiple sub-wavelength crystalline coating layers for the first time. Our new numerical method has the potential to enable fast exploration of realistic mirror configurations, and hence to guide the search for optimal mirror geometries, beam shapes and coating materials for gravitational-wave detectors

Adding Gravitational Memory to Waveform Catalogs using BMS Balance Laws

Accurate models of gravitational waves from merging binary black holes are crucial for detectors to measure events and extract new science. One important feature that is currently missing from the Simulating eXtreme Spacetimes (SXS) Collaboration's catalog of waveforms for merging black holes, and other waveform catalogs, is the gravitational memory effect: a persistent, physical change to spacetime that is induced by the passage of transient radiation. We find, however, that by exploiting the Bondi-Metzner-Sachs (BMS) balance laws, which come from the extended BMS transformations, we can correct the strain waveforms in the SXS catalog to include the missing displacement memory. Our results show that these corrected waveforms satisfy the BMS balance laws to a much higher degree of accuracy. Furthermore, we find that these corrected strain waveforms coincide especially well with the waveforms obtained from Cauchy-characteristic extraction (CCE) that already exhibit memory effects. These corrected strain waveforms also evade the transient junk effects that are currently present in CCE waveforms. Lastly, we make our code for computing these contributions to the BMS balance laws and memory publicly available as a part of the python package $\texttt{sxs}$, thus enabling anyone to evaluate the expected memory effects and violation of the BMS balance laws

Nonlinearities in black hole ringdowns

The gravitational wave strain emitted by a perturbed black hole (BH) ringing down is typically modeled analytically using first-order BH perturbation theory. In this Letter we show that second-order effects are necessary for modeling ringdowns from BH merger simulations. Focusing on the strain's $(\ell,m)=(4,4)$ angular harmonic, we show the presence of a quadratic effect across a range of binary BH mass ratios that agrees with theoretical expectations. We find that the quadratic $(4,4)$ mode amplitude exhibits quadratic scaling with the fundamental $(2,2)$ mode -- its parent mode. The nonlinear mode's amplitude is comparable to or even larger than that of the linear $(4,4)$ modes. Therefore correctly modeling ringdown -- improving mismatches by an order of magnitude -- requires the inclusion of nonlinear effects

High-accuracy waveforms for black hole-neutron star systems with spinning black holes

The availability of accurate numerical waveforms is an important requirement for the creation and calibration of reliable waveform models for gravitational wave astrophysics. For black hole-neutron star binaries, very few accurate waveforms are however publicly available. Most recent models are calibrated to a large number of older simulations with good parameter space coverage for low-spin non-precessing binaries but limited accuracy, and a much smaller number of longer, more recent simulations limited to non-spinning black holes. In this paper, we present long, accurate numerical waveforms for three new systems that include rapidly spinning black holes, and one precessing configuration. We study in detail the accuracy of the simulations, and in particular perform for the first time in the context of BHNS binaries a detailed comparison of waveform extrapolation methods to the results of Cauchy Characteristic Extraction. The new waveforms have $0.99$) for binaries seen face-on. For edge-on observations, particularly for precessing systems, disagreements between models and simulations increase, and models that include precession and/or higher-order modes start to perform better than BHNS models that currently lack these features