Showing 1–50 of 279 results for author: Griffin, M

Herschel Gould Belt Survey in Taurus. II: A census of dense cores and filaments in the TMC1 region

Authors: Jason Kirk, Derek Ward-Thompson, James Di Francesco, Philippe André, David Bresnahan, Vera Könyves, Kenneth Marsh, Matt Griffin, Nicola Schneider, A. Men'shchikov, Pedro Palmeirim, Sylvain Bontemps, Doris Arzoumanian, Milena Benedettini, Stefania Pezzuto

Abstract: We present a catalogue of dense cores and filaments in a $3.8^{\circ}\times2.4^{\circ}$ field around the TMC1 region of the Taurus Molecular Cloud. The catalogue was created using photometric data from the Herschel SPIRE and PACS instruments in the 70 $μ$m, 160 $μ$m, 250 $μ$m, 350 $μ$m, and 500 $μ$m continuum bands. Extended structure in the region was reconstructed from a Herschel column density… ▽ More We present a catalogue of dense cores and filaments in a $3.8^{\circ}\times2.4^{\circ}$ field around the TMC1 region of the Taurus Molecular Cloud. The catalogue was created using photometric data from the Herschel SPIRE and PACS instruments in the 70 $μ$m, 160 $μ$m, 250 $μ$m, 350 $μ$m, and 500 $μ$m continuum bands. Extended structure in the region was reconstructed from a Herschel column density map. Power spectra and PDFs of this structure are presented. The PDF splits into log-normal and power-law forms, with the high-density power-law component associated primarily with the central part of TMC1. The total mass in the mapped region is 2000 M$_\odot$, of which 34% is above an extinction of AV $\sim$ 3 mag -- a level that appears as a break in the PDF and as the minimum column density at which dense cores are found. A total of 35 dense filaments were extracted from the column density map. These have a characteristic FWHM width of 0.07 pc, but the TMC1 filament itself has a mean FWHM of $\sim$ 0.13 pc. The thermally supercritical filaments in the region are aligned orthogonal to the prevailing magnetic field direction. Derived properties for the supercritical TMC1 filament support the scenario of it being relatively young. A catalogue of 44 robust and candidate prestellar cores is created and is assessed to be complete down to 0.1 M$_\odot$. The combined prestellar CMF for the TMC1 and L1495 regions is well fit by a single log-normal distribution and is comparable to the standard IMF. △ Less

Submitted 1 July, 2024; originally announced July 2024.

Comments: 21 pages, 23 figures, accepted for publication in MNRAS. For data products see http://gouldbelt-herschel.cea.fr/archives

Discovery of T center-like quantum defects in silicon

Authors: Yihuang Xiong, Jiongzhi Zheng, Shay McBride, Xueyue Zhang, Sinéad M. Griffin, Geoffroy Hautier

Abstract: Quantum technologies would benefit from the development of high performance quantum defects acting as single-photon emitters or spin-photon interface. Finding such a quantum defect in silicon is especially appealing in view of its favorable spin bath and high processability. While some color centers in silicon have been emerging in quantum applications, there is still a need to search and develop… ▽ More Quantum technologies would benefit from the development of high performance quantum defects acting as single-photon emitters or spin-photon interface. Finding such a quantum defect in silicon is especially appealing in view of its favorable spin bath and high processability. While some color centers in silicon have been emerging in quantum applications, there is still a need to search and develop new high performance quantum emitters. Searching a high-throughput computational database of more than 22,000 charged complex defects in silicon, we identify a series of defects formed by a group III element combined with carbon ((A-C)$\rm _{Si}$ with A=B,Al,Ga,In,Tl) and substituting on a silicon site. These defects are analogous structurally, electronically and chemically to the well-known T center in silicon ((C-C-H)$\rm_{Si}$) and their optical properties are mainly driven by an unpaired electron in a carbon $p$ orbital. They all emit in the telecom and some of these color centers show improved properties compared to the T center in terms of computed radiative lifetime or emission efficiency. We also show that the synthesis of hydrogenated T center-like defects followed by a dehydrogenation annealing step could be an efficient way of synthesis. All the T center-like defects show a higher symmetry than the T center making them easier to align with magnetic fields. Our work motivates further studies on the synthesis and control of this new family of quantum defects, and also demonstrates the use of high-throughput computational screening to detect new complex quantum defects. △ Less

Submitted 8 May, 2024; originally announced May 2024.

Nematicity of a Magnetic Helix

Authors: R. Tumbleson, S. A. Morley, E. Hollingworth, A. Singh, T. Bayaraa, N. G. Burdet, A. U. Saleheen, M. R. McCarter, D. Raftrey, R. J. Pandolfi, V. Esposito, G. L. Dakovski, F. -J. Decker, A. H. Reid, T. A. Assefa, P. Fischer, S. M. Griffin, S. D. Kevan, F. Hellman, J. J. Turner, S. Roy

Abstract: A system that possesses translational symmetry but breaks orientational symmetry is known as a nematic phase. While there are many examples of nematic phases in a wide range of contexts, such as in liquid crystals, complex oxides, and superconductors, of particular interest is the magnetic analogue, where the spin, charge, and orbital degrees of freedom of the electron are intertwined. The difficu… ▽ More A system that possesses translational symmetry but breaks orientational symmetry is known as a nematic phase. While there are many examples of nematic phases in a wide range of contexts, such as in liquid crystals, complex oxides, and superconductors, of particular interest is the magnetic analogue, where the spin, charge, and orbital degrees of freedom of the electron are intertwined. The difficulty of spin nematics is the unambiguous realization and characterization of the phase. Here we present an entirely new type of magnetic nematic phase, which replaces the basis of individual spins with magnetic helices. The helical basis allows for the direct measurement of the order parameters with soft X-ray scattering and a thorough characterization of the nematic phase and its thermodynamic transitions. We discover two distinct nematic phases with unique spatio-temporal correlation signatures. Using coherent X-ray methods, we find that near the phase boundary between the two nematic phases, fluctuations coexist on the timescale of both seconds and sub-nanoseconds. Additionally, we have determined that the fluctuations occur simultaneously with a reorientation of the magnetic helices, indicating that there is spontaneous symmetry breaking and new degrees of freedom become available. Our results provide a novel framework for characterizing exotic phases and the phenomena presented can be mapped onto a broad class of physical systems. △ Less

Submitted 19 April, 2024; originally announced April 2024.

Simulations of Classical Three-Body Thermalization in One Dimension

Authors: M. Eltohfa, Xinghan Wang, Colton M. Griffin, F. Robicheaux

Abstract: One-dimensional systems, such as nanowires or electrons moving along strong magnetic field lines, have peculiar thermalization physics. The binary collision of point-like particles, typically the dominant process for reaching thermal equilibrium in higher dimensional systems, cannot thermalize a 1D system. We study how dilute classical 1D gases thermalize through three-body collisions. We consider… ▽ More One-dimensional systems, such as nanowires or electrons moving along strong magnetic field lines, have peculiar thermalization physics. The binary collision of point-like particles, typically the dominant process for reaching thermal equilibrium in higher dimensional systems, cannot thermalize a 1D system. We study how dilute classical 1D gases thermalize through three-body collisions. We consider a system of identical classical point particles with pairwise repulsive inverse power-law potential $V_{ij} \propto 1/|x_i-x_j|^n$ or the pairwise Lennard-Jones potential. Using Monte Carlo methods, we compute a collision kernel and use it in the Boltzmann equation to evolve a perturbed thermal state with temperature $T$ toward equilibrium. We explain the shape of the kernel and its dependence on the system parameters. Additionally, we implement molecular dynamics simulations of a many-body gas and show agreement with the Boltzmann evolution in the low density limit. For the inverse power-law potential, the rate of thermalization is proportional to $ρ^2 T^{\frac{1}{2}-\frac{1}{n}}$ where $ρ$ is the number density. The corresponding proportionality constant decreases with increasing $n$. △ Less

Submitted 9 July, 2024; v1 submitted 29 February, 2024; originally announced March 2024.

Comments: 13 pages, 12 figures

Compressive-Sensing-Enhanced First-Principles Calculation of Photoluminescence Spectra in Color Centers: A Comparison between Theory and Experiment for the G Center in Silicon

Authors: Jiongzhi Zheng, Lukasz Komza, Yihuang Xiong, Natalya Sheremetyeva, Changpeng Lin, Sinéad M. Griffin, Alp Sipahigil, Geoffroy Hautier

Abstract: Photoluminescence (PL) spectra are a versatile tool for exploring the electronic and optical properties of quantum defect systems. In this work, we investigate the PL spectra of the G center in silicon by combining first-principles computations with a machine-learned compressive-sensing technique and experiment. We show that the compressive-sensing technique provides a speed up of approximately 20… ▽ More Photoluminescence (PL) spectra are a versatile tool for exploring the electronic and optical properties of quantum defect systems. In this work, we investigate the PL spectra of the G center in silicon by combining first-principles computations with a machine-learned compressive-sensing technique and experiment. We show that the compressive-sensing technique provides a speed up of approximately 20 times compared with the finite-displacement method with similar numerical accuracy. We compare theory and experiment and show good agreement for the historically proposed configuration B of the G center. In particular, we attribute the experimentally observed E-line of the G center to a local vibration mode mainly involving two substitutional C atoms and one interstitial Si atom. Our theoretical results also well reproduce and explain the experimental E-line energy shifts originating from the carbon isotopic effect. In addition, our results demonstrate that some highly anharmonic modes that are apparent in computed spectra could be absent experimentally because of their short lifetime. Our work not only provides a deeper understanding of the G-center defect but also paves the way to accelerate the calculation of PL spectra for color centers. △ Less

Submitted 16 February, 2024; v1 submitted 12 February, 2024; originally announced February 2024.

Cool Gaseous Exoplanets: surveying the new frontier with Twinkle

Authors: Luke Booth, Subhajit Sarkar, Matt Griffin, Billy Edwards

Abstract: Cool gaseous exoplanets ($1.75\ R_\oplus < R_\text{p} < 3\ R_\text{J}$, $200$ K $<T_\text{eq} < 1000$~K) are an as-yet understudied population, with great potential to expand our understanding of planetary atmospheres and formation mechanisms. In this paper, we outline the basis for a homogeneous survey of cool gaseous planets with Twinkle, a 0.45-m diameter space telescope with simultaneous spect… ▽ More Cool gaseous exoplanets ($1.75\ R_\oplus < R_\text{p} < 3\ R_\text{J}$, $200$ K $<T_\text{eq} < 1000$~K) are an as-yet understudied population, with great potential to expand our understanding of planetary atmospheres and formation mechanisms. In this paper, we outline the basis for a homogeneous survey of cool gaseous planets with Twinkle, a 0.45-m diameter space telescope with simultaneous spectral coverage from 0.5-4.5~$μ$m, set to launch in 2025. We find that Twinkle has the potential to characterise the atmospheres of 36 known cool gaseous exoplanets (11~sub-Neptunian, 11~Neptunian, 14~Jovian) at an SNR $\geq$ 5 during its 3-year primary mission, with the capability of detecting most major molecules predicted by equilibrium chemistry to > $5σ$ significance. We find that an injected mass-metallicity trend is well-recovered, demonstrating Twinkle's ability to elucidate this fundamental relationship into cool regime. We also find Twinkle will be able to detect cloud layers at 3$σ$ or greater in all cool gaseous planets for clouds at $\leq$ 10 Pa pressure level, but will be insensitive to clouds deeper than $10^4$ Pa in all cases. With these results we demonstrate the capability of the Twinkle mission to greatly expand the current knowledge of cool gaseous planets, enabling key insights and constraints to be obtained for this poorly-charted region of exoplanet parameter space. △ Less

Submitted 9 February, 2024; originally announced February 2024.

Comments: 14 pages, 6 figures, 6 tables

Filamentary Network and Magnetic Field Structures Revealed with BISTRO in the High-Mass Star-Forming Region NGC2264 : Global Properties and Local Magnetogravitational Configurations

Authors: Jia-Wei Wang, Patrick M. Koch, Seamus D. Clarke, Gary Fuller, Nicolas Peretto, Ya-Wen Tang, Hsi-Wei Yen, Shih-Ping Lai, Nagayoshi Ohashi, Doris Arzoumanian, Doug Johnstone, Ray Furuya, Shu-ichiro Inutsuka, Chang Won Lee, Derek Ward-Thompson, Valentin J. M. Le Gouellec, Hong-Li Liu, Lapo Fanciullo, Jihye Hwang, Kate Pattle, Frédérick Poidevin, Mehrnoosh Tahani, Takashi Onaka, Mark G. Rawlings, Eun Jung Chung , et al. (132 additional authors not shown)

Abstract: We report 850 $μ$m continuum polarization observations toward the filamentary high-mass star-forming region NGC 2264, taken as part of the B-fields In STar forming Regions Observations (BISTRO) large program on the James Clerk Maxwell Telescope (JCMT). These data reveal a well-structured non-uniform magnetic field in the NGC 2264C and 2264D regions with a prevailing orientation around 30 deg from… ▽ More We report 850 $μ$m continuum polarization observations toward the filamentary high-mass star-forming region NGC 2264, taken as part of the B-fields In STar forming Regions Observations (BISTRO) large program on the James Clerk Maxwell Telescope (JCMT). These data reveal a well-structured non-uniform magnetic field in the NGC 2264C and 2264D regions with a prevailing orientation around 30 deg from north to east. Field strengths estimates and a virial analysis for the major clumps indicate that NGC 2264C is globally dominated by gravity while in 2264D magnetic, gravitational, and kinetic energies are roughly balanced. We present an analysis scheme that utilizes the locally resolved magnetic field structures, together with the locally measured gravitational vector field and the extracted filamentary network. From this, we infer statistical trends showing that this network consists of two main groups of filaments oriented approximately perpendicular to one another. Additionally, gravity shows one dominating converging direction that is roughly perpendicular to one of the filament orientations, which is suggestive of mass accretion along this direction. Beyond these statistical trends, we identify two types of filaments. The type-I filament is perpendicular to the magnetic field with local gravity transitioning from parallel to perpendicular to the magnetic field from the outside to the filament ridge. The type-II filament is parallel to the magnetic field and local gravity. We interpret these two types of filaments as originating from the competition between radial collapsing, driven by filament self-gravity, and the longitudinal collapsing, driven by the region's global gravity. △ Less

Submitted 23 January, 2024; originally announced January 2024.

Comments: Accepted for publication in the Astrophysical Journal. 43 pages, 32 figures, and 4 tables (including Appendix)

A note on odd partition numbers

Authors: Michael Griffin, Ken Ono

Abstract: Ramanujan's celebrated partition congruences modulo $\ell\in \{5, 7, 11\}$ assert that $$ p(\ell n+δ_{\ell})\equiv 0\pmod{\ell}, $$ where $0<δ_{\ell}<\ell$ satisfies $24δ_{\ell}\equiv 1\pmod{\ell}.$ By proving Subbarao's Conjecture, Radu showed that there are no such congruences when it comes to parity. There are infinitely many odd (resp. even) partition numbers in every arithmetic progression. F… ▽ More Ramanujan's celebrated partition congruences modulo $\ell\in \{5, 7, 11\}$ assert that $$ p(\ell n+δ_{\ell})\equiv 0\pmod{\ell}, $$ where $0<δ_{\ell}<\ell$ satisfies $24δ_{\ell}\equiv 1\pmod{\ell}.$ By proving Subbarao's Conjecture, Radu showed that there are no such congruences when it comes to parity. There are infinitely many odd (resp. even) partition numbers in every arithmetic progression. For primes $\ell \geq 5,$ we give a new proof of the conclusion that there are infinitely many $m$ for which $p(\ell m+δ_{\ell})$ is odd. This proof uses a generalization, due to the second author and Ramsey, of a result of Mazur in his classic paper on the Eisenstein ideal. We also refine a classical criterion of Sturm for modular form congruences, which allows us to show that the smallest such $m$ satisfies $m<(\ell^2-1)/24,$ representing a significant improvement to the previous bound. △ Less

Submitted 16 March, 2024; v1 submitted 1 January, 2024; originally announced January 2024.

Comments: Corrects minor typos found by the referees

MSC Class: 11P81; 11P83; 05A17

Edge-sharing quasi-one-dimensional cuprate fragments in optimally substituted Cu/Pb apatite

Authors: Katherine Inzani, John Vinson, Sinéad M. Griffin

Abstract: The flurry of theoretical and experimental studies following the report of room-temperature superconductivity at ambient pressure in Cu-substituted lead apatite Cu$_x$Pb$_{10-x}$(PO$_4$)$_6$O (`LK99') have explored whether and how this system might host strongly correlated physics including superconductivity. While first-principles calculations at low doping ($x\approx1$) have indicated a Cu-… ▽ More The flurry of theoretical and experimental studies following the report of room-temperature superconductivity at ambient pressure in Cu-substituted lead apatite Cu$_x$Pb$_{10-x}$(PO$_4$)$_6$O (`LK99') have explored whether and how this system might host strongly correlated physics including superconductivity. While first-principles calculations at low doping ($x\approx1$) have indicated a Cu-$d^{9}$ configuration coordinated with oxygen giving rise to isolated, correlated bands, its other structural, electronic, and magnetic properties diverge significantly from those of other known cuprate systems. Here we find that higher densities of ordered Cu substitutions can result in the formation of contiguous edge-sharing Cu-O chains, akin to those found in some members of the cuprate superconductor family. Interestingly, while such quasi-one-dimensional edge-sharing chains are typically ferromagnetically coupled along the chain, we find an antiferromagnetic ground-state magnetic order for our cuprate fragments which is in proximity to a ferromagnetic quantum critical point. This is a result of the elongated Cu-Cu distance in Cu-substituted apatite that leads to larger Cu-O-Cu angles supporting antiferromagnetism, which we demonstrate to be controllable by strain. Finally, our electronic structure calculations confirm the low-dimensional nature of the system and show that the bandwidth is driven by the Cu-O plaquette connectivity, resulting in an intermediate correlated regime. △ Less

Submitted 9 January, 2024; v1 submitted 21 December, 2023; originally announced December 2023.

Comments: 6 pages plus Supplemental. Additional citations in revised version

Electron-phonon coupling in copper-substituted lead phosphate apatite

Authors: Alexander C. Tyner, Sinéad M. Griffin, Alexander V. Balatsky

Abstract: Recent reports of room-temperature, ambient pressure superconductivity in copper-substituted lead phosphate apatite, commonly referred to as LK99, have prompted numerous theoretical and experimental studies into its properties. As the electron-phonon interaction is a common mechanism for superconductivity, the electron-phonon coupling strength is an important quantity to compute for LK99. In this… ▽ More Recent reports of room-temperature, ambient pressure superconductivity in copper-substituted lead phosphate apatite, commonly referred to as LK99, have prompted numerous theoretical and experimental studies into its properties. As the electron-phonon interaction is a common mechanism for superconductivity, the electron-phonon coupling strength is an important quantity to compute for LK99. In this work, we compare the electron-phonon coupling strength among the proposed compositions of LK99. The results of our study are in alignment with the conclusion that LK99 is not a likely candidate for room-temperature superconductivity if electron-phonon interaction is to serve as the mechanism. △ Less

Submitted 14 June, 2024; v1 submitted 19 December, 2023; originally announced December 2023.

Comments: 4 pages, 4 figures, Published Version

Journal ref: J Low Temp Phys (2024)

Ab initio amorphous spin Hamiltonian for the description of topological spin textures in FeGe

Authors: Temuujin Bayaraa, Sinéad M. Griffin

Abstract: Topological spin textures in magnetic materials such as skyrmions and hopfions are interesting manifestations of geometric structures in real materials, concurrently having potential applications as information carriers. In the crystalline systems, the formation of these topological spin textures is well understood as a result of the competition between interactions due to symmetry breaking and fr… ▽ More Topological spin textures in magnetic materials such as skyrmions and hopfions are interesting manifestations of geometric structures in real materials, concurrently having potential applications as information carriers. In the crystalline systems, the formation of these topological spin textures is well understood as a result of the competition between interactions due to symmetry breaking and frustration. However, in systems without translation symmetry such as amorphous materials, a fundamental understanding of the driving mechanisms of non-trivial spin structures is lacking owing to the structural and interaction complexity in these systems. In this work, we use a suite of first-principles-based calculations to propose an ab initio spin Hamiltonian that accurately represents the diversity of structural and magnetic properties in the exemplar amorphous FeGe. Monte Carlo simulations of our amorphous Hamiltonian find emergent skyrmions that are driven by frustrated geometric and magnetic exchange, consistent with those observed in experiment. Moreover, we find that the diversity of local structural motifs results in a large range of exchange interactions, far beyond those found in crystalline materials. Finally, we observe the formation of large-scale emergent structures in amorphous materials, far beyond the relevant interaction length-scale in the systems, suggesting a new route to emergent correlated phases beyond the crystalline limit. △ Less

Submitted 13 November, 2023; originally announced November 2023.

Linear magneto-conductivity as a DC probe of time-reversal symmetry breaking

Authors: Veronika Sunko, Chunxiao Liu, Marc Vila, Ilyoun Na, Yuchen Tang, Vladyslav Kozii, Sinéad M. Griffin, Joel E. Moore, Joseph Orenstein

Abstract: Several optical experiments have shown that in magnetic materials the principal axes of response tensors can rotate in a magnetic field. Here we offer a microscopic explanation of this effect, and propose a closely related DC transport phenomenon -- an off-diagonal \emph{symmetric} conductivity linear in a magnetic field, which we refer to as linear magneto-conductivity (LMC). Although LMC has the… ▽ More Several optical experiments have shown that in magnetic materials the principal axes of response tensors can rotate in a magnetic field. Here we offer a microscopic explanation of this effect, and propose a closely related DC transport phenomenon -- an off-diagonal \emph{symmetric} conductivity linear in a magnetic field, which we refer to as linear magneto-conductivity (LMC). Although LMC has the same functional dependence on magnetic field as the Hall effect, its origin is fundamentally different: LMC requires time-reversal symmetry to be broken even before a magnetic field is applied, and is therefore a sensitive probe of magnetism. We demonstrate LMC in three different ways: via a tight-binding toy model, density functional theory calculations on MnPSe$_3$, and a semiclassical calculation. The third approach additionally identifies two distinct mechanisms yielding LMC: momentum-dependent band magnetization and Berry curvature. Finally, we propose an experimental geometry suitable for detecting LMC, and demonstrate its applicability using Landauer-Büttiker simulations. Our results emphasize the importance of measuring the full conductivity tensor in magnetic materials, and introduce LMC as a new transport probe of symmetry. △ Less

Submitted 24 October, 2023; originally announced October 2023.

Comments: 6+8 pages, 4+3 figures

Chiral Majorana hinge modes on a curved surface with magnetic impurities

Authors: Ilyoun Na, James G. McHugh, Sinéad M. Griffin, Luca Chirolli

Abstract: Chiral Majorana one-dimensional modes have been proposed as they key component for topological quantum computing. In this study, we explore their potential realization as hinge modes in higher-order topological superconductors. To create such phases, we engineer a sign-changing, time-reversal symmetry-breaking mass term through an ensemble of magnetic impurities on the surface of a sphere. The mag… ▽ More Chiral Majorana one-dimensional modes have been proposed as they key component for topological quantum computing. In this study, we explore their potential realization as hinge modes in higher-order topological superconductors. To create such phases, we engineer a sign-changing, time-reversal symmetry-breaking mass term through an ensemble of magnetic impurities on the surface of a sphere. The magnetization of this ensemble arises from the competition between the external magnetic field and the Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction among the impurities, mediated by the surface Majorana modes. We determine the magnetic phase diagram and identify the optimal magnetic field to minimize orbital effects and induce a sign changing mass term. This term opens a gap in the surface spectrum, resulting in a gapless one-dimensional chiral Majorana mode along the nodal line of the mass term, thereby implementing a second-order topological superconductor. △ Less

Submitted 13 October, 2023; originally announced October 2023.

Controlling Topology through Targeted Symmetry Manipulation in Magnetic Systems

Authors: Ilyoun Na, Marc Vila, Sinéad M. Griffin

Abstract: The possibility of selecting magnetic space groups by orienting the magnetization direction or tuning magnetic orders offers a vast playground for engineering symmetry protected topological phases in magnetic materials. In this work, we study how selective tuning of symmetry and magnetism can influence and control the resulting topology in a 2D magnetic system, and illustrate such procedure in the… ▽ More The possibility of selecting magnetic space groups by orienting the magnetization direction or tuning magnetic orders offers a vast playground for engineering symmetry protected topological phases in magnetic materials. In this work, we study how selective tuning of symmetry and magnetism can influence and control the resulting topology in a 2D magnetic system, and illustrate such procedure in the ferromagnetic monolayer MnPSe$_3$. Density functional theory calculations reveals a symmetry-protected accidental semimetalic (SM) phase for out-of-plane magnetization which becomes an insulator when the magnetization is tilted in-plane, reaching band gap values close to $100$ meV. We identify an order-two composite antiunitary symmetry and threefold rotational symmetry that induce the band crossing and classify the possible topological phases using symmetry analysis, which we support with tight-binding and $\mathbf{k}\cdot\mathbf{p}$ models. Breaking of inversion symmetry opens a gap in the SM phase, giving rise to a Chern insulator. We demonstrate this explicitly in the isostructural Janus compound Mn$_2$P$_2$S$_3$Se$_3$, which naturally exhibits Rashba spin-orbit coupling that breaks inversion symmetry. Our results map out the phase space of topological properties of ferromagnetic transition metal phosphorus trichalcogenides and demonstrate the potential of the magnetization-dependent metal-to-insulator transition as a spin switch in integrated two-dimensional electronics. △ Less

Submitted 9 October, 2023; originally announced October 2023.

Realistic non-collinear ground states of solids with source-free exchange correlation functional

Authors: Guy C. Moore, Matthew K. Horton, Aaron D. Kaplan, Sinéad M. Griffin, Kristin A. Persson

Abstract: In this work, we extend the source-free (SF) exchange correlation (XC) functional developed by Sangeeta Sharma and co-workers to plane-wave density functional theory (DFT) based on the projector augmented wave (PAW) method. This constraint is implemented by the current authors within the VASP source code, using a fast Poisson solver that capitalizes on the parallel three-dimensional fast Fourier t… ▽ More In this work, we extend the source-free (SF) exchange correlation (XC) functional developed by Sangeeta Sharma and co-workers to plane-wave density functional theory (DFT) based on the projector augmented wave (PAW) method. This constraint is implemented by the current authors within the VASP source code, using a fast Poisson solver that capitalizes on the parallel three-dimensional fast Fourier transforms (FFTs) implemented in VASP. Using this modified XC functional, we explore the improved convergence behavior that results from applying this constraint to the GGA-PBE+$U$+$J$ functional. In the process, we compare the non-collinear magnetic ground state computed by each functional and their SF counterpart for a select number of magnetic materials in order to provide a metric for comparing with experimentally determined magnetic orderings. We observe significantly improved agreement with experimentally measured magnetic ground state structures after applying the source-free constraint. Furthermore, we explore the importance of considering probability current densities in spin polarized systems, even under no applied field. We analyze the XC torque as well, in order to provide theoretical and computational analyses of the net XC magnetic torque induced by the source-free constraint. Along these lines, we highlight the importance of properly considering the real-space integral of the source-free local magnetic XC field. Our analysis on probability currents, net torque, and constant terms draws additional links to the rich body of previous research on spin-current density functional theory (SCDFT), and paves the way for future extensions and corrections to the SF corrected XC functional. △ Less

Submitted 24 May, 2024; v1 submitted 29 September, 2023; originally announced October 2023.

A Simple Electrode Insulation and Channel Fabrication Technique for High-Electric Field Microfluidics

Authors: Gaurav Anand, Samira Safaripour, Jaynie Tercovich, Jenna Capozzi, Mark Griffin, Nathan Schin, Nicholas Mirra, Craig Snoeyink

Abstract: A simple and robust electrode insulation technique that can withstand a voltage as high as $\mathrm{1000~V}$, which is equivalent to an electric field strength of $\sim 1MV/m$ across a $\mathrm{10~μm}$ channel filled with an electrolyte of conductivity $\sim 0.1~S/m$, i.e., higher than sea water's conductivity, is introduced. A multi-dielectric layers approach is adopted to fabricate the blocked e… ▽ More A simple and robust electrode insulation technique that can withstand a voltage as high as $\mathrm{1000~V}$, which is equivalent to an electric field strength of $\sim 1MV/m$ across a $\mathrm{10~μm}$ channel filled with an electrolyte of conductivity $\sim 0.1~S/m$, i.e., higher than sea water's conductivity, is introduced. A multi-dielectric layers approach is adopted to fabricate the blocked electrodes, which helps reduce the number of material defects. Dielectric insulation with an exceptional breakdown electric field strength for an electrolyte confined between electrodes can have a wide range of applications in microfluidics, like high electric field strength-based dielectrophoresis. The voltage-current characteristics are studied for various concentrations of sodium chloride solution to estimate the insulation strength of the proposed materials, and the breakdown strength is calculated at the point where the electrical insulation failed. A detailed adhesion technique is also demonstrated, which will reduce the ambiguity around the fabrication of a sealed channel using SU-8. △ Less

Submitted 7 September, 2023; v1 submitted 1 August, 2023; originally announced August 2023.

Comments: 12 pages

Origin of correlated isolated flat bands in copper-substituted lead phosphate apatite

Authors: Sinéad M. Griffin

Abstract: A recent report of room temperature superconductivity at ambient pressure in Cu-substituted apatite (`LK99') has invigorated interest in the understanding of what materials and mechanisms can allow for high-temperature superconductivity. Here I perform density functional theory calculations on Cu-substituted lead phosphate apatite, identifying correlated isolated flat bands at the Fermi level, a c… ▽ More A recent report of room temperature superconductivity at ambient pressure in Cu-substituted apatite (`LK99') has invigorated interest in the understanding of what materials and mechanisms can allow for high-temperature superconductivity. Here I perform density functional theory calculations on Cu-substituted lead phosphate apatite, identifying correlated isolated flat bands at the Fermi level, a common signature of high transition temperatures in already established families of superconductors. I elucidate the origins of these isolated bands as arising from a structural distortion induced by the Cu ions and a chiral charge density wave from the Pb lone pairs. These results suggest that a minimal two-band model can encompass much of the low-energy physics in this system. Finally, I discuss the implications of my results on possible superconductivity in Cu-doped apatite. △ Less

Submitted 3 August, 2023; v1 submitted 31 July, 2023; originally announced July 2023.

ContriMix: Scalable stain color augmentation for domain generalization without domain labels in digital pathology

Authors: Tan H. Nguyen, Dinkar Juyal, Jin Li, Aaditya Prakash, Shima Nofallah, Chintan Shah, Sai Chowdary Gullapally, Limin Yu, Michael Griffin, Anand Sampat, John Abel, Justin Lee, Amaro Taylor-Weiner

Abstract: Differences in staining and imaging procedures can cause significant color variations in histopathology images, leading to poor generalization when deploying deep-learning models trained from a different data source. Various color augmentation methods have been proposed to generate synthetic images during training to make models more robust, eliminating the need for stain normalization during test… ▽ More Differences in staining and imaging procedures can cause significant color variations in histopathology images, leading to poor generalization when deploying deep-learning models trained from a different data source. Various color augmentation methods have been proposed to generate synthetic images during training to make models more robust, eliminating the need for stain normalization during test time. Many color augmentation methods leverage domain labels to generate synthetic images. This approach causes three significant challenges to scaling such a model. Firstly, incorporating data from a new domain into deep-learning models trained on existing domain labels is not straightforward. Secondly, dependency on domain labels prevents the use of pathology images without domain labels to improve model performance. Finally, implementation of these methods becomes complicated when multiple domain labels (e.g., patient identification, medical center, etc) are associated with a single image. We introduce ContriMix, a novel domain label free stain color augmentation method based on DRIT++, a style-transfer method. Contrimix leverages sample stain color variation within a training minibatch and random mixing to extract content and attribute information from pathology images. This information can be used by a trained ContriMix model to create synthetic images to improve the performance of existing classifiers. ContriMix outperforms competing methods on the Camelyon17-WILDS dataset. Its performance is consistent across different slides in the test set while being robust to the color variation from rare substances in pathology images. We make our code and trained ContriMix models available for research use. The code for ContriMix can be found at https://gitlab.com/huutan86/contrimix △ Less

Submitted 8 March, 2024; v1 submitted 7 June, 2023; originally announced June 2023.

The JCMT BISTRO Survey: Studying the Complex Magnetic Field of L43

Authors: Janik Karoly, Derek Ward-Thompson, Kate Pattle, David Berry, Anthony Whitworth, Jason Kirk, Pierre Bastien, Tao-Chung Ching, Simon Coude, Jihye Hwang, Woojin Kwon, Archana Soam, Jia-Wei Wang, Tetsuo Hasegawa, Shih-Ping Lai, Keping Qiu, Doris Arzoumanian, Tyler L. Bourke, Do-Young Byun, Huei-Ru Vivien Chen, Wen Ping Chen, Mike Chen, Zhiwei Chen, Jungyeon Cho, Minho Choi , et al. (133 additional authors not shown)

Abstract: We present observations of polarized dust emission at 850 $μ$m from the L43 molecular cloud which sits in the Ophiuchus cloud complex. The data were taken using SCUBA-2/POL-2 on the James Clerk Maxwell Telescope as a part of the BISTRO large program. L43 is a dense ($N_{\rm H_2}\sim 10^{22}$-10$^{23}$ cm$^{-2}$) complex molecular cloud with a submillimetre-bright starless core and two protostellar… ▽ More We present observations of polarized dust emission at 850 $μ$m from the L43 molecular cloud which sits in the Ophiuchus cloud complex. The data were taken using SCUBA-2/POL-2 on the James Clerk Maxwell Telescope as a part of the BISTRO large program. L43 is a dense ($N_{\rm H_2}\sim 10^{22}$-10$^{23}$ cm$^{-2}$) complex molecular cloud with a submillimetre-bright starless core and two protostellar sources. There appears to be an evolutionary gradient along the isolated filament that L43 is embedded within, with the most evolved source closest to the Sco OB2 association. One of the protostars drives a CO outflow that has created a cavity to the southeast. We see a magnetic field that appears to be aligned with the cavity walls of the outflow, suggesting interaction with the outflow. We also find a magnetic field strength of up to $\sim$160$\pm$30 $μ$G in the main starless core and up to $\sim$90$\pm$40 $μ$G in the more diffuse, extended region. These field strengths give magnetically super- and sub-critical values respectively and both are found to be roughly trans-Alfvénic. We also present a new method of data reduction for these denser but fainter objects like starless cores. △ Less

Submitted 22 May, 2023; v1 submitted 18 May, 2023; originally announced May 2023.

Comments: Accepted for publication in ApJ. 23 pages, 9 figures (7 main text, 2 appendix)

Intrinsic Origin and Enhancement of Topological Responses in Ferrimagnetic Antiperovskite Mn4N

Authors: Temuujin Bayaraa, Vsevolod Ivanov, Liang Z. Tan, Sinéad M. Griffin

Abstract: Using first-principles calculations we investigate the intrinsic origins of the anomalous Hall effect (AHE) and the anomalous Nernst effect (ANE) in antiperovskite ferrimagnet Mn4N. We predict that the AHE is significantly enhanced under both compressive and tensile strain, however, the ANE generally decreases under epitaxial strain, except for 1% compressive strain. We connect this behavior to th… ▽ More Using first-principles calculations we investigate the intrinsic origins of the anomalous Hall effect (AHE) and the anomalous Nernst effect (ANE) in antiperovskite ferrimagnet Mn4N. We predict that the AHE is significantly enhanced under both compressive and tensile strain, however, the ANE generally decreases under epitaxial strain, except for 1% compressive strain. We connect this behavior to the evolution of the Berry curvature with strain, suggesting similar strategies for achieving large AHE and ANE changes with modest amounts of strain. Finally, we find that the non-monotonic characteristics of the AHE and ANE stem from the formation and movement of new Weyl points at the periphery of the Brillouin Zone under compressive and tensile strains. △ Less

Submitted 5 April, 2023; originally announced April 2023.

On the effects of strain, defects, and interactions on the topological properties of HfTe5

Authors: Na Hyun Jo, Omar A. Ashour, Zhixue Shu, Chris Jozwiak, Aaron Bostwick, Sae Hee Ryu, Kai Sun, Tai Kong, Sinead M. Griffin, Eli Rotenberg

Abstract: Topological insulators are characterized by spin-momentum-locked massless surface states which are robust under various perturbations. Manipulating such surface states is a topic of vigorous research, as a possible route for the realization of emergent many-body physics in topological systems. Thus far, time-reversal symmetry breaking via Coulomb and magnetic perturbations has been a dominant appr… ▽ More Topological insulators are characterized by spin-momentum-locked massless surface states which are robust under various perturbations. Manipulating such surface states is a topic of vigorous research, as a possible route for the realization of emergent many-body physics in topological systems. Thus far, time-reversal symmetry breaking via Coulomb and magnetic perturbations has been a dominant approach for the tuning of topological states. However, the effect of the structural degrees of freedom on quasi-particle dynamics in topological materials remains elusive. In this work, we demonstrate a transition in HfTe5 between distinct topological phases as a function of either Te vacancy concentration or applied strain; these phases are characterized theoretically as a transition from strong to weak topological insulator and experimentally by a transition from sharp surface states and Dirac crossing to a Fermi-liquid-like quasiparticle state in which these surface-localized features are heavily suppressed. Although vacancies can result in various consequences such as scattering, doping, and structural distortions, we show that changes in the lattice constants play the foremost role in determining the electronic structure, self-energy, and topological states of HfTe5. Our results demonstrate the possibility of using both defect chemistry and strain as control parameters for topological phase transitions and associated many-body physics. △ Less

Submitted 19 March, 2023; originally announced March 2023.

Journal ref: Phys. Rev. B 109, 235122 (2024)

Local atomic stacking and symmetry in twisted graphene trilayers

Authors: Isaac M. Craig, Madeline Van Winkle, Catherine Groschner, Kaidi Zhang, Nikita Dowlatshahi, Ziyan Zhu, Takashi Taniguchi, Kenji Watanabe, Sinéad M. Griffin, D. Kwabena Bediako

Abstract: Moiré superlattices formed from twisting trilayers of graphene are an ideal model for studying electronic correlation, and offer several advantages over bilayer analogues, including more robust and tunable superconductivity and a wide range of twist angles associated with flat band formation. Atomic reconstruction, which strongly impacts the electronic structure of twisted graphene structures, has… ▽ More Moiré superlattices formed from twisting trilayers of graphene are an ideal model for studying electronic correlation, and offer several advantages over bilayer analogues, including more robust and tunable superconductivity and a wide range of twist angles associated with flat band formation. Atomic reconstruction, which strongly impacts the electronic structure of twisted graphene structures, has been suggested to play a major role in the relative versatility of superconductivity in trilayers. Here, we exploit an inteferometric 4D-STEM approach to image a wide range of trilayer graphene structures. Our results unveil a considerably different model for moiré lattice relaxation in trilayers than that proposed from previous measurements, informing a thorough understanding of how reconstruction modulates the atomic stacking symmetries crucial for establishing superconductivity and other correlated phases in twisted graphene trilayers. △ Less

Submitted 6 June, 2024; v1 submitted 16 March, 2023; originally announced March 2023.

Comments: 22 pages, 5 figures

Encoding multistate charge order and chirality in endotaxial heterostructures

Authors: S. Husremović, B. H. Goodge, M. Erodici, K. Inzani, A. Mier, S. M. Ribet, K. C. Bustillo, T. Taniguchi, K. Watanabe, C. Ophus, S. M. Griffin, D. K. Bediako

Abstract: High-density phase change memory (PCM) storage is proposed for materials with multiple intermediate resistance states, which have been observed in 1$T$-TaS$_2$ due to charge density wave (CDW) phase transitions. However, the metastability responsible for this behavior makes the presence of multistate switching unpredictable in TaS$_2$ devices. Here, we demonstrate the fabrication of nanothick vert… ▽ More High-density phase change memory (PCM) storage is proposed for materials with multiple intermediate resistance states, which have been observed in 1$T$-TaS$_2$ due to charge density wave (CDW) phase transitions. However, the metastability responsible for this behavior makes the presence of multistate switching unpredictable in TaS$_2$ devices. Here, we demonstrate the fabrication of nanothick verti-lateral $H$-TaS$_2$/1$T$-TaS$_2$ heterostructures in which the number of endotaxial metallic $H$-TaS$_2$ monolayers dictates the number of resistance transitions in 1$T$-TaS$_2$ lamellae near room temperature. Further, we also observe optically active heterochirality in the CDW superlattice structure, which is modulated in concert with the resistivity steps, and we show how strain engineering can be used to nucleate these polytype conversions. This work positions the principle of endotaxial heterostructures as a promising conceptual framework for reliable, non-volatile, and multi-level switching of structure, chirality, and resistance. △ Less

Submitted 3 September, 2023; v1 submitted 8 March, 2023; originally announced March 2023.

High-throughput identification of spin-photon interfaces in silicon

Authors: Yihuang Xiong, Céline Bourgois, Natalya Sheremetyeva, Wei Chen, Diana Dahliah, Hanbin Song, Sinéad M. Griffin, Alp Sipahigil, Geoffroy Hautier

Abstract: Color centers in host semiconductors are prime candidates for spin-photon interfaces that would enable numerous quantum applications. The discovery of an optimal spin-photon interface in silicon would move quantum information technologies towards a mature semiconductor technology. However, the space of possible charged defects in a host is very large, making the identification of promising quantum… ▽ More Color centers in host semiconductors are prime candidates for spin-photon interfaces that would enable numerous quantum applications. The discovery of an optimal spin-photon interface in silicon would move quantum information technologies towards a mature semiconductor technology. However, the space of possible charged defects in a host is very large, making the identification of promising quantum defects from experiments only extremely challenging. Here, we use high-throughput first principles computational screening to identify spin-photon interfaces among more than 1000 substitutional and interstitial charged defects in silicon. We evaluate the most promising defects by considering their optical properties, spin multiplicity, and formation energies. The use of a single-shot hybrid functional approach is critical in enabling the screening of a large number of defects with a reasonable accuracy in the calculated optical and electronic properties. We identify three new promising spin-photon interface as potential bright emitters in the telecom band: $\rm Ti_{i}^{+}$, $\rm Fe_{i}^{0}$, and $\rm Ru_{i}^{0}$. These candidates are excited through defect-bound excitons, stressing the importance of considering these type of defects in silicon if operations in the telecom band is targeted. Our work paves the way to further large scale computational screening for quantum defects in silicon and other hosts. △ Less

Submitted 2 March, 2023; originally announced March 2023.

Large ordered moment with strong easy-plane anisotropy and vortex-domain pattern in the kagome ferromagnet Fe$_3$Sn

Authors: Lilian Prodan, Donald M. Evans, Sinéad M. Griffin, Andreas Östlin, Markus Altthaler, Erik Lysne, Irina G. Filippova, Serghei Shova, Liviu Chioncel, Vladimir Tsurkan, István Kézsmárki

Abstract: We report the structural and magnetic properties of high-quality bulk single crystals of the kagome ferromagnet Fe$_3$Sn. The dependence of magnetisation on the magnitude and orientation of the external field reveals strong easy-plane type uniaxial magnetic anisotropy, which shows a monotonous increase from $K_1=-0.99\times 10^6 J/m^3$ at 300\,K to $-1.23\times10^6 J/m^3$ at 2\,K. Our \textit{ab i… ▽ More We report the structural and magnetic properties of high-quality bulk single crystals of the kagome ferromagnet Fe$_3$Sn. The dependence of magnetisation on the magnitude and orientation of the external field reveals strong easy-plane type uniaxial magnetic anisotropy, which shows a monotonous increase from $K_1=-0.99\times 10^6 J/m^3$ at 300\,K to $-1.23\times10^6 J/m^3$ at 2\,K. Our \textit{ab initio} electronic structure calculations yield the value of total magnetic moment of about 6.9 $μ_B$/f.u. and a magnetocrystalline anisotropy energy density of 0.406\,meV/f.u. ($1.16\times10^6 J/m^3$) both being in good agreement with the experimental values. The self-consistent DFT computations for the components of the spin/orbital moments indicate that the small difference between the saturation magnetisations measured along and perpendicular to the kagome layers results from the subtle balance between the Fe and Sn spin/orbital moments on the different sites. In zero field, magnetic force microscopy reveals micrometer-scale magnetic vortices with weakly pinned cores that vanish at $\sim$3\,T applied perpendicular to the kagome plane. Our micromagnetic simulations, using the experimentally determined value of anisotropy, well reproduce the observed vortex-domain structure. The present study, in comparison with the easy-axis ferromagnet Fe$_3$Sn$_2$, shows that varying the stacking of kagome layers provides an efficient control over magnetic anisotropy in this family of Fe-based kagome magnets. △ Less

Submitted 6 February, 2023; originally announced February 2023.

Comments: 10 pages, 5 figures

First BISTRO observations of the dark cloud Taurus L1495A-B10: the role of the magnetic field in the earliest stages of low-mass star formation

Authors: Derek Ward-Thompson, Janik Karoly, Kate Pattle, Anthony Whitworth, Jason Kirk, David Berry, Pierre Bastien, Tao-Chung Ching, Simon Coude, Jihye Hwang, Woojin Kwon, Archana Soam, Jia-Wei Wang, Tetsuo Hasegawa, Shih-Ping Lai, Keping Qiu, Doris Arzoumanian, Tyler L. Bourke, Do-Young Byun, Huei-Ru Vivien Chen, Wen Ping Chen, Mike Chen, Zhiwei Chen, Jungyeon Cho, Minho Choi , et al. (133 additional authors not shown)

Abstract: We present BISTRO Survey 850 μm dust emission polarisation observations of the L1495A-B10 region of the Taurus molecular cloud, taken at the JCMT. We observe a roughly triangular network of dense filaments. We detect 9 of the dense starless cores embedded within these filaments in polarisation, finding that the plane-of-sky orientation of the core-scale magnetic field lies roughly perpendicular to… ▽ More We present BISTRO Survey 850 μm dust emission polarisation observations of the L1495A-B10 region of the Taurus molecular cloud, taken at the JCMT. We observe a roughly triangular network of dense filaments. We detect 9 of the dense starless cores embedded within these filaments in polarisation, finding that the plane-of-sky orientation of the core-scale magnetic field lies roughly perpendicular to the filaments in almost all cases. We also find that the large-scale magnetic field orientation measured by Planck is not correlated with any of the core or filament structures, except in the case of the lowest-density core. We propose a scenario for early prestellar evolution that is both an extension to, and consistent with, previous models, introducing an additional evolutionary transitional stage between field-dominated and matter-dominated evolution, observed here for the first time. In this scenario, the cloud collapses first to a sheet-like structure. Uniquely, we appear to be seeing this sheet almost face-on. The sheet fragments into filaments, which in turn form cores. However, the material must reach a certain critical density before the evolution changes from being field-dominated to being matter-dominated. We measure the sheet surface density and the magnetic field strength at that transition for the first time and show consistency with an analytical prediction that had previously gone untested for over 50 years (Mestel 1965). △ Less

Submitted 23 February, 2023; originally announced February 2023.

Comments: 14 pages, 5 figures. ApJ accepted

Midgap state requirements for optically active quantum defects

Authors: Yihuang Xiong, Milena Mathew, Sinéad M. Griffin, Alp Sipahigil, Geoffroy Hautier

Abstract: Optically active quantum defects play an important role in quantum sensing, computing, and communication. The electronic structure and the single-particle energy levels of these quantum defects in the semiconducting host have been used to understand their opto-electronic properties. Optical excitations that are central for their initialization and readout are linked to transitions between occupied… ▽ More Optically active quantum defects play an important role in quantum sensing, computing, and communication. The electronic structure and the single-particle energy levels of these quantum defects in the semiconducting host have been used to understand their opto-electronic properties. Optical excitations that are central for their initialization and readout are linked to transitions between occupied and unoccupied single-particle states. It is commonly assumed that only quantum defects introducing levels well within the band gap and far from the band edges are of interest for quantum technologies as they mimic an isolated atom embedded in the host. In this perspective, we contradict this common assumption and show that optically active defects with energy levels close to the band edges can display similar properties. We highlight quantum defects that are excited through transitions to or from a band-like level (bound exciton), such as the T center and Se$\rm _{Si}^+$ in silicon. We also present how defects such as the silicon divacancy in diamond can involve transitions between localized levels that are above the conduction band or below the valence band. Loosening the commonly assumed requirement on the electronic structure of quantum defects offers opportunities in quantum defects design and discovery, especially in smaller band gap hosts such as silicon. We discuss the challenges in terms of operating temperature for photoluminescence or radiative lifetime in this regime. We also highlight how these alternative type of defects bring their own needs in terms of theoretical developments and fundamental understanding. This perspective clarifies the electronic structure requirement for quantum defects and will facilitate the identification and design of new color centers for quantum applications especially driven by first principles computations. △ Less

Submitted 21 February, 2023; originally announced February 2023.

A Simple Approach for Local and Global Variable Importance in Nonlinear Regression Models

Authors: Emily T. Winn-Nuñez, Maryclare Griffin, Lorin Crawford

Abstract: The ability to interpret machine learning models has become increasingly important as their usage in data science continues to rise. Most current interpretability methods are optimized to work on either (\textit{i}) a global scale, where the goal is to rank features based on their contributions to overall variation in an observed population, or (\textit{ii}) the local level, which aims to detail o… ▽ More The ability to interpret machine learning models has become increasingly important as their usage in data science continues to rise. Most current interpretability methods are optimized to work on either (\textit{i}) a global scale, where the goal is to rank features based on their contributions to overall variation in an observed population, or (\textit{ii}) the local level, which aims to detail on how important a feature is to a particular individual in the data set. In this work, a new operator is proposed called the "GlObal And Local Score" (GOALS): a simple \textit{post hoc} approach to simultaneously assess local and global feature variable importance in nonlinear models. Motivated by problems in biomedicine, the approach is demonstrated using Gaussian process regression where the task of understanding how genetic markers are associated with disease progression both within individuals and across populations is of high interest. Detailed simulations and real data analyses illustrate the flexible and efficient utility of GOALS over state-of-the-art variable importance strategies. △ Less

Submitted 10 August, 2023; v1 submitted 3 February, 2023; originally announced February 2023.

Broad-Range Directional Detection of Light Dark Matter in Cryogenic Ice

Authors: Nora Taufertshöfer, Maurice Garcia-Sciveres, Sinéad M. Griffin

Abstract: We propose hexagonal ice (H$_2$O) as a new target for light dark matter (DM) direct detection. Ice, a polar material, is suitable for single phonon detection through DM scattering for which we consider light dark photon and light scalar mediator models. We report a rate sensitivity down to a DM mass of $\sim$keV, constituting a broader mass range than other promising candidates. We find better sen… ▽ More We propose hexagonal ice (H$_2$O) as a new target for light dark matter (DM) direct detection. Ice, a polar material, is suitable for single phonon detection through DM scattering for which we consider light dark photon and light scalar mediator models. We report a rate sensitivity down to a DM mass of $\sim$keV, constituting a broader mass range than other promising candidates. We find better sensitivity for near-term experimental thresholds from the presence of high-frequency phonons. These advantages, and ice's availability, make it highly promising for single-phonon detection. △ Less

Submitted 13 April, 2023; v1 submitted 11 January, 2023; originally announced January 2023.

Comments: 4 figures

Structural spillage: an efficient method to identify non-crystalline topological materials

Authors: Daniel Muñoz-Segovia, Paul Corbae, Dániel Varjas, Frances Hellman, Sinéad M. Griffin, Adolfo G. Grushin

Abstract: While topological materials are not restricted to crystals, there is no efficient method to diagnose topology in non-crystalline solids such as amorphous materials. Here we introduce the structural spillage, a new indicator that predicts the unknown topological phase of a non-crystalline solid, which is compatible with first-principles calculations. We illustrate its potential with tight-binding a… ▽ More While topological materials are not restricted to crystals, there is no efficient method to diagnose topology in non-crystalline solids such as amorphous materials. Here we introduce the structural spillage, a new indicator that predicts the unknown topological phase of a non-crystalline solid, which is compatible with first-principles calculations. We illustrate its potential with tight-binding and first-principles calculations of amorphous bismuth, predicting a bilayer to be a new topologically nontrivial material. Our work opens up the efficient prediction of non-crystalline solids via first-principles and high-throughput searches. △ Less

Submitted 6 January, 2023; originally announced January 2023.

Comments: 6+18 pages, 4+8 figures

Journal ref: Physical Review Research 5, L042011 (2023)

JCMT BISTRO Observations: Magnetic Field Morphology of Bubbles Associated with NGC 6334

Authors: Mehrnoosh Tahani, Pierre Bastien, Ray S. Furuya, Kate Pattle, Doug Johnstone, Doris Arzoumanian, Yasuo Doi, Tetsuo Hasegawa, Shu-ichiro Inutsuka, Simon Coudé, Laura Fissel, Michael Chun-Yuan Chen, Frédérick Poidevin, Sarah Sadavoy, Rachel Friesen, Patrick M. Koch, James Di Francesco, Gerald H. Moriarty-Schieven, Zhiwei Chen, Eun Jung Chung, Chakali Eswaraiah, Lapo Fanciullo, Tim Gledhill, Valentin J. M. Le Gouellec, Thiem Hoang , et al. (120 additional authors not shown)

Abstract: We study the HII regions associated with the NGC 6334 molecular cloud observed in the sub-millimeter and taken as part of the B-fields In STar-forming Region Observations (BISTRO) Survey. In particular, we investigate the polarization patterns and magnetic field morphologies associated with these HII regions. Through polarization pattern and pressure calculation analyses, several of these bubbles… ▽ More We study the HII regions associated with the NGC 6334 molecular cloud observed in the sub-millimeter and taken as part of the B-fields In STar-forming Region Observations (BISTRO) Survey. In particular, we investigate the polarization patterns and magnetic field morphologies associated with these HII regions. Through polarization pattern and pressure calculation analyses, several of these bubbles indicate that the gas and magnetic field lines have been pushed away from the bubble, toward an almost tangential (to the bubble) magnetic field morphology. In the densest part of NGC 6334, where the magnetic field morphology is similar to an hourglass, the polarization observations do not exhibit observable impact from HII regions. We detect two nested radial polarization patterns in a bubble to the south of NGC 6334 that correspond to the previously observed bipolar structure in this bubble. Finally, using the results of this study, we present steps (incorporating computer vision; circular Hough Transform) that can be used in future studies to identify bubbles that have physically impacted magnetic field lines. △ Less

Submitted 21 December, 2022; originally announced December 2022.

Comments: Accepted for publication in Astrophysical Journal (ApJ)

Rotational and Dilational Reconstruction in Transition Metal Dichalcogenide Moiré Bilayers

Authors: Madeline Van Winkle, Isaac M. Craig, Stephen Carr, Medha Dandu, Karen C. Bustillo, Jim Ciston, Colin Ophus, Takashi Taniguchi, Kenji Watanabe, Archana Raja, Sinéad M. Griffin, D. Kwabena Bediako

Abstract: Lattice reconstruction and corresponding strain accumulation play a key role in defining the electronic structure of two-dimensional moiré superlattices, including those of transition metal dichalcogenides (TMDs). Imaging of TMD moirés has so far provided a qualitative understanding of this relaxation process in terms of interlayer stacking energy, while models of the underlying deformation mechan… ▽ More Lattice reconstruction and corresponding strain accumulation play a key role in defining the electronic structure of two-dimensional moiré superlattices, including those of transition metal dichalcogenides (TMDs). Imaging of TMD moirés has so far provided a qualitative understanding of this relaxation process in terms of interlayer stacking energy, while models of the underlying deformation mechanisms have relied on simulations. Here, we use interferometric four-dimensional scanning transmission electron microscopy to quantitatively map the mechanical deformations through which reconstruction occurs in small-angle twisted bilayer MoS2 and WSe2/MoS2 heterobilayers. We provide direct evidence that local rotations govern relaxation for twisted homobilayers, while local dilations are prominent in heterobilayers possessing a sufficiently large lattice mismatch. Encapsulation of the moiré layers in hBN further localizes and enhances these in-plane reconstruction pathways, suppressing out-of-plane corrugation. We also find that extrinsic uniaxial heterostrain, which introduces a lattice constant difference in twisted homobilayers, leads to accumulation and redistribution of reconstruction strain, demonstrating another route to modify the moiré potential. △ Less

Submitted 10 April, 2023; v1 submitted 13 December, 2022; originally announced December 2022.

Comments: 27 pages, 5 figures

Prediction of topological phases in metastable ferromagnetic MPX$_3$ monolayers

Authors: Natalya Sheremetyeva, Ilyoun Na, Anay Saraf, Sinéad M. Griffin, Geoffroy Hautier

Abstract: Density functional theory calculations are carried out to study the electronic and topological properties of $M$P$X_3$ ($M$ = Mn, Fe, Co, Ni, and $X$ = S, Se) monolayers in the ferromagnetic (FM) metastable magnetic state. We find that FM MnPSe$_3$ monolayers host topological semimetal signatures that are gapped out when spin-orbit coupling (SOC) is included. These findings are supported by explic… ▽ More Density functional theory calculations are carried out to study the electronic and topological properties of $M$P$X_3$ ($M$ = Mn, Fe, Co, Ni, and $X$ = S, Se) monolayers in the ferromagnetic (FM) metastable magnetic state. We find that FM MnPSe$_3$ monolayers host topological semimetal signatures that are gapped out when spin-orbit coupling (SOC) is included. These findings are supported by explicit calculations of the Berry curvature and the Chern number. The choice of the Hubbard-$U$ parameter to describe the $d$-electrons is thoroughly discussed, as well as the influence of using a hybrid-functional approach. The presence of band inversions and the associated topological features are found to be formalism-dependent. Nevertheless, routes to achieve the topological phase via the application of external biaxial strain are demonstrated. Within the hybrid-functional picture, topological band structures are recovered under a pressure of 15% (17 GPa). The present work provides a potential avenue for uncovering new topological phases in metastable ferromagnetic phases. △ Less

Submitted 7 December, 2022; v1 submitted 5 December, 2022; originally announced December 2022.

The JCMT BISTRO-2 Survey: Magnetic Fields of the Massive DR21 Filament

Authors: Tao-Chung Ching, Keping Qiu, Di Li, Zhiyuan Ren, Shih-Ping Lai, David Berry, Kate Pattle, Ray Furuya, Derek Ward-Thompson, Doug Johnstone, Patrick M. Koch, Chang Won Lee, Thiem Hoang, Tetsuo Hasegawa, Woojin Kwon, Pierre Bastien, Chakali Eswaraiah, Jia-Wei Wang, Kyoung Hee Kim, Jihye Hwang, Archana Soam, A-Ran Lyo, Junhao Liu, Valentin J. M. Le Gouellec, Doris Arzoumanian , et al. (132 additional authors not shown)

Abstract: We present 850 $μ$m dust polarization observations of the massive DR21 filament from the B-fields In STar-forming Region Observations (BISTRO) survey, using the POL-2 polarimeter and the SCUBA-2 camera on the James Clerk Maxwell Telescope. We detect ordered magnetic fields perpendicular to the parsec-scale ridge of the DR21 main filament. In the sub-filaments, the magnetic fields are mainly parall… ▽ More We present 850 $μ$m dust polarization observations of the massive DR21 filament from the B-fields In STar-forming Region Observations (BISTRO) survey, using the POL-2 polarimeter and the SCUBA-2 camera on the James Clerk Maxwell Telescope. We detect ordered magnetic fields perpendicular to the parsec-scale ridge of the DR21 main filament. In the sub-filaments, the magnetic fields are mainly parallel to the filamentary structures and smoothly connect to the magnetic fields of the main filament. We compare the POL-2 and Planck dust polarization observations to study the magnetic field structures of the DR21 filament on 0.1--10 pc scales. The magnetic fields revealed in the Planck data are well aligned with those of the POL-2 data, indicating a smooth variation of magnetic fields from large to small scales. The plane-of-sky magnetic field strengths derived from angular dispersion functions of dust polarization are 0.6--1.0 mG in the DR21 filament and $\sim$ 0.1 mG in the surrounding ambient gas. The mass-to-flux ratios are found to be magnetically supercritical in the filament and slightly subcritical to nearly critical in the ambient gas. The alignment between column density structures and magnetic fields changes from random alignment in the low-density ambient gas probed by Planck to mostly perpendicular in the high-density main filament probed by JCMT. The magnetic field structures of the DR21 filament are in agreement with MHD simulations of a strongly magnetized medium, suggesting that magnetic fields play an important role in shaping the DR21 main filament and sub-filaments. △ Less

Submitted 4 December, 2022; originally announced December 2022.

Comments: 26 pages, 13 figures, ApJ accepted

The JCMT BISTRO Survey: A Spiral Magnetic Field in a Hub-filament Structure, Monoceros R2

Authors: Jihye Hwang, Jongsoo Kim, Kate Pattle, Chang Won Lee, Patrick M. Koch, Doug Johnstone, Kohji Tomisaka, Anthony Whitworth, Ray S. Furuya, Ji-hyun Kang, A-Ran Lyo, Eun Jung Chung, Doris Arzoumanian, Geumsook Park, Woojin Kwon, Shinyoung Kim, Motohide Tamura, Jungmi Kwon, Archana Soam, Ilseung Han, Thiem Hoang, Kyoung Hee Kim, Takashi Onaka, Eswaraiah Chakali, Derek Ward-Thompson , et al. (135 additional authors not shown)

Abstract: We present and analyze observations of polarized dust emission at 850 $μ$m towards the central 1 pc $\times$ 1 pc hub-filament structure of Monoceros R2 (Mon R2). The data are obtained with SCUBA-2/POL-2 on the James Clerk Maxwell Telescope (JCMT) as part of the BISTRO (B-fields in Star-forming Region Observations) survey. The orientations of the magnetic field follow the spiral structure of Mon R… ▽ More We present and analyze observations of polarized dust emission at 850 $μ$m towards the central 1 pc $\times$ 1 pc hub-filament structure of Monoceros R2 (Mon R2). The data are obtained with SCUBA-2/POL-2 on the James Clerk Maxwell Telescope (JCMT) as part of the BISTRO (B-fields in Star-forming Region Observations) survey. The orientations of the magnetic field follow the spiral structure of Mon R2, which are well-described by an axisymmetric magnetic field model. We estimate the turbulent component of the magnetic field using the angle difference between our observations and the best-fit model of the underlying large-scale mean magnetic field. This estimate is used to calculate the magnetic field strength using the Davis-Chandrasekhar-Fermi method, for which we also obtain the distribution of volume density and velocity dispersion using a column density map derived from $Herschel$ data and the C$^{18}$O ($J$ = 3-2) data taken with HARP on the JCMT, respectively. We make maps of magnetic field strengths and mass-to-flux ratios, finding that magnetic field strengths vary from 0.02 to 3.64 mG with a mean value of 1.0 $\pm$ 0.06 mG, and the mean critical mass-to-flux ratio is 0.47 $\pm$ 0.02. Additionally, the mean Alfvén Mach number is 0.35 $\pm$ 0.01. This suggests that in Mon R2, magnetic fields provide resistance against large-scale gravitational collapse, and magnetic pressure exceeds turbulent pressure. We also investigate the properties of each filament in Mon R2. Most of the filaments are aligned along the magnetic field direction and are magnetically sub-critical. △ Less

Submitted 13 December, 2022; v1 submitted 12 October, 2022; originally announced October 2022.

Comments: This paper is accepted to the ApJ

Twinkle -- a small satellite spectroscopy mission for the next phase of exoplanet science

Authors: Ian Stotesbury, Billy Edwards, Jean-Francois Lavigne, Vasco Pesquita, James J. Veilleux, Philip Windred, Ahmed Al-Refaie, Lawrence Bradley, Sushuang Ma, Giorgio Savini, Giovanna Tinetti, Til Birnstiel, Sally Dodson-Robinson, Barbara Ercolano, Dax Feliz, Scott Gaudi, Nina Hernitschek, Daniel Holdsworth, Ing-Guey Jiang, Matt Griffin, Nataliea Lowson, Karan Molaverdikhani, Hilding Neilson, Caprice Phillips, Thomas Preibisch , et al. (13 additional authors not shown)

Abstract: With a focus on off-the-shelf components, Twinkle is the first in a series of cost competitive small satellites managed and financed by Blue Skies Space Ltd. The satellite is based on a high-heritage Airbus platform that will carry a 0.45 m telescope and a spectrometer which will provide simultaneous wavelength coverage from 0.5-4.5 $\rm{μm}$. The spacecraft prime is Airbus Stevenage while the tel… ▽ More With a focus on off-the-shelf components, Twinkle is the first in a series of cost competitive small satellites managed and financed by Blue Skies Space Ltd. The satellite is based on a high-heritage Airbus platform that will carry a 0.45 m telescope and a spectrometer which will provide simultaneous wavelength coverage from 0.5-4.5 $\rm{μm}$. The spacecraft prime is Airbus Stevenage while the telescope is being developed by Airbus Toulouse and the spectrometer by ABB Canada. Scheduled to begin scientific operations in 2025, Twinkle will sit in a thermally-stable, sun-synchronous, low-Earth orbit. The mission has a designed operation lifetime of at least seven years and, during the first three years of operation, will conduct two large-scale survey programmes: one focused on Solar System objects and the other dedicated to extrasolar targets. Here we present an overview of the architecture of the mission, refinements in the design approach, and some of the key science themes of the extrasolar survey. △ Less

Submitted 7 September, 2022; originally announced September 2022.

Comments: Presented at SPIE Astronomical Telescopes & Instrumentation 2022

Journal ref: SPIE Proceedings Volume 12180, Space Telescopes and Instrumentation 2022: Optical, Infrared, and Millimeter Wave; 1218033 (2022)

Understanding the chemical shifts of aqueous electrolyte species adsorbed in carbon nanopores

Authors: Anagha Sasikumar, John M. Griffin, Céline Merlet

Abstract: Interfaces between aqueous electrolytes and nanoporous carbons are involved in a number of technological applications such as energy storage and capacitive deionization. The disordered nature of the carbon materials makes it challenging to characterize ion adsorption and relationships between materials properties and performance. Nuclear magnetic spectroscopy can be very helpful in that respect th… ▽ More Interfaces between aqueous electrolytes and nanoporous carbons are involved in a number of technological applications such as energy storage and capacitive deionization. The disordered nature of the carbon materials makes it challenging to characterize ion adsorption and relationships between materials properties and performance. Nuclear magnetic spectroscopy can be very helpful in that respect thanks to its nuclei specificity and ability to distinguish between ions in the bulk and in pores. Nevertheless, several factors can affect the measured chemical shifts making it difficult to interpret experimental results. We use complementary methods, namely density functional theory calculations, molecular dynamics simulations and a mesoscopic model, to investigate various factors affecting the chemical shifts of aqueous electrolyte species. We show that the relative importance of these factors depend on the ion nature. In particular, the hydration shell has a much more pronounced effect for large polarizable ions such as Rb$^+$ and Cs$^+$ compared to Li$^+$. △ Less

Submitted 20 July, 2022; originally announced July 2022.

Nonunique fraction of Fock exchange for defects in two-dimensional materials

Authors: Wei Chen, Sinéad M. Griffin, Gian-Marco Rignanese, Geoffroy Hautier

Abstract: By investigating the vacancy and substitutional defects in monolayer WS$_2$ with hybrid functionals, we find that there is no unique amount of Fock exchange that concurrently satisfies the generalized Koopmans' condition and reproduces the band gap and band-edge positions. Fixing the mixing parameter of Fock exchange based upon the band gap can lead to qualitatively incorrect defect physics in two… ▽ More By investigating the vacancy and substitutional defects in monolayer WS$_2$ with hybrid functionals, we find that there is no unique amount of Fock exchange that concurrently satisfies the generalized Koopmans' condition and reproduces the band gap and band-edge positions. Fixing the mixing parameter of Fock exchange based upon the band gap can lead to qualitatively incorrect defect physics in two-dimensional materials. Instead, excellent agreement is achieved with experiment and many-body perturbation theory within $GW$ approximation once the mixing parameters are tuned individually for the defects and the band edges. We show the departure from a unique optimized mixing parameter is inherent to two-dimensional systems as the band edges experience a reduced screening whilst the localized defects are subject to bulklike screening. △ Less

Submitted 9 October, 2022; v1 submitted 1 June, 2022; originally announced June 2022.

Comments: Main manuscript (7 pages, 3 figures) with supplemental material (3 pages, 5 figures)

Computational needs of quantum mechanical calculations of materials for high-energy physics

Authors: Sinéad M. Griffin

Abstract: Searches for new physics in high-energy physics (HEP) experiments commonly rely on interactions with materials. A burgeoning direction is the accurate calculation and design of materials for HEP applications. In this Snowmass contribution, I briefly motivate the science need for quantum mechanical calculations of materials for HEP and outline the range of questions that such calculations can addre… ▽ More Searches for new physics in high-energy physics (HEP) experiments commonly rely on interactions with materials. A burgeoning direction is the accurate calculation and design of materials for HEP applications. In this Snowmass contribution, I briefly motivate the science need for quantum mechanical calculations of materials for HEP and outline the range of questions that such calculations can address. With this information, I assess the computational needs for ab initio calculations in HEP, the specific computational resources and workflows used by state-of-the-art methods, and finally identify promising future directions such as the use of machine learning and strongly-correlated quantum mechanical calculations moving towards materials calculations on quantum computers. △ Less

Submitted 21 May, 2022; originally announced May 2022.

Comments: Contribution to Snowmass 2021

Performance of the polarization leakage correction in the PILOT data

Authors: J-Ph. Bernard, A. Bernard, H. Roussel, I. Choubani, D. Alina, J. Aumont, A. Hughes, I. Ristorcelli, S. Stever, T. Matsumura S. Sugiyama, K. Komatsu, G. de Gasperis, K. Ferriere, V. Guillet, N. Ysard, P. Ade, P. de Bernardis, N. Bray, B. Crane, J. P. Dubois, M. Griffin, P. Hargrave, Y. Longval, S. Louvel, B. Maffei , et al. (11 additional authors not shown)

Abstract: The Polarized Instrument for Long-wavelength Observation of the Tenuous interstellar medium (PILOT) is a balloon-borne experiment that aims to measure the polarized emission of thermal dust at a wavelength of 240 um (1.2 THz). The PILOT experiment flew from Timmins, Ontario, Canada in 2015 and 2019 and from Alice Springs, Australia in April 2017. The in-flight performance of the instrument during… ▽ More The Polarized Instrument for Long-wavelength Observation of the Tenuous interstellar medium (PILOT) is a balloon-borne experiment that aims to measure the polarized emission of thermal dust at a wavelength of 240 um (1.2 THz). The PILOT experiment flew from Timmins, Ontario, Canada in 2015 and 2019 and from Alice Springs, Australia in April 2017. The in-flight performance of the instrument during the second flight was described in Mangilli et al. 2019. In this paper, we present data processing steps that were not presented in Mangilli et al. 2019 and that we have recently implemented to correct for several remaining instrumental effects. The additional data processing concerns corrections related to detector cross-talk and readout circuit memory effects, and leakage from total intensity to polarization. We illustrate the above effects and the performance of our corrections using data obtained during the third flight of PILOT, but the methods used to assess the impact of these effects on the final science-ready data, and our strategies for correcting them will be applied to all PILOT data. We show that the above corrections, and in particular that for the intensity to polarization leakage, which is most critical for accurate polarization measurements with PILOT, are accurate to better than 0.4 % as measured on Jupiter during flight#3. △ Less

Submitted 7 May, 2022; originally announced May 2022.

New Horizons: Scalar and Vector Ultralight Dark Matter

Authors: D. Antypas, A. Banerjee, C. Bartram, M. Baryakhtar, J. Betz, J. J. Bollinger, C. Boutan, D. Bowring, D. Budker, D. Carney, G. Carosi, S. Chaudhuri, S. Cheong, A. Chou, M. D. Chowdhury, R. T. Co, J. R. Crespo López-Urrutia, M. Demarteau, N. DePorzio, A. V. Derbin, T. Deshpande, M. D. Chowdhury, L. Di Luzio, A. Diaz-Morcillo, J. M. Doyle , et al. (104 additional authors not shown)

Abstract: The last decade has seen unprecedented effort in dark matter model building at all mass scales coupled with the design of numerous new detection strategies. Transformative advances in quantum technologies have led to a plethora of new high-precision quantum sensors and dark matter detection strategies for ultralight ($<10\,$eV) bosonic dark matter that can be described by an oscillating classical,… ▽ More The last decade has seen unprecedented effort in dark matter model building at all mass scales coupled with the design of numerous new detection strategies. Transformative advances in quantum technologies have led to a plethora of new high-precision quantum sensors and dark matter detection strategies for ultralight ($<10\,$eV) bosonic dark matter that can be described by an oscillating classical, largely coherent field. This white paper focuses on searches for wavelike scalar and vector dark matter candidates. △ Less

Submitted 28 March, 2022; originally announced March 2022.

Comments: Snowmass 2021 White Paper

Hard ferromagnetism down to the thinnest limit of iron-intercalated tantalum disulfide

Authors: Samra Husremović, Katherine Inzani, Catherine K. Groschner, Isaac M. Craig, Karen C. Bustillo, Peter Ercius, Nathanael P. Kazmierczak, Jacob Syndikus, Madeline Van Winkle, Shaul Aloni, Takashi Taniguchi, Kenji Watanabe, Sinéad M. Griffin, D. Kwabena Bediako

Abstract: Two-dimensional (2D) magnetic crystals hold promise for miniaturized and ultralow power electronic devices that exploit spin manipulation. In these materials, large, controllable magnetocrystalline anisotropy is a prerequisite for the stabilization and manipulation of long-range magnetic order. In known 2D magnetic crystals, relatively weak magnetocrystalline anisotropy results in typically soft f… ▽ More Two-dimensional (2D) magnetic crystals hold promise for miniaturized and ultralow power electronic devices that exploit spin manipulation. In these materials, large, controllable magnetocrystalline anisotropy is a prerequisite for the stabilization and manipulation of long-range magnetic order. In known 2D magnetic crystals, relatively weak magnetocrystalline anisotropy results in typically soft ferromagnetism. Here, we demonstrate that ferromagnetic order persists down to the thinnest limit of Fe$_x$TaS$_2$ (Fe-intercalated bilayer 2H-TaS$_2$) with giant coercivities up to 3 tesla. We prepare Fe-intercalated TaS$_2$ by chemical intercalation of van der Waals layered 2H-TaS$_2$ crystals and perform variable-temperature quantum transport, transmission electron microscopy, and confocal Raman spectroscopy measurements to shed new light on the coupled effects of dimensionality, degree of intercalation, and intercalant order/disorder on the hard ferromagnetic behavior of Fe$_x$TaS$_2$. More generally, we show that chemical intercalation gives access to a rich synthetic parameter space for low-dimensional magnets, in which magnetic properties can be tailored by the choice of the host material and intercalant identity/amount, in addition to the manifold distinctive degrees of freedom available in atomically thin, van der Waals crystals. △ Less

Submitted 25 July, 2022; v1 submitted 10 March, 2022; originally announced March 2022.

Journal ref: Journal of the American Chemical Society 2022 144, 12167-12176

Improved Pathwise Coordinate Descent for Power Penalties

Authors: Maryclare Griffin

Abstract: Pathwise coordinate descent algorithms have been used to compute entire solution paths for lasso and other penalized regression problems quickly with great success. They improve upon cold start algorithms by solving the problems that make up the solution path sequentially for an ordered set of tuning parameter values, instead of solving each problem separately. However, extending pathwise coordina… ▽ More Pathwise coordinate descent algorithms have been used to compute entire solution paths for lasso and other penalized regression problems quickly with great success. They improve upon cold start algorithms by solving the problems that make up the solution path sequentially for an ordered set of tuning parameter values, instead of solving each problem separately. However, extending pathwise coordinate descent algorithms to more the general bridge or power family of $\ell_q$ penalties is challenging. Faster algorithms for computing solution paths for these penalties are needed because $\ell_q$ penalized regression problems can be nonconvex and especially burdensome to solve. In this paper, we show that a reparameterization of $\ell_q$ penalized regression problems is more amenable to pathwise coordinate descent algorithms. This allows us to improve computation of the mode-thresholding function for $\ell_q$ penalized regression problems in practice and introduce two separate pathwise algorithms. We show that either pathwise algorithm is faster than the corresponding cold-start alternative, and demonstrate that different pathwise algorithms may be more likely to reach better solutions. △ Less

Submitted 14 August, 2023; v1 submitted 4 March, 2022; originally announced March 2022.

First principles study of the T-center in Silicon

Authors: Diana Dhaliah, Yihuang Xiong, Alp Sipahigil, Sinéad M. Griffin, Geoffroy Hautier

Abstract: The T-center in silicon is a well-known carbon-based color center that has been recently considered for quantum technology applications. Using first principles computations, we show that the excited state is formed by a defect-bound exciton made of a localized defect state occupied by an electron to which a hole is bound. The localized state is of strong carbon \textit{p} character and reminiscent… ▽ More The T-center in silicon is a well-known carbon-based color center that has been recently considered for quantum technology applications. Using first principles computations, we show that the excited state is formed by a defect-bound exciton made of a localized defect state occupied by an electron to which a hole is bound. The localized state is of strong carbon \textit{p} character and reminiscent of the localization of the unpaired electron in the ethyl radical molecule. The radiative lifetime for the defect-bound exciton is calculated to be on the order of $μ$s, much longer than other quantum defects such as the NV center in diamond and in agreement with experiments. The longer lifetime is associated with the small transition dipole moment as a result of the very different nature of the localized and delocalized states forming the defect-bound exciton. Finally, we use first principles calculations to assess the stability of the T-center. We find the T-center to be stable against decomposition into simpler defects when keeping the stoichiometry fixed. However, we identify that the T-center is easily prone to (de)hydrogenation and so requires very precise annealing conditions (temperature and atmosphere) to be efficiently formed. △ Less

Submitted 30 March, 2022; v1 submitted 8 February, 2022; originally announced February 2022.

Distributions of Hook lengths in integer partitions

Authors: Michael Griffin, Ken Ono, Wei-Lun Tsai

Abstract: Motivated by the many roles that hook lengths play in mathematics, we study the distribution of the number of $t$-hooks in the partitions of $n$. We prove that the limiting distribution is normal with mean $μ_t(n)\sim \frac{\sqrt{6n}}π-\frac{t}{2}$ and variance $σ_t^2(n)\sim \frac{(π^2-6)\sqrt{6n}}{2π^3}.$ Furthermore, we prove that the distribution of the number of hook lengths that are multiples… ▽ More Motivated by the many roles that hook lengths play in mathematics, we study the distribution of the number of $t$-hooks in the partitions of $n$. We prove that the limiting distribution is normal with mean $μ_t(n)\sim \frac{\sqrt{6n}}π-\frac{t}{2}$ and variance $σ_t^2(n)\sim \frac{(π^2-6)\sqrt{6n}}{2π^3}.$ Furthermore, we prove that the distribution of the number of hook lengths that are multiples of a fixed $t\geq 4$ in partitions of $n$ converge to a shifted Gamma distribution with parameter $k=(t-1)/2$ and scale $θ=\sqrt{2/(t-1)}.$ △ Less

Submitted 22 August, 2022; v1 submitted 17 January, 2022; originally announced January 2022.

Comments: Minor edits and a new first paragraph on previous related work (added at the suggestion of a referee)

MSC Class: 11P82; 05A17

B-fields in Star-Forming Region Observations (BISTRO): Magnetic Fields in the Filamentary Structures of Serpens Main

Authors: Woojin Kwon, Kate Pattle, Sarah Sadavoy, Charles L. H. Hull, Doug Johnstone, Derek Ward-Thompson, James Di Francesco, Patrick M. Koch, Ray Furuya, Yasuo Doi, Valentin J. M. Le Gouellec, Jihye Hwang, A-Ran Lyo, Archana Soam, Xindi Tang, Thiem Hoang, Florian Kirchschlager, Chakali Eswaraiah, Lapo Fanciullo, Kyoung Hee Kim, Takashi Onaka, Vera Könyves, Ji-hyun Kang, Chang Won Lee, Motohide Tamura , et al. (127 additional authors not shown)

Abstract: We present 850 $μ$m polarimetric observations toward the Serpens Main molecular cloud obtained using the POL-2 polarimeter on the James Clerk Maxwell Telescope (JCMT) as part of the B-fields In STar-forming Region Observations (BISTRO) survey. These observations probe the magnetic field morphology of the Serpens Main molecular cloud on about 6000 au scales, which consists of cores and six filament… ▽ More We present 850 $μ$m polarimetric observations toward the Serpens Main molecular cloud obtained using the POL-2 polarimeter on the James Clerk Maxwell Telescope (JCMT) as part of the B-fields In STar-forming Region Observations (BISTRO) survey. These observations probe the magnetic field morphology of the Serpens Main molecular cloud on about 6000 au scales, which consists of cores and six filaments with different physical properties such as density and star formation activity. Using the histogram of relative orientation (HRO) technique, we find that magnetic fields are parallel to filaments in less dense filamentary structures where $N_{H_2} < 0.93\times 10^{22}$ cm$^{-2}$ (magnetic fields perpendicular to density gradients), while being perpendicular to filaments (magnetic fields parallel to density gradients) in dense filamentary structures with star formation activity. Moreover, applying the HRO technique to denser core regions, we find that magnetic field orientations change to become perpendicular to density gradients again at $N_{H_2} \approx 4.6 \times 10^{22}$ cm$^{-2}$. This can be interpreted as a signature of core formation. At $N_{H_2} \approx 16 \times 10^{22}$ cm$^{-2}$ magnetic fields change back to being parallel to density gradients once again, which can be understood to be due to magnetic fields being dragged in by infalling material. In addition, we estimate the magnetic field strengths of the filaments ($B_{POS} = 60-300~μ$G)) using the Davis-Chandrasekhar-Fermi method and discuss whether the filaments are gravitationally unstable based on magnetic field and turbulence energy densities. △ Less

Submitted 13 January, 2022; originally announced January 2022.

Comments: 18 pages, accepted for publication in ApJ

Limiting Betti distributions of Hilbert schemes on $n$ points

Authors: Michael Griffin, Ken Ono, Larry Rolen, Wei-Lun Tsai

Abstract: Hausel and Rodriguez-Villegas recently observed that work of Göttsche, combined with a classical result of Erdős and Lehner on integer partitions, implies that the limiting Betti distribution for the Hilbert schemes $(\mathbb{C}^2)^{[n]}$ on $n$ points, as $n\rightarrow +\infty,$ is a \textit{Gumbel distribution}. In view of this example, they ask for further such Betti distributions. We answer th… ▽ More Hausel and Rodriguez-Villegas recently observed that work of Göttsche, combined with a classical result of Erdős and Lehner on integer partitions, implies that the limiting Betti distribution for the Hilbert schemes $(\mathbb{C}^2)^{[n]}$ on $n$ points, as $n\rightarrow +\infty,$ is a \textit{Gumbel distribution}. In view of this example, they ask for further such Betti distributions. We answer this question for the quasihomogeneous Hilbert schemes $((\mathbb{C}^2)^{[n]})^{T_{α,β}}$ that are cut out by torus actions. We prove that their limiting distributions are also of Gumbel type. To obtain this result, we combine work of Buryak, Feigin, and Nakajima on these Hilbert schemes with our generalization of the result of Erdős and Lehner, which gives the distribution of the number of parts in partitions that are multiples of a fixed integer $A\geq 2.$ Furthermore, if $p_k(A;n)$ denotes the number of partitions of $n$ with exactly $k$ parts that are multiples of $A$, then we obtain the asymptotic $$ p_k(A,n)\sim \frac{24^{\frac k2-\frac14}(n-Ak)^{\frac k2-\frac34}}{\sqrt2\left(1-\frac1A\right)^{\frac k2-\frac14}k!A^{k+\frac12}(2π)^k}e^{2π\sqrt{\frac1{6}\left(1-\frac1A\right)(n-Ak)}}, $$ a result which is of independent interest. △ Less

Submitted 5 April, 2022; v1 submitted 9 December, 2021; originally announced December 2021.

Comments: Minor revision addressing referee comments

MSC Class: 14C05; 14F99; 11P82

Manipulation of spin orientation via ferroelectric switching in Fe-doped Bi2WO6 from first principles

Authors: Katherine Inzani, Nabaraj Pokhrel, Nima Leclerc, Zachary Clemens, Sriram P. Ramkumar, Sinéad M. Griffin, Elizabeth A. Nowadnick

Abstract: Atomic-scale control of spins by electric fields is highly desirable for future technological applications. Magnetically-doped Aurivillius-phase oxides present one route to achieve this, with magnetic ions substituted into the ferroelectric structure at dilute concentrations, resulting in spin-charge coupling. However, there has been minimal exploration of the ferroelectric switching pathways in t… ▽ More Atomic-scale control of spins by electric fields is highly desirable for future technological applications. Magnetically-doped Aurivillius-phase oxides present one route to achieve this, with magnetic ions substituted into the ferroelectric structure at dilute concentrations, resulting in spin-charge coupling. However, there has been minimal exploration of the ferroelectric switching pathways in this materials class, limiting predictions of the influence of an electric field on the magnetic spins in the structure. Here, we determine the ferroelectric switching pathways of the end member of the Aurivilius phase family, Bi2WO6, using a combination of group theoretic analysis and density functional theory calculations. We find that in the ground state P21ab phase, a two-step switching pathway via C2 and Cm intermediate phases provides the lowest energy barrier. Considering iron substitutions on the W-site in Bi2WO6, we determine the spin easy axis. By tracking the change in spin directionality during ferroelectric switching, we find that a 90 degree switch in the polarization direction leads to a 112 degree reorientation of the spin easy axis. The low symmetry crystal-field environment of Bi2WO6 and magnetoelastic coupling on the magnetic dopant provide a route to spin control via and applied electric field. △ Less

Submitted 7 December, 2021; originally announced December 2021.

Comments: 11 pages including Appendices; Submitted to Physical Review B

Log-Gaussian Cox Process Modeling of Large Spatial Lightning Data using Spectral and Laplace Approximations

Authors: Megan L. Gelsinger, Maryclare Griffin, David S. Matteson, Joseph Guinness

Abstract: Lightning is a destructive and highly visible product of severe storms, yet there is still much to be learned about the conditions under which lightning is most likely to occur. The GOES-16 and GOES-17 satellites, launched in 2016 and 2018 by NOAA and NASA, collect a wealth of data regarding individual lightning strike occurrence and potentially related atmospheric variables. The acute nature and… ▽ More Lightning is a destructive and highly visible product of severe storms, yet there is still much to be learned about the conditions under which lightning is most likely to occur. The GOES-16 and GOES-17 satellites, launched in 2016 and 2018 by NOAA and NASA, collect a wealth of data regarding individual lightning strike occurrence and potentially related atmospheric variables. The acute nature and inherent spatial correlation in lightning data renders standard regression analyses inappropriate. Further, computational considerations are foregrounded by the desire to analyze the immense and rapidly increasing volume of lightning data. We present a new computationally feasible method that combines spectral and Laplace approximations in an EM algorithm, denoted SLEM, to fit the widely popular log-Gaussian Cox process model to large spatial point pattern datasets. In simulations, we find SLEM is competitive with contemporary techniques in terms of speed and accuracy. When applied to two lightning datasets, SLEM provides better out-of-sample prediction scores and quicker runtimes, suggesting its particular usefulness for analyzing lightning data, which tend to have sparse signals. △ Less

Submitted 30 November, 2021; originally announced November 2021.

Microscopic Theory of Magnetic Disorder-Induced Decoherence in Superconducting Nb Films

Authors: Evan Sheridan, Thomas F. Harrelson, Eric Sivonxay, Kristin A. Persson, M. Virginia P. Altoé, Irfan Siddiqi, D. Frank Ogletree, David I. Santiago, Sinéad M. Griffin

Abstract: The performance of superconducting qubits is orders of magnitude below what is expected from theoretical estimates based on the loss tangents of the constituent bulk materials. This has been attributed to the presence of uncontrolled surface oxides formed during fabrication which can introduce defects and impurities that create decoherence channels. Here, we develop an ab initio Shiba theory to in… ▽ More The performance of superconducting qubits is orders of magnitude below what is expected from theoretical estimates based on the loss tangents of the constituent bulk materials. This has been attributed to the presence of uncontrolled surface oxides formed during fabrication which can introduce defects and impurities that create decoherence channels. Here, we develop an ab initio Shiba theory to investigate the microscopic origin of magnetic-induced decoherence in niobium thin film superconductors and the formation of native oxides. Our ab initio calculations encompass the roles of structural disorder, stoichiometry, and strain on the formation of decoherence-inducing local spin moments. With parameters derived from these first-principles calculations we develop an effective quasi-classical model of magnetic-induced losses in the superconductor. We identify d-channel losses (associated with oxygen vacancies) as especially parasitic, resulting in a residual zero temperature surface impedance. This work provides a route to connecting atomic scale properties of superconducting materials and macroscopic decoherence channels affecting quantum systems. △ Less

Submitted 23 November, 2021; originally announced November 2021.

Comments: 10 pages plus supplemental materials