Showing 51–100 of 493 results for author: Andersen, M

Interplay of Hidden Orbital Order and Superconductivity in CeCoIn5

Authors: Weijiong Chen, Clara Neerup Breiø, Freek Massee, M. P. Allan, C. Petrovic, J. C. Séamus Davis, P. J. Hirschfeld, Brian M. Andersen, Andreas Kreisel

Abstract: Visualizing atomic-orbital degrees of freedom is a frontier challenge in scanned microscopy. Some types of orbital order are virtually imperceptible to normal scattering techniques because they do not reduce the overall crystal lattice symmetry. A good example is dxz/dyz (π,π) orbital order in tetragonal lattices. For enhanced detectability, here we consider the quasiparticle scattering interferen… ▽ More Visualizing atomic-orbital degrees of freedom is a frontier challenge in scanned microscopy. Some types of orbital order are virtually imperceptible to normal scattering techniques because they do not reduce the overall crystal lattice symmetry. A good example is dxz/dyz (π,π) orbital order in tetragonal lattices. For enhanced detectability, here we consider the quasiparticle scattering interference (QPI) signature of such (π,π) orbital order in both normal and superconducting phases. The theory reveals that sublattice-specific QPI signatures generated by the orbital order should emerge strongly in the superconducting phase. Sublattice-resolved QPI visualization in superconducting CeCoIn5 then reveals two orthogonal QPI patterns at lattice-substitutional impurity atoms. We analyze the energy dependence of these two orthogonal QPI patterns and find the intensity peaked near E=0, as predicted when such (π) orbital order is inte,πrtwined with d-wave superconductivity. Sublattice-resolved superconductive QPI techniques thus represent a new approach for study of hidden orbital order. △ Less

Submitted 9 May, 2023; v1 submitted 28 October, 2022; originally announced October 2022.

Comments: 19 pages, 5 figures

Report number: NBI CMT 2022

Journal ref: Nature Communications 14, 2984 (2023)

Chern-insulator phases and spontaneous spin and valley order in a moiré lattice model for magic-angle twisted bilayer graphene

Authors: Clara N. Breiø, Brian M. Andersen

Abstract: At a certain "magic" relative twist angle of two graphene sheets it remains a challenge to obtain a detailed description of the proliferation of correlated topological electronic phases and their filling-dependence. We perform a self-consistent real-space Hartree-Fock study of an effective moir{é} lattice model to map out the preferred ordered phases as a function of Coulomb interaction strength a… ▽ More At a certain "magic" relative twist angle of two graphene sheets it remains a challenge to obtain a detailed description of the proliferation of correlated topological electronic phases and their filling-dependence. We perform a self-consistent real-space Hartree-Fock study of an effective moir{é} lattice model to map out the preferred ordered phases as a function of Coulomb interaction strength and moir{é} flat-band filling factor. It is found that a quantum valley Hall phase, previously discovered at charge neutrality, is present at all integer fillings for sufficiently large interactions. However, except from charge neutrality additional spontaneous spin/valley polarization is present in the ground state at nonzero integer fillings, leading to Chern-insulator phases and anomalous quantum Hall effects at odd filling factors, thus constituting an example of interaction-driven nontrivial topology. At weaker interactions, all nonzero integer fillings feature metallic inhomogeneous spin/valley ordered phases which may also break additional point group symmetries of the system. We discuss these findings in the light of previous theoretical studies and recent experimental developments of magic-angle twisted bilayer graphene. △ Less

Submitted 26 October, 2022; originally announced October 2022.

Comments: 11 pages, 6 figures

Report number: CMT NBI 2022

Journal ref: Phys. Rev. B 107, 165114 (2023)

Euclid Near Infrared Spectrometer and Photometer instrument flight model presentation, performance and ground calibration results summary

Authors: T. Maciaszek, A. Ealet, W. Gillard, K. Jahnke, R. Barbier, E. Prieto, W. Bon, A. Bonnefoi, A. Caillat, M. Carle, A. Costille, F. Ducret, C. Fabron, B. Foulon, J. L. Gimenez, E. Grassi, M. Jaquet, D. Le Mignant, L. Martin, T. Pamplona, P. Sanchez, J. C. Clémens, L. Caillat, M. Niclas, A. Secroun , et al. (73 additional authors not shown)

Abstract: The NISP (Near Infrared Spectrometer and Photometer) is one of the two Euclid instruments. It operates in the near-IR spectral region (950-2020nm) as a photometer and spectrometer. The instrument is composed of: a cold (135 K) optomechanical subsystem consisting of a Silicon carbide structure, an optical assembly, a filter wheel mechanism, a grism wheel mechanism, a calibration unit, and a thermal… ▽ More The NISP (Near Infrared Spectrometer and Photometer) is one of the two Euclid instruments. It operates in the near-IR spectral region (950-2020nm) as a photometer and spectrometer. The instrument is composed of: a cold (135 K) optomechanical subsystem consisting of a Silicon carbide structure, an optical assembly, a filter wheel mechanism, a grism wheel mechanism, a calibration unit, and a thermal control system, a detection system based on a mosaic of 16 H2RG with their front-end readout electronic, and a warm electronic system (290 K) composed of a data processing / detector control unit and of an instrument control unit that interfaces with the spacecraft via a 1553 bus for command and control and via Spacewire links for science data. This paper presents: the final architecture of the flight model instrument and subsystems, and the performance and the ground calibration measurement done at NISP level and at Euclid Payload Module level at operational cold temperature. △ Less

Submitted 18 October, 2022; originally announced October 2022.

Comments: 18 pages, to appear in Proceedings of the SPIE

Journal ref: Proceedings of the SPIE, Volume 12180, id. 121801K 18 pp. (2022)

The Disk Population in a Distant Massive Protocluster

Authors: Yu Cheng, Jonathan C. Tan, John J. Tobin, Ruben Fedriani, Morten Andersen, Junfeng Wang

Abstract: The unprecedented angular resolution and sensitivity of ALMA makes it possible to unveil disk populations in distant ($>$2 kpc), embedded young cluster environments. We have conducted an observation towards the central region of the massive protocluster G286.21+0.16 at 1.3 mm. With a spatial resolution of 23 mas and a sensitivity of 15 $\rm μJy~beam^{-1}$, we detect a total of 38 protostellar disk… ▽ More The unprecedented angular resolution and sensitivity of ALMA makes it possible to unveil disk populations in distant ($>$2 kpc), embedded young cluster environments. We have conducted an observation towards the central region of the massive protocluster G286.21+0.16 at 1.3 mm. With a spatial resolution of 23 mas and a sensitivity of 15 $\rm μJy~beam^{-1}$, we detect a total of 38 protostellar disks. These disks have dust masses ranging from about 53 to 1825 $M_\oplus$, assuming a dust temperature of 20 K. This sample is not closely associated with previously identified dense cores, as would be expected for disks around Class 0 protostars. Thus, we expect our sample, being flux limited, to be mainly composed of Class I/flat-spectrum source disks, since these are typically more massive than Class II disks. Furthermore, we find that the distributions of disk masses and radii are statistically indistinguishable with those of the Class I/flat-spectrum objects in the Orion molecular cloud, indicating similar processes are operating in G286.21+0.16 to regulate disk formation and evolution. The cluster center appears to host a massive protostellar system composed of three sources within 1200 au, including a potential binary with 600 au projected separation. Relative to this center, there is no evidence for widespread mass segregation in the disk population. We do find a tentative trend of increasing disk radius versus distance from the cluster center, which may point to the influence of dynamical interactions being stronger in the central regions. △ Less

Submitted 18 October, 2022; originally announced October 2022.

Comments: 23 pages, 7 Figures, accepted by ApJ

Anatomy of Spin and Current Generation from Magnetization Gradients in Topological Insulators and Rashba Metals

Authors: Panagiotis Kotetes, Hano O. M. Sura, Brian M. Andersen

Abstract: We explore the spin density and charge currents arising on the surface of a topological insulator and in a 2D Rashba metal due to magnetization gradients. For topological insulators a single interconversion coefficient controls the generation of both quantities. This coefficient is quantized to a value proportional to the vorticity of the Dirac point which constitutes a hallmark of parity anomaly… ▽ More We explore the spin density and charge currents arising on the surface of a topological insulator and in a 2D Rashba metal due to magnetization gradients. For topological insulators a single interconversion coefficient controls the generation of both quantities. This coefficient is quantized to a value proportional to the vorticity of the Dirac point which constitutes a hallmark of parity anomaly at finite density. As such, it also unveils a robust route to disentangle and detect the protected states of a topological insulator on a given surface. In stark contrast, Rashba metals do not exhibit such anomalies since they contain an even number of helical branches. Nonetheless, also these are governed by quantized responses which, however, are not protected against weak disorder. Furthermore, we find that for Rashba metals the interconversion coefficients demonstrate discontinuities and a nontrivial interplay upon varying the chemical potential, the strength of the spin-orbit coupling, and a pairing gap. Our results have implications for the binding between magnetic skyrmions and superconducting vortices, the emergence of Majorana zero modes, and pave the way for superconducting diode effects mediated by out-of-plane magnetization gradients. △ Less

Submitted 16 October, 2023; v1 submitted 13 October, 2022; originally announced October 2022.

Comments: 23 pages, 6 figures and 4 appendices. Extended and improved discussion. Added new results, figures, one appendix, and a few references. To appear in PRB

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

Variability in electricity consumption by category of consumer: the impact on electricity load profiles

Authors: Philipp Andreas Gunkel, Henrik Klinge Jacobsen, Claire-Marie Bergaentzlé, Fabian Scheller, Frits Møller Andersen

Abstract: Residential electrification of transport and heat is changing consumption and its characteristics significantly. Previous studies have demonstrated the impact of socio-techno-economic determinants on residential consumption. However, they fail to capture the distributional characteristics of such consumer groups, which impact network planning and flexibility assessment. Using actual residential el… ▽ More Residential electrification of transport and heat is changing consumption and its characteristics significantly. Previous studies have demonstrated the impact of socio-techno-economic determinants on residential consumption. However, they fail to capture the distributional characteristics of such consumer groups, which impact network planning and flexibility assessment. Using actual residential electricity consumption profile data for 720,000 households in Denmark, we demonstrate that heat pumps are more likely to influence aggregated peak consumption than electric vehicles. At the same time, other socio-economic factors, such as occupancy, dwelling area and income, show little impact. Comparing the extrapolation of a comprehensive rollout of heat pumps or electric vehicles indicates that the most common consumer category deploying heat pumps has 14% more maximum consumption during peak load hours, 46% more average consumption and twice the higher median compared to households owning an electric vehicle. Electric vehicle show already flexibility with coincidence factors that ranges between 5-15% with a maximum of 17% whereas heat pumps are mostly baseload. The detailed and holistic outcomes of this study support flexibility assessment and grid planning in future studies but also the operation of flexible technologies. △ Less

Submitted 7 October, 2022; originally announced October 2022.

Comments: 37 pages, 18 figures, journal article

Bayesian Inference with Projected Densities

Authors: Jasper Marijn Everink, Yiqiu Dong, Martin Skovgaard Andersen

Abstract: Constraints are a natural choice for prior information in Bayesian inference. In various applications, the parameters of interest lie on the boundary of the constraint set. In this paper, we use a method that implicitly defines a constrained prior such that the posterior assigns positive probability to the boundary of the constraint set. We show that by projecting posterior mass onto the constrain… ▽ More Constraints are a natural choice for prior information in Bayesian inference. In various applications, the parameters of interest lie on the boundary of the constraint set. In this paper, we use a method that implicitly defines a constrained prior such that the posterior assigns positive probability to the boundary of the constraint set. We show that by projecting posterior mass onto the constraint set, we obtain a new posterior with a rich probabilistic structure on the boundary of that set. If the original posterior is a Gaussian, then such a projection can be done efficiently. We apply the method to Bayesian linear inverse problems, in which case samples can be obtained by repeatedly solving constrained least squares problems, similar to a MAP estimate, but with perturbations in the data. When combined into a Bayesian hierarchical model and the constraint set is a polyhedral cone, we can derive a Gibbs sampler to efficiently sample from the hierarchical model. To show the effect of projecting the posterior, we applied the method to deblurring and computed tomography examples. △ Less

Submitted 26 September, 2022; originally announced September 2022.

MSC Class: 62F15; 65C05; 90C25

Heat capacity double transitions in time-reversal symmetry broken superconductors

Authors: Henrik S. Røising, Glenn Wagner, Mercè Roig, Astrid T. Rømer, Brian M. Andersen

Abstract: Standard superconductors display a ubiquitous discontinuous jump in the electronic specific heat at the critical superconducting transition temperature. In a growing class of unconventional superconductors, however, a second order parameter component may get stabilized and produce a second heat capacity jump at a lower temperature, typically associated with the spontaneous breaking of time-reversa… ▽ More Standard superconductors display a ubiquitous discontinuous jump in the electronic specific heat at the critical superconducting transition temperature. In a growing class of unconventional superconductors, however, a second order parameter component may get stabilized and produce a second heat capacity jump at a lower temperature, typically associated with the spontaneous breaking of time-reversal symmetry. The splitting of the two specific heat discontinuities can be controlled by external perturbations such as chemical substitution, hydrostatic pressure, or uniaxial strain. We develop a theoretical quantitative multi-band framework to determine the ratio of the heat capacity jumps, given the band structure and the order parameter momentum structure. We discuss the conditions of the gap profile which determine the amplitude of the second jump. We apply our formalism to the case of Sr$_2$RuO$_4$, and using the gap functions from a microscopic random phase approximation calculation, we show that recently-proposed accidentally degenerate order parameters may exhibit a strongly suppressed second heat capacity jump. We discuss the origin of this result and consider also the role of spatial inhomogeneity on the specific heat. Our results provide a possible explanation of why a second heat capacity jump has so far evaded experimental detection in Sr$_2$RuO$_4$. △ Less

Submitted 30 November, 2022; v1 submitted 8 September, 2022; originally announced September 2022.

Comments: 13+8 pages, 8+8 figures

Report number: NBI CMT 2022

Journal ref: Phys. Rev. B 106, 174518 (2022)

Gemini North Adaptive Optics (GNAO) facility overview and status updates

Authors: Gaetano Sivo, Julia Scharwächter, Manuel Lazo, Célia Blain, Stephen Goodsell, Marcos van Dam, Martin Tschimmel, Henry Roe, Jennifer Lotz, Kim Tomassino-Reed, William Rambold, Courtney Raich, Ricardo Cardenes, Angelic Ebbers, Tim Gaggstatter, Pedro Gigoux, Thomas Schneider, Charles Cavedoni, Stacy Kang, Stanislas Karewicz, Heather Carr, Jesse Ball, Paul Hirst, Emmanuel Chirre, John White , et al. (32 additional authors not shown)

Abstract: The Gemini North Adaptive Optics (GNAO) facility is the upcoming AO facility for Gemini North providing a state-of-the-art AO system for surveys and time domain science in the era of JWST and Rubin operations. GNAO will be optimized to feed the Gemini infrared Multi Object Spectrograph (GIRMOS). While GIRMOS is the primary science driver for defining the capabilities of GNAO, any instrument oper… ▽ More The Gemini North Adaptive Optics (GNAO) facility is the upcoming AO facility for Gemini North providing a state-of-the-art AO system for surveys and time domain science in the era of JWST and Rubin operations. GNAO will be optimized to feed the Gemini infrared Multi Object Spectrograph (GIRMOS). While GIRMOS is the primary science driver for defining the capabilities of GNAO, any instrument operating with an f/32 beam can be deployed using GNAO. The GNAO project includes the development of a new laser guide star facility which will consist of four side-launched laser beams supporting the two primary AO modes of GNAO: a wide-field mode providing an improved image quality over natural seeing for a 2-arcminute circular field-of-view and a narrow-field mode providing near diffraction-limited performance over a 20x20 arcsecond square field-of-view. The GNAO wide field mode will enable GIRMOS's multi-IFU configuration in which the science beam to each individual IFU will be additionally corrected using multi-object AO within GIRMOS. The GNAO narrow field mode will feed the GIRMOS tiled IFU configuration in which all IFUs are combined into a "super"-IFU in the center of the field. GNAO also includes the development of a new Real Time Controller, a new GNAO Facility System Controller and finally the development of a new AO Bench. We present in this paper an overview of the GNAO facility and provide a status update of each product. △ Less

Submitted 30 August, 2022; originally announced August 2022.

Comments: SPIE conference 2022 Montreal

Superconductivity from repulsive interactions on the kagome lattice

Authors: Astrid T. Rømer, Shinibali Bhattacharyya, Roser Valentí, Morten H. Christensen, Brian M. Andersen

Abstract: The discovery of superconductivity in layered vanadium-based kagome metals $A$V$_3$Sb$_5$ ($A$: K, Rb, Cs) has added a new family of materials to the growing class of possible unconventional superconductors. However, the nature of the superconducting pairing in these materials remains elusive. We present a microscopic theoretical study of the leading superconducting instabilities on the kagome lat… ▽ More The discovery of superconductivity in layered vanadium-based kagome metals $A$V$_3$Sb$_5$ ($A$: K, Rb, Cs) has added a new family of materials to the growing class of possible unconventional superconductors. However, the nature of the superconducting pairing in these materials remains elusive. We present a microscopic theoretical study of the leading superconducting instabilities on the kagome lattice based on spin- and charge-fluctuation mediated Cooper pairing. The applied methodology includes effects of both on-site and nearest-neighbor repulsive Coulomb interactions. Near the upper van Hove filling -- relevant for the $A$V$_3$Sb$_5$ materials -- we find a rich phase diagram with several pairing symmetries being nearly degenerate. In particular, while a substantial fraction of the phase diagram is occupied by a spin-singlet order parameter transforming as a two-dimensional irreducible representation of the point group, several nodal spin-triplet pairing states remain competitive. We compute the band and interaction parameter-dependence of the hierarchy of the leading superconducting instabilities, and determine the detailed momentum dependence of the resulting preferred gap structures. Crucially, for moderate values of the interaction parameters, the individual pairing states depend strongly on momentum and exhibit multiple nodes on the Fermi surface. We discuss the properties of these superconducting gap structures in light of recent experimental developments of the $A$V$_3$Sb$_5$ materials. △ Less

Submitted 29 August, 2022; originally announced August 2022.

Comments: 14 pages, 7 figures

Report number: NBI CMT 2022

Journal ref: Phys. Rev. B 106, 174514 (2022)

Resolving the core of R136 in the optical

Authors: Venu M. Kalari, Elliott P. Horch, Ricardo Salinas, Jorick S. Vink, Morten Andersen, Joachim M. Bestenlehner, Monica Rubio

Abstract: The sharpest optical images of the R136 cluster in the Large Magellanic Cloud are presented, allowing for the first time to resolve members of the central core, including R136a1, the most massive star known. These data were taken using the Gemini speckle imager Zorro in medium-band filters with effective wavelengths similar to BVRI achieving angular resolutions between 30-40 mas. All stars previou… ▽ More The sharpest optical images of the R136 cluster in the Large Magellanic Cloud are presented, allowing for the first time to resolve members of the central core, including R136a1, the most massive star known. These data were taken using the Gemini speckle imager Zorro in medium-band filters with effective wavelengths similar to BVRI achieving angular resolutions between 30-40 mas. All stars previously known in the literature, having $V<$16 mag within the central $2''\times2''$ were recovered. Visual companions ($\geq$40 mas; 2000 au) were detected for the WN5h stars R136 a1, and a3. Photometry of the visual companion of a1 suggests it is of mid O spectral type. Based on new photometric luminosities using the resolved Zorro imaging, the masses of the individual WN5h stars are estimated to be between 150-200 $M_{\odot}$, lowering significantly the present-day masses of some of the most massive stars known. These mass estimates are critical anchor points for establishing the stellar upper-mass function. △ Less

Submitted 11 August, 2022; v1 submitted 26 July, 2022; originally announced July 2022.

Comments: Accepted for publication in ApJ; 14 pages, 8 figures, 3 tables

Loop currents in $A$V$_3$Sb$_5$ kagome metals: multipolar and toroidal magnetic orders

Authors: Morten H. Christensen, Turan Birol, Brian M. Andersen, Rafael M. Fernandes

Abstract: Experiments in the recently discovered vanadium-based kagome metals have suggested that their charge-ordered state displays not only bond distortions, characteristic of a ``real" charge density-wave (rCDW), but also time-reversal symmetry-breaking, typical of loop currents described by an ``imaginary" charge density-wave (iCDW). Here, we combine density-functional theory, group-theory, and phenome… ▽ More Experiments in the recently discovered vanadium-based kagome metals have suggested that their charge-ordered state displays not only bond distortions, characteristic of a ``real" charge density-wave (rCDW), but also time-reversal symmetry-breaking, typical of loop currents described by an ``imaginary" charge density-wave (iCDW). Here, we combine density-functional theory, group-theory, and phenomenological modeling to investigate the complex charge-ordered states that arise from interactions between the low-energy van Hove singularities present in the electronic structure of $A$V$_3$Sb$_5$. We find two broad classes of mixed iCDW-rCDW configurations: triple-$\mathbf{Q}$ iCDW, triple-$\mathbf{Q}$ rCDW order, dubbed $3\mathbf{Q}$-$3\mathbf{Q}$, and double-$\mathbf{Q}$ iCDW, single-$\mathbf{Q}$ rCDW order, dubbed $2\mathbf{Q}$-$1\mathbf{Q}$. Moreover, we identify seven different types of iCDW order, stemming from the different vanadium-orbital and kagome-sublattice structures of the two pairs of van Hove singularities present above and below the Fermi level. While the $2\mathbf{Q}$-$1\mathbf{Q}$ states trigger an orthorhombic distortion that breaks the threefold rotational symmetry of the kagome lattice, the $3\mathbf{Q}$-$3\mathbf{Q}$ states induce various types of subsidiary uniform magnetic orders, from conventional ferromagnetism to magnetic octupolar, magnetic toroidal, and even magnetic monopolar order. We show that these exotic orders display unique magneto-striction, magneto-electric, and magneto-electric-striction properties that can be probed experimentally to identify which iCDW state is realized in these compounds. We briefly discuss the impact of an out-of-plane modulation of the charge order and the interplay between these complex charge-ordered states and superconductivity. △ Less

Submitted 14 October, 2022; v1 submitted 26 July, 2022; originally announced July 2022.

Comments: 20 pages, 8 figures

Report number: NBI CMT 2022

Journal ref: Phys. Rev. B 106, 144504 (2022)

Predicting binding energies of astrochemically relevant molecules via machine learning

Authors: Torben Villadsen, Niels F. W. Ligterink, Mie Andersen

Abstract: The behaviour of molecules in space is to a large extent governed by where they freeze out or sublimate. The molecular binding energy is thus an important parameter for many astrochemical studies. This parameter is usually determined with time-consuming experiments, computationally expensive quantum chemical calculations, or the inexpensive, but inaccurate, linear addition method. In this work we… ▽ More The behaviour of molecules in space is to a large extent governed by where they freeze out or sublimate. The molecular binding energy is thus an important parameter for many astrochemical studies. This parameter is usually determined with time-consuming experiments, computationally expensive quantum chemical calculations, or the inexpensive, but inaccurate, linear addition method. In this work we propose a new method based on machine learning for predicting binding energies that is accurate, yet computationally inexpensive. A machine learning model based on Gaussian Process Regression is created and trained on a database of binding energies of molecules collected from laboratory experiments presented in the literature. The molecules in the database are categorized by their features, such as mono- or multilayer coverage, binding surface, functional groups, valence electrons, and H-bond acceptors and donors. The performance of the model is assessed with five-fold and leave-one-molecule-out cross validation. Predictions are generally accurate, with differences between predicted and literature binding energies values of less than $\pm$20\%. The validated model is used to predict the binding energies of twenty one molecules that have recently been detected in the interstellar medium, but for which binding energy values are not known. A simplified model is used to visualize where the snowlines of these molecules would be located in a protoplanetary disk. This work demonstrates that machine learning can be employed to accurately and rapidly predict binding energies of molecules. Machine learning complements current laboratory experiments and quantum chemical computational studies. The predicted binding energies will find use in the modelling of astrochemical and planet-forming environments. △ Less

Submitted 8 July, 2022; originally announced July 2022.

Comments: Accepted in astronomy and astrophysics

Journal ref: A&A 666, A45 (2022)

Micromagnetics of magnetic chemical modulations in soft-magnetic cylindrical nanowires

Authors: L. Álvaro-Gómez, S. Ruiz-Gómez, C. Fernández-González, M. Schöbitz, N. Mille, J. Hurst, D. Tiwari, A. De Riz, I. M. Andersen, J. Bachmann, L. Cagnon, M. Foerster, L. Aballe, R. Belkhou, J. C Toussaint, C. Thirion, A. Masseboeuf, D. Gusakova, L. Pérez, O. Fruchart

Abstract: We analyze the micromagnetics of short longitudinal modulations of a high-magnetization material in cylindrical nanowires made of a soft-magnetic material of lower magnetization such as permalloy, combining magnetic microscopy, analytical modeling, and micromagnetic simulations. The mismatch of magnetization induces curling of magnetization around the axis in the modulations, in an attempt to scre… ▽ More We analyze the micromagnetics of short longitudinal modulations of a high-magnetization material in cylindrical nanowires made of a soft-magnetic material of lower magnetization such as permalloy, combining magnetic microscopy, analytical modeling, and micromagnetic simulations. The mismatch of magnetization induces curling of magnetization around the axis in the modulations, in an attempt to screen the interfacial magnetic charges. The curling angle increases with modulation length, until a plateau is reached with nearly full charge screening for a specific length scale~$Δ_\mathrm{mod}$, larger than the dipolar exchange length of any of the two materials. The curling circulation can be switched by the Oersted field arising from a charge current with typical magnitude $10^{12} A/m^{2}$ for a diameter of $\sim$100 nm, and reaching a maximum for $Δ_\mathrm{mod}$. △ Less

Submitted 6 September, 2022; v1 submitted 13 May, 2022; originally announced May 2022.

Journal ref: Phys. Rev. B 106, 054433 (2022)

Superconductivity in multiorbital systems with repulsive interactions: Hund's pairing vs. spin-fluctuation pairing

Authors: Mercè Roig, Astrid T. Rømer, Andreas Kreisel, P. J. Hirschfeld, Brian M. Andersen

Abstract: Hund's pairing refers to Cooper pairing generated by onsite interactions that become attractive due to large Hund's exchange $J$. This is possible in multiorbital systems even when all local bare interactions are repulsive, since attraction in specific channels are given by certain linear combinations of interaction parameters. On the other hand, pairing processes such as the exchange of spin fluc… ▽ More Hund's pairing refers to Cooper pairing generated by onsite interactions that become attractive due to large Hund's exchange $J$. This is possible in multiorbital systems even when all local bare interactions are repulsive, since attraction in specific channels are given by certain linear combinations of interaction parameters. On the other hand, pairing processes such as the exchange of spin fluctuations, are also present. We compare these two mechanisms on an equal footing using electronic bands appropriate for different classes of multiorbital systems over a wide range of interaction parameters. We find that for systems without clear nesting features, the superconducting state generated by the Hund's mechanism agrees well with that from the full fluctuation exchange vertex when Hund's exchange and spin-orbit coupling are sufficiently large. On the other hand, for systems characterized by a peaked finite-momentum particle-hole susceptibility, spin-fluctuation pairing generally dominates over Hund's pairing. We conclude that Hund's pairing states are unlikely to be realized in systems like Sr$_2$RuO$_4$ and generic iron-based superconductors. △ Less

Submitted 12 September, 2022; v1 submitted 7 April, 2022; originally announced April 2022.

Comments: 6 pages, 2 figures, including supplemental material

Journal ref: Phys. Rev. B 106, L100501 (2022)

Leading superconducting instabilities in three-dimensional models for Sr2RuO4

Authors: Astrid T. Rømer, T. A. Maier, Andreas Kreisel, P. J. Hirschfeld, Brian M. Andersen

Abstract: The superconductor Sr2RuO4 has been the subject of enormous interest over more than two decades, but until now the form of its order parameter has not been determined. Since groundbreaking NMR experiments revealed that the pairs are of dominant spin-singlet character, attention has focused on time-reversal symmetry breaking linear combinations of $s$-, $d$- and $g$-wave one-dimensional (1D) irredu… ▽ More The superconductor Sr2RuO4 has been the subject of enormous interest over more than two decades, but until now the form of its order parameter has not been determined. Since groundbreaking NMR experiments revealed that the pairs are of dominant spin-singlet character, attention has focused on time-reversal symmetry breaking linear combinations of $s$-, $d$- and $g$-wave one-dimensional (1D) irreducible representations. However, a state of the form $d_{xz}+id_{yz}$ has also been proposed. We present a systematic study of the stability of various superconducting candidate states, assuming that pairing is driven by the fluctuation exchange mechanism, including a realistic three-dimensional Fermi surface, full treatment of both local and non-local spin-orbit couplings, and a wide range of interaction parameters $U,J,U',J'$. The leading superconducting instabilities are found to exhibit nodal even-parity $A_{1g} (s')$ or $B_{1g} (d_{x^2-y^2})$ symmetries, similar to the findings in two-dimensional models without longer-range Coulomb interaction which tends to favor $d_{xy}$ over $d_{x^2-y^2}$. Within the so-called Hund's coupling mean-field pairing scenario, the $E_g (d_{xz}/d_{yz})$ solution can be stabilized for large $J$ and specific forms of the spin-orbit coupling, but for all cases studied here the eigenvalues of other superconducting solutions are significantly larger when the full fluctuation exchange vertex is included in the pairing kernel. Additionally, we compute the spin susceptibility in relevant superconducting candidate phases and compare to recent neutron scattering and NMR Knight shift measurements. △ Less

Submitted 6 April, 2022; originally announced April 2022.

Comments: 11 pages, 5 figures

Report number: NBI CMT 2022

Journal ref: Phys. Rev. Research 4, 033011 (2022)

A Framework for Improving the Reliability of Black-box Variational Inference

Authors: Manushi Welandawe, Michael Riis Andersen, Aki Vehtari, Jonathan H. Huggins

Abstract: Black-box variational inference (BBVI) now sees widespread use in machine learning and statistics as a fast yet flexible alternative to Markov chain Monte Carlo methods for approximate Bayesian inference. However, stochastic optimization methods for BBVI remain unreliable and require substantial expertise and hand-tuning to apply effectively. In this paper, we propose Robust and Automated Black-bo… ▽ More Black-box variational inference (BBVI) now sees widespread use in machine learning and statistics as a fast yet flexible alternative to Markov chain Monte Carlo methods for approximate Bayesian inference. However, stochastic optimization methods for BBVI remain unreliable and require substantial expertise and hand-tuning to apply effectively. In this paper, we propose Robust and Automated Black-box VI (RABVI), a framework for improving the reliability of BBVI optimization. RABVI is based on rigorously justified automation techniques, includes just a small number of intuitive tuning parameters, and detects inaccurate estimates of the optimal variational approximation. RABVI adaptively decreases the learning rate by detecting convergence of the fixed--learning-rate iterates, then estimates the symmetrized Kullback--Leibler (KL) divergence between the current variational approximation and the optimal one. It also employs a novel optimization termination criterion that enables the user to balance desired accuracy against computational cost by comparing (i) the predicted relative decrease in the symmetrized KL divergence if a smaller learning were used and (ii) the predicted computation required to converge with the smaller learning rate. We validate the robustness and accuracy of RABVI through carefully designed simulation studies and on a diverse set of real-world model and data examples. △ Less

Submitted 16 May, 2024; v1 submitted 29 March, 2022; originally announced March 2022.

Selectivity Trends and Role of Adsorbate-Adsorbate Interactions in CO Hydrogenation on Rhodium Catalysts

Authors: Martin Deimel, Hector Prats, Michael Seibt, Karsten Reuter, Mie Andersen

Abstract: Predictive-quality computational modeling of heterogeneously catalyzed reactions has emerged as an important tool for the analysis and assessment of activity and activity trends. In contrast, more subtle selectivities and selectivity trends still pose a significant challenge to prevalent microkinetic modeling approaches that typically employ a mean-field approximation (MFA). Here, we focus on CO h… ▽ More Predictive-quality computational modeling of heterogeneously catalyzed reactions has emerged as an important tool for the analysis and assessment of activity and activity trends. In contrast, more subtle selectivities and selectivity trends still pose a significant challenge to prevalent microkinetic modeling approaches that typically employ a mean-field approximation (MFA). Here, we focus on CO hydrogenation on Rh catalysts with the possible products methane, acetaldehyde, ethanol and water. This reaction has already been subject to a number of experimental and theoretical studies with conflicting views on the factors controlling activity and selectivity towards the more valuable higher oxygenates. Using accelerated first-principles kinetic Monte Carlo (KMC) simulations and explicitly and systematically accounting for adsorbate-adsorbate interactions through a cluster expansion approach, we model the reaction on the low-index Rh(111) and stepped Rh(211) surfaces. We find that the Rh(111) facet is selective towards methane, while the Rh(211) facet exhibits a similar selectivity towards methane and acetaldehyde. This is consistent with the experimental selectivity observed for larger, predominantly (111)-exposing Rh nanoparticles and resolves the discrepancy to earlier first-principles MFA microkinetic work that found the Rh(111) facet to be selective towards acetaldehyde. While the latter work tried to approximately account for lateral interactions through coverage-dependent rate expressions, our analysis demonstrates that this fails to sufficiently capture concomitant correlations among the adsorbed reaction intermediates that crucially determine the overall selectivity. △ Less

Submitted 29 March, 2022; originally announced March 2022.

Journal ref: ACS Catal. 2022, 12, 13, 7907-7917

Low-rank flat-field correction for artifact reduction in spectral computed tomography

Authors: Katrine O Bangsgaard, Genoveva Burca, Evelina Ametova, Martin S Andersen, Jakob S Jørgensen

Abstract: Spectral computed tomography has received considerable interest in recent years since spectral measurements contain much richer information about the object of interest. In spectral computed tomography, we are interested in the energy channel-wise reconstructions of the object. However, such reconstructions suffer from low signal-to-noise ratio and share the challenges of conventional low-dose com… ▽ More Spectral computed tomography has received considerable interest in recent years since spectral measurements contain much richer information about the object of interest. In spectral computed tomography, we are interested in the energy channel-wise reconstructions of the object. However, such reconstructions suffer from low signal-to-noise ratio and share the challenges of conventional low-dose computed tomography such as ring artifacts. Ring artifacts arise from errors in the flat-field correction and can significantly degrade the quality of the reconstruction. We propose an extended flat-field model that exploits high correlation in the spectral flat-fields to reduce ring artifacts in the channel-wise reconstructions. The extended model relies on the assumption that the spectral flat-fields can be well-approximated by a low-rank matrix. Our proposed model works directly on the spectral flat-fields and can be combined with any existing reconstruction model, e.g., filtered back projection and iterative methods. The proposed model is validated on a neutron data set. The results show that our method successfully diminishes ring artifacts and improves the quality of the reconstructions. Moreover, the results indicate that our method is robust; it only needs a single spectral flat-field image, whereas existing methods need multiple spectral flat-field images to reach a similar level of ring reduction. △ Less

Submitted 3 March, 2022; originally announced March 2022.

Comments: 14 pages, 8 figures

MSC Class: 65R32; 65C20; 94A08; 65K10 ACM Class: G.3; G.1.6

A Bayesian Approach to CT Reconstruction with Uncertain Geometry

Authors: Frederik H. Pedersen, Jakob S. Jørgensen, Martin S. Andersen

Abstract: Computed tomography is a method for synthesizing volumetric or cross-sectional images of an object from a collection of projections. Popular reconstruction methods for computed tomography are based on idealized models and assumptions that may not be valid in practice. One such assumption is that the exact projection geometry is known. The projection geometry describes the relative location of the… ▽ More Computed tomography is a method for synthesizing volumetric or cross-sectional images of an object from a collection of projections. Popular reconstruction methods for computed tomography are based on idealized models and assumptions that may not be valid in practice. One such assumption is that the exact projection geometry is known. The projection geometry describes the relative location of the radiation source, object, and detector for each projection. However, in practice, the geometric parameters used to describe the position and orientation of the radiation source, object, and detector are estimated quantities with uncertainty. A failure to accurately estimate the geometry may lead to reconstructions with severe misalignment artifacts that significantly decrease their scientific or diagnostic value. We propose a novel reconstruction method that jointly estimates the reconstruction and the projection geometry. The reconstruction method is based on a Bayesian approach that yields a point estimate for the reconstruction and geometric parameters and, in addition, provides valuable information regarding their uncertainty. This is achieved by approximately sampling from the joint posterior distribution of the reconstruction and projection geometry using a hierarchical Gibbs sampler. Using real tomographic data, we demonstrate that the proposed reconstruction method significantly reduces misalignment artifacts. Compared with two commonly used alignment methods, our proposed method achieves comparable or better results under challenging conditions. △ Less

Submitted 2 March, 2022; originally announced March 2022.

Comments: 15 pages, 6 figures

MSC Class: 65R32; 65C20; 94A08; 65K10 ACM Class: G.3; G.1.6

Predicting binding motifs of complex adsorbates using machine learning with a physics-inspired graph representation

Authors: Wenbin Xu, Karsten Reuter, Mie Andersen

Abstract: Computational screening in heterogeneous catalysis relies increasingly on machine learning models for predicting key input parameters due to the high cost of computing these directly using first-principles methods. This becomes especially relevant when considering complex materials spaces, e.g. alloys, or complex reaction mechanisms with adsorbates that may exhibit bi- or higher-dentate adsorption… ▽ More Computational screening in heterogeneous catalysis relies increasingly on machine learning models for predicting key input parameters due to the high cost of computing these directly using first-principles methods. This becomes especially relevant when considering complex materials spaces, e.g. alloys, or complex reaction mechanisms with adsorbates that may exhibit bi- or higher-dentate adsorption motifs. Here we present a data-efficient approach to the prediction of binding motifs and associated adsorption enthalpies of complex adsorbates at transition metals (TMs) and their alloys based on a customized Wasserstein Weisfeiler-Lehman graph kernel and Gaussian Process Regression. The model shows good predictive performance, not only for the elemental TMs on which it was trained, but also for an alloy based on these TMs. Furthermore, incorporation of minimal new training data allows for predicting an out-of-domain TM. We believe the model may be useful in active learning approaches, for which we present an ensemble uncertainty estimation approach. △ Less

Submitted 26 July, 2022; v1 submitted 23 February, 2022; originally announced February 2022.

Comments: update open source repository and reference

Superconducting Instabilities in Strongly-Correlated Infinite-Layer Nickelates

Authors: Andreas Kreisel, Brian M. Andersen, Astrid T. Rømer, Ilya M. Eremin, Frank Lechermann

Abstract: The discovery of superconductivity in infinite-layer nickelates has added a new family of materials to the fascinating growing class of unconventional superconductors. By incorporating the strongly correlated multi-orbital nature of the low-energy electronic degrees of freedom, we compute the leading superconducting instability from magnetic fluctuations relevant for infinite-layer nickelates. Spe… ▽ More The discovery of superconductivity in infinite-layer nickelates has added a new family of materials to the fascinating growing class of unconventional superconductors. By incorporating the strongly correlated multi-orbital nature of the low-energy electronic degrees of freedom, we compute the leading superconducting instability from magnetic fluctuations relevant for infinite-layer nickelates. Specifically, by properly including the doping dependence of the Ni $d_{x^2-y^2}$ and $d_{z^2}$ orbitals as well as the self-doping band, we uncover a transition from $d$-wave pairing symmetry to nodal $s_\pm$ superconductivity, driven by strong fluctuations in the $d_{z^2}$-dominated orbital states. We discuss the properties of the resulting superconducting condensates in light of recent tunneling and penetration depth experiments probing the detailed superconducting gap structure of these materials. △ Less

Submitted 12 August, 2022; v1 submitted 22 February, 2022; originally announced February 2022.

Comments: 7 pages, 4 figures, supplement included

Report number: NBI CMT 2022

Journal ref: Phys. Rev. Lett. 129, 077002 (2022)

Classical model for survival resonances close to Talbot time

Authors: Mikkel F. Andersen, Sandro Wimberger

Abstract: We present a classical approximation for the peaks of survival resonances occurring when diffracting matter waves from absorption potentials. Generally our simplified model describes the absorption-diffraction process around the Talbot time very well. Classical treatments of this process are presently lacking. For purely imaginary potentials, the classical model duplicates quantum mechanical calcu… ▽ More We present a classical approximation for the peaks of survival resonances occurring when diffracting matter waves from absorption potentials. Generally our simplified model describes the absorption-diffraction process around the Talbot time very well. Classical treatments of this process are presently lacking. For purely imaginary potentials, the classical model duplicates quantum mechanical calculations. The classical model allows for simple evolution of phase-space probability densities, which in the limit of the effective Planck's constant going to zero allows for a compact analytical expression of the survival probability as a function of remaining parameters. Our work extends the range of processes that can be described through classical analogues. △ Less

Submitted 10 February, 2022; originally announced February 2022.

Journal ref: Phys. Rev. A 105, 013322 (2022)

An Isolated Stellar-Mass Black Hole Detected Through Astrometric Microlensing

Authors: Kailash C. Sahu, Jay Anderson, Stefano Casertano, Howard E. Bond, Andrzej Udalski, Martin Dominik, Annalisa Calamida, Andrea Bellini, Thomas M. Brown, Marina Rejkuba, Varun Bajaj, Noe Kains, Henry C. Ferguson, Chris L. Fryer, Philip Yock, Przemek Mroz, Szymon Kozlowski, Pawel Pietrukowicz, Radek Poleski, Jan Skowron, Igor Soszynski, Michael K. Szymanski, Krzysztof Ulaczyk, Lukasz Wyrzykowski, Richard Barry , et al. (68 additional authors not shown)

Abstract: We report the first unambiguous detection and mass measurement of an isolated stellar-mass black hole (BH). We used the Hubble Space Telescope (HST) to carry out precise astrometry of the source star of the long-duration (t_E~270 days), high-magnification microlensing event MOA-2011-BLG-191/OGLE-2011-BLG-0462 (hereafter designated as MOA-11-191/OGLE-11-462), in the direction of the Galactic bulge.… ▽ More We report the first unambiguous detection and mass measurement of an isolated stellar-mass black hole (BH). We used the Hubble Space Telescope (HST) to carry out precise astrometry of the source star of the long-duration (t_E~270 days), high-magnification microlensing event MOA-2011-BLG-191/OGLE-2011-BLG-0462 (hereafter designated as MOA-11-191/OGLE-11-462), in the direction of the Galactic bulge. HST imaging, conducted at eight epochs over an interval of six years, reveals a clear relativistic astrometric deflection of the background star's apparent position. Ground-based photometry of MOA-11-191/OGLE-11-462 shows a parallactic signature of the effect of the Earth's motion on the microlensing light curve. Combining the HST astrometry with the ground-based light curve and the derived parallax, we obtain a lens mass of 7.1 +/- 1.3 Msun and a distance of 1.58 +/- 0.18 kpc. We show that the lens emits no detectable light, which, along with having a mass higher than is possible for a white dwarf or neutron star, confirms its BH nature. Our analysis also provides an absolute proper motion for the BH. The proper motion is offset from the mean motion of Galactic-disk stars at similar distances by an amount corresponding to a transverse space velocity of ~45 km/s, suggesting that the BH received a 'natal kick' from its supernova explosion. Previous mass determinations for stellar-mass BHs have come from radial-velocity measurements of Galactic X-ray binaries, and from gravitational radiation emitted by merging BHs in binary systems in external galaxies. Our mass measurement is the first for an isolated stellar-mass BH using any technique. △ Less

Submitted 22 July, 2022; v1 submitted 31 January, 2022; originally announced January 2022.

Comments: 37 pages, Published in ApJ

Journal ref: ApJ, 933, 83 (2022)

Theory of spin-excitation anisotropy in the nematic phase of FeSe obtained from RIXS measurements

Authors: Andreas Kreisel, P. J. Hirschfeld, Brian M. Andersen

Abstract: Recent resonant inelastic x-ray scattering (RIXS) experiments have detected a significant high-energy spin-excitation anisotropy in the nematic phase of the enigmatic iron-based superconductor FeSe, whose origin remains controversial. We apply an itinerant model previously used to describe the spin-excitation anisotropy as measured by neutron scattering measurements, with magnetic fluctuations inc… ▽ More Recent resonant inelastic x-ray scattering (RIXS) experiments have detected a significant high-energy spin-excitation anisotropy in the nematic phase of the enigmatic iron-based superconductor FeSe, whose origin remains controversial. We apply an itinerant model previously used to describe the spin-excitation anisotropy as measured by neutron scattering measurements, with magnetic fluctuations included within the RPA approximation. The calculated RIXS cross section exhibits overall agreement with the RIXS data, including the high energy spin-excitation anisotropy. △ Less

Submitted 22 April, 2022; v1 submitted 28 January, 2022; originally announced January 2022.

Comments: 14 pages, 6 figures, published version

Report number: CMT NBI 2022

Journal ref: Frontiers in Physics 10, 859424 (2022)

Towards a secure API client generator for IoT devices

Authors: Anders Aaen Springborg, Martin Kaldahl Andersen, Kaare Holland Hattel, Michele Albano

Abstract: Given the success of IoT platforms, more developers and companies want to include the technology in their portfolio. However, in the case of single board microcontrollers, the support for networking operations is not ideal, and different IoT platforms allow access to the networking submodule via different libraries and system calls, leading to a steeper learning curve. Code generators for API clie… ▽ More Given the success of IoT platforms, more developers and companies want to include the technology in their portfolio. However, in the case of single board microcontrollers, the support for networking operations is not ideal, and different IoT platforms allow access to the networking submodule via different libraries and system calls, leading to a steeper learning curve. Code generators for API clients can enhance productivity, but they tend to generate universal purpose code, and on the other hand the networking primitives of IoT devices are platform specific, especially when security mechanisms such as Transport Layer Security are part of the picture. This paper presents \texttt{cpp-tiny-client}, an API client generator developed as a plugin for the OpenAPI Generator project, which can tailor the generated code based on the IoT platform specified by the user. Our work allows to generate correct code for API clients for IoT devices, and thus can empower a developer with more productivity and a faster time-to-market for its own applications. By combining together mainstream technologies only, \texttt{cpp-tiny-client} offers a gentle learning curve. Moreover, experiments show that the generated code has a reasonable footprint, at least with respect to the IoT devices that were used in the validation of the work. The code related to this work is available through the OpenAPI Generator project~\cite{OpenAPIGenerator}. This technical report is an extension of~\cite{acmsac22}, and it integrates the information presented at the ACM SAC 2022 conference. △ Less

Submitted 1 January, 2022; originally announced January 2022.

Comments: The work was accepted as 4-pages paper to ACM SAC 2022. Together with the TPC, I am providing a technical report here with the information I have to cut from the 4-pages version

ACM Class: C.2; D.2

Measured potential profile in a quantum anomalous Hall system suggests bulk-dominated current flow

Authors: Ilan T. Rosen, Molly P. Andersen, Linsey K. Rodenbach, Lixuan Tai, Peng Zhang, Kang L. Wang, M. A. Kastner, David Goldhaber-Gordon

Abstract: Ideally, quantum anomalous Hall systems should display zero longitudinal resistance. Yet in experimental quantum anomalous Hall systems elevated temperature can make the longitudinal resistance finite, indicating dissipative flow of electrons. Here, we show that the measured potentials at multiple locations within a device at elevated temperature are well-described by solution of Laplace's equatio… ▽ More Ideally, quantum anomalous Hall systems should display zero longitudinal resistance. Yet in experimental quantum anomalous Hall systems elevated temperature can make the longitudinal resistance finite, indicating dissipative flow of electrons. Here, we show that the measured potentials at multiple locations within a device at elevated temperature are well-described by solution of Laplace's equation, assuming spatially-uniform conductivity, suggesting non-equilibrium current flows through the two-dimensional bulk. Extrapolation suggests that at even lower temperatures current may still flow primarily through the bulk rather than, as had been assumed, through edge modes. An argument for bulk current flow previously applied to quantum Hall systems supports this picture. △ Less

Submitted 7 November, 2022; v1 submitted 24 December, 2021; originally announced December 2021.

Comments: 6 pages, 4 figures, plus supplemental materials

Report number: Phys. Rev. Lett. 129, 246602 (2022)

Evaluation of diagnostic test procedures for SARS-CoV-2 using latent class models: comparison of antigen test kits and sampling for PCR testing based on Danish national data registries

Authors: Jacob Stærk-Østergaard, Carsten Kirkeby, Lasse Engbo Christiansen, Michael Asger Andersen, Camilla Holten Møller, Marianne Voldstedlund, Matthew J. Denwood

Abstract: Antigen test kits have been used extensively as a screening tool during the worldwide pandemic of coronavirus (SARS-CoV-2). While it is generally expected that taking samples for analysis with PCR testing gives more reliable results than using antigen test kits, the overall sensitivity and specificity of the two protocols in the field have not yet been estimated without assuming that the PCR test… ▽ More Antigen test kits have been used extensively as a screening tool during the worldwide pandemic of coronavirus (SARS-CoV-2). While it is generally expected that taking samples for analysis with PCR testing gives more reliable results than using antigen test kits, the overall sensitivity and specificity of the two protocols in the field have not yet been estimated without assuming that the PCR test constitutes a gold standard. We use latent class models to estimate the in situ performance of both PCR and antigen testing, using data from the Danish national registries. The results are based on 240,000 paired tests results sub-selected from the 55 million test results that were obtained in Denmark during the period from February 2021 until June 2021. We found that the specificity of both tests is very high in our data sample (>99.7%), while the sensitivity of PCR sampling was estimated to be 95.7% (95% CI: 92.8-98.4%) and that of the antigen test kits used in Denmark over the study period was estimated at 53.8% (95% CI: 49.8-57.9%). Our findings can be used as supplementary information for consideration when implementing serial testing strategies that employ a confirmatory PCR sample following a positive result from an antigen test kit, such as the policy used in Denmark. We note that while this strategy reduces the number of false positives associated with antigen test screening, it also increases the false negatives. We demonstrate that the balance of trading false positives for false negatives only favours the use of serial testing when the expected true prevalence is low. Our results contain substantial uncertainty in the estimates for sensitivity due to the relatively small number of positive test results over this period: validation of our findings in a population with higher prevalence would therefore be highly relevant for future work. △ Less

Submitted 21 December, 2021; originally announced December 2021.

30.000 ways to reach 55% decarbonization of the European electricity sector

Authors: Tim T. Pedersen, Mikael Skou Andersen, Marta Victoria, Gorm B. Andresen

Abstract: Climate change mitigation is a global challenge that, however, needs to be resolved by national-level authorities, resembling a tragedy of the commons. This paradox is reflected at European scale, as climate commitments are made by the EU collectively, but implementation is the responsibility of individual Member States. Here, we investigate 30.000 near-optimal effort-sharing scenarios where the E… ▽ More Climate change mitigation is a global challenge that, however, needs to be resolved by national-level authorities, resembling a tragedy of the commons. This paradox is reflected at European scale, as climate commitments are made by the EU collectively, but implementation is the responsibility of individual Member States. Here, we investigate 30.000 near-optimal effort-sharing scenarios where the European electricity sector is decarbonized by at least 55% relative to 1990, in line with 2030 ambitions. Using a highly detailed brownfield electricity system optimization model, the optimal electricity system is simulated for a suite of effort-sharing scenarios. Results reveal large inequalities in the efforts required to decarbonize national electricity sectors, with some countries facing cost-optimal pathways to reach 55% emission reductions, while others are confronted with relatively high abatement costs. Specifically, we find that several countries with modest or low levels of GDP per capita will experience high abatement costs, and when passed over into electricity prices this may lead to increased energy poverty in certain parts of Europe △ Less

Submitted 30 November, 2022; v1 submitted 14 December, 2021; originally announced December 2021.

Finite rank kernel varieties: A variant of Hilbert's Nullstellensatz for graphons and applications to Hadamard matrices

Authors: Madelyn Andersen

Abstract: Graphons are symmetric measurable functions that arise from a sequence of graphs. A graphon variety is the a set of all graphons defined by a condition of the form $t(g, W) = 0$ for a fixed quantum graph $g$, where $t(.,.)$ is the homomorphism density and a quantum graph is a formal linear combination of multigraphs. Using a method of representing graphs as polynomials, we construct an epimorphism… ▽ More Graphons are symmetric measurable functions that arise from a sequence of graphs. A graphon variety is the a set of all graphons defined by a condition of the form $t(g, W) = 0$ for a fixed quantum graph $g$, where $t(.,.)$ is the homomorphism density and a quantum graph is a formal linear combination of multigraphs. Using a method of representing graphs as polynomials, we construct an epimorphism from the space of quantum graphs to a subring of the complex polynomial ring that is invariant under permutations of variables. When graphons are of finite rank, we demonstrate that an analog of the "ideal" inverse in Algebraic Geometry is an ideal in our polynomial representation. Defining an algebraic kernel set using kernel varieties, we demonstrate that we can call such sets closed under the Zariski Topology. We determine several ties to Algebraic Geometry as a result of utilizing finite rank kernels and discover that a weaker version of Hilbert's Nullstellensatz applies to kernel zero-sets with respect to homomorphism density. Throughout, we examine the connection between Algebraic Geometry and Graphon Theory. △ Less

Submitted 12 December, 2021; v1 submitted 7 December, 2021; originally announced December 2021.

Orbital-Selective High-Temperature Cooper Pairing Developed in the Two-Dimensional Limit

Authors: Chaofei Liu, Andreas Kreisel, Shan Zhong, Yu Li, Brian M. Andersen, P. J. Hirschfeld, Jian Wang

Abstract: The orbital multiplicity in multiband superconductors yields orbital differentiation in normal-state properties, and can lead to orbital-selective spin-fluctuation Cooper pairing. This phenomenon has become increasingly pivotal in clarifying the pairing 'enigma' particularly for multiband high-temperature superconductors. In one-unit-cell (1-UC) FeSe/SrTiO3, the thinnest and highest-Tc member of i… ▽ More The orbital multiplicity in multiband superconductors yields orbital differentiation in normal-state properties, and can lead to orbital-selective spin-fluctuation Cooper pairing. This phenomenon has become increasingly pivotal in clarifying the pairing 'enigma' particularly for multiband high-temperature superconductors. In one-unit-cell (1-UC) FeSe/SrTiO3, the thinnest and highest-Tc member of iron-based superconductors, the standard electron-hole Fermi pocket nesting scenario is apparently not applicable since the Gamma-centered hole pockets are absent, so the actual pairing mechanism is the subject of intense debate. Here, by measuring high-resolution Bogoliubov quasiparticle interference, we report observations of highly anisotropic magnetic Cooper pairing in 1-UC FeSe. From a theoretical point of view, it is important to incorporate effects of electronic correlations within a spin-fluctuation pairing calculation, where the dxy orbital becomes coherence-suppressed. The resulting pairing gap is compatible with the experimental findings, which suggests that high-Tc Cooper pairing with orbital selectivity applies to 1-UC FeSe. Our findings imply the general existence of orbital selectivity in iron-based superconductors and the universal importance of electron correlations in high-Tc superconductors. △ Less

Submitted 1 December, 2021; originally announced December 2021.

Journal ref: Nano Letters 22, 3245-3251 (2022)

Optical Vortex Braiding with Bessel Beams

Authors: Andrew A. Voitiv, Jasmine M. Andersen, Mark E. Siemens, Mark T. Lusk

Abstract: We propose the braiding of optical vortices in a laser beam with more than 2π rotation by superposing Bessel modes with a plane wave. We experimentally demonstrate this by using a Bessel-Gaussian beam and a coaxial Gaussian, and we present measurement of three complete braids. The amount of braiding is fundamentally limited only by the numerical aperture of the system and we discuss how braiding c… ▽ More We propose the braiding of optical vortices in a laser beam with more than 2π rotation by superposing Bessel modes with a plane wave. We experimentally demonstrate this by using a Bessel-Gaussian beam and a coaxial Gaussian, and we present measurement of three complete braids. The amount of braiding is fundamentally limited only by the numerical aperture of the system and we discuss how braiding can be controlled experimentally for any number of vortices. △ Less

Submitted 19 November, 2021; originally announced November 2021.

Journal ref: Optics Letters Vol. 45 No. 6 (2020), 1321-1324

Hydrodynamics Explanation for the Splitting of Higher-charge Optical Vortices

Authors: Andrew A. Voitiv, Jasmine M. Andersen, Patrick C. Ford, Mark T. Lusk, Mark E. Siemens

Abstract: We show that a two-dimensional hydrodynamics model provides a physical explanation for the splitting of higher-charge optical vortices under elliptical deformations. The model is applicable to laser light and quantum fluids alike. The study delineates vortex breakups from vortex unions under different forms of asymmetry in the beam, and it is also applied to explain the motion of intact higher-cha… ▽ More We show that a two-dimensional hydrodynamics model provides a physical explanation for the splitting of higher-charge optical vortices under elliptical deformations. The model is applicable to laser light and quantum fluids alike. The study delineates vortex breakups from vortex unions under different forms of asymmetry in the beam, and it is also applied to explain the motion of intact higher-charge vortices. △ Less

Submitted 9 March, 2022; v1 submitted 16 November, 2021; originally announced November 2021.

Journal ref: Optics Letters Vol. 47 No. 6 pp. 1391-1394 (2022)

Amplitude Structure of Optical Vortices Determines Annihilation Dynamics

Authors: Jasmine M. Andersen, Andrew A. Voitiv, Patrick C. Ford, Mark E. Siemens

Abstract: We show that annihilation dynamics between oppositely charged optical vortex pairs can be manipulated by the initial size of the vortex cores, consistent with hydrodynamics. When sufficiently close together, vortices with strongly overlapped cores annihilate more quickly than vortices with smaller cores that must wait for diffraction to cause meaningful core overlap. Numerical simulations and expe… ▽ More We show that annihilation dynamics between oppositely charged optical vortex pairs can be manipulated by the initial size of the vortex cores, consistent with hydrodynamics. When sufficiently close together, vortices with strongly overlapped cores annihilate more quickly than vortices with smaller cores that must wait for diffraction to cause meaningful core overlap. Numerical simulations and experimental measurements for vortices with hyperbolic tangent cores of various initial sizes show that hydrodynamics governs their motion, and reveal distinct phases of vortex recombination; decreasing the core size of an annihilating pair can prevent the annihilation event. △ Less

Submitted 5 September, 2023; v1 submitted 4 November, 2021; originally announced November 2021.

Journal ref: J. Opt. Soc. Am. A 40, 223-228 (2023)

No swan song for Sun-as-a-star helioseismology: performances of the Solar-SONG prototype for individual mode characterisation

Authors: S. N. Breton, P. L. Pallé, R. A. García, M. Fredslund Andersen, F. Grundahl, J. Christensen-Dalsgaard, H. Kjeldsen, S. Mathur

Abstract: The GOLF instrument on board SoHO has been in operation for almost 25 years but aging of the instrument has now strongly affected its performance, especially in the low-frequency p-mode region. At the end of the SoHO mission, the ground-based network BiSON will remain the only facility able to perform Sun-integrated helioseismic observations. Therefore, we want to assess the helioseismic performan… ▽ More The GOLF instrument on board SoHO has been in operation for almost 25 years but aging of the instrument has now strongly affected its performance, especially in the low-frequency p-mode region. At the end of the SoHO mission, the ground-based network BiSON will remain the only facility able to perform Sun-integrated helioseismic observations. Therefore, we want to assess the helioseismic performances of an échelle spectrograph like SONG. Indeed, the high precision of such an instrument and the quality of the data acquired for asteroseismic purpose calls for an evaluation of the instrument ability to perform global radial-velocity measurements of the solar disk. Data acquired during the Solar-SONG 2018 observation campaign at the Teide Observatory are used to study mid- and low-frequency p modes. A Solar-SONG time series of 30-day duration is reduced with a combination of the traditional IDL iSONG pipeline and a new Python pipeline described in this paper. A mode fitting method built around a Bayesian approach is then performed on the Solar-SONG and contemporaneous GOLF, BiSON, and HMI data. For this contemporaneous time series, Solar-SONG is able to characterise p modes at a lower frequency than BiSON and GOLF (1750μHz against 1946 and 2157 μHz respectively), while for HMI it is possible to characterise a mode at 1686 μHz. The decrease of GOLF sensitivity is then evaluated through the evolution of its low-frequency p-mode characterisation abilities over the years. [abridged] △ Less

Submitted 25 October, 2021; originally announced October 2021.

Comments: 24 pages, 17 figures. Accepted in A&A

Journal ref: A&A 658, A27 (2022)

Fractional AC Josephson effect in a topological insulator proximitized by a self-formed superconductor

Authors: Ilan T. Rosen, Christie J. Trimble, Molly P. Andersen, Evgeny Mikheev, Yanbin Li, Yunzhi Liu, Lixuan Tai, Peng Zhang, Kang L. Wang, Yi Cui, M. A. Kastner, James R. Williams, David Goldhaber-Gordon

Abstract: A lateral Josephson junction in which the surface of a 3D topological insulator serves as the weak link should support topologically protected excitations related to Majorana fermions. The resulting $4π$-periodic current-phase relationship could be detected under high-frequency excitation by the suppression of odd Shapiro steps. Here, we demonstrate such devices through the self-formation of a Pd-… ▽ More A lateral Josephson junction in which the surface of a 3D topological insulator serves as the weak link should support topologically protected excitations related to Majorana fermions. The resulting $4π$-periodic current-phase relationship could be detected under high-frequency excitation by the suppression of odd Shapiro steps. Here, we demonstrate such devices through the self-formation of a Pd-Te superconducting layer from a telluride topological insulator, and observe suppressed first and third Shapiro steps. Yet many of our devices, including ones where the Pd-Te layer is bolstered by an additional Al layer, show no suppression of Shapiro steps. Though we rule out the known trivial causes of suppressed Shapiro steps in our devices, our results serve as a warning that the non-equilibrium physics of these systems must be better understood before confidently claiming the observation of Majorana states. △ Less

Submitted 3 October, 2021; originally announced October 2021.

Comments: 8 pages, 5 figures

Supercurrents and spontaneous time-reversal symmetry breaking by nonmagnetic disorder in unconventional superconductors

Authors: Clara N. Breiø, P. J. Hirschfeld, Brian M. Andersen

Abstract: Recently, a theoretical study [Li {\it et al.},~npj Quantum Materials 6, 36 (2021)] investigated a model of a disordered $d$-wave superconductor, and reported local time-reversal symmetry breaking current loops for sufficiently high disorder levels. Since the pure $d$-wave superconducting state does not break time-reversal symmetry, it is surprising that such persistent currents arise purely from… ▽ More Recently, a theoretical study [Li {\it et al.},~npj Quantum Materials 6, 36 (2021)] investigated a model of a disordered $d$-wave superconductor, and reported local time-reversal symmetry breaking current loops for sufficiently high disorder levels. Since the pure $d$-wave superconducting state does not break time-reversal symmetry, it is surprising that such persistent currents arise purely from nonmagnetic disorder. Here we perform a detailed theoretical investigation of such disorder-induced orbital currents, and show that the occurrence of the currents can be traced to local extended $s$-wave pairing. We discuss the energetics leading to regions of $s\pm id$ order, which support spontaneous local currents in the presence of inhomogeneous density modulations. △ Less

Submitted 28 September, 2021; originally announced September 2021.

Comments: 8 pages, 5 figures

Report number: CMT NBI 2021

Journal ref: Phys. Rev. B 105, 014504 (2022)

Superconducting gap symmetry from Bogoliubov quasiparticle interference analysis on {Sr}$_2${RuO}$_4$

Authors: Shinibali Bhattacharyya, Andreas Kreisel, X. Kong, T. Berlijn, Astrid T. Rømer, Brian M. Andersen, P. J. Hirschfeld

Abstract: The nature of the superconducting order parameter in {Sr}$_2${RuO}$_4$ has generated intense interest in recent years. Since the superconducting gap is very small, high resolution methods such as scanning tunneling spectroscopy might be the best chance to directly resolve the gap symmetry. Recently, a Bogoliubov quasiparticle interference imaging (BQPI) experiment has suggested that the… ▽ More The nature of the superconducting order parameter in {Sr}$_2${RuO}$_4$ has generated intense interest in recent years. Since the superconducting gap is very small, high resolution methods such as scanning tunneling spectroscopy might be the best chance to directly resolve the gap symmetry. Recently, a Bogoliubov quasiparticle interference imaging (BQPI) experiment has suggested that the $d_{x^2-y^2}$ gap symmetry is appropriate for {Sr}$_2${RuO}$_4$. In this work, we use a material-specific theoretical approach based on Wannier functions of the surface of {Sr}$_2${RuO}$_4$ to calculate the continuum density of states as detected in scanning tunneling microscopy experiments. We examine several different proposed gap order parameters, and calculate the expected BQPI pattern for each case. Comparing to the available experimental data, our results suggest that a $s'+id_{xy}$ gap order parameter is the most probable state, but the measured BQPI patterns still display features unaccounted for by the theory for any of the states currently under discussion. △ Less

Submitted 18 June, 2023; v1 submitted 22 September, 2021; originally announced September 2021.

Comments: 14 pages, 9 figures

Journal ref: Phys. Rev. B 107, 144505 (2023)

Micro-architectural Analysis of a Learned Index

Authors: Mikkel Møller Andersen, Pınar Tözün

Abstract: Since the publication of The Case for Learned Index Structures in 2018, there has been a rise in research that focuses on learned indexes for different domains and with different functionalities. While the effectiveness of learned indexes as an alternative to traditional index structures such as B+Trees have already been demonstrated by several studies, previous work tend to focus on higher-level… ▽ More Since the publication of The Case for Learned Index Structures in 2018, there has been a rise in research that focuses on learned indexes for different domains and with different functionalities. While the effectiveness of learned indexes as an alternative to traditional index structures such as B+Trees have already been demonstrated by several studies, previous work tend to focus on higher-level performance metrics such as throughput and index size. In this paper, our goal is to dig deeper and investigate how learned indexes behave at a micro-architectural level compared to traditional indexes. More specifically, we focus on previously proposed learned index structure ALEX, which is a tree-based in-memory index structure that consists of a hierarchy of machine learned models. Unlike the original proposal for learned indexes, ALEX is designed from the ground up to allow updates and inserts. Therefore, it enables more dynamic workloads using learned indexes. In this work, we perform a micro-architectural analysis of ALEX and compare its behavior to the tree-based index structures that are not based on learned models, i.e., ART and B+Tree. Our results show that ALEX is bound by memory stalls, mainly stalls due to data misses from the last-level cache. Compared to ART and B+Tree, ALEX exhibits fewer stalls and a lower cycles-per-instruction value across different workloads. On the other hand, the amount of instructions required to handle out-of-bound inserts in ALEX can increase the instructions needed per request significantly (10X) for write-heavy workloads. However, the micro-architectural behavior shows that this increase in the instruction footprint exhibit high instruction-level parallelism, and, therefore, does not negatively impact the overall execution time. △ Less

Submitted 17 September, 2021; originally announced September 2021.

Comments: Under submission

Field tunable three-dimensional magnetic nanotextures in cobalt-nickel nanowires

Authors: I. M. Andersen, D. Wolf, L. A. Rodriguez, A. Lubk, D. Oliveros, C. Bran, T. Niermann, U. K. Rößler, M. Vazquez, C. Gatel, E. Snoeck

Abstract: Cylindrical magnetic nanowires with large transversal magnetocrystalline anisotropy have been shown to sustain non-trivial magnetic configurations resulting from the interplay of spatial confinement, exchange, and anisotropies. Exploiting these peculiar 3D spin configurations and their solitonic inhomogeneities are prospected to improve magnetization switching in future spintronics, such as power-… ▽ More Cylindrical magnetic nanowires with large transversal magnetocrystalline anisotropy have been shown to sustain non-trivial magnetic configurations resulting from the interplay of spatial confinement, exchange, and anisotropies. Exploiting these peculiar 3D spin configurations and their solitonic inhomogeneities are prospected to improve magnetization switching in future spintronics, such as power-saving magnetic memory and logic applications. Here we employ holographic vector field electron tomography to reconstruct the remanent magnetic states in CoNi nanowires with 10 nm resolution in 3D, with a particular focus on domain walls between remanent states and ubiquitous real-structure effects stemming from irregular morphology and anisotropy variations. By tuning the applied magnetic field direction, both longitudinal and transverse multi-vortex states of different chiralities and peculiar 3D features such as shifted vortex cores are stabilized. The chiral domain wall between the longitudinal vortices of opposite chiralities exhibits a complex 3D shape characterized by a push out of the central vortex line and a gain in exchange and anisotropy energy. A similar complex 3D texture, including bent vortex lines, forms at the domain boundary between transverse-vortex states and longitudinal configurations. Micromagnetic simulations allow an understanding of the origin of the observed complex magnetic states. △ Less

Submitted 28 July, 2021; originally announced July 2021.

Comments: 22 pages, 5 figures, preprent of original paper

Journal ref: Phys. Rev. Res. 3, (2021) 033085

Theory of the charge-density wave in $A$V$_3$Sb$_5$ kagome metals

Authors: Morten H. Christensen, Turan Birol, Brian M. Andersen, Rafael M. Fernandes

Abstract: The family of metallic kagome compounds $A$V$_3$Sb$_5$ ($A$=K, Rb, Cs) was recently discovered to exhibit both superconductivity and charge order. The nature of the charge-density wave (CDW) phase is presently unsettled, which complicates the interpretation of the superconducting ground state. In this paper, we use group-theory and density-functional theory (DFT) to derive and solve a phenomenolog… ▽ More The family of metallic kagome compounds $A$V$_3$Sb$_5$ ($A$=K, Rb, Cs) was recently discovered to exhibit both superconductivity and charge order. The nature of the charge-density wave (CDW) phase is presently unsettled, which complicates the interpretation of the superconducting ground state. In this paper, we use group-theory and density-functional theory (DFT) to derive and solve a phenomenological Landau model for this CDW state. The DFT results reveal three unstable phonon modes with the same in-plane momentum but different out-of-plane momenta, whose frequencies depend strongly on the electronic temperature. This is indicative of an electronically-driven CDW, stabilized by features of the in-plane electronic dispersion. Motivated by the DFT analysis, we construct a Landau free-energy expansion for coupled CDW order parameters with wave-vectors at the $M$ and $L$ points of the hexagonal Brillouin zone. We find an unusual trilinear term coupling these different order parameters, which can promote the simultaneous condensation of both CDWs even if the two modes are not nearly-degenerate. We classify the different types of coupled multi-$\bf{Q}$ CDW orders, focusing on those that break the sixfold rotational symmetry and lead to a unit-cell doubling along all three crystallographic directions, as suggested by experiments. We determine a region in parameter space, characterized by large nonlinear Landau coefficients, where these phases - dubbed staggered tri-hexagonal and staggered Star-of-David - are the leading instabilities of the system. Finally, we discuss the implications of our results for the kagome metals. △ Less

Submitted 12 July, 2021; v1 submitted 9 July, 2021; originally announced July 2021.

Comments: 17 pages, 11 figures

Report number: NBI CMT 2021

Journal ref: Phys. Rev. B 104, 214513 (2021)

Revealing the competition between charge-density wave and superconductivity in CsV$_3$Sb$_5$ through uniaxial strain

Authors: Tiema Qian, Morten H. Christensen, Chaowei Hu, Amartyajyoti Saha, Brian M. Andersen, Rafael M. Fernandes, Turan Birol, Ni Ni

Abstract: In this paper we report the impact of uniaxial strain $\varepsilon$ applied along the crystalline $a$ axis on the newly discovered kagome superconductor CsV$_3$Sb$_5$. At ambient conditions, CsV$_3$Sb$_5$ shows a charge-density wave (CDW) transition at $T_{\rm CDW}=94.5$ K and superconducts below $T_C = 3.34$ K. In our study, when the uniaxial strain $\varepsilon$ is varied from $-0.90\%$ to… ▽ More In this paper we report the impact of uniaxial strain $\varepsilon$ applied along the crystalline $a$ axis on the newly discovered kagome superconductor CsV$_3$Sb$_5$. At ambient conditions, CsV$_3$Sb$_5$ shows a charge-density wave (CDW) transition at $T_{\rm CDW}=94.5$ K and superconducts below $T_C = 3.34$ K. In our study, when the uniaxial strain $\varepsilon$ is varied from $-0.90\%$ to $0.90\%$, $T_C$ monotonically increases by $\sim 33\%$ from 3.0 K to 4.0 K, giving rise to the empirical relation $T_C (\varepsilon)=3.4+0.56\varepsilon+0.12\varepsilon^2$. On the other hand, for $\varepsilon$ changing from $-0.76\%$ to $1.26\%$, $T_{\rm CDW}$ decreases monotonically by $\sim 10\%$ from 97.5 K to 87.5 K with $T_{\rm CDW}(\varepsilon)=94.5-4.72\varepsilon-0.60\varepsilon^2$. The opposite response of $T_C$ and $T_{\rm CDW}$ to the uniaxial strain suggests strong competition between these two orders. Comparison with hydrostatic pressure measurements indicate that it is the change in the $c$-axis that is responsible for these behaviors of the CDW and superconducting transitions, and that the explicit breaking of the sixfold rotational symmetry by strain has a negligible effect. Combined with our first-principles calculations and phenomenological analysis, we conclude that the enhancement in $T_C$ with decreasing $c$ is caused primarily by the suppression of $T_{\rm CDW}$, rather than strain-induced modifications in the bare superconducting parameters. We propose that the sensitivity of $T_{\rm CDW}$ with respect to the changes in the $c$-axis arises from the impact of the latter on the trilinear coupling between the $M_1^+$ and $L_2^-$ phonon modes associated with the CDW. Overall, our work reveals that the $c$-axis lattice parameter, which can be controlled by both pressure and uniaxial strain, is a powerful tuning knob for the phase diagram of CsV$_3$Sb$_5$. △ Less

Submitted 12 July, 2021; v1 submitted 9 July, 2021; originally announced July 2021.

Comments: 8 pages, 5 figures

Journal ref: Phys. Rev. B 104 (2021), 144506

Assessing the forensic value of DNA evidence from Y chromosomes and mitogenomes

Authors: Mikkel M Andersen, David J Balding

Abstract: Y-chromosomal and mitochondrial DNA profiles have been used as evidence in courts for decades, yet the problem of evaluating the weight of evidence has not been adequately resolved. Both are lineage markers (inherited from just one parent), which presents different interpretation challenges compared with standard autosomal DNA profiles (inherited from both parents), for which recombination increas… ▽ More Y-chromosomal and mitochondrial DNA profiles have been used as evidence in courts for decades, yet the problem of evaluating the weight of evidence has not been adequately resolved. Both are lineage markers (inherited from just one parent), which presents different interpretation challenges compared with standard autosomal DNA profiles (inherited from both parents), for which recombination increases profile diversity and weakens the effects of relatedness. We review approaches to the evaluation of lineage marker profiles for forensic identification, focussing on the key roles of profile mutation rate and relatedness. Higher mutation rates imply fewer individuals matching the profile of an alleged contributor, but they will be more closely related. This makes it challenging to evaluate the possibility that one of these matching individuals could be the true source, because relatedness may make them more plausible alternative contributors than less-related individuals, and they may not be well mixed in the population. These issues reduce the usefulness of profile databases drawn from a broad population: the larger the population, the lower the profile relative frequency because of lower relatedness with the alleged contributor. Many evaluation methods do not adequately take account of relatedness, but its effects have become more pronounced with the latest generation of high-mutation-rate Y profiles. △ Less

Submitted 7 July, 2021; originally announced July 2021.

Mr. Plotter: Unifying Data Reduction Techniques in Storage and Visualization Systems

Authors: Sam Kumar, Michael P Andersen, David E. Culler

Abstract: As the rate of data collection continues to grow rapidly, developing visualization tools that scale to immense data sets is a serious and ever-increasing challenge. Existing approaches generally seek to decouple storage and visualization systems, performing just-in-time data reduction to transparently avoid overloading the visualizer. We present a new architecture in which the visualizer and data… ▽ More As the rate of data collection continues to grow rapidly, developing visualization tools that scale to immense data sets is a serious and ever-increasing challenge. Existing approaches generally seek to decouple storage and visualization systems, performing just-in-time data reduction to transparently avoid overloading the visualizer. We present a new architecture in which the visualizer and data store are tightly coupled. Unlike systems that read raw data from storage, the performance of our system scales linearly with the size of the final visualization, essentially independent of the size of the data. Thus, it scales to massive data sets while supporting interactive performance (sub-100 ms query latency). This enables a new class of visualization clients that automatically manage data, quickly and transparently requesting data from the underlying database without requiring the user to explicitly initiate queries. It lays a groundwork for supporting truly interactive exploration of big data and opens new directions for research on scalable information visualization systems. △ Less

Submitted 23 June, 2021; originally announced June 2021.

Comments: 14 pages; Originally published in May 2018 as a technical report in the UC Berkeley EECS Technical Report Series (see https://www2.eecs.berkeley.edu/Pubs/TechRpts/2018/EECS-2018-85.html)

Report number: Technical Report No. UCB/EECS-2018-85

Buckling and yield strength estimation of architected materials under arbitrary loads

Authors: Morten N. Andersen, Yiqiang Wang, Fengwen Wang, Ole Sigmund

Abstract: Buckling strength estimation of architected materials has mainly been restricted to load cases oriented along symmetry axes. However, realistic load scenarios normally exhibit more general stress distributions. In this paper we propose a simple yet accurate method to estimate the buckling strength of stretch-dominated lattice structures based on individual member analysis. As an integral part of t… ▽ More Buckling strength estimation of architected materials has mainly been restricted to load cases oriented along symmetry axes. However, realistic load scenarios normally exhibit more general stress distributions. In this paper we propose a simple yet accurate method to estimate the buckling strength of stretch-dominated lattice structures based on individual member analysis. As an integral part of the method, the yield strength is also determined. This simplified model is verified by rigorous numerical analysis. In particular, we efficiently compute the complete buckling strength surfaces of an orthotropic bulk modulus optimal plate lattice structure and isotropic stiffness optimal plate and truss lattice structures subjected to rotated uni-axial loads, where the ratio between the highest and lowest buckling strength is found to be 1.77, 2.11 and 2.41, respectively. For comparison, we also provide their yield strength surfaces, where the corresponding ratios are 1.84, 1.16 and 1.79. Furthermore, we use the knowledge gained from the simplified model to create a new configuration of the isotropic plate lattice structure with a more isotropic buckling strength surface and buckling strength ratio of 1.24, without deterioration of the stiffness or yield strength. The proposed method provides a valuable tool to quickly estimate the microstructural buckling strength of stretch-dominated lattice structures, especially for applications where the stress state is non-uniform such as infill in additive manufacturing. △ Less

Submitted 13 May, 2021; originally announced May 2021.

Comments: 19 pages with 5 figures

TOI-1431b/MASCARA-5b: A Highly Irradiated Ultra-Hot Jupiter Orbiting One of the Hottest & Brightest Known Exoplanet Host Stars

Authors: Brett Christopher Addison, Emil Knudstrup, Ian Wong, Guillaume Hebrard, Patrick Dorval, Ignas Snellen, Simon Albrecht, Aaron Bello-Arufe, Jose-Manuel Almenara, Isabelle Boisse, Xavier Bonfils, Shweta Dalal, Olivier Demangeon, Sergio Hoyer, Flavien Kiefer, N. C. Santos, Grzegorz Nowak, Rafael Luque, Monika Stangret, Enric Palle, Rene Tronsgaard, Victoria Antoci, Lars A. Buchhave, Maximilian N. Gunther, Tansu Daylan , et al. (48 additional authors not shown)

Abstract: We present the discovery of a highly irradiated and moderately inflated ultra-hot Jupiter, TOI-1431b/MASCARA-5b (HD 201033b), first detected by NASA's Transiting Exoplanet Survey Satellite mission (TESS) and the Multi-site All-Sky CAmeRA (MASCARA). The signal was established to be of planetary origin through radial velocity measurements obtained using SONG, SOPHIE, FIES, NRES, and EXPRES, which sh… ▽ More We present the discovery of a highly irradiated and moderately inflated ultra-hot Jupiter, TOI-1431b/MASCARA-5b (HD 201033b), first detected by NASA's Transiting Exoplanet Survey Satellite mission (TESS) and the Multi-site All-Sky CAmeRA (MASCARA). The signal was established to be of planetary origin through radial velocity measurements obtained using SONG, SOPHIE, FIES, NRES, and EXPRES, which show a reflex motion of $K=294.1\pm1.1$ m s$^{-1}$. A joint analysis of the TESS and ground-based photometry and radial velocity measurements reveals that TOI-1431b has a mass of $M_{p}=3.12\pm0.18$ $\rm{M_J}$ ($990\pm60$ M$_{\oplus}$), an inflated radius of $R_{p}=1.49\pm0.05$ $\rm{R_J}$ ($16.7\pm0.6$ R$_{\oplus}$), and an orbital period of $P=2.650237\pm0.000003$ d. Analysis of the spectral energy distribution of the host star reveals that the planet orbits a bright ($\mathrm{V}=8.049$ mag) and young ($0.29^{+0.32}_{-0.19}$ Gyr) Am type star with $T_{\rm eff}=7690^{+400}_{-250}$ $\rm{K}$, resulting in a highly irradiated planet with an incident flux of $\langle F \rangle=7.24^{+0.68}_{-0.64}\times$10$^9$ erg s$^{-1}$ cm$^{-2}$ ($5300^{+500}_{-470}\mathrm{S_{\oplus}}$) and an equilibrium temperature of $T_{eq}=2370\pm70$ K. TESS photometry also reveals a secondary eclipse with a depth of $127^{+4}_{-5}$ppm as well as the full phase curve of the planet's thermal emission in the red-optical. This has allowed us to measure the dayside and nightside temperature of its atmosphere as $T_\mathrm{day}=3004\pm64$ K and $T_\mathrm{night}=2583\pm63$ K, the second hottest measured nightside temperature. The planet's low day/night temperature contrast ($\sim$420 K) suggests very efficient heat transport between the dayside and nightside hemispheres. △ Less

Submitted 23 September, 2021; v1 submitted 25 April, 2021; originally announced April 2021.

Comments: Accepted for publication in the Astronomical Journal. 39 pages, 18 figures, and 4 tables

Computer Algebra in R with caracas

Authors: Mikkel Meyer Andersen, Søren Højsgaard

Abstract: The capability of R to do symbolic mathematics is enhanced by the caracas package. This package uses the Python computer algebra library SymPy as a back-end but caracas is tightly integrated in the R environment, thereby enabling the R user with symbolic mathematics within R. Key components of the caracas package are illustrated in this paper. Examples are taken from statistics and mathematics. Th… ▽ More The capability of R to do symbolic mathematics is enhanced by the caracas package. This package uses the Python computer algebra library SymPy as a back-end but caracas is tightly integrated in the R environment, thereby enabling the R user with symbolic mathematics within R. Key components of the caracas package are illustrated in this paper. Examples are taken from statistics and mathematics. The caracas package integrates well with e.g. Rmarkdown, and as such creation of scientific reports and teaching is supported. △ Less

Submitted 19 April, 2021; v1 submitted 12 April, 2021; originally announced April 2021.

Regularization by Denoising Sub-sampled Newton Method for Spectral CT Multi-Material Decomposition

Authors: Alessandro Perelli, Martin S. Andersen

Abstract: Spectral Computed Tomography (CT) is an emerging technology that enables to estimate the concentration of basis materials within a scanned object by exploiting different photon energy spectra. In this work, we aim at efficiently solving a model-based maximum-a-posterior problem to reconstruct multi-materials images with application to spectral CT. In particular, we propose to solve a regularized o… ▽ More Spectral Computed Tomography (CT) is an emerging technology that enables to estimate the concentration of basis materials within a scanned object by exploiting different photon energy spectra. In this work, we aim at efficiently solving a model-based maximum-a-posterior problem to reconstruct multi-materials images with application to spectral CT. In particular, we propose to solve a regularized optimization problem based on a plug-in image-denoising function using a randomized second order method. By approximating the Newton step using a sketching of the Hessian of the likelihood function, it is possible to reduce the complexity while retaining the complex prior structure given by the data-driven regularizer. We exploit a non-uniform block sub-sampling of the Hessian with inexact but efficient Conjugate gradient updates that require only Jacobian-vector products for denoising term. Finally, we show numerical and experimental results for spectral CT materials decomposition. △ Less

Submitted 25 March, 2021; originally announced March 2021.

Comments: Accepted in Philosophical Transactions A, issue "Synergistic tomographic image reconstruction (Part 1)"

Interface between graphene and liquid Cu from molecular dynamics simulations

Authors: Juan Santiago Cingolani, Martin Deimel, Simone Köcher, Christoph Scheurer, Karsten Reuter, Mie Andersen

Abstract: Controllable synthesis of defect-free graphene is crucial for applications since the properties of graphene are highly sensitive to any deviations from the crystalline lattice. We focus here on the emerging use of liquid Cu catalysts, which has high potential for fast and efficient industrial-scale production of high-quality graphene. The interface between graphene and liquid Cu is studied using f… ▽ More Controllable synthesis of defect-free graphene is crucial for applications since the properties of graphene are highly sensitive to any deviations from the crystalline lattice. We focus here on the emerging use of liquid Cu catalysts, which has high potential for fast and efficient industrial-scale production of high-quality graphene. The interface between graphene and liquid Cu is studied using force field and ab initio molecular dynamics, revealing a complete or partial embedding of finite-sized flakes. By analyzing flakes of different sizes we find that the size-dependence of the embedding can be rationalized based on the energy cost of embedding versus bending the graphene flake. The embedding itself is driven by the formation of covalent bonds between the under-coordinated edge C atoms and the liquid Cu surface, which is accompanied by a significant charge transfer. In contrast, the central flake atoms are located around or slightly above 3 Å from the liquid Cu surface and exhibit weak vdW-bonding and much lower charge transfer. The structural and electronic properties of the embedded state revealed in our work provides the atomic-scale information needed to develop effective models to explain the special growth observed in experiments where various interesting phenomena such as flake self-assembly and rotational alignment, high growth speeds and low defect densities in the final graphene product have been observed. △ Less

Submitted 5 March, 2021; originally announced March 2021.

Comments: This article may be downloaded for personal use only. Any other use requires prior permission of the author and AIP Publishing. This article appeared in J. Chem. Phys. 153, 074702 (2020) and may be found at https://doi.org/10.1063/5.0020126

Journal ref: J. Chem. Phys. 153, 074702 (2020)

Challenges and Opportunities in High-dimensional Variational Inference

Authors: Akash Kumar Dhaka, Alejandro Catalina, Manushi Welandawe, Michael Riis Andersen, Jonathan Huggins, Aki Vehtari

Abstract: Current black-box variational inference (BBVI) methods require the user to make numerous design choices -- such as the selection of variational objective and approximating family -- yet there is little principled guidance on how to do so. We develop a conceptual framework and set of experimental tools to understand the effects of these choices, which we leverage to propose best practices for maxim… ▽ More Current black-box variational inference (BBVI) methods require the user to make numerous design choices -- such as the selection of variational objective and approximating family -- yet there is little principled guidance on how to do so. We develop a conceptual framework and set of experimental tools to understand the effects of these choices, which we leverage to propose best practices for maximizing posterior approximation accuracy. Our approach is based on studying the pre-asymptotic tail behavior of the density ratios between the joint distribution and the variational approximation, then exploiting insights and tools from the importance sampling literature. Our framework and supporting experiments help to distinguish between the behavior of BBVI methods for approximating low-dimensional versus moderate-to-high-dimensional posteriors. In the latter case, we show that mass-covering variational objectives are difficult to optimize and do not improve accuracy, but flexible variational families can improve accuracy and the effectiveness of importance sampling -- at the cost of additional optimization challenges. Therefore, for moderate-to-high-dimensional posteriors we recommend using the (mode-seeking) exclusive KL divergence since it is the easiest to optimize, and improving the variational family or using model parameter transformations to make the posterior and optimal variational approximation more similar. On the other hand, in low-dimensional settings, we show that heavy-tailed variational families and mass-covering divergences are effective and can increase the chances that the approximation can be improved by importance sampling. △ Less

Submitted 30 June, 2021; v1 submitted 1 March, 2021; originally announced March 2021.