Showing 1–50 of 68 results for author: Gibbs, A

Two Combinatorial Interpretations of Rascal Numbers

Authors: Amelia Gibbs, Brian K. Miceli

Abstract: The main goal of this paper is to assign two combinatorial interpretations to the elements of the Rascal Triangle defined by Angorro et al. The first interpretation involves counting ascents in binary words, while the second interpretation involves pattern avoidance in the ascent sequences studied by Duncan and Steingrímsson, and related directly to a more recent paper by Baxter and Pudwell. The main goal of this paper is to assign two combinatorial interpretations to the elements of the Rascal Triangle defined by Angorro et al. The first interpretation involves counting ascents in binary words, while the second interpretation involves pattern avoidance in the ascent sequences studied by Duncan and Steingrímsson, and related directly to a more recent paper by Baxter and Pudwell. △ Less

Submitted 17 May, 2024; originally announced May 2024.

MSC Class: 05

Metastable ordered states induced by low temperature annealing of δ-Ag2/3V2O5

Authors: T. Kubo, K. Kojima, N. Katayama, T. Runčevski, R. E. Dinnebier, A. S. Gibbs, M. Isobe, H. Sawa

Abstract: In δ-Ag2/3V2O5 with charge degrees of freedom in V, it is known that the charge ordering state and physical properties of V that appear at low temperatures depend strongly on the ordering state of Ag. In this study, we focused on the Ag ions in the interlayer and studied the structure using synchrotron radiation powder diffraction in dependence on temperature. We found that when the sample is slow… ▽ More In δ-Ag2/3V2O5 with charge degrees of freedom in V, it is known that the charge ordering state and physical properties of V that appear at low temperatures depend strongly on the ordering state of Ag. In this study, we focused on the Ag ions in the interlayer and studied the structure using synchrotron radiation powder diffraction in dependence on temperature. We found that when the sample is slowly cooled from room temperature and ordering occurs at the Ag sites, V4+/V5+ charge ordering of V and subsequent V4+-V4+ structural dimers are produced. Although quenching the sample from room temperature suppresses the ordering of Ag, annealing at around 160 K promotes partial ordering of Ag and allows a metastable phase to be realized. This metastable phase is maintained even when the temperature is lowered again, producing a remarkable change in low-temperature properties. These results indicate that the ordered state of Ag, which is the key to control the charge-ordered state and physical properties, can be controlled by low-temperature annealing. The results of this study may provide a methodology for the realization of metastable states in a wide range of material groups of vanadium compounds, where competition among various charge ordered states underlies the physical properties. △ Less

Submitted 11 March, 2024; originally announced March 2024.

Comments: 10 pages, 6 figures

Direct-imaging Discovery of a Substellar Companion Orbiting the Accelerating Variable Star, HIP 39017

Authors: Taylor L. Tobin, Thayne Currie, Yiting Li, Jeffrey Chilcote, Timothy D. Brandt, Brianna Lacy, Masayuki Kuzuhara, Maria Vincent, Mona El Morsy, Vincent Deo, Jonathan P. Williams, Olivier Guyon, Julien Lozi, Sebastien Vievard, Nour Skaf, Kyohoon Ahn, Tyler Groff, N. Jeremy Kasdin, Taichi Uyama, Motohide Tamura, Aidan Gibbs, Briley L. Lewis, Rachel Bowens-Rubin, Maïssa Salama, Qier An , et al. (1 additional authors not shown)

Abstract: We present the direct-imaging discovery of a substellar companion (a massive planet or low-mass brown dwarf) to the young, $γ$ Doradus-type variable star, HIP 39017 (HD 65526). The companion's SCExAO/CHARIS JHK ($1.1-2.4μ$m) spectrum and Keck/NIRC2 L$^{\prime}$ photometry indicate that it is an L/T transition object. A comparison of the JHK+L$^{\prime}$ spectrum to several atmospheric model grids… ▽ More We present the direct-imaging discovery of a substellar companion (a massive planet or low-mass brown dwarf) to the young, $γ$ Doradus-type variable star, HIP 39017 (HD 65526). The companion's SCExAO/CHARIS JHK ($1.1-2.4μ$m) spectrum and Keck/NIRC2 L$^{\prime}$ photometry indicate that it is an L/T transition object. A comparison of the JHK+L$^{\prime}$ spectrum to several atmospheric model grids finds a significantly better fit to cloudy models than cloudless models. Orbit modeling with relative astrometry and precision stellar astrometry from Hipparcos and Gaia yields a semi-major axis of $23.8^{+8.7}_{-6.1}$ au, a dynamical companion mass of $30^{+31}_{-12}$~M$_J$, and a mass ratio of $\sim$1.9\%, properties most consistent with low-mass brown dwarfs. However, its mass estimated from luminosity models is a lower $\sim$13.8 $M_{\rm J}$ due to an estimated young age ($\lesssim$ 115 Myr); using a weighted posterior distribution informed by conservative mass constraints from luminosity evolutionary models yields a lower dynamical mass of $23.6_{-7.4}^{+9.1}$~M$_J$ and a mass ratio of $\sim$1.4\%. Analysis of the host star's multi-frequency $γ$ Dor-type pulsations, astrometric monitoring of HIP 39017b, and Gaia Data Release 4 astrometry of the star will clarify the system age and better constrain the mass and orbit of the companion. This discovery further reinforces the improved efficiency of targeted direct-imaging campaigns informed by long-baseline, precision stellar astrometry. △ Less

Submitted 15 May, 2024; v1 submitted 6 March, 2024; originally announced March 2024.

Comments: Published in AJ: April 9, 2024

Journal ref: 2024, AJ, 167, 205

Using theory and experiments of spheres moving near boundaries to optimize the method of images for regularized Stokeslets

Authors: Hoa Nguyen, Amelia Gibbs, Frank Healy, Orrin Shindell, Ricardo Cortez, Kate Brown, Jonathan McCoy, Bruce Rodenborn

Abstract: The general system of images for regularized Stokeslets (GSIRS) developed by Cortez and Varela (2015) is used extensively to model Stokes flow phenomena such as microorganisms swimming near a boundary. Our collaborative team uses dynamically similar scaled macroscopic experiments to test theories for forces and torques on spheres moving near a boundary and use these data and the method of regulari… ▽ More The general system of images for regularized Stokeslets (GSIRS) developed by Cortez and Varela (2015) is used extensively to model Stokes flow phenomena such as microorganisms swimming near a boundary. Our collaborative team uses dynamically similar scaled macroscopic experiments to test theories for forces and torques on spheres moving near a boundary and use these data and the method of regularized Stokeslets (MRS) created by Cortez et al. (2015) to calibrate the GSIRS. We find excellent agreement between theory and experiments, which provides the first experimental validation of such series solutions for spheres moving near an infinite plane boundary. We test two surface discretization methods commonly used in the literature: the six-patch method and the spherical centroidal Voronoi tessellation (SCVT) method. Our data show that the SCVT method provides the most accurate results when the motional symmetry is broken by the presence of a boundary. We use theory and the MRS to find optimal values for the regularization parameter in free space and show that the optimal regularization parameter values can be fit with simple formulae when using the SCVT method, so that other researchers have an easy reference. We also present a regularization function with higher order accuracy when compared with the regularization function previously introduced by Cortez et al. (2005). The simulated force and torque values compare very well with experiments and theory for a wide range of boundary distances. But, we find the simulations lose accuracy when the gap between the edge of the sphere and the wall is smaller than the average distance between grid points using the SCVT method. Our computational parameters and MATLAB and PYTHON implementations of the Lee and Leal (1980) theory provide researchers with important resources to optimize the numerical simulations of spheres moving near boundaries. △ Less

Submitted 11 July, 2024; v1 submitted 29 January, 2024; originally announced January 2024.

Comments: 38 pages, 10 figures

An efficient frequency-independent numerical method for computing the far-field pattern induced by polygonal obstacles

Authors: A. Gibbs, S. Langdon

Abstract: For problems of time-harmonic scattering by rational polygonal obstacles, embedding formulae express the far-field pattern induced by any incident plane wave in terms of the far-field patterns for a relatively small (frequency-independent) set of canonical incident angles. Although these remarkable formulae are exact in theory, here we demonstrate that: (i) they are highly sensitive to numerical e… ▽ More For problems of time-harmonic scattering by rational polygonal obstacles, embedding formulae express the far-field pattern induced by any incident plane wave in terms of the far-field patterns for a relatively small (frequency-independent) set of canonical incident angles. Although these remarkable formulae are exact in theory, here we demonstrate that: (i) they are highly sensitive to numerical errors in practice, and (ii) direct calculation of the coefficients in these formulae may be impossible for particular sets of canonical incident angles, even in exact arithmetic. Only by overcoming these practical issues can embedding formulae provide a highly efficient approach to computing the far-field pattern induced by a large number of incident angles. Here we address challenges (i) and (ii), supporting our theory with numerical experiments. Challenge (i) is solved using techniques from computational complex analysis: we reformulate the embedding formula as a complex contour integral and prove that this is much less sensitive to numerical errors. In practice, this contour integral can be efficiently evaluated by residue calculus. Challenge (ii) is addressed using techniques from numerical linear algebra: we oversample, considering more canonical incident angles than are necessary, thus expanding the set of valid coefficient vectors. The coefficient vector can then be selected using either a least squares approach or column subset selection. △ Less

Submitted 5 March, 2024; v1 submitted 26 October, 2023; originally announced October 2023.

MSC Class: 35J05; 78A45; 30E20; 65F20

A caustic terminating at an inflection point

Authors: J. R. Ockendon, H. Ockendon, R. H. Tew, D. P. Hewett, A. Gibbs

Abstract: We present an asymptotic and numerical study of the evolution of an incoming wavefield which has a caustic close to a curve with an inflection point. Our results reveal the emergence of a wavefield which resembles that of a shadow boundary but has a maximum amplitude along the tangent at the inflection point. We present an asymptotic and numerical study of the evolution of an incoming wavefield which has a caustic close to a curve with an inflection point. Our results reveal the emergence of a wavefield which resembles that of a shadow boundary but has a maximum amplitude along the tangent at the inflection point. △ Less

Submitted 6 December, 2023; v1 submitted 20 October, 2023; originally announced October 2023.

SAGUARO: Time-domain Infrastructure for the Fourth Gravitational-wave Observing Run and Beyond

Authors: Griffin Hosseinzadeh, Kerry Paterson, Jillian C. Rastinejad, Manisha Shrestha, Philip N. Daly, Michael J. Lundquist, David J. Sand, Wen-fai Fong, K. Azalee Bostroem, Saarah Hall, Samuel D. Wyatt, Alex R. Gibbs, Eric Christensen, William Lindstrom, Jonathan Nation, Joseph Chatelain, Curtis McCully

Abstract: We present upgraded infrastructure for Searches after Gravitational Waves Using ARizona Observatories (SAGUARO) during LIGO, Virgo, and KAGRA's fourth gravitational-wave (GW) observing run (O4). These upgrades implement many of the lessons we learned after a comprehensive analysis of potential electromagnetic counterparts to the GWs discovered during the previous observing run. We have developed a… ▽ More We present upgraded infrastructure for Searches after Gravitational Waves Using ARizona Observatories (SAGUARO) during LIGO, Virgo, and KAGRA's fourth gravitational-wave (GW) observing run (O4). These upgrades implement many of the lessons we learned after a comprehensive analysis of potential electromagnetic counterparts to the GWs discovered during the previous observing run. We have developed a new web-based target and observation manager (TOM) that allows us to coordinate sky surveys, vet potential counterparts, and trigger follow-up observations from one centralized portal. The TOM includes software that aggregates all publicly available information on the light curves and possible host galaxies of targets, allowing us to rule out potential contaminants like active galactic nuclei, variable stars, solar-system objects, and preexisting supernovae, as well as to assess the viability of any plausible counterparts. We have also upgraded our image-subtraction pipeline by assembling deeper reference images and training a new neural network-based real-bogus classifier. These infrastructure upgrades will aid coordination by enabling the prompt reporting of observations, discoveries, and analysis to the GW follow-up community, and put SAGUARO in an advantageous position to discover kilonovae in the remainder of O4 and beyond. Many elements of our open-source software stack have broad utility beyond multimessenger astronomy, and will be particularly relevant in the "big data" era of transient discoveries by the Vera C. Rubin Observatory. △ Less

Submitted 14 March, 2024; v1 submitted 12 October, 2023; originally announced October 2023.

Comments: updated to match accepted version

Journal ref: The Astrophysical Journal, 964:35 (18pp), 2024 March 20

The development of HISPEC for Keck and MODHIS for TMT: science cases and predicted sensitivities

Authors: Quinn M. Konopacky, Ashley D. Baker, Dimitri Mawet, Michael P. Fitzgerald, Nemanja Jovanovic, Charles Beichman, Garreth Ruane, Rob Bertz, Hiroshi Terada, Richard Dekany, Larry Lingvay, Marc Kassis, David Anderson, Motohide Tamura, Bjorn Benneke, Thomas Beatty, Tuan Do, Shogo Nishiyama, Peter Plavchan, Jason Wang, Ji Wang, Adam Burgasser, Jean-Baptiste Ruffio, Huihao Zhang, Aaron Brown , et al. (50 additional authors not shown)

Abstract: HISPEC is a new, high-resolution near-infrared spectrograph being designed for the W.M. Keck II telescope. By offering single-shot, R=100,000 between 0.98 - 2.5 um, HISPEC will enable spectroscopy of transiting and non-transiting exoplanets in close orbits, direct high-contrast detection and spectroscopy of spatially separated substellar companions, and exoplanet dynamical mass and orbit measureme… ▽ More HISPEC is a new, high-resolution near-infrared spectrograph being designed for the W.M. Keck II telescope. By offering single-shot, R=100,000 between 0.98 - 2.5 um, HISPEC will enable spectroscopy of transiting and non-transiting exoplanets in close orbits, direct high-contrast detection and spectroscopy of spatially separated substellar companions, and exoplanet dynamical mass and orbit measurements using precision radial velocity monitoring calibrated with a suite of state-of-the-art absolute and relative wavelength references. MODHIS is the counterpart to HISPEC for the Thirty Meter Telescope and is being developed in parallel with similar scientific goals. In this proceeding, we provide a brief overview of the current design of both instruments, and the requirements for the two spectrographs as guided by the scientific goals for each. We then outline the current science case for HISPEC and MODHIS, with focuses on the science enabled for exoplanet discovery and characterization. We also provide updated sensitivity curves for both instruments, in terms of both signal-to-noise ratio and predicted radial velocity precision. △ Less

Submitted 19 September, 2023; originally announced September 2023.

Comments: 25 pages, 9 figures. To appear in the Proceedings of SPIE: Techniques and Instrumentation for Detection of Exoplanets XI, vol. 12680 (2023)

Integral equation methods for acoustic scattering by fractals

Authors: A. M. Caetano, S. N. Chandler-Wilde, X. Claeys, A. Gibbs, D. P. Hewett, A. Moiola

Abstract: We study sound-soft time-harmonic acoustic scattering by general scatterers, including fractal scatterers, in 2D and 3D space. For an arbitrary compact scatterer $Γ$ we reformulate the Dirichlet boundary value problem for the Helmholtz equation as a first kind integral equation (IE) on $Γ$ involving the Newton potential. The IE is well-posed, except possibly at a countable set of frequencies, and… ▽ More We study sound-soft time-harmonic acoustic scattering by general scatterers, including fractal scatterers, in 2D and 3D space. For an arbitrary compact scatterer $Γ$ we reformulate the Dirichlet boundary value problem for the Helmholtz equation as a first kind integral equation (IE) on $Γ$ involving the Newton potential. The IE is well-posed, except possibly at a countable set of frequencies, and reduces to existing single-layer boundary IEs when $Γ$ is the boundary of a bounded Lipschitz open set, a screen, or a multi-screen. When $Γ$ is uniformly of $d$-dimensional Hausdorff dimension in a sense we make precise (a $d$-set), the operator in our equation is an integral operator on $Γ$ with respect to $d$-dimensional Hausdorff measure, with kernel the Helmholtz fundamental solution, and we propose a piecewise-constant Galerkin discretization of the IE, which converges in the limit of vanishing mesh width. When $Γ$ is the fractal attractor of an iterated function system of contracting similarities we prove convergence rates under assumptions on $Γ$ and the IE solution, and describe a fully discrete implementation using recently proposed quadrature rules for singular integrals on fractals. We present numerical results for a range of examples and make our software available as a Julia code. △ Less

Submitted 30 April, 2024; v1 submitted 5 September, 2023; originally announced September 2023.

MSC Class: 28A80; 65R20; 78A45

Anisotropic Seebeck coefficient of $\mathrm{Sr}_2\mathrm{Ru}\mathrm{O}_4$ in the incoherent regime

Authors: Ramzy Daou, Sylvie Hébert, Gaël Grissonnanche, Elena Hassinger, Louis Taillefer, Haruka Taniguchi, Yoshiteru Maeno, Alexandra S. Gibbs, Andrew P. Mackenzie

Abstract: Intuitive entropic interpretations of the thermoelectric effect in metals predict an isotropic Seebeck coefficient at high temperatures in the incoherent regime even in anisotropic metals since entropy is not directional. $\mathrm{Sr}_2\mathrm{Ru}\mathrm{O}_4$ is an enigmatic material known for a well characterised anisotropic normal state and unconventional superconductivity. Recent ab-initio tra… ▽ More Intuitive entropic interpretations of the thermoelectric effect in metals predict an isotropic Seebeck coefficient at high temperatures in the incoherent regime even in anisotropic metals since entropy is not directional. $\mathrm{Sr}_2\mathrm{Ru}\mathrm{O}_4$ is an enigmatic material known for a well characterised anisotropic normal state and unconventional superconductivity. Recent ab-initio transport calculations of $\mathrm{Sr}_2\mathrm{Ru}\mathrm{O}_4$ that include the effect of strong electronic correlations predicted an enhanced high-temperature anisotropy of the Seebeck coefficient at temperatures above 300 K, but experimental evidence is missing. From measurements on clean $\mathrm{Sr}_2\mathrm{Ru}\mathrm{O}_4$ single crystals along both crystallographic directions, we find that the Seebeck coefficient becomes increasingly isotropic upon heating towards room temperature as generally expected. Above 300 K, however, $S$ acquires a new anisotropy which rises up to the highest temperatures measured (750 K), in qualitative agreement with calculations. This is a challenge to entropic interpretations and highlights the lack of an intuitive framework to understand the anisotropy of thermopower at high temperatures. △ Less

Submitted 4 September, 2023; originally announced September 2023.

Comments: 5 pages, 2 figures

Journal ref: Phys. Rev. B 108, L121106 (2023)

Numerical evaluation of oscillatory integrals via automated steepest descent contour deformation

Authors: A. Gibbs, D. P. Hewett, D. Huybrechs

Abstract: Steepest descent methods combining complex contour deformation with numerical quadrature provide an efficient and accurate approach for the evaluation of highly oscillatory integrals. However, unless the phase function governing the oscillation is particularly simple, their application requires a significant amount of a priori analysis and expert user input, to determine the appropriate contour de… ▽ More Steepest descent methods combining complex contour deformation with numerical quadrature provide an efficient and accurate approach for the evaluation of highly oscillatory integrals. However, unless the phase function governing the oscillation is particularly simple, their application requires a significant amount of a priori analysis and expert user input, to determine the appropriate contour deformation, and to deal with the non-uniformity in the accuracy of standard quadrature techniques associated with the coalescence of stationary points (saddle points) with each other, or with the endpoints of the original integration contour. In this paper we present a novel algorithm for the numerical evaluation of oscillatory integrals with general polynomial phase functions, which automates the contour deformation process and avoids the difficulties typically encountered with coalescing stationary points and endpoints. The inputs to the algorithm are simply the phase and amplitude functions, the endpoints and orientation of the original integration contour, and a small number of numerical parameters. By a series of numerical experiments we demonstrate that the algorithm is accurate and efficient over a large range of frequencies, even for examples with a large number of coalescing stationary points and with endpoints at infinity. As a particular application, we use our algorithm to evaluate cuspoid canonical integrals from scattering theory. A Matlab implementation of the algorithm is made available and is called PathFinder. △ Less

Submitted 6 December, 2023; v1 submitted 14 July, 2023; originally announced July 2023.

MSC Class: 65D30; 30E20; 41A60

Magnetic ground states and excitations in Zn-doped averieite -- a family of oxide-based $S=1/2$ kagome antiferromagnets

Authors: M. Georgopoulou, D. Boldrin, B. Fåk, P. Manuel, A. Gibbs, J. Ollivier, E. Suard, A. S. Wills

Abstract: Spin-1/2 kagome materials have recently attracted a resurgence of interest as they are considered an ideal host of the quantum spin liquid (QSL) state, which can underpin functionality such as superconductivity. Here we report the first synthesis and characterization of a new oxide-based distorted $S=1/2$ kagome antiferromagnet (KAFM) in the Zn$_x$Cu$_{5-x}$(VO$_4$)$_{2}$O$_2$CsCl (termed Zn… ▽ More Spin-1/2 kagome materials have recently attracted a resurgence of interest as they are considered an ideal host of the quantum spin liquid (QSL) state, which can underpin functionality such as superconductivity. Here we report the first synthesis and characterization of a new oxide-based distorted $S=1/2$ kagome antiferromagnet (KAFM) in the Zn$_x$Cu$_{5-x}$(VO$_4$)$_{2}$O$_2$CsCl (termed Zn$\mathbf {_x}$) series, namely Zn$_2$-averievite, Zn$\mathbf {_2}$ ($x=2$). Using magnetometry, synchrotron diffraction and neutron scattering we demonstrate an evolution of ground states with $x$ in Zn$\mathbf {_x}$; from long-range magnetic order in averievite ($x=0$), via a spin-glass-like ground state in Zn$\mathbf {_1}$, to a quantum spin liquid (QSL) in Zn$\mathbf {_2}$ for which inelastic neutron scattering reveals a gapless continuum of excitations. Similar to archetypal $S = 1/2$ KAFMs herbertsmithite and SrCr$_{8.19}$Ga$_{3.81}$O$_{19}$ (SCGO), the dynamic magnetic susceptibility of Zn$\mathbf {_2}$ shows scaling behavior consistent with proximity to a quantum critical point. The results demonstrate that the new Zn$\mathbf {_2}$ material is an excellent test bed for achieving the elusive goal of charge carrier doping in QSL states of $S = 1/2$ KAFMs, in-line with previous theoretical studies. △ Less

Submitted 26 June, 2023; originally announced June 2023.

Comments: 15 pages, 7 figures, Supplemental Material

The crystal and magnetic structure of cesium superoxide

Authors: R. A. Ewings, M. Reehuis, F. Orlandi, P. Manuel, D. D. Khalyavin, A. S. Gibbs, A. D. Fortes, A. Hoser, A. J. Princep, M. Jansen

Abstract: CsO2 is a member of the family of alkali superoxides (formula AO2 with A= Na, K, Rb and Cs) that exhibit magnetic behavior arising from open $p$-shell electrons residing on O2- molecules. We use neutron diffraction to solve the crystal and magnetic structures of CsO2, and observe a complex series of structures on cooling from room temperature to 1.6 K. These include an incommensurate modulation al… ▽ More CsO2 is a member of the family of alkali superoxides (formula AO2 with A= Na, K, Rb and Cs) that exhibit magnetic behavior arising from open $p$-shell electrons residing on O2- molecules. We use neutron diffraction to solve the crystal and magnetic structures of CsO2, and observe a complex series of structures on cooling from room temperature to 1.6 K. These include an incommensurate modulation along the a-axis of the structure at intermediate temperatures, which then locks into a commensurate modulation that doubles the unit cell compared to the previously supposed orthorhombic unit cell. In both incommensurate and commensurate phases our structural solution involves a staggering of the cesium ion positions along the b-axis, in contrast to studies of other alkali superoxides in which staggered tilts of the O2- dimers relative to the c-axis are seen. Below T ~ 10 K we observe magnetic Bragg reflections arising from an antiferromagnetically ordered structure with a wavevector of k = (0,0,0) (relative to the doubled crystallographic unit cell), with moments that point predominantly along the b-axis with a small component along the a-axis that hints at possible anisotropic exchange coupling (consistent with the crystal structure). Measurements of the magnetic Bragg reflections in an applied magnetic field suggest a spin-flop transition takes place between 2 T and 4 T in which moments likely flop to point along the crystallographic a-axis. Our measurements indicate that CsO2 is an interesting example of magnetic properties being inherently linked to the crystal structure, in that the staggered displacement of the cesium ions activates antisymmetric exchange which then permits the observed spin canting. △ Less

Submitted 30 June, 2023; v1 submitted 26 June, 2023; originally announced June 2023.

Comments: 13 pages, 11 figures, 6 tables

Numerical evaluation of singular integrals on non-disjoint self-similar fractal sets

Authors: Andrew Gibbs, David P. Hewett, Botond Major

Abstract: We consider the numerical evaluation of a class of double integrals with respect to a pair of self-similar measures over a self-similar fractal set (the attractor of an iterated function system), with a weakly singular integrand of logarithmic or algebraic type. In a recent paper [Gibbs, Hewett and Moiola, Numer. Alg., 2023] it was shown that when the fractal set is "disjoint" in a certain sense (… ▽ More We consider the numerical evaluation of a class of double integrals with respect to a pair of self-similar measures over a self-similar fractal set (the attractor of an iterated function system), with a weakly singular integrand of logarithmic or algebraic type. In a recent paper [Gibbs, Hewett and Moiola, Numer. Alg., 2023] it was shown that when the fractal set is "disjoint" in a certain sense (an example being the Cantor set), the self-similarity of the measures, combined with the homogeneity properties of the integrand, can be exploited to express the singular integral exactly in terms of regular integrals, which can be readily approximated numerically. In this paper we present a methodology for extending these results to cases where the fractal is non-disjoint but non-overlapping (in the sense that the open set condition holds). Our approach applies to many well-known examples including the Sierpinski triangle, the Vicsek fractal, the Sierpinski carpet, and the Koch snowflake. △ Less

Submitted 6 September, 2023; v1 submitted 23 March, 2023; originally announced March 2023.

MSC Class: 65D30; 28A80

Structure and dynamics of the negative thermal expansion material Cd(CN)$_2$ under pressure

Authors: Chloe S. Coates, Mia Baise, Johanthan M. Bulled, Emma H. Wolpert, Joshua W. Makepeace, Helen C. Walker, Alexandra S. Gibbs, A. Dominic Fortes, Ben Slater, Andrew L. Goodwin

Abstract: We use a combination of variable-temperature / variable-pressure neutron powder diffraction, variable-pressure inelastic neutron scattering, and quantum chemical calculations to interrogate the behaviour of the negative thermal expansion (NTE) material $^{114}$Cd(CN)$_2$ under hydrostatic pressure. We determine the equation of state of the ambient-pressure phase, and discover the so-called `warm h… ▽ More We use a combination of variable-temperature / variable-pressure neutron powder diffraction, variable-pressure inelastic neutron scattering, and quantum chemical calculations to interrogate the behaviour of the negative thermal expansion (NTE) material $^{114}$Cd(CN)$_2$ under hydrostatic pressure. We determine the equation of state of the ambient-pressure phase, and discover the so-called `warm hardening' effect whereby the material becomes elastically stiffer as it is heated. We also identify a number of high-pressure phases, and map out the phase behaviour of Cd(CN)$_2$ over the range $0\leq p\leq0.5$\,GPa, $100\leq T\leq300$\,K. As expected for an NTE material, the low-energy phonon frequencies are found to soften under pressure, and we determine an effective Gr{ü}neisen parameter for these modes. Finally, we show that the elastic behaviour of Cd(CN)$_2$ is sensitive to the local Cd coordination environment, which suggests an interplay between short- (phononic) and long-timescale (cyanide flips) fluctuations in Cd(CN)$_2$. △ Less

Submitted 20 February, 2023; originally announced February 2023.

EDEN Survey: Small Transiting Planet Detection Limits and Constraints on the Occurrence Rates for Late M Dwarfs within 15 pc

Authors: Jeremy Dietrich, Dániel Apai, Martin Schlecker, Kevin K. Hardegree-Ullman, Benjamin V. Rackham, Nicolas Kurtovic, Karan Molaverdikhani, Paul Gabor, Thomas Henning, Wen-Ping Chen, Luigi Mancini, Alex Bixel, Aidan Gibbs, Richard P. Boyle, Samantha Brown-Sevilla, Remo Burn, Timmy N. Delage, Lizxandra Flores-Rivera, Riccardo Franceschi, Gabriele Pichierri, Sofia Savvidou, Jonas Syed, Ivan Bruni, Wing-Huen Ip, Chow-Choong Ngeow , et al. (7 additional authors not shown)

Abstract: Earth-sized exoplanets that transit nearby, late spectral type red dwarfs will be prime targets for atmospheric characterization in the coming decade. Such systems, however, are difficult to find via wide-field transit surveys like Kepler or TESS. Consequently, the presence of such transiting planets is unexplored and the occurrence rates of short-period Earth-sized planets around late M dwarfs re… ▽ More Earth-sized exoplanets that transit nearby, late spectral type red dwarfs will be prime targets for atmospheric characterization in the coming decade. Such systems, however, are difficult to find via wide-field transit surveys like Kepler or TESS. Consequently, the presence of such transiting planets is unexplored and the occurrence rates of short-period Earth-sized planets around late M dwarfs remain poorly constrained. Here, we present the deepest photometric monitoring campaign of 22 nearby late M dwarf stars, using data from over 500 nights on seven 1-2 meter class telescopes. Our survey includes all known single quiescent northern late M dwarfs within 15 pc. We use transit-injection-and-recovery tests to quantify the completeness of our survey, successfully identify most ($>80\%$) transiting short-period (0.5-1 d) super-Earths ($R > 1.9 R_\oplus$), and are sensitive ($\sim50\%$) to transiting Earth-sized planets ($1.0-1.2 R_\oplus$). Our high sensitivity to transits with a near-zero false positive rate demonstrates an efficient survey strategy. Our survey does not yield a transiting planet detection, yet it provides the most sensitive upper limits on transiting planets orbiting our target stars. Finally, we explore multiple hypotheses about the occurrence rates of short-period planets (from Earth-sized planets to giant planets) around late M dwarfs. We show, for example, that giant planets at short periods ($<1$ day) are uncommon around our target stars. Our dataset provides some insight into occurrence rates of short-period planets around TRAPPIST-1-like stars, and our results can help test planetary formation and system evolution models, as well as guide future observations of nearby late M dwarfs. △ Less

Submitted 8 February, 2023; originally announced February 2023.

Comments: 27 pages, 11 figures

Surveying Nearby Brown Dwarfs with HGCA: Direct Imaging Discovery of a Faint, High-Mass Brown Dwarf Orbiting HD 176535 A

Authors: Yiting Li, Timothy D. Brandt, G. Mirek Brandt, Qier An, Kyle Franson, Trent J. Dupuy, Minghan Chen, Rachel Bowens-Rubin, Briley L. Lewis, Brendan P. Bowler, Aidan Gibbs, Rocio Kiman, Jacqueline Faherty, Thayne Currie, Rebecca Jensen-Clem, Hengyue Zhang Ezequiel Contreras-Martinez, Michael P. Fitzgerald, Benjamin A. Mazin, Maxwell Millar-Blanchaer

Abstract: Brown dwarfs with well-measured masses, ages and luminosities provide direct benchmark tests of substellar formation and evolutionary models. We report the first results from a direct imaging survey aiming to find and characterize substellar companions to nearby accelerating stars with the assistance of the Hipparcos-Gaia Catalog of Accelerations (HGCA). In this paper, we present a joint high-cont… ▽ More Brown dwarfs with well-measured masses, ages and luminosities provide direct benchmark tests of substellar formation and evolutionary models. We report the first results from a direct imaging survey aiming to find and characterize substellar companions to nearby accelerating stars with the assistance of the Hipparcos-Gaia Catalog of Accelerations (HGCA). In this paper, we present a joint high-contrast imaging and astrometric discovery of a substellar companion to HD 176535 A, a K3.5V main-sequence star aged approximately $3.59_{-1.15}^{+0.87}$ Gyrs at a distance of $36.99 \pm 0.03$ pc. In advance of our high-contrast imaging observations, we combined precision HARPS RVs and HGCA astrometry to predict the potential companion's location and mass. We thereafter acquired two nights of KeckAO/NIRC2 direct imaging observations in the $L'$ band, which revealed a companion with a contrast of $ΔL'_p = 9.20\pm0.06$ mag at a projected separation of $\approx$0.$\!\!''35$ ($\approx$13 AU) from the host star. We revise our orbital fit by incorporating our dual-epoch relative astrometry using the open-source MCMC orbit fitting code $\tt orvara$. HD 176535 B is a new benchmark dwarf useful for constraining the evolutionary and atmospheric models of high-mass brown dwarfs. We found a luminosity of $\rm log(L_{bol}/L_{\odot}) = -5.26\pm0.07$ and a model-dependent effective temperature of $980 \pm 35$ K for HD 176535 B. Our dynamical mass suggests that some substellar evolutionary models may be underestimating luminosity for high-mass T dwarfs. Given its angular separation and luminosity, HD 176535 B would make a promising candidate for Aperture Masking Interferometry with JWST and GRAVITY/KPIC, and further spectroscopic characterization with instruments like the CHARIS/SCExAO/Subaru integral field spectrograph. △ Less

Submitted 16 May, 2023; v1 submitted 25 January, 2023; originally announced January 2023.

A Hausdorff-measure boundary element method for acoustic scattering by fractal screens

Authors: António M. Caetano, Simon N. Chandler-Wilde, Andrew Gibbs, David P. Hewett, Andrea Moiola

Abstract: Sound-soft fractal screens can scatter acoustic waves even when they have zero surface measure. To solve such scattering problems we make what appears to be the first application of the boundary element method (BEM) where each BEM basis function is supported in a fractal set, and the integration involved in the formation of the BEM matrix is with respect to a non-integer order Hausdorff measure ra… ▽ More Sound-soft fractal screens can scatter acoustic waves even when they have zero surface measure. To solve such scattering problems we make what appears to be the first application of the boundary element method (BEM) where each BEM basis function is supported in a fractal set, and the integration involved in the formation of the BEM matrix is with respect to a non-integer order Hausdorff measure rather than the usual (Lebesgue) surface measure. Using recent results on function spaces on fractals, we prove convergence of the Galerkin formulation of this ``Hausdorff BEM'' for acoustic scattering in $\mathbb{R}^{n+1}$ ($n=1,2$) when the scatterer, assumed to be a compact subset of $\mathbb{R}^n\times\{0\}$, is a $d$-set for some $d\in (n-1,n]$, so that, in particular, the scatterer has Hausdorff dimension $d$. For a class of fractals that are attractors of iterated function systems, we prove convergence rates for the Hausdorff BEM and superconvergence for smooth antilinear functionals, under certain natural regularity assumptions on the solution of the underlying boundary integral equation. We also propose numerical quadrature routines for the implementation of our Hausdorff BEM, along with a fully discrete convergence analysis, via numerical (Hausdorff measure) integration estimates and inverse estimates on fractals, estimating the discrete condition numbers. Finally, we show numerical experiments that support the sharpness of our theoretical results, and our solution regularity assumptions, including results for scattering in $\mathbb{R}^2$ by Cantor sets, and in $\mathbb{R}^3$ by Cantor dusts. △ Less

Submitted 7 February, 2024; v1 submitted 13 December, 2022; originally announced December 2022.

Comments: 61 pages

MSC Class: 28A80; 65R20; 65N38; 78A45; 78M15

Journal ref: Numer. Math., 156, 2024, 463-532

Creating and controlling Dirac fermions, Weyl fermions, and nodal lines in the magnetic antiperovskite Eu$_3$PbO

Authors: Moritz M. Hirschmann, Alexandra S. Gibbs, Fabio Orlandi, Dmitry Khalyavin, Pascal Manuel, Vahideh Abdolazimi, Alexander Yaresko, Jürgen Nuss, H. Takagi, Andreas P. Schnyder, Andreas W. Rost

Abstract: The band topology of magnetic semimetals is of interest both from the fundamental science point of view and with respect to potential spintronics and memory applications. Unfortunately, only a handful of suitable topological semimetals with magnetic order have been discovered so far. One such family that hosts these characteristics is the antiperovskites, A$_3$BO, a family of 3D Dirac semimetals.… ▽ More The band topology of magnetic semimetals is of interest both from the fundamental science point of view and with respect to potential spintronics and memory applications. Unfortunately, only a handful of suitable topological semimetals with magnetic order have been discovered so far. One such family that hosts these characteristics is the antiperovskites, A$_3$BO, a family of 3D Dirac semimetals. The A=Eu$^{2+}$ compounds magnetically order with multiple phases as a function of applied magnetic field. Here, by combining band structure calculations with neutron diffraction and magnetic measurements, we establish the antiperovskite Eu$_3$PbO as a new topological magnetic semimetal. This topological material exhibits a multitude of different topological phases with ordered Eu moments which can be easily controlled by an external magnetic field. The topological phase diagram of Eu$_3$PbO includes an antiferromagnetic Dirac phase, as well as ferro- and ferrimagnetic phases with both Weyl points and nodal lines. For each of these phases, we determine the bulk band dispersions, the surface states, and the topological invariants by means of $\textit{ab-initio}$ and tight-binding calculations. Our discovery of these topological phases introduces Eu$_3$PbO as a new platform to study and manipulate the interplay of band topology, magnetism, and transport. △ Less

Submitted 25 July, 2022; originally announced July 2022.

Limits on the Auroral Generation of H$_3^+$ in Brown Dwarf and Extrasolar Giant Planet Atmospheres with Keck/NIRSPEC

Authors: Aidan Gibbs, Michael Fitzgerald

Abstract: The molecular ion H$_3^+$ is a potentially powerful tracer of the ionospheres and thermal structures of Jovian planets, but has never been detected in a planetary mass object outside of the solar system. Models predict that H$_3^+$ emission driven by EUV flux and solar wind on hot Jupiters, or by powerful aurorae on brown dwarfs, will be between $10^2$ and $10^5\times$ more intense than that of Ju… ▽ More The molecular ion H$_3^+$ is a potentially powerful tracer of the ionospheres and thermal structures of Jovian planets, but has never been detected in a planetary mass object outside of the solar system. Models predict that H$_3^+$ emission driven by EUV flux and solar wind on hot Jupiters, or by powerful aurorae on brown dwarfs, will be between $10^2$ and $10^5\times$ more intense than that of Jupiter. If optimal conditions for the production of emission do exist, the emission may be detectable by current ground-based instruments or in the near future. We present the first search for H$_3^+$ line emission in brown dwarfs with Keck/NIRSPEC $L^\prime$ high-resolution spectroscopy. Additionally, we survey stars hosting giant planets at semi-major axes near $0.1-0.2$ au, which models suggest may be the best planetary targets. No candidate H$_3^+$ emission is found. The limits we place on the emission of H$_3^+$ from brown dwarfs indicates that auroral generation of H$_3^+$ in these environments likely does not linearly scale from the processes found on Jupiter, plausibly due to deeper atmospheric penetration by precipitating auroral electrons. Detection of H$_3^+$ emission in brown dwarfs may be possible with the James Webb Space Telescope (JWST), or future thirty-meter class telescopes. △ Less

Submitted 7 June, 2022; originally announced June 2022.

Comments: 22 pages, 12 figures

Competing spin-orbital singlet states in the 4$d^4$ honeycomb ruthenate Ag$_3$LiRu$_2$O$_6$

Authors: T. Takayama, M. Blankenhorn, J. Bertinshaw, D. Haskel, N. A. Bogdanov, K. Kitagawa, A. N. Yaresko, A. Krajewska, S. Bette, G. McNally, A. S. Gibbs, Y. Matsumoto, D. P. Sari, I. Watanabe, G. Fabbris, W. Bi, T. I. Larkin, K. S. Rabinovich, A. V. Boris, H. Ishii, H. Yamaoka, T. Irifune, R. Bewley, C. J. Ridley, C. L. Bull , et al. (3 additional authors not shown)

Abstract: When spin-orbit-entangled $d$-electrons reside on a honeycomb lattice, rich quantum states are anticipated to emerge, as exemplified by the $d^5$ Kitaev materials. Distinct yet equally intriguing physics may be realized with a $d$-electron count other than $d^5$. We found that the layered ruthenate Ag$_3$LiRu$_2$O$_6$ with $d^4$ Ru$^{4+}$ ions at ambient pressure forms a honeycomb lattice of spin-… ▽ More When spin-orbit-entangled $d$-electrons reside on a honeycomb lattice, rich quantum states are anticipated to emerge, as exemplified by the $d^5$ Kitaev materials. Distinct yet equally intriguing physics may be realized with a $d$-electron count other than $d^5$. We found that the layered ruthenate Ag$_3$LiRu$_2$O$_6$ with $d^4$ Ru$^{4+}$ ions at ambient pressure forms a honeycomb lattice of spin-orbit-entangled singlets, which is a playground for frustrated excitonic magnetism. Under pressure, the singlet state does not develop the expected excitonic magnetism but experiences two successive transitions to other nonmagnetic phases, first to an intermediate phase with moderate distortion of honeycomb lattice, and eventually to a high-pressure phase with very short Ru-Ru dimer bonds. While the strong dimerization in the high-pressure phase originates from a molecular orbital formation as in the sister compound Li$_2$RuO$_3$, the intermediate phase represents a spin-orbit-coupled $J$-dimer state which is stabilized by the admixture of upper-lying $J_{\rm eff} = 1$-derived states. We argue that the $J$-dimer state is induced by a pseudo-Jahn-Teller effect associated with the low-lying spin-orbital excited states and is unique to spin-orbit-entangled $d^4$ systems. The discovery of competing singlet phases demonstrates rich spin-orbital physics of $d^4$ honeycomb compounds and paves the way for realization of unconventional magnetism. △ Less

Submitted 24 May, 2022; originally announced May 2022.

Comments: 12 pages, 5 figures

Report number: Phys. Rev. Res. 4, 043079 (2022)

Numerical Quadrature for Singular Integrals on Fractals

Authors: A. Gibbs, D. P. Hewett, A. Moiola

Abstract: We present and analyse numerical quadrature rules for evaluating regular and singular integrals on self-similar fractal sets. The integration domain $\mathbb{R}^n$ is assumed to be the compact attractor of an iterated function system of contracting similarities satisfying the open set condition. Integration is with respect to any ``invariant'' (also known as ``balanced'' or ``self-similar'') measu… ▽ More We present and analyse numerical quadrature rules for evaluating regular and singular integrals on self-similar fractal sets. The integration domain $\mathbb{R}^n$ is assumed to be the compact attractor of an iterated function system of contracting similarities satisfying the open set condition. Integration is with respect to any ``invariant'' (also known as ``balanced'' or ``self-similar'') measure supported on $Γ$, including in particular the Hausdorff measure $\mathcal{H}^d$ restricted to $Γ$, where $d$ is the Hausdorff dimension of $Γ$. Both single and double integrals are considered. Our focus is on composite quadrature rules in which integrals over $Γ$ are decomposed into sums of integrals over suitable partitions of $Γ$ into self-similar subsets. For certain singular integrands of logarithmic or algebraic type we show how in the context of such a partitioning the invariance property of the measure can be exploited to express the singular integral exactly in terms of regular integrals. For the evaluation of these regular integrals we adopt a composite barycentre rule, which for sufficiently regular integrands exhibits second-order convergence with respect to the maximum diameter of the subsets. As an application we show how this approach, combined with a singularity-subtraction technique, can be used to accurately evaluate the singular double integrals that arise in Hausdorff-measure Galerkin boundary element methods for acoustic wave scattering by fractal screens. △ Less

Submitted 14 July, 2022; v1 submitted 22 December, 2021; originally announced December 2021.

MSC Class: 65D32 (Primary); 28A78; 28A80 (Secondary)

Journal ref: Numer. Algorithms, 92, 2023, pp. 2071-2124

A Systematic Exploration of Kilonova Candidates from Neutron Star Mergers During the Third Gravitational Wave Observing Run

Authors: J. Rastinejad, K. Paterson, W. Fong, D. J. Sand, M. J. Lundquist, G. Hosseinzadeh, E. Christensen, P. N. Daly, A. R. Gibbs, S. Hall, F. Shelly, S. Yang

Abstract: We present a comprehensive analysis of 653 optical candidate counterparts reported during the third gravitational wave (GW) observing run. Our sample concentrates on candidates from the 15 events (published in GWTC-2, GWTC-3 or not retracted on GraceDB) that had a >1% chance of including a neutron star in order to assess their viability as true kilonovae. In particular, we leverage tools available… ▽ More We present a comprehensive analysis of 653 optical candidate counterparts reported during the third gravitational wave (GW) observing run. Our sample concentrates on candidates from the 15 events (published in GWTC-2, GWTC-3 or not retracted on GraceDB) that had a >1% chance of including a neutron star in order to assess their viability as true kilonovae. In particular, we leverage tools available in real time, including pre-merger detections and cross-matching with catalogs (i.e. point source, variable star, quasar and host galaxy redshift datasets), to eliminate 65% of candidates in our sample. We further employ spectroscopic classifications, late-time detections and light curve behavior analyses, and conclude that 66 candidates remain viable kilonovae. These candidates lack sufficient information to determine their classifications, and the majority would require luminosities greater than that of AT2017gfo. Pre-merger detections in public photometric survey data and comparison of catalogued host galaxy redshifts with the GW event distances are critical to incorporate into vetting procedures, as these tools eliminated >20% and >30% of candidates, respectively. We expect that such tools which leverage archival information will significantly reduce the strain on spectroscopic and photometric follow-up resources in future observing runs. Finally, we discuss the critical role prompt updates from GW astronomers to the EM community play in reducing the number of candidates requiring vetting. △ Less

Submitted 1 February, 2022; v1 submitted 17 December, 2021; originally announced December 2021.

Comments: Accepted to ApJ

Formation of orbital molecules on a pyrochlore lattice induced by A-O bond covalency

Authors: A. Krajewska, A. Yaresko, J. Nuss, A. S. Gibbs, S. Bette, M. Blankenhorn, R. E. Dinnebier, D. P. Sari, I. Watanabe, J. Bertinshaw, H. Gretarsson, K. Ishii, D. Matsumura, T. Tsuji, M. Isobe, B. Keimer, H. Takagi, T. Takayama

Abstract: The pyrochlore ruthenate In$_2$Ru$_2$O$_7$ displays a subtle competition between spin-orbital entanglement and molecular orbital formation. At room temperature, a spin-orbit-entangled singlet state was identified. With decreasing temperature, In$_2$Ru$_2$O$_7$ undergoes multiple structural transitions and eventually forms a nonmagnetic ground state with semi-isolated Ru$_2$O units on the pyrochlor… ▽ More The pyrochlore ruthenate In$_2$Ru$_2$O$_7$ displays a subtle competition between spin-orbital entanglement and molecular orbital formation. At room temperature, a spin-orbit-entangled singlet state was identified. With decreasing temperature, In$_2$Ru$_2$O$_7$ undergoes multiple structural transitions and eventually forms a nonmagnetic ground state with semi-isolated Ru$_2$O units on the pyrochlore lattice. The dominant hopping through the Ru-O-Ru linkage leads to molecular orbital formation within the Ru$_2$O units. This molecular orbital formation is unique in that it involves the O$^{2-}$ anions, unlike the transition-metal dimers observed in systems with edge-sharing octahedra. We argue that the covalent character of In-O bonds plays a pivotal role in the structural transitions and molecular orbital formation and such bonding character of "$A$-site" ions is an important ingredient for electronic phase competition in complex transition metal oxides. △ Less

Submitted 9 November, 2021; originally announced November 2021.

Site-selective $d^{10}$/$d^0$ substitution in a $S = 1/2$ spin ladder Ba$_2$CuTe$_{1-x}$W$_x$O$_6$ ($0 \leq x \leq 0.3$)

Authors: Charlotte E. Pughe, Otto H. J. Mustonen, Alexandra S. Gibbs, Martin Etter, Cheng Liu, Siân E. Dutton, Aidan Friskney, Neil C. Hyatt, Gavin B. G. Stenning, Heather M. Mutch, Fiona C. Coomer, Edmund J. Cussen

Abstract: Remarkably, doping isovalent $d^{10}$ and $d^0$ cations onto the $B$'' site in $A_2B$'$B$''O$_6$ double perovskites has the power to direct the magnetic interactions between magnetic $B$' cations. This is due to changes in orbital hybridization, which favors different superexchange pathways, and leads to the formation of alternative magnetic structures depending on whether $B$'' is $d^{10}$ or… ▽ More Remarkably, doping isovalent $d^{10}$ and $d^0$ cations onto the $B$'' site in $A_2B$'$B$''O$_6$ double perovskites has the power to direct the magnetic interactions between magnetic $B$' cations. This is due to changes in orbital hybridization, which favors different superexchange pathways, and leads to the formation of alternative magnetic structures depending on whether $B$'' is $d^{10}$ or $d^0$. Furthermore, the competition generated by introducing mixtures of $d^{10}$ and $d^0$ cations can drive the material into the realms of exotic quantum magnetism. Here, a W$^{6+}$ $d^0$ dopant was introduced to a $d^{10}$ hexagonal perovskite Ba$_2$CuTeO$_6$, which possesses a spin ladder geometry of Cu$^{2+}$ cations, creating a Ba$_2$CuTe$_{1-x}$W$_x$O$_6$ solid solution ($x$ = 0 - 0.3). Neutron and synchrotron X-ray diffraction show that W$^{6+}$ is almost exclusively substituted for Te$^{6+}$ on the corner-sharing site within the spin ladder, in preference to the face-sharing site between ladders. This means the intra-ladder interactions are selectively tuned by the $d^0$ cations. Bulk magnetic measurements suggest this suppresses magnetic ordering in a similar manner to that observed for the spin-liquid like material Sr$_2$CuTe$_{1-x}$W$_x$O$_6$. This further demonstrates the utility of $d^{10}$ and $d^0$ dopants as a tool for tuning magnetic ground states in a wide range of perovskites and perovskite-derived structures. △ Less

Submitted 5 August, 2021; originally announced August 2021.

Comments: 26+24 pages, 7+20 figures

Journal ref: Inorganic Chemistry 61, 4033-4045 (2022)

A large sub-Neptune transiting the thick-disk M4V TOI-2406

Authors: R. D. Wells, B. V. Rackham, N. Schanche, R. Petrucci, Y. Gomez Maqueo Chew, B. -O. Demory, A. J. Burgasser, R. Burn, F. J. Pozuelos, M. N. Gunther, L. Sabin, U. Schroffenegger, M. A. Gomez-Munoz, K. G. Stassun, V. Van Grootel, S. B. Howell, D. Sebastian, A. H. M. J. Triaud, D. Apai, I. Plauchu-Frayn, C. A. Guerrero, P. F. Guillen, A. Landa, G. Melgoza, F. Montalvo , et al. (49 additional authors not shown)

Abstract: Large sub-Neptunes are uncommon around the coolest stars in the Galaxy and are rarer still around those that are metal-poor. However, owing to the large planet-to-star radius ratio, these planets are highly suitable for atmospheric study via transmission spectroscopy in the infrared, such as with JWST. Here we report the discovery and validation of a sub-Neptune orbiting the thick-disk, mid-M dwar… ▽ More Large sub-Neptunes are uncommon around the coolest stars in the Galaxy and are rarer still around those that are metal-poor. However, owing to the large planet-to-star radius ratio, these planets are highly suitable for atmospheric study via transmission spectroscopy in the infrared, such as with JWST. Here we report the discovery and validation of a sub-Neptune orbiting the thick-disk, mid-M dwarf star TOI-2406. We first infer properties of the host star by analysing the star's near-infrared spectrum, spectral energy distribution, and Gaia parallax. We use multi-band photometry to confirm that the transit event is on-target and achromatic, and we statistically validate the TESS signal as a transiting exoplanet. We then determine physical properties of the planet through global transit modelling of the TESS and ground-based time-series data. We determine the host to be a metal-poor M4V star, located at a distance of 56 pc, with a sub-solar metallicity $(\mathrm{[Fe/H] = -0.38 \pm 0.07})$, and a member of the thick disk. The planet is a relatively large sub-Neptune for the M-dwarf planet population, with $\mathrm{R_p = 2.94 \pm 0.17} \mathrm{R_\oplus}$ and $\mathrm{P = 3.077}$ d, producing transits of 2% depth. We note the orbit has a non-zero eccentricity to 3$\mathrmσ$, prompting questions about the dynamical history of the system. This system is an interesting outcome of planet formation and presents a benchmark for large-planet formation around metal-poor, low-mass stars. The system warrants further study, in particular radial velocity follow-up to determine the planet mass and constrain possible bound companions. Furthermore, TOI-2406 b is a good target for future atmospheric study through transmission spectroscopy, particularly in the category of warm sub-Neptunes. △ Less

Submitted 29 July, 2021; originally announced July 2021.

Comments: Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 653, A97 (2021)

Discovery of Superconductivity in (Ba,K)SbO$_{3}$

Authors: Minu Kim, Graham M. McNally, Hun-Ho Kim, Mohamed Oudah, Alexandra Gibbs, Pascal Manuel, Robert Green, Tomohiro Takayama, Alexander Yaresko, Ulrich Wedig, Masahiko Isobe, Reinhard K. Kremer, D. A. Bonn, Bernhard Keimer, Hidenori Takagi

Abstract: Superconducting bismuthates (Ba,K)BiO$_{3}$ (BKBO) constitute an interesting class of superconductors in that superconductivity with a remarkably high $T_\mathrm{c}$ of 30 K arises in proximity to charge density wave (CDW) order. Prior understanding on the driving mechanism of the CDW and superconductivity emphasizes the role of either bismuth (negative $U$ model) or oxygen ions (ligand hole model… ▽ More Superconducting bismuthates (Ba,K)BiO$_{3}$ (BKBO) constitute an interesting class of superconductors in that superconductivity with a remarkably high $T_\mathrm{c}$ of 30 K arises in proximity to charge density wave (CDW) order. Prior understanding on the driving mechanism of the CDW and superconductivity emphasizes the role of either bismuth (negative $U$ model) or oxygen ions (ligand hole model). While holes in BKBO presumably reside on oxygen owing to their negative charge transfer energy, so far there has been no other comparative material studied. Here, we introduce (Ba,K)SbO$_{3}$ (BKSO) in which the Sb 5$s$ orbital energy is higher than that of the Bi 6$s$ orbitals enabling tuning of the charge transfer energy from negative to slightly positive. The parent compound BaSbO$_{3-δ}$ shows a larger CDW gap compared to the undoped bismuthate BaBiO$_{3}$. As the CDW order is suppressed via potassium substitution up to 65 %, superconductivity emerges, rising up to $T_\mathrm{c}$ = 15 K. This value is lower than the maximum $T_\mathrm{c}$ of BKBO, but higher by more than a factor of two at comparable potassium concentrations. The discovery of an enhanced CDW gap and superconductivity in BKSO indicates that the sign of the charge transfer energy may not be crucial, but instead strong metal-oxygen covalency plays the essential role in constituting a CDW and high-$T_\mathrm{c}$ superconductivity in the main-group perovskite oxides. △ Less

Submitted 14 July, 2021; originally announced July 2021.

Journal ref: Nat. Mater. 21, 627-633 (2022)

The impact and recovery of asteroid 2018 LA

Authors: Peter Jenniskens, Mohutsiwa Gabadirwe, Qing-Zhu Yin, Alexander Proyer, Oliver Moses, Tomas Kohout, Fulvio Franchi, Roger L. Gibson, Richard Kowalski, Eric J. Christensen, Alex R. Gibbs, Aren Heinze, Larry Denneau, Davide Farnocchia, Paul W. Chodas, William Gray, Marco Micheli, Nick Moskovitz, Christopher A. Onken, Christian Wolf, Hadrien A. R. Devillepoix, Quanzhi Ye, Darrel K. Robertson, Peter Brown, Esko Lyytinen , et al. (41 additional authors not shown)

Abstract: The June 2, 2018, impact of asteroid 2018 LA over Botswana is only the second asteroid detected in space prior to impacting over land. Here, we report on the successful recovery of meteorites. Additional astrometric data refine the approach orbit and define the spin period and shape of the asteroid. Video observations of the fireball constrain the asteroid's position in its orbit and were used to… ▽ More The June 2, 2018, impact of asteroid 2018 LA over Botswana is only the second asteroid detected in space prior to impacting over land. Here, we report on the successful recovery of meteorites. Additional astrometric data refine the approach orbit and define the spin period and shape of the asteroid. Video observations of the fireball constrain the asteroid's position in its orbit and were used to triangulate the location of the fireball's main flare over the Central Kalahari Game Reserve. 23 meteorites were recovered. A consortium study of eight of these classifies Motopi Pan as a HED polymict breccia derived from howardite, cumulate and basaltic eucrite, and diogenite lithologies. Before impact, 2018 LA was a solid rock of about 156 cm diameter with high bulk density about 2.85 g/cm3, a relatively low albedo pV about 0.25, no significant opposition effect on the asteroid brightness, and an impact kinetic energy of about 0.2 kt. The orbit of 2018 LA is consistent with an origin at Vesta (or its Vestoids) and delivery into an Earth-impacting orbit via the nu_6 resonance. The impact that ejected 2018 LA in an orbit towards Earth occurred 22.8 +/- 3.8 Ma ago. Zircons record a concordant U-Pb age of 4563 +/- 11 Ma and a consistent 207Pb/206Pb age of 4563 +/- 6 Ma. A much younger Pb-Pb phosphate resetting age of 4234 +/- 41 Ma was found. From this impact chronology, we discuss what is the possible source crater of Motopi Pan and the age of Vesta's Veneneia impact basin. △ Less

Submitted 12 May, 2021; originally announced May 2021.

Comments: Meteoritics & Planetary Science (2021)

High sensitivity heat capacity measurements on Sr2RuO4 under uniaxial pressure

Authors: You-Sheng Li, Naoki Kikugawa, Dmitry A. Sokolov, Fabian Jerzembeck, Alexandra S. Gibbs, Yoshiteru Maeno, Clifford W. Hicks, Jörg Schmalian, Michael Nicklas, Andrew P. Mackenzie

Abstract: A key question regarding the unconventional superconductivity of Sr$_2$RuO$_4$ remains whether the order parameter is single- or two-component. Under a hypothesis of two-component superconductivity, uniaxial pressure is expected to lift their degeneracy, resulting in a split transition. The most direct and fundamental probe of a split transition is heat capacity. Here, we report measurement of hea… ▽ More A key question regarding the unconventional superconductivity of Sr$_2$RuO$_4$ remains whether the order parameter is single- or two-component. Under a hypothesis of two-component superconductivity, uniaxial pressure is expected to lift their degeneracy, resulting in a split transition. The most direct and fundamental probe of a split transition is heat capacity. Here, we report measurement of heat capacity of samples subject to large and highly homogeneous uniaxial pressure. We place an upper limit on the heat-capacity signature of any second transition of a few per cent of that of the primary superconducting transition. The normalized jump in heat capacity, $ΔC/C$, grows smoothly as a function of uniaxial pressure, favouring order parameters which are allowed to maximize in the same part of the Brillouin zone as the well-studied van Hove singularity. Thanks to the high precision of our measurements, these findings place stringent constraints on theories of the superconductivity of Sr$_2$RuO$_4$. △ Less

Submitted 3 March, 2021; originally announced March 2021.

Comments: This manuscript is an extended version of the scientific part of arXiv:1906.07597, which is being split into two separate papers. An extended discussion on the experimental methods can be found at arXiv:2009.03125. It includes supplemental material

Journal ref: PNAS 118, e2020492118 (2021)

Searches after Gravitational Waves Using ARizona Observatories (SAGUARO): Observations and Analysis from Advanced LIGO/Virgo's Third Observing Run

Authors: K. Paterson, M. J. Lundquist, J. C. Rastinejad, W. Fong, D. J. Sand, J. E. Andrews, R. C. Amaro, O. Eskandari, S. Wyatt, P. N. Daly, H. Bradley, S. Zhou-Wright, S. Valenti, S. Yang, E. Christensen, A. R. Gibbs, F. Shelly, C. Bilinski, L. Chomiuk, A. Corsi, M. R. Drout, R. J. Foley, P. Gabor, P. Garnavich, C. J. Grier , et al. (12 additional authors not shown)

Abstract: With the conclusion of the third observing run for Advanced LIGO/Virgo (O3), we present a detailed analysis of both triggered and serendipitous observations of 17 gravitational wave (GW) events (7 triggered and 10 purely serendipitous) from the Searches After Gravitational-waves Using ARizona Observatories (SAGUARO) program. We searched a total of 4935 deg$^2$ down to a median 5$σ$ transient detec… ▽ More With the conclusion of the third observing run for Advanced LIGO/Virgo (O3), we present a detailed analysis of both triggered and serendipitous observations of 17 gravitational wave (GW) events (7 triggered and 10 purely serendipitous) from the Searches After Gravitational-waves Using ARizona Observatories (SAGUARO) program. We searched a total of 4935 deg$^2$ down to a median 5$σ$ transient detection depth of 21.1 AB mag using the Mt Lemmon 1.5 m telescope, the discovery engine for SAGUARO. In addition to triggered events within 24~hours, our transient search encompassed a time interval following GW events of $<120$~hrs, providing observations on $\sim$ 1/2 of the events accessible to the Mt Lemmon 1.5 m telescope. We covered 2.1--86\% of the LVC total probability ($P_{\rm total}$) for individual events, with a median $P_{\rm total} \approx 8\%$ within $<120$~hours. Following improvements to our pipeline and the addition of serendipitous observations, we find a total of 7 new optical candidates across 5 GW events which we are unable to rule out after searching for additional information and comparing to kilonova models. Using both publicly available and our own late-time data, we investigated a total of 252 optical candidates for these 17 events, finding only 65\% were followed up in some capacity by the community. Of the total 252 candidates, we are able to rule out an additional 12 previously reported counterpart candidates. In light of these results, we discuss lessons learned from the SAGUARO GW counterpart search. We discuss how community coordination of observations and candidate follow-up, as well as the role of archival data, are crucial to improving the efficiency of follow-up efforts and preventing unnecessary duplication of effort with limited EM resources. △ Less

Submitted 21 December, 2020; originally announced December 2020.

Comments: 44 pages, 21 figures, 2 tables, submitted to ApJ

Computational Skills by Stealth in Secondary School Data Science

Authors: Wesley Burr, Fanny Chevalier, Christopher Collins, Alison L Gibbs, Raymond Ng, Chris Wild

Abstract: The unprecedented growth in the availability of data of all types and qualities and the emergence of the field of data science has provided an impetus to finally realizing the implementation of the full breadth of the Nolan and Temple Lang proposed integration of computing concepts into statistics curricula at all levels in statistics and new data science programs and courses. Moreover, data scien… ▽ More The unprecedented growth in the availability of data of all types and qualities and the emergence of the field of data science has provided an impetus to finally realizing the implementation of the full breadth of the Nolan and Temple Lang proposed integration of computing concepts into statistics curricula at all levels in statistics and new data science programs and courses. Moreover, data science, implemented carefully, opens accessible pathways to stem for students for whom neither mathematics nor computer science are natural affinities, and who would traditionally be excluded. We discuss a proposal for the stealth development of computational skills in students' first exposure to data science through careful, scaffolded exposure to computation and its power. The intent of this approach is to support students, regardless of interest and self-efficacy in coding, in becoming data-driven learners, who are capable of asking complex questions about the world around them, and then answering those questions through the use of data-driven inquiry. This discussion is presented in the context of the International Data Science in Schools Project which recently published computer science and statistics consensus curriculum frameworks for a two-year secondary school data science program, designed to make data science accessible to all. △ Less

Submitted 8 October, 2020; originally announced October 2020.

Comments: 38 pages, 8 figures

Acoustic scattering by impedance screens/cracks with fractal boundary: well-posedness analysis and boundary element approximation

Authors: J. Bannister, A. Gibbs, D. P. Hewett

Abstract: We study time-harmonic scattering in $\mathbb{R}^n$ ($n=2,3$) by a planar screen (a "crack" in the context of linear elasticity), assumed to be a non-empty bounded relatively open subset $Γ$ of the hyperplane $\mathbb{R}^{n-1}\times \{0\}$, on which impedance (Robin) boundary conditions are imposed. In contrast to previous studies, $Γ$ can have arbitrarily rough (possibly fractal) boundary. To obt… ▽ More We study time-harmonic scattering in $\mathbb{R}^n$ ($n=2,3$) by a planar screen (a "crack" in the context of linear elasticity), assumed to be a non-empty bounded relatively open subset $Γ$ of the hyperplane $\mathbb{R}^{n-1}\times \{0\}$, on which impedance (Robin) boundary conditions are imposed. In contrast to previous studies, $Γ$ can have arbitrarily rough (possibly fractal) boundary. To obtain well-posedness for such $Γ$ we show how the standard impedance boundary value problem and its associated system of boundary integral equations must be supplemented with additional solution regularity conditions, which hold automatically when $\partialΓ$ is smooth. We show that the associated system of boundary integral operators is compactly perturbed coercive in an appropriate function space setting, strengthening previous results. This permits the use of Mosco convergence to prove convergence of boundary element approximations on smoother "prefractal" screens to the limiting solution on a fractal screen. We present accompanying numerical results, validating our theoretical convergence results, for three-dimensional scattering by a Koch snowflake and a square snowflake. △ Less

Submitted 8 March, 2022; v1 submitted 28 September, 2020; originally announced September 2020.

Realising square and diamond lattice $S=1/2$ Heisenberg antiferromagnet models in the $α$ and $β$ phases of the coordination framework, KTi(C$_2$O$_4$)$_2\cdot$\textit{x}H$_2$O

Authors: Aly H. Abdeldaim, Teng Li, Lewis Farrar, Alexander A. Tsirlin, Wenjiao Yao, Alexandra S. Gibbs, Pascal Manuel, Philip Lightfoot, Gøran J. Nilsen, Lucy Clark

Abstract: We report the crystal structures and magnetic properties of two psuedo-polymorphs of the $S=1/2$ Ti$^{3+}$ coordination framework, KTi(C$_2$O$_4$)$_2\cdot$xH$_2$O. Single-crystal X-ray and powder neutron diffraction measurements on $α$-KTi(C$_2$O$_4$)$_2\cdot$xH$_2$O confirm its structure in the tetragonal $I4/mcm$ space group with a square planar arrangement of Ti$^{3+}$ ions. Magnetometry and sp… ▽ More We report the crystal structures and magnetic properties of two psuedo-polymorphs of the $S=1/2$ Ti$^{3+}$ coordination framework, KTi(C$_2$O$_4$)$_2\cdot$xH$_2$O. Single-crystal X-ray and powder neutron diffraction measurements on $α$-KTi(C$_2$O$_4$)$_2\cdot$xH$_2$O confirm its structure in the tetragonal $I4/mcm$ space group with a square planar arrangement of Ti$^{3+}$ ions. Magnetometry and specific heat measurements reveal weak antiferromagnetic interactions, with $J_1\approx7$ K and $J_2/J_1=0.11$ indicating a slight frustration of nearest- and next-nearest-neighbor interactions. Below $1.8$ K, $α$ undergoes a transition to G-type antiferromagnetic order with magnetic moments aligned along the $c$ axis of the tetragonal structure. The estimated ordered moment of Ti$^{3+}$ in $α$ is suppressed from its spin-only value to $0.62(3)~μ_B$, thus verifying the two-dimensional nature of the magnetic interactions within the system. $β$-KTi(C$_2$O$_4$)$_2\cdot$2H$_2$O, on the other hand, realises a three-dimensional diamond-like magnetic network of Ti$^{3+}$ moments within a hexagonal $P6_222$ structure. An antiferromagnetic exchange coupling of $J\approx54$ K -- an order of magnitude larger than in $α$ -- is extracted from magnetometry and specific heat data. $β$ undergoes Néel ordering at $T_N=28$ K, with the magnetic moments aligned within the $ab$ plane and a slightly reduced ordered moment of $0.79~μ_B$ per Ti$^{3+}$. Through density-functional theory calculations, we address the origin of the large difference in the exchange parameters between the $α$ and $β$ psuedo-polymorphs. Given their observed magnetic behaviors, we propose $α$-KTi(C$_2$O$_4$)$_2\cdot$xH$_2$O and $β$-KTi(C$_2$O$_4$)$_2\cdot$2H$_2$O as close to ideal model $S=1/2$ Heisenberg square and diamond lattice antiferromagnets, respectively. △ Less

Submitted 7 October, 2020; v1 submitted 7 September, 2020; originally announced September 2020.

Comments: 12 pages, 7 figures

EDEN: Sensitivity Analysis and Transiting Planet Detection Limits for Nearby Late Red Dwarfs

Authors: Aidan Gibbs, Alex Bixel, Benjamin Rackham, Daniel Apai, Martin Schlecker, Nestor Espinoza, Luigi Mancini, Wen-Ping Chen, Thomas Henning, Paul Gabor, Richard Boyle, Jose Perez Chavez, Allie Mousseau, Jeremy Dietrich, Quentin Jay Socia, Wing Ip, Chow-Choong Ngeow, Anli Tsai, Asmita Bhandare, Victor Marian, Hans Baehr, Samantha Brown, Maximilian Haberle, Miriam Keppler, Karan Molaverdikhani , et al. (1 additional authors not shown)

Abstract: Small planets are common around late-M dwarfs and can be detected through highly precise photometry by the transit method. Planets orbiting nearby stars are particularly important as they are often the best-suited for future follow-up studies. We present observations of three nearby M-dwarfs referred to as EIC-1, EIC-2, and EIC-3, and use them to search for transits and set limits on the presence… ▽ More Small planets are common around late-M dwarfs and can be detected through highly precise photometry by the transit method. Planets orbiting nearby stars are particularly important as they are often the best-suited for future follow-up studies. We present observations of three nearby M-dwarfs referred to as EIC-1, EIC-2, and EIC-3, and use them to search for transits and set limits on the presence of planets. On most nights our observations are sensitive to Earth-sized transiting planets, and photometric precision is similar to or better than TESS for faint late-M dwarfs of the same magnitude (I=15 mag). We present our photometry and transit search pipeline, which utilizes simple median detrending in combination with transit least squares based transit detection (Hippke & Heller 2019).For these targets, and transiting planets between one and two Earth radii, we achieve an average transit detection probability of 60% between periods of 0.5 and 2 days, 30% between 2 and 5 days,and 10% between 5 and 10 days. These sensitivities are conservative compared to visual searches. △ Less

Submitted 23 February, 2020; originally announced February 2020.

Comments: Accepted to AJ

Fast hybrid numerical-asymptotic boundary element methods for high frequency screen and aperture problems based on least-squares collocation

Authors: Andrew Gibbs, David Hewett, Daan Huybrechs, Emile Parolin

Abstract: We present a hybrid numerical-asymptotic (HNA) boundary element method (BEM) for high frequency scattering by two-dimensional screens and apertures, whose computational cost to achieve any prescribed accuracy remains bounded with increasing frequency. Our method is a collocation implementation of the high order hp HNA approximation space of Hewett et al. IMA J. Numer. Anal. 35 (2015), pp.1698- 172… ▽ More We present a hybrid numerical-asymptotic (HNA) boundary element method (BEM) for high frequency scattering by two-dimensional screens and apertures, whose computational cost to achieve any prescribed accuracy remains bounded with increasing frequency. Our method is a collocation implementation of the high order hp HNA approximation space of Hewett et al. IMA J. Numer. Anal. 35 (2015), pp.1698- 1728, where a Galerkin implementation was studied. An advantage of the current collocation scheme is that the one-dimensional highly oscillatory singular integrals appearing in the BEM matrix entries are significantly easier to evaluate than the two-dimensional integrals appearing in the Galerkin case, which leads to much faster computation times. Here we compute the required integrals at frequency-independent cost using the numerical method of steepest descent, which involves complex contour deformation. The change from Galerkin to collocation is nontrivial because naive collocation implementations based on square linear systems suffer from severe numerical instabilities associated with the numerical redundancy of the HNA basis, which produces highly ill-conditioned BEM matrices. In this paper we show how these instabilities can be removed by oversampling, and solving the resulting overdetermined collocation system in a weighted least-squares sense using a truncated singular value decomposition. On the basis of our numerical experiments, the amount of oversampling required to stabilise the method is modest (around 25% typically suffices) and independent of frequency. As an application of our method we present numerical results for high frequency scattering by prefractal approximations to the middle-third Cantor set. △ Less

Submitted 2 April, 2020; v1 submitted 20 December, 2019; originally announced December 2019.

MSC Class: 65N38; 65R20; 78A45; 78M15; 78M35

Large Easy Axis Anisotropy in the One-Dimensional Magnet BaMo(PO$_4$)$_2$

Authors: Aly H. Abdeldaim, Danis I. Badrtdinov, Alexandra S. Gibbs, Pascal Manuel, Helen C. Walker, Manh Duc Le, Chien Hung Wu, Dariusz Wardecki, Sten-Gunnar Eriksson, Yaroslav O. Kvashnin, Alexander A. Tsirlin, Gøran J. Nilsen

Abstract: We present an extensive experimental and theoretical study on the low-temperature magnetic properties of the monoclinic anhydrous alum compound BaMo(PO$_4$)$_2$. The magnetic susceptibility reveals strong antiferromagnetic interactions $θ_{CW} = -167$~K and long-range magnetic order at $T_N=22$~K, in agreement with a recent report. Powder neutron diffraction furthermore shows that the order is col… ▽ More We present an extensive experimental and theoretical study on the low-temperature magnetic properties of the monoclinic anhydrous alum compound BaMo(PO$_4$)$_2$. The magnetic susceptibility reveals strong antiferromagnetic interactions $θ_{CW} = -167$~K and long-range magnetic order at $T_N=22$~K, in agreement with a recent report. Powder neutron diffraction furthermore shows that the order is collinear, with the moments near the $ac$ plane. Neutron spectroscopy reveals a large excitation gap $Δ= 15$~meV in the low-temperature ordered phase, suggesting a much larger easy-axis spin anisotropy than anticipated. However, the large anisotropy justifies the relatively high ordered moment, Néel temperature, and collinear order observed experimentally, and is furthermore reproduced in a first principles calculations using a new computational scheme. We therefore propose BaMo(PO$_4$)$_2$ to host $S=1$ antiferromagnetic chains with large easy-axis anisotropy, which has been theoretically predicted to realize novel excitation continua. △ Less

Submitted 9 December, 2019; originally announced December 2019.

Comments: 11 pages, 7 figures

High sensitivity heat capacity measurements on Sr2RuO4 under uniaxial pressure

Authors: Y. -S. Li, N. Kikugawa, D. A. Sokolov, F. Jerzembeck, A. S. Gibbs, Y. Maeno, C. W. Hicks, M. Nicklas, A. P. Mackenzie

Abstract: A key question regarding the unconventional superconductivity of Sr$_2$RuO$_4$ remains whether the order parameter is single- or two-component. Under a hypothesis of two-component superconductivity, uniaxial pressure is expected to lift their degeneracy, resulting in a split transition. The most direct and fundamental probe of a split transition is heat capacity. Here, we report development of new… ▽ More A key question regarding the unconventional superconductivity of Sr$_2$RuO$_4$ remains whether the order parameter is single- or two-component. Under a hypothesis of two-component superconductivity, uniaxial pressure is expected to lift their degeneracy, resulting in a split transition. The most direct and fundamental probe of a split transition is heat capacity. Here, we report development of new high-frequency methodology for measurement of heat capacity of samples subject to large and highly homogeneous uniaxial pressure. We place an upper limit on the heat capacity signature of any second transition of a few per cent of the primary superconducting transition. The normalized jump in heat capacity, $ΔC/C$, grows smoothly as a function of uniaxial pressure, but we find no qualitative evidence of a pressure-induced order parameter transition. Thanks to the high precision of our measurements, these findings place stringent constraints on theories of the superconductivity of Sr$_2$RuO$_4$. △ Less

Submitted 30 April, 2020; v1 submitted 18 June, 2019; originally announced June 2019.

Comments: 9 pages, 5 figures, supplementary material included as separate PDF file

Searches After Gravitational-waves Using ARizona Observatories (SAGUARO): System Overview and First Results from Advanced LIGO/Virgo's Third Observing Run

Authors: M. J. Lundquist, K. Paterson, W. Fong, D. J. Sand, J. E. Andrews, I. Shivaei, P. N. Daly, S. Valenti, S. Yang, E. Christensen, A. R. Gibbs, F. Shelly, S. Wyatt, O. Kuhn, R. C. Amaro, I. Arcavi, P. Behroozi, N. Butler, L. Chomiuk, A. Corsi, M. R. Drout, E. Egami, X. Fan, R. J. Foley, B. Frye , et al. (31 additional authors not shown)

Abstract: We present Searches After Gravitational-waves Using ARizona Observatories (SAGUARO), a comprehensive effort dedicated to the discovery and characterization of optical counterparts to gravitational wave (GW) events. SAGUARO utilizes ground-based facilities ranging from 1.5m to 10m in diameter, located primarily in the Northern Hemisphere. We provide an overview of SAGUARO's telescopic resources, pi… ▽ More We present Searches After Gravitational-waves Using ARizona Observatories (SAGUARO), a comprehensive effort dedicated to the discovery and characterization of optical counterparts to gravitational wave (GW) events. SAGUARO utilizes ground-based facilities ranging from 1.5m to 10m in diameter, located primarily in the Northern Hemisphere. We provide an overview of SAGUARO's telescopic resources, pipeline for transient detection, and database for candidate visualization. We describe SAGUARO's discovery component, which utilizes the $5$~deg$^2$ field-of-view optical imager on the Mt. Lemmon 1.5m telescope, reaching limits of $\approx 21.3$~AB mag while rapidly tiling large areas. We also describe the follow-up component of SAGUARO, used for rapid vetting and monitoring of optical candidates. With the onset of Advanced LIGO/Virgo's third observing run, we present results from the first three SAGUARO searches following the GW events S190408an, S190425z and S190426c, which serve as a valuable proof-of-concept of SAGUARO. We triggered and searched 15, 60 and 60 deg$^{2}$ respectively, 17.6, 1.4 and 41.8 hrs after the initial GW alerts. We covered 7.8, 3.0 and 5.1\% of the total probability within the GW event localizations, reaching 3$σ$ limits of 19.8, 21.3 and 20.8 AB mag, respectively. Although no viable counterparts associated with these events were found, we recovered 6 known transients and ruled out 5 potential candidates. We also present Large Binocular Telescope spectroscopy of PS19eq/SN2019ebq, a promising kilonova candidate that was later determined to be a supernova. With the ability to tile large areas and conduct detailed follow-up, SAGUARO represents a significant addition to GW counterpart searches. △ Less

Submitted 17 July, 2019; v1 submitted 14 June, 2019; originally announced June 2019.

Comments: 16 pages, 7 figures, 1 table. Accepted to ApJL

Magnetic frustration and spontaneous rotational symmetry breaking in PdCrO2

Authors: Dan Sun, Dmitry Sokolov, Jack Bartlett, Jhuma Sannigrahi, Seunghyun Khim, Pallavi Kushwaha, Dmitry D. Khalyavin, Pascal Manuel, Alexandra S. Gibbs, Andrew P. Mackenzie, Clifford W. Hicks

Abstract: In the triangular layered magnet PdCrO2 the intralayer magnetic interactions are strong, however the lattice structure frustrates interlayer interactions. In spite of this, long-range, 120$^\circ$ antiferromagnetic order condenses at $T_N = 38$~K. We show here through neutron scattering measurements under in-plane uniaxial stress and in-plane magnetic field that this occurs through a spontaneous l… ▽ More In the triangular layered magnet PdCrO2 the intralayer magnetic interactions are strong, however the lattice structure frustrates interlayer interactions. In spite of this, long-range, 120$^\circ$ antiferromagnetic order condenses at $T_N = 38$~K. We show here through neutron scattering measurements under in-plane uniaxial stress and in-plane magnetic field that this occurs through a spontaneous lifting of the three-fold rotational symmetry of the nonmagnetic lattice, which relieves the interlayer frustration. We also show through resistivity measurements that uniaxial stress can suppress thermal magnetic disorder within the antiferromagnetic phase. △ Less

Submitted 16 May, 2019; v1 submitted 26 April, 2019; originally announced April 2019.

Comments: 9 pages, 9 figures

Journal ref: Phys. Rev. B 100, 094414 (2019)

A high frequency boundary element method for scattering by a class of multiple obstacles

Authors: Andrew Gibbs, Simon Chandler-Wilde, Stephen Langdon, Andrea Moiola

Abstract: We propose a boundary element method for problems of time-harmonic acoustic scattering by multiple obstacles in two dimensions, at least one of which is a convex polygon. By combining a Hybrid Numerical Asymptotic (HNA) approximation space on the convex polygon with standard polynomial-based approximation spaces on each of the other obstacles, we show that the number of degrees of freedom required… ▽ More We propose a boundary element method for problems of time-harmonic acoustic scattering by multiple obstacles in two dimensions, at least one of which is a convex polygon. By combining a Hybrid Numerical Asymptotic (HNA) approximation space on the convex polygon with standard polynomial-based approximation spaces on each of the other obstacles, we show that the number of degrees of freedom required in the HNA space to maintain a given accuracy needs to grow only logarithmically with respect to the frequency, as opposed to the (at least) linear growth required by standard polynomial-based schemes. This method is thus most effective when the convex polygon is many wavelengths in diameter and the small obstacles have a combined perimeter comparable to the problem wavelength. △ Less

Submitted 26 February, 2020; v1 submitted 11 March, 2019; originally announced March 2019.

MSC Class: 35J05; 65N38; 65R20

Intermultiplet transitions and magnetic long-range order in Sm-based pyrochlores

Authors: Viviane Pecanha-Antonio, Erxi Feng, Xiao Sun, Devashibhai Adroja, Helen C. Walker, Alexandra Gibbs, Fabio Orlandi, Yixi Su, Thomas Brückel

Abstract: We present bulk and neutron scattering measurements performed on the isotopically enriched $^{154}\mathrm{Sm_2Ti_2O_7}$ and $^{154}\mathrm{Sm_2Sn_2O_7}$ samples. Both compounds display sharp heat capacity anomalies, at 350 mK and 440 mK, respectively. Inelastic neutron scattering measurements are employed to determine the crystalline electric field (CEF) level scheme, which includes transitions be… ▽ More We present bulk and neutron scattering measurements performed on the isotopically enriched $^{154}\mathrm{Sm_2Ti_2O_7}$ and $^{154}\mathrm{Sm_2Sn_2O_7}$ samples. Both compounds display sharp heat capacity anomalies, at 350 mK and 440 mK, respectively. Inelastic neutron scattering measurements are employed to determine the crystalline electric field (CEF) level scheme, which includes transitions between the ground-state and first excited $J$ multiplets of the $\mathrm{Sm}^{3+}$ ion. To further validate those results, the single-ion magnetic susceptibility of the compounds is calculated and compared with the experimental DC-susceptibility measured in low applied magnetic fields. It is demonstrated that the inclusion of intermultiplet transitions in the CEF analysis is fundamental to the understanding of the intermediate and, more importantly, low temperature magnetic behaviour of the Sm-based pyrochlores. Finally, the heat capacity anomaly is shown to correspond to the onset of an all-in-all-out long-range order in the stannate sample, while in the titanate a dipolar long-range order can be only indirectly inferred. △ Less

Submitted 3 January, 2019; originally announced January 2019.

Comments: 13 pages, 10 Figures

Journal ref: Phys. Rev. B 99, 134415 (2019)

High-resolution photoemission on Sr2RuO4 reveals correlation-enhanced effective spin-orbit coupling and dominantly local self-energies

Authors: A. Tamai, M. Zingl, E. Rozbicki, E. Cappelli, S. Ricco, A. de la Torre, S. McKeown Walker, F. Y. Bruno, P. D. C. King, W. Meevasana, M. Shi, M. Radovic, N. C. Plumb, A. S. Gibbs, A. P. Mackenzie, C. Berthod, H. Strand, M. Kim, A. Georges, F. Baumberger

Abstract: We explore the interplay of electron-electron correlations and spin-orbit coupling in the model Fermi liquid Sr2RuO4 using laser-based angle-resolved photoemission spectroscopy. Our precise measurement of the Fermi surface confirms the importance of spin-orbit coupling in this material and reveals that its effective value is enhanced by a factor of about two, due to electronic correlations. The se… ▽ More We explore the interplay of electron-electron correlations and spin-orbit coupling in the model Fermi liquid Sr2RuO4 using laser-based angle-resolved photoemission spectroscopy. Our precise measurement of the Fermi surface confirms the importance of spin-orbit coupling in this material and reveals that its effective value is enhanced by a factor of about two, due to electronic correlations. The self-energies for the $β$ and $γ$ sheets are found to display significant angular dependence. By taking into account the multi-orbital composition of quasiparticle states, we determine self-energies associated with each orbital component directly from the experimental data. This analysis demonstrates that the perceived angular dependence does not imply momentum-dependent many-body effects, but arises from a substantial orbital mixing induced by spin-orbit coupling. A comparison to single-site dynamical mean-field theory further supports the notion of dominantly local orbital self-energies, and provides strong evidence for an electronic origin of the observed non-linear frequency dependence of the self-energies, leading to `kinks' in the quasiparticle dispersion of Sr2RuO4. △ Less

Submitted 7 May, 2019; v1 submitted 16 December, 2018; originally announced December 2018.

Comments: PRX accepted

Journal ref: Phys. Rev. X 9, 021048 (2019)

VLT/SPHERE Multi-Wavelength High-Contrast Imaging of the HD 115600 Debris Disk: New Constraints on the Dust Geometry and the Presence of Young Giant Planets

Authors: Aidan Gibbs, Kevin Wagner, Daniel Apai, Attila Mo'or, Thayne Currie, Mickael Bonnefoy, Maud Langlois, Carey Lisse

Abstract: Young and dynamically active planetary systems can form disks of debris that are easier to image than the planets themselves. The morphology and evolution of these disks can help to infer the properties of the putative planets responsible for generating and shaping the debris structures. We present integral field spectroscopy and dual-band imaging from VLT/SPHERE ($1.0 - 1.7 μm$) of the debris dis… ▽ More Young and dynamically active planetary systems can form disks of debris that are easier to image than the planets themselves. The morphology and evolution of these disks can help to infer the properties of the putative planets responsible for generating and shaping the debris structures. We present integral field spectroscopy and dual-band imaging from VLT/SPHERE ($1.0 - 1.7 μm$) of the debris disk around the young F2V/F3V star HD 115600. We aim to 1) characterize the geometry and composition of the debris ring, 2) search for thermal emission of young giant planets, and 3) in the absence of detected planets, to refine the inferred properties of plausible planets around HD 115600 to prepare future attempts to detect them. Using a different dust scattering model (ZODIPIC) than in the discovery paper (Henyey-Greenstein) to model the disk geometry, we find $a_0 = 46 \pm 2$ au for the disk's central radius and offsets $Δα$, $Δδ$ = $ -1.0 \pm 0.5, 0.5 \pm 0.5$ au. This offset is smaller than previously found, suggesting that unseen planets of lower masses could be sculpting the disk. Spectroscopy of the disk in Y-J bands with SPHERE shows reddish color, which becomes neutral or slightly blue in H band seen with GPI, broadly consistent with a mixed bulk disk composition of processed organics and water ice. While our observed field contains numerous background objects at wide separations, no exoplanet has been directly observed to a mass sensitivity limit of $2-3(5-7) M_{\rm J}$ between a projected separation of 40 and 200 au for hot (cold)-start models. △ Less

Submitted 21 November, 2018; originally announced November 2018.

Comments: 7 figures

Pressure-induced collapse of spin-orbital Mott state in the hyperhoneycomb iridate $β$-Li$_2$IrO$_3$

Authors: T. Takayama, A. Krajewska, A. S. Gibbs, A. N. Yaresko, H. Ishii, H. Yamaoka, K. Ishii, N. Hiraoka, N. P. Funnell, C. L. Bull, H. Takagi

Abstract: Hyperhoneycomb iridate $β$-Li$_2$IrO$_3$ is a three-dimensional analogue of two-dimensional honeycomb iridates, such as $α$-Li$_2$IrO$_3$, which recently appeared as another playground for the physics of Kitaev-type spin liquid. $β$-Li$_2$IrO$_3$ shows a non-collinear spiral ordering of spin-orbital-entangled $J_{\rm eff}$ = 1/2 moments at low temperature, which is known to be suppressed under a p… ▽ More Hyperhoneycomb iridate $β$-Li$_2$IrO$_3$ is a three-dimensional analogue of two-dimensional honeycomb iridates, such as $α$-Li$_2$IrO$_3$, which recently appeared as another playground for the physics of Kitaev-type spin liquid. $β$-Li$_2$IrO$_3$ shows a non-collinear spiral ordering of spin-orbital-entangled $J_{\rm eff}$ = 1/2 moments at low temperature, which is known to be suppressed under a pressure of $\sim$2 GPa. With further increase of pressure, a structural transition is observed at $P_{\rm S}$ $\sim$ 4 GPa at room temperature. Using the neutron powder diffraction technique, the crystal structure in the high-pressure phase of $β$-Li$_2$IrO$_3$ above $P_{\rm S}$ was refined, which indicates the formation of Ir$_2$ dimers on the zig-zag chains, with the Ir-Ir distance even shorter than that of metallic Ir. We argue that the strong dimerization stabilizes the bonding molecular orbital state comprising the two local $d_{zx}$-orbitals on the Ir-O$_2$-Ir bond plane, which conflicts with the equal superposition of $d_{xy}$-, $d_{yz}$- and $d_{zx}$- orbitals in the $J_{\rm eff}$ = 1/2 wave function produced by strong spin-orbit coupling. The results of resonant inelastic x-ray scattering (RIXS) measurements and the electronic structure calculations are fully consistent with the collapse of the $J_{\rm eff}$ = 1/2 state. A subtle competition of various electronic phases is universal in honeycomb-based Kitaev materials. △ Less

Submitted 16 August, 2018; originally announced August 2018.

Journal ref: Phys. Rev. B 99, 125127 (2019)

Magellan Adaptive Optics Imaging of PDS 70: Measuring the Mass Accretion Rate of a Young Giant Planet within a Gapped Disk

Authors: Kevin Wagner, Katherine B. Follette, Laird M. Close, Dániel Apai, Aidan Gibbs, Miriam Keppler, André Müller, Thomas Henning, Markus Kasper, Ya-Lin Wu, Joseph Long, Jared Males, Katie Morzinski, Melissa McClure

Abstract: PDS 70b is a recently discovered and directly imaged exoplanet within the wide ($\gtrsim$40 au) cavity around PDS 70 (Keppler et al. 2018, Müller et al. 2018). Ongoing accretion onto the central star suggests that accretion onto PDS 70b may also be ongoing. We present the first high contrast images at H$α$ (656 nm) and nearby continuum (643 nm) of PDS 70 utilizing the MagAO system. The combination… ▽ More PDS 70b is a recently discovered and directly imaged exoplanet within the wide ($\gtrsim$40 au) cavity around PDS 70 (Keppler et al. 2018, Müller et al. 2018). Ongoing accretion onto the central star suggests that accretion onto PDS 70b may also be ongoing. We present the first high contrast images at H$α$ (656 nm) and nearby continuum (643 nm) of PDS 70 utilizing the MagAO system. The combination of these filters allows for the accretion rate of the young planet to be inferred, as hot infalling hydrogen gas will emit strongly at H$α$ over the optical continuum. We detected a source in H$α$ at the position of PDS 70b on two sequential nights in May 2018, for which we establish a false positive probability of $<$0.1%. We conclude that PDS 70b is a young, actively accreting planet. We utilize the H$α$ line luminosity to derive a mass accretion rate of $\dot M= 10^{-8\pm1}$ M$_{Jup}/yr$, where the large uncertainty is primarily due to the unknown amount of optical extinction from the circumstellar and circumplanetary disks. PDS 70b represents the second case of an accreting planet interior to a disk gap, and is among the early examples of a planet observed during its formation. △ Less

Submitted 15 December, 2020; v1 submitted 27 July, 2018; originally announced July 2018.

Comments: 6 pages, 5 figures, accepted to ApJL

Micron-scale measurements of low anisotropic strain response of local $T_c$ in Sr$_2$RuO$_4$

Authors: Christopher A. Watson, Alexandra S. Gibbs, Andrew P. Mackenzie, Clifford W. Hicks, Kathryn A. Moler

Abstract: Strontium ruthenate (Sr$_2$RuO$_4$) is a multiband superconductor that displays evidence of topological superconductivity, although a model of the order parameter that is consistent with all experiments remains elusive. We integrated a piezoelectric-based strain apparatus with a scanning superconducting quantum interference device (SQUID) microscope to map the diamagnetic response of single-crysta… ▽ More Strontium ruthenate (Sr$_2$RuO$_4$) is a multiband superconductor that displays evidence of topological superconductivity, although a model of the order parameter that is consistent with all experiments remains elusive. We integrated a piezoelectric-based strain apparatus with a scanning superconducting quantum interference device (SQUID) microscope to map the diamagnetic response of single-crystal Sr$_2$RuO$_4$ as a function of temperature, uniaxial pressure, and position with micron-scale spatial resolution. We thereby obtained local measurements of the superconducting transition temperature $T_c$ vs. anisotropic strain $ε$ with sufficient sensitivity for comparison to theoretical models that assume a uniform $p_x\pm ip_y$ order parameter. We found that $T_c$ varies with position and that the locally measured $T_c$ vs. $ε$ curves are quadratic ($T_c\proptoε^2$), as allowed by the C$_4$ symmetry of the crystal lattice. We did not observe the low-strain linear cusp ($T_c\propto \left| ε\right|$) that would be expected for a two-component order parameter such as $p_x\pm ip_y$. These results provide new input for models of the order parameter of Sr$_2$RuO$_4$. △ Less

Submitted 18 July, 2018; originally announced July 2018.

Comments: 18 pages, 6 figures, 2 tables

Journal ref: Phys. Rev. B 98, 094521 (2018)

Timekeeping infrastructure for the Catalina Sky Survey

Authors: Robert L. Seaman, Alex R. Gibbs

Abstract: Time domain science forms an increasing fraction of astronomical programs at many facilities. Synoptic and targeted observing modes of transient, varying, and moving sources rely on precise clocks to provide the underlying time tags. Often precision is mistaken for accuracy, or the precise time signals never reach the instrumentation in the first place. We will discuss issues of deploying a stable… ▽ More Time domain science forms an increasing fraction of astronomical programs at many facilities. Synoptic and targeted observing modes of transient, varying, and moving sources rely on precise clocks to provide the underlying time tags. Often precision is mistaken for accuracy, or the precise time signals never reach the instrumentation in the first place. We will discuss issues of deploying a stable high-precision GNSS clock on a remote mountaintop, and of conveying the resulting time signals to a computer in a way that permits hardware timestamping of the camera shutter (or equivalent) rather than the arbitrary delays encountered with non-real-time data acquisition software. Strengths and limitations of the Network Time Protocol will be reviewed. Timekeeping infrastructure deployed for the Catalina Sky Survey will serve as an example. △ Less

Submitted 3 July, 2018; originally announced July 2018.

Comments: 12 pages, SPIE proceedings, AS107-93

Spin dynamics of coupled spin ladders near quantum criticality in Ba2CuTeO6

Authors: David Macdougal, Alexandra S. Gibbs, Tao Ying, Stefan Wessel, Helen C. Walker, David Voneshen, Frederic Mila, Hidenori Takagi, Radu Coldea

Abstract: We report inelastic neutron scattering measurements of the magnetic excitations in Ba2CuTeO6, proposed by ab initio calculations to magnetically realize weakly coupled antiferromagnetic two-leg spin-1/2 ladders. Isolated ladders are expected to have a singlet ground state protected by a spin gap. Ba2CuTeO6 orders magnetically, but with a small Neel temperature relative to the exchange strength, su… ▽ More We report inelastic neutron scattering measurements of the magnetic excitations in Ba2CuTeO6, proposed by ab initio calculations to magnetically realize weakly coupled antiferromagnetic two-leg spin-1/2 ladders. Isolated ladders are expected to have a singlet ground state protected by a spin gap. Ba2CuTeO6 orders magnetically, but with a small Neel temperature relative to the exchange strength, suggesting that the interladder couplings are relatively small and only just able to stabilize magnetic order, placing Ba2CuTeO6 close in parameter space to the critical point separating the gapped phase and Neel order. Through comparison of the observed spin dynamics with linear spin wave theory and quantum Monte Carlo calculations, we propose values for all relevant intra- and interladder exchange parameters, which place the system on the ordered side of the phase diagram in proximity to the critical point. We also compare high field magnetization data with quantum Monte Carlo predictions for the proposed model of coupled ladders. △ Less

Submitted 7 December, 2018; v1 submitted 11 June, 2018; originally announced June 2018.

Comments: 14 pages, 12 figures

Journal ref: Phys. Rev. B 98, 174410 (2018)

Numerically stable computation of embedding formulae for scattering by polygons

Authors: Andrew Gibbs, Stephen Langdon, Andrea Moiola

Abstract: For problems of time-harmonic scattering by polygonal obstacles, embedding formulae provide a useful means of computing the far-field coefficient induced by any incident plane wave, given the far-field coefficient of a relatively small set of canonical problems. The number of such problems to be solved depends only on the geometry of the scatterer. Whilst the formulae themselves are exact in theor… ▽ More For problems of time-harmonic scattering by polygonal obstacles, embedding formulae provide a useful means of computing the far-field coefficient induced by any incident plane wave, given the far-field coefficient of a relatively small set of canonical problems. The number of such problems to be solved depends only on the geometry of the scatterer. Whilst the formulae themselves are exact in theory, any implementation will inherit numerical error from the method used to solve the canonical problems. This error can lead to numerical instabilities. Here, we present an effective approach to identify and regulate these instabilities. This approach is subsequently extended to the case where the incident wave is a Herglotz wave function, and we suggest how this could potentially remove frequency dependence of a T-matrix method. △ Less

Submitted 23 May, 2018; originally announced May 2018.

MSC Class: 35J05; 74J20; 76Q05; 78A45

Magnetic interactions in PdCrO$_2$ and their effects on its magnetic structure

Authors: Manh Duc Le, Seyyoung Jeon, A. I. Kolesnikov, D. J. Voneshen, A. S. Gibbs, Jun Sung Kim, Jinwon Jeong, Han-Jin Noh, Changhwi Park, Jaejun Yu, T. G. Perring, Je-Geun Park

Abstract: We report a neutron scattering study of the metallic triangular lattice antiferromagnet PdCrO$_2$. Powder neutron diffraction measurements confirm that the crystalline space group symmetry remains $R\bar{3}m$ below $T_N$. This implies that magnetic interactions consistent with the crystal symmetry do not stabilise the non-coplanar magnetic structure which was one of two structures previously propo… ▽ More We report a neutron scattering study of the metallic triangular lattice antiferromagnet PdCrO$_2$. Powder neutron diffraction measurements confirm that the crystalline space group symmetry remains $R\bar{3}m$ below $T_N$. This implies that magnetic interactions consistent with the crystal symmetry do not stabilise the non-coplanar magnetic structure which was one of two structures previously proposed on the basis of single crystal neutron diffraction measurements. Inelastic neutron scattering measurements find two gaps at low energies which can be explained as arising from a dipolar-type exchange interaction. This symmetric anisotropic interaction also stabilises a magnetic structure very similar to the coplanar magnetic structure which was also suggested by the single crystal diffraction study. The higher energy magnon dispersion can be modelled by linear spin wave theory with exchange interactions up to sixth nearest-neighbors, but discrepancies remain which hint at additional effects unexplained by the linear theory. △ Less

Submitted 16 April, 2018; originally announced April 2018.

Comments: 10 pages, 5 figures, submitted to Phys. Rev. B

Journal ref: Phys. Rev. B 98, 024429 (2018)