Showing 101–150 of 493 results for author: Andersen, M

Hydrodynamics of noncircular vortices in beams of light and other two-dimensional fluids

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

Abstract: The motion of noncircular two-dimensional vortices is shown to depend on a form of coupling between vortex ellipticity and the gradient of fluid density. The approach is based on the perspective that an elliptic vortex can be described as the projection of a virtual construct, a circular vortex with a symmetry axis that is tilted with respect to the direction of propagation. The resulting kinetic… ▽ More The motion of noncircular two-dimensional vortices is shown to depend on a form of coupling between vortex ellipticity and the gradient of fluid density. The approach is based on the perspective that an elliptic vortex can be described as the projection of a virtual construct, a circular vortex with a symmetry axis that is tilted with respect to the direction of propagation. The resulting kinetic equation offers insights into how tilt and vortex velocity coevolve in few-body nonequilibrium settings such as vortex pair nucleation and annihilation. The model is developed and applied in association with optical vortices, and optical experiments are used to verify its predictive power. It is valid for quantum fluids and classical hydrodynamics settings as well. △ Less

Submitted 20 September, 2021; v1 submitted 16 February, 2021; originally announced February 2021.

Journal ref: Phys. Rev. A 104, 033520 (2021)

The stellar content of H72.97-69.39, a potential super star cluster in the making

Authors: M. Andersen, H. Zinnecker, A. S. Hirschauer, O. Nayak, M. Meixner

Abstract: Young Massive Clusters (YMCs) and Super Star Clusters (SSCs) represent an extreme mode of star formation. Far-infrared imaging of the Magellanic Clouds has identified one potential embedded SSC, HSO BMHERICC J72.971176-69.391112 (HH in short), in the south-west outskirts of the Large Magellanic Cloud. We present Gemini Flamingos 2 and GSAOI near-infrared imaging of a 3'x3' region around HH in orde… ▽ More Young Massive Clusters (YMCs) and Super Star Clusters (SSCs) represent an extreme mode of star formation. Far-infrared imaging of the Magellanic Clouds has identified one potential embedded SSC, HSO BMHERICC J72.971176-69.391112 (HH in short), in the south-west outskirts of the Large Magellanic Cloud. We present Gemini Flamingos 2 and GSAOI near-infrared imaging of a 3'x3' region around HH in order to characterize the stellar content of the cluster. The stellar content is probed down to 1.5 Msun. We find substantial dust extinction across the cluster region, extending up to A_K of 3. Deeply embedded stars are associated with ALMA-detected molecular gas suggesting that star formation is ongoing. The high spatial resolution of the GSAOI data allows identification of the central massive object associated with the 13CO ALMA observations and to detect fainter low-mass stars around the H30alpha ALMA source. The morphology of the molecular gas and the nebulosity from adjacent star formation suggest they have interacted covering a region of several pc. The total stellar content in the cluster is estimated from the intermediate- and high-mass stellar content to be at least 10000 Msun, less than R136 with up to 100 000 Msun within 4.7 pc radius, but places it in the regime of a super star cluster. Based on the extinction determination of individual stars we estimate a molecular gas mass in the vicinity of HH of 6600 Msun, suggesting more star formation can be expected. △ Less

Submitted 12 February, 2021; originally announced February 2021.

Comments: Acceepted, AJ

Precise radial velocities of giant stars XV. Mysterious nearly periodic radial velocity variations in the eccentric binary $ε$ Cygni

Authors: Paul Heeren, Sabine Reffert, Trifon Trifonov, Ka Ho Wong, Man Hoi Lee, Jorge Lillo-Box, Andreas Quirrenbach, Torben Arentoft, Simon Albrecht, Frank Grundahl, Mads Fredslund Andersen, Victoria Antoci, Pere L. Pallé

Abstract: Using the Hamilton Echelle Spectrograph at Lick Observatory, we have obtained precise radial velocities (RVs) of a sample of 373 G- and K-giant stars over more than 12 years, leading to the discovery of several single and multiple planetary systems. The RVs of the long-period (~53 years) spectroscopic binary $ε$ Cyg (HIP 102488) are found to exhibit additional regular variations with a much shorte… ▽ More Using the Hamilton Echelle Spectrograph at Lick Observatory, we have obtained precise radial velocities (RVs) of a sample of 373 G- and K-giant stars over more than 12 years, leading to the discovery of several single and multiple planetary systems. The RVs of the long-period (~53 years) spectroscopic binary $ε$ Cyg (HIP 102488) are found to exhibit additional regular variations with a much shorter period (~291 days). We intend to improve the orbital solution of the $ε$ Cyg system and attempt to identify the cause of the nearly periodic shorter period variations, which might be due to an additional substellar companion. We used precise RV measurements of the K-giant star $ε$ Cyg from Lick Observatory, in combination with a large set of RVs collected more recently with the SONG telescope, as well as archival data sets. Our Keplerian model to the RVs characterizes the orbit of the spectroscopic binary to higher precision than achieved previously, resulting in a semi-major axis of $a = 15.8 \mathrm{AU}$, an eccentricity of $e = 0.93$, and a minimum mass of the secondary of $m \sin i = 0.265 M_\odot$. Additional short-period RV variations closely resemble the signal of a Jupiter-mass planet orbiting the evolved primary component with a period of $291 \mathrm{d}$, but the period and amplitude of the putative orbit change strongly over time. Furthermore, in our stability analysis of the system, no stable orbits could be found in a large region around the best fit. Both of these findings deem a planetary cause of the RV variations unlikely. Most of the investigated alternative scenarios, such as an hierarchical triple or stellar spots, also fail to explain the observed variability convincingly. Due to its very eccentric binary orbit, it seems possible, however, that $ε$ Cyg could be an extreme example of a heartbeat system. △ Less

Submitted 3 February, 2021; originally announced February 2021.

Comments: 17 pages, 13 figures, accepted to A&A

Journal ref: A&A 647, A160 (2021)

A more probable explanation for a continuum flash in the direction of a redshift $\approx$ 11 galaxy

Authors: Charles L. Steinhardt, Michael I. Andersen, Gabriel B. Brammer, Lise Christensen, Johan P. U. Fynbo, Bo Milvang-Jensen, Pascal A. Oesch, Sune Toft

Abstract: Recent work reported the discovery of a gamma-ray burst (GRB) associated with the galaxy GN-z11 at $z\sim 11$. The extreme improbability of the transient source being a GRB in the very early Universe requires robust elimination of all plausible alternative hypotheses. We identify numerous examples of similar transient signals in separate archival MOSFIRE observations and argue that Solar system ob… ▽ More Recent work reported the discovery of a gamma-ray burst (GRB) associated with the galaxy GN-z11 at $z\sim 11$. The extreme improbability of the transient source being a GRB in the very early Universe requires robust elimination of all plausible alternative hypotheses. We identify numerous examples of similar transient signals in separate archival MOSFIRE observations and argue that Solar system objects -- natural or artificial -- are a far more probable explanation for these phenomena. An appendix has been added in response to additional points raised in Jiang et al. (2021), which do not change the conclusion. △ Less

Submitted 22 February, 2021; v1 submitted 29 January, 2021; originally announced January 2021.

Superconducting state of Sr$_2$RuO$_4$ in the presence of longer-range Coulomb interactions

Authors: Astrid T. Rømer, P. J. Hirschfeld, Brian M. Andersen

Abstract: The symmetry of the superconducting condensate in Sr$_2$RuO$_4$ remains controversial after nuclear magnetic resonance (NMR) experiments recently overturned the dominant chiral $p$-wave paradigm. Several theoretical proposals have been put forward to account for existing experiments, including a $d+ig$-wave admixture, conjectured to be stabilized by longer-range Coulomb interactions. We perform a… ▽ More The symmetry of the superconducting condensate in Sr$_2$RuO$_4$ remains controversial after nuclear magnetic resonance (NMR) experiments recently overturned the dominant chiral $p$-wave paradigm. Several theoretical proposals have been put forward to account for existing experiments, including a $d+ig$-wave admixture, conjectured to be stabilized by longer-range Coulomb interactions. We perform a material-specific microscopic theoretical study of pairing by spin- and charge-fluctuations in Sr$_2$RuO$_4$, including the effects of spin-orbit coupling, and both local and longer-range Coulomb repulsion. The latter has important consequences for Sr$_2$RuO$_4$ due to the near-degeneracy of symmetry-distinct pairing states in this material. We find that both the $g$- and $d_{x^2-y^2}$-wave channels remain noncompetitive compared to leading nodal $s'$, $d_{xy}$, and helical ($p$) solutions. This suggests nodal time-reversal symmetry broken $s'+id_{xy}$ or $s'+ip$ phases, promoted by longer-range Coulomb repulsion, as the most favorable candidates for Sr$_2$RuO$_4$. We analyse the properties of these states, and show that the $s'+id_{xy}$ solution agrees with the bulk of available experimental data, including recent discoveries from NMR, muon spin relaxation ($μ$SR), and ultrasound measurements. △ Less

Submitted 18 January, 2021; originally announced January 2021.

Comments: 4 pages, 2 figures + supplementary material

Report number: CMT NBI 2021

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

Observing the changing surface structures of the active K giant sigma Gem with SONG

Authors: H. Korhonen, R. M. Roettenbacher, S. Gu, F. Grundahl, M. F. Andersen, G. W. Henry, J. Jessen-Hansen, V. Antoci, P. L. Pallé

Abstract: Aims: We aim to study the spot evolution and differential rotation in the magnetically active cool K-type giant star sigma Gem from broadband photometry and continuous spectroscopic observations that span 150 nights. Methods: We use high-resolution, high signal-to-noise ratio spectra obtained with the Hertzsprung SONG telescope to reconstruct surface (photospheric) temperature maps with Doppler im… ▽ More Aims: We aim to study the spot evolution and differential rotation in the magnetically active cool K-type giant star sigma Gem from broadband photometry and continuous spectroscopic observations that span 150 nights. Methods: We use high-resolution, high signal-to-noise ratio spectra obtained with the Hertzsprung SONG telescope to reconstruct surface (photospheric) temperature maps with Doppler imaging techniques. The 303 observations span 150 nights and allow for a detailed analysis of the spot evolution and surface differential rotation. The Doppler imaging results are compared to simultaneous broadband photometry from the Tennessee State University T3 0.4 m Automated Photometric Telescope. The activity from the stellar chromosphere, which is higher in the stellar atmosphere, is also studied using SONG observations of Balmer H alpha line profiles and correlated with the photospheric activity. Results: The temperature maps obtained during eight consecutive stellar rotations show mainly high-latitude or polar spots, with the main spot concentrations above latitude 45 deg. The spots concentrate around phase 0.25 near the beginning of our observations and around phase 0.75 towards the end. The photometric observations confirm a small jump in spot phases that occurred in February 2016. The cross-correlation of the temperature maps reveals rather strong solar-like differential rotation, giving a relative surface differential rotation coefficient of $α$ = 0.10 +/- 0.02. There is a weak correlation between the locations of starspots and enhanced emission in the chromosphere at some epochs. △ Less

Submitted 30 December, 2020; originally announced December 2020.

Comments: 19 pages, 10 figures, accepted for publication in Astronomy & Astrophysics, in the published version the Tables B1 and B2 will only be available via CDS

Journal ref: A&A 646, A6 (2021)

Topological Superconductivity Induced by Magnetic Texture Crystals

Authors: Daniel Steffensen, Morten H. Christensen, Brian M. Andersen, Panagiotis Kotetes

Abstract: We present a detailed investigation of the topological phases and Majorana fermion (MF) excitations that arise from the bulk interplay between (un)conventional one/two-band spin-singlet superconductivity and a number of magnetic texture crystals. The latter define inhomogeneous magnetization profiles which consist of a periodically-repeating primitive cell. Here we focus on magnetic texture crysta… ▽ More We present a detailed investigation of the topological phases and Majorana fermion (MF) excitations that arise from the bulk interplay between (un)conventional one/two-band spin-singlet superconductivity and a number of magnetic texture crystals. The latter define inhomogeneous magnetization profiles which consist of a periodically-repeating primitive cell. Here we focus on magnetic texture crystals with a primitive cell of the helix, whirl, and skyrmion types, which feature distinct symmetry properties. We identify a multitude of accessible topological phases which harbor flat, uni- or bi-directional, (quasi-)helical, or chiral MF edge modes. This rich variety originates from the interplay between topological phases with gapped and nodal bulk energy spectra. The types of the emerging topological superconducting phases and the features of the arising MFs are solely determined by the properties and compatibility of the so-called magnetic and pairing point/space groups. Our analysis is general and does not rely on specific parameters of the models employed here to exemplify the topological scenarios which become accessible. Therefore, our results can be extended to systems with multiple bands, are relevant for a wide range of layered materials and hybrid devices, and provide predictions for strong, weak and crystalline topological phases. △ Less

Submitted 19 March, 2022; v1 submitted 17 December, 2020; originally announced December 2020.

Comments: 46 pages, 25 panels of figures, 5 tables, 4 appendices; Version close to the one appearing soon in Physical Review Research

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

Inter-orbital nematicity and the origin of a single electron Fermi pocket in FeSe

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

Abstract: The electronic structure of the enigmatic iron-based superconductor FeSe has puzzled researchers since spectroscopic probes failed to observe the expected electron pocket at the $Y$ point in the 1-Fe Brillouin zone. It has been speculated that this pocket, essential for an understanding of the superconducting state, is either absent or incoherent. Here, we perform a theoretical study of the prefer… ▽ More The electronic structure of the enigmatic iron-based superconductor FeSe has puzzled researchers since spectroscopic probes failed to observe the expected electron pocket at the $Y$ point in the 1-Fe Brillouin zone. It has been speculated that this pocket, essential for an understanding of the superconducting state, is either absent or incoherent. Here, we perform a theoretical study of the preferred nematic order originating from nearest-neighbor Coulomb interactions in an electronic model relevant for FeSe. We find that at low temperatures the dominating nematic components are of inter-orbital $d_{xz}-d_{xy}$ and $d_{yz}-d_{xy}$ character, with spontaneously broken amplitudes for these two components. This inter-orbital nematic order naturally leads to distinct hybridization gaps at the $X$ and $Y$ points of the 1-Fe Brillouin zone, and may thereby produce highly anisotropic Fermi surfaces with only a single electron pocket at one of these momentum-space locations. The associated superconducting gap structure obtained with the generated low-energy electronic band structure from spin-fluctuation mediated pairing agrees well with that measured experimentally. Finally, from a comparison of the computed spin susceptibility to available neutron scattering data, we discuss the necessity of additional self-energy effects, and explore the role of orbital-dependent quasiparticle weights as a minimal means to include them. △ Less

Submitted 13 December, 2020; originally announced December 2020.

Comments: 8 pages, 5 figures + Appedices

Report number: CMT NBI 2020

Journal ref: Phys. Rev. B 103, 054505 (2021)

The Real-Time Correlation Function of Floquet Conformal Fields

Authors: Malthe Andersen, Frederik Nørfjand, Nikolaj Thomas Zinner

Abstract: Conformal field theories (CFT) play an important role in quantum field theories as they allow for exact studies that are hard to access without the conformal symmetries in place. A remarkable recent development is the application of conformal field theory to periodically driven systems, the so-called Floquet-driven CFTs. A number of recent works have found that there are classes of periodically dr… ▽ More Conformal field theories (CFT) play an important role in quantum field theories as they allow for exact studies that are hard to access without the conformal symmetries in place. A remarkable recent development is the application of conformal field theory to periodically driven systems, the so-called Floquet-driven CFTs. A number of recent works have found that there are classes of periodically driven (1+1)-dimensional CFTs that can be analytically solved, and that show both heating and non-heating phases, as well as unusual properties in terms of information scrambling. Hitherto, most results have focused on the dynamics of such models at stroboscopic times, i.e. samples taken at each period of the time evolution. In this paper, we consider time-evolution at arbitrary real times and derive formulas for the exact time-propagation of heating and non-heating classes of Floquet-driven CFTs. This requires a subtle treatment of non-stroboscopic time that leads to different interpretations in terms of curved space-time dynamics in these systems compared to previous results in the literature. △ Less

Submitted 17 November, 2020; originally announced November 2020.

Comments: 10 pages, 10 figures

Journal ref: Phys. Rev. D 103, 056005 (2021)

Role of anion in the pairing interaction of iron-based superconductivity

Authors: J. -X. Yin, Y. Y. Zhao, Zheng Wu, X. X. Wu, A. Kreisel, B. M. Andersen, Gennevieve Macam, Sen Zhou, Rui Wu, Limin Liu, Hanbin Deng, Changjiang Zhu, Yuan Li, Yingkai Sun, Zhi-Quan Huang, Feng-Chuan Chuang, Hsin Lin, C. -S. Ting, J. -P. Hu, Z. Q. Wang, P. C. Dai, H. Ding, S. H. Pan

Abstract: High-temperature iron-based superconductivity develops in a structure with unusual lattice-orbital geometry, based on a planar layer of Fe atoms with 3d orbitals and tetrahedrally coordinated by anions. Here we elucidate the electronic role of anions in the iron-based superconductors utilizing state-of-the-art scanning tunneling microscopy. By measuring the local electronic structure, we find that… ▽ More High-temperature iron-based superconductivity develops in a structure with unusual lattice-orbital geometry, based on a planar layer of Fe atoms with 3d orbitals and tetrahedrally coordinated by anions. Here we elucidate the electronic role of anions in the iron-based superconductors utilizing state-of-the-art scanning tunneling microscopy. By measuring the local electronic structure, we find that As anion in Ba0.4K0.6Fe2As2 has a striking impact on the electron pairing. The superconducting electronic feature can be switched off/on by removing/restoring As atoms on Fe layer at the atomic scale. Our analysis shows that this remarkable atomic switch effect is related to the geometrical cooperation between anion mediated hopping and unconventional pairing interaction. Our results uncover that the local Fe-anion coupling is fundamental for the pairing interaction of iron-based superconductivity, and promise the potential of bottom-up engineering of electron pairing. △ Less

Submitted 15 November, 2020; originally announced November 2020.

Orbital-dependent self-energy effects and consequences for the superconducting gap structure in multi-orbital correlated electron systems

Authors: Kristofer Björnson, Andreas Kreisel, Astrid T. Rømer, Brian M. Andersen

Abstract: We perform a theoretical study of the effects of electronic correlations on the superconducting gap structure of multi-band superconductors. In particular, by comparing standard RPA-based spin-fluctuation mediated gap structures to those obtained within the FLEX formalism for an iron-based superconductor, we obtain directly the feedback effects from electron-electron interactions on the momentum-s… ▽ More We perform a theoretical study of the effects of electronic correlations on the superconducting gap structure of multi-band superconductors. In particular, by comparing standard RPA-based spin-fluctuation mediated gap structures to those obtained within the FLEX formalism for an iron-based superconductor, we obtain directly the feedback effects from electron-electron interactions on the momentum-space gap structure. We show how self-energy effects can lead to an orbital inversion of the orbital-resolved spin susceptibility, and thereby invert the hierarchy of the most important orbitals channels for superconducting pairing. This effect has important consequences for the detailed gap variations on the Fermi surface. We expect such self-energy feedback on the pairing gap to be generally relevant for superconductivity in strongly correlated multi-orbital systems. △ Less

Submitted 13 January, 2021; v1 submitted 9 November, 2020; originally announced November 2020.

Comments: 8 pages, 5 figures

Report number: NBI CMT 2020

Journal ref: Phys. Rev. B 103, 024508 (2021)

On the competition for ultimately stiff and strong architected materials

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

Abstract: Advances in manufacturing techniques may now realize virtually any imaginable microstructures, paving the way for architected materials with properties beyond those found in nature. This has lead to a quest for closing gaps in property-space by carefully designed metamaterials. Development of mechanical metamaterials has gone from open truss lattice structures to closed plate lattice structures wi… ▽ More Advances in manufacturing techniques may now realize virtually any imaginable microstructures, paving the way for architected materials with properties beyond those found in nature. This has lead to a quest for closing gaps in property-space by carefully designed metamaterials. Development of mechanical metamaterials has gone from open truss lattice structures to closed plate lattice structures with stiffness close to theoretical bounds. However, the quest for optimally stiff and strong materials is complex. Plate lattice structures have higher stiffness and (yield) strength but are prone to buckling at low volume fractions. Hence here, truss lattice structures may still be optimal. To make things more complicated, hollow trusses or structural hierarchy bring closed-walled microstructures back in the competition. Based on analytical and numerical studies of common microstructures from the literature, we provide higher order interpolation schemes for their effective stiffness and (buckling) strength. Furthermore, we provide a case study based on multi-property Ashby charts for weight-optimal porous beams under bending, that demonstrates the intricate interplay between structure and microarchitecture that plays the key role in the design of ultimate load carrying structures. The provided interpolation schemes may also be used to account for microstructural yield and buckling in multiscale design optimization schemes. △ Less

Submitted 8 November, 2020; originally announced November 2020.

Comments: 18 pages main manuscript with 6 figures and 10 pages appendix with 2 figures

MSC Class: 74Q15

Properties of the Hyades, the eclipsing binary HD27130, and the oscillating red giant $ε$ Tau

Authors: K. Brogaard, E. Pakštienė, F. Grundahl, Š. Mikolaitis, G. Tautvaišienė, D. Slumstrup, G. J. J. Talens, D. A. VandenBerg, A. Miglio, T. Arentoft, H. Kjeldsen, R. Janulis, A. Drazdauskas, A. Marchini, R. Minkevičiūtė, E. Stonkutė, V. Bagdonas, M. Fredslund Andersen, J. Jessen-Hansen, P. L. Pallé, P. Dorval, I. A. G. Snellen, G. P. P. L. Otten, T. R. White

Abstract: Eclipsing binary stars allow derivation of accurate and precise masses and radii. When they reside in star clusters, properties of even higher precision, along with additional information, can be extracted. Asteroseismology of solar-like oscillations offers similar possibilities for single stars. We improve the previously established properties of the Hyades eclipsing binary HD27130 and re-assess… ▽ More Eclipsing binary stars allow derivation of accurate and precise masses and radii. When they reside in star clusters, properties of even higher precision, along with additional information, can be extracted. Asteroseismology of solar-like oscillations offers similar possibilities for single stars. We improve the previously established properties of the Hyades eclipsing binary HD27130 and re-assess the asteroseismic properties of the giant star $ε$ Tau. The physical properties of these members of the Hyades are then used to constrain the helium content and age of the cluster. New multi-colour light curves were combined with multi-epoch radial velocities to yield masses and radii of HD27130. $T_{\rm eff}$ was derived from spectroscopy and photometry, and verified using the Gaia parallax. We estimate the cluster age from re-evaluated asteroseismic properties of $ε$ Tau while using HD27130 to constrain the helium content. The masses and radii, and $T_{\rm eff}$ of HD 27130 were found to be $M=1.0245\pm0.0024 M_{\odot}$, $R=0.9226\pm0.015 R_{\odot}$, $T_{\rm eff}=5650\pm50$ K for the primary, and $M=0.7426\pm0.0016 M_{\odot}$, $R=0.7388\pm0.026 R_{\odot}$, $T_{\rm eff}=4300\pm100$ K for the secondary component. Our re-evaluation of $ε$ Tau suggests that the previous literature estimates are trustworthy, and that the Hipparcos parallax is more reliable than the Gaia DR2 parallax. The helium content of HD27130 and thus of the Hyades is found to be $Y=0.27$ but with significant model dependence. Correlations with the adopted metallicity results in a robust helium enrichment law with $\frac{ΔY}{ΔZ}$ close to 1.2. We estimate the age of the Hyades to be 0.9 $\pm$ 0.1 (stat) $\pm$ 0.1 (sys) Gyr in slight tension with recent age estimates based on the cluster white dwarfs. (abridged) △ Less

Submitted 4 November, 2020; originally announced November 2020.

Comments: 14 pages, 8 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 645, A25 (2021)

Singular magnetic anisotropy in the nematic phase of FeSe

Authors: R. Zhou, D. D. Scherer, H. Mayaffre, P. Toulemonde, M. Ma, Y. Li, B. M. Andersen, M. -H. Julien

Abstract: FeSe is arguably the simplest, yet the most enigmatic, iron-based superconductor. Its nematic but non-magnetic ground state is unprecedented in this class of materials and stands out as a current puzzle. Here, our NMR measurements in the nematic state of mechanically detwinned FeSe reveal that both the Knight shift and the spin-lattice relaxation rate 1/T_1 possess an in-plane anisotropy opposite… ▽ More FeSe is arguably the simplest, yet the most enigmatic, iron-based superconductor. Its nematic but non-magnetic ground state is unprecedented in this class of materials and stands out as a current puzzle. Here, our NMR measurements in the nematic state of mechanically detwinned FeSe reveal that both the Knight shift and the spin-lattice relaxation rate 1/T_1 possess an in-plane anisotropy opposite to that of the iron pnictides LaFeAsO and BaFe2As2. Using a microscopic electron model that includes spin-orbit coupling, our calculations show that an opposite quasiparticle weight ratio between the d_xz and d_yz orbitals leads to an opposite anisotropy of the orbital magnetic susceptibility, which explains our Knight shift results. We attribute this property to a different nature of nematic order in the two compounds, predominantly bond-type in FeSe and onsite ferro-orbital in pnictides. The T_1 anisotropy is found to be inconsistent with existing neutron scattering data in FeSe, showing that the spin fluctuation spectrum reveals surprises at low energy, possibly from fluctuations that do not break C_4 symmetry. Therefore, our results reveal that important information is hidden in these anisotropies and they place stringent constraints on the low-energy spin correlations as well as on the nature of nematicity in FeSe. △ Less

Submitted 2 October, 2020; originally announced October 2020.

Comments: Submitted version. More recent version and Supplementary Information (24 pages) available upon request

Journal ref: npj Quantum Mater. 5, 93 (2020)

Robust, Accurate Stochastic Optimization for Variational Inference

Authors: Akash Kumar Dhaka, Alejandro Catalina, Michael Riis Andersen, Måns Magnusson, Jonathan H. Huggins, Aki Vehtari

Abstract: We consider the problem of fitting variational posterior approximations using stochastic optimization methods. The performance of these approximations depends on (1) how well the variational family matches the true posterior distribution,(2) the choice of divergence, and (3) the optimization of the variational objective. We show that even in the best-case scenario when the exact posterior belongs… ▽ More We consider the problem of fitting variational posterior approximations using stochastic optimization methods. The performance of these approximations depends on (1) how well the variational family matches the true posterior distribution,(2) the choice of divergence, and (3) the optimization of the variational objective. We show that even in the best-case scenario when the exact posterior belongs to the assumed variational family, common stochastic optimization methods lead to poor variational approximations if the problem dimension is moderately large. We also demonstrate that these methods are not robust across diverse model types. Motivated by these findings, we develop a more robust and accurate stochastic optimization framework by viewing the underlying optimization algorithm as producing a Markov chain. Our approach is theoretically motivated and includes a diagnostic for convergence and a novel stopping rule, both of which are robust to noisy evaluations of the objective function. We show empirically that the proposed framework works well on a diverse set of models: it can automatically detect stochastic optimization failure or inaccurate variational approximation △ Less

Submitted 3 September, 2020; v1 submitted 1 September, 2020; originally announced September 2020.

Journal ref: NeurIPS 2020

Pair Correlations and Photoassociation Dynamics of Two Atoms in an Optical Tweezer

Authors: M. Weyland, S. S. Szigeti, R. A. B. Hobbs, P. Ruksasakchai, L. Sanchez, M. F. Andersen

Abstract: We investigate the photoassociation dynamics of exactly two laser-cooled $^{85}$Rb atoms in an optical tweezer and reveal fundamentally different behavior to photoassociation in many-atom ensembles. We observe non-exponential decay in our two-atom experiment that cannot be described by a single rate coefficient and find its origin in our system's pair correlation. This is in stark contrast to many… ▽ More We investigate the photoassociation dynamics of exactly two laser-cooled $^{85}$Rb atoms in an optical tweezer and reveal fundamentally different behavior to photoassociation in many-atom ensembles. We observe non-exponential decay in our two-atom experiment that cannot be described by a single rate coefficient and find its origin in our system's pair correlation. This is in stark contrast to many-atom photoassociation dynamics, which are governed by exponential decay with a single rate coefficient. We also investigate photoassociation in a three-atom system, thereby probing the transition from two-atom dynamics to many-atom dynamics. Our experiments reveal additional reaction dynamics that are only accessible through the control of single atoms and suggest photoassociation could measure pair correlations in few-atom systems. It further showcases our complete control over the quantum state of individual atoms and molecules, which provides information unobtainable from many-atom experiments. △ Less

Submitted 25 January, 2021; v1 submitted 31 August, 2020; originally announced September 2020.

Comments: main text: 5 pages, 4 figures supplemental: 5 pages, 7 figures

Journal ref: Phys. Rev. Lett. 126, 083401 (2021)

Trapping Majorana Zero Modes in Vortices of Magnetic Texture Crystals Coupled to Nodal Superconductors

Authors: Daniel Steffensen, Brian M. Andersen, Panagiotis Kotetes

Abstract: We propose a mechanism for opening a full bulk energy gap and inducing vortex Majorana zero modes (MZMs) in nodal superconductors (SCs). We show that this becomes possible by coupling the nodal SC of interest to a magnetic texture crystal. The latter consists of superpositions of magnetic textures which repeat periodically in space according to suitable wave vectors that enable spin-flip scatterin… ▽ More We propose a mechanism for opening a full bulk energy gap and inducing vortex Majorana zero modes (MZMs) in nodal superconductors (SCs). We show that this becomes possible by coupling the nodal SC of interest to a magnetic texture crystal. The latter consists of superpositions of magnetic textures which repeat periodically in space according to suitable wave vectors that enable spin-flip scattering between all pairs of nodes of the SC, and thus open a full gap in its bulk energy spectrum. In this event, MZMs can be trapped in spin or shift vortices introduced in the magnetic texture crystal. Our approach is generic and applies to nodal SCs of spin-singlet, -triplet, or -mixed type of pairing. Therefore, it promises to find application in a variety of nodal SCs, where the magnetic textures appear either spontaneously due to electron-electron interactions, or are imposed by nanomagnets, or become induced by coupling to a lattice of localized magnetic moments. △ Less

Submitted 24 October, 2021; v1 submitted 24 August, 2020; originally announced August 2020.

Comments: 17 pages, 5 Figures. Modified title and text. Extended version with new results on the experimental realization. Applications regarding hybrid devices will be discussed elsewhere

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

The PDR structure and kinematics around the compact HII regions S235A and S235C with [CII], [13CII], [OI] and HCO+ line profiles

Authors: M. S. Kirsanova, V. Ossenkopf-Okada, L. D. Anderson, P. A. Boley, J. H. Bieging, Ya. N. Pavlyuchenkov, M. Luisi, N. Schneider, M. Andersen, M. R. Samal, A. M. Sobolev, C. Buchbender, R. Aladro, Y. Okada

Abstract: The aim of the present work is to study structure and gas kinematics in the photodissociation regions (PDRs) around the compact HII regions S235A and S235C. We observe the [CII], [13CII] and [OI] line emission, using SOFIA/upGREAT and complement them by data of HCO+ and CO. We use the [13CII] line to measure the optical depth of the [CII] emission, and find that the [CII] line profiles are influen… ▽ More The aim of the present work is to study structure and gas kinematics in the photodissociation regions (PDRs) around the compact HII regions S235A and S235C. We observe the [CII], [13CII] and [OI] line emission, using SOFIA/upGREAT and complement them by data of HCO+ and CO. We use the [13CII] line to measure the optical depth of the [CII] emission, and find that the [CII] line profiles are influenced by self-absorption, while the [13CII] line remains unaffected by these effects. Hence, for dense PDRs, [13CII] emission is a better tracer of gas kinematics. The optical depth of the [CII] line is up to 10 in S235A. We find an expanding motion of the [CII]-emitting layer of the PDRs into the front molecular layer in both regions. Comparison of the gas and dust columns shows that gas components visible neither in the [CII] nor in low-J CO lines may contribute to the total column across S235A. We test whether the observed properties of the PDRs match the predictions of spherical models of expanding HII region + PDR + molecular cloud. Integrated intensities of the [13CII], [CII] and [OI] lines are well-represented by the model, but the models do not reproduce the double-peaked [CII] line profiles due to an insufficient column density of C+. The model predicts that the [OI] line could be a more reliable tracer of gas kinematics, but the foreground self-absorbing material does not allow using it in the considered regions. △ Less

Submitted 30 July, 2020; originally announced July 2020.

Comments: accepted to MNRAS

Spin-orbit quantum impurity and quantization in a topological magnet

Authors: Jia-Xin Yin, Nana Shumiya, Yuxiao Jiang, Huibin Zhou, Gennevieve Macam, Hano Omar Mohammad Sura, Songtian S. Zhang, Zijia Cheng, Zurab Guguchia, Yangmu Li, Qi Wang, Maksim Litskevich, Ilya Belopolski, Xian Yang, Tyler A. Cochran, Guoqing Chang, Qi Zhang, Zhi-Quan Huang, Feng-Chuan Chuang, Hsin Lin, Hechang Lei, Brian M. Andersen, Ziqiang Wang, Shuang Jia, M. Zahid Hasan

Abstract: Quantum states induced by single-atomic impurities are at the frontier of physics and material science. While such states have been reported in high-temperature superconductors and dilute magnetic semiconductors, they are unexplored in topological magnets which can feature spin-orbit tunability. Here we use spin-polarized scanning tunneling microscopy/spectroscopy (STM/S) to study the engineered q… ▽ More Quantum states induced by single-atomic impurities are at the frontier of physics and material science. While such states have been reported in high-temperature superconductors and dilute magnetic semiconductors, they are unexplored in topological magnets which can feature spin-orbit tunability. Here we use spin-polarized scanning tunneling microscopy/spectroscopy (STM/S) to study the engineered quantum impurity in a topological magnet Co3Sn2S2. We find that each substituted In impurity introduces a striking localized bound state. Our systematic magnetization-polarized probe reveals that this bound state is spin-down polarized, in lock with a negative orbital magnetization. Moreover, the magnetic bound states of neighboring impurities interact to form quantized orbitals, exhibiting an intriguing spin-orbit splitting, analogous to the splitting of the topological fermion line. Our work collectively demonstrates the strong spin-orbit effect of the single-atomic impurity at the quantum level, suggesting that a nonmagnetic impurity can introduce spin-orbit coupled magnetic resonance in topological magnets. △ Less

Submitted 26 July, 2020; originally announced July 2020.

Comments: To appear in Nature Communications (2020), both experiment and theory included. arXiv admin note: text overlap with arXiv:2002.11783

Journal ref: Nature Communications 11, 4415 (2020)

State Space Expectation Propagation: Efficient Inference Schemes for Temporal Gaussian Processes

Authors: William J. Wilkinson, Paul E. Chang, Michael Riis Andersen, Arno Solin

Abstract: We formulate approximate Bayesian inference in non-conjugate temporal and spatio-temporal Gaussian process models as a simple parameter update rule applied during Kalman smoothing. This viewpoint encompasses most inference schemes, including expectation propagation (EP), the classical (Extended, Unscented, etc.) Kalman smoothers, and variational inference. We provide a unifying perspective on thes… ▽ More We formulate approximate Bayesian inference in non-conjugate temporal and spatio-temporal Gaussian process models as a simple parameter update rule applied during Kalman smoothing. This viewpoint encompasses most inference schemes, including expectation propagation (EP), the classical (Extended, Unscented, etc.) Kalman smoothers, and variational inference. We provide a unifying perspective on these algorithms, showing how replacing the power EP moment matching step with linearisation recovers the classical smoothers. EP provides some benefits over the traditional methods via introduction of the so-called cavity distribution, and we combine these benefits with the computational efficiency of linearisation, providing extensive empirical analysis demonstrating the efficacy of various algorithms under this unifying framework. We provide a fast implementation of all methods in JAX. △ Less

Submitted 12 July, 2020; originally announced July 2020.

Comments: Accepted to International Conference on Machine Learning (ICML) 2020

On the Remarkable Superconductivity of FeSe and its Close Cousins

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

Abstract: Emergent electronic phenomena in iron-based superconductors have been at the forefront of condensed matter physics for more than a decade. Much has been learned about the origin and intertwined roles of ordered phases, including nematicity, magnetism, and superconductivity, in this fascinating class of materials. In recent years, focus has been centered on the peculiar and highly unusual propertie… ▽ More Emergent electronic phenomena in iron-based superconductors have been at the forefront of condensed matter physics for more than a decade. Much has been learned about the origin and intertwined roles of ordered phases, including nematicity, magnetism, and superconductivity, in this fascinating class of materials. In recent years, focus has been centered on the peculiar and highly unusual properties of FeSe and its close cousins. This family of materials has attracted considerable attention due to the discovery of unexpected superconducting gap structures, a wide range of superconducting critical temperatures, and evidence for nontrivial band topology, including associated spin-helical surface states and vortex-induced Majorana bound states. Here, we review superconductivity in iron chalcogenide superconductors, including bulk FeSe, doped bulk FeSe, FeTe$_{1-x}$Se$_x$, intercalated FeSe materials, and monolayer FeSe and FeTe$_{1-x}$Se$_x$ on SrTiO$_3$. We focus on the superconducting properties, including a survey of the relevant experimental studies, and a discussion of the different proposed theoretical pairing scenarios. In the last part of the paper, we review the growing recent evidence for nontrivial topological effects in FeSe-related materials, focusing again on interesting implications for superconductivity. △ Less

Submitted 10 September, 2020; v1 submitted 6 July, 2020; originally announced July 2020.

Comments: review article, 57 pages, 23 figures

Report number: NBI CMT 2020

Journal ref: Symmetry 12, 1402 (2020)

Spatially dispersing Yu-Shiba-Rusinov states in the unconventional superconductor $\mathrm{FeTe}_{0.55}\mathrm{Se}_{0.45}$

Authors: Damianos Chatzopoulos, Doohee Cho, Koen M. Bastiaans, Gorm O. Steffensen, Damian Bouwmeester, Alireza Akbari, Genda Gu, Jens Paaske, Brian M. Andersen, Milan P. Allan

Abstract: By using scanning tunneling microscopy (STM) we find and characterize dispersive, energy-symmetric in-gap states in the iron-based superconductor $\mathrm{FeTe}_{0.55}\mathrm{Se}_{0.45}$, a material that exhibits signatures of topological superconductivity, and Majorana bound states at vortex cores or at impurity locations. We use a superconducting STM tip for enhanced energy resolution, which ena… ▽ More By using scanning tunneling microscopy (STM) we find and characterize dispersive, energy-symmetric in-gap states in the iron-based superconductor $\mathrm{FeTe}_{0.55}\mathrm{Se}_{0.45}$, a material that exhibits signatures of topological superconductivity, and Majorana bound states at vortex cores or at impurity locations. We use a superconducting STM tip for enhanced energy resolution, which enables us to show that impurity states can be tuned through the Fermi level with varying tip-sample distance. We find that the impurity state is of the Yu-Shiba-Rusinov (YSR) type, and argue that the energy shift is caused by the low superfluid density in $\mathrm{FeTe}_{0.55}\mathrm{Se}_{0.45}$, which allows the electric field of the tip to slightly penetrate the sample. We model the newly introduced tip-gating scenario within the single-impurity Anderson model and find good agreement to the experimental data. △ Less

Submitted 23 June, 2020; originally announced June 2020.

Report number: NBI QDEV CMT 2021

Journal ref: Nature Communications 12, 298 (2021)

Asteroseismic masses of four evolved planet-hosting stars using SONG and TESS: resolving the retired A-star mass controversy

Authors: Sai Prathyusha Malla, Dennis Stello, Daniel Huber, Benjamin T. Montet, Timothy R. Bedding, Mads Fredslund Andersen, Frank Grundahl, Jens Jessen-Hansen, Daniel R. Hey, Pere L. Palle, Licai Deng, Chunguang Zhang, Xiaodian Chen, James Lloyd, Victoria Antoci

Abstract: The study of planet occurrence as a function of stellar mass is important for a better understanding of planet formation. Estimating stellar mass, especially in the red giant regime, is difficult. In particular, stellar masses of a sample of evolved planet-hosting stars based on spectroscopy and grid-based modelling have been put to question over the past decade with claims they were overestimated… ▽ More The study of planet occurrence as a function of stellar mass is important for a better understanding of planet formation. Estimating stellar mass, especially in the red giant regime, is difficult. In particular, stellar masses of a sample of evolved planet-hosting stars based on spectroscopy and grid-based modelling have been put to question over the past decade with claims they were overestimated. Although efforts have been made in the past to reconcile this dispute using asteroseismology, results were inconclusive. In an attempt to resolve this controversy, we study four more evolved planet-hosting stars in this paper using asteroseismology, and we revisit previous results to make an informed study of the whole ensemble in a self-consistent way. For the four new stars, we measure their masses by locating their characteristic oscillation frequency, $\mathrmν_{\mathrm{max}}$, from their radial velocity time series observed by SONG. For two stars, we are also able to measure the large frequency separation, $\mathrm{Δν}$, helped by extended SONG single-site and dual-site observations and new TESS observations. We establish the robustness of the $\mathrmν_{\mathrm{max}}$-only-based results by determining the stellar mass from $\mathrm{Δν}$, and from both $\mathrm{Δν}$ and $\mathrmν_{\mathrm{max}}$. We then compare the seismic masses of the full ensemble of 16 stars with the spectroscopic masses from three different literature sources. We find an offset between the seismic and spectroscopic mass scales that is mass-dependent, suggesting that the previously claimed overestimation of spectroscopic masses only affects stars more massive than about 1.6 M$_\mathrm{\odot}$. △ Less

Submitted 22 June, 2020; v1 submitted 13 June, 2020; originally announced June 2020.

Comments: 13 pages, 12 figures, Accepted for publication in the Main Journal of MNRAS

Journal ref: MNRAS 496 2020 5423-5435

Sun-like Oscillations in the Population II giant HD 122563

Authors: Orlagh Creevey, Frédéric Thévenin, Frank Grundahl, Enrico Corsaro, Mads F. Andersen, Victoria Antoci, Lionel Bigot, Remo Collet, Pere L. Pallé, Bernard Pichon, David Salabert

Abstract: We have been monitoring the metal-poor Population II giant, HD 122563, for radial velocity variations since 2016 using the SONG telescope on Tenerife. We have detected the global seismic quantity, numax, which provides information related to the stellar parameters. By combining these data with complementary data, we derive a new precise surface gravity, radius and distance to the star. Our results… ▽ More We have been monitoring the metal-poor Population II giant, HD 122563, for radial velocity variations since 2016 using the SONG telescope on Tenerife. We have detected the global seismic quantity, numax, which provides information related to the stellar parameters. By combining these data with complementary data, we derive a new precise surface gravity, radius and distance to the star. Our results are corroborated by using the parallax from Gaia DR2. We present these results and some of their implications. △ Less

Submitted 29 May, 2020; v1 submitted 28 May, 2020; originally announced May 2020.

Comments: Proceedings of the French Astronomy Society Meeting, SF2A, June 2019; v2: typo in reported numax in Table, authors corrected

Journal ref: 2019sf2a.conf..425C

Stellar Variability in a Forming Massive Star Cluster

Authors: Yu Cheng, Morten Andersen, Jonathan C. Tan

Abstract: We present a near-infrared (NIR) variability analysis for an 6\arcmin $\times$ 6\arcmin region, which encompasses the massive protocluster G286.21+0.17. The total sample comprises more than 5000 objects, of which 562 show signs of a circumstellar disk based on their infrared colors. The data includes HST observations taken in two epochs separated by 3 years in the F110W and F160W bands. 363 object… ▽ More We present a near-infrared (NIR) variability analysis for an 6\arcmin $\times$ 6\arcmin region, which encompasses the massive protocluster G286.21+0.17. The total sample comprises more than 5000 objects, of which 562 show signs of a circumstellar disk based on their infrared colors. The data includes HST observations taken in two epochs separated by 3 years in the F110W and F160W bands. 363 objects (7% of the sample) exhibit NIR variability at a significant level (Stetson index >1.7), and a higher variability fraction (14%) is found for the young stellar objects (YSOs) with disk excesses. We identified 4 high amplitude (>0.6 mag) variables seen in both NIR bands. Follow up and archival observations of the most variable object in this survey (G286.2032+0.1740) reveal a rising light curve over 8 years from 2011 to 2019, with a K band brightening of 3.5 mag. Overall the temporal behavior of G286.2032+0.1740 resembles that of typical FU Ori objects, however its pre-burst luminosity indicates it has a very low mass ($<0.12\:M_\odot$), making it an extreme case of an outburst event that is still ongoing. △ Less

Submitted 17 May, 2020; originally announced May 2020.

Comments: 13 pages, 8 figures, accepted by ApJ

Gas Kinematics of the Massive Protocluster G286.21+0.17 Revealed by ALMA

Authors: Yu Cheng, Jonathan C. Tan, Mengyao Liu, Wanggi Lim, Morten Andersen

Abstract: We study the gas kinematics and dynamics of the massive protocluster G286.21+0.17 with the Atacama Large Millimeter/submillimeter Array using spectral lines of $C^{18}O$(2-1), $N_2D^+$(3-2), $DCO^+$(3-2) and $DCN$(3-2). On the parsec clump scale, $C^{18}O$ emission appears highly filamentary around the systemic velocity. $N_2D^+$ and $DCO^+$ are more closely associated with the dust continuum.… ▽ More We study the gas kinematics and dynamics of the massive protocluster G286.21+0.17 with the Atacama Large Millimeter/submillimeter Array using spectral lines of $C^{18}O$(2-1), $N_2D^+$(3-2), $DCO^+$(3-2) and $DCN$(3-2). On the parsec clump scale, $C^{18}O$ emission appears highly filamentary around the systemic velocity. $N_2D^+$ and $DCO^+$ are more closely associated with the dust continuum. $DCN$ is strongly concentrated towards the protocluster center, where no or only weak detection is seen for $N_2D^+$ and $DCO^+$, possibly due to this region being at a relatively evolved evolutionary stage. Spectra of 76 continuum defined dense cores, typically a few 1000 AU in size, are analysed to measure their centroid velocities and internal velocity dispersions. There are no statistically significant velocity offsets of the cores among the different dense gas tracers. Furthermore, the majority (71\%) of the dense cores have subthermal velocity offsets with respect to their surrounding, lower density $C^{18}O$ emitting gas. Within the uncertainties, the dense cores in G286 show internal kinematics that are consistent with being in virial equilibrium. On clumps scales, the core to core velocity dispersion is also similar to that required for virial equilibrium in the protocluster potential. However, the distribution in velocity of the cores is largely composed of two spatially distinct groups, which indicates that the dense molecular gas has not yet relaxed to virial equilibrium, perhaps due to there being recent/continuous infall into the system. △ Less

Submitted 17 May, 2020; originally announced May 2020.

Comments: 26 pages, 17 figures, accepted by ApJ

Journal ref: The Astrophysical Journal, Volume 894, 2020, Number 2, Page 87

Doppler imaging and differential rotation of $σ^{2}$ Coronae Borealis using SONG

Authors: Yue Xiang, Shenghong Gu, A. Collier Cameron, J. R. Barnes, J. Christensen-Dalsgaard, F. Grundahl, V. Antoci, M. F. Andersen, P. L. Pallé

Abstract: We present new Doppler images of both components of the double-lined binary $σ^{2}$ CrB, based on the high-resolution spectroscopic data collected during 11 nights in 2015 March--April. The observed spectra form two independent data sets with sufficient phase coverage. We apply the least-squares deconvolution to all observed spectra to obtain high signal-to-noise mean profiles, from which we deriv… ▽ More We present new Doppler images of both components of the double-lined binary $σ^{2}$ CrB, based on the high-resolution spectroscopic data collected during 11 nights in 2015 March--April. The observed spectra form two independent data sets with sufficient phase coverage. We apply the least-squares deconvolution to all observed spectra to obtain high signal-to-noise mean profiles, from which we derive the Doppler images of both components of $σ^{2}$ CrB simultaneously. The surfaces of both F9 and G0 components are dominated by pronounced polar spots. The F9 component exhibits a weak spot at latitude 30$^{\circ}$ and its mid-to-low latitudes are relatively featureless. The G0 star shows an extended spot structure at latitude 30$^{\circ}$, and its surface spot coverage is larger than that of the F9 star, which suggests a higher level of magnetic activity. With the cross-correlation method, we derive a solar-like surface differential rotation on the G0 star of $σ^{2}$ CrB for the first time, and the surface shear rate is $ΔΩ= 0.180 \pm 0.004$ rad d$^{-1}$ and $α= ΔΩ/ Ω_{eq} = 0.032 \pm 0.001$. We do not obtain a clear surface shear law for the F9 star due to the lack of mid-to-low latitude features, but detect a systematic longitude shift of high-latitude spots, which indicates a slower rotation with respect to the co-rotating frame. △ Less

Submitted 6 May, 2020; originally announced May 2020.

Comments: 9 figures, 2 tables, accepted for publication in ApJ

HST survey of the Orion Nebula Cluster in the H$_2$O 1.4 $μ$m absorption band: I. A census of substellar and planetary mass objects

Authors: Massimo Robberto, Mario Gennaro, Maria Giulia Ubeira Gabellini, Lynne A. Hillenbrand, Camilla Pacifici, Leonardo Ubeda, Morten Andersen, Travis Barman, Andrea Bellini, Nicola Da Rio, Selma E. de Mink, Giuseppe Lodato, Carlo Felice Manara, Imants Platais, Laurent Pueyo, Giovanni M. Strampelli, Jonathan C. Tan, Leonardo Testi

Abstract: In order to obtain a complete census of the stellar and sub-stellar population, down to a few M$_{Jup}$ in the $\sim1$ Myr old Orion Nebula Cluster, we used the infrared channel of the Wide Field Camera 3 of the Hubble Space Telescope with the F139M and F130N filters. These bandpasses correspond to the $1.4 μ$m H$_2$O absorption feature and an adjacent line-free continuum region. Out of $4,504$ de… ▽ More In order to obtain a complete census of the stellar and sub-stellar population, down to a few M$_{Jup}$ in the $\sim1$ Myr old Orion Nebula Cluster, we used the infrared channel of the Wide Field Camera 3 of the Hubble Space Telescope with the F139M and F130N filters. These bandpasses correspond to the $1.4 μ$m H$_2$O absorption feature and an adjacent line-free continuum region. Out of $4,504$ detected sources, $3,352$ (about $75\%$) appear fainter than m$_{130}=14$ (Vega mag) in the F130N filter, a brightness corresponding to the hydrogen-burning limit mass (M$\simeq 0.072 M_\odot$) at $\sim 1$ Myr. Of these, however, only $742$ sources have a negative F130M-139N color index, indicative of the presence of H$_2$O vapor in absorption, and can therefore be classified as bona-fide M and L dwarfs, with effective temperatures T$\lesssim 2850$ K at an assumed $1$ Myr cluster age. On our color-magnitude diagram, this population of sources with H$_2$O absorption appears clearly distinct from the larger background population of highly reddened stars and galaxies with positive F130M-F139N color index, and can be traced down to the sensitivity limit of our survey, m$_{130}\simeq 21.5$, corresponding to a $1$ Myr old $\simeq 3 $M$_{Jup}$, planetary mass object under about 2 magnitudes of visual extinction. Theoretical models of the BT-Settl family predicting substellar isochrones of $1, 2$ and $3$ Myr (down to $\sim 1 $M$_{Jup}$) fail to reproduce the observed H$_2$O color index at M$\lesssim 20 $M$_{Jup}$. We perform a Bayesian analysis to determine extinction, mass and effective temperature of each sub-stellar member of our sample, together with its membership probability. △ Less

Submitted 28 April, 2020; originally announced April 2020.

Comments: Accepted for publication in the Astrophysical Journal. The resolution of several figures has been downgraded to comply with the size limit of arXiv submissions

Correlation-induced insulating topological phases at charge neutrality in twisted bilayer graphene

Authors: Yuan Da Liao, Jian Kang, Clara N. Breiø, Xiao Yan Xu, Han-Qing Wu, Brian M. Andersen, Rafael M. Fernandes, Zi Yang Meng

Abstract: Twisted bilayer graphene (TBG) provides a unique framework to elucidate the interplay between strong correlations and topological phenomena in two-dimensional systems. The existence of multiple electronic degrees of freedom -- charge, spin, and valley -- gives rise to a plethora of possible ordered states and instabilities. Identifying which of them are realized in the regime of strong correlation… ▽ More Twisted bilayer graphene (TBG) provides a unique framework to elucidate the interplay between strong correlations and topological phenomena in two-dimensional systems. The existence of multiple electronic degrees of freedom -- charge, spin, and valley -- gives rise to a plethora of possible ordered states and instabilities. Identifying which of them are realized in the regime of strong correlations is fundamental to shed light on the nature of the superconducting and correlated insulating states observed in the TBG experiments. Here, we use unbiased, sign-problem-free quantum Monte Carlo simulations to solve an effective interacting lattice model for TBG at charge neutrality. Besides the usual cluster Hubbard-like repulsion, this model also contains an assisted hopping interaction that emerges due to the non-trivial topological properties of TBG. Such a non-local interaction fundamentally alters the phase diagram at charge neutrality, gapping the Dirac cones even for infinitesimally small interaction. As the interaction strength increases, a sequence of different correlated insulating phases emerge, including a quantum valley Hall state with topological edge states, an intervalley-coherent insulator, and a valence bond solid. The charge-neutrality correlated insulating phases discovered here provide the sought-after reference states needed for a comprehensive understanding of the insulating states at integer fillings and the proximate superconducting states of TBG. △ Less

Submitted 23 January, 2021; v1 submitted 26 April, 2020; originally announced April 2020.

Comments: 15 pages, 9 figures, 2 tables

Journal ref: Phys. Rev. X 11, 011014 (2021)

Practical Hilbert space approximate Bayesian Gaussian processes for probabilistic programming

Authors: Gabriel Riutort-Mayol, Paul-Christian Bürkner, Michael R. Andersen, Arno Solin, Aki Vehtari

Abstract: Gaussian processes are powerful non-parametric probabilistic models for stochastic functions. However, the direct implementation entails a complexity that is computationally intractable when the number of observations is large, especially when estimated with fully Bayesian methods such as Markov chain Monte Carlo. In this paper, we focus on a low-rank approximate Bayesian Gaussian processes, based… ▽ More Gaussian processes are powerful non-parametric probabilistic models for stochastic functions. However, the direct implementation entails a complexity that is computationally intractable when the number of observations is large, especially when estimated with fully Bayesian methods such as Markov chain Monte Carlo. In this paper, we focus on a low-rank approximate Bayesian Gaussian processes, based on a basis function approximation via Laplace eigenfunctions for stationary covariance functions. The main contribution of this paper is a detailed analysis of the performance, and practical recommendations for how to select the number of basis functions and the boundary factor. Intuitive visualizations and recommendations, make it easier for users to improve approximation accuracy and computational performance. We also propose diagnostics for checking that the number of basis functions and the boundary factor are adequate given the data. The approach is simple and exhibits an attractive computational complexity due to its linear structure, and it is easy to implement in probabilistic programming frameworks. Several illustrative examples of the performance and applicability of the method in the probabilistic programming language Stan are presented together with the underlying Stan model code. △ Less

Submitted 22 March, 2022; v1 submitted 23 April, 2020; originally announced April 2020.

Comments: 27 pages, 18 figures

Theory of Strain-Induced Magnetic Order and Splitting of $T_c$ and $T_{\rm TRSB}$ in Sr$_2$RuO$_4$

Authors: Astrid T. Rømer, Andreas Kreisel, Marvin A. Müller, P. J. Hirschfeld, Ilya M. Eremin, Brian M. Andersen

Abstract: The internal structure of the superconducting state in Sr$_2$RuO$_4$ remains elusive at present, and exhibits evidence for time-reversal symmetry breaking. Recent muon spin relaxation measurements under uniaxial strain have revealed an increasing splitting between the superconducting critical temperature $T_c$ and the onset of time-reversal symmetry breaking $T_{\rm TRSB}$ with applied strain [Gri… ▽ More The internal structure of the superconducting state in Sr$_2$RuO$_4$ remains elusive at present, and exhibits evidence for time-reversal symmetry breaking. Recent muon spin relaxation measurements under uniaxial strain have revealed an increasing splitting between the superconducting critical temperature $T_c$ and the onset of time-reversal symmetry breaking $T_{\rm TRSB}$ with applied strain [Grinenko et al., arXiv:2001.08152]. In addition, static magnetic order is induced by the uniaxial strain beyond $\sim$1 GPa, indicating that unstrained Sr$_2$RuO$_4$ is close to a magnetic quantum critical point. Here, we perform a theoretical study of the magnetic susceptibility and the associated pairing structure as a function of uniaxial strain. It is found that the recent muon relaxation data can be qualitatively explained from the perspective of spin-fluctuation mediated pairing and the associated strain-dependence of accidentally degenerate pair states in unstrained Sr$_2$RuO$_4$. In addition, while unstrained Sr$_2$RuO$_4$ features mainly $(2π/3,2π/3)$ magnetic fluctuations, uniaxial strain promotes $(π,\pmπ/2)$ magnetism. △ Less

Submitted 30 March, 2020; originally announced March 2020.

Comments: 8 pages, 4 figures

Report number: CMT NBI 2020

Journal ref: Phys. Rev. B 102, 054506 (2020)

Preferential Batch Bayesian Optimization

Authors: Eero Siivola, Akash Kumar Dhaka, Michael Riis Andersen, Javier Gonzalez, Pablo Garcia Moreno, Aki Vehtari

Abstract: Most research in Bayesian optimization (BO) has focused on \emph{direct feedback} scenarios, where one has access to exact values of some expensive-to-evaluate objective. This direction has been mainly driven by the use of BO in machine learning hyper-parameter configuration problems. However, in domains such as modelling human preferences, A/B tests, or recommender systems, there is a need for me… ▽ More Most research in Bayesian optimization (BO) has focused on \emph{direct feedback} scenarios, where one has access to exact values of some expensive-to-evaluate objective. This direction has been mainly driven by the use of BO in machine learning hyper-parameter configuration problems. However, in domains such as modelling human preferences, A/B tests, or recommender systems, there is a need for methods that can replace direct feedback with \emph{preferential feedback}, obtained via rankings or pairwise comparisons. In this work, we present preferential batch Bayesian optimization (PBBO), a new framework that allows finding the optimum of a latent function of interest, given any type of parallel preferential feedback for a group of two or more points. We do so by using a Gaussian process model with a likelihood specially designed to enable parallel and efficient data collection mechanisms, which are key in modern machine learning. We show how the acquisitions developed under this framework generalize and augment previous approaches in Bayesian optimization, expanding the use of these techniques to a wider range of domains. An extensive simulation study shows the benefits of this approach, both with simulated functions and four real data sets. △ Less

Submitted 31 August, 2021; v1 submitted 25 March, 2020; originally announced March 2020.

Comments: 6 pages + 7 pages in supplementary material

The TESS-Keck Survey I: A Warm Sub-Saturn-mass Planet and a Caution about Stray Light in TESS Cameras

Authors: Paul A. Dalba, Arvind F. Gupta, Joseph E. Rodriguez, Diana Dragomir, Chelsea X. Huang, Stephen R. Kane, Samuel N. Quinn, Allyson Bieryla, Gilbert A. Esquerdo, Benjamin J. Fulton, Nicholas Scarsdale, Natalie M. Batalha, Corey Beard, Aida Behmard, Ashley Chontos, Ian J. M. Crossfield, Courtney D. Dressing, Steven Giacalone, Michelle L. Hill, Lea A. Hirsch, Andrew W. Howard, Daniel Huber, Howard Isaacson, Molly Kosiarek, Jack Lubin , et al. (41 additional authors not shown)

Abstract: We report the detection of a Saturn-size exoplanet orbiting HD 332231 (TOI 1456) in light curves from the Transiting Exoplanet Survey Satellite (TESS). HD 332231, an F8 dwarf star with a V-band magnitude of 8.56, was observed by TESS in Sectors 14 and 15. We detect a single-transit event in the Sector 15 presearch data conditioning (PDC) light curve. We obtain spectroscopic follow-up observations… ▽ More We report the detection of a Saturn-size exoplanet orbiting HD 332231 (TOI 1456) in light curves from the Transiting Exoplanet Survey Satellite (TESS). HD 332231, an F8 dwarf star with a V-band magnitude of 8.56, was observed by TESS in Sectors 14 and 15. We detect a single-transit event in the Sector 15 presearch data conditioning (PDC) light curve. We obtain spectroscopic follow-up observations of HD 332231 with the Automated Planet Finder, Keck I, and SONG telescopes. The orbital period we infer from the radial velocity (RV) observations leads to the discovery of another transit in Sector 14 that was masked by PDC due to scattered light contamination. A joint analysis of the transit and RV data confirms the planetary nature of HD 332231 b, a Saturn-size ($0.867^{+0.027}_{-0.025} \; R_{\rm J}$), sub-Saturn-mass ($0.244\pm0.021 \; M_{\rm J}$) exoplanet on a 18.71 day circular orbit. The low surface gravity of HD 332231 b and the relatively low stellar flux it receives make it a compelling target for transmission spectroscopy. Also, the stellar obliquity is likely measurable via the Rossiter-McLaughlin effect, an exciting prospect given the 0.14 au orbital separation of HD 332231 b. The spectroscopic observations do not provide substantial evidence for any additional planets in the HD 332231 system, but continued RV monitoring is needed to further characterize this system. We also predict that the frequency and duration of masked data in the PDC light curves for TESS Sectors 14-16 could hide transits of some exoplanets with orbital periods between 10.5 and 17.5 days. △ Less

Submitted 29 April, 2020; v1 submitted 23 March, 2020; originally announced March 2020.

Comments: 17 pages, 9 figures, 4 tables. Published in the Astronomical Journal

Outflow from Outer-arm Starburst in a Grazing Collision between Galaxies

Authors: Michele Kaufman, Bruce G. Elmegreen, Morten Andersen, Debra Meloy Elmegreen, Curtis Struck, Frederic Bournaud, Elias Brinks, James C. McGarry

Abstract: New ALMA CO, HCO+, and 100 GHz continuum observations and Gemini NIFS K-band spectra are combined with previous radio to X-ray data to study Feature i, a highly luminous starburst clump on an outer arm of the interacting galaxy NGC 2207. This clump has an optically-opaque dust cone extending out of its 170 pc core. The measured CO accounts for the dark cone extinction if almost all the gas and dus… ▽ More New ALMA CO, HCO+, and 100 GHz continuum observations and Gemini NIFS K-band spectra are combined with previous radio to X-ray data to study Feature i, a highly luminous starburst clump on an outer arm of the interacting galaxy NGC 2207. This clump has an optically-opaque dust cone extending out of its 170 pc core. The measured CO accounts for the dark cone extinction if almost all the gas and dust there is in front of the star clusters. An associated approaching CO outflow has v_z ~ 16 km/s, an estimated molecular mass 8 x 10^6 Msun, and rises to heights ~ 0.9 kpc. A receding CO outflow on the far side with v_z ~ 28 km/s is less extensive. The observed star formation in the core over 10 Myr can supply the dark cone kinetic energy via supernovae and stellar winds. Other signs of intense activity are variable radio continuum suggesting an embedded supernova or other outburst, X-ray emission, and Br-gamma and He I lines with 82 km/s line widths. According to previous models, the retrograde encounter suffered by NGC 2207 caused loss of angular momentum. This compressed its outer disk. We suggest that the resulting inward crashing gas stream hit a massive HI cloud as it was being compressed by a spiral density wave and triggered the observed starburst. △ Less

Submitted 19 March, 2020; originally announced March 2020.

Comments: accepted by AJ

Non-local correlations in Iron Pnictides and Chalcogenides

Authors: Shinibali Bhattacharyya, Kristofer Björnson, Karim Zantout, Daniel Steffensen, Laura Fanfarillo, Andreas Kreisel, Roser Valentí, Brian M. Andersen, P. J. Hirschfeld

Abstract: Deviations of low-energy electronic structure of iron-based superconductors from density functional theory predictions have been parametrized in terms of band- and orbital-dependent mass renormalizations and energy shifts. The former have typically been described in terms of a local self-energy within the framework of dynamical mean field theory, while the latter appears to require non-local effec… ▽ More Deviations of low-energy electronic structure of iron-based superconductors from density functional theory predictions have been parametrized in terms of band- and orbital-dependent mass renormalizations and energy shifts. The former have typically been described in terms of a local self-energy within the framework of dynamical mean field theory, while the latter appears to require non-local effects due to interband scattering. By calculating the renormalized bandstructure in both random phase approximation (RPA) and the two-particle self-consistent approximation (TPSC), we show that correlations in pnictide systems like LaFeAsO and LiFeAs can be described rather well by a non-local self-energy. In particular, Fermi pocket shrinkage as seen in experiment occurs due to repulsive interband finite-energy scattering. For the canonical iron chalcogenide system FeSe in its bulk tetragonal phase, the situation is however more complex since even including momentum-dependent band renormalizations cannot explain experimental findings. We propose that the long-range Coulomb interaction may play an important role in band-structure renormalization in FeSe. We further compare our evaluations of non-local quasiparticle scattering lifetime within RPA and TPSC with experimental data for LiFeAs. △ Less

Submitted 17 July, 2020; v1 submitted 3 March, 2020; originally announced March 2020.

Comments: 12 Pages, 8 figures

Journal ref: Phys. Rev. B 102, 035109 (2020)

Spin-orbit quantum impurity in a topological kagome magnet

Authors: Jia-Xin Yin, Nana Shumiya, Yuxiao Jiang, Huibin Zhou, Gennevieve Macam, Songtian S. Zhang, Hano Omar Mohammad Sura, Zijia Cheng, Zurab Guguchia, Yangmu Li, Qi Wang, Maksim Litskevich, Ilya Belopolski, Xian Yang, Tyler A. Cochran, Guoqing Chang, Qi Zhang, Brian M. Andersen, Zhi-Quan Huang, Feng-Chuan Chuang, Hsin Lin, Hechang Lei, Ziqiang Wang, Shuang Jia, M. Zahid Hasan

Abstract: Quantum states induced by single-atomic-impurities are the current frontier of material and information science. Recently the spin-orbit coupled correlated kagome magnets are emerging as a new class of topological quantum materials, although the effect of single-atomic impurities remains unexplored. Here we use state-of-the-art scanning tunneling microscopy/spectroscopy (STM/S) to study the atomic… ▽ More Quantum states induced by single-atomic-impurities are the current frontier of material and information science. Recently the spin-orbit coupled correlated kagome magnets are emerging as a new class of topological quantum materials, although the effect of single-atomic impurities remains unexplored. Here we use state-of-the-art scanning tunneling microscopy/spectroscopy (STM/S) to study the atomic indium impurity in a topological kagome magnet Co3Sn2S2, which is designed to support the spin-orbit quantum state. We find each impurity features a strongly localized bound state. Our systematic magnetization-polarized tunneling probe reveals its spin-down polarized nature with an unusual moment of -5uB, indicative of additional orbital magnetization. As the separation between two impurities progressively shrinks, their respective bound states interact and form quantized molecular orbital states. The molecular orbital of three neighboring impurities further exhibits an intriguing splitting owing to the combination of geometry, magnetism, and spin-orbit coupling, analogous to the splitting of the topological Weyl fermion line12,19. Our work demonstrates the quantum-level interplay between magnetism and spin-orbit coupling at an individual atomic impurity, which provides insights into the emergent impurity behavior in a topological kagome magnet and the potential of spin-orbit quantum impurities for information science. △ Less

Submitted 26 February, 2020; originally announced February 2020.

Fluctuation-driven superconductivity in Sr$_2$RuO$_4$ from weak repulsive interactions

Authors: Astrid T. Rømer, Brian M. Andersen

Abstract: We provide results for the leading superconducting instabilities for a model pertaining to Sr$_2$RuO$_4$ obtained within spin-fluctuation mediated superconductivity in the very weak-coupling limit. The theory incorporates spin-orbit coupling (SOC) effects both in the band structure and in the pairing kernel in the form of associated magnetic anisotropies. The leading superconducting phase is found… ▽ More We provide results for the leading superconducting instabilities for a model pertaining to Sr$_2$RuO$_4$ obtained within spin-fluctuation mediated superconductivity in the very weak-coupling limit. The theory incorporates spin-orbit coupling (SOC) effects both in the band structure and in the pairing kernel in the form of associated magnetic anisotropies. The leading superconducting phase is found to be $d_{x^2-y^2}$ and a nodal $s$-wave state. However, the odd-parity helical solution can become leading either for small SOC and Hund's coupling $J$ in the weak $U$-limit, or in the opposite limit with large SOC and $J$ at larger values of the Hubbard-$U$. The odd-parity chiral solution is never found to be leading. Finally we discuss the form of the resulting superconducting spectral gaps in the different explored parameter regimes. △ Less

Submitted 30 January, 2020; originally announced January 2020.

Comments: 5 pages, 5 figures

Report number: NBI CMT 2020

Journal ref: Modern Physics Letters B, Vol. 34, No. 19n20, 2040052 (2020)

Reducing bias in dual flip angle T1-mapping in human brain at 7T

Authors: H. Olsson, M. Andersen, R. Wirestam, J. Lätt, G. Helms

Abstract: Purpose: To address the systematic bias in whole-brain dual flip angle (DFA) T1-mapping at 7T by optimizing the flip angle pair and carefully selecting RF pulse shape and duration. Theory and Methods: Spoiled gradient echoes can be used to estimate whole-brain maps of T1. This can be accomplished by using only two acquisitions with different flip angles, i.e., a DFA-based approach. Although DFA-ba… ▽ More Purpose: To address the systematic bias in whole-brain dual flip angle (DFA) T1-mapping at 7T by optimizing the flip angle pair and carefully selecting RF pulse shape and duration. Theory and Methods: Spoiled gradient echoes can be used to estimate whole-brain maps of T1. This can be accomplished by using only two acquisitions with different flip angles, i.e., a DFA-based approach. Although DFA-based T1-mapping is seemingly straightforward to implement, it is sensitive to bias caused by incomplete spoiling and incidental magnetization transfer (MT) effects. Further bias is introduced by the increased B0 and B1+ inhomogeneities at 7T. Experiments were performed to determine the optimal flip angle pair and appropriate RF pulse shape and duration. Obtained T1 estimates were validated using inversion recovery prepared EPI and compared to literature values. A multi-echo readout was used to increase SNR, enabling quantification of R2* and susceptibility, X. Results: Incomplete spoiling was observed above a local flip angle of approximately 20 degrees. An asymmetric gauss-filtered sinc pulse with a constant duration of 700 us showed a sufficiently flat frequency response profile to avoid incomplete excitation in areas with high B0 offsets. A pulse duration of 700 us minimized effects from incidental MT. Conclusion: When performing DFA-based T1-mapping one should (i) limit the higher flip angle to avoid incomplete spoiling, (ii) use a RF pulse shape insensitive to B0 inhomogeneities and (iii) apply a constant RF pulse duration, balanced to minimize incidental MT. △ Less

Submitted 24 January, 2020; originally announced January 2020.

Comments: 24 pages, 8 figures

Direct Measurements of Collisional Dynamics in Cold Atom Triads

Authors: L. A. Reynolds, E. Schwartz, U. Ebling, M. Weyland, J. Brand, M. F. Andersen

Abstract: The introduction of optical tweezers for trapping atoms has opened remarkable opportunities for manipulating few-body systems. Here, we present the first bottom-up assembly of atom triads. We directly observe atom loss through inelastic collisions at the single event level, overcoming the substantial challenge in many-atom experiments of distinguishing one-, two-, and three-particle processes. We… ▽ More The introduction of optical tweezers for trapping atoms has opened remarkable opportunities for manipulating few-body systems. Here, we present the first bottom-up assembly of atom triads. We directly observe atom loss through inelastic collisions at the single event level, overcoming the substantial challenge in many-atom experiments of distinguishing one-, two-, and three-particle processes. We measure a strong suppression of three-body loss, which is not fully explained by the presently availably theory for three-body processes. The suppression of losses could indicate the presence of local anti-correlations due to the interplay of attractive short range interactions and low dimensional confinement. Our methodology opens a promising pathway in experimental few-body dynamics. △ Less

Submitted 15 January, 2020; originally announced January 2020.

Comments: 11 pages, 4 figures, includes Supplemental Material

Journal ref: Phys. Rev. Lett. 124, 073401 (2020)

Leave-One-Out Cross-Validation for Bayesian Model Comparison in Large Data

Authors: Måns Magnusson, Michael Riis Andersen, Johan Jonasson, Aki Vehtari

Abstract: Recently, new methods for model assessment, based on subsampling and posterior approximations, have been proposed for scaling leave-one-out cross-validation (LOO) to large datasets. Although these methods work well for estimating predictive performance for individual models, they are less powerful in model comparison. We propose an efficient method for estimating differences in predictive performa… ▽ More Recently, new methods for model assessment, based on subsampling and posterior approximations, have been proposed for scaling leave-one-out cross-validation (LOO) to large datasets. Although these methods work well for estimating predictive performance for individual models, they are less powerful in model comparison. We propose an efficient method for estimating differences in predictive performance by combining fast approximate LOO surrogates with exact LOO subsampling using the difference estimator and supply proofs with regards to scaling characteristics. The resulting approach can be orders of magnitude more efficient than previous approaches, as well as being better suited to model comparison. △ Less

Submitted 3 January, 2020; originally announced January 2020.

Journal ref: Proceedings of the 23rd International Conference on Artificial Intelligence and Statistics (AISTATS), PMLR 108:341-351, 2020

Gaussian process with derivative information for the analysis of the sunlight adverse effects on color of rock art paintings

Authors: Gabriel Riutort-Mayol, Michael Riis Andersen, Aki Vehtari, José Luis Lerma

Abstract: Microfading Spectrometry (MFS) is a method for assessing light sensitivity color (spectral) variations of cultural heritage objects. The MFS technique provides measurements of the surface under study, where each point of the surface gives rise to a time-series that represents potential spectral (color) changes due to sunlight exposition over time. Color fading is expected to be non-decreasing as a… ▽ More Microfading Spectrometry (MFS) is a method for assessing light sensitivity color (spectral) variations of cultural heritage objects. The MFS technique provides measurements of the surface under study, where each point of the surface gives rise to a time-series that represents potential spectral (color) changes due to sunlight exposition over time. Color fading is expected to be non-decreasing as a function of time and stabilize eventually. These properties can be expressed in terms of the partial derivatives of the functions. We propose a spatio-temporal model that takes this information into account by jointly modeling the spatio-temporal process and its derivative process using Gaussian processes (GPs). We fitted the proposed model to MFS data collected from the surface of prehistoric rock art paintings. A multivariate covariance function in a GP allows modeling trichromatic image color variables jointly with spatial distances and time points variables as inputs to evaluate the covariance structure of the data. We demonstrated that the colorimetric variables are useful for predicting the color fading time-series for new unobserved spatial locations. Furthermore, constraining the model using derivative sign observations for monotonicity was shown to be beneficial in terms of both predictive performance and application-specific interpretability. △ Less

Submitted 7 November, 2019; originally announced November 2019.

Comments: arXiv admin note: substantial text overlap with arXiv:1910.12575

Machine Learning-Based Analysis of Sperm Videos and Participant Data for Male Fertility Prediction

Authors: Steven A. Hicks, Jorunn M. Andersen, Oliwia Witczak, Vajira Thambawita, Påll Halvorsen, Hugo L. Hammer, Trine B. Haugen, Michael A. Riegler

Abstract: Methods for automatic analysis of clinical data are usually targeted towards a specific modality and do not make use of all relevant data available. In the field of male human reproduction, clinical and biological data are not used to its fullest potential. Manual evaluation of a semen sample using a microscope is time-consuming and requires extensive training. Furthermore, the validity of manual… ▽ More Methods for automatic analysis of clinical data are usually targeted towards a specific modality and do not make use of all relevant data available. In the field of male human reproduction, clinical and biological data are not used to its fullest potential. Manual evaluation of a semen sample using a microscope is time-consuming and requires extensive training. Furthermore, the validity of manual semen analysis has been questioned due to limited reproducibility, and often high inter-personnel variation. The existing computer-aided sperm analyzer systems are not recommended for routine clinical use due to methodological challenges caused by the consistency of the semen sample. Thus, there is a need for an improved methodology. We use modern and classical machine learning techniques together with a dataset consisting of 85 videos of human semen samples and related participant data to automatically predict sperm motility. Used techniques include simple linear regression and more sophisticated methods using convolutional neural networks. Our results indicate that sperm motility prediction based on deep learning using sperm motility videos is rapid to perform and consistent. The algorithms performed worse when participant data was added. In conclusion, machine learning-based automatic analysis may become a valuable tool in male infertility investigation and research. △ Less

Submitted 29 October, 2019; originally announced October 2019.

Comments: Preprint, accepted by Nature Scientific Reports for publication 24.10.2019

Quantum phase transition of correlated iron-based superconductivity in LiFe$_{1-x}$Co$_x$As

Authors: Jia-Xin Yin, Songtian S. Zhang, Guangyang Dai, Yuanyuan Zhao, Andreas Kreisel, Gennevieve Macam, Xianxin Wu, Hu Miao, Zhi-Quan Huang, Johannes H. J. Martiny, Brian M. Andersen, Nana Shumiya, Daniel Multer, Maksim Litskevich, Zijia Cheng, Xian Yang, Tyler A. Cochran, Guoqing Chang, Ilya Belopolski, Lingyi Xing, Xiancheng Wang, Yi Gao, Feng-Chuan Chuang, Hsin Lin, Ziqiang Wang , et al. (3 additional authors not shown)

Abstract: The interplay between unconventional Cooper pairing and quantum states associated with atomic scale defects is a frontier of research with many open questions. So far, only a few of the high-temperature superconductors allow this intricate physics to be studied in a widely tunable way. We use scanning tunneling microscopy (STM) to image the electronic impact of Co atoms on the ground state of the… ▽ More The interplay between unconventional Cooper pairing and quantum states associated with atomic scale defects is a frontier of research with many open questions. So far, only a few of the high-temperature superconductors allow this intricate physics to be studied in a widely tunable way. We use scanning tunneling microscopy (STM) to image the electronic impact of Co atoms on the ground state of the LiFe$_{1-x}$Co$_x$As system. We observe that impurities progressively suppress the global superconducting gap and introduce low energy states near the gap edge, with the superconductivity remaining in the strong-coupling limit. Unexpectedly, the fully opened gap evolves into a nodal state before the Cooper pair coherence is fully destroyed. Our systematic theoretical analysis shows that these new observations can be quantitatively understood by the nonmagnetic Born-limit scattering effect in a s$\pm$-wave superconductor, unveiling the driving force of the superconductor to metal quantum phase transition. △ Less

Submitted 24 October, 2019; originally announced October 2019.

Comments: 22 pages, 12 figures, includes Supplementary Materials. To appear in Phys. Rev. Lett

Journal ref: Phys. Rev. Lett. 123, 217004 (2019)

Uncertainty-aware Sensitivity Analysis Using Rényi Divergences

Authors: Topi Paananen, Michael Riis Andersen, Aki Vehtari

Abstract: For nonlinear supervised learning models, assessing the importance of predictor variables or their interactions is not straightforward because it can vary in the domain of the variables. Importance can be assessed locally with sensitivity analysis using general methods that rely on the model's predictions or their derivatives. In this work, we extend derivative based sensitivity analysis to a Baye… ▽ More For nonlinear supervised learning models, assessing the importance of predictor variables or their interactions is not straightforward because it can vary in the domain of the variables. Importance can be assessed locally with sensitivity analysis using general methods that rely on the model's predictions or their derivatives. In this work, we extend derivative based sensitivity analysis to a Bayesian setting by differentiating the Rényi divergence of a model's predictive distribution. By utilising the predictive distribution instead of a point prediction, the model uncertainty is taken into account in a principled way. Our empirical results on simulated and real data sets demonstrate accurate and reliable identification of important variables and interaction effects compared to alternative methods. △ Less

Submitted 9 March, 2021; v1 submitted 17 October, 2019; originally announced October 2019.

Journal ref: Proceedings of the Thirty-Seventh Conference on Uncertainty in Artificial Intelligence, PMLR 161:1185-1194, 2021

Precise radial velocities of giant stars XIII. A second Jupiter orbiting in 4:3 resonance in the 7 CMa system

Authors: R. Luque, T. Trifonov, S. Reffert, A. Quirrenbach, M. H. Lee, S. Albrecht, M. Fredslund Andersen, V. Antoci, F. Grundahl, C. Schwab, V. Wolthoff

Abstract: We report the discovery of a second planet orbiting the K giant star 7 CMa based on 166 high-precision radial velocities obtained with Lick, HARPS, UCLES and SONG. The periodogram analysis reveals two periodic signals of approximately 745 and 980 d, associated to planetary companions. A double-Keplerian orbital fit of the data reveals two Jupiter-like planets with minimum masses… ▽ More We report the discovery of a second planet orbiting the K giant star 7 CMa based on 166 high-precision radial velocities obtained with Lick, HARPS, UCLES and SONG. The periodogram analysis reveals two periodic signals of approximately 745 and 980 d, associated to planetary companions. A double-Keplerian orbital fit of the data reveals two Jupiter-like planets with minimum masses $m_b\sin i \sim 1.9 \,\mathrm{M_{J}}$ and $m_c\sin i \sim 0.9 \,\mathrm{M_{J}}$, orbiting at semi-major axes of $a_b \sim 1.75\,\mathrm{au}$ and $a_c \sim 2.15\,\mathrm{au}$, respectively. Given the small orbital separation and the large minimum masses of the planets close encounters may occur within the time baseline of the observations, thus, a more accurate N-body dynamical modeling of the available data is performed. The dynamical best-fit solution leads to collision of the planets and we explore the long-term stable configuration of the system in a Bayesian framework, confirming that 13% of the posterior samples are stable for at least 10 Myr. The result from the stability analysis indicates that the two-planets are trapped in a low-eccentricity 4:3 mean-motion resonance. This is only the third discovered system to be inside a 4:3 resonance, making it very valuable for planet formation and orbital evolution models. △ Less

Submitted 13 October, 2019; originally announced October 2019.

Comments: Accepted in A&A

Journal ref: A&A 631, A136 (2019)

Variation of shear moduli across superconducting phase transitions

Authors: Dimitri Labat, Panagiotis Kotetes, Brian M. Andersen, Indranil Paul

Abstract: We study how shear moduli of a correlated metal change across superconducting phase transitions. Using a microscopic theory we explain why for most classes of superconductors this change is small. The Fe-based and the A15 systems are notable exceptions where the change is boosted by five orders of magnitude. We show that this boost is a consequence of enhanced nematic correlation. The theory expla… ▽ More We study how shear moduli of a correlated metal change across superconducting phase transitions. Using a microscopic theory we explain why for most classes of superconductors this change is small. The Fe-based and the A15 systems are notable exceptions where the change is boosted by five orders of magnitude. We show that this boost is a consequence of enhanced nematic correlation. The theory explains the unusual temperature dependence of the orthorhombic shear and the back-bending of the nematic transition line in the superconducting phase of the Fe-based systems. △ Less

Submitted 7 October, 2019; originally announced October 2019.

Comments: 5 pages, 3 Figures

Report number: CMT NBI 2019

Journal ref: Phys. Rev. B 101, 144502 (2020)

Visualization of local magnetic moments emerging from impurities in the Hund's metal states of FeSe

Authors: Sang Yong Song, J. H. J. Martiny, A. Kreisel, B. M. Andersen, Jungpil Seo

Abstract: Understanding the origin of the magnetism of high temperature superconductors is crucial for establishing their unconventional pairing mechanism. Recently, theory predicts that FeSe is close to a magnetic quantum critical point, and thus weak perturbations such as impurities could induce local magnetic moments. To elucidate such quantum instability, we have employed scanning tunneling microscopy a… ▽ More Understanding the origin of the magnetism of high temperature superconductors is crucial for establishing their unconventional pairing mechanism. Recently, theory predicts that FeSe is close to a magnetic quantum critical point, and thus weak perturbations such as impurities could induce local magnetic moments. To elucidate such quantum instability, we have employed scanning tunneling microscopy and spectroscopy. In particular, we have grown FeSe film on superconducting Pb(111) using molecular beam epitaxy and investigated magnetic excitation caused by impurities in the proximity-induced superconducting gap of FeSe. Our study provides a deep insight into the origin of the magnetic ordering of FeSe by showing the way local magnetic moments develop in response to impurities near the magnetic quantum critical point. △ Less

Submitted 17 March, 2020; v1 submitted 30 September, 2019; originally announced September 2019.

Comments: 35 pages including main text and supplemental material

Report number: CMT NBI 2019

Journal ref: Phys. Rev. Lett. 124, 117001 (2020)

Effects of spin-orbit coupling on spin-fluctuation induced pairing in iron-based superconductors

Authors: Daniel D. Scherer, Brian M. Andersen

Abstract: We perform a theoretical study of the leading pairing instabilities and the associated superconducting gap functions within the spin-fluctuation mediated pairing scenario in the presence of spin-orbit coupling (SOC). Focussing on iron-based superconductors (FeSCs), our model Hamiltonian consists of a realistic density functional theory (DFT)-derived ten-band hopping term, spin-orbit coupling, and… ▽ More We perform a theoretical study of the leading pairing instabilities and the associated superconducting gap functions within the spin-fluctuation mediated pairing scenario in the presence of spin-orbit coupling (SOC). Focussing on iron-based superconductors (FeSCs), our model Hamiltonian consists of a realistic density functional theory (DFT)-derived ten-band hopping term, spin-orbit coupling, and electron-electron interactions included via the multi-orbital Hubbard-Hund Hamiltonian. We perform an extensive parameter sweep and investigate different doping regimes including cases with only hole- or only electron Fermi pockets. In addition, we explore two different bandstructures: a rather generic band derived for LaFeAsO but known to represent standard DFT-obtained bands for iron-based superconductors, and a band specifically tailored for FeSe which exhibits a notably different Fermi surface compared to the generic case. It is found that for the generic FeSCs band, even rather large SOC has negligible effect on the resulting gap structure; the $s_{+-}$ (pseudo-)spin singlet pairing remains strongly favored and SOC does not lead to any SOC-characteristic gap oscillations along the various Fermi surfaces. By contrast in the strongly hole-doped case featuring only hole-pockets around the $Γ$-point, the leading solution is $d$-wave pseudo-spin singlet, but with a notable SOC-driven tendency towards helical pseudo-spin triplet pairing, which may even become the leading instability. In the heavily electron doped situation, featuring only electron pockets centered around the $M$-point, the leading superconducting instabilities are pseudo-spin singlets with SOC favoring the $s$-wave case as compared to $d$-wave pairing, which is the favored gap symmetry for vanishing SOC. △ Less

Submitted 3 September, 2019; originally announced September 2019.

Comments: 19+15 pages, 21+16 figures

Report number: NBI CMT 2019

Pairing in the Two-Dimensional Hubbard Model from Weak to Strong Coupling

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

Abstract: The Hubbard model is the simplest model that is believed to exhibit superconductivity arising from purely repulsive interactions, and has been extensively applied to explore a variety of unconventional superconducting systems. Here we study the evolution of the leading superconducting instabilities of the single-orbital Hubbard model on a two-dimensional square lattice as a function of onsite Coul… ▽ More The Hubbard model is the simplest model that is believed to exhibit superconductivity arising from purely repulsive interactions, and has been extensively applied to explore a variety of unconventional superconducting systems. Here we study the evolution of the leading superconducting instabilities of the single-orbital Hubbard model on a two-dimensional square lattice as a function of onsite Coulomb repulsion $U$ and band filling by calculating the irreducible particle-particle scattering vertex obtained from dynamical cluster approximation (DCA) calculations, and compare the results to both perturbative Kohn-Luttinger (KL) theory as well as the widely used random phase approximation (RPA) spin-fluctuation pairing scheme. Near half-filling we find remarkable agreement of the hierarchy of the leading pairing states between these three methods, implying adiabatic continuity between weak- and strong-coupling pairing solutions of the Hubbard model. The $d_{x^2-y^2}$-wave instability is robust to increasing $U$ near half-filling as expected. Away from half filling, the predictions of KL and RPA at small $U$ for transitions to other pair states agree with DCA at intermediate $U$ as well as recent diagrammatic Monte Carlo calculations. RPA results fail only in the very dilute limit, where it yields a $d_{xy}$ ground state instead of a $p$-wave state established by diagrammatic Monte Carlo and low-order perturbative methods, as well as our DCA calculations. We discuss the origins of this discrepancy, highlighting the crucial role of the vertex corrections neglected in the RPA approach. Overall, comparison of the various methods over the entire phase diagram strongly suggests a smooth crossover of the superconducting interaction generated by local Hubbard interactions between weak and strong coupling. △ Less

Submitted 2 September, 2019; originally announced September 2019.

Comments: 9 pages, 5 figures

Report number: CMT NBI 2019

Journal ref: Phys. Rev. Research 2, 013108 (2020)

OGLE-2015-BLG-1649Lb: A gas giant planet around a low-mass dwarf

Authors: Masayuki Nagakane, Chien-Hsiu Lee, Naoki Koshimoto, Daisuke Suzuki, Andrzej Udalski, Jean-Philippe Beaulieu, Takahiro Sumi, David Bennett, Ian A. Bond, Nicholas J. Rattenbury, Etienne Bachelet, Martin Dominik, Fumio Abe, Richard Barry, Aparna Bhattacharya, Martin Donachie, H. Fujii, Akihiko Fukui, Yuki Hirao, Yoshitaka Itow, Y. Kamei, Iona Kondo, Man Cheung Alex Li, Y. Matsubara, Taro Matsuo , et al. (44 additional authors not shown)

Abstract: We report the discovery of an exoplanet in microlensing event OGLE-2015-BLG-1649. The planet/host-star mass ratio is $q =7.2 \times 10^{-3}$ and the projected separation normalized by the Einstein radius is $s = 0.9$. The upper limit of the lens flux is obtained from adaptive optics observations by IRCS/Subaru, which excludes the probability of a G-dwarf or more massive host star and helps to put… ▽ More We report the discovery of an exoplanet in microlensing event OGLE-2015-BLG-1649. The planet/host-star mass ratio is $q =7.2 \times 10^{-3}$ and the projected separation normalized by the Einstein radius is $s = 0.9$. The upper limit of the lens flux is obtained from adaptive optics observations by IRCS/Subaru, which excludes the probability of a G-dwarf or more massive host star and helps to put a tighter constraint on the lens mass as well as commenting on the formation scenarios of giant planets orbiting low-mass stars. We conduct a Bayesian analysis including constraints on the lens flux to derive the probability distribution of the physical parameters of the lens system. We thereby find that the masses of the host star and planet are $M_{L} = 0.34 \pm 0.19 M_{\odot}$ and $M_{p} = 2.5^{+1.5}_{-1.4} M_{Jup}$, respectively. The distance to the system is $D_{L} = 4.23^{+1.51}_{-1.64}$kpc. The projected star-planet separation is $a_{\perp} = 2.07^{+0.65}_{-0.77}$AU. The lens-source relative proper motion of the event is quite high, at $\sim 7.1 \, {\rm mas/yr}$. Therefore, we may be able to determine the lens physical parameters uniquely or place much stronger constraints on them by measuring the color-dependent centroid shift and/or the image elongation with additional high resolution imaging already a few years from now. △ Less

Submitted 2 October, 2019; v1 submitted 25 July, 2019; originally announced July 2019.

Comments: 24 pages, 5 figures, 3 tables. Accepted for publication in AJ