Showing 1–50 of 220 results for author: Madsen, K

Calculating the energy profile of an enzymatic reaction on a quantum computer

Authors: Patrick Ettenhuber, Mads Bøttger Hansen, Irfansha Shaik, Stig Elkjær Rasmussen, Pier Paolo Poier, Niels Kristian Madsen, Marco Majland, Frank Jensen, Lars Olsen, Nikolaj Thomas Zinner

Abstract: Quantum computing (QC) provides a promising avenue toward enabling quantum chemistry calculations, which are classically impossible due to a computational complexity that increases exponentially with system size. As fully fault-tolerant algorithms and hardware, for which an exponential speedup is predicted, are currently out of reach, recent research efforts are dedicated to developing and scaling… ▽ More Quantum computing (QC) provides a promising avenue toward enabling quantum chemistry calculations, which are classically impossible due to a computational complexity that increases exponentially with system size. As fully fault-tolerant algorithms and hardware, for which an exponential speedup is predicted, are currently out of reach, recent research efforts are dedicated to developing and scaling algorithms for Noisy Intermediate-Scale Quantum (NISQ) devices to showcase the practical utility of such machines. To demonstrate the utility of NISQ devices in the field of chemistry, we apply our recently developed FAST-VQE algorithm and a novel quantum gate reduction strategy based on propositional satisfiability together with standard optimization tools for the simulation of the rate-determining proton transfer step for CO2 hydration catalysed by carbonic anhydrase resulting in the first application of a quantum computing device for the simulation of an enzymatic reaction. To this end, we have combined classical force field simulations with quantum mechanical methods on classical and quantum computers in a hybrid calculation approach. The presented technique significantly enhances the accuracy and capabilities of QC-based molecular modeling and finally pushes it into compelling and realistic applications. The framework is general and can be applied beyond the case of computational enzymology. △ Less

Submitted 20 August, 2024; originally announced August 2024.

Neural-network-enabled molecular dynamics study of HfO$_2$ phase transitions

Authors: Sebastian Bichelmaier, Jesús Carrete, Georg K. H. Madsen

Abstract: The advances of machine-learned force fields have opened up molecular dynamics (MD) simulations for compounds for which ab-initio MD is too resource-intensive and phenomena for which classical force fields are insufficient. Here we describe a neural-network force field parametrized to reproduce the r2SCAN potential energy landscape of HfO$_2$. Based on an automatic differentiable implementation of… ▽ More The advances of machine-learned force fields have opened up molecular dynamics (MD) simulations for compounds for which ab-initio MD is too resource-intensive and phenomena for which classical force fields are insufficient. Here we describe a neural-network force field parametrized to reproduce the r2SCAN potential energy landscape of HfO$_2$. Based on an automatic differentiable implementation of the isothermal-isobaric (NPT) ensemble with flexible cell fluctuations, we study the phase space of HfO$_2$. We find excellent predictive capabilities regarding the lattice constants and experimental X-ray diffraction data. The phase transition away from monoclinic is clearly visible at a temperature around 2000 K, in agreement with available experimental data and previous calculations. Another abrupt change in lattice constants occurs around 3000 K. While the resulting lattice constants are closer to cubic, they exhibit a small tetragonal distortion, and there is no associated change in volume. We show that this high-temperature structure is in agreement with the available high-temperature diffraction data. △ Less

Submitted 5 August, 2024; originally announced August 2024.

International Astrophysical Consortium for High-energy Calibration: Summary of the 15th IACHEC Workshop

Authors: K. K. Madsen, V. Burwitz, K. Forster, C. E. Grant, M. Guainazzi, V. Kashyap, H. L. Marshall, E. D. Miller, L. Natalucci, P. P. Plucinsky, Y. Terada

Abstract: In this report, we summarize the activities of the International Astronomical Consortium for High Energy Calibration (IACHEC) from the 15th IACHEC Workshop in Pelham, Germany. Sixty scientists directly involved in the calibration of operational and future high-energy missions gathered for 3.5 days to discuss the status of the cross-calibration between the current international complement of X-ray… ▽ More In this report, we summarize the activities of the International Astronomical Consortium for High Energy Calibration (IACHEC) from the 15th IACHEC Workshop in Pelham, Germany. Sixty scientists directly involved in the calibration of operational and future high-energy missions gathered for 3.5 days to discuss the status of the cross-calibration between the current international complement of X-ray observatories and the possibilities to improve it. This summary consists of reports from the Working Groups with topics ranging across the identification and characterization of standard calibration sources, multi-observatory cross-calibration campaigns, appropriate and new statistical techniques, calibration of instruments and characterization of background, preservation of knowledge, and results for the benefit of the astronomical community. △ Less

Submitted 12 July, 2024; originally announced July 2024.

Comments: 10 pages, 1 figure. arXiv admin note: text overlap with arXiv:2111.01613

Digging its own Site: Linear Coordination Stabilizes a Pt1/Fe2O3 Single-Atom Catalyst

Authors: Ali Rafsanjani-Abbasi, Florian Buchner, Faith J. Lewis, Lena Puntscher, Florian Kraushofer, Panukorn Sombut, Moritz Eder, Jiri Pavelec, Erik Rheinfrank, Giada Franceschi, Viktor C. Birschitzky, Michele Riva, Cesare Franchini, Michael Schmid, Ulrike Diebold, Matthias Meier, Georg K. H. Madsen, Gareth S. Parkinson

Abstract: Determining the local coordination of the active site is a pre-requisite for the reliable modeling of single-atom catalysts (SACs). Obtaining such information is difficult on powder-based systems, so much emphasis is placed on density functional theory-based computations based on idealized low-index surfaces of the support. In this work, we investigate how Pt atoms bind to the (1-102) facet of Fe2… ▽ More Determining the local coordination of the active site is a pre-requisite for the reliable modeling of single-atom catalysts (SACs). Obtaining such information is difficult on powder-based systems, so much emphasis is placed on density functional theory-based computations based on idealized low-index surfaces of the support. In this work, we investigate how Pt atoms bind to the (1-102) facet of Fe2O3, a common support material in SAC. Using a combination of scanning tunneling microscopy (STM), x-ray photoelectron spectroscopy (XPS), and an extensive computational evolutionary search, we find that Pt atoms significantly reconfigure the support lattice to facilitate a pseudo-linear coordination to surface oxygen atoms. Despite breaking three surface Fe-O bonds, this geometry is favored by 0.84 eV over the best configuration involving an unperturbed support. We suggest that the linear O-Pt-O configuration is common in reactive Pt-based SAC systems because it balances thermal stability with the ability to adsorb reactants from the gas phase, and that extensive structural searches are likely necessary to determine realistic active site geometry in single-atom catalysis. △ Less

Submitted 26 June, 2024; originally announced June 2024.

A neural-network-backed effective harmonic potential study of the ambient pressure phases of hafnia

Authors: Sebastian Bichelmaier, Jesús Carrete, Ralf Wanzenböck, Florian Buchner, Georg K. H. Madsen

Abstract: Phonon-based approaches and molecular dynamics are widely established methods for gaining access to a temperature-dependent description of material properties. However, when a compound's phase space is vast, density-functional-theory-backed studies quickly reach prohibitive levels of computational expense. Here, we explore the complex phase structure of HfO2 using effective harmonic potentials bas… ▽ More Phonon-based approaches and molecular dynamics are widely established methods for gaining access to a temperature-dependent description of material properties. However, when a compound's phase space is vast, density-functional-theory-backed studies quickly reach prohibitive levels of computational expense. Here, we explore the complex phase structure of HfO2 using effective harmonic potentials based on a neural-network force field (NNFF) as a surrogate model. We detail the data acquisition and training strategy that enable the NNFF to provide almost ab-initio accuracy at a significantly reduced cost and present a recipe for automation. We demonstrate how the NNFF can generalize beyond its training data and that it is transferable between several phases of hafnia. We find that the thermal expansion of the low-symmetry phases agrees well with experimental results and we determine the P-43m phase to be the favorable (stoichiometric) cubic phase over the established Fm-3m. In contrast, the experimental lattice constants of the cubic phases are substantially larger than what is calculated for the corresponding stoichiometric phases. Furthermore, we show that the stoichiometric cubic phases are unlikely to be thermodynamically stable compared to the tetragonal and monoclinic phases, and hypothesize that they only exist in defect-stabilized forms. △ Less

Submitted 15 June, 2024; originally announced June 2024.

End-to-End Delivery in LEO Mega-constellations and the Reordering Problem

Authors: Rasmus Sibbern Frederiksen, Thomas Gundgaard Mulvad, Israel Leyva-Mayorga, Tatiana Kozlova Madsen, Federico Chiariotti

Abstract: Low Earth orbit (LEO) satellite mega-constellations with hundreds or thousands of satellites and inter-satellite links (ISLs) have the potential to provide global end-to-end connectivity. Furthermore, if the physical distance between source and destination is sufficiently long, end-to-end routing over the LEO constellation can provide lower latency when compared to the terrestrial infrastructure d… ▽ More Low Earth orbit (LEO) satellite mega-constellations with hundreds or thousands of satellites and inter-satellite links (ISLs) have the potential to provide global end-to-end connectivity. Furthermore, if the physical distance between source and destination is sufficiently long, end-to-end routing over the LEO constellation can provide lower latency when compared to the terrestrial infrastructure due to the faster propagation of electromagnetic waves in space than in optic fiber. However, the frequent route changes due to the movement of the satellites result in the out-of-order delivery of packets, causing sudden changes to the Round-Trip Time (RTT) that can be misinterpreted as congestion by congestion control algorithms. In this paper, the performance of three widely used congestion control algorithms, Cubic, Reno, and BBR, is evaluated in an emulated LEO satellite constellation with Free-Space Optical (FSO) ISLs. Furthermore, we perform a sensitivity analysis for Cubic by changing the satellite constellation parameters, length of the routes, and the positions of the source and destination to identify problematic routing scenarios. The results show that route changes can have profound transient effects on the goodput of the connection, posing problems for typical broadband applications. △ Less

Submitted 13 May, 2024; originally announced May 2024.

Comments: Submitted to IEEE PIMRC Workshops 2024

SYNCS: Synthetic Data and Contrastive Self-Supervised Training for Central Sulcus Segmentation

Authors: Vladyslav Zalevskyi, Kristoffer Hougaard Madsen

Abstract: Bipolar disorder (BD) and schizophrenia (SZ) are severe mental disorders with profound societal impact. Identifying risk markers early is crucial for understanding disease progression and enabling preventive measures. The Danish High Risk and Resilience Study (VIA) focuses on understanding early disease processes, particularly in children with familial high risk (FHR). Understanding structural bra… ▽ More Bipolar disorder (BD) and schizophrenia (SZ) are severe mental disorders with profound societal impact. Identifying risk markers early is crucial for understanding disease progression and enabling preventive measures. The Danish High Risk and Resilience Study (VIA) focuses on understanding early disease processes, particularly in children with familial high risk (FHR). Understanding structural brain changes associated with these diseases during early stages is essential for effective interventions. The central sulcus (CS) is a prominent brain landmark related to brain regions involved in motor and sensory processing. Analyzing CS morphology can provide valuable insights into neurodevelopmental abnormalities in the FHR group. However, segmenting the central sulcus (CS) presents challenges due to its variability, especially in adolescents. This study introduces two novel approaches to improve CS segmentation: synthetic data generation to model CS variability and self-supervised pre-training with multi-task learning to adapt models to new cohorts. These methods aim to enhance segmentation performance across diverse populations, eliminating the need for extensive preprocessing. △ Less

Submitted 22 March, 2024; originally announced March 2024.

The High Energy X-ray Probe (HEX-P): Instrument and Mission Profile

Authors: Kristin K. Madsen, Javier A. García, Daniel Stern, Rashied Armini, Stefano Basso, Diogo Coutinho, Brian W. Grefenstette, Steven Kenyon, Alberto Moretti, Patrick Morrisey, Kirpal Nandra, Giovanni Pareschi, Peter Predehl, Arne Rau, Daniele Spiga, Jörn Willms, William W. Zhang

Abstract: The High Energy X-ray Probe is a proposed NASA probe-class mission that combines the power of high angular resolution with a broad X-ray bandpass to provide the necessary leap in capabilities to address the important astrophysical questions of the next decade. HEX-P achieves breakthrough performance by combining technologies developed by experienced international partners. HEX-P will be launched i… ▽ More The High Energy X-ray Probe is a proposed NASA probe-class mission that combines the power of high angular resolution with a broad X-ray bandpass to provide the necessary leap in capabilities to address the important astrophysical questions of the next decade. HEX-P achieves breakthrough performance by combining technologies developed by experienced international partners. HEX-P will be launched into L1 to enable high observing efficiency. To meet the science goals, the payload consists of a suite of co-aligned X-ray telescopes designed to cover the 0.2 - 80 keV bandpass. The High Energy Telescope (HET) has an effective bandpass of 2 - 80 keV, and the Low Energy Telescope (LET) has an effective bandpass of 0.2 - 20 keV. The combination of bandpass and high observing efficiency delivers a powerful platform for broad science to serve a wide community. The baseline mission is five years, with 30% of the observing time dedicated to the PI-led program and 70% to a General Observer (GO) program. The GO program will be executed along with the PI-led program. △ Less

Submitted 7 December, 2023; originally announced December 2023.

Comments: 10 pages, 8 figures, to be submitted to Frontiers

The High Energy X-ray Probe (HEX-P): The Future of Hard X-ray Dual AGN Science

Authors: Ryan W. Pfeifle, Peter G. Boorman, Kimberly A. Weaver, Johannes Buchner, Francesca Civano, Kristin Madsen, Daniel Stern, Núria Torres-Albà, Emanuele Nardini, Claudio Ricci, Stefano Marchesi, D. R. Ballantyne, Dominic Sicilian, Chien-Ting Chen, Elias Kammoun, Ryan C. Hickox, Javier A. García, Labani Mallick

Abstract: A fundamental goal of modern-day astrophysics is to understand the connection between supermassive black hole (SMBH) growth and galaxy evolution. Merging galaxies offer one of the most dramatic channels for galaxy evolution known, capable of driving inflows of gas into galactic nuclei, potentially fueling both star formation and central SMBH activity. Dual active galactic nuclei (dual AGNs) in lat… ▽ More A fundamental goal of modern-day astrophysics is to understand the connection between supermassive black hole (SMBH) growth and galaxy evolution. Merging galaxies offer one of the most dramatic channels for galaxy evolution known, capable of driving inflows of gas into galactic nuclei, potentially fueling both star formation and central SMBH activity. Dual active galactic nuclei (dual AGNs) in late-stage mergers with nuclear pair separations $<10$ kpc are thus ideal candidates to study SMBH growth along the merger sequence since they coincide with the most transformative period for galaxies. However, dual AGNs can be extremely difficult to confirm and study. Hard X-ray ($>10$ keV) studies offer a relatively contamination-free tool for probing the dense obscuring environments predicted to surround the majority of dual AGN in late-stage mergers. To date, only a handful of the brightest and closest systems have been studied at these energies due to the demanding instrumental requirements involved. We demonstrate the unique capabilities of HEX-P to spatially resolve the soft and - for the first time - hard X-ray counterparts of closely-separated ($\sim2''-5''$) dual AGNs in the local Universe. By incorporating state-of-the-art physical torus models, we reproduce realistic broadband X-ray spectra expected for deeply embedded accreting SMBHs. Hard X-ray spatially resolved observations of dual AGNs - accessible only to HEX-P - will hence transform our understanding of dual AGN in the nearby Universe. △ Less

Submitted 9 November, 2023; originally announced November 2023.

Comments: Under review for Frontiers in Astronomy and Space Sciences

The High Energy X-ray Probe (HEX-P): Supernova remnants, pulsar wind nebulae, and nuclear astrophysics

Authors: Stephen Reynolds, Hongjun An, Moaz Abdelmaguid, Jason Alford, Chris L. Fryer, Kaya Mori, Melania Nynka, Jaegeun Park, Yukikatsu Terada, Jooyun Woo, Aya Bamba, Priyadarshini Bangale, Rebecca Diesing, Jordan Eagle, Stefano Gabici, Joseph Gelfand, Brian Grefenstette, Javier Garcia, Chanho Kim, Sajan Kumar, Brydyn Mac Intyre, Kristin Madsen, Silvia Manconi, Yugo Motogami, Hayato Ohsumi , et al. (7 additional authors not shown)

Abstract: HEX-P is a probe-class mission concept that will combine high spatial resolution X-ray imaging ($<10"$ full width at half maximum) and broad spectral coverage (0.2--80 keV) with an effective area far superior to current facilities (including XMM-Newton and NuSTAR) to enable revolutionary new insights into a variety of important astrophysical problems. HEX-P is ideally suited to address important p… ▽ More HEX-P is a probe-class mission concept that will combine high spatial resolution X-ray imaging ($<10"$ full width at half maximum) and broad spectral coverage (0.2--80 keV) with an effective area far superior to current facilities (including XMM-Newton and NuSTAR) to enable revolutionary new insights into a variety of important astrophysical problems. HEX-P is ideally suited to address important problems in the physics and astrophysics of supernova remnants (SNRs) and pulsar-wind nebulae (PWNe). For shell SNRs, HEX-P can greatly improve our understanding via more accurate spectral characterization and localization of non-thermal X-ray emission from both non-thermal-dominated SNRs and those containing both thermal and non-thermal components, and can discover previously unknown non-thermal components in SNRs. Multi-epoch HEX-P observations of several young SNRs (e.g., Cas A and Tycho) are expected to detect year-scale variabilities of X-ray filaments and knots, thus enabling us to determine fundamental parameters related to diffusive shock acceleration, such as local magnetic field strengths and maximum electron energies. For PWNe, HEX-P will provide spatially-resolved, broadband X-ray spectral data separately from their pulsar emission, allowing us to study how particle acceleration, cooling, and propagation operate in different evolution stages of PWNe. HEX-P is also poised to make unique and significant contributions to nuclear astrophysics of Galactic radioactive sources by improving detections of, or limits on, $^{44}$Ti in the youngest SNRs and by potentially discovering rare nuclear lines as evidence of double neutron star mergers. Throughout the paper, we present simulations of each class of objects, demonstrating the power of both the imaging and spectral capabilities of HEX-P to advance our knowledge of SNRs, PWNe, and nuclear astrophysics. △ Less

Submitted 8 November, 2023; originally announced November 2023.

Comments: 34 pages, 15 figures; part of a suite of papers describing the HEX-P hard X-ray mission concept

The High Energy X-ray Probe (HEX-P): The Circum-nuclear Environment of Growing Supermassive Black Holes

Authors: P. G. Boorman, N. Torres-Albà, A. Annuar, S. Marchesi, R. Pfeifle, D. Stern, F. Civano, M. Baloković, J. Buchner, C. Ricci, D. M. Alexander, W. N. Brandt, M. Brightman, C. T. Chen, S. Creech, P. Gandhi, J. A. García, F. Harrison, R. Hickox, E. Kammoun, S. LaMassa, G. Lanzuisi, L. Marcotulli, K. Madsen, G. Matt , et al. (10 additional authors not shown)

Abstract: Ever since the discovery of the first Active Galactic Nuclei (AGN), substantial observational and theoretical effort has been invested into understanding how massive black holes have evolved across cosmic time. Circum-nuclear obscuration is now established as a crucial component, with almost every AGN observed known to display signatures of some level of obscuration in their X-ray spectra. But des… ▽ More Ever since the discovery of the first Active Galactic Nuclei (AGN), substantial observational and theoretical effort has been invested into understanding how massive black holes have evolved across cosmic time. Circum-nuclear obscuration is now established as a crucial component, with almost every AGN observed known to display signatures of some level of obscuration in their X-ray spectra. But despite more than six decades of effort, substantial open questions remain: How does the accretion power impact the structure of the circum-nuclear obscurer? What are the dynamical properties of the obscurer? Can dense circum-nuclear obscuration exist around intrinsically weak AGN? How many intermediate mass black holes occupy the centers of dwarf galaxies? In this paper, we showcase a number of next-generation prospects attainable with the High Energy X-ray Probe (https://hexp.org) to contribute towards solving these questions in the 2030s. The uniquely broad (0.2--80 keV) and strictly simultaneous X-ray passband of HEX-P makes it ideally suited for studying the temporal co-evolution between the central engine and circum-nuclear obscurer. Improved sensitivities and reduced background will enable the development of spectroscopic models complemented by current and future multi-wavelength observations. We show that the angular resolution of HEX-P both below and above 10 keV will enable the discovery and confirmation of accreting massive black holes at both low accretion power and low black hole masses even when concealed by thick obscuration. In combination with other next-generation observations of the dusty hearts of nearby galaxies, HEX-P will hence be pivotal in paving the way towards a complete picture of black hole growth and galaxy co-evolution. △ Less

Submitted 8 November, 2023; originally announced November 2023.

Comments: To be submitted for peer review to a HEX-P Collection in Frontiers of Astronomy and Space Science. Constructive comments welcome!

The High Energy X-ray Probe (HEX-P): Sensitive broadband X-ray observations of transient phenomena in the 2030s

Authors: Murray Brightman, Raffaella Margutti, Ava Polzin, Amruta Jaodand, Kenta Hotokezaka, Jason A. J. Alford, Gregg Hallinan, Elias Kammoun, Kunal Mooley, Megan Masterson, Lea Marcotulli, Arne Rau, George A. Younes, Daniel Stern, Javier A. García, Kristin Madsen

Abstract: HEX-P will launch at a time when the sky is being routinely scanned for transient gravitational wave, electromagnetic and neutrino phenomena that will require the capabilities of a sensitive, broadband X-ray telescope for follow up studies. These include the merger of compact objects such as neutron stars and black holes, stellar explosions, and the birth of new compact objects. \hexp\ will probe… ▽ More HEX-P will launch at a time when the sky is being routinely scanned for transient gravitational wave, electromagnetic and neutrino phenomena that will require the capabilities of a sensitive, broadband X-ray telescope for follow up studies. These include the merger of compact objects such as neutron stars and black holes, stellar explosions, and the birth of new compact objects. \hexp\ will probe the accretion and ejecta from these transient phenomena through the study of relativistic outflows and reprocessed emission, provide unique capabilities for understanding jet physics, and potentially revealing the nature of the central engine. △ Less

Submitted 8 November, 2023; originally announced November 2023.

Comments: Submitted to Frontiers in Astronomy and Space Sciences

The High Energy X-ray Probe (HEX-P): resolving the nature of Sgr A* flares, compact object binaries and diffuse X-ray emission in the Galactic Center and beyond

Authors: Kaya Mori, Gabriele Ponti, Matteo Bachetti, Arash Bodaghee, Jonathan Grindlay, Jaesub Hong, Roman Krivonos, Ekaterina Kuznetsova, Shifra Mandel, Antonio Rodriguez, Giovanni Stel, Shuo Zhang, Tong Bao, Franz Bauer, Maica Clavel, Benjamin Coughenour, Javier A. Garcia, Julian Gerber, Brian Grefenstette, Amruta Jaodand, Bret Lehmer, Kristin Madsen, Melania Nynka, Peter Predehl, Ciro Salcedo , et al. (2 additional authors not shown)

Abstract: HEX-P is a probe-class mission concept that will combine high spatial resolution X-ray imaging ($<10"$ FWHM) and broad spectral coverage (0.2-80 keV) with an effective area far superior to current facilities' (including XMM-Newton and NuSTAR). These capabilities will enable revolutionary new insights into a variety of important astrophysical problems. We present scientific objectives and simulatio… ▽ More HEX-P is a probe-class mission concept that will combine high spatial resolution X-ray imaging ($<10"$ FWHM) and broad spectral coverage (0.2-80 keV) with an effective area far superior to current facilities' (including XMM-Newton and NuSTAR). These capabilities will enable revolutionary new insights into a variety of important astrophysical problems. We present scientific objectives and simulations of HEX-P observations of the Galactic Center (GC) and Bulge. We demonstrate the unique and powerful capabilities of the HEX-P observatory for studying both X-ray point sources and diffuse X-ray emission. HEX-P will be uniquely equipped to explore a variety of major topics in Galactic astrophysics, allowing us to (1) investigate broad-band properties of X-ray flares emitted from the supermassive black hole (BH) at Sgr A* and probe the associated particle acceleration and emission mechanisms; (2) identify hard X-ray sources detected by NuSTAR and determine X-ray point source populations in different regions and luminosity ranges; (3) determine the distribution of compact object binaries in the nuclear star cluster and the composition of the Galactic Ridge X-ray emission; (4) identify X-ray transients and measure fundamental parameters such as BH spin; (5) find hidden pulsars in the GC; (6) search for BH-OB binaries and hard X-ray flares from young stellar objects in young massive clusters; (7) measure white dwarf (WD) masses of magnetic CVs to deepen our understanding of CV evolution and the origin of WD magnetic fields; (8) explore primary particle accelerators in the GC in synergy with future TeV and neutrino observatories; (9) map out cosmic-ray distributions by observing non-thermal X-ray filaments; (10) explore past X-ray outbursts from Sgr A* through X-ray reflection components from giant molecular clouds. △ Less

Submitted 8 November, 2023; originally announced November 2023.

The High Energy X-ray Probe (HEX-P): Galactic PeVatrons, star clusters, superbubbles, microquasar jets, and gamma-ray binaries

Authors: Kaya Mori, Stephen Reynolds, Hongjun An, Aya Bamba, Roman Krivonos, Naomi Tsuji, Moaz Abdelmaguid, Jason Alford, Priyadarshini Bangale, Silvia Celli, Rebecca Diesing, Jordan Eagle, Chris L. Fryer, Stefano Gabici, Joseph Gelfand, Brian Grefenstette, Javier Garcia, Chanho Kim, Sajan Kumar, Ekaterina Kuznetsova, Brydyn Mac Intyre, Kristin Madsen, Silvia Manconi, Yugo Motogami, Hayato Ohsumi , et al. (10 additional authors not shown)

Abstract: HEX-P is a probe-class mission concept that will combine high spatial resolution X-ray imaging (<10" FWHM) and broad spectral coverage (0.2-80 keV) with an effective area far superior to current facilities (including XMM-Newton and NuSTAR) to enable revolutionary new insights into a variety of important astrophysical problems. With the recent discoveries of over 40 ultra-high-energy gamma-ray sour… ▽ More HEX-P is a probe-class mission concept that will combine high spatial resolution X-ray imaging (<10" FWHM) and broad spectral coverage (0.2-80 keV) with an effective area far superior to current facilities (including XMM-Newton and NuSTAR) to enable revolutionary new insights into a variety of important astrophysical problems. With the recent discoveries of over 40 ultra-high-energy gamma-ray sources (detected above 100 TeV) and neutrino emission in the Galactic Plane, we have entered a new era of multi-messenger astrophysics facing the exciting reality of Galactic PeVatrons. In the next decade, as more Galactic PeVatrons and TeV gamma-ray sources are expected to be discovered, the identification of their acceleration and emission mechanisms will be the most pressing issue in both particle and high-energy astrophysics. In this paper, along with its companion papers (Reynolds et al. 2023, Mori et al. 2023), we will present that HEX-P is uniquely suited to address important problems in various cosmic-ray accelerators, including Galactic PeVatrons, through investigating synchrotron X-ray emission of TeV-PeV electrons produced by both leptonic and hadronic processes. △ Less

Submitted 8 November, 2023; originally announced November 2023.

Comments: 31 pages, 12 figures, submitted to FrASS

The High Energy X-ray Probe (HEX-P): Bringing the Cosmic X-ray Background into Focus

Authors: Francesca Civano, Xiurui Zhao, Peter Boorman, Stefano Marchesi, Tonima Ananna, Samantha Creech, Chien-Ting Chen, Ryan Hickox, Daniel Stern, Kristin Madsen, Javier A. García, Ross Silver, James Aird, David M. Alexander, Mislav Balokovíc, William N. Brandt, Johannes Buchner, Poshak Gandhi, Elias Kammoun, Stephanie LaMassa, Giorgio Lanzuisi, Andrea Merloni, Alberto Moretti, Kirpal Nandra, Emanuele Nardini , et al. (7 additional authors not shown)

Abstract: Since the discovery of the cosmic X-ray background (CXB), astronomers have strived to understand the accreting supermassive black holes (SMBHs) contributing to its peak in the 10-40 keV band. Existing soft X-ray telescopes could study this population up to only 10 keV, and, while NuSTAR (focusing on 3--24 keV) made great progress, it also left significant uncertainties in characterizing the hard X… ▽ More Since the discovery of the cosmic X-ray background (CXB), astronomers have strived to understand the accreting supermassive black holes (SMBHs) contributing to its peak in the 10-40 keV band. Existing soft X-ray telescopes could study this population up to only 10 keV, and, while NuSTAR (focusing on 3--24 keV) made great progress, it also left significant uncertainties in characterizing the hard X-ray population, crucial for calibrating current population synthesis models. This paper presents an in-depth analysis of simulations of two extragalactic surveys (deep and wide) with the High-Energy X-ray Probe (HEX-P), each observed for 2 Ms. Applying established source detection techniques, we show that HEX-P surveys will reach a flux of $\sim$10$^{-15}$ erg s$^{-1}$ cm$^{-2}$ in the 10-40 keV band, an order of magnitude fainter than current NuSTAR surveys. With the large sample of new hard X-ray detected sources ($\sim2000$), we showcase HEX-P's ability to resolve more than 80% of the CXB up to 40 keV into individual sources. The expected precision of HEX-P's resolved background measurement will allow us to distinguish between population synthesis models of SMBH growth. HEX-P leverages accurate broadband (0.5-40 keV) spectral analysis and the combination of soft and hard X-ray colors to provide obscuration constraints even for the fainter sources, with the overall objective of measuring the Compton-thick fraction. With unprecedented sensitivity in the 10--40 keV band, HEX-P will explore the hard X-ray emission from AGN to flux limits never reached before, thus expanding the parameter space for serendipitous discoveries. Consequently, it is plausible that new models will be needed to capture the population HEX-P will unveil. △ Less

Submitted 16 November, 2023; v1 submitted 8 November, 2023; originally announced November 2023.

Comments: 28 pages, 13 Figures, 2 tables; submitted to Frontiers in Astronomy and Space Sciences; minor edits in the text, abstract and figure labels

The High Energy X-ray Probe (HEX-P): the most powerful jets through the lens of a superb X-ray eye

Authors: Lea Marcotulli, Marco Ajello, Markus Böttcher, Paolo Coppi, Luigi Costamante, Laura Di Gesu, Manel Errando, Javier A. García, Andrea Gokus, Ioannis Liodakis, Greg Madejski, Kristin Madsen, Alberto Moretti, Riccardo Middei, Felicia McBride, Maria Petropoulou, Bindu Rani, Tullia Sbarrato, Daniel Stern, Georgios Vasilopoulos, Michael Zacharias, Haocheng Zhang, the HEX-P Collaboration

Abstract: A fraction of the active supermassive black holes at the centers of galaxies in our Universe are capable of launching extreme kiloparsec-long relativistic jets. These jets are known multiband (radio to $γ$-ray) and multimessenger (neutrino) emitters, and some of them have been monitored over several decades at all accessible wavelengths. However, many open questions remain unanswered about the pro… ▽ More A fraction of the active supermassive black holes at the centers of galaxies in our Universe are capable of launching extreme kiloparsec-long relativistic jets. These jets are known multiband (radio to $γ$-ray) and multimessenger (neutrino) emitters, and some of them have been monitored over several decades at all accessible wavelengths. However, many open questions remain unanswered about the processes powering these highly energetic phenomena. These jets intrinsically produce soft-to-hard X-ray emission that extends from $E\sim0.1\,\rm keV$ up to $E>100\,\rm keV$. Simultaneous broadband X-ray coverage, combined with excellent timing and imaging capabilities, is required to uncover the physics of jets. Indeed, truly simultaneous soft-to-hard X-ray coverage, in synergy with current and upcoming high-energy facilities (such as IXPE, COSI, CTAO, etc.) and neutrino detectors (e.g., IceCube), would enable us to disentangle the particle population responsible for the high-energy radiation from these jets. A sensitive hard X-ray survey ($F_{8-24\,\rm keV}<10^{-15}\,\rm erg~cm^{-2}~s^{-1}$) could unveil the bulk of their population in the early Universe. Acceleration and radiative processes responsible for the majority of their X-ray emission would be pinned down by microsecond timing capabilities at both soft and hard X-rays. Furthermore, imaging jet structures for the first time in the hard X-ray regime could unravel the origin of their high-energy emission. The proposed Probe-class mission concept High Energy X-ray Probe (HEX-P) combines all these required capabilities, making it the crucial next-generation X-ray telescope in the multi-messenger, time-domain era. HEX-P will be the ideal mission to unravel the science behind the most powerful accelerators in the universe. △ Less

Submitted 8 November, 2023; originally announced November 2023.

Comments: 30 pages, 12 figures, under review from Frontiers in Astronomy and Space Sciences

The High Energy X-ray Probe (HEX-P): Probing Accretion onto Stellar Mass Black Holes

Authors: Riley Connors, John Tomsick, Paul Draghis, Benjamin Coughenour, Aarran Shaw, Javier Garcia, Dominic Walton, Kristin Madsen, Daniel Stern, Nicole Cavero Rodriguez, Thomas Dauser, Melania Del Santo, Jiachen Jiang, Henric Krawczynski, Honghui Liu, Joseph Neilsen, Michael Nowak, Sean Pike, Andrea Santangelo, Navin Sridhar, Andrew West, Joern Wilms, the HEX-P Team

Abstract: Accretion is a universal astrophysical process that plays a key role in cosmic history, from the epoch of reionization to galaxy and stellar formation and evolution. Accreting stellar-mass black holes in X-ray binaries are one of the best laboratories to study the accretion process and probe strong gravity -- and most importantly, to measure the angular momentum, or spin, of black holes, and its r… ▽ More Accretion is a universal astrophysical process that plays a key role in cosmic history, from the epoch of reionization to galaxy and stellar formation and evolution. Accreting stellar-mass black holes in X-ray binaries are one of the best laboratories to study the accretion process and probe strong gravity -- and most importantly, to measure the angular momentum, or spin, of black holes, and its role as a powering mechanism for relativistic astrophysical phenomena. Comprehensive characterization of the disk-corona system of accreting black holes, and their co-evolution, is fundamental to measurements of black hole spin. Here, we use simulated data to demonstrate how key unanswered questions in the study of accreting stellar-mass black holes will be addressed by the {\it High Energy X-ray Probe} (\hexp). \hexp\ is a probe-class mission concept that will combine high spatial resolution X-ray imaging and broad spectral coverage ($0.2\mbox{--}80$keV) with a sensitivity superior to current facilities (including \xmm\ and \nustar) to enable revolutionary new insights into a variety of important astrophysical problems. We illustrate the capability of \hexp\ to: 1) measure the evolving structures of black hole binary accretion flows down to low ($\lesssim0.1\%$) Eddington-scaled luminosities via detailed X-ray reflection spectroscopy; 2) provide unprecedented spectral observations of the coronal plasma, probing its elusive geometry and energetics; 3) perform detailed broadband studies of stellar mass black holes in nearby galaxies, thus expanding the repertoire of sources we can use to study accretion physics and determine the fundamental nature of black holes; and 4) act as a complementary observatory to a range of future ground and space-based astronomical observatories, thus providing key spectral measurements of the multi-component emission from the inner accretion flows of BH-XRBs. △ Less

Submitted 8 November, 2023; originally announced November 2023.

Comments: 31 pages, 7 figures

The High Energy X-ray Probe (HEX-P): Constraining Supermassive Black Hole Growth with Population Spin Measurements

Authors: J. M. Piotrowska, J. A. García, D. J. Walton, R. S. Beckmann, D. Stern, D. R. Ballantyne, D. R. Wilkins, S. Bianchi, P. G. Boorman, J. Buchner, C. -T. Chen, P. Coppi, T. Dauser, A. C. Fabian, E. Kammoun, K. Madsen, L. Mallick, G. Matt, G. Matzeu, E. Nardini, A. Pizzetti, S. Puccetti, C. Ricci, F. Tombesi, N. Torres-Albà , et al. (2 additional authors not shown)

Abstract: Constraining the primary growth channel of supermassive black holes (SMBH) remains one the most actively debated questions in the context of cosmological structure formation. Owing to the expected connection between SMBH spin parameter evolution and the accretion and merger history of individual black holes, population spin measurements offer a rare observational window into the SMBH cosmic growth… ▽ More Constraining the primary growth channel of supermassive black holes (SMBH) remains one the most actively debated questions in the context of cosmological structure formation. Owing to the expected connection between SMBH spin parameter evolution and the accretion and merger history of individual black holes, population spin measurements offer a rare observational window into the SMBH cosmic growth. As of today, the most common method for estimating SMBH spin relies on modeling the relativistically broaden atomic profiles in the reflection spectrum observed in X-rays. In this paper, we study the observational requirements needed to confidently distinguish between the primary SMBH growth channels, based on their distinct spin-mass distributions predicted by the Horizon-AGN cosmological simulation. In doing so, we characterize outstanding limitations associated with the existing measurements and discuss the landscape of future observational campaigns, which can be planned and executed with future X-ray observatories. We focus our attention on the High-Energy X-ray Probe (HEX-P), a concept probe-class mission aimed to serve the high-energy community in the 2030s. △ Less

Submitted 8 November, 2023; originally announced November 2023.

Comments: 26 pages, 11 figures, submitted for publication in Frontiers in Astronomy and Space Sciences

The High Energy X-ray Probe: Resolved X-ray Populations in Extragalactic Environments

Authors: Bret D. Lehmer, Kristen Garofali, Breanna A. Binder, Francesca Fornasini, Neven Vulic, Andreas Zezas, Ann Hornschemeier, Margaret Lazzarini, Hannah Moon, Toni Venters, Daniel Wik, Mihoko Yukita, Matteo Bachetti, Javier A. García, Brian Grefenstette, Kristin Madsen, Kaya Mori, Daniel Stern

Abstract: We construct simulated galaxy data sets based on the High Energy X-ray Probe (HEX-P) mission concept to demonstrate the significant advances in galaxy science that will be yielded by the HEX-P observatory. The combination of high spatial resolution imaging ($<$20 arcsec FWHM), broad spectral coverage (0.2-80 keV), and sensitivity superior to current facilities (e.g., XMM-Newton and NuSTAR) will en… ▽ More We construct simulated galaxy data sets based on the High Energy X-ray Probe (HEX-P) mission concept to demonstrate the significant advances in galaxy science that will be yielded by the HEX-P observatory. The combination of high spatial resolution imaging ($<$20 arcsec FWHM), broad spectral coverage (0.2-80 keV), and sensitivity superior to current facilities (e.g., XMM-Newton and NuSTAR) will enable HEX-P to detect hard (4-25 keV) X-ray emission from resolved point-source populations within $\sim$800 galaxies and integrated emission from $\sim$6000 galaxies out to 100 Mpc. These galaxies cover wide ranges of galaxy types (e.g., normal, starburst, and passive galaxies) and properties (e.g., metallicities and star-formation histories). In such galaxies, HEX-P will: (1) provide unique information about X-ray binary populations, including accretor demographics (black hole and neutron stars), distributions of accretion states and state transition cadences; (2) place order-of-magnitude more stringent constraints on inverse Compton emission associated with particle acceleration in starburst environments; and (3) put into clear context the contributions from X-ray emitting populations to both ionizing the surrounding interstellar medium in low-metallicity galaxies and heating the intergalactic medium in the $z > 8$ Universe. △ Less

Submitted 8 November, 2023; originally announced November 2023.

Comments: Accepted for publication in Frontiers in Astronomy and Space Sciences

The High Energy X-ray Probe (HEX-P): Studying Extreme Accretion with Ultraluminous X-ray Sources

Authors: Matteo Bachetti, Matthew J. Middleton, Ciro Pinto, Andrés Gúrpide, Dominic J. Walton, Murray Brightman, Bret Lehmer, Timothy P. Roberts, Georgios Vasilopoulos, Jason Alford, Roberta Amato, Elena Ambrosi, Lixin Dai, Hannah P. Earnshaw, Hamza El Byad, Javier A. García, Gian Luca Israel, Amruta Jaodand, Kristin Madsen, Chandreyee Maitra, Shifra Mandel, Kaya Mori, Fabio Pintore, Ken Ohsuga, Maura Pilia , et al. (3 additional authors not shown)

Abstract: Ultraluminous X-ray sources (ULXs) represent an extreme class of accreting compact objects: from the identification of some of the accretors as neutron stars to the detection of powerful winds travelling at 0.1-0.2 c, the increasing evidence points towards ULXs harbouring stellar-mass compact objects undergoing highly super-Eddington accretion. Measuring their intrinsic properties, such as the acc… ▽ More Ultraluminous X-ray sources (ULXs) represent an extreme class of accreting compact objects: from the identification of some of the accretors as neutron stars to the detection of powerful winds travelling at 0.1-0.2 c, the increasing evidence points towards ULXs harbouring stellar-mass compact objects undergoing highly super-Eddington accretion. Measuring their intrinsic properties, such as the accretion rate onto the compact object, the outflow rate, the masses of accretor/companion -- hence their progenitors, lifetimes, and future evolution -- is challenging due to ULXs being mostly extragalactic and in crowded fields. Yet ULXs represent our best opportunity to understand super-Eddington accretion physics and the paths through binary evolution to eventual double compact object binaries and gravitational wave sources. Through a combination of end-to-end and single-source simulations, we investigate the ability of HEX-P to study ULXs in the context of their host galaxies and compare it to XMM-Newton and NuSTAR, the current instruments with the most similar capabilities. HEX-P's higher sensitivity, which is driven by its narrow point-spread function and low background, allows it to detect pulsations and broad spectral features from ULXs better than XMM-Newton and NuSTAR. We describe the value of HEX-P in understanding ULXs and their associated key physics, through a combination of broadband sensitivity, timing resolution, and angular resolution, which make the mission ideal for pulsation detection and low-background, broadband spectral studies. △ Less

Submitted 10 November, 2023; v1 submitted 8 November, 2023; originally announced November 2023.

Comments: 24 pages, 9 figures. Accepted for publication in Frontiers of Astronomy and Space Science. Minor corrections included

The High Energy X-ray Probe (HEX-P): A New Window into Neutron Star Accretion

Authors: R. M. Ludlam, C. Malacaria, E. Sokolova-Lapa, F. Fuerst, P. Pradhan, A. W. Shaw, K. Pottschmidt, S. Pike, G. Vasilopoulos, J. Wilms, J. A. García, K. Madsen, D. Stern, C. Maitra, M. Del Santo, D. J. Walton, M. C. Brumback, J. van den Eijnden

Abstract: Accreting neutron stars (NSs) represent a unique laboratory for probing the physics of accretion in the presence of strong magnetic fields ($B\gtrsim 10^8$ G). Additionally, the matter inside the NS itself exists in an ultra-dense, cold state that cannot be reproduced in Earth-based laboratories. Hence, observational studies of these objects are a way to probe the most extreme physical regimes. He… ▽ More Accreting neutron stars (NSs) represent a unique laboratory for probing the physics of accretion in the presence of strong magnetic fields ($B\gtrsim 10^8$ G). Additionally, the matter inside the NS itself exists in an ultra-dense, cold state that cannot be reproduced in Earth-based laboratories. Hence, observational studies of these objects are a way to probe the most extreme physical regimes. Here we present an overview of the field and discuss the most important outstanding problems related to NS accretion. We show how these open questions regarding accreting NSs in both low-mass and high-mass X-ray binary systems can be addressed with the High-Energy X-ray Probe (HEX-P) via simulated data. In particular, with the broad X-ray passband and improved sensitivity afforded by a low X-ray background, HEX-P will be able to 1) distinguish between competing continuum emission models; 2) provide tighter upper limits on NS radii via reflection modeling techniques that are independent and complementary to other existing methods; 3) constrain magnetic field geometry, plasma parameters, and accretion column emission patterns by characterizing fundamental and harmonic cyclotron lines and exploring their behavior with pulse phase; 4) directly measure the surface magnetic field strength of highly magnetized NSs at the lowest accretion luminosities; as well as 5) detect cyclotron line features in extragalactic sources and probe their dependence on luminosity in the super-Eddington regime in order to distinguish between geometrical evolution and accretion-induced decay of the magnetic field. In these ways HEX-P will provide an essential new tool for exploring the physics of NSs, their magnetic fields, and the physics of extreme accretion. △ Less

Submitted 8 November, 2023; originally announced November 2023.

Comments: 37 pages, 13 figures, 4 tables, To be published in a special topical issue by Frontiers in Astronomy and Space Sciences on probe-class mission concept HEX-P

The High Energy X-ray Probe (HEX-P): Probing the physics of the X-ray corona in active galactic nuclei

Authors: E. Kammoun, A. M. Lohfink, M. Masterson, D. R. Wilkins, X. Zhao, M. Baloković, P. G. Boorman, R. M. T. Connors, P. Coppi, A. C. Fabian, J. A. García, K. K. Madsen, N. Rodriguez Cavero, N. Sridhar, D. Stern, J. Tomsick, T. Wevers, D. J. Walton, S. Bianchi, J. Buchner, F. Civano, G. Lanzuisi, L. Mallick, G. Matt, A. Merloni , et al. (6 additional authors not shown)

Abstract: The hard X-ray emission in active galactic nuclei (AGN) and black hole X-ray binaries is thought to be produced by a hot cloud of electrons referred to as the corona. This emission, commonly described by a power law with a high-energy cutoff, is suggestive of Comptonization by thermal electrons. While several hypotheses have been proposed to explain the origin, geometry, and composition of the cor… ▽ More The hard X-ray emission in active galactic nuclei (AGN) and black hole X-ray binaries is thought to be produced by a hot cloud of electrons referred to as the corona. This emission, commonly described by a power law with a high-energy cutoff, is suggestive of Comptonization by thermal electrons. While several hypotheses have been proposed to explain the origin, geometry, and composition of the corona, we still lack a clear understanding of this fundamental component. NuSTAR has been playing a key role improving our knowledge of X-ray coronae thanks to its unprecedented sensitivity above 10 keV. However, these constraints are limited to bright, nearby sources. The High Energy X-ray Probe (HEX-P) is a probe-class mission concept combining high spatial resolution X-ray imaging and broad spectral coverage (0.2-80 keV) with a sensitivity superior to current facilities. In this paper, we highlight the major role that HEX-P will play in further advancing our insights of X-ray coronae, notably in AGN. We demonstrate how HEX-P will measure key properties and track the temporal evolution of coronae in unobscured AGN. This will allow us to determine their electron distribution and test the dominant emission mechanisms. Furthermore, we show how HEX-P will accurately estimate the coronal properties of obscured AGN in the local Universe, helping address fundamental questions about AGN unification. In addition, HEX-P will characterize coronae in a large sample of luminous quasars at cosmological redshifts for the first time and track the evolution of coronae in transient systems in real time. We also demonstrate how HEX-P will enable estimating the coronal geometry using spectral-timing techniques. HEX-P will thus be essential to understand the evolution and growth of black holes over a broad range of mass, distance, and luminosity, and will help uncover the black holes' role in shaping the Universe. △ Less

Submitted 8 November, 2023; originally announced November 2023.

Comments: Submitted to Frontiers in Astronomy and Space Sciences after reviewers reports

Wavelet-based image decomposition method for NuSTAR stray light background studies

Authors: Andrey Mukhin, Roman Krivonos, Alexey Vikhlinin, Brian Grefenstette, Kristin Madsen, Daniel Wik

Abstract: The large side aperture of the NuSTAR telescope for unfocused photons (so-called stray light) is a known source of rich astrophysical information. To support many studies based on the NuSTAR stray light data, we present a fully automatic method for determining detector area suitable for background analysis and free from any kind of focused X-ray flux. The method's main idea is `a trous' wavelet im… ▽ More The large side aperture of the NuSTAR telescope for unfocused photons (so-called stray light) is a known source of rich astrophysical information. To support many studies based on the NuSTAR stray light data, we present a fully automatic method for determining detector area suitable for background analysis and free from any kind of focused X-ray flux. The method's main idea is `a trous' wavelet image decomposition, capable of detecting structures of any spatial scale and shape, which makes the method of general use. Applied to the NuSTAR data, the method provides a detector image region with the highest possible statistical quality, suitable for the NuSTAR stray light studies. We developed an open-source Python nuwavdet package, which implements the presented method. The package contains subroutines to generate detector image region for further stray light analysis and/or to produce a list of detector bad-flagged pixels for processing in the NuSTAR Data Analysis Software for conventional X-ray analysis. △ Less

Submitted 27 October, 2023; v1 submitted 16 October, 2023; originally announced October 2023.

Comments: 9 pages, 10 figures. Published in Journal of Astronomical Telescopes, Instruments, and Systems (JATIS)

Journal ref: Journal of Astronomical Telescopes, Instruments, and Systems 9(4), 048001 (26 October 2023)

Quantitative predictions of the thermal conductivity in transition metal dichalcogenides: The impact of point defects in MoS$_2$ and WS$_2$ monolayers

Authors: Srinivisan Mahendran, Jesús Carrete, Andreas Isacsson, Georg K. H. Madsen, Paul Erhart

Abstract: Transition metal dichalcogenides are investigated for various applications at the nanoscale thanks to their unique combination of properties and dimensionality. For many of the anticipated applications, heat conduction plays an important role. At the same time, these materials often contain relatively large amounts of point defects. Here, we provide a systematic analysis of the impact of intrinsic… ▽ More Transition metal dichalcogenides are investigated for various applications at the nanoscale thanks to their unique combination of properties and dimensionality. For many of the anticipated applications, heat conduction plays an important role. At the same time, these materials often contain relatively large amounts of point defects. Here, we provide a systematic analysis of the impact of intrinsic and selected extrinsic defects on the lattice thermal conductivity of MoS$_2$ and WS$_2$ monolayers. We combine Boltzmann transport theory and the Green's function-based T-matrix approach for the calculation of scattering rates. The force constants for the defect configurations are obtained from density functional theory calculations via a regression approach, which allows us to sample a rather large number of defects at a moderate computational cost and to systematically enforce both the translational and rotational acoustic sum rules. The calculated lattice thermal conductivity is in quantitative agreement with experimental data for heat transport and defect concentrations for both MoS$_2$ and WS$_2$. Crucially, this demonstrates that the strong deviation from a 1/T-temperature dependence of the lattice thermal conductivity observed experimentally, can be fully explained by the presence of point defects. We furthermore predict the scattering strengths of the intrinsic defects to decrease in the sequence $V_{Mo}\approx V_{2S}^=>V_{2S}^\perp>V_S>S_{ad}$ in both materials, while the scattering rates for the extrinsic (adatom) defects decrease with increasing mass such that Li$_{ad}$>Na$_{ad}$>K$_{ad}$. Compared to earlier work, we find that both intrinsic and extrinsic adatoms are relatively weak scatterers. We attribute this difference to the treatment of the translational and rotational acoustic sum rules, which if not enforced can lead to spurious contributions in the zero-frequency limit. △ Less

Submitted 13 October, 2023; originally announced October 2023.

Comments: 8 pages, 4 figures

Journal ref: Journal of Physical Chemistry C 128, 1709 (2024)

Efficient time-dependent vibrational coupled cluster computations with time-dependent basis sets at the two-mode coupling level: full and hybrid TDMVCC[2]

Authors: Andreas Buchgraitz Jensen, Mads Greisen Højlund, Alberto Zoccante, Niels Kristian Madsen, Ove Christiansen

Abstract: The computation of the nuclear quantum dynamics of molecules is challenging, requiring both accuracy and efficiency to be applicable to systems of interest. Recently, theories have been developed for employing time-dependent basis functions (denoted modals) with vibrational coupled cluster theory (TDMVCC). The TDMVCC method was introduced along with a pilot implementation, which illustrated good a… ▽ More The computation of the nuclear quantum dynamics of molecules is challenging, requiring both accuracy and efficiency to be applicable to systems of interest. Recently, theories have been developed for employing time-dependent basis functions (denoted modals) with vibrational coupled cluster theory (TDMVCC). The TDMVCC method was introduced along with a pilot implementation, which illustrated good accuracy in benchmark computations. In this paper we report an efficient implementation of TDMVCC covering the case where the wave function and Hamiltonian contain up to two-mode couplings. After a careful regrouping of terms, the wave function can be propagated with a cubic computational scaling with respect to the number of degrees of freedom. We discuss the use of a restricted set of active one-mode basis functions for each mode, as well as two interesting limits: i) the use of a full active basis where the variational modal determination amounts essentially to the variational determination of a time-dependent reference state for the cluster expansion; and ii) the use of a single function as active basis for some degrees of freedom. The latter case defines a hybrid TDMVCC/TDH approach which can obtain an even lower computational scaling. The resulting computational scaling for hybrid and full TDMVCC[2] are illustrated for polyaromatic hydrocarbons (PAHs) with up to 264 modes. Finally, computations on the internal vibrational redistribution on benzoic acid (39 modes) are used to show the faster convergence of TDMVCC/TDH hybrid computations towards TDMVCC compared to simple neglect of some degrees of freedom. △ Less

Submitted 29 August, 2023; originally announced August 2023.

Comments: Will be submitted to Journal of Chemical Physics

Neural-Network Force Field Backed Nested Sampling: Study of the Silicon p-T Phase Diagram

Authors: N. Unglert, J. Carrete, L. B. Pártay, G. K. H. Madsen

Abstract: Nested sampling is a promising method for calculating phase diagrams of materials, however, the computational cost limits its applicability if ab-initio accuracy is required. In the present work, we report on the efficient use of a neural-network force field in conjunction with the nested-sampling algorithm. We train our force fields on a recently reported database of silicon structures and demons… ▽ More Nested sampling is a promising method for calculating phase diagrams of materials, however, the computational cost limits its applicability if ab-initio accuracy is required. In the present work, we report on the efficient use of a neural-network force field in conjunction with the nested-sampling algorithm. We train our force fields on a recently reported database of silicon structures and demonstrate our approach on the low-pressure region of the silicon pressure-temperature phase diagram between 0 and \SI{16}{GPa}. The simulated phase diagram shows a good agreement with experimental results, closely reproducing the melting line. Furthermore, all of the experimentally stable structures within the investigated pressure range are also observed in our simulations. We point out the importance of the choice of exchange-correlation functional for the training data and show how the meta-GGA r2SCAN plays a pivotal role in achieving accurate thermodynamic behaviour using nested-sampling. We furthermore perform a detailed analysis of the exploration of the potential energy surface and highlight the critical role of a diverse training data set. △ Less

Submitted 22 August, 2023; originally announced August 2023.

Electron-induced non-monotonic pressure dependence of the lattice thermal conductivity of θ-TaN

Authors: Ashis Kundu, Yani Chen, Xiaolong Yang, Fanchen Meng, Jesús Carrete, Mukul Kabir, Georg K. H. Madsen, Wu Li

Abstract: Recent theoretical and experimental research suggests that $θ$-TaN is a semimetal with high thermal conductivity ($κ$), primarily due to the contribution of phonons ($κ_\texttt{ph}$). By using first-principles calculations, we show a non-monotonic pressure dependence of the $κ$ of $θ$-TaN. $κ_\texttt{ph}$ first increases until it reaches a maximum at around 60~GPa, and then decreases. This anomalo… ▽ More Recent theoretical and experimental research suggests that $θ$-TaN is a semimetal with high thermal conductivity ($κ$), primarily due to the contribution of phonons ($κ_\texttt{ph}$). By using first-principles calculations, we show a non-monotonic pressure dependence of the $κ$ of $θ$-TaN. $κ_\texttt{ph}$ first increases until it reaches a maximum at around 60~GPa, and then decreases. This anomalous behaviour is a consequence of the competing pressure responses of phonon-phonon and phonon-electron interactions, in contrast to the known materials BAs and BP, where the non-monotonic pressure dependence is caused by the interplay between different phonon-phonon scattering channels. Although TaN has phonon dispersion features similar to BAs at ambient pressure, its response to pressure is different and an overall stiffening of the phonon branches takes place. Consequently, the relevant phonon-phonon scattering weakens as pressure increases. However, the increased electronic density of states near the Fermi level, and specifically the emergence of additional pockets of the Fermi surface at the high-symmetry L point in the Brillouin zone, leads to a substantial increase in phonon-electron scattering at high pressures, driving a decrease in $κ_{\mathrm{ph}}$. At intermediate pressures ($\sim$~20$-$70~GPa), the $κ$ of TaN surpasses that of BAs. Our work provides deeper insight into phonon transport in semimetals and metals where phonon-electron scattering is relevant. △ Less

Submitted 19 March, 2024; v1 submitted 28 July, 2023; originally announced July 2023.

Comments: 6 pages, 4 figures

Journal ref: Phys. Rev. Lett. 132, 116301 (2024)

How to verify the precision of density-functional-theory implementations via reproducible and universal workflows

Authors: Emanuele Bosoni, Louis Beal, Marnik Bercx, Peter Blaha, Stefan Blügel, Jens Bröder, Martin Callsen, Stefaan Cottenier, Augustin Degomme, Vladimir Dikan, Kristjan Eimre, Espen Flage-Larsen, Marco Fornari, Alberto Garcia, Luigi Genovese, Matteo Giantomassi, Sebastiaan P. Huber, Henning Janssen, Georg Kastlunger, Matthias Krack, Georg Kresse, Thomas D. Kühne, Kurt Lejaeghere, Georg K. H. Madsen, Martijn Marsman , et al. (20 additional authors not shown)

Abstract: In the past decades many density-functional theory methods and codes adopting periodic boundary conditions have been developed and are now extensively used in condensed matter physics and materials science research. Only in 2016, however, their precision (i.e., to which extent properties computed with different codes agree among each other) was systematically assessed on elemental crystals: a firs… ▽ More In the past decades many density-functional theory methods and codes adopting periodic boundary conditions have been developed and are now extensively used in condensed matter physics and materials science research. Only in 2016, however, their precision (i.e., to which extent properties computed with different codes agree among each other) was systematically assessed on elemental crystals: a first crucial step to evaluate the reliability of such computations. We discuss here general recommendations for verification studies aiming at further testing precision and transferability of density-functional-theory computational approaches and codes. We illustrate such recommendations using a greatly expanded protocol covering the whole periodic table from Z=1 to 96 and characterizing 10 prototypical cubic compounds for each element: 4 unaries and 6 oxides, spanning a wide range of coordination numbers and oxidation states. The primary outcome is a reference dataset of 960 equations of state cross-checked between two all-electron codes, then used to verify and improve nine pseudopotential-based approaches. Such effort is facilitated by deploying AiiDA common workflows that perform automatic input parameter selection, provide identical input/output interfaces across codes, and ensure full reproducibility. Finally, we discuss the extent to which the current results for total energies can be reused for different goals (e.g., obtaining formation energies). △ Less

Submitted 26 May, 2023; originally announced May 2023.

Comments: Main text: 23 pages, 4 figures. Supplementary: 68 pages. Nature Review Physics 2023

Journal ref: Nat. Rev. Phys. 6, 45 (2024)

Coronal Heating as Determined by the Solar Flare Frequency Distribution Obtained by Aggregating Case Studies

Authors: James Paul Mason, Alexandra Werth, Colin G. West, Allison A. Youngblood, Donald L. Woodraska, Courtney Peck, Kevin Lacjak, Florian G. Frick, Moutamen Gabir, Reema A. Alsinan, Thomas Jacobsen, Mohammad Alrubaie, Kayla M. Chizmar, Benjamin P. Lau, Lizbeth Montoya Dominguez, David Price, Dylan R. Butler, Connor J. Biron, Nikita Feoktistov, Kai Dewey, N. E. Loomis, Michal Bodzianowski, Connor Kuybus, Henry Dietrick, Aubrey M. Wolfe , et al. (977 additional authors not shown)

Abstract: Flare frequency distributions represent a key approach to addressing one of the largest problems in solar and stellar physics: determining the mechanism that counter-intuitively heats coronae to temperatures that are orders of magnitude hotter than the corresponding photospheres. It is widely accepted that the magnetic field is responsible for the heating, but there are two competing mechanisms th… ▽ More Flare frequency distributions represent a key approach to addressing one of the largest problems in solar and stellar physics: determining the mechanism that counter-intuitively heats coronae to temperatures that are orders of magnitude hotter than the corresponding photospheres. It is widely accepted that the magnetic field is responsible for the heating, but there are two competing mechanisms that could explain it: nanoflares or Alfvén waves. To date, neither can be directly observed. Nanoflares are, by definition, extremely small, but their aggregate energy release could represent a substantial heating mechanism, presuming they are sufficiently abundant. One way to test this presumption is via the flare frequency distribution, which describes how often flares of various energies occur. If the slope of the power law fitting the flare frequency distribution is above a critical threshold, $α=2$ as established in prior literature, then there should be a sufficient abundance of nanoflares to explain coronal heating. We performed $>$600 case studies of solar flares, made possible by an unprecedented number of data analysts via three semesters of an undergraduate physics laboratory course. This allowed us to include two crucial, but nontrivial, analysis methods: pre-flare baseline subtraction and computation of the flare energy, which requires determining flare start and stop times. We aggregated the results of these analyses into a statistical study to determine that $α= 1.63 \pm 0.03$. This is below the critical threshold, suggesting that Alfvén waves are an important driver of coronal heating. △ Less

Submitted 9 May, 2023; originally announced May 2023.

Comments: 1,002 authors, 14 pages, 4 figures, 3 tables, published by The Astrophysical Journal on 2023-05-09, volume 948, page 71

Measuring the Cosmic X-ray Background in 3-20keV with Straylight from NuSTAR

Authors: Steven Rossland, Daniel Wik, Brian Grefenstette, Nico Cappelluti, Francesca Civano, Fabio Gastaldello, Roberto Gilli, Fiona Harrison, Ann Hornschemeier, Ryan Hickox, Roman Krivonos, Kristin Madsen, Silvano Molendi, Andrew Ptak, Daniel Stern, Andreas Zoglauer

Abstract: By characterizing the contribution of stray light to large datasets from the NuSTAR X-ray observatory collected over 2012--2017, we report a measurement of the cosmic X-ray background in the 3--20 keV energy range. These data represent $\sim20\%$ sky coverage while avoiding Galactic Ridge X-ray emission and are less weighted by deep, survey fields than previous measurements with NuSTAR. Images in… ▽ More By characterizing the contribution of stray light to large datasets from the NuSTAR X-ray observatory collected over 2012--2017, we report a measurement of the cosmic X-ray background in the 3--20 keV energy range. These data represent $\sim20\%$ sky coverage while avoiding Galactic Ridge X-ray emission and are less weighted by deep, survey fields than previous measurements with NuSTAR. Images in narrow energy bands are stacked in detector space and spatially fit with a model representing the stray light and uniform pattern expected from the cosmic X-ray background and the instrumental background, respectively. We establish baseline flux values from Earth-occulted data and validate the fitting method on stray light observations of the Crab, which further serve to calibrate the resulting spectra. We present independent spectra of the cosmic X-ray background with the FPMA and FPMB detector arrays, which are in excellent agreement with the canonical characterization by HEAO 1 and are $10\%$ lower than most subsequent measurements; $F_{\rm{3-20~keV}}^{FPMA} = 2.63 \times 10^{-11}~\rm{erg~s^{-1}~cm^{-2}~deg^{-2}}$ and $F_{\rm{3-20~keV}}^{FPMB} = 2.58 \times 10^{-11}~\rm{erg~s^{-1}~cm^{-2}~deg^{-2}}$. We discuss these results in light of previous measurements of the cosmic X-ray background and consider the impact of systematic uncertainties on our spectra. △ Less

Submitted 19 April, 2023; v1 submitted 16 April, 2023; originally announced April 2023.

Comments: 23 pages, 15 figures, 4 tables

Higher-order Organization in the Human Brain from Matrix-Based Rényi's Entropy

Authors: Qiang Li, Shujian Yu, Kristoffer H Madsen, Vince D Calhoun, Armin Iraji

Abstract: Pairwise metrics are often employed to estimate statistical dependencies between brain regions, however they do not capture higher-order information interactions. It is critical to explore higher-order interactions that go beyond paired brain areas in order to better understand information processing in the human brain. To address this problem, we applied multivariate mutual information, specifica… ▽ More Pairwise metrics are often employed to estimate statistical dependencies between brain regions, however they do not capture higher-order information interactions. It is critical to explore higher-order interactions that go beyond paired brain areas in order to better understand information processing in the human brain. To address this problem, we applied multivariate mutual information, specifically, Total Correlation and Dual Total Correlation to reveal higher-order information in the brain. In this paper, we estimate these metrics using matrix-based Rényi's entropy, which offers a direct and easily interpretable approach that is not limited by direct assumptions about probability distribution functions of multivariate time series. We applied these metrics to resting-state fMRI data in order to examine higher-order interactions in the brain. Our results showed that the higher-order information interactions captured increase gradually as the interaction order increases. Furthermore, we observed a gradual increase in the correlation between the Total Correlation and Dual Total Correlation as the interaction order increased. In addition, the significance of Dual Total Correlation values compared to Total Correlation values also indicate that the human brain exhibits synergy dominance during the resting state. △ Less

Submitted 25 April, 2023; v1 submitted 21 March, 2023; originally announced March 2023.

Comments: 5 pages, 3 figures; Accepted to Data Science and Learning Workshop: Unraveling the Brain. A satellite workshop of ICASSP 2023

Journal ref: 2023 IEEE International Conference on Acoustics, Speech, and Signal Processing Workshops (ICASSPW)

Towards Phytoplankton Parasite Detection Using Autoencoders

Authors: Simon Bilik, Daniel Batrakhanov, Tuomas Eerola, Lumi Haraguchi, Kaisa Kraft, Silke Van den Wyngaert, Jonna Kangas, Conny Sjöqvist, Karin Madsen, Lasse Lensu, Heikki Kälviäinen, Karel Horak

Abstract: Phytoplankton parasites are largely understudied microbial components with a potentially significant ecological impact on phytoplankton bloom dynamics. To better understand their impact, we need improved detection methods to integrate phytoplankton parasite interactions in monitoring aquatic ecosystems. Automated imaging devices usually produce high amount of phytoplankton image data, while the oc… ▽ More Phytoplankton parasites are largely understudied microbial components with a potentially significant ecological impact on phytoplankton bloom dynamics. To better understand their impact, we need improved detection methods to integrate phytoplankton parasite interactions in monitoring aquatic ecosystems. Automated imaging devices usually produce high amount of phytoplankton image data, while the occurrence of anomalous phytoplankton data is rare. Thus, we propose an unsupervised anomaly detection system based on the similarity of the original and autoencoder-reconstructed samples. With this approach, we were able to reach an overall F1 score of 0.75 in nine phytoplankton species, which could be further improved by species-specific fine-tuning. The proposed unsupervised approach was further compared with the supervised Faster R-CNN based object detector. With this supervised approach and the model trained on plankton species and anomalies, we were able to reach the highest F1 score of 0.86. However, the unsupervised approach is expected to be more universal as it can detect also unknown anomalies and it does not require any annotated anomalous data that may not be always available in sufficient quantities. Although other studies have dealt with plankton anomaly detection in terms of non-plankton particles, or air bubble detection, our paper is according to our best knowledge the first one which focuses on automated anomaly detection considering putative phytoplankton parasites or infections. △ Less

Submitted 17 August, 2023; v1 submitted 15 March, 2023; originally announced March 2023.

A broadband X-ray imaging spectroscopy in the 2030s: the FORCE mission

Authors: Koji Mori, Takeshi G. Tsuru, Kazuhiro Nakazawa, Yoshihiro Ueda, Shin Watanabe, Takaaki Tanaka, Manabu Ishida, Hironori Matsumoto, Hisamitsu Awaki, Hiroshi Murakami, Masayoshi Nobukawa, Ayaki Takeda, Yasushi Fukazawa, Hiroshi Tsunemi, Tadayuki Takahashi, Ann Hornschemeier, Takashi Okajima, William W. Zhang, Brian J. Williams, Tonia Venters, Kristin Madsen, Mihoko Yukita, Hiroki Akamatsu, Aya Bamba, Teruaki Enoto , et al. (27 additional authors not shown)

Abstract: In this multi-messenger astronomy era, all the observational probes are improving their sensitivities and overall performance. The Focusing on Relativistic universe and Cosmic Evolution (FORCE) mission, the product of a JAXA/NASA collaboration, will reach a 10 times higher sensitivity in the hard X-ray band ($E >$ 10~keV) in comparison with any previous hard X-ray missions, and provide simultaneou… ▽ More In this multi-messenger astronomy era, all the observational probes are improving their sensitivities and overall performance. The Focusing on Relativistic universe and Cosmic Evolution (FORCE) mission, the product of a JAXA/NASA collaboration, will reach a 10 times higher sensitivity in the hard X-ray band ($E >$ 10~keV) in comparison with any previous hard X-ray missions, and provide simultaneous soft X-ray coverage. FORCE aims to be launched in the early 2030s, providing a perfect hard X-ray complement to the ESA flagship mission Athena. FORCE will be the most powerful X-ray probe for discovering obscured/hidden black holes and studying high energy particle acceleration in our Universe and will address how relativistic processes in the universe are realized and how these affect cosmic evolution. FORCE, which will operate over 1--79 keV, is equipped with two identical pairs of supermirrors and wideband X-ray imagers. The mirror and imager are connected by a high mechanical stiffness extensible optical bench with alignment monitor systems with a focal length of 12~m. A light-weight silicon mirror with multi-layer coating realizes a high angular resolution of $<15''$ in half-power diameter in the broad bandpass. The imager is a hybrid of a brand-new SOI-CMOS silicon-pixel detector and a CdTe detector responsible for the softer and harder energy bands, respectively. FORCE will play an essential role in the multi-messenger astronomy in the 2030s with its broadband X-ray sensitivity. △ Less

Submitted 13 March, 2023; originally announced March 2023.

Comments: 12 pages, 8 figures. Proceedings of SPIE Astronomical Telescopes and Instrumentation 2022

Journal ref: Proc. of SPIE vol. 12181, 1218122 (2022)

Deep Ensembles vs. Committees for Uncertainty Estimation in Neural-Network Force Fields: Comparison and Application to Active Learning

Authors: Jesús Carrete, Hadrián Montes-Campos, Ralf Wanzenböck, Esther Heid, Georg K. H. Madsen

Abstract: A reliable uncertainty estimator is a key ingredient in the successful use of machine-learning force fields for predictive calculations. Important considerations are correlation with error, overhead during training and inference, and efficient workflows to systematically improve the force field. However, in the case of neural-network force fields, simple committees are often the only option consid… ▽ More A reliable uncertainty estimator is a key ingredient in the successful use of machine-learning force fields for predictive calculations. Important considerations are correlation with error, overhead during training and inference, and efficient workflows to systematically improve the force field. However, in the case of neural-network force fields, simple committees are often the only option considered due to their easy implementation. Here we present a generalization of the deep-ensemble design, based on multiheaded neural networks and a heteroscedastic loss, that can efficiently deal with uncertainties in both the energy and the forces. We compare uncertainty metrics based on deep ensembles, committees and bootstrap-aggregation ensembles using data for an ionic liquid and a perovskite surface. We demonstrate an adversarial approach to active learning to efficiently and progressively refine the force fields. That active learning workflow is realistically possible thanks to exceptionally fast training enabled by residual learning and a nonlinear learned optimizer. △ Less

Submitted 17 February, 2023; originally announced February 2023.

Managing Expert Disagreement for the Policy Process and Beyond

Authors: Ulrike Hahn, Jens Koed Madsen, Chris Reed

Abstract: In this paper, we outline a new proposal for communicating scientific debate to policymakers and other stakeholders in circumstances where there is substantial disagreement within the scientific literature. In those circumstances, it seems important to provide policy makers both with a useful, balanced summary that is representative of opinion in the field large, and to transparently communicate t… ▽ More In this paper, we outline a new proposal for communicating scientific debate to policymakers and other stakeholders in circumstances where there is substantial disagreement within the scientific literature. In those circumstances, it seems important to provide policy makers both with a useful, balanced summary that is representative of opinion in the field large, and to transparently communicate the actual evidence-base. To this end, we propose the compilation of argument maps through a collective intelligence process; these maps are then given to a wide sample of the relevant research community for evaluation and summary opinion in an IGM style IGM style poll (see igmchicago.org), which provides a representative view of opinion on the issue at stake within the wider scientific community. Policymakers then receive these two artefacts (map and poll) as their expert advice. Such a process would help overcome the resource limitations of the traditional expert advice process, while also providing greater balance by drawing on the expertise of researchers beyond the leading proponents of particular theories within a field. And, the actual evidence base would be transparent. In this paper, we present a pilot project stepping through the map building component of such a policy advice scheme. We detail process, products, and issues encountered by implementing in the OVA (Online Visualisation of Argument tool, ova.arg-tech.org) an argument map with sample evidence from the behavioural literature on communicating probabilities, as a central issue within pandemic. △ Less

Submitted 28 December, 2022; originally announced December 2022.

Comments: 12 pages, 1 Figure

NuSTAR spectral analysis beyond 79 keV with stray light

Authors: G. Mastroserio, B. W. Grefenstette, P. Thalhammer, D. J. K. Buisson, M. C. Brumback, R. M. Ludlam, R. M. T. Connors, J. A. Garcıa, V. Grinberg, K. K. Madsen, H. Miyasaka, J. A. Tomsick, J. Wilms

Abstract: Due to the structure of the NuSTAR telescope, photons at large off-axis (> 1deg) can reach the detectors directly (stray light), without passing through the instrument optics. At these off-axis angles NuSTAR essentially turns into a collimated instrument and the spectrum can extend to energies above the Pt k-edge (79 keV) of the multi-layers, which limits the effective area bandpass of the optics.… ▽ More Due to the structure of the NuSTAR telescope, photons at large off-axis (> 1deg) can reach the detectors directly (stray light), without passing through the instrument optics. At these off-axis angles NuSTAR essentially turns into a collimated instrument and the spectrum can extend to energies above the Pt k-edge (79 keV) of the multi-layers, which limits the effective area bandpass of the optics. We present the first scientific spectral analysis beyond 79 keV using a Cygnus X-1 observation in StrayCats, the catalog of stray light observations. This serendipitous stray light observation occurred simultaneously with an INTEGRAL observation. When the spectra are modeled together in the 30-120 keV energy band, we find that the NuSTAR stray light flux is well calibrated and constrained to be consistent with the INTEGRAL flux at the 90% confidence level. Furthermore, we explain how to treat the background of the stray light spectral analysis, which is especially important at high energies. △ Less

Submitted 28 September, 2022; originally announced September 2022.

Comments: Accepted for publication in ApJ

AT2019wxt: An ultra-stripped supernova candidate discovered in electromagnetic follow-up of a gravitational wave trigger

Authors: Hinna Shivkumar, Amruta D. Jaodand, Arvind Balasubramanian, Christoffer Fremling, Alessandra Corsi, Anastasios Tzanidakis, Samaya Nissanke, Mansi Kasliwal, Murray Brightman, Geert Raaijmakers, Kristin Kruse Madsen, Fiona Harrison, Dario Carbone, Nayana A. J., Jean-Michel Désert, Igor Andreoni

Abstract: We present optical, radio and X-ray observations of a rapidly-evolving transient AT2019wxt (PS19hgw), discovered during the search for an electromagnetic (EM) counterpart to the gravitational-wave (GW) trigger S191213g (LIGO Scientific Collaboration & Virgo Collaboration 2019a). Although S191213g was not confirmed as a significant GW event in the off-line analysis of LIGO-Virgo data, AT2019wxt rem… ▽ More We present optical, radio and X-ray observations of a rapidly-evolving transient AT2019wxt (PS19hgw), discovered during the search for an electromagnetic (EM) counterpart to the gravitational-wave (GW) trigger S191213g (LIGO Scientific Collaboration & Virgo Collaboration 2019a). Although S191213g was not confirmed as a significant GW event in the off-line analysis of LIGO-Virgo data, AT2019wxt remained an interesting transient due its peculiar nature. The optical/NIR light curve of AT2019wxt displayed a double-peaked structure evolving rapidly in a manner analogous to currently know ultra-stripped supernovae (USSNe) candidates. This double-peaked structure suggests presence of an extended envelope around the progenitor, best modelled with two-components: i) early-time shock-cooling emission and ii) late-time radioactive $^{56}$Ni decay. We constrain the ejecta mass of AT2019wxt at $M_{ej} \approx{0.20 M_{\odot}}$ which indicates a significantly stripped progenitor that was possibly in a binary system. We also followed-up AT2019wxt with long-term Chandra and Jansky Very Large Array observations spanning $\sim$260 days. We detected no definitive counterparts at the location of AT2019wxt in these long-term X-ray and radio observational campaigns. We establish the X-ray upper limit at $9.93\times10^{-17}$ erg cm$^{-2}$ s$^{-1}$ and detect an excess radio emission from the region of AT2019wxt. However, there is little evidence for SN1993J- or GW170817-like variability of the radio flux over the course of our observations. A substantial host galaxy contribution to the measured radio flux is likely. The discovery and early-time peak capture of AT2019wxt in optical/NIR observation during EMGW follow-up observations highlights the need of dedicated early, multi-band photometric observations to identify USSNe. △ Less

Submitted 18 August, 2022; originally announced August 2022.

Comments: 20 pages, 14 figures, Submitted to ApJ

The Low Temperature Corona in ESO 511$-$G030 Revealed by NuSTAR and XMM-Newton

Authors: Zuobin Zhang, Jiachen Jiang, Honghui Liu, Cosimo Bambi, Christopher S. Reynolds, Andrew C. Fabian, Thomas Dauser, Kristin Madsen, Andrew Young, Luigi Gallo, Zhibo Yu, John Tomsick

Abstract: We present the results from a coordinated XMM-Newton $+$ NuSTAR observation of the Seyfert 1 Galaxy ESO 511$-$G030. With this joint monitoring programme, we conduct a detailed variability and spectral analysis. The source remained in a low flux and very stable state throughout the observation period, although there are slight fluctuations of flux over long timescales. The broadband (0.3-78~keV) sp… ▽ More We present the results from a coordinated XMM-Newton $+$ NuSTAR observation of the Seyfert 1 Galaxy ESO 511$-$G030. With this joint monitoring programme, we conduct a detailed variability and spectral analysis. The source remained in a low flux and very stable state throughout the observation period, although there are slight fluctuations of flux over long timescales. The broadband (0.3-78~keV) spectrum shows the presence of a power-law continuum with a soft excess below 2~keV, a relatively narrow iron K$α$ emission ($\sim$6.4~keV), and an obvious cutoff at high energies. We find that the soft excess can be modeled by two different possible scenarios: a warm ($kT_{\rm e} \sim$ 0.19~keV) and optically thick ($τ- 18\sim25$) Comptonizing corona or a relativistic reflection from a high-density ($\log [n_{\rm e}/{\rm cm}^{-3}]=17.1 \sim 18.5$) inner disc. All models require a low temperature ($kT_{\rm e} \sim$ 13~keV) for the hot corona. △ Less

Submitted 18 May, 2023; v1 submitted 2 August, 2022; originally announced August 2022.

Comments: 18 pages, 15 figures. v2: refereed version

Journal ref: Astrophys.J. 949: 4 (2023)

StrayCats II: An Updated Catalog of NuSTAR Stray Light Observations

Authors: R. M. Ludlam, B. W. Grefenstette, M. C. Brumback, J. A. Tomsick, D. J. K. Buisson, B. M. Coughenour, G. Mastroserio, D. Wik, R. Krivonos, A. D. Jaodand, K. K. Madsen

Abstract: We present an updated catalog of StrayCats (a catalog of NuSTAR stray light observations of X-ray sources) that includes nearly 18 additional months of observations. StrayCats v2 has an added 53 sequence IDs, 106 rows, and 3 new identified stray light (SL) sources in comparison to the original catalog. The total catalog now has 489 unique sequence IDs, 862 entries, and 83 confirmed StrayCats sourc… ▽ More We present an updated catalog of StrayCats (a catalog of NuSTAR stray light observations of X-ray sources) that includes nearly 18 additional months of observations. StrayCats v2 has an added 53 sequence IDs, 106 rows, and 3 new identified stray light (SL) sources in comparison to the original catalog. The total catalog now has 489 unique sequence IDs, 862 entries, and 83 confirmed StrayCats sources. Additionally, we provide new resources for the community to gauge the utility and spectral state of the source in a given observation. We have created long term light curves for each identified SL source using MAXI and Swift/BAT data when available. Further, source extraction regions for 632 identified SL observations were created and are available to the public. In this paper we present an overview of the updated catalog and new resources for each identified StrayCats SL source. △ Less

Submitted 20 June, 2022; originally announced June 2022.

Comments: Accepted for publication in ApJ, 6 pages, 5 figures, 1 table

Measuring the Evolution of the NuSTAR Detector Gains

Authors: Brian Grefenstette, Murray Brightman, Hannah P. Earnshaw, Karl Forster, Kristin K. Madsen, Hiromasa Miyasaka

Abstract: The memo describes the methods used to track the long-term gain variations in the NuSTAR detectors. It builds on the analysis presented in Madsen et al. (2015) using the deployable calibration source to measure the gain drift in the NuSTAR CdZnTe detectors. This is intended to be a live document that is periodically updated as new entries are required in the NuSTAR gain CALDB files. This document… ▽ More The memo describes the methods used to track the long-term gain variations in the NuSTAR detectors. It builds on the analysis presented in Madsen et al. (2015) using the deployable calibration source to measure the gain drift in the NuSTAR CdZnTe detectors. This is intended to be a live document that is periodically updated as new entries are required in the NuSTAR gain CALDB files. This document covers analysis up through early-2022 and the gain v011 CALDB file released in version 20240226. △ Less

Submitted 27 February, 2024; v1 submitted 8 June, 2022; originally announced June 2022.

Comments: 11 page, 7 figures. Intended as a living, easy-to-find document. No intention of submitting this to a journal

Extending the baseline for SMC X-1's spin and orbital behavior with NuSTAR stray light

Authors: McKinley C. Brumback, B. W. Grefenstette, D. J. K. Buisson, M. Bachetti, R. Connors, J. A. Garcia, A. Jaodand, R. Krivonos, R. Ludlam, K. K. Madsen, G. Mastroserio, J. A. Tomsick, D. Wik

Abstract: StrayCats, the catalog of NuSTAR stray light observations, contains data from bright X-ray sources that fall within crowded source regions. These observations offer unique additional data with which to monitor sources like X-ray binaries that show variable timing behavior. In this work, we present a timing analysis of stray light data of the high mass X-ray binary SMC X-1, the first scientific ana… ▽ More StrayCats, the catalog of NuSTAR stray light observations, contains data from bright X-ray sources that fall within crowded source regions. These observations offer unique additional data with which to monitor sources like X-ray binaries that show variable timing behavior. In this work, we present a timing analysis of stray light data of the high mass X-ray binary SMC X-1, the first scientific analysis of a single source from the StrayCats project. We describe the process of screening stray light data for scientific analysis, verify the orbital ephemeris, and create both time and energy resolved pulse profiles. We find that the orbital ephemeris of SMC X-1 is unchanged and confirm a long-term spin up rate of $\dotν=(2.52\pm0.03)\times10^{-11}$ Hz s$^{-1}$. We also note that the shape of SMC X-1's pulse profile, while remaining double-peaked, varies significantly with time and only slightly with energy. △ Less

Submitted 22 February, 2022; originally announced February 2022.

Comments: 10 pages, 8 figures

Reconstruction of the NuSTAR point spread function using single-laser metrology

Authors: Hannah P. Earnshaw, Kristin K. Madsen, Karl Forster, Brian W. Grefenstette, Murray Brightman, Andreas Zoglauer, Fiona Harrison

Abstract: This paper describes a method by which the metrology system of the Nuclear Spectroscopic Telescope Array (NuSTAR) X-ray space observatory, which uses two lasers to characterize the relative motion of the optics and focal plane benches, can be approximated should one laser fail. The two benches are separated by a ten-meter-long rigid mast that undergoes small amounts of thermal flexing which need t… ▽ More This paper describes a method by which the metrology system of the Nuclear Spectroscopic Telescope Array (NuSTAR) X-ray space observatory, which uses two lasers to characterize the relative motion of the optics and focal plane benches, can be approximated should one laser fail. The two benches are separated by a ten-meter-long rigid mast that undergoes small amounts of thermal flexing which need to be compensated for in order to produce a non-blurred image. We analyze the trends of mast motion by archival observation parameters in order to discover whether the mast motion in future observations can be predicted. We find that, by using the solar aspect angle (SAA), observation date, and orbital phase, we can simulate the motion of one laser by translating the track produced by the other and applying modifications to the resulting mast aspect solution, allowing the reconstruction of a minimally distorted point spread function in most cases. We will implement the generation of simulated mast files alongside the usual NuSTAR data reduction pipeline for contingency purposes. This work has implications for reducing the risk of implementing laser metrology systems on future missions that use deployable masts to achieve the long focal lengths required in high-energy astronomy by mitigating the impact of a metrology laser failure in the extended phase of a mission. △ Less

Submitted 9 February, 2022; originally announced February 2022.

Comments: 36 pages, 20 figures, 1 table. Accepted for publication in JATIS

Electric Field Models of Transcranial Magnetic Stimulation Coils with Arbitrary Geometries: Reconstruction from Incomplete Magnetic Field Measurements

Authors: Kristoffer Hougaard Madsen, Maria Drakaki, Axel Thielscher

Abstract: Background: Calculation of the electric field induced by transcranial magnetic stimulation (TMS) in the brain requires accurate models of the stimulation coils. Reconstructing models from the measured magnetic fields of coils so far worked only for flat coil geometries and required data of their full magnetic field distribution. Objective: To reconstruct models of coils with arbitrary winding ge… ▽ More Background: Calculation of the electric field induced by transcranial magnetic stimulation (TMS) in the brain requires accurate models of the stimulation coils. Reconstructing models from the measured magnetic fields of coils so far worked only for flat coil geometries and required data of their full magnetic field distribution. Objective: To reconstruct models of coils with arbitrary winding geometries from spatially incomplete magnetic field measurements. Methods: Dipole approximation via minimum norm estimation with a cross validation procedure simultaneously assessing predictability and reproducibility of the field approximation. The methods were validated on both simulated and acquired magnetic flux density data. Furthermore, reconstruction of coil models from sparsely sampled data is investigated and a procedure for obtaining sparse dipole coil models based on orthogonal matching pursuit is suggested. Results: Given that measured data from regions around the coil are available, the magnetic vector potential can accurately be represented by fitted dipole models even from sparsely sampled data. Considerable sparsification of dipole models is possible while retaining accurate representation of the field also close to the coil. Conclusion: Our versatile approach removes an important hurdle for constructing coil models from measurement data and enables a flexible trade-off between the accuracy and computational efficiency of the reconstructed dipole models. △ Less

Submitted 29 December, 2021; originally announced December 2021.

Orbital decay in M82 X-2

Authors: Matteo Bachetti, Marianne Heida, Thomas Maccarone, Daniela Huppenkothen, Gian Luca Israel, Didier Barret, Murray Brightman, McKinley Brumback, Hannah P. Earnshaw, Karl Forster, Felix Fürst, Brian W. Grefenstette, Fiona A. Harrison, Amruta D. Jaodand, Kristin K. Madsen, Matthew Middleton, Sean N. Pike, Maura Pilia, Juri Poutanen, Daniel Stern, John A. Tomsick, Dominic J. Walton, Natalie Webb, Jörn Wilms

Abstract: M82 X-2 is the first pulsating ultraluminous X-ray source (PULX) discovered. The luminosity of these extreme pulsars, if isotropic, implies an extreme mass transfer rate. An alternative is to assume a much lower mass transfer rate, but with an apparent luminosity boosted by geometrical beaming. Only an independent measurement of the mass transfer rate can help discriminate between these two scenar… ▽ More M82 X-2 is the first pulsating ultraluminous X-ray source (PULX) discovered. The luminosity of these extreme pulsars, if isotropic, implies an extreme mass transfer rate. An alternative is to assume a much lower mass transfer rate, but with an apparent luminosity boosted by geometrical beaming. Only an independent measurement of the mass transfer rate can help discriminate between these two scenarios. In this Paper, we follow the orbit of the neutron star for seven years, measure the decay of the orbit ($\dot{P}_{orb}/{P}_{orb}\approx-8\cdot10^{-6}\mathrm{yr}^{-1}$), and argue that this orbital decay is driven by extreme mass transfer of more than 150 times the mass transfer limit set by the Eddington luminosity. If this is true, the mass available to the accretor is more than enough to justify its luminosity, with no need for beaming. This also strongly favors models where the accretor is a highly-magnetized neutron star. △ Less

Submitted 29 August, 2022; v1 submitted 1 December, 2021; originally announced December 2021.

Comments: 18 pages, 6 figures, accepted for publication in ApJ

Journal ref: ApJ 937, 125 (2022)

IACHEC 2020/2021 Pandemic Report

Authors: K. K. Madsen, V. Burwitz, K. Forster, C. E. Grant, M. Guainazzi, V. Kashyap, H. L. Marshall, E. D. Miller, L. Natalucci, P. P. Plucinsky, Y. Terada

Abstract: In this report we summarize the activities of the International Astronomical Consortium for High Energy Calibration (IACHEC) and the work done since the last in-person meeting in Japan (Shonan Village Center), May 2019, through two virtual meetings that were held in November 2020 and May 2021. The on-line only meetings divided the contents of the usual in-person workshop between mission updates an… ▽ More In this report we summarize the activities of the International Astronomical Consortium for High Energy Calibration (IACHEC) and the work done since the last in-person meeting in Japan (Shonan Village Center), May 2019, through two virtual meetings that were held in November 2020 and May 2021. The on-line only meetings divided the contents of the usual in-person workshop between mission updates and working group updates. The November meeting was dedicated to mission calibration updates and the current status of the cross-calibration between NuSTAR, Swift, and NICER, which frequently join together in observations of bright transients, and a review of the XMM-Newton and Chandra cross-calibration. Results between \nustar\ and \swift\ overall show good agreement, but issues persist in the overlap region 3--5 keV for bright source with large dust scattering halos. The NICER cross-calibration is still progressing and evolving, while for the XMM-Newton and Chandra cross-calibration systematic differences both in the absolute flux and spectral shape determination still exists on different classes of sources. The meeting in May was focused on the Working Group progress and reports summarized here. △ Less

Submitted 2 November, 2021; originally announced November 2021.

Comments: Yearly activity report of the IACHEC, 16 pages, 2 figures

2021 Effective Area calibration of the Nuclear Spectroscopic Telescope ARray (NuSTAR)

Authors: Kristin K. Madsen, Karl Forster, Brian W. Grefenstette, Fiona A. Harrison, Hiromasa Miyasaka

Abstract: We present here the updated calibration of The Nuclear Spectroscopic Telescope ARray NuSTAR, which was performed using data on the Crab accumulated over the last 9 years in orbit. The basis for this new calibration contains over 250ks of focused Crab (imaged through the optics) and over 500ks of stray-light Crab (not imaged through optics). We measured an epoch averaged Crab spectrum of the stray-… ▽ More We present here the updated calibration of The Nuclear Spectroscopic Telescope ARray NuSTAR, which was performed using data on the Crab accumulated over the last 9 years in orbit. The basis for this new calibration contains over 250ks of focused Crab (imaged through the optics) and over 500ks of stray-light Crab (not imaged through optics). We measured an epoch averaged Crab spectrum of the stray-light Crab data and define a canonical Crab spectrum of Gamma = 2.103 +- 0.001 and N = 9.69 +- 0.02 keV-1 cm-2 s-1 at 1 keV, which we use as our calibration standard. The new calibration, released in the CALDB update 20211020, provides significant updates to: 1) the detector absorption component, 2) the detector response function, and 3) the effective area vignetting function. The calibration improves agreement between FPMA and FPMB across detectors with a standard deviation of 1.7% for repeat observations between off-axis angles of 1-4 arcmin, and the measured flux has increased by 5-15%, with 5% below 1 arcmin off-axis angle, 10% between 1-2 arcmin, and 15 above 4arcmin. △ Less

Submitted 21 October, 2021; originally announced October 2021.

Comments: 26 pages, 18 figures

Accurate effective harmonic potential treatment of the high-temperature cubic phase of Hafnia

Authors: Sebastian Bichelmaier, Jesús Carrete, Michael Nelhiebel, Georg K. H. Madsen

Abstract: HfO$_2$ is an important high-$κ$ dielectric and ferroelectric, exhibiting a complex potential energy landscape with several phases close in energy. It is, however, a strongly anharmonic solid, and thus describing its temperature-dependent behavior is methodologically challenging. We propose an approach based on self-consistent, effective harmonic potentials and higher-order corrections to study th… ▽ More HfO$_2$ is an important high-$κ$ dielectric and ferroelectric, exhibiting a complex potential energy landscape with several phases close in energy. It is, however, a strongly anharmonic solid, and thus describing its temperature-dependent behavior is methodologically challenging. We propose an approach based on self-consistent, effective harmonic potentials and higher-order corrections to study the potential energy surface of anharmonic materials. The introduction of a reweighting procedure enables the usage of unregularized regression methods and efficiently harnesses the information contained in every data point obtained from density functional theory. This renders the approach highly efficient and a promising candidate for large-scale studies of materials and phase transitions. We detail the approach and test it on the example of the high-temperature cubic phase of HfO$_2$. Our results for the thermal expansion coefficient, $α_V \approx 3.3\times 10^{-5}$ K$^{-1}$, are in agreement with existing experimental ($α_V \approx 4\pm1\times 10^{-5}$ K$^{-1}$) and theoretical ($α_V \approx 5\pm1\times 10^{-5}$ K$^{-1}$) work. Likewise, the bulk modulus agrees well with experiment. We show the detailed temperature dependence of these quantities. △ Less

Submitted 10 October, 2021; originally announced October 2021.

Signaling Design for Cooperative Resource Allocation and its Impact to Reliability

Authors: Rasmus Liborius Bruun, C. Santiago Morejón García, Troels B. Sørensen, Nuno K. Pratas, Tatiana Kozlova Madsen, Preben Mogensen

Abstract: Decentralized cooperative resource allocation schemes for robotic swarms are essential to enable high reliability in high throughput data exchanges. These cooperative schemes require control signaling with the aim to avoid half-duplex problems at the receiver and mitigate interference. We propose two cooperative resource allocation schemes, device sequential and group scheduling, and introduce a c… ▽ More Decentralized cooperative resource allocation schemes for robotic swarms are essential to enable high reliability in high throughput data exchanges. These cooperative schemes require control signaling with the aim to avoid half-duplex problems at the receiver and mitigate interference. We propose two cooperative resource allocation schemes, device sequential and group scheduling, and introduce a control signaling design. We observe that failure in the reception of these control signals leads to non-cooperative behavior and to significant performance degradation. The cause of these failures are identified and specific countermeasures are proposed and evaluated. We compare the proposed resource allocation schemes against the NR sidelink mode 2 resource allocation and show that even though signaling has an important impact on the resource allocation performance, our proposed device sequential and group scheduling resource allocation schemes improve reliability by an order of magnitude compared to sidelink mode 2. △ Less

Submitted 15 September, 2022; v1 submitted 15 September, 2021; originally announced September 2021.

Magnetic correlations in triangular antiferromagnet FeGa$_2$S$_4$

Authors: K. Guratinder, M. Schmidt, H. C. Walker, R. Bewley, M. Wörle, D. Cabra, S. A. Osorio, M. Villalba, A. K. Madsen, L. Keller, A. Wildes, P. Puphal, A. Cervellino, Ch. Rüegg, O. Zaharko

Abstract: The crystal structure and magnetic correlations in triangular antiferromagnet FeGa$_2$S$_4$ are studied by x-ray diffraction, magnetic susceptibility, neutron diffraction and neutron inelastic scattering. We report significant mixing at the cation sites and disentangle magnetic properties dominated by major and minor magnetic sites. The magnetic short-range correlations at 0.77 Å$^{-1}$ correspond… ▽ More The crystal structure and magnetic correlations in triangular antiferromagnet FeGa$_2$S$_4$ are studied by x-ray diffraction, magnetic susceptibility, neutron diffraction and neutron inelastic scattering. We report significant mixing at the cation sites and disentangle magnetic properties dominated by major and minor magnetic sites. The magnetic short-range correlations at 0.77 Å$^{-1}$ correspond to the major sites and being static at base temperature they evolve into dynamic correlations around 30 - 50 K. The minor sites contribute to the magnetic peak at 0.6 Å$^{-1}$, which vanishes at 5.5 K. Our analytical studies of triangular lattice models with bilinear and biquadratic terms provide the ratios between exchanges for the proposed ordering vectors. The modelling of the inelastic neutron spectrum within linear spin wave theory results in the set of exchange couplings $J_1=1.7$\,meV, $J_2=0.9$\,meV, $J_3=0.8$\,meV for the bilinear Heisenberg Hamiltonian. However, not all features of the excitation spectrum are explained with this model. △ Less

Submitted 14 July, 2021; originally announced July 2021.

Comments: Accepted in Phys. Rev. B

A Differentiable Neural-Network Force Field for Ionic Liquids

Authors: Hadrián Montes-Campos, Jesús Carrete, Sebastian Bichelmaier, Luis M. Varela, Georg K. H. Madsen

Abstract: We present NeuralIL, a model for the potential energy of an ionic liquid that accurately reproduces first-principles results with orders-of-magnitude savings in computational cost. Based on a multilayer perceptron and spherical Bessel descriptors of the atomic environments, NeuralIL is implemented in such a way as to be fully automatically differentiable. It can thus be trained on ab-initio forces… ▽ More We present NeuralIL, a model for the potential energy of an ionic liquid that accurately reproduces first-principles results with orders-of-magnitude savings in computational cost. Based on a multilayer perceptron and spherical Bessel descriptors of the atomic environments, NeuralIL is implemented in such a way as to be fully automatically differentiable. It can thus be trained on ab-initio forces instead of just energies, to make the most out of the available data, and can efficiently predict arbitrary derivatives of the potential energy. Using ethylammonium nitrate as the test system, we obtain out-of-sample accuracies better than 2 meV/atom (<0.05 kcal/mol) in the energies and 70 meV/Å in the forces. We show that encoding the element specific density in the spherical Bessel descriptors is key to achieving this. Harnessing the information provided by the forces drastically reduces the amount of atomic configurations required to train a neural network force field based on atom-centered descriptors. We choose the Swish-1 activation function and discuss the role of this choice in keeping the neural network differentiable. Furthermore, the possibility of training on small data sets allows for an ensemble-learning approach to the detection of extrapolation. Finally, we find that a separate treatment of long-range interactions is not required to achieve a high-quality representation of the potential energy surface of these dense ionic systems. △ Less

Submitted 11 November, 2021; v1 submitted 30 June, 2021; originally announced June 2021.