-
Path Optimization Sheaves
Authors:
Michael Moy,
Robert Cardona,
Robert Green,
Jacob Cleveland,
Alan Hylton,
Robert Short
Abstract:
Motivated by efforts to incorporate sheaves into networking, we seek to reinterpret pathfinding algorithms in terms of cellular sheaves, using Dijkstra's algorithm as an example. We construct sheaves on a graph with distinguished source and sink vertices, in which paths are represented by sections. The first sheaf is a very general construction that can be applied to other algorithms, while the se…
▽ More
Motivated by efforts to incorporate sheaves into networking, we seek to reinterpret pathfinding algorithms in terms of cellular sheaves, using Dijkstra's algorithm as an example. We construct sheaves on a graph with distinguished source and sink vertices, in which paths are represented by sections. The first sheaf is a very general construction that can be applied to other algorithms, while the second is created specifically to capture the decision making of Dijkstra's algorithm. In both cases, Dijkstra's algorithm can be described as a systematic process of extending local sections to global sections. We discuss the relationship between the two sheaves and summarize how other pathfinding algorithms can be interpreted in a similar way. While the sheaves presented here address paths and pathfinding algorithms, we suggest that future work could explore connections to other concepts from graph theory and other networking algorithms. This work was supported by the NASA Internship Project and SCaN Internship Project during the summer of 2020.
△ Less
Submitted 10 December, 2020;
originally announced December 2020.
-
On the one dimensional Dirac equation with potential
Authors:
Burak Erdogan,
William R. Green
Abstract:
We investigate $L^1\to L^\infty$ dispersive estimates for the one dimensional Dirac equation with a potential. In particular, we show that the Dirac evolution satisfies the natural $t^{-\frac12}$ decay rate, which may be improved to $t^{-\frac32}$ at the cost of spatial weights when the thresholds are regular. We classify the structure of threshold obstructions, showing that there is at most a one…
▽ More
We investigate $L^1\to L^\infty$ dispersive estimates for the one dimensional Dirac equation with a potential. In particular, we show that the Dirac evolution satisfies the natural $t^{-\frac12}$ decay rate, which may be improved to $t^{-\frac32}$ at the cost of spatial weights when the thresholds are regular. We classify the structure of threshold obstructions, showing that there is at most a one dimensional space at each threshold. We show that, in the presence of a threshold resonance, the Dirac evolution satisfies the natural decay rate, and satisfies the faster weighted bound except for a piece of rank at most two, one per threshold. Further, we prove high energy dispersive bounds that are near optimal with respect to the required smoothness of the initial data. To do so we use a variant of a high energy argument that was originally developed to study Kato smoothing estimates for magnetic Schrödinger operators. This method has never been used before to obtain $L^1 \to L^\infty$ estimates. As a consequence of our analysis we prove a uniform limiting absorption principle, Strichartz estimates and prove the existence of an eigenvalue free region for the one dimensional Dirac operator with a non-self-adjoint potential.
△ Less
Submitted 18 July, 2023; v1 submitted 26 November, 2020;
originally announced November 2020.
-
Topological waves in passive and active fluids on curved surfaces: a unified picture
Authors:
Richard Green,
Jay Armas,
Jan de Boer,
Luca Giomi
Abstract:
We investigate the occurrence of topologically protected waves in classical fluids confined on curved surfaces. Using a combination of topological band theory and real space analysis, we demonstrate the existence of a system-independent mechanism behind topological protection in two-dimensional passive and active fluids. This allows us to formulate an index theorem linking the number of modes, det…
▽ More
We investigate the occurrence of topologically protected waves in classical fluids confined on curved surfaces. Using a combination of topological band theory and real space analysis, we demonstrate the existence of a system-independent mechanism behind topological protection in two-dimensional passive and active fluids. This allows us to formulate an index theorem linking the number of modes, determined by the topology of Fourier space, to the real space topology of the surface on which they are hosted. With this framework in hand, we review two examples of topological waves in two-dimensional fluids, namely oceanic shallow-water waves propagating on the Earth's rotating surface and momentum waves in active polar fluids spontaneously "flocking" on substrates endowed with a ${\rm U}(1)$ isometry (e.g. surfaces of revolution). Our work suggests some simple rules to engineer topological modes on surfaces in passive and active soft matter systems.
△ Less
Submitted 24 November, 2020;
originally announced November 2020.
-
Placing High-Redshift Quasars in Perspective: a Catalog of Spectroscopic Properties from the Gemini Near Infrared Spectrograph -- Distant Quasar Survey
Authors:
Brandon M. Matthews,
Ohad Shemmer,
Cooper Dix,
Michael S. Brotherton,
Adam D. Myers,
I. Andruchow,
W. N. Brandt,
Gabriel A. Ferrero,
S. C. Gallagher,
Richard Green,
Paulina Lira,
Richard M. Plotkin,
Gordon T. Richards,
Jessie C. Runnoe,
Donald P. Schneider,
Yue Shen,
Michael A. Strauss,
Beverley J. Wills
Abstract:
We present spectroscopic measurements for 226 sources from the Gemini Near Infrared Spectrograph - Distant Quasar Survey (GNIRS-DQS). Being the largest uniform, homogeneous survey of its kind, it represents a flux-limited sample ($m_{i}$ ${\lesssim}$ 19.0 mag, $H$ ${\lesssim}$ 16.5 mag) of Sloan Digital Sky Survey (SDSS) quasars at 1.5 ${\lesssim}$ $z$ ${\lesssim}$ 3.5 with a monochromatic luminos…
▽ More
We present spectroscopic measurements for 226 sources from the Gemini Near Infrared Spectrograph - Distant Quasar Survey (GNIRS-DQS). Being the largest uniform, homogeneous survey of its kind, it represents a flux-limited sample ($m_{i}$ ${\lesssim}$ 19.0 mag, $H$ ${\lesssim}$ 16.5 mag) of Sloan Digital Sky Survey (SDSS) quasars at 1.5 ${\lesssim}$ $z$ ${\lesssim}$ 3.5 with a monochromatic luminosity ($λL_λ$) at 5100 ${\unicode{xC5}}$ in the range of $10^{44} - 10^{46}$ erg $\rm{s}^{-1}$. A combination of the GNIRS and SDSS spectra covers principal quasar diagnostic features, chiefly the C IV $λ$1549, Mg II $λλ$2798, 2803, H$β$ $λ$4861, and [O III] $λλ$4959, 5007 emission lines, in each source. The spectral inventory will be utilized primarily to develop prescriptions for obtaining more accurate and precise redshifts, black hole masses, and accretion rates for all quasars. Additionally, the measurements will facilitate an understanding of the dependence of rest-frame ultraviolet-optical spectral properties of quasars on redshift, luminosity, and Eddington ratio, and test whether the physical properties of the quasar central engine evolve over cosmic time.
△ Less
Submitted 21 November, 2020;
originally announced November 2020.
-
Strain induced stabilization of a static Jahn-Teller distortion in the O$^*$-phase of La$_{7/8}$Sr$_{1/8}$MnO$_3$
Authors:
Michael Dettbarn,
Volodymyr B. Zabolotnyy,
Abdul V. Tcakaev,
Ronny Sutarto,
Feizhou He,
Yunzhong Chen,
Robert J. Green,
Vladimir Hinkov
Abstract:
At room temperature, bulk La$_{7/8}$Sr$_{1/8}$MnO$_3$ is in the dynamic Jahn--Teller O$^*$ phase, but undergoes a transition to a static, magnetically ordered Jahn--Teller phase at lower temperatures. Here we study a $6$ unit cells thin film of this compound grown on SrTiO$_3$, resulting in small compressive strain due to a lattice mismatch of $\lesssim 0.2\%$. We combine X-ray absorption spectros…
▽ More
At room temperature, bulk La$_{7/8}$Sr$_{1/8}$MnO$_3$ is in the dynamic Jahn--Teller O$^*$ phase, but undergoes a transition to a static, magnetically ordered Jahn--Teller phase at lower temperatures. Here we study a $6$ unit cells thin film of this compound grown on SrTiO$_3$, resulting in small compressive strain due to a lattice mismatch of $\lesssim 0.2\%$. We combine X-ray absorption spectroscopy with multiplet ligand field theory to study the local electronic and magnetic properties of Mn in the film. We determine the Mn $d_{3z^2-r^2}$ orbital to be $0.13\;\text{eV}$ lower in energy than the $d_{x^2-y^2}$, which is a disproportionately large splitting given the small degree of compressive strain. We interpret this as resulting from the strain providing a preferential orientation for the MnO$_6$ octahedra, which are strongly susceptible to such a deformation in the vicinity of the phase transition. Hence, they collectively elongate along the $c$ axis into a static Jahn--Teller arrangement. Furthermore, we demonstrate the strongly covalent character of La$_{7/8}$Sr$_{1/8}$MnO$_3$, with a contribution of nearly $50\%$ of the one-ligand-hole configuration $d^{5} \underline{L}^1$ to the ground state wavefunction. Finally, we find the system to be in a high-spin configuration, with the projection of the local magnetic moment on the quantization axis being about $3.7\;μ_{\text{B}}/\text{Mn}$. We show, however, that the system is close to a high-spin--low-spin transition, which might be triggered by crystal field effects.
△ Less
Submitted 13 November, 2020;
originally announced November 2020.
-
Evidence for bond-disproportionation in LiNiO2 from x-ray absorption spectroscopy
Authors:
R. J. Green,
H. Wadati,
T. Z. Regier,
A. J. Achkar,
C. McMahon,
J. P. Clancy,
H. A. Dabkowska,
B. D. Gaulin,
G. A. Sawatzky,
D. G. Hawthorn
Abstract:
The electronic structure of LiNiO$_2$, a promising Li-ion battery cathode material, has remained a challenge to understand due to its highly covalent yet correlated nature. Here we elucidate the electronic structure in LiNiO$_2$ and the related compound NaNiO$_2$ using x-ray absorption spectra (XAS) and quantum many-body calculations. Notably, we use inverse partial fluorescence yield to correctly…
▽ More
The electronic structure of LiNiO$_2$, a promising Li-ion battery cathode material, has remained a challenge to understand due to its highly covalent yet correlated nature. Here we elucidate the electronic structure in LiNiO$_2$ and the related compound NaNiO$_2$ using x-ray absorption spectra (XAS) and quantum many-body calculations. Notably, we use inverse partial fluorescence yield to correctly measure the Ni $L$-edge XAS, which is inaccurate using conventional methods. We show that the XAS are indicative of a strong Jahn-Teller effect in NaNiO$_2$ and a bond disproportionated state in LiNiO$_2$, supporting a theory of a high-entropy, glassy disproportionated state that stabilizes charging cycles in LiNiO$_2$.
△ Less
Submitted 12 November, 2020;
originally announced November 2020.
-
Lattice continuum-limit study of nucleon quasi-PDFs
Authors:
Constantia Alexandrou,
Krzysztof Cichy,
Martha Constantinou,
Jeremy R. Green,
Kyriakos Hadjiyiannakou,
Karl Jansen,
Floriano Manigrasso,
Aurora Scapellato,
Fernanda Steffens
Abstract:
The quasi-PDF approach provides a path to computing parton distribution functions (PDFs) using lattice QCD. This approach requires matrix elements of a power-divergent operator in a nucleon at high momentum and one generically expects discretization effects starting at first order in the lattice spacing $a$. Therefore, it is important to demonstrate that the continuum limit can be reliably taken a…
▽ More
The quasi-PDF approach provides a path to computing parton distribution functions (PDFs) using lattice QCD. This approach requires matrix elements of a power-divergent operator in a nucleon at high momentum and one generically expects discretization effects starting at first order in the lattice spacing $a$. Therefore, it is important to demonstrate that the continuum limit can be reliably taken and to understand the size and shape of lattice artifacts. In this work, we report a calculation of isovector unpolarized and helicity PDFs using lattice ensembles with $N_f=2+1+1$ Wilson twisted mass fermions, a pion mass of approximately 370 MeV, and three different lattice spacings. Our results show a significant dependence on $a$, and the continuum extrapolation produces a better agreement with phenomenology. The latter is particularly true for the antiquark distribution at small momentum fraction $x$, where the extrapolation changes its sign.
△ Less
Submitted 13 April, 2021; v1 submitted 2 November, 2020;
originally announced November 2020.
-
Search for Gravitational Waves Associated with Gamma-Ray Bursts Detected by Fermi and Swift During the LIGO-Virgo Run O3a
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
R. Abbott,
T. D. Abbott,
S. Abraham,
F. Acernese,
K. Ackley,
C. Adams,
R. X. Adhikari,
V. B. Adya,
C. Affeldt,
M. Agathos,
K. Agatsuma,
N. Aggarwal,
O. D. Aguiar,
A. Aich,
L. Aiello,
A. Ain,
P. Ajith,
G. Allen,
A. Allocca,
P. A. Altin,
A. Amato,
S. Anand,
A. Ananyeva
, et al. (1228 additional authors not shown)
Abstract:
We search for gravitational-wave transients associated with gamma-ray bursts detected by the Fermi and Swift satellites during the first part of the third observing run of Advanced LIGO and Advanced Virgo (1 April 2019 15:00 UTC - 1 October 2019 15:00 UTC). 105 gamma-ray bursts were analyzed using a search for generic gravitational-wave transients; 32 gamma-ray bursts were analyzed with a search t…
▽ More
We search for gravitational-wave transients associated with gamma-ray bursts detected by the Fermi and Swift satellites during the first part of the third observing run of Advanced LIGO and Advanced Virgo (1 April 2019 15:00 UTC - 1 October 2019 15:00 UTC). 105 gamma-ray bursts were analyzed using a search for generic gravitational-wave transients; 32 gamma-ray bursts were analyzed with a search that specifically targets neutron star binary mergers as short gamma-ray burst progenitors. We describe a method to calculate the probability that triggers from the binary merger targeted search are astrophysical and apply that method to the most significant gamma-ray bursts in that search. We find no significant evidence for gravitational-wave signals associated with the gamma-ray bursts that we followed up, nor for a population of unidentified subthreshold signals. We consider several source types and signal morphologies, and report for these lower bounds on the distance to each gamma-ray burst.
△ Less
Submitted 20 August, 2021; v1 submitted 27 October, 2020;
originally announced October 2020.
-
Population Properties of Compact Objects from the Second LIGO-Virgo Gravitational-Wave Transient Catalog
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
R. Abbott,
T. D. Abbott,
S. Abraham,
F. Acernese,
K. Ackley,
A. Adams,
C. Adams,
R. X. Adhikari,
V. B. Adya,
C. Affeldt,
M. Agathos,
K. Agatsuma,
N. Aggarwal,
O. D. Aguiar,
L. Aiello,
A. Ain,
P. Ajith,
G. Allen,
A. Allocca,
P. A. Altin,
A. Amato,
S. Anand,
A. Ananyeva
, et al. (1316 additional authors not shown)
Abstract:
We report on the population of the 47 compact binary mergers detected with a false-alarm rate 1/yr in the second LIGO--Virgo Gravitational-Wave Transient Catalog, GWTC-2. We observe several characteristics of the merging binary black hole (BBH) population not discernible until now. First, we find that the primary mass spectrum contains structure beyond a power-law with a sharp high-mass cut-off; i…
▽ More
We report on the population of the 47 compact binary mergers detected with a false-alarm rate 1/yr in the second LIGO--Virgo Gravitational-Wave Transient Catalog, GWTC-2. We observe several characteristics of the merging binary black hole (BBH) population not discernible until now. First, we find that the primary mass spectrum contains structure beyond a power-law with a sharp high-mass cut-off; it is more consistent with a broken power law with a break at $39.7^{+20.3}_{-9.1}\,M_\odot$, or a power law with a Gaussian feature peaking at $33.1^{+4.0}_{-5.6}\,M_\odot$ (90\% credible interval). While the primary mass distribution must extend to $\sim65\,M_\odot$ or beyond, only $2.9^{+3.5}_{1.7}\%$ of systems have primary masses greater than $45\,M_\odot$. Second, we find that a fraction of BBH systems have component spins misaligned with the orbital angular momentum, giving rise to precession of the orbital plane. Moreover, 12% to 44% of BBH systems have spins tilted by more than $90^\circ$, giving rise to a negative effective inspiral spin parameter $χ_\mathrm{eff}$. Under the assumption that such systems can only be formed by dynamical interactions, we infer that between 25% and 93% of BBH with non-vanishing $|χ_\mathrm{eff}| > 0.01$ are dynamically assembled. Third, we estimate merger rates, finding $\mathcal{R}_\text{BBH} = 23.9^{+14.3}_{8.6}$ Gpc$^{-3}$ yr$^{-1}$ for BBH and $\mathcal{R}_\text{BNS}= 320^{+490}_{-240}$ Gpc$^{-3}$ yr$^{-1}$ for binary neutron stars. We find that the BBH rate likely increases with redshift ($85\%$ credibility), but not faster than the star-formation rate ($86\%$ credibility). Additionally, we examine recent exceptional events in the context of our population models, finding that the asymmetric masses of GW190412 and the high component masses of GW190521 are consistent with our models, but the low secondary mass of GW190814 makes it an outlier.
△ Less
Submitted 25 February, 2021; v1 submitted 27 October, 2020;
originally announced October 2020.
-
Tests of General Relativity with Binary Black Holes from the second LIGO-Virgo Gravitational-Wave Transient Catalog
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
R. Abbott,
T. D. Abbott,
S. Abraham,
F. Acernese,
K. Ackley,
A. Adams,
C. Adams,
R. X. Adhikari,
V. B. Adya,
C. Affeldt,
M. Agathos,
K. Agatsuma,
N. Aggarwal,
O. D. Aguiar,
L. Aiello,
A. Ain,
P. Ajith,
G. Allen,
A. Allocca,
P. A. Altin,
A. Amato,
S. Anand,
A. Ananyeva
, et al. (1322 additional authors not shown)
Abstract:
Gravitational waves enable tests of general relativity in the highly dynamical and strong-field regime. Using events detected by LIGO-Virgo up to 1 October 2019, we evaluate the consistency of the data with predictions from the theory. We first establish that residuals from the best-fit waveform are consistent with detector noise, and that the low- and high-frequency parts of the signals are in ag…
▽ More
Gravitational waves enable tests of general relativity in the highly dynamical and strong-field regime. Using events detected by LIGO-Virgo up to 1 October 2019, we evaluate the consistency of the data with predictions from the theory. We first establish that residuals from the best-fit waveform are consistent with detector noise, and that the low- and high-frequency parts of the signals are in agreement. We then consider parametrized modifications to the waveform by varying post-Newtonian and phenomenological coefficients, improving past constraints by factors of ${\sim}2$; we also find consistency with Kerr black holes when we specifically target signatures of the spin-induced quadrupole moment. Looking for gravitational-wave dispersion, we tighten constraints on Lorentz-violating coefficients by a factor of ${\sim}2.6$ and bound the mass of the graviton to $m_g \leq 1.76 \times 10^{-23} \mathrm{eV}/c^2$ with 90% credibility. We also analyze the properties of the merger remnants by measuring ringdown frequencies and damping times, constraining fractional deviations away from the Kerr frequency to $δ\hat{f}_{220} = 0.03^{+0.38}_{-0.35}$ for the fundamental quadrupolar mode, and $δ\hat{f}_{221} = 0.04^{+0.27}_{-0.32}$ for the first overtone; additionally, we find no evidence for postmerger echoes. Finally, we determine that our data are consistent with tensorial polarizations through a template-independent method. When possible, we assess the validity of general relativity based on collections of events analyzed jointly. We find no evidence for new physics beyond general relativity, for black hole mimickers, or for any unaccounted systematics.
△ Less
Submitted 16 June, 2021; v1 submitted 27 October, 2020;
originally announced October 2020.
-
GWTC-2: Compact Binary Coalescences Observed by LIGO and Virgo During the First Half of the Third Observing Run
Authors:
R. Abbott,
T. D. Abbott,
S. Abraham,
F. Acernese,
K. Ackley,
A. Adams,
C. Adams,
R. X. Adhikari,
V. B. Adya,
C. Affeldt,
M. Agathos,
K. Agatsuma,
N. Aggarwal,
O. D. Aguiar,
L. Aiello,
A. Ain,
P. Ajith,
S. Akcay,
G. Allen,
A. Allocca,
P. A. Altin,
A. Amato,
S. Anand,
A. Ananyeva,
S. B. Anderson
, et al. (1327 additional authors not shown)
Abstract:
We report on gravitational wave discoveries from compact binary coalescences detected by Advanced LIGO and Advanced Virgo in the first half of the third observing run (O3a) between 1 April 2019 15:00 UTC and 1 October 2019 15:00. By imposing a false-alarm-rate threshold of two per year in each of the four search pipelines that constitute our search, we present 39 candidate gravitational wave event…
▽ More
We report on gravitational wave discoveries from compact binary coalescences detected by Advanced LIGO and Advanced Virgo in the first half of the third observing run (O3a) between 1 April 2019 15:00 UTC and 1 October 2019 15:00. By imposing a false-alarm-rate threshold of two per year in each of the four search pipelines that constitute our search, we present 39 candidate gravitational wave events. At this threshold, we expect a contamination fraction of less than 10%. Of these, 26 candidate events were reported previously in near real-time through GCN Notices and Circulars; 13 are reported here for the first time. The catalog contains events whose sources are black hole binary mergers up to a redshift of ~0.8, as well as events whose components could not be unambiguously identified as black holes or neutron stars. For the latter group, we are unable to determine the nature based on estimates of the component masses and spins from gravitational wave data alone. The range of candidate events which are unambiguously identified as binary black holes (both objects $\geq 3~M_\odot$) is increased compared to GWTC-1, with total masses from $\sim 14~M_\odot$ for GW190924_021846 to $\sim 150~M_\odot$ for GW190521. For the first time, this catalog includes binary systems with significantly asymmetric mass ratios, which had not been observed in data taken before April 2019. We also find that 11 of the 39 events detected since April 2019 have positive effective inspiral spins under our default prior (at 90% credibility), while none exhibit negative effective inspiral spin. Given the increased sensitivity of Advanced LIGO and Advanced Virgo, the detection of 39 candidate events in ~26 weeks of data (~1.5 per week) is consistent with GWTC-1.
△ Less
Submitted 8 March, 2021; v1 submitted 27 October, 2020;
originally announced October 2020.
-
Sudden collapse of magnetic order in oxygen deficient nickelate films
Authors:
Jiarui Li,
Robert J. Green,
Zhen Zhang,
Ronny Sutarto,
Jerzy T. Sadowski,
Zhihai Zhu,
Grace Zhang,
Da Zhou,
Yifei Sun,
Feizhou He,
Shriram Ramanathan,
Riccardo Comin
Abstract:
Oxygen vacancies play a crucial role in the control of the electronic, magnetic, ionic, and transport properties of functional oxide perovskites. Rare earth nickelates (RENiO$_{3-x}$) have emerged over the years as a rich platform to study the interplay between the lattice, the electronic structure, and ordered magnetism. In this study, we investigate the evolution of the electronic and magnetic s…
▽ More
Oxygen vacancies play a crucial role in the control of the electronic, magnetic, ionic, and transport properties of functional oxide perovskites. Rare earth nickelates (RENiO$_{3-x}$) have emerged over the years as a rich platform to study the interplay between the lattice, the electronic structure, and ordered magnetism. In this study, we investigate the evolution of the electronic and magnetic structure in thin films of RENiO$_{3-x}$, using a combination of X-ray absorption spectroscopy and imaging, resonant X-ray scattering, and extended multiplet ligand field theory modeling. We find that oxygen vacancies modify the electronic configuration within the Ni-O orbital manifolds, leading to a dramatic evolution of long-range electronic transport pathways despite the absence of nanoscale phase separation. Remarkably, magnetism is robust to substantial levels of carrier doping, and only a moderate weakening of the $(1/4, 1/4, 1/4)_{pc}$ antiferromagnetic order parameter is observed, whereas the magnetic transition temperature is largely unchanged. Only at a certain point long-range magnetism is abruptly erased without an accompanying structural transition. We propose the progressive disruption of the 3D magnetic superexchange pathways upon introduction of point defects as the mechanism behind the sudden collapse of magnetic order in oxygen-deficient nickelates. Our work demonstrates that, unlike most other oxides, ordered magnetism in RENiO$_{3-x}$ is mostly insensitive to carrier doping. The sudden collapse of ordered magnetism upon oxygen removal may provide a new mechanism for solid-state magneto-ionic switching and new applications in antiferromagnetic spintronics.
△ Less
Submitted 16 October, 2020;
originally announced October 2020.
-
Identifying when Precession can be Measured in Gravitational Waveforms
Authors:
Rhys Green,
Charlie Hoy,
Stephen Fairhurst,
Mark Hannam,
Francesco Pannarale,
Cory Thomas
Abstract:
In binary-black-hole systems where the black-hole spins are misaligned with the orbital angular momentum, precession effects leave characteristic modulations in the emitted gravitational waveform. Here, we investigate where in the parameter space we will be able to accurately identify precession, for likely observations over coming LIGO-Virgo-KAGRA observing runs. Despite the large number of param…
▽ More
In binary-black-hole systems where the black-hole spins are misaligned with the orbital angular momentum, precession effects leave characteristic modulations in the emitted gravitational waveform. Here, we investigate where in the parameter space we will be able to accurately identify precession, for likely observations over coming LIGO-Virgo-KAGRA observing runs. Despite the large number of parameters that characterise a precessing binary, we perform a large scale systematic study to identify the impact of each source parameter on the measurement of precession. We simulate a fiducial binary at moderate mass-ratio, signal-to-noise ratio (SNR), and spins, such that precession will be clearly identifiable, then successively vary each parameter while holding the remaining parameters fixed. As expected, evidence for precession increases with signal-to noise-ratio (SNR), higher in-plane spins, more unequal component masses, and higher inclination, but our study provides a quantitative illustration of each of these effects, and informs our intuition on parameter dependencies that have not yet been studied in detail, for example, the effect of varying the relative strength of the two polarisations, the total mass, and the aligned-spin components. We also measure the "precession SNR" $ρ_p$, which was introduced in Refs[1,2] to quantify the signal power associated with precession. By comparing $ρ_p$ with both Bayes factors and the recovered posterior distributions, we find it is a reliable metric for measurability that accurately predicts when the detected signal contains evidence for precession.
△ Less
Submitted 12 May, 2021; v1 submitted 8 October, 2020;
originally announced October 2020.
-
Properties and astrophysical implications of the 150 Msun binary black hole merger GW190521
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
R. Abbott,
T. D. Abbott,
S. Abraham,
F. Acernese,
K. Ackley,
C. Adams,
R. X. Adhikari,
V. B. Adya,
C. Affeldt,
M. Agathos,
K. Agatsuma,
N. Aggarwal,
O. D. Aguiar,
A. Aich,
L. Aiello,
A. Ain,
P. Ajith,
S. Akcay,
G. Allen,
A. Allocca,
P. A. Altin,
A. Amato,
S. Anand
, et al. (1233 additional authors not shown)
Abstract:
The gravitational-wave signal GW190521 is consistent with a binary black hole merger source at redshift 0.8 with unusually high component masses, $85^{+21}_{-14}\,M_{\odot}$ and $66^{+17}_{-18}\,M_{\odot}$, compared to previously reported events, and shows mild evidence for spin-induced orbital precession. The primary falls in the mass gap predicted by (pulsational) pair-instability supernova theo…
▽ More
The gravitational-wave signal GW190521 is consistent with a binary black hole merger source at redshift 0.8 with unusually high component masses, $85^{+21}_{-14}\,M_{\odot}$ and $66^{+17}_{-18}\,M_{\odot}$, compared to previously reported events, and shows mild evidence for spin-induced orbital precession. The primary falls in the mass gap predicted by (pulsational) pair-instability supernova theory, in the approximate range $65 - 120\,M_{\odot}$. The probability that at least one of the black holes in GW190521 is in that range is 99.0%. The final mass of the merger $(142^{+28}_{-16}\,M_{\odot})$ classifies it as an intermediate-mass black hole. Under the assumption of a quasi-circular binary black hole coalescence, we detail the physical properties of GW190521's source binary and its post-merger remnant, including component masses and spin vectors. Three different waveform models, as well as direct comparison to numerical solutions of general relativity, yield consistent estimates of these properties. Tests of strong-field general relativity targeting the merger-ringdown stages of coalescence indicate consistency of the observed signal with theoretical predictions. We estimate the merger rate of similar systems to be $0.13^{+0.30}_{-0.11}\,{\rm Gpc}^{-3}\,\rm{yr}^{-1}$. We discuss the astrophysical implications of GW190521 for stellar collapse, and for the possible formation of black holes in the pair-instability mass gap through various channels: via (multiple) stellar coalescence, or via hierarchical merger of lower-mass black holes in star clusters or in active galactic nuclei. We find it to be unlikely that GW190521 is a strongly lensed signal of a lower-mass black hole binary merger. We also discuss more exotic possible sources for GW190521, including a highly eccentric black hole binary, or a primordial black hole binary.
△ Less
Submitted 2 September, 2020;
originally announced September 2020.
-
GW190521: A Binary Black Hole Merger with a Total Mass of $150 ~ M_{\odot}$
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
R. Abbott,
T. D. Abbott,
S. Abraham,
F. Acernese,
K. Ackley,
C. Adams,
R. X. Adhikari,
V. B. Adya,
C. Affeldt,
M. Agathos,
K. Agatsuma,
N. Aggarwal,
O. D. Aguiar,
A. Aich,
L. Aiello,
A. Ain,
P. Ajith,
S. Akcay,
G. Allen,
A. Allocca,
P. A. Altin,
A. Amato,
S. Anand
, et al. (1232 additional authors not shown)
Abstract:
On May 21, 2019 at 03:02:29 UTC Advanced LIGO and Advanced Virgo observed a short duration gravitational-wave signal, GW190521, with a three-detector network signal-to-noise ratio of 14.7, and an estimated false-alarm rate of 1 in 4900 yr using a search sensitive to generic transients. If GW190521 is from a quasicircular binary inspiral, then the detected signal is consistent with the merger of tw…
▽ More
On May 21, 2019 at 03:02:29 UTC Advanced LIGO and Advanced Virgo observed a short duration gravitational-wave signal, GW190521, with a three-detector network signal-to-noise ratio of 14.7, and an estimated false-alarm rate of 1 in 4900 yr using a search sensitive to generic transients. If GW190521 is from a quasicircular binary inspiral, then the detected signal is consistent with the merger of two black holes with masses of $85^{+21}_{-14} M_{\odot}$ and $66^{+17}_{-18} M_{\odot}$ (90 % credible intervals). We infer that the primary black hole mass lies within the gap produced by (pulsational) pair-instability supernova processes, and has only a 0.32 % probability of being below $65 M_{\odot}$. We calculate the mass of the remnant to be $142^{+28}_{-16} M_{\odot}$, which can be considered an intermediate mass black hole (IMBH). The luminosity distance of the source is $5.3^{+2.4}_{-2.6}$ Gpc, corresponding to a redshift of $0.82^{+0.28}_{-0.34}$. The inferred rate of mergers similar to GW190521 is $0.13^{+0.30}_{-0.11}\,\mathrm{Gpc}^{-3}\,\mathrm{yr}^{-1}$.
△ Less
Submitted 2 September, 2020;
originally announced September 2020.
-
Gravitational-wave surrogate models powered by artificial neural networks: The ANN-Sur for waveform generation
Authors:
Sebastian Khan,
Rhys Green
Abstract:
Inferring the properties of black holes and neutron stars is a key science goal of gravitational-wave (GW) astronomy. To extract as much information as possible from GW observations we must develop methods to reduce the cost of Bayesian inference. In this paper, we use artificial neural networks (ANNs) and the parallelisation power of graphics processing units (GPUs) to improve the surrogate model…
▽ More
Inferring the properties of black holes and neutron stars is a key science goal of gravitational-wave (GW) astronomy. To extract as much information as possible from GW observations we must develop methods to reduce the cost of Bayesian inference. In this paper, we use artificial neural networks (ANNs) and the parallelisation power of graphics processing units (GPUs) to improve the surrogate modelling method, which can produce accelerated versions of existing models. As a first application of our method, ANN-Sur, we build a time-domain surrogate model of the spin-aligned binary black hole (BBH) waveform model SEOBNRv4. We achieve median mismatches of 2e-5 and mismatches no worse than 2e-3. For a typical BBH waveform generated from 12 Hz with a total mass of $60 M_\odot$ the original SEOBNRv4 model takes 1812 ms. Existing bespoke code optimisations (SEOBNRv4opt) reduced this to 91.6 ms and the interpolation based, frequency-domain surrogate SEOBNRv4ROM can generate this waveform in 6.9 ms. Our ANN-Sur model, when run on a CPU takes 2.7 ms and just 0.4 ms when run on a GPU. ANN-Sur can also generate large batches of waveforms simultaneously. We find that batches of up to 10^4 waveforms can be evaluated on a GPU in just 163 ms, corresponding to a time per waveform of 0.016 ms. This method is a promising way to utilise the parallelisation power of GPUs to drastically increase the computational efficiency of Bayesian parameter estimation.
△ Less
Submitted 5 November, 2020; v1 submitted 29 August, 2020;
originally announced August 2020.
-
Non-normality and non-monotonic dynamics in complex reaction networks
Authors:
Zachary G. Nicolaou,
Takashi Nishikawa,
Schuyler B. Nicholson,
Jason R. Green,
Adilson E. Motter
Abstract:
Complex chemical reaction networks, which underlie many industrial and biological processes, often exhibit non-monotonic changes in chemical species concentrations, typically described using nonlinear models. Such non-monotonic dynamics are in principle possible even in linear models if the matrices defining the models are non-normal, as characterized by a necessarily non-orthogonal set of eigenve…
▽ More
Complex chemical reaction networks, which underlie many industrial and biological processes, often exhibit non-monotonic changes in chemical species concentrations, typically described using nonlinear models. Such non-monotonic dynamics are in principle possible even in linear models if the matrices defining the models are non-normal, as characterized by a necessarily non-orthogonal set of eigenvectors. However, the extent to which non-normality is responsible for non-monotonic behavior remains an open question. Here, using a master equation to model the reaction dynamics, we derive a general condition for observing non-monotonic dynamics of individual species, establishing that non-normality promotes non-monotonicity but is not a requirement for it. In contrast, we show that non-normality is a requirement for non-monotonic dynamics to be observed in the Rényi entropy. Using hydrogen combustion as an example application, we demonstrate that non-monotonic dynamics under experimental conditions are supported by a linear chain of connected components, in contrast with the dominance of a single giant component observed in typical random reaction networks. The exact linearity of the master equation enables development of rigorous theory and simulations for dynamical networks of unprecedented size (approaching $10^5$ dynamical variables, even for a network of only 20 reactions and involving less than 100 atoms). Our conclusions are expected to hold for other combustion processes, and the general theory we develop is applicable to all chemical reaction networks, including biological ones.
△ Less
Submitted 21 August, 2020;
originally announced August 2020.
-
On the $L^p$ boundedness of the Wave Operators for fourth order Schrödinger operators
Authors:
Michael Goldberg,
William R. Green
Abstract:
We consider the fourth order Schrödinger operator $H=Δ^2+V(x)$ in three dimensions with real-valued potential $V$. Let $H_0=Δ^2$, if $V$ decays sufficiently and there are no eigenvalues or resonances in the absolutely continuous spectrum of $H$ then the wave operators $W_{\pm}= s\,-\,\lim_{t\to \pm \infty} e^{itH}e^{-itH_0}$ extend to bounded operators on $L^p(\mathbb R^3)$ for all $1<p<\infty$.
We consider the fourth order Schrödinger operator $H=Δ^2+V(x)$ in three dimensions with real-valued potential $V$. Let $H_0=Δ^2$, if $V$ decays sufficiently and there are no eigenvalues or resonances in the absolutely continuous spectrum of $H$ then the wave operators $W_{\pm}= s\,-\,\lim_{t\to \pm \infty} e^{itH}e^{-itH_0}$ extend to bounded operators on $L^p(\mathbb R^3)$ for all $1<p<\infty$.
△ Less
Submitted 12 March, 2021; v1 submitted 17 August, 2020;
originally announced August 2020.
-
From Ji to Jaffe-Manohar orbital angular momentum in Lattice QCD using a direct derivative method
Authors:
M. Engelhardt,
J. R. Green,
N. Hasan,
S. Krieg,
S. Meinel,
J. Negele,
A. Pochinsky,
S. Syritsyn
Abstract:
A Lattice QCD approach to quark orbital angular momentum in the proton based on generalized transverse momentum-dependent parton distributions (GTMDs) is enhanced methodologically by incorporating a direct derivative technique. This improvement removes a significant numerical bias that had been seen to afflict results of a previous study. In particular, the value obtained for Ji quark orbital angu…
▽ More
A Lattice QCD approach to quark orbital angular momentum in the proton based on generalized transverse momentum-dependent parton distributions (GTMDs) is enhanced methodologically by incorporating a direct derivative technique. This improvement removes a significant numerical bias that had been seen to afflict results of a previous study. In particular, the value obtained for Ji quark orbital angular momentum is reconciled with the one obtained independently via Ji's sum rule, validating the GMTD approach. Since GTMDs simultaneously contain information about the quark impact parameter and transverse momentum, they permit a direct evaluation of the cross product of the latter. They are defined through proton matrix elements of a quark bilocal operator containing a Wilson line; the choice in Wilson line path allows one to continuously interpolate from Ji to Jaffe-Manohar quark orbital angular momentum. The latter is seen to be significantly enhanced in magnitude compared to Ji quark orbital angular momentum, confirming previous results.
△ Less
Submitted 9 August, 2020;
originally announced August 2020.
-
Complete parameter inference for GW150914 using deep learning
Authors:
Stephen R. Green,
Jonathan Gair
Abstract:
The LIGO and Virgo gravitational-wave observatories have detected many exciting events over the past five years. As the rate of detections grows with detector sensitivity, this poses a growing computational challenge for data analysis. With this in mind, in this work we apply deep learning techniques to perform fast likelihood-free Bayesian inference for gravitational waves. We train a neural-netw…
▽ More
The LIGO and Virgo gravitational-wave observatories have detected many exciting events over the past five years. As the rate of detections grows with detector sensitivity, this poses a growing computational challenge for data analysis. With this in mind, in this work we apply deep learning techniques to perform fast likelihood-free Bayesian inference for gravitational waves. We train a neural-network conditional density estimator to model posterior probability distributions over the full 15-dimensional space of binary black hole system parameters, given detector strain data from multiple detectors. We use the method of normalizing flows---specifically, a neural spline normalizing flow---which allows for rapid sampling and density estimation. Training the network is likelihood-free, requiring samples from the data generative process, but no likelihood evaluations. Through training, the network learns a global set of posteriors: it can generate thousands of independent posterior samples per second for any strain data consistent with the prior and detector noise characteristics used for training. By training with the detector noise power spectral density estimated at the time of GW150914, and conditioning on the event strain data, we use the neural network to generate accurate posterior samples consistent with analyses using conventional sampling techniques.
△ Less
Submitted 7 August, 2020;
originally announced August 2020.
-
Gravitational-wave constraints on the equatorial ellipticity of millisecond pulsars
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
R. Abbott,
T. D. Abbott,
S. Abraham,
F. Acernese,
K. Ackley,
A. Adams,
C. Adams,
R. X. Adhikari,
V. B. Adya,
C. Affeldt,
M. Agathos,
K. Agatsuma,
N. Aggarwal,
O. D. Aguiar,
L. Aiello,
A. Ain,
P. Ajith,
G. Allen,
A. Allocca,
P. A. Altin,
A. Amato,
S. Anand,
A. Ananyeva
, et al. (1311 additional authors not shown)
Abstract:
We present a search for continuous gravitational waves from five radio pulsars, comprising three recycled pulsars (PSR J0437-4715, PSR J0711-6830, and PSR J0737-3039A) and two young pulsars: the Crab pulsar (J0534+2200) and the Vela pulsar (J0835-4510). We use data from the third observing run of Advanced LIGO and Virgo combined with data from their first and second observing runs. For the first t…
▽ More
We present a search for continuous gravitational waves from five radio pulsars, comprising three recycled pulsars (PSR J0437-4715, PSR J0711-6830, and PSR J0737-3039A) and two young pulsars: the Crab pulsar (J0534+2200) and the Vela pulsar (J0835-4510). We use data from the third observing run of Advanced LIGO and Virgo combined with data from their first and second observing runs. For the first time we are able to match (for PSR J0437-4715) or surpass (for PSR J0711-6830) the indirect limits on gravitational-wave emission from recycled pulsars inferred from their observed spin-downs, and constrain their equatorial ellipticities to be less than $10^{-8}$. For each of the five pulsars, we perform targeted searches that assume a tight coupling between the gravitational-wave and electromagnetic signal phase evolution. We also present constraints on PSR J0711-6830, the Crab pulsar and the Vela pulsar from a search that relaxes this assumption, allowing the gravitational-wave signal to vary from the electromagnetic expectation within a narrow band of frequencies and frequency derivatives.
△ Less
Submitted 13 October, 2020; v1 submitted 28 July, 2020;
originally announced July 2020.
-
Time integrable weighted dispersive estimates for the fourth order Schrödinger equation in three dimensions
Authors:
Michael Goldberg,
William R. Green
Abstract:
We consider the fourth order Schrödinger operator $H=Δ^2+V$ and show that if there are no eigenvalues or resonances in the absolutely continuous spectrum of $H$ that the solution operator $e^{-itH}$ satisfies a large time integrable $|t|^{-\frac54}$ decay rate between weighted spaces. This bound improves what is possible for the free case in two directions; both better time decay and smaller spati…
▽ More
We consider the fourth order Schrödinger operator $H=Δ^2+V$ and show that if there are no eigenvalues or resonances in the absolutely continuous spectrum of $H$ that the solution operator $e^{-itH}$ satisfies a large time integrable $|t|^{-\frac54}$ decay rate between weighted spaces. This bound improves what is possible for the free case in two directions; both better time decay and smaller spatial weights. In the case of a mild resonance at zero energy, we derive the operator-valued expansion $e^{-itH}P_{ac}(H)=t^{-\frac34} A_0+t^{-\frac54}A_1$ where $A_0:L^1\to L^\infty$ is an operator of rank at most four and $A_1$ maps between polynomially weighted spaces.
△ Less
Submitted 2 June, 2021; v1 submitted 9 July, 2020;
originally announced July 2020.
-
Hadronic light-by-light contribution to $(g-2)_μ$ from lattice QCD with SU(3) flavor symmetry
Authors:
En-Hung Chao,
Antoine Gérardin,
Jeremy R. Green,
Renwick J. Hudspith,
Harvey B. Meyer
Abstract:
We perform a lattice QCD calculation of the hadronic light-by-light contribution to $(g-2)_μ$ at the SU(3) flavor-symmetric point $m_π=m_K\simeq 420\,$MeV. The representation used is based on coordinate-space perturbation theory, with all QED elements of the relevant Feynman diagrams implemented in continuum, infinite Euclidean space. As a consequence, the effect of using finite lattices to evalua…
▽ More
We perform a lattice QCD calculation of the hadronic light-by-light contribution to $(g-2)_μ$ at the SU(3) flavor-symmetric point $m_π=m_K\simeq 420\,$MeV. The representation used is based on coordinate-space perturbation theory, with all QED elements of the relevant Feynman diagrams implemented in continuum, infinite Euclidean space. As a consequence, the effect of using finite lattices to evaluate the QCD four-point function of the electromagnetic current is exponentially suppressed. Thanks to the SU(3)-flavor symmetry, only two topologies of diagrams contribute, the fully connected and the leading disconnected. We show the equivalence in the continuum limit of two methods of computing the connected contribution, and introduce a sparse-grid technique for computing the disconnected contribution. Thanks to our previous calculation of the pion transition form factor, we are able to correct for the residual finite-size effects and extend the tail of the integrand. We test our understanding of finite-size effects by using gauge ensembles differing only by their volume. After a continuum extrapolation based on four lattice spacings, we obtain $a_μ^{\rm hlbl} = (65.4\pm 4.9 \pm 6.6)\times 10^{-11}$, where the first error results from the uncertainties on the individual gauge ensembles and the second is the systematic error of the continuum extrapolation. Finally, we estimate how this value will change as the light-quark masses are lowered to their physical values.
△ Less
Submitted 13 July, 2020; v1 submitted 29 June, 2020;
originally announced June 2020.
-
3D Mobility Models and Analysis for UAVs
Authors:
Peter J. Smith,
Pawel A. Dmochowski,
Ikram Singh,
Richard Green,
Carl P. Dettmann,
Justin P. Coon
Abstract:
We present a flexible family of 3D mobility models suitable for unmanned aerial vehicles (UAV). Based on stochastic differential equations, the models offer a unique property of explicitly incorporating the mobility control mechanism and environmental perturbation, while enabling tractable steady state solutions for properties such as position and connectivity. Specifically, motivated by UAV fligh…
▽ More
We present a flexible family of 3D mobility models suitable for unmanned aerial vehicles (UAV). Based on stochastic differential equations, the models offer a unique property of explicitly incorporating the mobility control mechanism and environmental perturbation, while enabling tractable steady state solutions for properties such as position and connectivity. Specifically, motivated by UAV flight data, for a symmetric mobility model with an arbitrary control mechanism, we derive the steady state distribution of the distance from the target position. We provide closed form expressions for the special cases of the Ornstein-Uhlenbeck (OU) process and on-off control (OC). We extend the model to incorporate imperfect positioning and asymmetric control. For a practically relevant scenario of partial symmetry (such as in the x-y plane), we present steady state position results for the OU control. Building on these results, we derive UAV connectivity probability results based on a SNR criterion in a Rayleigh fading environment.
△ Less
Submitted 27 June, 2020;
originally announced June 2020.
-
Pōniuā'ena: A Luminous $z=7.5$ Quasar Hosting a 1.5 Billion Solar Mass Black Hole
Authors:
Jinyi Yang,
Feige Wang,
Xiaohui Fan,
Joseph F. Hennawi,
Frederick B. Davies,
Minghao Yue,
Eduardo Banados,
Xue-Bing Wu,
Bram Venemans,
Aaron J. Barth,
Fuyan Bian,
Konstantina Boutsia,
Roberto Decarli,
Emanuele Paolo Farina,
Richard Green,
Linhua Jiang,
Jiang-Tao Li,
Chiara Mazzucchelli,
Fabian Walter
Abstract:
We report the discovery of a luminous quasar, J1007+2115 at $z=7.515$ ("Pōniuā'ena"), from our wide-field reionization-era quasar survey. J1007+2115 is the second quasar now known at $z>7.5$, deep into the reionization epoch. The quasar is powered by a $(1.5\pm0.2)\times10^9$ $M_{\odot}$ supermassive black hole (SMBH), based on its broad MgII emission-line profile from Gemini and Keck near-IR spec…
▽ More
We report the discovery of a luminous quasar, J1007+2115 at $z=7.515$ ("Pōniuā'ena"), from our wide-field reionization-era quasar survey. J1007+2115 is the second quasar now known at $z>7.5$, deep into the reionization epoch. The quasar is powered by a $(1.5\pm0.2)\times10^9$ $M_{\odot}$ supermassive black hole (SMBH), based on its broad MgII emission-line profile from Gemini and Keck near-IR spectroscopy. The SMBH in J1007+2115 is twice as massive as that in quasar J1342+0928 at $z=7.54$, the current quasar redshift record holder. The existence of such a massive SMBH just 700 million years after the Big Bang significantly challenges models of the earliest SMBH growth. Model assumptions of Eddington-limited accretion and a radiative efficiency of 0.1 require a seed black hole of $\gtrsim 10^{4}$ $M_{\odot}$ at $z=30$. This requirement suggests either a massive black hole seed as a result of direct collapse or earlier periods of rapid black hole growth with hyper-Eddington accretion and/or a low radiative efficiency. We measure the damping wing signature imprinted by neutral hydrogen absorption in the intergalactic medium (IGM) on J1007+2115's Ly$α$ line profile, and find that it is weaker than that of J1342+0928 and two other $z\gtrsim7$ quasars. We estimate an IGM volume-averaged neutral fraction $\langle x\rm_{HI}\rangle=0.39^{+0.22}_{-0.13}$. This range of values suggests a patchy reionization history toward different IGM sightlines. We detect the 158 $μ$m [C II] emission line in J1007+2115 with ALMA; this line centroid yields a systemic redshift of $z=7.5149\pm0.0004$ and indicates a star formation rate of $\sim210$ $M_{\odot}$ yr$^{-1}$ in its host galaxy.
△ Less
Submitted 23 June, 2020;
originally announced June 2020.
-
Distance Correlation Sure Independence Screening for Accelerated Feature Selection in Parkinson's Disease Vocal Data
Authors:
Dan Schellhas,
Bishal Neupane,
Deepak Thammineni,
Bhargav Kanumuri,
Robert C. Green II
Abstract:
With the abundance of machine learning methods available and the temptation of using them all in an ensemble method, having a model-agnostic method of feature selection is incredibly alluring. Principal component analysis was developed in 1901 and has been a strong contender in this role since, but in the end is an unsupervised method. It offers no guarantee that the features that are selected hav…
▽ More
With the abundance of machine learning methods available and the temptation of using them all in an ensemble method, having a model-agnostic method of feature selection is incredibly alluring. Principal component analysis was developed in 1901 and has been a strong contender in this role since, but in the end is an unsupervised method. It offers no guarantee that the features that are selected have good predictive power because it does not know what is being predicted. To this end, Peng et al. developed the minimum redundancy-maximum relevance (mRMR) method in 2005. It uses the mutual information not only between predictors but also includes the mutual information with the response in its calculation. Estimating mutual information and entropy tend to be expensive and problematic endeavors, which leads to excessive processing times even for dataset that is approximately 750 by 750 in a Leave-One-Subject-Out jackknife situation. To remedy this, we use a method from 2012 called Distance Correlation Sure Independence Screening (DC-SIS) which uses the distance correlation measure of Székely et al. to select features that have the greatest dependence with the response. We show that this method produces statistically indistinguishable results to the mRMR selection method on Parkinson's Disease vocal diagnosis data 90 times faster.
△ Less
Submitted 23 June, 2020;
originally announced June 2020.
-
GW190814: Gravitational Waves from the Coalescence of a 23 M$_\odot$ Black Hole with a 2.6 M$_\odot$ Compact Object
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
R. Abbott,
T. D. Abbott,
S. Abraham,
F. Acernese,
K. Ackley,
C. Adams,
R. X. Adhikari,
V. B. Adya,
C. Affeldt,
M. Agathos,
K. Agatsuma,
N. Aggarwal,
O. D. Aguiar,
A. Aich,
L. Aiello,
A. Ain,
P. Ajith,
S. Akcay,
G. Allen,
A. Allocca,
P. A. Altin,
A. Amato,
S. Anand
, et al. (1232 additional authors not shown)
Abstract:
We report the observation of a compact binary coalescence involving a 22.2 - 24.3 $M_{\odot}$ black hole and a compact object with a mass of 2.50 - 2.67 $M_{\odot}$ (all measurements quoted at the 90$\%$ credible level). The gravitational-wave signal, GW190814, was observed during LIGO's and Virgo's third observing run on August 14, 2019 at 21:10:39 UTC and has a signal-to-noise ratio of 25 in the…
▽ More
We report the observation of a compact binary coalescence involving a 22.2 - 24.3 $M_{\odot}$ black hole and a compact object with a mass of 2.50 - 2.67 $M_{\odot}$ (all measurements quoted at the 90$\%$ credible level). The gravitational-wave signal, GW190814, was observed during LIGO's and Virgo's third observing run on August 14, 2019 at 21:10:39 UTC and has a signal-to-noise ratio of 25 in the three-detector network. The source was localized to 18.5 deg$^2$ at a distance of $241^{+41}_{-45}$ Mpc; no electromagnetic counterpart has been confirmed to date. The source has the most unequal mass ratio yet measured with gravitational waves, $0.112^{+0.008}_{-0.009}$, and its secondary component is either the lightest black hole or the heaviest neutron star ever discovered in a double compact-object system. The dimensionless spin of the primary black hole is tightly constrained to $\leq 0.07$. Tests of general relativity reveal no measurable deviations from the theory, and its prediction of higher-multipole emission is confirmed at high confidence. We estimate a merger rate density of 1-23 Gpc$^{-3}$ yr$^{-1}$ for the new class of binary coalescence sources that GW190814 represents. Astrophysical models predict that binaries with mass ratios similar to this event can form through several channels, but are unlikely to have formed in globular clusters. However, the combination of mass ratio, component masses, and the inferred merger rate for this event challenges all current models for the formation and mass distribution of compact-object binaries.
△ Less
Submitted 22 June, 2020;
originally announced June 2020.
-
The anomalous magnetic moment of the muon in the Standard Model
Authors:
T. Aoyama,
N. Asmussen,
M. Benayoun,
J. Bijnens,
T. Blum,
M. Bruno,
I. Caprini,
C. M. Carloni Calame,
M. Cè,
G. Colangelo,
F. Curciarello,
H. Czyż,
I. Danilkin,
M. Davier,
C. T. H. Davies,
M. Della Morte,
S. I. Eidelman,
A. X. El-Khadra,
A. Gérardin,
D. Giusti,
M. Golterman,
Steven Gottlieb,
V. Gülpers,
F. Hagelstein,
M. Hayakawa
, et al. (107 additional authors not shown)
Abstract:
We review the present status of the Standard Model calculation of the anomalous magnetic moment of the muon. This is performed in a perturbative expansion in the fine-structure constant $α$ and is broken down into pure QED, electroweak, and hadronic contributions. The pure QED contribution is by far the largest and has been evaluated up to and including $\mathcal{O}(α^5)$ with negligible numerical…
▽ More
We review the present status of the Standard Model calculation of the anomalous magnetic moment of the muon. This is performed in a perturbative expansion in the fine-structure constant $α$ and is broken down into pure QED, electroweak, and hadronic contributions. The pure QED contribution is by far the largest and has been evaluated up to and including $\mathcal{O}(α^5)$ with negligible numerical uncertainty. The electroweak contribution is suppressed by $(m_μ/M_W)^2$ and only shows up at the level of the seventh significant digit. It has been evaluated up to two loops and is known to better than one percent. Hadronic contributions are the most difficult to calculate and are responsible for almost all of the theoretical uncertainty. The leading hadronic contribution appears at $\mathcal{O}(α^2)$ and is due to hadronic vacuum polarization, whereas at $\mathcal{O}(α^3)$ the hadronic light-by-light scattering contribution appears. Given the low characteristic scale of this observable, these contributions have to be calculated with nonperturbative methods, in particular, dispersion relations and the lattice approach to QCD. The largest part of this review is dedicated to a detailed account of recent efforts to improve the calculation of these two contributions with either a data-driven, dispersive approach, or a first-principle, lattice-QCD approach. The final result reads $a_μ^\text{SM}=116\,591\,810(43)\times 10^{-11}$ and is smaller than the Brookhaven measurement by 3.7$σ$. The experimental uncertainty will soon be reduced by up to a factor four by the new experiment currently running at Fermilab, and also by the future J-PARC experiment. This and the prospects to further reduce the theoretical uncertainty in the near future-which are also discussed here-make this quantity one of the most promising places to look for evidence of new physics.
△ Less
Submitted 13 November, 2020; v1 submitted 8 June, 2020;
originally announced June 2020.
-
Bayesian inference for compact binary coalescences with BILBY: Validation and application to the first LIGO--Virgo gravitational-wave transient catalogue
Authors:
I. M. Romero-Shaw,
C. Talbot,
S. Biscoveanu,
V. D'Emilio,
G. Ashton,
C. P. L. Berry,
S. Coughlin,
S. Galaudage,
C. Hoy,
M. Huebner,
K. S. Phukon,
M. Pitkin,
M. Rizzo,
N. Sarin,
R. Smith,
S. Stevenson,
A. Vajpeyi,
M. Arene,
K. Athar,
S. Banagiri,
N. Bose,
M. Carney,
K. Chatziioannou,
J. A. Clark,
M. Colleoni
, et al. (34 additional authors not shown)
Abstract:
Gravitational waves provide a unique tool for observational astronomy. While the first LIGO--Virgo catalogue of gravitational-wave transients (GWTC-1) contains eleven signals from black hole and neutron star binaries, the number of observations is increasing rapidly as detector sensitivity improves. To extract information from the observed signals, it is imperative to have fast, flexible, and scal…
▽ More
Gravitational waves provide a unique tool for observational astronomy. While the first LIGO--Virgo catalogue of gravitational-wave transients (GWTC-1) contains eleven signals from black hole and neutron star binaries, the number of observations is increasing rapidly as detector sensitivity improves. To extract information from the observed signals, it is imperative to have fast, flexible, and scalable inference techniques. In a previous paper, we introduced BILBY: a modular and user-friendly Bayesian inference library adapted to address the needs of gravitational-wave inference. In this work, we demonstrate that BILBY produces reliable results for simulated gravitational-wave signals from compact binary mergers, and verify that it accurately reproduces results reported for the eleven GWTC-1 signals. Additionally, we provide configuration and output files for all analyses to allow for easy reproduction, modification, and future use. This work establishes that BILBY is primed and ready to analyse the rapidly growing population of compact binary coalescence gravitational-wave signals.
△ Less
Submitted 7 June, 2021; v1 submitted 1 June, 2020;
originally announced June 2020.
-
CIV emission line properties and uncertainties in black hole mass estimates of z ~ 3.5 quasars
Authors:
Wenwen Zuo,
Xue-Bing Wu,
Xiaohui Fan,
Richard Green,
Weiming Yi,
Andreas Schulze,
Ran Wang,
Fuyan Bian
Abstract:
Using a high luminosity ($L_{\rm bol} \sim 10^{47.5} - 10^{48.3}$ erg s$^{-1}$), high redshift ($3.2 < z < 3.8$) quasar sample of 19 quasars with optical and near-infrared spectroscopy, we investigate the reliability of the CIV-based black hole mass estimates ($M_{\rm BH}$). The median logarithm of the CIV- and H$β$-based $M_{\rm BH}$ ratios is 0.110 dex with the scatter of 0.647 dex. The CIV-to-H…
▽ More
Using a high luminosity ($L_{\rm bol} \sim 10^{47.5} - 10^{48.3}$ erg s$^{-1}$), high redshift ($3.2 < z < 3.8$) quasar sample of 19 quasars with optical and near-infrared spectroscopy, we investigate the reliability of the CIV-based black hole mass estimates ($M_{\rm BH}$). The median logarithm of the CIV- and H$β$-based $M_{\rm BH}$ ratios is 0.110 dex with the scatter of 0.647 dex. The CIV-to-H$β$ BH mass differences are significantly correlated with the CIV FWHMs, blueshifts and asymmetries. Corrections of the CIV FWHM using the blueshift and asymmetry reduce the scatter of the mass differences by $\sim$ 0.04-0.2 dex. Quasars in our sample accrete at the Eddington ratio $R_{\rm EDD}>0.3$ and cover a considerable range of blueshifts, with 18/19 of the quasars showing CIV blueshifts (with the median value of 1126 km s$^{-1}$) and 14/19 of the quasars showing CIV blueshifts larger than 500 km s$^{-1}$. It suggests that not all quasars with high Eddington ratios show large blueshifts. The Baldwin effect between the CIV rest-frame equivalent width (REW) and the continuum luminosity at 1350 $\overset{\lower.5em\circ}{\mathrm{A}}$ is not seen, likely due to the limited luminosity range of our sample. We find a lack of flux in the red wing of the composite spectrum with larger CIV blueshift, and detect a higher ratio of [OIII] quasars with REW$_{\rm [OIII]}>5$ $\overset{\lower.5em\circ}{\mathrm{A}}$ in the subsample with lower CIV blueshift. It is more likely that they are caused by the combination of the Eddington ratio and the orientation effect.
△ Less
Submitted 8 May, 2020;
originally announced May 2020.
-
GW190412: Observation of a Binary-Black-Hole Coalescence with Asymmetric Masses
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
R. Abbott,
T. D. Abbott,
S. Abraham,
F. Acernese,
K. Ackley,
C. Adams,
R. X. Adhikari,
V. B. Adya,
C. Affeldt,
M. Agathos,
K. Agatsuma,
N. Aggarwal,
O. D. Aguiar,
A. Aich,
L. Aiello,
A. Ain,
P. Ajith,
S. Akcay,
G. Allen,
A. Allocca,
P. A. Altin,
A. Amato,
S. Anand
, et al. (1232 additional authors not shown)
Abstract:
We report the observation of gravitational waves from a binary-black-hole coalescence during the first two weeks of LIGO's and Virgo's third observing run. The signal was recorded on April 12, 2019 at 05:30:44 UTC with a network signal-to-noise ratio of 19. The binary is different from observations during the first two observing runs most notably due to its asymmetric masses: a ~30 solar mass blac…
▽ More
We report the observation of gravitational waves from a binary-black-hole coalescence during the first two weeks of LIGO's and Virgo's third observing run. The signal was recorded on April 12, 2019 at 05:30:44 UTC with a network signal-to-noise ratio of 19. The binary is different from observations during the first two observing runs most notably due to its asymmetric masses: a ~30 solar mass black hole merged with a ~8 solar mass black hole companion. The more massive black hole rotated with a dimensionless spin magnitude between 0.22 and 0.60 (90% probability). Asymmetric systems are predicted to emit gravitational waves with stronger contributions from higher multipoles, and indeed we find strong evidence for gravitational radiation beyond the leading quadrupolar order in the observed signal. A suite of tests performed on GW190412 indicates consistency with Einstein's general theory of relativity. While the mass ratio of this system differs from all previous detections, we show that it is consistent with the population model of stellar binary black holes inferred from the first two observing runs.
△ Less
Submitted 24 August, 2020; v1 submitted 17 April, 2020;
originally announced April 2020.
-
Improvement, generalization, and scheme conversion of Wilson-line operators on the lattice in the auxiliary field approach
Authors:
Jeremy R. Green,
Karl Jansen,
Fernanda Steffens
Abstract:
Nonlocal quark bilinear operators connected by link paths are used for studying parton distribution functions (PDFs) and transverse momentum-dependent PDFs of hadrons using lattice QCD. The nonlocality makes it difficult to understand the renormalization and improvement of these operators using standard methods. In previous work, we showed that by introducing an auxiliary field on the lattice, one…
▽ More
Nonlocal quark bilinear operators connected by link paths are used for studying parton distribution functions (PDFs) and transverse momentum-dependent PDFs of hadrons using lattice QCD. The nonlocality makes it difficult to understand the renormalization and improvement of these operators using standard methods. In previous work, we showed that by introducing an auxiliary field on the lattice, one can understand an on-axis Wilson-line operator as the product of two local operators in an extended theory. In this paper, we provide details about the calculation in perturbation theory of the factor for conversion from our lattice-suitable renormalization scheme to the MS-bar scheme. Extending our work, we study Symanzik improvement of the extended theory to understand the pattern of discretization effects linear in the lattice spacing, $a$, which are present even if the lattice fermion action exactly preserves chiral symmetry. This provides a prospect for an eventual $O(a)$ improvement of lattice calculations of PDFs. We also generalize our approach to apply to Wilson lines along lattice diagonals and to piecewise-straight link paths.
△ Less
Submitted 16 April, 2020; v1 submitted 21 February, 2020;
originally announced February 2020.
-
Prescriptions for Correcting Ultraviolet-Based Redshifts for Luminous Quasars at High Redshift
Authors:
Cooper Dix,
Ohad Shemmer,
Michael S. Brotherton,
Richard F. Green,
Michelle Mason,
Adam D. Myers
Abstract:
High-redshift quasars typically have their redshift determined from rest-frame ultraviolet (UV) emission lines. However, these lines, and more specifically the prominent C IV $λ1549$ emission line, are typically blueshifted yielding highly uncertain redshift estimates compared to redshifts determined from rest-frame optical emission lines. We present near-infrared spectroscopy of 18 luminous quasa…
▽ More
High-redshift quasars typically have their redshift determined from rest-frame ultraviolet (UV) emission lines. However, these lines, and more specifically the prominent C IV $λ1549$ emission line, are typically blueshifted yielding highly uncertain redshift estimates compared to redshifts determined from rest-frame optical emission lines. We present near-infrared spectroscopy of 18 luminous quasars at $2.15 < z < 3.70$ that allows us to obtain reliable systemic redshifts for these sources. Together with near-infrared spectroscopy of an archival sample of 44 quasars with comparable luminosities and redshifts, we provide prescriptions for correcting UV-based redshifts. Our prescriptions reduce velocity offsets with respect to the systemic redshifts by $\sim140$ km s$^{-1}$ and reduce the uncertainty on the UV-based redshift by $\sim25\%$ with respect to the best method currently used for determining such values. We also find that the redshifts determined from the Sloan Digital Sky Survey Pipeline for our sources suffer from significant uncertainties, which cannot be easily mitigated. We discuss the potential of our prescriptions to improve UV-based redshift corrections given a much larger sample of high redshift quasars with near-infrared spectra.
△ Less
Submitted 19 February, 2020;
originally announced February 2020.
-
Gravitational-wave parameter estimation with autoregressive neural network flows
Authors:
Stephen R. Green,
Christine Simpson,
Jonathan Gair
Abstract:
We introduce the use of autoregressive normalizing flows for rapid likelihood-free inference of binary black hole system parameters from gravitational-wave data with deep neural networks. A normalizing flow is an invertible mapping on a sample space that can be used to induce a transformation from a simple probability distribution to a more complex one: if the simple distribution can be rapidly sa…
▽ More
We introduce the use of autoregressive normalizing flows for rapid likelihood-free inference of binary black hole system parameters from gravitational-wave data with deep neural networks. A normalizing flow is an invertible mapping on a sample space that can be used to induce a transformation from a simple probability distribution to a more complex one: if the simple distribution can be rapidly sampled and its density evaluated, then so can the complex distribution. Our first application to gravitational waves uses an autoregressive flow, conditioned on detector strain data, to map a multivariate standard normal distribution into the posterior distribution over system parameters. We train the model on artificial strain data consisting of IMRPhenomPv2 waveforms drawn from a five-parameter $(m_1, m_2, φ_0, t_c, d_L)$ prior and stationary Gaussian noise realizations with a fixed power spectral density. This gives performance comparable to current best deep-learning approaches to gravitational-wave parameter estimation. We then build a more powerful latent variable model by incorporating autoregressive flows within the variational autoencoder framework. This model has performance comparable to Markov chain Monte Carlo and, in particular, successfully models the multimodal $φ_0$ posterior. Finally, we train the autoregressive latent variable model on an expanded parameter space, including also aligned spins $(χ_{1z}, χ_{2z})$ and binary inclination $θ_{JN}$, and show that all parameters and degeneracies are well-recovered. In all cases, sampling is extremely fast, requiring less than two seconds to draw $10^4$ posterior samples.
△ Less
Submitted 18 February, 2020;
originally announced February 2020.
-
Landau Level Phases in Bilayer Graphene under Pressure at Charge Neutrality
Authors:
Brett R. Green,
Jorge O. Sofo
Abstract:
Bilayer graphene in a magnetic field hosts a variety of ordered phases built from eight Landau levels close in energy to the neutrality point. These levels are characterized by orbital $n=0,1$, valley $ξ=+,-$ and spin $σ=\uparrow,\downarrow$; their relative energies depend strongly on the Coulomb interaction, magnetic field, and interlayer bias. We treat interactions at the Hartree-Fock level, inc…
▽ More
Bilayer graphene in a magnetic field hosts a variety of ordered phases built from eight Landau levels close in energy to the neutrality point. These levels are characterized by orbital $n=0,1$, valley $ξ=+,-$ and spin $σ=\uparrow,\downarrow$; their relative energies depend strongly on the Coulomb interaction, magnetic field, and interlayer bias. We treat interactions at the Hartree-Fock level, including the effects of metallic gates, layer separation, spatial extent of the $p_z$ orbitals, all Slonczewski-Weiss-McClure tight-binding parameters, and pressure. We obtain the ground state as function of the applied magnetic field, bias, and pressure. The gates, layer separation and extent of the $p_z$ orbitals weaken the Coulomb interaction at different length scales; these effects distort the phase diagram but do not change its topology. However, previously-predicted continuous transitions become discontinuous when all tight-binding parameters are included nonperturbatively. We find that pressure increases the importance of the noninteracting scale with respect to the Coulomb energy, which drives phase transitions to occur at lower fields. This brings two orbitally polarized states not yet predicted or observed into the experimentally accessible region of the phase diagram, in addition to previously-identified valley-, spin-, and partially orbitally polarized states.
△ Less
Submitted 13 February, 2020;
originally announced February 2020.
-
Specht module branching rules for wreath products of symmetric groups
Authors:
Reuben Green
Abstract:
We review a class of modules for the wreath product S(m) wr S(n) of two symmetric groups which are analogous to the Specht modules of the symmetric group, and prove a pair of branching rules for this family of modules. These branching rules describe the behaviour of these wreath product Specht modules under restriction to the wreath products S(m-1) wr S(n) and S(m) wr S(n-1). In particular, we see…
▽ More
We review a class of modules for the wreath product S(m) wr S(n) of two symmetric groups which are analogous to the Specht modules of the symmetric group, and prove a pair of branching rules for this family of modules. These branching rules describe the behaviour of these wreath product Specht modules under restriction to the wreath products S(m-1) wr S(n) and S(m) wr S(n-1). In particular, we see that these restrictions of wreath product Specht modules have Specht module filtrations, and we obtain combinatorial interpretations of the multiplicities in these filtrations.
△ Less
Submitted 19 January, 2020;
originally announced January 2020.
-
Ubiquitous impact of localised impurity states on the exchange coupling mechanism in magnetic topological insulators
Authors:
Thiago R. F. Peixoto,
Hendrik Bentmann,
Philipp Rüßmann,
Abdul-Vakhab Tcakaev,
Martin Winnerlein,
Steffen Schreyeck,
Sonja Schatz,
Raphael Crespo Vidal,
Fabian Stier,
Volodymyr Zabolotnyy,
Robert J. Green,
Chul Hee Min,
Celso I. Fornari,
Henriette Maaß,
Hari Babu Vasili,
Pierluigi Gargiani,
Manuel Valvidares,
Alessandro Barla,
Jens Buck,
Moritz Hoesch,
Florian Diekmann,
Sebastian Rohlf,
Matthias Kalläne,
Kai Rossnagel,
Charles Gould
, et al. (5 additional authors not shown)
Abstract:
Since the discovery of the quantum anomalous Hall effect in the magnetically doped topological insulators (MTI) Cr:(Bi,Sb)$_2$Te$_3$ and V:(Bi,Sb)$_2$Te$_3$, the search for the exchange coupling mechanisms underlying the onset of ferromagnetism has been a central issue, and a variety of different scenarios have been put forward. By combining resonant photoemission, X-ray magnetic dichroism and mul…
▽ More
Since the discovery of the quantum anomalous Hall effect in the magnetically doped topological insulators (MTI) Cr:(Bi,Sb)$_2$Te$_3$ and V:(Bi,Sb)$_2$Te$_3$, the search for the exchange coupling mechanisms underlying the onset of ferromagnetism has been a central issue, and a variety of different scenarios have been put forward. By combining resonant photoemission, X-ray magnetic dichroism and multiplet ligand field theory, we determine the local electronic and magnetic configurations of V and Cr impurities in (Bi,Sb)$_2$Te$_3$. While strong pd hybridisation is found for both dopant types, their 3d densities of states show pronounced differences. State-of-the-art first-principles calculations show how these impurity states mediate characteristic short-range pd exchange interactions, whose strength sensitively varies with the position of the 3d states relative to the Fermi level. Measurements on films with varying host stoichiometry support this trend. Our results establish the essential role of impurity-state mediated exchange interactions in the magnetism of MTI.
△ Less
Submitted 15 January, 2020;
originally announced January 2020.
-
Time-information uncertainty relations in thermodynamics
Authors:
Schuyler B. Nicholson,
Luis Pedro Garcia-Pintos,
Adolfo del Campo,
Jason R. Green
Abstract:
Physical systems that power motion and create structure in a fixed amount of time dissipate energy and produce entropy. Whether living or synthetic, systems performing these dynamic functions must balance dissipation and speed. Here, we show that rates of energy and entropy exchange are subject to a speed limit -- a time-information uncertainty relation -- imposed by the rates of change in the inf…
▽ More
Physical systems that power motion and create structure in a fixed amount of time dissipate energy and produce entropy. Whether living or synthetic, systems performing these dynamic functions must balance dissipation and speed. Here, we show that rates of energy and entropy exchange are subject to a speed limit -- a time-information uncertainty relation -- imposed by the rates of change in the information content of the system. This uncertainty relation bounds the time that elapses before the change in a thermodynamic quantity has the same magnitude as its initial standard deviation. From this general bound, we establish a family of speed limits for heat, work, entropy production, and entropy flow depending on the experimental constraints on the system. In all of these inequalities, the time scale of transient dynamical fluctuations is universally bounded by the Fisher information. Moreover, they all have a mathematical form that mirrors the Mandelstam-Tamm version of the time-energy uncertainty relation in quantum mechanics. These bounds on the speed of arbitrary observables apply to transient systems away from thermodynamic equilibrium, independent of the physical assumptions about the stochastic dynamics or their function.
△ Less
Submitted 15 January, 2020;
originally announced January 2020.
-
GW190425: Observation of a Compact Binary Coalescence with Total Mass $\sim 3.4 M_{\odot}$
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
B. P. Abbott,
R. Abbott,
T. D. Abbott,
S. Abraham,
F. Acernese,
K. Ackley,
C. Adams,
R. X. Adhikari,
V. B. Adya,
C. Affeldt,
M. Agathos,
K. Agatsuma,
N. Aggarwal,
O. D. Aguiar,
L. Aiello,
A. Ain,
P. Ajith,
G. Allen,
A. Allocca,
M. A. Aloy,
P. A. Altin,
A. Amato,
S. Anand
, et al. (1177 additional authors not shown)
Abstract:
On 2019 April 25, the LIGO Livingston detector observed a compact binary coalescence with signal-to-noise ratio 12.9. The Virgo detector was also taking data that did not contribute to detection due to a low signal-to-noise ratio, but were used for subsequent parameter estimation. The 90% credible intervals for the component masses range from 1.12 to 2.52 $M_{\odot}$ (1.45 to 1.88 $M_{\odot}$ if w…
▽ More
On 2019 April 25, the LIGO Livingston detector observed a compact binary coalescence with signal-to-noise ratio 12.9. The Virgo detector was also taking data that did not contribute to detection due to a low signal-to-noise ratio, but were used for subsequent parameter estimation. The 90% credible intervals for the component masses range from 1.12 to 2.52 $M_{\odot}$ (1.45 to 1.88 $M_{\odot}$ if we restrict the dimensionless component spin magnitudes to be smaller than 0.05). These mass parameters are consistent with the individual binary components being neutron stars. However, both the source-frame chirp mass $1.44^{+0.02}_{-0.02} M_{\odot}$ and the total mass $3.4^{+0.3}_{-0.1}\,M_{\odot}$ of this system are significantly larger than those of any other known binary neutron star system. The possibility that one or both binary components of the system are black holes cannot be ruled out from gravitational-wave data. We discuss possible origins of the system based on its inconsistency with the known Galactic binary neutron star population. Under the assumption that the signal was produced by a binary neutron star coalescence, the local rate of neutron star mergers is updated to $250-2810 \text{Gpc}^{-3}\text{yr}^{-1}$.
△ Less
Submitted 7 April, 2020; v1 submitted 6 January, 2020;
originally announced January 2020.
-
A Joint Fermi-GBM and LIGO/Virgo Analysis of Compact Binary Mergers From the First and Second Gravitational-wave Observing Runs
Authors:
The Fermi Gamma-ray Burst Monitor Team,
the LIGO Scientific Collaboration,
the Virgo Collaboration,
:,
R. Hamburg,
C. Fletcher,
E. Burns,
A. Goldstein,
E. Bissaldi,
M. S. Briggs,
W. H. Cleveland,
M. M. Giles,
C. M. Hui,
D. Kocevski,
S. Lesage,
B. Mailyan,
C. Malacaria,
S. Poolakkil,
R. Preece,
O. J. Roberts,
P. Veres,
A. von Kienlin,
C. A. Wilson-Hodge,
J. Wood,
R. Abbott
, et al. (1241 additional authors not shown)
Abstract:
We present results from offline searches of Fermi Gamma-ray Burst Monitor (GBM) data for gamma-ray transients coincident with the compact binary coalescences observed by the gravitational-wave (GW) detectors Advanced LIGO and Advanced Virgo during their first and second observing runs. In particular, we perform follow-up for both confirmed events and low significance candidates reported in the LIG…
▽ More
We present results from offline searches of Fermi Gamma-ray Burst Monitor (GBM) data for gamma-ray transients coincident with the compact binary coalescences observed by the gravitational-wave (GW) detectors Advanced LIGO and Advanced Virgo during their first and second observing runs. In particular, we perform follow-up for both confirmed events and low significance candidates reported in the LIGO/Virgo catalog GWTC-1. We search for temporal coincidences between these GW signals and GBM triggered gamma-ray bursts (GRBs). We also use the GBM Untargeted and Targeted subthreshold searches to find coincident gamma-rays below the on-board triggering threshold. This work implements a refined statistical approach by incorporating GW astrophysical source probabilities and GBM visibilities of LIGO/Virgo sky localizations to search for cumulative signatures of coincident subthreshold gamma-rays. All search methods recover the short gamma-ray burst GRB 170817A occurring ~1.7 s after the binary neutron star merger GW170817. We also present results from a new search seeking GBM counterparts to LIGO single-interferometer triggers. This search finds a candidate joint event, but given the nature of the GBM signal and localization, as well as the high joint false alarm rate of $1.1 \times 10^{-6}$ Hz, we do not consider it an astrophysical association. We find no additional coincidences.
△ Less
Submitted 24 February, 2020; v1 submitted 3 January, 2020;
originally announced January 2020.
-
Strain-Driven Disproportionation at a Correlated Oxide Metal-Insulator Transition
Authors:
T. H. Kim,
T. R. Paudel,
R. J. Green,
K. Song,
H. -S. Lee,
S. -Y. Choi,
J. Irwin,
B. Noesges,
L. J. Brillson,
M. S. Rzchowski,
G. A. Sawatzky,
E. Y. Tsymbal,
C. B. Eom
Abstract:
Metal-to-insulator phase transitions in complex oxide thin films are exciting phenomena which may be useful for device applications, but in many cases the physical mechanism responsible for the transition is not fully understood. Here we demonstrate that epitaxial strain generates local disproportionation of the NiO6 octahedra, driven through changes in the oxygen stoichiometry, and that this dire…
▽ More
Metal-to-insulator phase transitions in complex oxide thin films are exciting phenomena which may be useful for device applications, but in many cases the physical mechanism responsible for the transition is not fully understood. Here we demonstrate that epitaxial strain generates local disproportionation of the NiO6 octahedra, driven through changes in the oxygen stoichiometry, and that this directly modifies the metal-to-insulator phase transition in epitaxial (001) NdNiO3 thin films. Theoretically, we predict that the Ni-O-Ni bond angle decreases, while octahedral tilt and local disproportionation of the NiO6 octahedra increases resulting in a small band gap in otherwise metallic system. This is driven by an increase in oxygen vacancy concentration in the rare-earth nickelates with increasing in-plane biaxial tensile strain. Experimentally, we find an increase in pseudocubic unit-cell volume and resistivity with increasing biaxial tensile strain, corroborating our theoretical predictions. With electron energy loss spectroscopy and x-ray absorption, we find a reduction of the Ni valence with increasing tensile strain. These results indicate that epitaxial strain modifies the oxygen stoichiometry of rare-earth perovskite thin films and through this mechanism affect the metal-to-insulator phase transition in these compounds.
△ Less
Submitted 29 December, 2019;
originally announced December 2019.
-
Open data from the first and second observing runs of Advanced LIGO and Advanced Virgo
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
R. Abbott,
T. D. Abbott,
S. Abraham,
F. Acernese,
K. Ackley,
C. Adams,
R. X. Adhikari,
V. B. Adya,
C. Affeldt,
M. Agathos,
K. Agatsuma,
N. Aggarwal,
O. D. Aguiar,
A. Aich,
L. Aiello,
A. Ain,
P. Ajith,
G. Allen,
A. Allocca,
P. A. Altin,
A. Amato,
S. Anand,
A. Ananyeva
, et al. (1223 additional authors not shown)
Abstract:
Advanced LIGO and Advanced Virgo are actively monitoring the sky and collecting gravitational-wave strain data with sufficient sensitivity to detect signals routinely. In this paper we describe the data recorded by these instruments during their first and second observing runs. The main data products are the gravitational-wave strain arrays, released as time series sampled at 16384 Hz. The dataset…
▽ More
Advanced LIGO and Advanced Virgo are actively monitoring the sky and collecting gravitational-wave strain data with sufficient sensitivity to detect signals routinely. In this paper we describe the data recorded by these instruments during their first and second observing runs. The main data products are the gravitational-wave strain arrays, released as time series sampled at 16384 Hz. The datasets that include this strain measurement can be freely accessed through the Gravitational Wave Open Science Center at http://gw-openscience.org, together with data-quality information essential for the analysis of LIGO and Virgo data, documentation, tutorials, and supporting software.
△ Less
Submitted 25 January, 2021; v1 submitted 25 December, 2019;
originally announced December 2019.
-
On theoretical upper limits for valid timesteps of implicit ODE methods
Authors:
K. R. Green,
G. W. Patrick,
R. J. Spiteri
Abstract:
Implicit methods for the numerical solution of initial-value problems may admit multiple solutions at any given time step. Accordingly, their nonlinear solvers may converge to any of these solutions. Below a critical timestep, exactly one of the solutions (the consistent solution) occurs on a solution branch (the principal branch) that can be continuously and monotonically continued back to zero t…
▽ More
Implicit methods for the numerical solution of initial-value problems may admit multiple solutions at any given time step. Accordingly, their nonlinear solvers may converge to any of these solutions. Below a critical timestep, exactly one of the solutions (the consistent solution) occurs on a solution branch (the principal branch) that can be continuously and monotonically continued back to zero timestep.
Standard step-size control can promote convergence to consistent solutions by adjusting the timestep to maintain an error estimate below a given tolerance. However, simulations for symplectic systems or large physical systems are often run with constant timesteps and are thus more susceptible to convergence to inconsistent solutions. Because simulations cannot be reliably continued from inconsistent solutions, the critical timestep is a theoretical upper bound for valid timesteps.
△ Less
Submitted 18 December, 2019;
originally announced December 2019.
-
Excited hairy black holes: dynamical construction and level transitions
Authors:
Pablo Bosch,
Stephen R. Green,
Luis Lehner,
Hugo Roussille
Abstract:
We study the dynamics of unstable Reissner-Nordström anti-de Sitter black holes under charged scalar field perturbations in spherical symmetry. We unravel their general behavior and approach to the final equlibrium state. In the first part of this work, we present a numerical analysis of massive charged scalar field quasinormal modes. We identify the known mode families - superradiant modes, zero-…
▽ More
We study the dynamics of unstable Reissner-Nordström anti-de Sitter black holes under charged scalar field perturbations in spherical symmetry. We unravel their general behavior and approach to the final equlibrium state. In the first part of this work, we present a numerical analysis of massive charged scalar field quasinormal modes. We identify the known mode families - superradiant modes, zero-damped modes, AdS modes, and the near-horizon mode - and we track their migration under variation of the black hole and field parameters. We show that the zero-damped modes become superradiantly unstable for large RNAdS with large gauge coupling; the leading unstable mode is identified with the near-horizon condensation instability. In the second part, we present results of numerical simulations of perturbed large RNAdS, showing the nonlinear development of these unstable modes. For generic initial conditions, charge and mass are ransferred from the black hole to the scalar field, until an equilibrium solution with a scalar condensate is reached. We use results from the linear analysis, however, to select special initial data corresponding to an unstable overtone mode. We find that these data evolve to produce a new equilibrium state - an excited hairy black hole with the scalar condensate in an overtone configuration. This state is, however, unstable, and the black hole eventually decays to the generic end state. Nevertheless, this demonstrates the potential relevance of overtone modes as transients in black hole dynamics.
△ Less
Submitted 11 December, 2019;
originally announced December 2019.
-
Developments in the position-space approach to the HLbL contribution to the muon $g-2$ on the lattice
Authors:
Nils Asmussen,
En-Hung Chao,
Antoine Gérardin,
Jeremy R. Green,
Renwick J. Hudspith,
Harvey B. Meyer,
Andreas Nyffeler
Abstract:
The measurement of the anomalous magnetic moment of the muon and its prediction allow for a high-precision test of the Standard Model (SM). In this proceedings article we present ongoing work combining lattice QCD and continuum QED in order to determine an important SM contribution to the magnetic moment, the hadronic light-by-light contribution. We compute the quark-connected contribution in the…
▽ More
The measurement of the anomalous magnetic moment of the muon and its prediction allow for a high-precision test of the Standard Model (SM). In this proceedings article we present ongoing work combining lattice QCD and continuum QED in order to determine an important SM contribution to the magnetic moment, the hadronic light-by-light contribution. We compute the quark-connected contribution in the Mainz position-space approach and investigate the long-distance part of our data using calculations of the $π^0$-pole and charged pion loop contributions.
△ Less
Submitted 13 November, 2019;
originally announced November 2019.
-
GRAVITAS: A Model Checking Based Planning and Goal Reasoning Framework for Autonomous Systems
Authors:
Hadrien Bride,
Jin Song Dong,
Ryan Green,
Zhe Hou,
Brendan Mahony,
Martin Oxenham
Abstract:
While AI techniques have found many successful applications in autonomous systems, many of them permit behaviours that are difficult to interpret and may lead to uncertain results. We follow the "verification as planning" paradigm and propose to use model checking techniques to solve planning and goal reasoning problems for autonomous systems. We give a new formulation of Goal Task Network (GTN) t…
▽ More
While AI techniques have found many successful applications in autonomous systems, many of them permit behaviours that are difficult to interpret and may lead to uncertain results. We follow the "verification as planning" paradigm and propose to use model checking techniques to solve planning and goal reasoning problems for autonomous systems. We give a new formulation of Goal Task Network (GTN) that is tailored for our model checking based framework. We then provide a systematic method that models GTNs in the model checker Process Analysis Toolkit (PAT). We present our planning and goal reasoning system as a framework called Goal Reasoning And Verification for Independent Trusted Autonomous Systems (GRAVITAS) and discuss how it helps provide trustworthy plans in an uncertain environment. Finally, we demonstrate the proposed ideas in an experiment that simulates a survey mission performed by the REMUS-100 autonomous underwater vehicle.
△ Less
Submitted 3 October, 2019;
originally announced October 2019.
-
Rho resonance, timelike pion form factor, and implications for lattice studies of the hadronic vacuum polarisation
Authors:
Felix Erben,
Jeremy R. Green,
Daniel Mohler,
Hartmut Wittig
Abstract:
We study isospin-1 P-wave $ππ$ scattering in lattice QCD with two flavours of O($a$) improved Wilson fermions. For pion masses ranging from $m_π=265$ MeV to $m_π=437$ MeV, we determine the energy spectrum in the centre-of-mass frame and in three moving frames. We obtain the scattering phase shifts using Lüscher's finite-volume quantisation condition. Fitting the dependence of the phase shifts on t…
▽ More
We study isospin-1 P-wave $ππ$ scattering in lattice QCD with two flavours of O($a$) improved Wilson fermions. For pion masses ranging from $m_π=265$ MeV to $m_π=437$ MeV, we determine the energy spectrum in the centre-of-mass frame and in three moving frames. We obtain the scattering phase shifts using Lüscher's finite-volume quantisation condition. Fitting the dependence of the phase shifts on the scattering momentum to a Breit-Wigner form allows us to determine the resonance parameters $m_ρ$ and $g_{ρππ}$. By combining the scattering phase shifts with the decay matrix element of the vector current, we calculate the timelike pion form factor, $F_π$, and compare the results to the Gounaris-Sakurai representation of the form factor in terms of the resonance parameters. In addition, we fit our data for the form factor to the functional form suggested by the Omnès representation, which allows for the extraction of the charge radius of the pion. As a further application, we discuss the long-distance behaviour of the vector correlator, which is dominated by the two-pion channel. We reconstruct the long-distance part in two ways: one based on the finite-volume energies and matrix elements and the other based on $F_π$. It is shown that this part can be accurately constrained using the reconstructions, which has important consequences for lattice calculations of the hadronic vacuum polarisation contribution to the muon anomalous magnetic moment.
△ Less
Submitted 2 October, 2019;
originally announced October 2019.
-
Strategic Scientific Plan for Gemini Observatory
Authors:
J. P. Blakeslee,
A. Adamson,
C. Davis,
R. Díaz,
B. Miller,
A. Peck,
R. Rutten,
G. Sivo,
J. Thomas-Osip,
T. Boroson,
R. Carrasco,
E. Dennihy,
M. Díaz,
L. Ferrarese,
R. Green,
P. Hirst,
N. Hwang,
I. Jørgensen,
H. Kim,
S. Kleinman,
K. Labrie,
T. Lee,
J. Lotz,
S. Leggett,
L. Medina
, et al. (8 additional authors not shown)
Abstract:
We present the Strategic Scientific Plan (SSP) for the direction and activities of the Gemini Observatory in the 2020s. The overarching goal is to ensure that Gemini best uses the available resources to serve the needs of its international user community throughout the coming decade. The actionable items fall into three general categories: (1) preserving Gemini's current facilities and strengths;…
▽ More
We present the Strategic Scientific Plan (SSP) for the direction and activities of the Gemini Observatory in the 2020s. The overarching goal is to ensure that Gemini best uses the available resources to serve the needs of its international user community throughout the coming decade. The actionable items fall into three general categories: (1) preserving Gemini's current facilities and strengths; (2) developing instrumentation and software systems, including data pipelines, to enable new scientific capabilities that build on those strengths; (3) strategizing how visiting instruments can deliver additional valuable capabilities. We provide a high-level timeline (schematically illustrated in one figure) for the main developments discussed in this SSP. The schedule is ambitious, but in light of the recent Gemini in the Era of Multi-Messenger Astronomy (GEMMA) award from the NSF, the plan becomes achievable. Lists of milestones are given for gauging progress. As these milestones are reached and new instruments become available, some current instruments will need to be retired; we make recommendations in this regard. The final section concludes by reemphasizing the importance of a strong partnership committed to the needs of all members.
△ Less
Submitted 19 September, 2019;
originally announced September 2019.
-
A guide to LIGO-Virgo detector noise and extraction of transient gravitational-wave signals
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
B. P. Abbott,
R. Abbott,
T. D. Abbott,
S. Abraham,
F. Acernese,
K. Ackley,
C. Adams,
V. B. Adya,
C. Affeldt,
M. Agathos,
K. Agatsuma,
N. Aggarwal,
O. D. Aguiar,
L. Aiello,
A. Ain,
P. Ajith,
T. Alford,
G. Allen,
A. Allocca,
M. A. Aloy,
P. A. Altin,
A. Amato,
A. Ananyeva
, et al. (1113 additional authors not shown)
Abstract:
The LIGO Scientific Collaboration and the Virgo Collaboration have cataloged eleven confidently detected gravitational-wave events during the first two observing runs of the advanced detector era. All eleven events were consistent with being from well-modeled mergers between compact stellar-mass objects: black holes or neutron stars. The data around the time of each of these events have been made…
▽ More
The LIGO Scientific Collaboration and the Virgo Collaboration have cataloged eleven confidently detected gravitational-wave events during the first two observing runs of the advanced detector era. All eleven events were consistent with being from well-modeled mergers between compact stellar-mass objects: black holes or neutron stars. The data around the time of each of these events have been made publicly available through the gravitational-wave open science center. The entirety of the gravitational-wave strain data from the first and second observing runs have also now been made publicly available. There is considerable interest among the broad scientific community in understanding the data and methods used in the analyses. In this paper, we provide an overview of the detector noise properties and the data analysis techniques used to detect gravitational-wave signals and infer the source properties. We describe some of the checks that are performed to validate the analyses and results from the observations of gravitational-wave events. We also address concerns that have been raised about various properties of LIGO-Virgo detector noise and the correctness of our analyses as applied to the resulting data.
△ Less
Submitted 10 February, 2020; v1 submitted 29 August, 2019;
originally announced August 2019.
-
Teukolsky formalism for nonlinear Kerr perturbations
Authors:
Stephen R. Green,
Stefan Hollands,
Peter Zimmerman
Abstract:
We develop a formalism to treat higher order (nonlinear) metric perturbations of the Kerr spacetime in a Teukolsky framework. We first show that solutions to the linearized Einstein equation with nonvanishing stress tensor can be decomposed into a pure gauge part plus a zero mode (infinitesimal perturbation of the mass and spin) plus a perturbation arising from a certain scalar ("Debye-Hertz") pot…
▽ More
We develop a formalism to treat higher order (nonlinear) metric perturbations of the Kerr spacetime in a Teukolsky framework. We first show that solutions to the linearized Einstein equation with nonvanishing stress tensor can be decomposed into a pure gauge part plus a zero mode (infinitesimal perturbation of the mass and spin) plus a perturbation arising from a certain scalar ("Debye-Hertz") potential, plus a so-called "corrector tensor." The scalar potential is a solution to the spin $-2$ Teukolsky equation with a source. This source, as well as the tetrad components of the corrector tensor, are obtained by solving certain decoupled ordinary differential equations involving the stress tensor. As we show, solving these ordinary differential equations reduces simply to integrations in the coordinate $r$ in outgoing Kerr-Newman coordinates, so in this sense, the problem is reduced to the Teukolsky equation with source, which can be treated by a separation of variables ansatz. Since higher order perturbations are subject to a linearized Einstein equation with a stress tensor obtained from the lower order perturbations, our method also applies iteratively to the higher order metric perturbations, and could thus be used to analyze the nonlinear coupling of perturbations in the near-extremal Kerr spacetime, where weakly turbulent behavior has been conjectured to occur. Our method could also be applied to the study of perturbations generated by a pointlike body traveling on a timelike geodesic in Kerr, which is relevant to the extreme mass ratio inspiral problem.
△ Less
Submitted 14 September, 2020; v1 submitted 24 August, 2019;
originally announced August 2019.