Showing 1–50 of 299 results for author: Kraemer, S

Determining the Extents, Geometries, and Kinematics of Narrow-Line Region Outflows in Nearby Seyfert Galaxies

Authors: Garrett E. Polack, Mitchell Revalski, D. Michael Crenshaw, Travis C. Fischer, Henrique R. Schmitt, Steven B. Kraemer, Beena Meena, Marc Rafelski

Abstract: Outflowing gas from supermassive black holes in the centers of active galaxies has been postulated as a major contributor to galactic evolution. To explore the interaction between narrow-line region (NLR) outflows and their host galaxies, we use Hubble Space Telescope (HST) Space Telescope Imaging Spectrograph (STIS) spectra and Wide Field Camera 3 (WFC3) images of 15 nearby (z < 0.02) active gala… ▽ More Outflowing gas from supermassive black holes in the centers of active galaxies has been postulated as a major contributor to galactic evolution. To explore the interaction between narrow-line region (NLR) outflows and their host galaxies, we use Hubble Space Telescope (HST) Space Telescope Imaging Spectrograph (STIS) spectra and Wide Field Camera 3 (WFC3) images of 15 nearby (z < 0.02) active galactic nuclei (AGN) to determine the extents and geometries of their NLRs. We combine new HST WFC3 continuum and [O III] $λ$5007A images of 11 AGN with 4 archival AGN to match existing spectra from HST STIS. For the 6 AGN with suitable long-slit coverage of their NLRs, we use isophotal fitting of ground-based images, continuum-subtracted [O III] images, and the STIS spectra, to resolve, measure, and de-project the gas kinematics to the plane of the host galaxy disk and distinguish NLR outflows from galaxy rotation and/or kinematically disturbed gas. We find an average [O III] extent of $\sim$680pc with a correlation between gas extent and [O III] luminosity of R$_\mathrm{[O III]}$ $\propto$ L$_{\text{[O III]}}^{0.39}$. The measured extents depend strongly on the depth of the [O III] images, highlighting the importance of adopting uniform thresholds when analyzing scaling relationships. The outflows reach from 39-88% of the full NLR extents, and we find that all 6 of the AGN with STIS coverage of their entire NLRs show strong kinematic evidence for outflows, despite previous uncertainty for these AGN. This suggests that NLR outflows are ubiquitous in moderate luminosity AGN and that standard criteria for kinematic modeling are essential for identifying outflows. △ Less

Submitted 24 July, 2024; originally announced July 2024.

Comments: Accepted for publication in ApJ on July 9, 2024. The paper has 21 pages, 10 figures, and 6 tables

10 Years of Fair Representations: Challenges and Opportunities

Authors: Mattia Cerrato, Marius Köppel, Philipp Wolf, Stefan Kramer

Abstract: Fair Representation Learning (FRL) is a broad set of techniques, mostly based on neural networks, that seeks to learn new representations of data in which sensitive or undesired information has been removed. Methodologically, FRL was pioneered by Richard Zemel et al. about ten years ago. The basic concepts, objectives and evaluation strategies for FRL methodologies remain unchanged to this day. In… ▽ More Fair Representation Learning (FRL) is a broad set of techniques, mostly based on neural networks, that seeks to learn new representations of data in which sensitive or undesired information has been removed. Methodologically, FRL was pioneered by Richard Zemel et al. about ten years ago. The basic concepts, objectives and evaluation strategies for FRL methodologies remain unchanged to this day. In this paper, we look back at the first ten years of FRL by i) revisiting its theoretical standing in light of recent work in deep learning theory that shows the hardness of removing information in neural network representations and ii) presenting the results of a massive experimentation (225.000 model fits and 110.000 AutoML fits) we conducted with the objective of improving on the common evaluation scenario for FRL. More specifically, we use automated machine learning (AutoML) to adversarially "mine" sensitive information from supposedly fair representations. Our theoretical and experimental analysis suggests that deterministic, unquantized FRL methodologies have serious issues in removing sensitive information, which is especially troubling as they might seem "fair" at first glance. △ Less

Submitted 4 July, 2024; originally announced July 2024.

Towards Universal Mesh Movement Networks

Authors: Mingrui Zhang, Chunyang Wang, Stephan Kramer, Joseph G. Wallwork, Siyi Li, Jiancheng Liu, Xiang Chen, Matthew D. Piggott

Abstract: Solving complex Partial Differential Equations (PDEs) accurately and efficiently is an essential and challenging problem in all scientific and engineering disciplines. Mesh movement methods provide the capability to improve the accuracy of the numerical solution without increasing the overall mesh degree of freedom count. Conventional sophisticated mesh movement methods are extremely expensive and… ▽ More Solving complex Partial Differential Equations (PDEs) accurately and efficiently is an essential and challenging problem in all scientific and engineering disciplines. Mesh movement methods provide the capability to improve the accuracy of the numerical solution without increasing the overall mesh degree of freedom count. Conventional sophisticated mesh movement methods are extremely expensive and struggle to handle scenarios with complex boundary geometries. However, existing learning-based methods require re-training from scratch given a different PDE type or boundary geometry, which limits their applicability, and also often suffer from robustness issues in the form of inverted elements. In this paper, we introduce the Universal Mesh Movement Network (UM2N), which -- once trained -- can be applied in a non-intrusive, zero-shot manner to move meshes with different size distributions and structures, for solvers applicable to different PDE types and boundary geometries. UM2N consists of a Graph Transformer (GT) encoder for extracting features and a Graph Attention Network (GAT) based decoder for moving the mesh. We evaluate our method on advection and Navier-Stokes based examples, as well as a real-world tsunami simulation case. Our method outperforms existing learning-based mesh movement methods in terms of the benchmarks described above. In comparison to the conventional sophisticated Monge-Ampère PDE-solver based method, our approach not only significantly accelerates mesh movement, but also proves effective in scenarios where the conventional method fails. Our project page is at https://erizmr.github.io/UM2N/. △ Less

Submitted 1 July, 2024; v1 submitted 29 June, 2024; originally announced July 2024.

Multiscale simulation of spatially correlated microstructure via a latent space representation

Authors: Reese E. Jones, Craig M. Hamel, Dan Bolintineanu, Kyle Johnson, Robert Buarque de Macedo, Jan Fuhg, Nikolaos Bouklas, Sharlotte Kramer

Abstract: When deformation gradients act on the scale of the microstructure of a part due to geometry and loading, spatial correlations and finite-size effects in simulation cells cannot be neglected. We propose a multiscale method that accounts for these effects using a variational autoencoder to encode the structure-property map of the stochastic volume elements making up the statistical description of th… ▽ More When deformation gradients act on the scale of the microstructure of a part due to geometry and loading, spatial correlations and finite-size effects in simulation cells cannot be neglected. We propose a multiscale method that accounts for these effects using a variational autoencoder to encode the structure-property map of the stochastic volume elements making up the statistical description of the part. In this paradigm the autoencoder can be used to directly encode the microstructure or, alternatively, its latent space can be sampled to provide likely realizations. We demonstrate the method on three examples using the common additively manufactured material AlSi10Mg in: (a) a comparison with direct numerical simulation of the part microstructure, (b) a push forward of microstructural uncertainty to performance quantities of interest, and (c) a simulation of functional gradation of a part with stochastic microstructure. △ Less

Submitted 29 May, 2024; originally announced May 2024.

Comments: 23 pages, 25 figures

A practical guide to light-sheet microscopy for nanoscale imaging: Looking beyond the cell

Authors: Stephanie N. Kramer, Jeanpun Antarasen, Cole R. Reinholt, Lydia Kisley

Abstract: We present a comprehensive guide to light-sheet microscopy (LSM) to assist scientists in navigating the practical implementation of this microscopy technique. Emphasizing the applicability of LSM to image both static microscale and nanoscale features, as well as diffusion dynamics, we present the fundamental concepts of microscopy, progressing through beam profile considerations, to image reconstr… ▽ More We present a comprehensive guide to light-sheet microscopy (LSM) to assist scientists in navigating the practical implementation of this microscopy technique. Emphasizing the applicability of LSM to image both static microscale and nanoscale features, as well as diffusion dynamics, we present the fundamental concepts of microscopy, progressing through beam profile considerations, to image reconstruction. We outline key practical decisions in constructing a home-built system and provide insight into the alignment and calibration processes. We briefly discuss the conditions necessary for constructing a continuous 3D image and introduce our home-built code for data analysis. By providing this guide, we aim to alleviate the challenges associated with designing and constructing LSM systems and offer scientists new to LSM a valuable resource in navigating this complex field. △ Less

Submitted 14 May, 2024; originally announced May 2024.

Comments: 65 pages, 16 figures, submitted to Journal of Applied Physics

The Canadian VirusSeq Data Portal & Duotang: open resources for SARS-CoV-2 viral sequences and genomic epidemiology

Authors: Erin E. Gill, Baofeng Jia, Carmen Lia Murall, Raphaël Poujol, Muhammad Zohaib Anwar, Nithu Sara John, Justin Richardsson, Ashley Hobb, Abayomi S. Olabode, Alexandru Lepsa, Ana T. Duggan, Andrea D. Tyler, Arnaud N'Guessan, Atul Kachru, Brandon Chan, Catherine Yoshida, Christina K. Yung, David Bujold, Dusan Andric, Edmund Su, Emma J. Griffiths, Gary Van Domselaar, Gordon W. Jolly, Heather K. E. Ward, Henrich Feher , et al. (45 additional authors not shown)

Abstract: The COVID-19 pandemic led to a large global effort to sequence SARS-CoV-2 genomes from patient samples to track viral evolution and inform public health response. Millions of SARS-CoV-2 genome sequences have been deposited in global public repositories. The Canadian COVID-19 Genomics Network (CanCOGeN - VirusSeq), a consortium tasked with coordinating expanded sequencing of SARS-CoV-2 genomes acro… ▽ More The COVID-19 pandemic led to a large global effort to sequence SARS-CoV-2 genomes from patient samples to track viral evolution and inform public health response. Millions of SARS-CoV-2 genome sequences have been deposited in global public repositories. The Canadian COVID-19 Genomics Network (CanCOGeN - VirusSeq), a consortium tasked with coordinating expanded sequencing of SARS-CoV-2 genomes across Canada early in the pandemic, created the Canadian VirusSeq Data Portal, with associated data pipelines and procedures, to support these efforts. The goal of VirusSeq was to allow open access to Canadian SARS-CoV-2 genomic sequences and enhanced, standardized contextual data that were unavailable in other repositories and that meet FAIR standards (Findable, Accessible, Interoperable and Reusable). The Portal data submission pipeline contains data quality checking procedures and appropriate acknowledgement of data generators that encourages collaboration. Here we also highlight Duotang, a web platform that presents genomic epidemiology and modeling analyses on circulating and emerging SARS-CoV-2 variants in Canada. Duotang presents dynamic changes in variant composition of SARS-CoV-2 in Canada and by province, estimates variant growth, and displays complementary interactive visualizations, with a text overview of the current situation. The VirusSeq Data Portal and Duotang resources, alongside additional analyses and resources computed from the Portal (COVID-MVP, CoVizu), are all open-source and freely available. Together, they provide an updated picture of SARS-CoV-2 evolution to spur scientific discussions, inform public discourse, and support communication with and within public health authorities. They also serve as a framework for other jurisdictions interested in open, collaborative sequence data sharing and analyses. △ Less

Submitted 7 May, 2024; originally announced May 2024.

Star Formation Rates in [Ne V] 3426 Å Selected Active Galactic Nuclei: Evidence for a Decrease along the Main Sequence?

Authors: Léa Feuillet, Marcio Meléndez, Steve Kraemer, Henrique Schmitt, Travis Fischer, James Reeves

Abstract: Studying the behavior along the galaxy main sequence is key in furthering our understanding of the possible connection between AGN activity and star formation. We select a sample of 1215 AGN from the catalog of SDSS galaxy properties from the Portsmouth group by detection of the high-ionization [Ne V] 3426 Å emission line. Our sample extends from 10$^{40}$ to 10$^{42.5}$ erg/s in [Ne V] luminosity… ▽ More Studying the behavior along the galaxy main sequence is key in furthering our understanding of the possible connection between AGN activity and star formation. We select a sample of 1215 AGN from the catalog of SDSS galaxy properties from the Portsmouth group by detection of the high-ionization [Ne V] 3426 Å emission line. Our sample extends from 10$^{40}$ to 10$^{42.5}$ erg/s in [Ne V] luminosity in a redshift range z = 0.17 to 0.57. We compare the specific star formation rates (sSFRs, SFR scaled by galaxy mass) obtained from the corrected [O II] and Hα luminosities, and the SED-determined values from Portsmouth. We find that the emission-line-based sSFR values are unreliable for the [Ne V] sample due to the AGN contribution, and proceed with the SED sSFRs for our study of the main sequence. We find evidence for a decrease in sSFR along the main sequence in the [Ne V] sample which is consistent with results from the hard X-ray BAT AGN sample, which extends to lower redshifts than our [Ne V] sample. Although we do not find evidence that the concurrent AGN activity is suppressing star formation, our results are consistent with a lower gas fraction in the host galaxies of the AGN as compared to that of the star forming galaxies. If the evacuation of gas, and therefore suppression of star formation is due to AGN activity, it must have occurred in a previous epoch. △ Less

Submitted 3 May, 2024; originally announced May 2024.

Comments: 10 pages, 8 figures, 1 table, accepted for publication in ApJ

Testing Tachyon-Dominated Cosmology with Type Ia Supernovae

Authors: Samuel H. Kramer, Ian H. Redmount

Abstract: An open or hyperbolic Friedmann-Robertson-Walker spacetime dominated by tachyonic dark matter can exhibit an ``inflected'' expansion -- initially decelerating, later accelerating -- similar but not identical to that of now-standard $Λ$CDM models dominated by dark energy. The features of the tachyonic model can be extracted by fitting the redshift-distance relation of the model to data obtained by… ▽ More An open or hyperbolic Friedmann-Robertson-Walker spacetime dominated by tachyonic dark matter can exhibit an ``inflected'' expansion -- initially decelerating, later accelerating -- similar but not identical to that of now-standard $Λ$CDM models dominated by dark energy. The features of the tachyonic model can be extracted by fitting the redshift-distance relation of the model to data obtained by treating Type Ia supernovae as standard candles. Here such a model is fitted to samples of 186 and 1048 Type Ia supernovae from the literature. The fits yield values of $H_0=(66.6\pm1.5)~\hbox{km/s/Mpc}$ and $H_0=(69.6\pm0.4)~\hbox{km/s/Mpc}$, respectively, for the current-time Hubble parameter, and $t_0=(8.35\pm0.68)~\hbox{Gyr}$ and $t_0=(8.15\pm0.36)~\hbox{Gyr}$, respectively, for the comoving-time age of the Universe. Tests of the model against other observations will be undertaken in subsequent works. △ Less

Submitted 19 March, 2024; originally announced March 2024.

Comments: 20 pages, 8 figures

Supermassive Black Hole Winds in X-rays -- SUBWAYS. III. A population study on ultra-fast outflows

Authors: V. E. Gianolli, S. Bianchi, P-O Petrucci, M. Brusa, G. Chartas, G. Lanzuisi, G. A. Matzeu, M. Parra, F. Ursini, E. Behar, M. Bischetti, A. Comastri, E. Costantini, G. Cresci, M. Dadina, B. De Marco, A. De Rosa, F. Fiore, M. Gaspari, R. Gilli, M. Giustini, M. Guainazzi, A. R. King, S. Kraemer, G. Kriss , et al. (22 additional authors not shown)

Abstract: The detection of blue-shifted absorption lines likely associated with ionized Iron K-shell transitions in the X-ray spectra of many Active Galactic Nuclei (AGN) suggests the presence of a highly ionized gas outflowing with mildly relativistic velocities (0.03c-0.6c), named Ultra-Fast Outflow (UFO). Within the SUBWAYS project we characterized these winds starting from a sample of 22 radio-quiet qua… ▽ More The detection of blue-shifted absorption lines likely associated with ionized Iron K-shell transitions in the X-ray spectra of many Active Galactic Nuclei (AGN) suggests the presence of a highly ionized gas outflowing with mildly relativistic velocities (0.03c-0.6c), named Ultra-Fast Outflow (UFO). Within the SUBWAYS project we characterized these winds starting from a sample of 22 radio-quiet quasars at 0.1 < z < 0.4, and compared the results with similar studies in the literature on samples of 42 local radio-quiet Seyfert galaxies and 14 high redshift radio-quiet quasars. The scope of our work is a statistical study of UFO parameters and incidence, considering key physical properties of the sources, e.g. supermassive black hole (SMBH) mass, bolometric luminosity, accretion rates and Spectral Energy Distribution, with the aim of gaining new insights into the UFO launching mechanisms. We find indications that highly luminous AGN with steeper X-ray/UV ratio, are more likely to host UFO. The presence of UFO is not significantly related to any other AGN property in our sample. These findings suggest that the UFO phenomenon may be transient. Focusing on AGN with UFO, other important results are: (1) faster UFO have larger ionization parameters and column densities; (2) X-ray radiation plays a more crucial role in driving highly ionized winds compared to UV; (3) the correlation between outflow velocity and luminosity is significantly flatter than what expected for radiatively driven winds; (4) more massive BH experience higher wind mass-losses, suppressing accretion of matter onto the BH; (5) the UFO launching radius is positively correlated with the Eddington ratio. Furthermore, our analysis suggest the involvement of multiple launching mechanisms, including radiation pressure and magneto-hydrodynamic processes, rather than pointing to a single, universally applicable mechanism. △ Less

Submitted 11 April, 2024; v1 submitted 14 March, 2024; originally announced March 2024.

Comments: 52 pages, 154 figures and 7 tables. Accepted in Astronomy & Astrophysics, in press

Journal ref: A&A 687, A235 (2024)

Resolving Dual Active Galactic Nuclei with ~100 pc separation in MCG-03-34-64

Authors: Anna Trindade Falcao, T. J. Turner, S. B. Kraemer, V. Braito, J. Reeves, H. R. Schmitt, L. Feuillet

Abstract: We report the serendipitous multiwavelength discovery of a candidate dual black hole system with a separation of ~100 pc, in the gas-rich luminous infrared galaxy MCG-03-34-64 (z=0.016). Hubble Space Telescope/ACS observations show three distinct optical centroids in the [O III] narrow-band and F814W images. Subsequent analysis of Chandra/ACIS data shows two spatially-resolved peaks of equal inten… ▽ More We report the serendipitous multiwavelength discovery of a candidate dual black hole system with a separation of ~100 pc, in the gas-rich luminous infrared galaxy MCG-03-34-64 (z=0.016). Hubble Space Telescope/ACS observations show three distinct optical centroids in the [O III] narrow-band and F814W images. Subsequent analysis of Chandra/ACIS data shows two spatially-resolved peaks of equal intensity in the neutral Fe Ka (6.2-6.6 keV) band, while high-resolution radio continuum observations with the Very Large Array at 8.46 GHz (3.6 cm band) show two spatially-coincident radio peaks. Fast shocks as the ionizing source seem unlikely, given the energies required for production of Fe Ka. If confirmed, the separation of ~100 pc would represent the closest dual AGN reported to date with spatially-resolved, multiwavelength observations. △ Less

Submitted 20 July, 2024; v1 submitted 12 March, 2024; originally announced March 2024.

Comments: Submitted to ApJ

How causal inference concepts can guide research into the effects of climate on infectious diseases

Authors: Laura Andrea Barrero Guevara, Sarah C Kramer, Tobias Kurth, Matthieu Domenech de Cellès

Abstract: A pressing question resulting from global warming is how infectious diseases will be affected by climate change. Answering this question requires research into the effects of weather on the population dynamics of transmission and infection; elucidating these effects, however, has proven difficult due to the challenges of assessing causality from the predominantly observational data available in ep… ▽ More A pressing question resulting from global warming is how infectious diseases will be affected by climate change. Answering this question requires research into the effects of weather on the population dynamics of transmission and infection; elucidating these effects, however, has proven difficult due to the challenges of assessing causality from the predominantly observational data available in epidemiological research. Here, we show how concepts from causal inference -- the sub-field of statistics aiming at inferring causality from data -- can guide that research. Through a series of case studies, we illustrate how such concepts can help assess study design and strategically choose a study's location, evaluate and reduce the risk of bias, and interpret the multifaceted effects of meteorological variables on transmission. More broadly, we argue that interdisciplinary approaches based on explicit causal frameworks are crucial for reliably estimating the effect of weather and accurately predicting the consequences of climate change. △ Less

Submitted 19 February, 2024; originally announced February 2024.

Amplifying Exploration in Monte-Carlo Tree Search by Focusing on the Unknown

Authors: Cedric Derstroff, Jannis Brugger, Jannis Blüml, Mira Mezini, Stefan Kramer, Kristian Kersting

Abstract: Monte-Carlo tree search (MCTS) is an effective anytime algorithm with a vast amount of applications. It strategically allocates computational resources to focus on promising segments of the search tree, making it a very attractive search algorithm in large search spaces. However, it often expends its limited resources on reevaluating previously explored regions when they remain the most promising… ▽ More Monte-Carlo tree search (MCTS) is an effective anytime algorithm with a vast amount of applications. It strategically allocates computational resources to focus on promising segments of the search tree, making it a very attractive search algorithm in large search spaces. However, it often expends its limited resources on reevaluating previously explored regions when they remain the most promising path. Our proposed methodology, denoted as AmEx-MCTS, solves this problem by introducing a novel MCTS formulation. Central to AmEx-MCTS is the decoupling of value updates, visit count updates, and the selected path during the tree search, thereby enabling the exclusion of already explored subtrees or leaves. This segregation preserves the utility of visit counts for both exploration-exploitation balancing and quality metrics within MCTS. The resultant augmentation facilitates in a considerably broader search using identical computational resources, preserving the essential characteristics of MCTS. The expanded coverage not only yields more precise estimations but also proves instrumental in larger and more complex problems. Our empirical evaluation demonstrates the superior performance of AmEx-MCTS, surpassing classical MCTS and related approaches by a substantial margin. △ Less

Submitted 13 February, 2024; originally announced February 2024.

Comments: 10 pages, 7 figures

Deconstructing the spin susceptibility of a cuprate superconductor

Authors: R. Zhou, I. Vinograd, M. Hirata, T. Wu, H. Mayaffre, S. Krämer, W. N. Hardy, R. Liang, D. A. Bonn, T. Loew, J. Porras, B. Keimer, M. -H. Julien

Abstract: A major obstacle to understanding high-Tc cuprates is that superconductivity precludes observing normal-state properties at low temperatures. One prime example is the normal-state spin susceptibility: although its decrease upon cooling far above Tc typifies pseudogap behavior, its behavior at low temperatures is generally unknown. Here, our measurements in high magnetic fields expose the spin susc… ▽ More A major obstacle to understanding high-Tc cuprates is that superconductivity precludes observing normal-state properties at low temperatures. One prime example is the normal-state spin susceptibility: although its decrease upon cooling far above Tc typifies pseudogap behavior, its behavior at low temperatures is generally unknown. Here, our measurements in high magnetic fields expose the spin susceptibility of YBa2Cu3Oy down to low temperatures. Even though superconductivity is suppressed by the field, we uncover two thermally-activated contributions alongside a residual susceptibility at T=0 due to gapless excitations. We relate these two distinct gaps to short-range charge-density waves and to the formation of spin singlets similar to those found in certain quantum spin systems. These phenomena thus collectively contribute to the pseudogap in the spin susceptibility at low temperature, supplementing short-lived antiferromagnetism known to initiate pseudogap behavior at high temperatures. We therefore propose that the pseudogap should be regarded as a composite property. △ Less

Submitted 4 February, 2024; originally announced February 2024.

Classifying Intermediate Redshift Galaxies in SDSS: Alternative Diagnostic Diagrams

Authors: Léa Feuillet, Marcio Meléndez, Steve Kraemer, Henrique Schmitt, Travis Fischer, James Reeves

Abstract: We select a sample of 1,437 active galactic nuclei (AGN) from the catalog of the Sloan Digital Sky Survey (SDSS) galaxy properties from the Portsmouth group by detection of the high-ionization [Ne V] 3426 Å emission line. We compare the fluxes of [Ne III] 3869 Å, [O III] 5007 Å, [O II] 3726, 3728 Å, and [O I] 6300 Å to that of [Ne V]. All four lines show a strong linear correlation with [Ne V], al… ▽ More We select a sample of 1,437 active galactic nuclei (AGN) from the catalog of the Sloan Digital Sky Survey (SDSS) galaxy properties from the Portsmouth group by detection of the high-ionization [Ne V] 3426 Å emission line. We compare the fluxes of [Ne III] 3869 Å, [O III] 5007 Å, [O II] 3726, 3728 Å, and [O I] 6300 Å to that of [Ne V]. All four lines show a strong linear correlation with [Ne V], although lines from ions with lower ionization potentials have a lower correlation coefficient. We investigate the use of two forbidden-line ratio (FLR) diagnostic diagrams that do not rely on H$α$ in order to classify high redshift galaxies. These use the [Ne III]/[O II] line ratio plotted against [O III]/[O I] and [O III]/[O II] respectively. We use photo-ionization modeling to characterize the behavior of the narrow-line region in AGN and star-forming regions and test the validity of our diagnostic diagrams. We also use a luminosity cutoff of log L[OIII] [erg/s] = 42, which lowers the contamination of the AGN region by star-forming galaxies down to 10% but does not remove Green Pea and Purple Grape galaxies from the AGN region. We also investigate the OHNO diagram which uses [Ne III]/[O II] plotted against [O III]/H$β$. Using our new diagnostic diagrams, we are able to reliably classify AGN up to a redshift of z $\leq$ 1.06, and add more than 822 new AGN to the [Ne V]-selected AGN sample. △ Less

Submitted 9 April, 2024; v1 submitted 28 December, 2023; originally announced December 2023.

Comments: Published February 12 2024 in ApJ, 15 pages, 10 figures

Journal ref: ApJ 962 104 (2024)

Peer Learning: Learning Complex Policies in Groups from Scratch via Action Recommendations

Authors: Cedric Derstroff, Mattia Cerrato, Jannis Brugger, Jan Peters, Stefan Kramer

Abstract: Peer learning is a novel high-level reinforcement learning framework for agents learning in groups. While standard reinforcement learning trains an individual agent in trial-and-error fashion, all on its own, peer learning addresses a related setting in which a group of agents, i.e., peers, learns to master a task simultaneously together from scratch. Peers are allowed to communicate only about th… ▽ More Peer learning is a novel high-level reinforcement learning framework for agents learning in groups. While standard reinforcement learning trains an individual agent in trial-and-error fashion, all on its own, peer learning addresses a related setting in which a group of agents, i.e., peers, learns to master a task simultaneously together from scratch. Peers are allowed to communicate only about their own states and actions recommended by others: "What would you do in my situation?". Our motivation is to study the learning behavior of these agents. We formalize the teacher selection process in the action advice setting as a multi-armed bandit problem and therefore highlight the need for exploration. Eventually, we analyze the learning behavior of the peers and observe their ability to rank the agents' performance within the study group and understand which agents give reliable advice. Further, we compare peer learning with single agent learning and a state-of-the-art action advice baseline. We show that peer learning is able to outperform single-agent learning and the baseline in several challenging discrete and continuous OpenAI Gym domains. Doing so, we also show that within such a framework complex policies from action recommendations beyond discrete action spaces can evolve. △ Less

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

Comments: 9 pages, 7 figures, AAAI-24

Journal ref: AAAI, vol. 38, no. 10, pp. 11766-11774, Mar. 2024

A Generalized Bandsplit Neural Network for Cinematic Audio Source Separation

Authors: Karn N. Watcharasupat, Chih-Wei Wu, Yiwei Ding, Iroro Orife, Aaron J. Hipple, Phillip A. Williams, Scott Kramer, Alexander Lerch, William Wolcott

Abstract: Cinematic audio source separation is a relatively new subtask of audio source separation, with the aim of extracting the dialogue, music, and effects stems from their mixture. In this work, we developed a model generalizing the Bandsplit RNN for any complete or overcomplete partitions of the frequency axis. Psychoacoustically motivated frequency scales were used to inform the band definitions whic… ▽ More Cinematic audio source separation is a relatively new subtask of audio source separation, with the aim of extracting the dialogue, music, and effects stems from their mixture. In this work, we developed a model generalizing the Bandsplit RNN for any complete or overcomplete partitions of the frequency axis. Psychoacoustically motivated frequency scales were used to inform the band definitions which are now defined with redundancy for more reliable feature extraction. A loss function motivated by the signal-to-noise ratio and the sparsity-promoting property of the 1-norm was proposed. We additionally exploit the information-sharing property of a common-encoder setup to reduce computational complexity during both training and inference, improve separation performance for hard-to-generalize classes of sounds, and allow flexibility during inference time with detachable decoders. Our best model sets the state of the art on the Divide and Remaster dataset with performance above the ideal ratio mask for the dialogue stem. △ Less

Submitted 1 December, 2023; v1 submitted 5 September, 2023; originally announced September 2023.

Comments: Accepted to the IEEE Open Journal of Signal Processing (ICASSP 2024 Track)

Longitudinal spin fluctuations driving field-reinforced superconductivity in UTe$_2$

Authors: Yo Tokunaga, Hironori Sakai, Shinsaku Kambe, Petr Opletal, Yoshifumi Tokiwa, Yoshinori Haga, Shunsaku Kitagawa, Kenji Ishida, Dai Aoki, Georg Knebel, Gerard Lapertot, Steffen Krämer, Mladen Horvatić

Abstract: Our measurements of $^{125}$Te NMR relaxations reveal an enhancement of electronic spin fluctuations above $μ_0H^*\sim15$ T, leading to their divergence in the vicinity of the metamagnetic transition at $μ_0H_m\approx35$ T, below which field-reinforced superconductivity appears when a magnetic field ($H$) is applied along the crystallographic $b$ axis. The NMR data evidence that these fluctuations… ▽ More Our measurements of $^{125}$Te NMR relaxations reveal an enhancement of electronic spin fluctuations above $μ_0H^*\sim15$ T, leading to their divergence in the vicinity of the metamagnetic transition at $μ_0H_m\approx35$ T, below which field-reinforced superconductivity appears when a magnetic field ($H$) is applied along the crystallographic $b$ axis. The NMR data evidence that these fluctuations are dominantly longitudinal, providing a key to understanding the peculiar superconducting phase diagram in $H\|b$, where such fluctuations enhance the pairing interactions. △ Less

Submitted 30 November, 2023; v1 submitted 20 July, 2023; originally announced July 2023.

Comments: 7 pages, 4 figures

Journal ref: Phys. Rev. Lett. 131, 226503 (2023)

No Small Scale Radio Jets Here: Multi-Epoch Observations of Radio Continuum Structures in NGC 1068 with the VLBA

Authors: Travis C. Fischer, Megan C. Johnson, Nathan J. Secrest, D. Michael Crenshaw, Steven B. Kraemer

Abstract: We present recent Very Long Baseline Array (VLBA) 5 GHz radio observations of the nearby, luminous Seyfert 2 galaxy NGC 1068 for comparison to similar VLBA observations made on 1997 April 26. By cross-correlating the positions of emitting regions across both epochs, we find that spatially-resolved extra-nuclear radio knots in this system have sub-relativistic transverse speeds (v < 0.1c). We discu… ▽ More We present recent Very Long Baseline Array (VLBA) 5 GHz radio observations of the nearby, luminous Seyfert 2 galaxy NGC 1068 for comparison to similar VLBA observations made on 1997 April 26. By cross-correlating the positions of emitting regions across both epochs, we find that spatially-resolved extra-nuclear radio knots in this system have sub-relativistic transverse speeds (v < 0.1c). We discuss sources of the observed knots and how the radio emission relates to additional phases of gas in the central ~150 pcs of this system. We suggest that the most likely explanation for the observed emission is synchrotron radiation formed by shocked host media via interactions between AGN winds and the host environment. △ Less

Submitted 26 June, 2023; originally announced June 2023.

Comments: 13 pages, 4 figures, accepted for publication in ApJ

Automated Scientific Discovery: From Equation Discovery to Autonomous Discovery Systems

Authors: Stefan Kramer, Mattia Cerrato, Sašo Džeroski, Ross King

Abstract: The paper surveys automated scientific discovery, from equation discovery and symbolic regression to autonomous discovery systems and agents. It discusses the individual approaches from a "big picture" perspective and in context, but also discusses open issues and recent topics like the various roles of deep neural networks in this area, aiding in the discovery of human-interpretable knowledge. Fu… ▽ More The paper surveys automated scientific discovery, from equation discovery and symbolic regression to autonomous discovery systems and agents. It discusses the individual approaches from a "big picture" perspective and in context, but also discusses open issues and recent topics like the various roles of deep neural networks in this area, aiding in the discovery of human-interpretable knowledge. Further, we will present closed-loop scientific discovery systems, starting with the pioneering work on the Adam system up to current efforts in fields from material science to astronomy. Finally, we will elaborate on autonomy from a machine learning perspective, but also in analogy to the autonomy levels in autonomous driving. The maximal level, level five, is defined to require no human intervention at all in the production of scientific knowledge. Achieving this is one step towards solving the Nobel Turing Grand Challenge to develop AI Scientists: AI systems capable of making Nobel-quality scientific discoveries highly autonomously at a level comparable, and possibly superior, to the best human scientists by 2050. △ Less

Submitted 3 May, 2023; originally announced May 2023.

Supermassive Black Hole Winds in X-rays: SUBWAYS. II. HST UV spectroscopy of winds at intermediate redshifts

Authors: M. Mehdipour, G. A. Kriss, M. Brusa, G. A. Matzeu, M. Gaspari, S. B. Kraemer, S. Mathur, E. Behar, S. Bianchi, M. Cappi, G. Chartas, E. Costantini, G. Cresci, M. Dadina, B. De Marco, A. De Rosa, J. P. Dunn, V. E. Gianolli, M. Giustini, J. S. Kaastra, A. R. King, Y. Krongold, F. La Franca, G. Lanzuisi, A. L. Longinotti , et al. (13 additional authors not shown)

Abstract: We present a UV spectroscopic study of ionized outflows in 21 active galactic nuclei (AGN), observed with the HST. The targets of the SUBWAYS sample were selected with the aim to probe the parameter space of the underexplored AGN between the local Seyfert galaxies and the luminous quasars at high redshifts. Our targets, spanning redshifts of 0.1-0.4 and bolometric luminosities (L_bol) of 10^45-10^… ▽ More We present a UV spectroscopic study of ionized outflows in 21 active galactic nuclei (AGN), observed with the HST. The targets of the SUBWAYS sample were selected with the aim to probe the parameter space of the underexplored AGN between the local Seyfert galaxies and the luminous quasars at high redshifts. Our targets, spanning redshifts of 0.1-0.4 and bolometric luminosities (L_bol) of 10^45-10^46 erg/s, have been observed with a large multi-wavelength campaign. Here, we model the UV spectra and look for different types of AGN outflows. We find that 60% of our targets show a presence of outflowing H I absorption, while 40% exhibit ionized outflows seen as absorption by either C IV, N V, or O VI. This is comparable to the occurrence of ionized outflows seen in the local Seyfert galaxies. All UV absorption lines in the sample are relatively narrow, with outflow velocities reaching up to -3300 km/s. We did not detect any UV counterparts to the X-ray ultra-fast outflows (UFOs), most likely due to their being too highly ionized. However, all SUBWAYS targets with an X-ray UFO demonstrate the presence of UV outflows at lower velocities. We find significant correlations between the column density (N) of the UV ions and L_bol of the AGN, with N of H I decreasing with L_bol, while N of O VI is increasing with L_bol. This is likely to be a photoionization effect, where toward higher AGN luminosities, the wind becomes more ionized, resulting in less absorption by neutral or low-ionization ions and more absorption by high-ionization ions. In addition, we find that N of the UV ions decreases as their outflow velocity increases. This may be explained by a mechanical power that is evacuating the UV-absorbing medium. Our observed relations are consistent with multiphase AGN feeding and feedback simulations indicating that a combination of both radiative and mechanical processes are in play. △ Less

Submitted 6 December, 2022; originally announced December 2022.

Comments: Accepted for publication in Astronomy & Astrophysics (A&A), 22 pages, 12 figures

Supermassive Black Hole Winds in X-rays -- SUBWAYS. I. Ultra-fast outflows in QSOs beyond the local Universe

Authors: G. A. Matzeu, M. Brusa, G. Lanzuisi, M. Dadina, S. Bianchi, G. Kriss, M. Mehdipour, E. Nardini, G. Chartas, R. Middei, E. Piconcelli, V. Gianolli, A. Comastri, A. L. Longinotti, Y. Krongold, F. Ricci, P. O. Petrucci, F. Tombesi, A. Luminari, L. Zappacosta, G. Miniutti, M. Gaspari, E. Behar, M. Bischetti, S. Mathur , et al. (26 additional authors not shown)

Abstract: We present a new X-ray spectroscopic study of $22$ luminous ($2\times10^{45}\lesssim L_{\rm bol}\rm /erg\,s^{-1} \lesssim 2\times10^{46}$) active galactic nuclei (AGNs) at intermediate-redshift ($0.1 \lesssim z \lesssim 0.4$), as part of the SUpermassive Black hole Winds in the x-rAYS (SUBWAYS) sample, mostly composed of quasars (QSOs) and type\,1 AGN. Here, 17 targets were observed with \textit{X… ▽ More We present a new X-ray spectroscopic study of $22$ luminous ($2\times10^{45}\lesssim L_{\rm bol}\rm /erg\,s^{-1} \lesssim 2\times10^{46}$) active galactic nuclei (AGNs) at intermediate-redshift ($0.1 \lesssim z \lesssim 0.4$), as part of the SUpermassive Black hole Winds in the x-rAYS (SUBWAYS) sample, mostly composed of quasars (QSOs) and type\,1 AGN. Here, 17 targets were observed with \textit{XMM-Newton} between 2019--2020 and the remaining 5 are from previous observations. The aim of this large campaign ($1.45\,\rm Ms$ duration) is to characterise the various manifestations of winds in the X-rays driven from supermassive black holes in AGN. In this paper we focus on the search and characterization of ultra-fast outflows (UFOs), which are typically detected through blueshifted absorption troughs in the Fe\,K band ($E>7\,\rm keV$). By following Monte Carlo procedures, we confirm the detection of absorption lines corresponding to highly ionised iron (e.g., Fe\,\textsc{xxv}\,H$α$, Fe\,\textsc{xxvi}\,Ly$α$) in 7/22 sources at the $\gtrsim95\%$ confidence level (for each individual line). The global combined probability of such absorption features in the sample is $>99.9\%$. The SUBWAYS campaign extends at higher luminosity and redshifts than previous local studies on Seyferts, obtained using \xmm and \suzaku observations. We find a UFO detection fraction of $\sim30\%$ on the total sample that is in agreement with the previous findings. This work independently provides further support for the existence of highly-ionised matter propagating at mildly relativistic speed ($\gtrsim0.1c$) in a considerable fraction of AGN over a broad range of luminosities, which is expected to play a key role in the self-regulated AGN feeding-feedback cycle, as also supported by hydrodynamical multiphase simulations. △ Less

Submitted 6 December, 2022; originally announced December 2022.

Comments: 30 pages, 10 figures, accepted for publication in A&A

Investigating the Narrow Line Region Dynamics in Nearby Active Galaxies

Authors: Beena Meena, D. Michael Crenshaw, Henrique R. Schmitt, Mitchell Revalski, Zo Chapman, Travis C. Fischer, Steven B. Kraemer, Justin H. Robinson, Julia Falcone, Garrett E. Polack

Abstract: We present dynamical models of the narrow line region (NLR) outflows in the nearby Seyfert galaxies Mrk 3, Mrk 78, NGC 1068, and NGC 4151 using observations from the Hubble Space Telescope and Apache Point Observatory. We employ long-slit spectroscopy to map the spatially-resolved outflow and rotational velocities of the ionized gas. We also perform surface brightness decompositions of host galaxy… ▽ More We present dynamical models of the narrow line region (NLR) outflows in the nearby Seyfert galaxies Mrk 3, Mrk 78, NGC 1068, and NGC 4151 using observations from the Hubble Space Telescope and Apache Point Observatory. We employ long-slit spectroscopy to map the spatially-resolved outflow and rotational velocities of the ionized gas. We also perform surface brightness decompositions of host galaxy images to constrain the enclosed stellar mass distributions as functions of distance from the supermassive black holes (SMBHs). Assuming that the NLR gas is accelerated by AGN radiation pressure, and subsequently decelerated by the host galaxy and SMBH gravitational potentials, we derive outflow velocity profiles where the gas is launched in situ at multiple distances from the SMBH. We find a strong correlation between the turnover (from acceleration to deceleration) radii from our models, with the turnovers seen in the observed velocities and spatially-resolved mass outflow rates for the AGN with bolometric luminosities $>$ 10$^{44}$ erg sec$^{-1}$. This consistency indicates that radiation pressure is the dominant driving mechanism behind the NLR outflows in these moderate-luminosity AGN, with a force multiplier $\sim$500 yielding the best agreement between the modeled and observed turnover radii. However, in Meena2021 we found that this trend may not hold at lower luminosities, where our modeled turnover distance for NGC 4051 is much smaller than in the observed kinematics. This result may indicate that either additional force(s) are responsible for accelerating the NLR outflows in low-luminosity AGN, or higher spatial resolution observations are required to quantify their turnover radii. △ Less

Submitted 5 December, 2022; originally announced December 2022.

Comments: 36 pages, 10 Figures, 5 Tables, 7 supplementary Figures, Accepted for publication in ApJ

Neural RELAGGS

Authors: Lukas Pensel, Stefan Kramer

Abstract: Multi-relational databases are the basis of most consolidated data collections in science and industry today. Most learning and mining algorithms, however, require data to be represented in a propositional form. While there is a variety of specialized machine learning algorithms that can operate directly on multi-relational data sets, propositionalization algorithms transform multi-relational data… ▽ More Multi-relational databases are the basis of most consolidated data collections in science and industry today. Most learning and mining algorithms, however, require data to be represented in a propositional form. While there is a variety of specialized machine learning algorithms that can operate directly on multi-relational data sets, propositionalization algorithms transform multi-relational databases into propositional data sets, thereby allowing the application of traditional machine learning and data mining algorithms without their modification. One prominent propositionalization algorithm is RELAGGS by Krogel and Wrobel, which transforms the data by nested aggregations. We propose a new neural network based algorithm in the spirit of RELAGGS that employs trainable composite aggregate functions instead of the static aggregate functions used in the original approach. In this way, we can jointly train the propositionalization with the prediction model, or, alternatively, use the learned aggegrations as embeddings in other algorithms. We demonstrate the increased predictive performance by comparing N-RELAGGS with RELAGGS and multiple other state-of-the-art algorithms. △ Less

Submitted 4 November, 2022; originally announced November 2022.

Comments: Submitted to Machine Learning Journal

A UFO Seen Edge-On: Resolving Ultrafast Outflow Emission on $\sim$200-pc Scales with $Chandra$ in the Active Nucleus of Mrk 34

Authors: W. Peter Maksym, Martin Elvis, Giuseppina Fabbiano, Anna Trindade-Falcão, Steven B. Kraemer, Travis C. Fischer, D. Michael Crenshaw, Thaisa Storchi-Bergmann

Abstract: We present $Chandra$ ACIS imaging spectroscopy of the nucleus of the Seyfert 2 Galaxy Mrk 34. We identify spatially and spectrally resolved features in the band that includes Fe K$α$, Fe XXV and Fe XXVI. These features indicate high-velocity ($\gtrsim15,000\,\rm{km\,s}^{-1}$ line-of-sight) material separated spanning $\sim$0.5 arcsec, within $\sim200$ pc of the nucleus. This outflow could have dep… ▽ More We present $Chandra$ ACIS imaging spectroscopy of the nucleus of the Seyfert 2 Galaxy Mrk 34. We identify spatially and spectrally resolved features in the band that includes Fe K$α$, Fe XXV and Fe XXVI. These features indicate high-velocity ($\gtrsim15,000\,\rm{km\,s}^{-1}$ line-of-sight) material separated spanning $\sim$0.5 arcsec, within $\sim200$ pc of the nucleus. This outflow could have deprojected velocities $\sim12-28\times$ greater than the [O III] emitting outflows, and could potentially dominate the kinetic power in the outflow. This emission may point to the origins of the optical and X-ray winds observed at larger radii, and could indicate a link between ultra-fast outflows and AGN feedback on $\gtrsim$kpc scales. △ Less

Submitted 20 March, 2023; v1 submitted 28 September, 2022; originally announced September 2022.

Comments: 13 pages, 7 figures. Resubmitted to ApJ with revisions

Design of experiments for the calibration of history-dependent models via deep reinforcement learning and an enhanced Kalman filter

Authors: Ruben Villarreal, Nikolaos N. Vlassis, Nhon N. Phan, Tommie A. Catanach, Reese E. Jones, Nathaniel A. Trask, Sharlotte L. B. Kramer, WaiChing Sun

Abstract: Experimental data is costly to obtain, which makes it difficult to calibrate complex models. For many models an experimental design that produces the best calibration given a limited experimental budget is not obvious. This paper introduces a deep reinforcement learning (RL) algorithm for design of experiments that maximizes the information gain measured by Kullback-Leibler (KL) divergence obtaine… ▽ More Experimental data is costly to obtain, which makes it difficult to calibrate complex models. For many models an experimental design that produces the best calibration given a limited experimental budget is not obvious. This paper introduces a deep reinforcement learning (RL) algorithm for design of experiments that maximizes the information gain measured by Kullback-Leibler (KL) divergence obtained via the Kalman filter (KF). This combination enables experimental design for rapid online experiments where traditional methods are too costly. We formulate possible configurations of experiments as a decision tree and a Markov decision process (MDP), where a finite choice of actions is available at each incremental step. Once an action is taken, a variety of measurements are used to update the state of the experiment. This new data leads to a Bayesian update of the parameters by the KF, which is used to enhance the state representation. In contrast to the Nash-Sutcliffe efficiency (NSE) index, which requires additional sampling to test hypotheses for forward predictions, the KF can lower the cost of experiments by directly estimating the values of new data acquired through additional actions. In this work our applications focus on mechanical testing of materials. Numerical experiments with complex, history-dependent models are used to verify the implementation and benchmark the performance of the RL-designed experiments. △ Less

Submitted 26 September, 2022; originally announced September 2022.

Comments: 40 pages, 20 figures

Observation of the radiative decay of the ${}^{229}\mathrm{Th}$ nuclear clock isomer

Authors: Sandro Kraemer, Janni Moens, Michail Athanasakis-Kaklamanakis, Silvia Bara, Kjeld Beeks, Premaditya Chhetri, Katerina Chrysalidis, Arno Claessens, Thomas E. Cocolios, João M. Correia, Hilde De Witte, Rafael Ferrer, Sarina Geldhof, Reinhard Heinke, Niyusha Hosseini, Mark Huyse, Ulli Köster, Yuri Kudryavtsev, Mustapha Laatiaoui, Razvan Lica, Goele Magchiels, Vladimir Manea, Clement Merckling, Lino M. C. Pereira, Sebastian Raeder , et al. (10 additional authors not shown)

Abstract: The nucleus of the radioisotope thorium-229 (${}^{229}$Th) features an isomer with an exceptionally low excitation energy that enables direct laser manipulation of nuclear states. For this reason, it is a leading candidate for use in next-generation optical clocks. This nuclear clock will be a unique tool, amongst others, for tests of fundamental physics. While first indirect experimental evidence… ▽ More The nucleus of the radioisotope thorium-229 (${}^{229}$Th) features an isomer with an exceptionally low excitation energy that enables direct laser manipulation of nuclear states. For this reason, it is a leading candidate for use in next-generation optical clocks. This nuclear clock will be a unique tool, amongst others, for tests of fundamental physics. While first indirect experimental evidence for the existence of such an extraordinary nuclear state is significantly older, the proof of existence has been delivered only recently by observing the isomer's electron conversion decay and its hyperfine structure in a laser spectroscopy study, revealing information on the isomer's excitation energy, nuclear spin and electromagnetic moments. Further studies reported the electron conversion lifetime and refined the isomer's energy. In spite of recent progress, the isomer's radiative decay, a key ingredient for the development of a nuclear clock, remained unobserved. In this Letter, we report the detection of the radiative decay of this low-energy isomer in thorium-229 (${}^{229\mathrm{m}}$Th). By performing vacuum-ultraviolet spectroscopy of ${}^{229\mathrm{m}}$Th incorporated into large-bandgap CaF${}_2$ and MgF${}_2$ crystals at the ISOLDE facility at CERN, the photon vacuum wavelength of the isomer's decay is measured as 148.71(42) nm, corresponding to an excitation energy of 8.338(24) eV. This value is in agreement with recent measurements, and decreases the uncertainty by a factor of seven. The half-life of ${}^{229\mathrm{m}}$Th embedded in MgF${}_2$ is determined to be 670(102) s. The observation of the radiative decay in a large-bandgap crystal has important consequences for the design of a future nuclear clock and the improved uncertainty of the energy eases the search for direct laser excitation of the atomic nucleus. △ Less

Submitted 21 September, 2022; originally announced September 2022.

A Fair Experimental Comparison of Neural Network Architectures for Latent Representations of Multi-Omics for Drug Response Prediction

Authors: Tony Hauptmann, Stefan Kramer

Abstract: Recent years have seen a surge of novel neural network architectures for the integration of multi-omics data for prediction. Most of the architectures include either encoders alone or encoders and decoders, i.e., autoencoders of various sorts, to transform multi-omics data into latent representations. One important parameter is the depth of integration: the point at which the latent representation… ▽ More Recent years have seen a surge of novel neural network architectures for the integration of multi-omics data for prediction. Most of the architectures include either encoders alone or encoders and decoders, i.e., autoencoders of various sorts, to transform multi-omics data into latent representations. One important parameter is the depth of integration: the point at which the latent representations are computed or merged, which can be either early, intermediate, or late. The literature on integration methods is growing steadily, however, close to nothing is known about the relative performance of these methods under fair experimental conditions and under consideration of different use cases. We developed a comparison framework that trains and optimizes multi-omics integration methods under equal conditions. We incorporated early integration and four recently published deep learning methods: MOLI, Super.FELT, OmiEmbed, and MOMA. Further, we devised a novel method, Omics Stacking, that combines the advantages of intermediate and late integration. Experiments were conducted on a public drug response data set with multiple omics data (somatic point mutations, somatic copy number profiles and gene expression profiles) that was obtained from cell lines, patient-derived xenografts, and patient samples. Our experiments confirmed that early integration has the lowest predictive performance. Overall, architectures that integrate triplet loss achieved the best results. Statistical differences can, overall, rarely be observed, however, in terms of the average ranks of methods, Super.FELT is consistently performing best in a cross-validation setting and Omics Stacking best in an external test set setting. The source code of all experiments is available under \url{https://github.com/kramerlab/Multi-Omics_analysis} △ Less

Submitted 31 August, 2022; originally announced August 2022.

Invariant Representations with Stochastically Quantized Neural Networks

Authors: Mattia Cerrato, Marius Köppel, Roberto Esposito, Stefan Kramer

Abstract: Representation learning algorithms offer the opportunity to learn invariant representations of the input data with regard to nuisance factors. Many authors have leveraged such strategies to learn fair representations, i.e., vectors where information about sensitive attributes is removed. These methods are attractive as they may be interpreted as minimizing the mutual information between a neural l… ▽ More Representation learning algorithms offer the opportunity to learn invariant representations of the input data with regard to nuisance factors. Many authors have leveraged such strategies to learn fair representations, i.e., vectors where information about sensitive attributes is removed. These methods are attractive as they may be interpreted as minimizing the mutual information between a neural layer's activations and a sensitive attribute. However, the theoretical grounding of such methods relies either on the computation of infinitely accurate adversaries or on minimizing a variational upper bound of a mutual information estimate. In this paper, we propose a methodology for direct computation of the mutual information between a neural layer and a sensitive attribute. We employ stochastically-activated binary neural networks, which lets us treat neurons as random variables. We are then able to compute (not bound) the mutual information between a layer and a sensitive attribute and use this information as a regularization factor during gradient descent. We show that this method compares favorably with the state of the art in fair representation learning and that the learned representations display a higher level of invariance compared to full-precision neural networks. △ Less

Submitted 2 December, 2022; v1 submitted 4 August, 2022; originally announced August 2022.

Comments: To appear in AAAI23

Competition between spin ordering and superconductivity near the pseudogap boundary in La2-xSrxCuO4: insights from NMR

Authors: I. Vinograd, R. Zhou, H. Mayaffre, S. Krämer, S. K. Ramakrishna, A. P. Reyes, T. Kurosawa, N. Momono, M. Oda, S. Komiya, S. Ono, M. Horio, J. Chang, M. -H. Julien

Abstract: When superconductivity is suppressed by high magnetic fields in La2-xSrxCuO4, striped antiferromagnetic (AFM) order becomes the magnetic ground state of the entire pseudogap regime, up to its end at the doping p* [M. Frachet, I. Vinograd et al., Nat. Phys. 16, 1064 (2020)]. Glass-like freezing of this state is detected in 139La NMR measurements of the spin-lattice relaxation rate 1/T1. Here, we pr… ▽ More When superconductivity is suppressed by high magnetic fields in La2-xSrxCuO4, striped antiferromagnetic (AFM) order becomes the magnetic ground state of the entire pseudogap regime, up to its end at the doping p* [M. Frachet, I. Vinograd et al., Nat. Phys. 16, 1064 (2020)]. Glass-like freezing of this state is detected in 139La NMR measurements of the spin-lattice relaxation rate 1/T1. Here, we present a quantitative analysis of 1/T1 data in the hole-doping range p=x=0.12-0.171, based on the Bloembergen-Purcell-Pound (BPP) theory, modified to include statistical distribution of parameters arising from strong spatial inhomogeneity. We observe spin fluctuations to slow down at temperatures T near the onset of static charge order and, overall, the effect of the field B may be seen as equivalent to strengthening stripe order by approaching p=0.12 doping. In details however, our analysis reveals significant departure from usual field-induced magnetic transitions. The continuous growth of the amplitude of the fluctuating moment with increasing B suggests a nearly-critical state in the B->0 limit, with very weak quasi-static moments possibly confined in small areas like vortex cores. Further, the nucleation of spin order in the vortex cores is shown to account quantitatively for both the value and the p dependence of a field scale characterizing bulk spin freezing. The correlation time of the fluctuating moment appears to depend exponentially on B/T (over the investigated range). This explains the timescale dependence of various experimental manifestations, including why, for transport measurements, the AFM moments may be considered static over a considerable range of B and T. These results make the high-field magnetic ground state up to p* an integral part of the discussion on putative quantum criticality. △ Less

Submitted 10 July, 2022; originally announced July 2022.

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

Dynamic neighbors: a proposal of a tool to characterize phase transitions

Authors: L. Aramis de Icaza Astiz, Atahualpa S. Kraemer, Gerardo Odriozola, Mariano López de Haro

Abstract: For molecular dynamics simulations of hard particles, we define dynamic neighbors as the distinct particles that collide with a given reference one during a specific time interval. This definition allows us to determine the distribution of the number of dynamic neighbors, its average, and its standard deviation. We will show that regardless of the time window used to identify dynamic neighbors, th… ▽ More For molecular dynamics simulations of hard particles, we define dynamic neighbors as the distinct particles that collide with a given reference one during a specific time interval. This definition allows us to determine the distribution of the number of dynamic neighbors, its average, and its standard deviation. We will show that regardless of the time window used to identify dynamic neighbors, their distribution is correlated with diffusion coefficients, structure, and configurational entropy. Thus, it is likely that the distribution of the number of dynamic neighbors may be employed as another tool to gain insights into the dynamic behavior of hard systems. We tested this approach on 2D and 3D systems consisting of monodisperse and binary mixtures of hard disks and spheres. Results show that implementing dynamic neighbors to define order parameters can sharpen the signals where transitions take place. △ Less

Submitted 24 September, 2023; v1 submitted 3 July, 2022; originally announced July 2022.

Comments: 11 pages, 9 figures

Competing magnetic phases in the frustrated spin-1/2 chain compound $β$-TeVO$_4$ probed by NMR

Authors: M. Pregelj, A. Zorko, D. Arčon, M. Klanjšek, N. Janša, P. Jeglič, O. Zaharko, S. Krämer, M. Horvatić, A. Prokofiev

Abstract: In frustrated spin-1/2 chains the competition between the nearest- and next-nearest-neighbor exchange interactions leads to a rich phase diagram that becomes even richer in the presence of perturbations in their material realizations. These effects are still largely unexplored, so that new insight into static and dynamic magnetism, in particular by sensitive local probes, is highly desired. Here w… ▽ More In frustrated spin-1/2 chains the competition between the nearest- and next-nearest-neighbor exchange interactions leads to a rich phase diagram that becomes even richer in the presence of perturbations in their material realizations. These effects are still largely unexplored, so that new insight into static and dynamic magnetism, in particular by sensitive local probes, is highly desired. Here we present a comprehensive $^{17}$O nuclear magnetic resonance study of $β$-TeVO$_4$, where the anisotropy of the main exchange interactions and additional weak interchain exchange interactions complement the theoretical phase diagram. Our results confirm the dynamical nature of the intriguing spin-stripe phase that has been reported in previous studies. In addition, we find that the magnetic order in the high-field phase, which develops just below the magnetization saturation, is consistent with an unusual type of spin-density-wave (SDW) order with different alignments of the magnetic moments on the neighboring chains. This is reminiscent of the ordering in the SDW phase, realized in the absence of the magnetic field, and is thus most likely stabilized by magnetic anisotropy. △ Less

Submitted 24 May, 2022; originally announced May 2022.

Comments: 6 pages, 5 figures

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

The interactions of SARS-CoV-2 with co-circulating pathogens: Epidemiological implications and current knowledge gaps

Authors: Anabelle Wong, Laura Barrero, Elizabeth Goult, Michael Briga, Sarah C. Kramer, Aleksandra Kovacevic, Lulla Opatowski, Matthieu Domenech de Cellès

Abstract: Despite the availability of effective vaccines, the persistence of SARS-CoV-2 suggests that co-circulation with other pathogens and resulting multi-epidemics may become increasingly frequent. To better forecast and control the risk of such multi-epidemics, it is essential to elucidate the potential interactions of SARS-CoV-2 with other pathogens; these interactions, however, remain poorly defined.… ▽ More Despite the availability of effective vaccines, the persistence of SARS-CoV-2 suggests that co-circulation with other pathogens and resulting multi-epidemics may become increasingly frequent. To better forecast and control the risk of such multi-epidemics, it is essential to elucidate the potential interactions of SARS-CoV-2 with other pathogens; these interactions, however, remain poorly defined. Here, we aimed to review the current body of evidence about SARS-CoV-2 interactions. To study pathogen interactions in a systematic way, we first developed a general framework to capture their major components: sign, strength, symmetry, duration, and mechanism. We then reviewed the experimental evidence from animal models about SARS-CoV-2 interactions. Of the 14 studies identified, 11 focused on the outcomes of co-infection with non-attenuated influenza A viruses and generally demonstrated that co-infection increased disease severity compared with either mono-infection. By contrast, the effect of co-infection on the viral load of either virus was variable and inconsistent across studies. Next, we reviewed the epidemiological evidence about SARS-CoV-2 interactions in human populations. Although numerous studies were identified, only few were specifically designed to infer interaction and many were prone to multiple biases, including confounding. Nevertheless, their results suggested that influenza and pneumococcal conjugate vaccinations were associated with reduced riskof SARS-CoV-2 infection. Finally, we formulated simple transmission models of SARS-CoV-2 co-circulation with a viral or a bacterial pathogen, showing how they can naturally incorporate the proposed framework. More generally, we argue that such models, when designed with an integrative and multidisciplinary perspective, will be invaluable tools to resolve the substantial uncertainties that remain about SARS-CoV-2 interactions. △ Less

Submitted 16 November, 2022; v1 submitted 16 May, 2022; originally announced May 2022.

Comments: Main text 35 pages including 5 figures. Appendices 1-5 (Table S1-S4, Model details). Reason for changes: Revision 1, updated literature search (as of 22 Aug 2022)

Deep-Learned Generators of Porosity Distributions Produced During Metal Additive Manufacturing

Authors: Francis Ogoke, Kyle Johnson, Michael Glinsky, Chris Laursen, Sharlotte Kramer, Amir Barati Farimani

Abstract: Laser Powder Bed Fusion has become a widely adopted method for metal Additive Manufacturing (AM) due to its ability to mass produce complex parts with increased local control. However, AM produced parts can be subject to undesirable porosity, negatively influencing the properties of printed components. Thus, controlling porosity is integral for creating effective parts. A precise understanding of… ▽ More Laser Powder Bed Fusion has become a widely adopted method for metal Additive Manufacturing (AM) due to its ability to mass produce complex parts with increased local control. However, AM produced parts can be subject to undesirable porosity, negatively influencing the properties of printed components. Thus, controlling porosity is integral for creating effective parts. A precise understanding of the porosity distribution is crucial for accurately simulating potential fatigue and failure zones. Previous research on generating synthetic porous microstructures have succeeded in generating parts with high density, isotropic porosity distributions but are often inapplicable to cases with sparser, boundary-dependent pore distributions. Our work bridges this gap by providing a method that considers these constraints by deconstructing the generation problem into its constitutive parts. A framework is introduced that combines Generative Adversarial Networks with Mallat Scattering Transform-based autocorrelation methods to construct novel realizations of the individual pore geometries and surface roughness, then stochastically reconstruct them to form realizations of a porous printed part. The generated parts are compared to the existing experimental porosity distributions based on statistical and dimensional metrics, such as nearest neighbor distances, pore volumes, pore anisotropies and scattering transform based auto-correlations. △ Less

Submitted 11 May, 2022; originally announced May 2022.

Efficient optimization of a regional water elevation model with an automatically generated adjoint

Authors: Tuomas Kärnä, Joseph G. Wallwork, Stephan C. Kramer

Abstract: Calibration of unknown model parameters is a common task in many ocean model applications. We present an adjoint-based optimization of an unstructured mesh shallow water model for the Baltic Sea. Spatially varying bottom friction parameter is tuned to minimize the misfit with respect to tide gauge sea surface height (SSH) observations. A key benefit of adjoint-based optimization is that computatio… ▽ More Calibration of unknown model parameters is a common task in many ocean model applications. We present an adjoint-based optimization of an unstructured mesh shallow water model for the Baltic Sea. Spatially varying bottom friction parameter is tuned to minimize the misfit with respect to tide gauge sea surface height (SSH) observations. A key benefit of adjoint-based optimization is that computational cost does not depend on the number of unknown variables. Adjoint models are, however, typically very laborious to implement. In this work, we leverage a domain specific language framework in which the discrete adjoint model can be obtained automatically. The adjoint model is both exactly compatible with the discrete forward model and computationally efficient. A gradient-based quasi-Newton method is used to minimize the misfit. Optimizing spatially-variable parameters is typically an under-determined problem and can lead to over-fitting. We employ Hessian-based regularization to penalize the spatial curvature of the friction field to overcome this problem. The SSH dynamics in the Baltic Sea are simulated for a 3-month period. Optimization of the bottom friction parameter results in significant improvement of the model performance. The results are especially encouraging in the complex Danish Straits region, highlighting the benefit of unstructured meshes. Domain specific language frameworks enable automated model analysis and provide easy access to adjoint modeling. Our application shows that this capability can be enabled with few efforts, and the optimization procedure is robust and computationally efficient. △ Less

Submitted 6 October, 2023; v1 submitted 3 May, 2022; originally announced May 2022.

MSC Class: 76M10; 37N10; 65M32 ACM Class: I.6; J.2

Calibrating constitutive models with full-field data via physics informed neural networks

Authors: Craig M. Hamel, Kevin N. Long, Sharlotte L. B. Kramer

Abstract: The calibration of solid constitutive models with full-field experimental data is a long-standing challenge, especially in materials which undergo large deformation. In this paper, we propose a physics-informed deep-learning framework for the discovery of constitutive model parameterizations given full-field displacement data and global force-displacement data. Contrary to the majority of recent l… ▽ More The calibration of solid constitutive models with full-field experimental data is a long-standing challenge, especially in materials which undergo large deformation. In this paper, we propose a physics-informed deep-learning framework for the discovery of constitutive model parameterizations given full-field displacement data and global force-displacement data. Contrary to the majority of recent literature in this field, we work with the weak form of the governing equations rather than the strong form to impose physical constraints upon the neural network predictions. The approach presented in this paper is computationally efficient, suitable for irregular geometric domains, and readily ingests displacement data without the need for interpolation onto a computational grid. A selection of canonical hyperelastic materials models suitable for different material classes is considered including the Neo-Hookean, Gent, and Blatz-Ko constitutive models as exemplars for general hyperelastic behavior, polymer behavior with lock-up, and compressible foam behavior respectively. We demonstrate that physics informed machine learning is an enabling technology and may shift the paradigm of how full-field experimental data is utilized to calibrate constitutive models under finite deformations. △ Less

Submitted 30 March, 2022; originally announced March 2022.

Quantifying Feedback from Narrow Line Region Outflows in Nearby Active Galaxies. IV. The Effects of Different Density Estimates on the Ionized Gas Masses and Outflow Rates

Authors: Mitchell Revalski, D. Michael Crenshaw, Marc Rafelski, Steven B. Kraemer, Garrett E. Polack, Anna Trindade Falcão, Travis C. Fischer, Beena Meena, Francisco Martinez, Henrique R. Schmitt, Nicholas R. Collins, Julia Falcone

Abstract: Active galactic nuclei (AGN) can launch outflows of ionized gas that may influence galaxy evolution, and quantifying their full impact requires spatially resolved measurements of the gas masses, velocities, and radial extents. We previously reported these quantities for the ionized narrow-line region (NLR) outflows in six low-redshift AGN, where the gas velocities and extents were determined from… ▽ More Active galactic nuclei (AGN) can launch outflows of ionized gas that may influence galaxy evolution, and quantifying their full impact requires spatially resolved measurements of the gas masses, velocities, and radial extents. We previously reported these quantities for the ionized narrow-line region (NLR) outflows in six low-redshift AGN, where the gas velocities and extents were determined from Hubble Space Telescope long-slit spectroscopy. However, calculating the gas masses required multi-component photoionization models to account for radial variations in the gas densities, which span $\sim$6 orders of magnitude. In order to simplify this method for larger samples with less spectral coverage, we compare these gas masses with those calculated from techniques in the literature. First, we use a recombination equation with three different estimates for the radial density profiles. These include constant densities, those derived from [S II], and power-law profiles based on constant values of the ionization parameter ($U$). Second, we use single-component photoionization models with power-law density profiles based on constant $U$, and allow $U$ to vary with radius based on the [O III]/H$β$ ratios. We find that assuming a constant density of $n_\mathrm{H} =$ 10$^2$ cm$^{-3}$ overestimates the gas masses for all six outflows, particularly at small radii where the outflow rates peak. The use of [S II] marginally matches the total gas masses, but also overestimates at small radii. Overall, single-component photoionization models where $U$ varies with radius are able to best match the gas mass and outflow rate profiles when there are insufficient emission lines to construct detailed models. △ Less

Submitted 14 June, 2022; v1 submitted 14 March, 2022; originally announced March 2022.

Comments: Accepted for Publication in ApJ on March 11, 2022. The paper has 20 pages and 8 figures, with results tabulated in the Appendix. Version two includes minor corrections to match the journal publication

Fair Interpretable Representation Learning with Correction Vectors

Authors: Mattia Cerrato, Alesia Vallenas Coronel, Marius Köppel, Alexander Segner, Roberto Esposito, Stefan Kramer

Abstract: Neural network architectures have been extensively employed in the fair representation learning setting, where the objective is to learn a new representation for a given vector which is independent of sensitive information. Various representation debiasing techniques have been proposed in the literature. However, as neural networks are inherently opaque, these methods are hard to comprehend, which… ▽ More Neural network architectures have been extensively employed in the fair representation learning setting, where the objective is to learn a new representation for a given vector which is independent of sensitive information. Various representation debiasing techniques have been proposed in the literature. However, as neural networks are inherently opaque, these methods are hard to comprehend, which limits their usefulness. We propose a new framework for fair representation learning that is centered around the learning of "correction vectors", which have the same dimensionality as the given data vectors. Correction vectors may be computed either explicitly via architectural constraints or implicitly by training an invertible model based on Normalizing Flows. We show experimentally that several fair representation learning models constrained in such a way do not exhibit losses in ranking or classification performance. Furthermore, we demonstrate that state-of-the-art results can be achieved by the invertible model. Finally, we discuss the law standing of our methodology in light of recent legislation in the European Union. △ Less

Submitted 7 February, 2022; originally announced February 2022.

Fair Interpretable Learning via Correction Vectors

Authors: Mattia Cerrato, Marius Köppel, Alexander Segner, Stefan Kramer

Abstract: Neural network architectures have been extensively employed in the fair representation learning setting, where the objective is to learn a new representation for a given vector which is independent of sensitive information. Various "representation debiasing" techniques have been proposed in the literature. However, as neural networks are inherently opaque, these methods are hard to comprehend, whi… ▽ More Neural network architectures have been extensively employed in the fair representation learning setting, where the objective is to learn a new representation for a given vector which is independent of sensitive information. Various "representation debiasing" techniques have been proposed in the literature. However, as neural networks are inherently opaque, these methods are hard to comprehend, which limits their usefulness. We propose a new framework for fair representation learning which is centered around the learning of "correction vectors", which have the same dimensionality as the given data vectors. The corrections are then simply summed up to the original features, and can therefore be analyzed as an explicit penalty or bonus to each feature. We show experimentally that a fair representation learning problem constrained in such a way does not impact performance. △ Less

Submitted 17 January, 2022; originally announced January 2022.

Comments: ICLR-21 Workshop on Responsible AI

Fair Group-Shared Representations with Normalizing Flows

Authors: Mattia Cerrato, Marius Köppel, Alexander Segner, Stefan Kramer

Abstract: The issue of fairness in machine learning stems from the fact that historical data often displays biases against specific groups of people which have been underprivileged in the recent past, or still are. In this context, one of the possible approaches is to employ fair representation learning algorithms which are able to remove biases from data, making groups statistically indistinguishable. In t… ▽ More The issue of fairness in machine learning stems from the fact that historical data often displays biases against specific groups of people which have been underprivileged in the recent past, or still are. In this context, one of the possible approaches is to employ fair representation learning algorithms which are able to remove biases from data, making groups statistically indistinguishable. In this paper, we instead develop a fair representation learning algorithm which is able to map individuals belonging to different groups in a single group. This is made possible by training a pair of Normalizing Flow models and constraining them to not remove information about the ground truth by training a ranking or classification model on top of them. The overall, ``chained'' model is invertible and has a tractable Jacobian, which allows to relate together the probability densities for different groups and ``translate'' individuals from one group to another. We show experimentally that our methodology is competitive with other fair representation learning algorithms. Furthermore, our algorithm achieves stronger invariance w.r.t. the sensitive attribute. △ Less

Submitted 17 January, 2022; originally announced January 2022.

Comments: ICLR-21 Workshop on Responsible AI

Efficient algorithm for simulating particles in real quasiperiodic environments

Authors: Alan Rodrigo Mendoza Sosa, Atahualpa S. Kraemer

Abstract: We introduce an algorithm based on Generalized Dual Method (GDM) to efficiently study the dynamics of a particle in quasiperiodic environments without the need to use periodic approximations or to save the information of the vertices that make up the quasiperiodic lattice. We show that the computation time and the memory required to find the tile in which a particle is located as a function of the… ▽ More We introduce an algorithm based on Generalized Dual Method (GDM) to efficiently study the dynamics of a particle in quasiperiodic environments without the need to use periodic approximations or to save the information of the vertices that make up the quasiperiodic lattice. We show that the computation time and the memory required to find the tile in which a particle is located as a function of the distance $R$ to the center of symmetry remains constant in our algorithm, while using the GDM directly both quantities go like $R^2$.This allows us to perform realistic simulations with low consumption of computational resources. The algorithm can be used to study any quasiperiodic lattice that can be produced by the cut-and-project method. Using this algorithm, we have calculated the free path length distribution in quasiperiodic Lorentz gases reproducing previous results and for systems with high symmetries at the Boltzmann-Grad limit. We have found for the Boltzmann-Grad limit, that the distribution of free paths depends on the rank $r$ of the quasiperiodic system and not on its symmetry. The distribution as a function of the free path length $l$ appears to be a combination of exponential decay and a power-law behavior. The latter seems to become important only for probabilities less than $(2^{r-2} r (r+1))^{-1}$, showing an exponential decaying free-path length distribution for $r \rightarrow \infty$, similar to what is observed in disordered systems. △ Less

Submitted 26 April, 2022; v1 submitted 15 November, 2021; originally announced November 2021.

Comments: 16 pages, 9 figures, paper for the algorithm in https://github.com/AlanRodrigoMendozaSosa/Quasiperiodic-Tiles

Tracking X-ray Outflows with Optical/IR Footprint Lines

Authors: Anna Trindade Falcao, S. B. Kraemer, D. M. Crenshaw, M. Melendez, M. Revalski, T. C. Fischer, H. R. Schmitt, T. J. Turner

Abstract: We use Cloudy photoionisation models to predict the flux profiles for optical/IR emission lines that trace the footprint of X-ray gas, such as [Fe X] 6375A and [Si X] 1.43$μ$m. These are a subset of coronal lines, from ions with ionisation potential $\geq$ that of O VII, i.e., 138eV. The footprint lines are formed in gas over the same range in ionisation state as the H and He-like of O and Ne ions… ▽ More We use Cloudy photoionisation models to predict the flux profiles for optical/IR emission lines that trace the footprint of X-ray gas, such as [Fe X] 6375A and [Si X] 1.43$μ$m. These are a subset of coronal lines, from ions with ionisation potential $\geq$ that of O VII, i.e., 138eV. The footprint lines are formed in gas over the same range in ionisation state as the H and He-like of O and Ne ions, which are also the source of X-ray emission lines. The footprint lines can be detected with optical and IR telescopes, such as the Hubble Space Telescope/STIS and James Webb Space Telescope/NIRSpec, and can potentially be used to measure the kinematics of the extended X-ray emission gas. As a test case, we use the footprints to quantify the properties of the X-ray outflow in the Seyfert 1 galaxy NGC 4151. To confirm the accuracy of our method, we compare our model predictions to the measured flux from archival STIS spectra and previous ground-based studies, and the results are in good agreement. We also use our X-ray footprint method to predict the mass profile for the X-ray emission-line gas in NGC 4151 and derive a total spatially-integrated X-ray mass of $7.8(\pm 2.1) \times 10^{5}~M_{\odot}$, in comparison to $5.4(\pm 1.1) \times 10^{5}~M_{\odot}$ measured from a Chandra X-ray analysis. Our results indicate that high-ionisation footprint emission lines in the optical and near-infrared can be used to accurately trace the kinematics and physical conditions of AGN ionised, X-ray emission-line gas. △ Less

Submitted 18 January, 2022; v1 submitted 21 October, 2021; originally announced October 2021.

Asymptotic Log-Det Sum-of-Ranks Minimization via Tensor (Alternating) Iteratively Reweighted Least Squares

Authors: Sebastian Krämer

Abstract: Affine sum-of-ranks minimization (ASRM) generalizes the affine rank minimization (ARM) problem from matrices to tensors. Here, the interest lies in the ranks of a family $\mathcal{K}$ of different matricizations. Transferring our priorly discussed results on asymptotic log-det rank minimization, we show that iteratively reweighted least squares with weight strength $p = 0$ remains a, theoretically… ▽ More Affine sum-of-ranks minimization (ASRM) generalizes the affine rank minimization (ARM) problem from matrices to tensors. Here, the interest lies in the ranks of a family $\mathcal{K}$ of different matricizations. Transferring our priorly discussed results on asymptotic log-det rank minimization, we show that iteratively reweighted least squares with weight strength $p = 0$ remains a, theoretically and practically, particularly viable method denoted as $\mathrm{IRLS}$-$0\mathcal{K}$. As in the matrix case, we prove global convergence of asymptotic minimizers of the log-det sum-of-ranks function to desired solutions. Further, we show local convergence of $\mathrm{IRLS}$-$0\mathcal{K}$ in dependence of the rate of decline of the therein appearing regularization parameter $γ\searrow 0$. For hierarchical families $\mathcal{K}$, we show how an alternating version ($\mathrm{AIRLS}$-$0\mathcal{K}$, related to prior work under the name $\mathrm{SALSA}$) can be evaluated solely through tensor tree network based operations. The method can thereby be applied to high dimensions through the avoidance of exponential computational complexity. Further, the otherwise crucial rank adaption process becomes essentially superfluous even for completion problems. In numerical experiments, we show that the therefor required subspace restrictions and relaxation of the affine constraint cause only a marginal loss of approximation quality. On the other hand, we demonstrate that $\mathrm{IRLS}$-$0\mathcal{K}$ allows to observe the theoretical phase transition also for generic tensor recoverability in practice. Concludingly, we apply $\mathrm{AIRLS}$-$0\mathcal{K}$ to larger scale problems. △ Less

Submitted 29 June, 2021; originally announced June 2021.

Comments: 29 pages + 20 pages supplementary material

MSC Class: 15A03; 15A29; 65J20; 90C31; 90C26

Asymptotic Log-Det Rank Minimization via (Alternating) Iteratively Reweighted Least Squares

Authors: Sebastian Krämer

Abstract: The affine rank minimization (ARM) problem is well known for both its applications and the fact that it is NP-hard. One of the most successful approaches, yet arguably underrepresented, is iteratively reweighted least squares (IRLS), more specifically $\mathrm{IRLS}$-$0$. Despite comprehensive empirical evidence that it overall outperforms nuclear norm minimization and related methods, it is still… ▽ More The affine rank minimization (ARM) problem is well known for both its applications and the fact that it is NP-hard. One of the most successful approaches, yet arguably underrepresented, is iteratively reweighted least squares (IRLS), more specifically $\mathrm{IRLS}$-$0$. Despite comprehensive empirical evidence that it overall outperforms nuclear norm minimization and related methods, it is still not understood to a satisfying degree. In particular, the significance of a slow decrease of the therein appearing regularization parameter denoted $γ$ poses interesting questions. While commonly equated to matrix recovery, we here consider the ARM independently. We investigate the particular structure and global convergence property behind the asymptotic minimization of the log-det objective function on which $\mathrm{IRLS}$-$0$ is based. We expand on local convergence theorems, now with an emphasis on the decline of $γ$, and provide representative examples as well as counterexamples such as a diverging $\mathrm{IRLS}$-$0$ sequence that clarify theoretical limits. We present a data sparse, alternating realization $\mathrm{AIRLS}$-$p$ (related to prior work under the name $\mathrm{SALSA}$) that, along with the rest of this work, serves as basis and introduction to the more general tensor setting. In conclusion, numerical sensitivity experiments are carried out that reconfirm the success of $\mathrm{IRLS}$-$0$ and demonstrate that in surprisingly many cases, a slower decay of $γ$ will yet lead to a solution of the ARM problem, up to the point that the exact theoretical phase transition for generic recoverability can be observed. Likewise, this suggests that non-convexity is less substantial and problematic for the log-det approach than it might initially appear. △ Less

Submitted 28 June, 2021; originally announced June 2021.

Comments: 26 pages + 16 pages supplementary material

MSC Class: 15A03; 15A29; 65J20; 90C31; 90C26

Hubble Space Telescope [O III] Emission-Line Kinematics in Two Nearby QSO2s: A Case for X-ray Feedback

Authors: Anna Trindade Falcão, S. B. Kraemer, T. C. Fischer, D. M. Crenshaw, M. Revalski, H. R. Schmitt, W. P. Maksym, M. Vestergaard, M. Elvis, C. M. Gaskell, F. Hamann, L. C. Ho, J. Hutchings, R. Mushotzky, H. Netzer, T. Storchi-Bergmann, T. J. Turner, M. J. Ward

Abstract: We present a dynamical study of the narrow-line regions in two nearby QSO2s. We construct dynamical models based on detailed photoionization models of the emission-line gas, including the effects of internal dust, to apply to observations of large-scale outflows from these AGNs. We use Mrk 477 and Mrk 34 in order to test our models against recent HST STIS observations of [O III] emission-line kine… ▽ More We present a dynamical study of the narrow-line regions in two nearby QSO2s. We construct dynamical models based on detailed photoionization models of the emission-line gas, including the effects of internal dust, to apply to observations of large-scale outflows from these AGNs. We use Mrk 477 and Mrk 34 in order to test our models against recent HST STIS observations of [O III] emission-line kinematics since these AGNs possess more energetic outflows than found in Seyfert galaxies. We find that the outflows within 500 pc are consistent with radiative acceleration of dusty gas, however, the outflows in Mrk 34 are significantly more extended and may not be directly accelerated by radiation. We characterize the properties of X-ray winds found from the expansion of [O III]-emitting gas close to the black hole. We show that such winds possess the kinetic energy density to disturb [O III] gas at 1.8 kpc, and have sufficient energy to entrain the [O III] clouds at 1.2 kpc. Assuming that the X-ray wind possesses the same radial mass distribution as the [O III] gas, we find that the peak kinetic luminosity for this wind is 2% of Mrk 34's bolometric luminosity, which is in the 0.5% - 5% range required by some models for efficient feedback. Our work shows that, although the kinetic luminosity as measured from [O III]-emitting gas is frequently low, X-ray winds may provide more than one order of magnitude higher kinetic power. △ Less

Submitted 25 May, 2021; originally announced May 2021.

Fast Private Parameter Learning and Inference for Sum-Product Networks

Authors: Ernst Althaus, Mohammad Sadeq Dousti, Stefan Kramer, Nick Johannes Peter Rassau

Abstract: A sum-product network (SPN) is a graphical model that allows several types of inferences to be drawn efficiently. There are two types of learning for SPNs: Learning the architecture of the model, and learning the parameters. In this paper, we tackle the second problem: We show how to learn the weights for the sum nodes, assuming the architecture is fixed, and the data is horizontally partitioned b… ▽ More A sum-product network (SPN) is a graphical model that allows several types of inferences to be drawn efficiently. There are two types of learning for SPNs: Learning the architecture of the model, and learning the parameters. In this paper, we tackle the second problem: We show how to learn the weights for the sum nodes, assuming the architecture is fixed, and the data is horizontally partitioned between multiple parties. The computations will preserve the privacy of each participant. Furthermore, we will use secret sharing instead of (homomorphic) encryption, which allows fast computations and requires little computational resources. To this end, we use a novel integer division to compute approximate real divisions. We also show how simple and private inferences can be performed using the learned SPN. △ Less

Submitted 14 October, 2021; v1 submitted 15 April, 2021; originally announced April 2021.

Gauging the effect of Supermassive Black Holes feedback on Quasar host galaxies

Authors: B. Dall'Agnol de Oliveira, T. Storchi-Bergmann, S. B. Kraemer, M. Villar Martín, A. Schnorr-Müller, H. R. Schmitt, D. Ruschel-Dutra, D. M. Crenshaw, T. C. Fischer

Abstract: In order to gauge the role that active galactic nuclei (AGN) play in the evolution of galaxies via the effect of kinetic feedback in nearby QSO$\,$2's ($z\sim0.3$), we observed eight such objects with bolometric luminosities $L_{bol} \sim 10^{46}\rm{erg\,s^{-1}}$ using Gemini GMOS-IFU's. The emission lines were fitted with at least two Gaussian curves, the broadest of which we attributed to gas ki… ▽ More In order to gauge the role that active galactic nuclei (AGN) play in the evolution of galaxies via the effect of kinetic feedback in nearby QSO$\,$2's ($z\sim0.3$), we observed eight such objects with bolometric luminosities $L_{bol} \sim 10^{46}\rm{erg\,s^{-1}}$ using Gemini GMOS-IFU's. The emission lines were fitted with at least two Gaussian curves, the broadest of which we attributed to gas kinetically disturbed by an outflow. We found that the maximum extent of the outflow ranges from $\sim$1 to 8 kpc, being ${\sim}\,0.5\,{\pm}\,0.3$ times the extent of the [O$\,$III] ionized gas region. Our `${\tt default}$' assumptions for the gas density (obtained from the {[S$\,$II] doublet) and outflow velocities resulted in peak mass outflow rates of $\dot{M}_{out}^{\tt defa}{\sim}\,3\,{-}\,30\,\rm{M_{\odot}}\,yr^{-1}$ and outflow power of $\dot{E}_{out}^{\tt defa}\sim\,10^{41}\,{-}\,10^{43}\,\mathrm{erg\,s^{-1}}$. The corresponding kinetic coupling efficiencies are $\varepsilon_f^{\tt defa}=\dot{E}_{out}^{\tt defa}/L_{bol}\,\sim7\times10^{-4}\,{-}\,0.5\,\%$, with the average efficiency being only $0.06\,\%$ ($0.01\,\%$ median), implying little feedback powers from ionized gas outflows in the host galaxies. We investigated the effects of varying assumptions and calculations on $\dot{M}_{out}$ and $\dot{E}_{out}$ regarding the ionized gas densities, velocities, masses and inclinations of the outflow relative to the plane of the sky, resulting in average uncertainties of one dex. In particular, we found that better indicators of the [O$\,$III] emitting gas density than the default [S$\,$II] line ratio, such as the [Ar$\,$IV]$λλ$4711,40 line ratio, result in almost an order of magnitude decrease in the $\varepsilon_f$. △ Less

Submitted 13 April, 2021; originally announced April 2021.

Comments: Accepted: 13 april 2021. Online supplementary material added at the end of the pdf

Tunneling spectroscopy of few-monolayer NbSe$_2$ in high magnetic field: Ising protection and triplet superconductivity

Authors: M. Kuzmanović, T. Dvir, D. LeBoeuf, S. Ilić, M. Haim, D. Möckli, S. Kraemer, M. Khodas, M. Houzet, J. S. Meyer, M. Aprili, H. Steinberg, C. H. L. Quay

Abstract: In conventional Bardeen-Cooper-Scrieffer (BCS) superconductors, Cooper pairs of electrons of opposite spin (i.e. singlet structure) form the ground state. Equal spin triplet pairs (ESTPs), as in superfluid $^3$He, are of great interest for superconducting spintronics and topological superconductivity, yet remain elusive. Recently, odd-parity ESTPs were predicted to arise in (few-)monolayer superco… ▽ More In conventional Bardeen-Cooper-Scrieffer (BCS) superconductors, Cooper pairs of electrons of opposite spin (i.e. singlet structure) form the ground state. Equal spin triplet pairs (ESTPs), as in superfluid $^3$He, are of great interest for superconducting spintronics and topological superconductivity, yet remain elusive. Recently, odd-parity ESTPs were predicted to arise in (few-)monolayer superconducting NbSe$_2$, from the non-colinearity between the out-of-plane Ising spin-orbit field (due to the lack of inversion symmetry in monolayer NbSe$_2$) and an applied in-plane magnetic field. These ESTPs couple to the singlet order parameter at finite field. Using van der Waals tunnel junctions, we perform spectroscopy of superconducting NbSe$_2$ flakes, of 2--25 monolayer thickness, measuring the quasiparticle density of states (DOS) as a function of applied in-plane magnetic field up to 33T. In flakes $\lesssim$ 15 monolayers thick the DOS has a single superconducting gap. In these thin samples, the magnetic field acts primarily on the spin (vs orbital) degree of freedom of the electrons, and superconductivity is further protected by the Ising field. The superconducting energy gap, extracted from our tunnelling spectra, decreases as a function of the applied magnetic field. However, in bilayer NbSe$_2$, close to the critical field (up to 30T, much larger than the Pauli limit), superconductivity appears to be more robust than expected from Ising protection alone. Our data can be explained by the above-mentioned ESTPs. △ Less

Submitted 2 August, 2022; v1 submitted 1 April, 2021; originally announced April 2021.

Locally commensurate charge-density wave with three-unit-cell periodicity in YBCO

Authors: I. Vinograd, R. Zhou, M. Hirata, T. Wu, H. Mayaffre, S. Krämer, R. Liang, W. N. Hardy, D. A. Bonn, M. -H. Julien

Abstract: In order to identify the mechanism responsible for the formation of charge-density waves (CDW) in cuprate superconductors, it is important to understand which aspects of the CDW's microscopic structure are generic and which are material-dependent. Here, we show that, at the local scale probed by NMR, long-range CDW order in YBa2Cu3Oy is unidirectional with a commensurate period of three unit cells… ▽ More In order to identify the mechanism responsible for the formation of charge-density waves (CDW) in cuprate superconductors, it is important to understand which aspects of the CDW's microscopic structure are generic and which are material-dependent. Here, we show that, at the local scale probed by NMR, long-range CDW order in YBa2Cu3Oy is unidirectional with a commensurate period of three unit cells (lambda = 3b), implying that the incommensurability found in X-ray scattering is ensured by phase slips (discommensurations). Furthermore, NMR spectra reveal a predominant oxygen character of the CDW with an out-of-phase relationship between certain lattice sites but no specific signature of a secondary CDW with lambda = 6b associated with a putative pair-density wave. These results shed light on universal aspects of the cuprate CDW. In particular, its spatial profile appears to generically result from the interplay between an incommensurate tendency at long length scales, possibly related to properties of the Fermi surface, and local commensuration effects, due to electron-electron interactions or lock-in to the lattice. △ Less

Submitted 24 March, 2021; originally announced March 2021.

Comments: Original submission (revised version available upon request)

Journal ref: Nat Commun 12, 3274 (2021)

Radiative Driving of the AGN Outflows in the Narrow-Line Seyfert 1 Galaxy NGC 4051

Authors: Beena Meena, D. Michael Crenshaw, Henrique R. Schmitt, Mitchell Revalski, Travis C. Fischer, Garrett E. Polack, Steven B. Kraemer, Dzhuliya Dashtamirova

Abstract: We explore the properties of ionized gas in the nuclear and circumnuclear environment of the narrow-line Seyfert 1 galaxy NGC 4051 using spectroscopic and imaging observations from the Hubble Space Telescope (HST) and Apache Point Observatory (APO)'s ARC 3.5m Telescope. We identify an unresolved moderate-density intermediate width component and a high-density broad component in the optical emissio… ▽ More We explore the properties of ionized gas in the nuclear and circumnuclear environment of the narrow-line Seyfert 1 galaxy NGC 4051 using spectroscopic and imaging observations from the Hubble Space Telescope (HST) and Apache Point Observatory (APO)'s ARC 3.5m Telescope. We identify an unresolved moderate-density intermediate width component and a high-density broad component in the optical emission lines from the active nucleus, as well as spatially-resolved emission extending up to $\sim$1 kpc in the AGN ionized narrow-line region (NLR) and $\sim$8 kpc in the stellar ionized host galaxy. The HST narrow-band image reveals a distinct conical structure in [O III] emission towards the NE, and the ionized gas kinematics shows up to two blueshifted velocity components, indicating outflows along the edges of a cone. We introduce an improved model of biconical outflow, with our line of sight passing through the wall of the cone, which suggests that the large number of outflowing UV absorbers seen in NGC 4051 are NLR clouds in absorption. Using the de-projection factors from the biconical geometry, we measure true outflow velocities up to 680 km s$^{-1}$ at a distance of $\sim$350 pc, however, we do not find any rotational signature inside a projected distance $\leq$ 10'' ($\sim$800 pc) from the nucleus. We compare the gas kinematics with analytical models based on a radiation-gravity formalism, which show that most of the observed NLR outflows are launched within $\sim$0.5 pc of the nucleus and can travel up to $\sim$1 kpc from this low-luminosity AGN. △ Less

Submitted 18 May, 2021; v1 submitted 22 March, 2021; originally announced March 2021.

Comments: 29 Pages, 16 Figures, 4 Tables. Accepted for publication in ApJ

Revealing three-dimensional quantum criticality by Sr-substitution in Han Purple

Authors: Stephan Allenspach, Pascal Puphal, Joosep Link, Ivo Heinmaa, Ekaterina Pomjakushina, Cornelius Krellner, Jakob Lass, Gregory S. Tucker, Christof Niedermayer, Shusaku Imajo, Yoshimitsu Kohama, Koichi Kindo, Steffen Krämer, Mladen Horvatić, Marcelo Jaime, Alexander Madsen, Antonietta Mira, Nicolas Laflorencie, Frédéric Mila, Bruce Normand, Christian Rüegg, Raivo Stern, Franziska Weickert

Abstract: Classical and quantum phase transitions (QPTs), with their accompanying concepts of criticality and universality, are a cornerstone of statistical thermodynamics. An exemplary controlled QPT is the field-induced magnetic ordering of a gapped quantum magnet. Although numerous "quasi-one-dimensional" coupled spin-chain and -ladder materials are known whose ordering transition is three-dimensional (3… ▽ More Classical and quantum phase transitions (QPTs), with their accompanying concepts of criticality and universality, are a cornerstone of statistical thermodynamics. An exemplary controlled QPT is the field-induced magnetic ordering of a gapped quantum magnet. Although numerous "quasi-one-dimensional" coupled spin-chain and -ladder materials are known whose ordering transition is three-dimensional (3D), quasi-2D systems are special for several physical reasons. Motivated by the ancient pigment Han Purple (BaCuSi$_{2}$O$_{6}$), a quasi-2D material displaying anomalous critical properties, we present a complete analysis of Ba$_{0.9}$Sr$_{0.1}$CuSi$_{2}$O$_{6}$. We measure the zero-field magnetic excitations by neutron spectroscopy and deduce the magnetic Hamiltonian. We probe the field-induced transition by combining magnetization, specific-heat, torque and magnetocalorimetric measurements with low-temperature nuclear magnetic resonance studies near the QPT. By a Bayesian statistical analysis and large-scale Quantum Monte Carlo simulations, we demonstrate unambiguously that observable 3D quantum critical scaling is restored by the structural simplification arising from light Sr-substitution in Han Purple. △ Less

Submitted 9 June, 2021; v1 submitted 11 March, 2021; originally announced March 2021.

Comments: 19 pages, 5 figures

Journal ref: Phys. Rev. Research 3, 023177 (2021)