Showing 1–50 of 72 results for author: Rasmussen, M

PARAFAC2-based Coupled Matrix and Tensor Factorizations with Constraints

Authors: Carla Schenker, Xiulin Wang, David Horner, Morten A. Rasmussen, Evrim Acar

Abstract: Data fusion models based on Coupled Matrix and Tensor Factorizations (CMTF) have been effective tools for joint analysis of data from multiple sources. While the vast majority of CMTF models are based on the strictly multilinear CANDECOMP/PARAFAC (CP) tensor model, recently also the more flexible PARAFAC2 model has been integrated into CMTF models. PARAFAC2 tensor models can handle irregular/ragge… ▽ More Data fusion models based on Coupled Matrix and Tensor Factorizations (CMTF) have been effective tools for joint analysis of data from multiple sources. While the vast majority of CMTF models are based on the strictly multilinear CANDECOMP/PARAFAC (CP) tensor model, recently also the more flexible PARAFAC2 model has been integrated into CMTF models. PARAFAC2 tensor models can handle irregular/ragged tensors and have shown to be especially useful for modelling dynamic data with unaligned or irregular time profiles. However, existing PARAFAC2-based CMTF models have limitations in terms of possible regularizations on the factors and/or types of coupling between datasets. To address these limitations, in this paper we introduce a flexible algorithmic framework that fits PARAFAC2-based CMTF models using Alternating Optimization (AO) and the Alternating Direction Method of Multipliers (ADMM). The proposed framework allows to impose various constraints on all modes and linear couplings to other matrix-, CP- or PARAFAC2-models. Experiments on various simulated and a real dataset demonstrate the utility and versatility of the proposed framework as well as its benefits in terms of accuracy and efficiency in comparison with state-of-the-art methods. △ Less

Submitted 18 June, 2024; originally announced June 2024.

Comments: 15 pages, 15 figures,1 table

Segment anything model for head and neck tumor segmentation with CT, PET and MRI multi-modality images

Authors: Jintao Ren, Mathis Rasmussen, Jasper Nijkamp, Jesper Grau Eriksen, Stine Korreman

Abstract: Deep learning presents novel opportunities for the auto-segmentation of gross tumor volume (GTV) in head and neck cancer (HNC), yet fully automatic methods usually necessitate significant manual refinement. This study investigates the Segment Anything Model (SAM), recognized for requiring minimal human prompting and its zero-shot generalization ability across natural images. We specifically examin… ▽ More Deep learning presents novel opportunities for the auto-segmentation of gross tumor volume (GTV) in head and neck cancer (HNC), yet fully automatic methods usually necessitate significant manual refinement. This study investigates the Segment Anything Model (SAM), recognized for requiring minimal human prompting and its zero-shot generalization ability across natural images. We specifically examine MedSAM, a version of SAM fine-tuned with large-scale public medical images. Despite its progress, the integration of multi-modality images (CT, PET, MRI) for effective GTV delineation remains a challenge. Focusing on SAM's application in HNC GTV segmentation, we assess its performance in both zero-shot and fine-tuned scenarios using single (CT-only) and fused multi-modality images. Our study demonstrates that fine-tuning SAM significantly enhances its segmentation accuracy, building upon the already effective zero-shot results achieved with bounding box prompts. These findings open a promising avenue for semi-automatic HNC GTV segmentation. △ Less

Submitted 27 February, 2024; originally announced February 2024.

Attractors of Linear Maps with Bounded Noise

Authors: Jeroen S. W. Lamb, Martin Rasmussen, Wei Hao Tey

Abstract: We consider invertible linear maps with additive spherical bounded noise. We show that minimal attractors of such random dynamical systems are unique, strictly convex and have a continuously differentiable boundary. Moreover, we present an auxiliary finite-dimensional deterministic boundary map for which the unit normal bundle of this boundary is globally attracting. We consider invertible linear maps with additive spherical bounded noise. We show that minimal attractors of such random dynamical systems are unique, strictly convex and have a continuously differentiable boundary. Moreover, we present an auxiliary finite-dimensional deterministic boundary map for which the unit normal bundle of this boundary is globally attracting. △ Less

Submitted 5 October, 2023; originally announced October 2023.

Comments: 8 pages, 1 figure

MSC Class: 37C05; 37H10; 37B35

Sensitivity of Simulations of Double Detonation Type Ia Supernova to Integration Methodology

Authors: Michael Zingale, Zhi Chen, Melissa Rasmussen, Abigail Polin, Max Katz, Alexander Smith Clark, Eric T. Johnson

Abstract: We study the coupling of hydrodynamics and reactions in simulations of the double detonation model for Type Ia supernovae. When assessing the convergence of simulations, the focus is usually on spatial resolution; however, the method of coupling the physics together as well as the tolerances used in integrating a reaction network also play an important role. In this paper, we explore how the choic… ▽ More We study the coupling of hydrodynamics and reactions in simulations of the double detonation model for Type Ia supernovae. When assessing the convergence of simulations, the focus is usually on spatial resolution; however, the method of coupling the physics together as well as the tolerances used in integrating a reaction network also play an important role. In this paper, we explore how the choices made in both coupling and integrating the reaction portion of a simulation (operator / Strang splitting vs.\ the simplified spectral deferred corrections method we introduced previously) influences the accuracy, efficiency, and the nucleosynthesis of simulations of double detonations. We find no need to limit reaction rates or reduce the simulation timestep to the reaction timescale. The entire simulation methodology used here is GPU-accelerated and made freely available as part of the Castro simulation code. △ Less

Submitted 13 March, 2024; v1 submitted 4 September, 2023; originally announced September 2023.

Comments: accepted to Astrophysical Journal this version has a new appendix (on shock burning) compared to previous in response to the referee

MAPSeg: Unified Unsupervised Domain Adaptation for Heterogeneous Medical Image Segmentation Based on 3D Masked Autoencoding and Pseudo-Labeling

Authors: Xuzhe Zhang, Yuhao Wu, Elsa Angelini, Ang Li, Jia Guo, Jerod M. Rasmussen, Thomas G. O'Connor, Pathik D. Wadhwa, Andrea Parolin Jackowski, Hai Li, Jonathan Posner, Andrew F. Laine, Yun Wang

Abstract: Robust segmentation is critical for deriving quantitative measures from large-scale, multi-center, and longitudinal medical scans. Manually annotating medical scans, however, is expensive and labor-intensive and may not always be available in every domain. Unsupervised domain adaptation (UDA) is a well-studied technique that alleviates this label-scarcity problem by leveraging available labels fro… ▽ More Robust segmentation is critical for deriving quantitative measures from large-scale, multi-center, and longitudinal medical scans. Manually annotating medical scans, however, is expensive and labor-intensive and may not always be available in every domain. Unsupervised domain adaptation (UDA) is a well-studied technique that alleviates this label-scarcity problem by leveraging available labels from another domain. In this study, we introduce Masked Autoencoding and Pseudo-Labeling Segmentation (MAPSeg), a $\textbf{unified}$ UDA framework with great versatility and superior performance for heterogeneous and volumetric medical image segmentation. To the best of our knowledge, this is the first study that systematically reviews and develops a framework to tackle four different domain shifts in medical image segmentation. More importantly, MAPSeg is the first framework that can be applied to $\textbf{centralized}$, $\textbf{federated}$, and $\textbf{test-time}$ UDA while maintaining comparable performance. We compare MAPSeg with previous state-of-the-art methods on a private infant brain MRI dataset and a public cardiac CT-MRI dataset, and MAPSeg outperforms others by a large margin (10.5 Dice improvement on the private MRI dataset and 5.7 on the public CT-MRI dataset). MAPSeg poses great practical value and can be applied to real-world problems. GitHub: https://github.com/XuzheZ/MAPSeg/. △ Less

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

Comments: CVPR 2024 camera-ready (8 pages, 3 figures) with the supplemental materials (5 pages, 4 figures). Xuzhe Zhang and Yuhao Wu are co-first authors. Andrew F. Laine and Yun Wang are co-senior supervising authors

Persistence of Minimal Invariant Sets for Certain Set-Valued Dynamical Systems: A Boundary Mapping Approach

Authors: K. Kourliouros, J. S. W. Lamb, M. Rasmussen, W. H. Tey, K. G. Timperi, D. Turaev

Abstract: We study the problem of persistence of minimal invariant sets with smooth boundary for a certain class of discrete-time set-valued dynamical systems, naturally arising in the context of random dynamical systems with bounded noise. In particular, we introduce a single-valued map, the so-called boundary map, which has the property that a certain class of invariant submanifolds for this map is in one… ▽ More We study the problem of persistence of minimal invariant sets with smooth boundary for a certain class of discrete-time set-valued dynamical systems, naturally arising in the context of random dynamical systems with bounded noise. In particular, we introduce a single-valued map, the so-called boundary map, which has the property that a certain class of invariant submanifolds for this map is in one-to-one correspondence with invariant sets for the corresponding set-valued map. We show that minimal invariant sets with smooth boundary persist under small perturbations of the set-valued map, provided that the associated boundary map is normally hyperbolic at the unit normal bundle of the boundary. △ Less

Submitted 3 March, 2023; originally announced March 2023.

MSC Class: 37H30; 37C05; 37D05; 37G35; 53D10

A Two-Stage Deep Representation Learning-Based Speech Enhancement Method Using Variational Autoencoder and Adversarial Training

Authors: Yang Xiang, Jesper Lisby Højvang, Morten Højfeldt Rasmussen, Mads Græsbøll Christensen

Abstract: This paper focuses on leveraging deep representation learning (DRL) for speech enhancement (SE). In general, the performance of the deep neural network (DNN) is heavily dependent on the learning of data representation. However, the DRL's importance is often ignored in many DNN-based SE algorithms. To obtain a higher quality enhanced speech, we propose a two-stage DRL-based SE method through advers… ▽ More This paper focuses on leveraging deep representation learning (DRL) for speech enhancement (SE). In general, the performance of the deep neural network (DNN) is heavily dependent on the learning of data representation. However, the DRL's importance is often ignored in many DNN-based SE algorithms. To obtain a higher quality enhanced speech, we propose a two-stage DRL-based SE method through adversarial training. In the first stage, we disentangle different latent variables because disentangled representations can help DNN generate a better enhanced speech. Specifically, we use the $β$-variational autoencoder (VAE) algorithm to obtain the speech and noise posterior estimations and related representations from the observed signal. However, since the posteriors and representations are intractable and we can only apply a conditional assumption to estimate them, it is difficult to ensure that these estimations are always pretty accurate, which may potentially degrade the final accuracy of the signal estimation. To further improve the quality of enhanced speech, in the second stage, we introduce adversarial training to reduce the effect of the inaccurate posterior towards signal reconstruction and improve the signal estimation accuracy, making our algorithm more robust for the potentially inaccurate posterior estimations. As a result, better SE performance can be achieved. The experimental results indicate that the proposed strategy can help similar DNN-based SE algorithms achieve higher short-time objective intelligibility (STOI), perceptual evaluation of speech quality (PESQ), and scale-invariant signal-to-distortion ratio (SI-SDR) scores. Moreover, the proposed algorithm can also outperform recent competitive SE algorithms. △ Less

Submitted 27 September, 2023; v1 submitted 16 November, 2022; originally announced November 2022.

Comments: Accepted by IEEE/ACM Transactions on Audio, Speech and Language Processing

On the quasi-ergodicity of absorbing Markov chains with unbounded transition densities, including random logistic maps with escape

Authors: Matheus M. Castro, Vincent P. H. Goverse, Jeroen S. W. Lamb, Martin Rasmussen

Abstract: In this paper, we consider absorbing Markov chains $X_n$ admitting a quasi-stationary measure $μ$ on $M$ where the transition kernel $\mathcal P$ admits an eigenfunction $0\leq η\in L^1(M,μ)$. We find conditions on the transition densities of $\mathcal P$ with respect to $μ$ which ensure that $η(x) μ(\mathrm d x)$ is a quasi-ergodic measure for $X_n$ and that the Yaglom limit converges to the quas… ▽ More In this paper, we consider absorbing Markov chains $X_n$ admitting a quasi-stationary measure $μ$ on $M$ where the transition kernel $\mathcal P$ admits an eigenfunction $0\leq η\in L^1(M,μ)$. We find conditions on the transition densities of $\mathcal P$ with respect to $μ$ which ensure that $η(x) μ(\mathrm d x)$ is a quasi-ergodic measure for $X_n$ and that the Yaglom limit converges to the quasi-stationary measure $μ$-almost surely. We apply this result to the random logistic map $X_{n+1} = ω_n X_n (1-X_n)$ absorbed at $\mathbb R \setminus [0,1],$ where $ω_n$ is an i.i.d sequence of random variables uniformly distributed in $[a,b],$ for $1\leq a <4$ and $b>4.$ △ Less

Submitted 26 March, 2024; v1 submitted 2 September, 2022; originally announced September 2022.

Comments: 35 pages, 0 figures

MSC Class: 37H05; 47B65; 60J05

An Improved Method for Coupling Hydrodynamics with Astrophysical Reaction Networks

Authors: M. Zingale, M. P. Katz, A. Nonaka, M. Rasmussen

Abstract: Reacting astrophysical flows can be challenging to model because of the difficulty in accurately coupling hydrodynamics and reactions. This can be particularly acute during explosive burning or at high temperatures where nuclear statistical equilibrium is established. We develop a new approach based on the ideas of spectral deferred corrections (SDC) coupling of explicit hydrodynamics and stiff re… ▽ More Reacting astrophysical flows can be challenging to model because of the difficulty in accurately coupling hydrodynamics and reactions. This can be particularly acute during explosive burning or at high temperatures where nuclear statistical equilibrium is established. We develop a new approach based on the ideas of spectral deferred corrections (SDC) coupling of explicit hydrodynamics and stiff reaction sources as an alternative to operator splitting that is simpler than the more comprehensive SDC approach we demonstrated previously. We apply the new method to a double detonation problem with a moderately-sized astrophysical nuclear reaction network and explore the timestep size and reaction network tolerances to show that the simplified-SDC approach provides improved coupling with decreased computational expense compared to traditional Strang operator splitting. This is all done in the framework of the Castro hydrodynamics code, and all algorithm implementations are freely available. △ Less

Submitted 26 July, 2022; v1 submitted 2 June, 2022; originally announced June 2022.

Comments: accepted to ApJ. Castro is available at https://github.com/amrex-astro/Castro -- all code for the results here is in the github repo

A deep representation learning speech enhancement method using $β$-VAE

Authors: Yang Xiang, Jesper Lisby Højvang, Morten Højfeldt Rasmussen, Mads Græsbøll Christensen

Abstract: In previous work, we proposed a variational autoencoder-based (VAE) Bayesian permutation training speech enhancement (SE) method (PVAE) which indicated that the SE performance of the traditional deep neural network-based (DNN) method could be improved by deep representation learning (DRL). Based on our previous work, we in this paper propose to use $β$-VAE to further improve PVAE's ability of repr… ▽ More In previous work, we proposed a variational autoencoder-based (VAE) Bayesian permutation training speech enhancement (SE) method (PVAE) which indicated that the SE performance of the traditional deep neural network-based (DNN) method could be improved by deep representation learning (DRL). Based on our previous work, we in this paper propose to use $β$-VAE to further improve PVAE's ability of representation learning. More specifically, our $β$-VAE can improve PVAE's capacity of disentangling different latent variables from the observed signal without the trade-off problem between disentanglement and signal reconstruction. This trade-off problem widely exists in previous $β$-VAE algorithms. Unlike the previous $β$-VAE algorithms, the proposed $β$-VAE strategy can also be used to optimize the DNN's structure. This means that the proposed method can not only improve PVAE's SE performance but also reduce the number of PVAE training parameters. The experimental results show that the proposed method can acquire better speech and noise latent representation than PVAE. Meanwhile, it also obtains a higher scale-invariant signal-to-distortion ratio, speech quality, and speech intelligibility. △ Less

Submitted 11 May, 2022; originally announced May 2022.

Comments: Submitted to Eurosipco

The Lyapunov spectrum for conditioned random dynamical systems

Authors: Matheus M. Castro, Dennis Chemnitz, Hugo Chu, Maximilian Engel, Jeroen S. W. Lamb, Martin Rasmussen

Abstract: We establish the existence of a full spectrum of Lyapunov exponents for memoryless random dynamical systems with absorption. To this end, we crucially embed the process conditioned to never being absorbed, the Q-process, into the framework of random dynamical systems, allowing us to study multiplicative ergodic properties. We show that the finite-time Lyapunov exponents converge in conditioned pro… ▽ More We establish the existence of a full spectrum of Lyapunov exponents for memoryless random dynamical systems with absorption. To this end, we crucially embed the process conditioned to never being absorbed, the Q-process, into the framework of random dynamical systems, allowing us to study multiplicative ergodic properties. We show that the finite-time Lyapunov exponents converge in conditioned probability and apply our results to iterated function systems and stochastic differential equations. △ Less

Submitted 8 April, 2022; originally announced April 2022.

Comments: 36 pages, 0 figures

MSC Class: 37H05; 37H15; 47D07; 60J05; 60J25

A Bayesian Permutation training deep representation learning method for speech enhancement with variational autoencoder

Authors: Yang Xiang, Jesper Lisby Højvang, Morten Højfeldt Rasmussen, Mads Græsbøll Christensen

Abstract: Recently, variational autoencoder (VAE), a deep representation learning (DRL) model, has been used to perform speech enhancement (SE). However, to the best of our knowledge, current VAE-based SE methods only apply VAE to the model speech signal, while noise is modeled using the traditional non-negative matrix factorization (NMF) model. One of the most important reasons for using NMF is that these… ▽ More Recently, variational autoencoder (VAE), a deep representation learning (DRL) model, has been used to perform speech enhancement (SE). However, to the best of our knowledge, current VAE-based SE methods only apply VAE to the model speech signal, while noise is modeled using the traditional non-negative matrix factorization (NMF) model. One of the most important reasons for using NMF is that these VAE-based methods cannot disentangle the speech and noise latent variables from the observed signal. Based on Bayesian theory, this paper derives a novel variational lower bound for VAE, which ensures that VAE can be trained in supervision, and can disentangle speech and noise latent variables from the observed signal. This means that the proposed method can apply the VAE to model both speech and noise signals, which is totally different from the previous VAE-based SE works. More specifically, the proposed DRL method can learn to impose speech and noise signal priors to different sets of latent variables for SE. The experimental results show that the proposed method can not only disentangle speech and noise latent variables from the observed signal but also obtain a higher scale-invariant signal-to-distortion ratio and speech quality score than the similar deep neural network-based (DNN) SE method. △ Less

Submitted 24 January, 2022; originally announced January 2022.

Comments: Accepted by ICASSP 2022

Existence and uniqueness of quasi-stationary and quasi-ergodic measures for absorbing Markov chains: a Banach lattice approach

Authors: Matheus M. Castro, Jeroen S. W. Lamb, Guillermo Olicón-Méndez, Martin Rasmussen

Abstract: We establish the existence and uniqueness of quasi-stationary and quasi-ergodic measures for almost surely absorbed discrete-time Markov chains under weak conditions. We obtain our results by exploiting Banach lattice properties of transition functions under natural regularity assumptions. We establish the existence and uniqueness of quasi-stationary and quasi-ergodic measures for almost surely absorbed discrete-time Markov chains under weak conditions. We obtain our results by exploiting Banach lattice properties of transition functions under natural regularity assumptions. △ Less

Submitted 28 February, 2024; v1 submitted 26 November, 2021; originally announced November 2021.

Comments: 24 pages, 0 figures

MSC Class: 37H05; 47B65; 60J05

Optimal Decremental Connectivity in Non-Sparse Graphs

Authors: Anders Aaman, Adam Karczmarz, Jakub Łącki, Nikos Parotsidis, Peter M. R. Rasmussen, Mikkel Thorup

Abstract: We present a dynamic algorithm for maintaining the connected and 2-edge-connected components in an undirected graph subject to edge deletions. The algorithm is Monte-Carlo randomized and processes any sequence of edge deletions in $O(m + n \operatorname{polylog} n)$ total time. Interspersed with the deletions, it can answer queries to whether any two given vertices currently belong to the same (2-… ▽ More We present a dynamic algorithm for maintaining the connected and 2-edge-connected components in an undirected graph subject to edge deletions. The algorithm is Monte-Carlo randomized and processes any sequence of edge deletions in $O(m + n \operatorname{polylog} n)$ total time. Interspersed with the deletions, it can answer queries to whether any two given vertices currently belong to the same (2-edge-)connected component in constant time. Our result is based on a general Monte-Carlo randomized reduction from decremental $c$-edge-connectivity to a variant of fully-dynamic $c$-edge-connectivity on a sparse graph. While being Monte-Carlo, our reduction supports a certain final self-check that can be used in Las Vegas algorithms for static problems such as Unique Perfect Matching. For non-sparse graphs with $Ω(n \operatorname{polylog} n)$ edges, our connectivity and $2$-edge-connectivity algorithms handle all deletions in optimal linear total time, using existing algorithms for the respective fully-dynamic problems. This improves upon an $O(m \log (n^2 / m) + n \operatorname{polylog} n)$-time algorithm of Thorup [J.Alg. 1999], which runs in linear time only for graphs with $Ω(n^2)$ edges. Our constant amortized cost for edge deletions in decremental connectivity in non-sparse graphs should be contrasted with an $Ω(\log n/\log\log n)$ worst-case time lower bound in the decremental setting [Alstrup, Thore Husfeldt, FOCS'98] as well as an $Ω(\log n)$ amortized time lower-bound in the fully-dynamic setting [Patrascu and Demaine STOC'04]. △ Less

Submitted 19 November, 2021; v1 submitted 17 November, 2021; originally announced November 2021.

COMAP Early Science: VII. Prospects for CO Intensity Mapping at Reionization

Authors: Patrick C. Breysse, Dongwoo T. Chung, Kieran A. Cleary, Håvard T. Ihle, Hamsa Padmanabhan, Marta B. Silva, J. Richard Bond, Jowita Borowska, Morgan Catha, Sarah E. Church, Delaney A. Dunne, Hans Kristian Eriksen, Marie Kristine Foss, Todd Gaier, Joshua Ott Gundersen, Andrew I. Harris, Richard Hobbs, Laura Keating, James W. Lamb, Charles R. Lawrence, Jonas G. S. Lunde, Norman Murray, Timothy J. Pearson, Liju Philip, Maren Rasmussen , et al. (7 additional authors not shown)

Abstract: We introduce COMAP-EoR, the next generation of the Carbon Monoxide Mapping Array Project aimed at extending CO intensity mapping to the Epoch of Reionization. COMAP-EoR supplements the existing 30 GHz COMAP Pathfinder with two additional 30 GHz instruments and a new 16 GHz receiver. This combination of frequencies will be able to simultaneously map CO(1--0) and CO(2--1) at reionization redshifts (… ▽ More We introduce COMAP-EoR, the next generation of the Carbon Monoxide Mapping Array Project aimed at extending CO intensity mapping to the Epoch of Reionization. COMAP-EoR supplements the existing 30 GHz COMAP Pathfinder with two additional 30 GHz instruments and a new 16 GHz receiver. This combination of frequencies will be able to simultaneously map CO(1--0) and CO(2--1) at reionization redshifts ($z\sim5-8$) in addition to providing a significant boost to the $z\sim3$ sensitivity of the Pathfinder. We examine a set of existing models of the EoR CO signal, and find power spectra spanning several orders of magnitude, highlighting our extreme ignorance about this period of cosmic history and the value of the COMAP-EoR measurement. We carry out the most detailed forecast to date of an intensity mapping cross-correlation, and find that five out of the six models we consider yield signal to noise ratios (S/N) $\gtrsim20$ for COMAP-EoR, with the brightest reaching a S/N above 400. We show that, for these models, COMAP-EoR can make a detailed measurement of the cosmic molecular gas history from $z\sim2-8$, as well as probe the population of faint, star-forming galaxies predicted by these models to be undetectable by traditional surveys. We show that, for the single model that does not predict numerous faint emitters, a COMAP-EoR-type measurement is required to rule out their existence. We briefly explore prospects for a third-generation Expanded Reionization Array (COMAP-ERA) capable of detecting the faintest models and characterizing the brightest signals in extreme detail. △ Less

Submitted 12 November, 2021; v1 submitted 10 November, 2021; originally announced November 2021.

Comments: Paper 7 of 7 in series. 19 pages, 10 figures, to be submitted to ApJ

COMAP Early Science: VI. A First Look at the COMAP Galactic Plane Survey

Authors: Thomas J. Rennie, Stuart E. Harper, Clive Dickinson, Liju Philip, Kieran A. Cleary, Richard J. Bond, Jowita Borowska, Patrick C. Breysse, Morgan Catha, Roke Cepeda-Arroita, Dongwoo T. Chung, Sarah E. Church, Delaney A. Dunne, Hans Kristian Eriksen, Marie Kristine Foss, Todd Gaier, Joshua Ott Gunderson, Andrew I. Harris, Brandon Hensley, Richard Hobbs, Håvard T. Ihle, James W. Lamb, Charles R. Lawrence, Jonas G. S. Lunde, Roberta Paladini , et al. (7 additional authors not shown)

Abstract: We present early results from the COMAP Galactic Plane Survey conducted between June 2019 and April 2021, spanning $20^\circ<\ell<40^\circ$ in Galactic longitude and $|b|<1.\!\!^{\circ}5$ in Galactic latitude with an angular resolution of $4.5^{\prime}$. The full survey will span $\ell \sim 20^{\circ}$- $220^{\circ}$ and will be the first large-scale radio continuum survey at $30$ GHz with sub-deg… ▽ More We present early results from the COMAP Galactic Plane Survey conducted between June 2019 and April 2021, spanning $20^\circ<\ell<40^\circ$ in Galactic longitude and $|b|<1.\!\!^{\circ}5$ in Galactic latitude with an angular resolution of $4.5^{\prime}$. The full survey will span $\ell \sim 20^{\circ}$- $220^{\circ}$ and will be the first large-scale radio continuum survey at $30$ GHz with sub-degree resolution. We present initial results from the first part of the survey, including diffuse emission and spectral energy distributions (SEDs) of HII regions and supernova remnants. Using low and high frequency surveys to constrain free-free and thermal dust emission contributions, we find evidence of excess flux density at $30\,$GHz in six regions that we interpret as anomalous microwave emission. Furthermore we model UCHII contributions using data from the $5\,$GHz CORNISH catalogue and reject this as the cause of the $30\,$GHz excess. Six known supernova remnants (SNR) are detected at $30\,$GHz, and we measure spectral indices consistent with the literature or show evidence of steepening. The flux density of the SNR W44 at $30\,$GHz is consistent with a power-law extrapolation from lower frequencies with no indication of spectral steepening in contrast with recent results from the Sardinia Radio Telescope. We also extract five hydrogen radio recombination lines to map the warm ionized gas, which can be used to estimate electron temperatures or to constrain continuum free-free emission. The full COMAP Galactic plane survey, to be released in 2023/2024, will be an invaluable resource for Galactic astrophysics. △ Less

Submitted 21 March, 2022; v1 submitted 10 November, 2021; originally announced November 2021.

Comments: Paper 6 of 7 in series. 28 pages, 10 figures, submitted to ApJ

COMAP Early Science: V. Constraints and Forecasts at $z \sim 3$

Authors: Dongwoo T. Chung, Patrick C. Breysse, Kieran A. Cleary, Håvard T. Ihle, Hamsa Padmanabhan, Marta B. Silva, J. Richard Bond, Jowita Borowska, Morgan Catha, Sarah E. Church, Delaney A. Dunne, Hans Kristian Eriksen, Marie Kristine Foss, Todd Gaier, Joshua Ott Gundersen, Stuart E. Harper, Andrew I. Harris, Brandon Hensley, Richard Hobbs, Laura C. Keating, Junhan Kim, James W. Lamb, Charles R. Lawrence, Jonas Gahr Sturtzel Lunde, Norman Murray , et al. (12 additional authors not shown)

Abstract: We present the current state of models for the $z\sim3$ carbon monoxide (CO) line-intensity signal targeted by the CO Mapping Array Project (COMAP) Pathfinder in the context of its early science results. Our fiducial model, relating dark matter halo properties to CO luminosities, informs parameter priors with empirical models of the galaxy-halo connection and previous CO(1-0) observations. The Pat… ▽ More We present the current state of models for the $z\sim3$ carbon monoxide (CO) line-intensity signal targeted by the CO Mapping Array Project (COMAP) Pathfinder in the context of its early science results. Our fiducial model, relating dark matter halo properties to CO luminosities, informs parameter priors with empirical models of the galaxy-halo connection and previous CO(1-0) observations. The Pathfinder early science data spanning wavenumbers $k=0.051$-$0.62\,$Mpc$^{-1}$ represent the first direct 3D constraint on the clustering component of the CO(1-0) power spectrum. Our 95% upper limit on the redshift-space clustering amplitude $A_{\rm clust}\lesssim70\,μ$K$^2$ greatly improves on the indirect upper limit of $420\,μ$K$^2$ reported from the CO Power Spectrum Survey (COPSS) measurement at $k\sim1\,$Mpc$^{-1}$. The COMAP limit excludes a subset of models from previous literature, and constrains interpretation of the COPSS results, demonstrating the complementary nature of COMAP and interferometric CO surveys. Using line bias expectations from our priors, we also constrain the squared mean line intensity-bias product, $\langle{Tb}\rangle^2\lesssim50\,μ$K$^2$, and the cosmic molecular gas density, $ρ_\text{H2}<2.5\times10^8\,M_\odot\,$Mpc$^{-3}$ (95% upper limits). Based on early instrument performance and our current CO signal estimates, we forecast that the five-year Pathfinder campaign will detect the CO power spectrum with overall signal-to-noise of 9-17. Between then and now, we also expect to detect the CO-galaxy cross-spectrum using overlapping galaxy survey data, enabling enhanced inferences of cosmic star-formation and galaxy-evolution history. △ Less

Submitted 4 March, 2022; v1 submitted 10 November, 2021; originally announced November 2021.

Comments: Paper 5 of 7 in series. 17 pages + appendix and bibliography (30 pages total); 15 figures, 6 tables; accepted for publication in ApJ; v3 reflects the accepted version with minor changes and additions to text

Journal ref: ApJ, 933, 186 (2022)

COMAP Early Science: IV. Power Spectrum Methodology and Results

Authors: Håvard T. Ihle, Jowita Borowska, Kieran A. Cleary, Hans Kristian Eriksen, Marie K. Foss, Stuart E. Harper, Junhan Kim, Jonas G. S. Lunde, Liju Philip, Maren Rasmussen, Nils-Ole Stutzer, Bade D. Uzgil, Duncan J. Watts, Ingunn Kathrine Wehus, J. Richard Bond, Patrick C. Breysse, Morgan Catha, Sarah E. Church, Dongwoo T. Chung, Clive Dickinson, Delaney A. Dunne, Todd Gaier, Joshua Ott Gundersen, Andrew I. Harris, Richard Hobbs , et al. (8 additional authors not shown)

Abstract: We present the power spectrum methodology used for the first-season COMAP analysis, and assess the quality of the current data set. The main results are derived through the Feed-feed Pseudo-Cross-Spectrum (FPXS) method, which is a robust estimator with respect to both noise modeling errors and experimental systematics. We use effective transfer functions to take into account the effects of instrum… ▽ More We present the power spectrum methodology used for the first-season COMAP analysis, and assess the quality of the current data set. The main results are derived through the Feed-feed Pseudo-Cross-Spectrum (FPXS) method, which is a robust estimator with respect to both noise modeling errors and experimental systematics. We use effective transfer functions to take into account the effects of instrumental beam smoothing and various filter operations applied during the low-level data processing. The power spectra estimated in this way have allowed us to identify a systematic error associated with one of our two scanning strategies, believed to be due to residual ground or atmospheric contamination. We omit these data from our analysis and no longer use this scanning technique for observations. We present the power spectra from our first season of observing and demonstrate that the uncertainties are integrating as expected for uncorrelated noise, with any residual systematics suppressed to a level below the noise. Using the FPXS method, and combining data on scales $k=0.051-0.62 \,\mathrm{Mpc}^{-1}$ we estimate $P_\mathrm{CO}(k) = -2.7 \pm 1.7 \times 10^4μ\textrm{K}^2\mathrm{Mpc}^3$, the first direct 3D constraint on the clustering component of the CO(1-0) power spectrum in the literature. △ Less

Submitted 6 April, 2022; v1 submitted 10 November, 2021; originally announced November 2021.

Comments: Paper 4 of 7 in series. 18 pages, 11 figures, as accepted in ApJ

COMAP Early Science: III. CO Data Processing

Authors: Marie K. Foss, Håvard T. Ihle, Jowita Borowska, Kieran A. Cleary, Hans Kristian Eriksen, Stuart E. Harper, Junhan Kim, James W. Lamb, Jonas G. S. Lunde, Liju Philip, Maren Rasmussen, Nils-Ole Stutzer, Bade D. Uzgil, Duncan J. Watts, Ingunn K. Wehus, David P. Woody, J. Richard Bond, Patrick C. Breysse, Morgan Catha, Sarah E. Church, Dongwoo T. Chung, Clive Dickinson, Delaney A. Dunne, Todd Gaier, Joshua Ott Gundersen , et al. (8 additional authors not shown)

Abstract: We describe the first season COMAP analysis pipeline that converts raw detector readouts to calibrated sky maps. This pipeline implements four main steps: gain calibration, filtering, data selection, and map-making. Absolute gain calibration relies on a combination of instrumental and astrophysical sources, while relative gain calibration exploits real-time total-power variations. High efficiency… ▽ More We describe the first season COMAP analysis pipeline that converts raw detector readouts to calibrated sky maps. This pipeline implements four main steps: gain calibration, filtering, data selection, and map-making. Absolute gain calibration relies on a combination of instrumental and astrophysical sources, while relative gain calibration exploits real-time total-power variations. High efficiency filtering is achieved through spectroscopic common-mode rejection within and across receivers, resulting in nearly uncorrelated white noise within single-frequency channels. Consequently, near-optimal but biased maps are produced by binning the filtered time stream into pixelized maps; the corresponding signal bias transfer function is estimated through simulations. Data selection is performed automatically through a series of goodness-of-fit statistics, including $χ^2$ and multi-scale correlation tests. Applying this pipeline to the first-season COMAP data, we produce a dataset with very low levels of correlated noise. We find that one of our two scanning strategies (the Lissajous type) is sensitive to residual instrumental systematics. As a result, we no longer use this type of scan and exclude data taken this way from our Season 1 power spectrum estimates. We perform a careful analysis of our data processing and observing efficiencies and take account of planned improvements to estimate our future performance. Power spectrum results derived from the first-season COMAP maps are presented and discussed in companion papers. △ Less

Submitted 30 November, 2021; v1 submitted 10 November, 2021; originally announced November 2021.

Comments: Paper 3 of 7 in series. 26 pages, 23 figures, submitted to ApJ

COMAP Early Science: II. Pathfinder Instrument

Authors: James W. Lamb, Kieran A. Cleary, David P. Woody, Morgan Catha, Dongwoo T. Chung, Joshua Ott Gundersen, Stuart E. Harper, Andrew I. Harris, Richard Hobbs, Håvard T. Ihle, Jonathon Kocz, Timothy J. Pearson, Liju Philip, Travis W. Powell, Lilian Basoalto, J. Richard Bond, Jowita Borowska, Patrick C. Breysse, Sarah E. Church, Clive Dickinson, Delaney A. Dunne, Hans Kristian Eriksen, Marie Kristine Foss, Todd Gaier, Junhan Kim , et al. (10 additional authors not shown)

Abstract: Line intensity mapping (LIM) is a new technique for tracing the global properties of galaxies over cosmic time. Detection of the very faint signals from redshifted carbon monoxide (CO), a tracer of star formation, pushes the limits of what is feasible with a total-power instrument. The CO Mapping Project (COMAP) Pathfinder is a first-generation instrument aiming to prove the concept and develop th… ▽ More Line intensity mapping (LIM) is a new technique for tracing the global properties of galaxies over cosmic time. Detection of the very faint signals from redshifted carbon monoxide (CO), a tracer of star formation, pushes the limits of what is feasible with a total-power instrument. The CO Mapping Project (COMAP) Pathfinder is a first-generation instrument aiming to prove the concept and develop the technology for future experiments, as well as delivering early science products. With 19 receiver channels in a hexagonal focal plane arrangement on a 10.4 m antenna, and an instantaneous 26-34 GHz frequency range with 2 MHz resolution, it is ideally suited to measuring CO($J$=1-0) from $z\sim3$. In this paper we discuss strategies for designing and building the Pathfinder and the challenges that were encountered. The design of the instrument prioritized LIM requirements over those of ancillary science. After a couple of years of operation, the instrument is well understood, and the first year of data is already yielding useful science results. Experience with this Pathfinder will drive the design of the next generations of experiments. △ Less

Submitted 29 November, 2021; v1 submitted 10 November, 2021; originally announced November 2021.

Comments: Paper 2 of 7 in series. 27 pages, 28 figures, submitted to ApJ

COMAP Early Science: I. Overview

Authors: Kieran A. Cleary, Jowita Borowska, Patrick C. Breysse, Morgan Catha, Dongwoo T. Chung, Sarah E. Church, Clive Dickinson, Hans Kristian Eriksen, Marie Kristine Foss, Joshua Ott Gundersen, Stuart E. Harper, Andrew I. Harris, Richard Hobbs, Håvard, T. Ihle, Junhan Kim, Jonathon Kocz, James W. Lamb, Jonas G. S. Lunde, Hamsa Padmanabhan, Timothy J. Pearson, Liju Philip, Travis W. Powell, Maren Rasmussen, Anthony C. S. Readhead , et al. (18 additional authors not shown)

Abstract: The CO Mapping Array Project (COMAP) aims to use line intensity mapping of carbon monoxide (CO) to trace the distribution and global properties of galaxies over cosmic time, back to the Epoch of Reionization (EoR). To validate the technologies and techniques needed for this goal, a Pathfinder instrument has been constructed and fielded. Sensitive to CO(1-0) emission from $z=2.4$-$3.4$ and a fainte… ▽ More The CO Mapping Array Project (COMAP) aims to use line intensity mapping of carbon monoxide (CO) to trace the distribution and global properties of galaxies over cosmic time, back to the Epoch of Reionization (EoR). To validate the technologies and techniques needed for this goal, a Pathfinder instrument has been constructed and fielded. Sensitive to CO(1-0) emission from $z=2.4$-$3.4$ and a fainter contribution from CO(2-1) at $z=6$-8, the Pathfinder is surveying $12$ deg$^2$ in a 5-year observing campaign to detect the CO signal from $z\sim3$. Using data from the first 13 months of observing, we estimate $P_\mathrm{CO}(k) = -2.7 \pm 1.7 \times 10^4μ\mathrm{K}^2 \mathrm{Mpc}^3$ on scales $k=0.051-0.62 \mathrm{Mpc}^{-1}$ - the first direct 3D constraint on the clustering component of the CO(1-0) power spectrum. Based on these observations alone, we obtain a constraint on the amplitude of the clustering component (the squared mean CO line temperature-bias product) of $\langle Tb\rangle^2<49$ $μ$K$^2$ - nearly an order-of-magnitude improvement on the previous best measurement. These constraints allow us to rule out two models from the literature. We forecast a detection of the power spectrum after 5 years with signal-to-noise ratio (S/N) 9-17. Cross-correlation with an overlapping galaxy survey will yield a detection of the CO-galaxy power spectrum with S/N of 19. We are also conducting a 30 GHz survey of the Galactic plane and present a preliminary map. Looking to the future of COMAP, we examine the prospects for future phases of the experiment to detect and characterize the CO signal from the EoR. △ Less

Submitted 29 November, 2021; v1 submitted 10 November, 2021; originally announced November 2021.

Comments: Paper 1 of 7 in series. 18 pages, 16 figures, submitted to ApJ

Rate-induced tipping and saddle-node bifurcation for a class of quadratic differential equations with nonautonomous asymptotic dynamics

Authors: Iacopo P. Longo, Carmen Núñez, Rafael Obaya, Martin Rasmussen

Abstract: An in-depth analysis of nonautonomous bifurcations of saddle-node type for scalar differential equations $x'=-x^2+q(t)\,x+p(t)$, where $q\colon\mathbb{R}\to\mathbb{R}$ and $p\colon\mathbb{R}\to\mathbb{R}$ are bounded and uniformly continuous, is fundamental to explain the absence or occurrence of rate-induced tipping for the differential equation $y' =(y-(2/π)\arctan(ct))^2+p(t)$ as the rate $c$ v… ▽ More An in-depth analysis of nonautonomous bifurcations of saddle-node type for scalar differential equations $x'=-x^2+q(t)\,x+p(t)$, where $q\colon\mathbb{R}\to\mathbb{R}$ and $p\colon\mathbb{R}\to\mathbb{R}$ are bounded and uniformly continuous, is fundamental to explain the absence or occurrence of rate-induced tipping for the differential equation $y' =(y-(2/π)\arctan(ct))^2+p(t)$ as the rate $c$ varies on $[0,\infty)$. A classical attractor-repeller pair, whose existence for $c=0$ is assumed, may persist for any $c>0$, or disappear for a certain critical rate $c=c_0$, giving rise to rate-induced tipping. A suitable example demonstrates that one can have more than one critical rate, and the existence of the classical attractor-repeller pair may return when $c$ increases. △ Less

Submitted 6 October, 2021; originally announced October 2021.

Comments: 40 pages, 10 figures

Journal ref: SIAM Journal on Applied Dynamical Systems, 20(1) (2021), 500-540

On Boundaries of $\varepsilon$-neighbourhoods of Planar Sets, Part II: Global Structure and Curvature

Authors: Jeroen S. W. Lamb, Martin Rasmussen, Kalle Timperi

Abstract: We study the global topological structure and smoothness of the boundaries of $\varepsilon$-neighbourhoods $E_\varepsilon = \{x \in \mathbb{R}^2 \, : \, \textrm{dist}(x, E) \leq \varepsilon \}$ of planar sets $E \subset \mathbb{R}^2$. We show that for a compact set $E$ and $\varepsilon > 0$ the boundary $\partial E_\varepsilon$ can be expressed as a disjoint union of an at most countably infinite… ▽ More We study the global topological structure and smoothness of the boundaries of $\varepsilon$-neighbourhoods $E_\varepsilon = \{x \in \mathbb{R}^2 \, : \, \textrm{dist}(x, E) \leq \varepsilon \}$ of planar sets $E \subset \mathbb{R}^2$. We show that for a compact set $E$ and $\varepsilon > 0$ the boundary $\partial E_\varepsilon$ can be expressed as a disjoint union of an at most countably infinite union of Jordan curves and a possibly uncountable, totally disconnected set of singularities. We also show that curvature is defined almost everywhere on the Jordan curve subsets of the boundary. △ Less

Submitted 15 July, 2021; originally announced July 2021.

Comments: 31 pages, 11 figures

MSC Class: 51F30 (Primary) 57K20; 54C50; 51M15; 58C06 (Secondary)

Learning Safe and Optimal Control Strategies for Storm Water Detention Ponds

Authors: Martijn A. Goorden, Kim G. Larsen, Jesper E. Nielsen, Thomas D. Nielsen, Michael R. Rasmussen, Jiri Srba

Abstract: Storm water detention ponds are used to manage the discharge of rainfall runoff from urban areas to nearby streams. Their purpose is to reduce the hydraulic impact and sediment loads of the receiving waters. Detention ponds are currently designed based on static controls: the output flow of a pond is capped at a fixed value. This is not optimal with respect to the current infrastructure capacity a… ▽ More Storm water detention ponds are used to manage the discharge of rainfall runoff from urban areas to nearby streams. Their purpose is to reduce the hydraulic impact and sediment loads of the receiving waters. Detention ponds are currently designed based on static controls: the output flow of a pond is capped at a fixed value. This is not optimal with respect to the current infrastructure capacity and for some detention ponds it might even violate current regulations set by the European Water Framework Directive. We apply formal methods to synthesize (i.e., derive automatically) a safe and optimal active controller. We model the storm water detention pond, including the urban catchment area and the rain forecasts, as a hybrid Markov decision process. Subsequently, we use the tool Uppaal Stratego to synthesize a control strategy minimizing the cost related to pollution (optimality) while guaranteeing no emergency overflow of the detention pond (safety). Simulation results for an existing pond show that Uppaal Stratego can learn optimal strategies that prevent emergency overflows, where the current static control is not always able to prevent it. At the same time, our approach can improve sedimentation during low rain periods. △ Less

Submitted 26 April, 2021; originally announced April 2021.

Comments: Accepted for the IFAC Conference on Analysis and Design of Hybrid Systems 2021. Full paper version

On Boundaries of $\varepsilon$-neighbourhoods of Planar Sets, Part I: Singularities

Authors: Jeroen S. W. Lamb, Martin Rasmussen, Kalle Timperi

Abstract: We study geometric and topological properties of singularities on the boundaries of $\varepsilon$-neighbourhoods $E_\varepsilon = \{x \in \mathbb{R}^2 : \textrm{dist}(x, E) \leq \varepsilon \}$ of planar sets $E \subset \mathbb{R}^2$. We develop a novel technique for analysing the boundary and obtain, for a compact set $E$ and $\varepsilon > 0$, a classification of singularities (i.e. non-smooth p… ▽ More We study geometric and topological properties of singularities on the boundaries of $\varepsilon$-neighbourhoods $E_\varepsilon = \{x \in \mathbb{R}^2 : \textrm{dist}(x, E) \leq \varepsilon \}$ of planar sets $E \subset \mathbb{R}^2$. We develop a novel technique for analysing the boundary and obtain, for a compact set $E$ and $\varepsilon > 0$, a classification of singularities (i.e. non-smooth points) on $\partial E_\varepsilon$ into eight categories. We also show that the set of singularities is either countable or the disjoint union of a countable set and a closed, totally disconnected, nowhere dense set. △ Less

Submitted 25 December, 2020; originally announced December 2020.

Comments: 40 pages, 9 figures

MSC Class: 51F30 (Primary) 57K20; 54C50; 51M15; 58C06 (Secondary)

Occupancy times for time-dependent stage-structured models

Authors: George Chappelle, Alan Hastings, Martin Rasmussen

Abstract: During their lifetimes, individuals in populations pass through different states, and the notion of an occupancy time describes the amount of time an individual spends in a given set of states. Questions related to this idea were studied in a recent paper by Roth and Caswell for cases where the environmental conditions are constant. However, it is truly important to consider the case where environ… ▽ More During their lifetimes, individuals in populations pass through different states, and the notion of an occupancy time describes the amount of time an individual spends in a given set of states. Questions related to this idea were studied in a recent paper by Roth and Caswell for cases where the environmental conditions are constant. However, it is truly important to consider the case where environments are changing randomly or in directional way through time, so the transition probabilities between different states change over time, motivating the use of time-dependent stage-structured models. Using absorbing inhomogenous Markov chains and the discrete-time McKendrick--von F{ö}rster equation, we derive explicit formulas for the occupancy time, its expectation, and its higher-order moments for stage-structured models with time-dependent transition rates. We apply our approach to study a time-dependent model of the Southern Fulmar, and obtain insights into how the number of breeding attempts depends on external conditions that vary through time. △ Less

Submitted 30 November, 2020; originally announced December 2020.

Comments: 16 pages, 10 figures

Tiling with Squares and Packing Dominos in Polynomial Time

Authors: Anders Aamand, Mikkel Abrahamsen, Thomas D. Ahle, Peter M. R. Rasmussen

Abstract: A polyomino is a polygonal region with axis parallel edges and corners of integral coordinates, which may have holes. In this paper, we consider planar tiling and packing problems with polyomino pieces and a polyomino container $P$. We give two polynomial time algorithms, one for deciding if $P$ can be tiled with $k\times k$ squares for any fixed $k$ which can be part of the input (that is, decidi… ▽ More A polyomino is a polygonal region with axis parallel edges and corners of integral coordinates, which may have holes. In this paper, we consider planar tiling and packing problems with polyomino pieces and a polyomino container $P$. We give two polynomial time algorithms, one for deciding if $P$ can be tiled with $k\times k$ squares for any fixed $k$ which can be part of the input (that is, deciding if $P$ is the union of a set of non-overlapping $k\times k$ squares) and one for packing $P$ with a maximum number of non-overlapping and axis-parallel $2\times 1$ dominos, allowing rotations by $90^\circ$. As packing is more general than tiling, the latter algorithm can also be used to decide if $P$ can be tiled by $2\times 1$ dominos. These are classical problems with important applications in VLSI design, and the related problem of finding a maximum packing of $2\times 2$ squares is known to be NP-Hard [J. Algorithms 1990]. For our three problems there are known pseudo-polynomial time algorithms, that is, algorithms with running times polynomial in the area of $P$. However, the standard, compact way to represent a polygon is by listing the coordinates of the corners in binary. We use this representation, and thus present the first polynomial time algorithms for the problems. Concretely, we give a simple $O(n\log n)$ algorithm for tiling with squares, and a more involved $O(n^3\,\text{polylog}\, n)$ algorithm for packing and tiling with dominos, where $n$ is the number of corners of $P$. △ Less

Submitted 9 August, 2021; v1 submitted 22 November, 2020; originally announced November 2020.

Comments: Compared to the first version, running times for domino packing have been improved and a simpler algorithm has been described

Modeling all alternative solutions for highly renewable energy systems

Authors: Tim T. Pedersen, Marta Victoria, Morten G. Rasmussen, Gorm B. Andresen

Abstract: As the world is transitioning towards highly renewable energy systems, advanced tools are needed to analyze such complex networks. Energy system design is, however, challenged by real-world objective functions consisting of a blurry mix of technical and socioeconomic agendas, with limitations that cannot always be clearly stated. As a result, it is highly likely that solutions which are techno-eco… ▽ More As the world is transitioning towards highly renewable energy systems, advanced tools are needed to analyze such complex networks. Energy system design is, however, challenged by real-world objective functions consisting of a blurry mix of technical and socioeconomic agendas, with limitations that cannot always be clearly stated. As a result, it is highly likely that solutions which are techno-economically suboptimal will be preferable. Here, we present a method capable of determining the continuum containing all techno-economically near-optimal solutions, moving the field of energy system modeling from discrete solutions to a new era where continuous solution ranges are available. The presented method is applied to study a range of technical and socioeconomic metrics on a model of the European electricity system. The near-optimal region is found to be relatively flat allowing for solutions that are slightly more expensive than the optimum but better in terms of equality, land use, and implementation time. △ Less

Submitted 29 June, 2021; v1 submitted 2 October, 2020; originally announced October 2020.

Comments: 25 pages, 7 figures, also available as preprint at: https://papers.ssrn.com/sol3/papers.cfm?abstract_id=3682045

Journal ref: Tim T. Pedersen, Marta Victoria, Morten G. Rasmussen, Gorm B. Andresen, Modeling all alternative solutions for highly renewable energy systems, Energy, Volume 234, 2021, 121294, ISSN 0360-5442,

Support of Closed Walks and Second Eigenvalue Multiplicity of the Normalized Adjacency Matrix

Authors: Theo McKenzie, Peter M. R. Rasmussen, Nikhil Srivastava

Abstract: We show that the multiplicity of the second normalized adjacency matrix eigenvalue of any connected graph of maximum degree $Δ$ is bounded by $O(n Δ^{7/5}/\log^{1/5-o(1)}n)$ for any $Δ$, and by $O(n\log^{1/2}d/\log^{1/4-o(1)}n)$ for simple $d$-regular graphs when $d\ge \log^{1/4}n$. In fact, the same bounds hold for the number of eigenvalues in any interval of width $λ_2/\log_Δ^{1-o(1)}n$ containi… ▽ More We show that the multiplicity of the second normalized adjacency matrix eigenvalue of any connected graph of maximum degree $Δ$ is bounded by $O(n Δ^{7/5}/\log^{1/5-o(1)}n)$ for any $Δ$, and by $O(n\log^{1/2}d/\log^{1/4-o(1)}n)$ for simple $d$-regular graphs when $d\ge \log^{1/4}n$. In fact, the same bounds hold for the number of eigenvalues in any interval of width $λ_2/\log_Δ^{1-o(1)}n$ containing the second eigenvalue $λ_2$. The main ingredient in the proof is a polynomial (in $k$) lower bound on the typical support of a closed random walk of length $2k$ in any connected graph, which in turn relies on new lower bounds for the entries of the Perron eigenvector of submatrices of the normalized adjacency matrix. △ Less

Submitted 16 June, 2023; v1 submitted 24 July, 2020; originally announced July 2020.

Comments: A previous version of this paper proved the main result for d-regular graphs. The current version proves a more general result for the normalized adjacency matrix of bounded degree graphs. New version fixes an incorrect citation in the proof of Proposition 5.2. 24pp, 3 figures

Journal ref: STOC 2021

A Speech Enhancement Algorithm based on Non-negative Hidden Markov Model and Kullback-Leibler Divergence

Authors: Yang Xiang, Liming Shi, Jesper Lisby Højvang, Morten Højfeldt Rasmussen, Mads Græsbøll Christensen

Abstract: In this paper, we propose a novel supervised single-channel speech enhancement method combing the the Kullback-Leibler divergence-based non-negative matrix factorization (NMF) and hidden Markov model (NMF-HMM). With the application of HMM, the temporal dynamics information of speech signals can be taken into account. In the training stage, the sum of Poisson, leading to the KL divergence measure,… ▽ More In this paper, we propose a novel supervised single-channel speech enhancement method combing the the Kullback-Leibler divergence-based non-negative matrix factorization (NMF) and hidden Markov model (NMF-HMM). With the application of HMM, the temporal dynamics information of speech signals can be taken into account. In the training stage, the sum of Poisson, leading to the KL divergence measure, is used as the observation model for each state of HMM. This ensures that a computationally efficient multiplicative update can be used for the parameter update of the proposed model. In the online enhancement stage, we propose a novel minimum mean-square error (MMSE) estimator for the proposed NMF-HMM. This estimator can be implemented using parallel computing, saving the time complexity. The performance of the proposed algorithm is verified by objective measures. The experimental results show that the proposed strategy achieves better speech enhancement performance than state-of-the-art speech enhancement methods. More specifically, compared with the traditional NMF-based speech enhancement methods, our proposed algorithm achieves a 5\% improvement for short-time objective intelligibility (STOI) and 0.18 improvement for perceptual evaluation of speech quality (PESQ). △ Less

Submitted 30 June, 2020; originally announced June 2020.

Quasi-ergodic limits for finite absorbing Markov chains

Authors: Fritz Colonius, Martin Rasmussen

Abstract: We present formulas for quasi-ergodic limits of finite absorbing Markov chains. Since the irreducible case has been solved in 1965 by Darroch and Seneta, we focus on the reducible case, and our results are based on a very precise asymptotic analysis of the (exponential and polynomial) growth behaviour along admissible paths We present formulas for quasi-ergodic limits of finite absorbing Markov chains. Since the irreducible case has been solved in 1965 by Darroch and Seneta, we focus on the reducible case, and our results are based on a very precise asymptotic analysis of the (exponential and polynomial) growth behaviour along admissible paths △ Less

Submitted 26 May, 2020; originally announced May 2020.

Comments: 30 pages

MSC Class: 15B51; 37A25; 60J10

Journal ref: Linear Algebra and its Applications 609 (2021), pp. 253-288

No Repetition: Fast Streaming with Highly Concentrated Hashing

Authors: Anders Aamand, Debarati Das, Evangelos Kipouridis, Jakob B. T. Knudsen, Peter M. R. Rasmussen, Mikkel Thorup

Abstract: To get estimators that work within a certain error bound with high probability, a common strategy is to design one that works with constant probability, and then boost the probability using independent repetitions. Important examples of this approach are small space algorithms for estimating the number of distinct elements in a stream, or estimating the set similarity between large sets. Using sta… ▽ More To get estimators that work within a certain error bound with high probability, a common strategy is to design one that works with constant probability, and then boost the probability using independent repetitions. Important examples of this approach are small space algorithms for estimating the number of distinct elements in a stream, or estimating the set similarity between large sets. Using standard strongly universal hashing to process each element, we get a sketch based estimator where the probability of a too large error is, say, 1/4. By performing $r$ independent repetitions and taking the median of the estimators, the error probability falls exponentially in $r$. However, running $r$ independent experiments increases the processing time by a factor $r$. Here we make the point that if we have a hash function with strong concentration bounds, then we get the same high probability bounds without any need for repetitions. Instead of $r$ independent sketches, we have a single sketch that is $r$ times bigger, so the total space is the same. However, we only apply a single hash function, so we save a factor $r$ in time, and the overall algorithms just get simpler. Fast practical hash functions with strong concentration bounds were recently proposed by Aamand em et al. (to appear in STOC 2020). Using their hashing schemes, the algorithms thus become very fast and practical, suitable for online processing of high volume data streams. △ Less

Submitted 2 April, 2020; originally announced April 2020.

Comments: 10 pages

Automatic image-based identification and biomass estimation of invertebrates

Authors: Johanna Ärje, Claus Melvad, Mads Rosenhøj Jeppesen, Sigurd Agerskov Madsen, Jenni Raitoharju, Maria Strandgård Rasmussen, Alexandros Iosifidis, Ville Tirronen, Kristian Meissner, Moncef Gabbouj, Toke Thomas Høye

Abstract: Understanding how biological communities respond to environmental changes is a key challenge in ecology and ecosystem management. The apparent decline of insect populations necessitates more biomonitoring but the time-consuming sorting and identification of taxa pose strong limitations on how many insect samples can be processed. In turn, this affects the scale of efforts to map invertebrate diver… ▽ More Understanding how biological communities respond to environmental changes is a key challenge in ecology and ecosystem management. The apparent decline of insect populations necessitates more biomonitoring but the time-consuming sorting and identification of taxa pose strong limitations on how many insect samples can be processed. In turn, this affects the scale of efforts to map invertebrate diversity altogether. Given recent advances in computer vision, we propose to replace the standard manual approach of human expert-based sorting and identification with an automatic image-based technology. We describe a robot-enabled image-based identification machine, which can automate the process of invertebrate identification, biomass estimation and sample sorting. We use the imaging device to generate a comprehensive image database of terrestrial arthropod species. We use this database to test the classification accuracy i.e. how well the species identity of a specimen can be predicted from images taken by the machine. We also test sensitivity of the classification accuracy to the camera settings (aperture and exposure time) in order to move forward with the best possible image quality. We use state-of-the-art Resnet-50 and InceptionV3 CNNs for the classification task. The results for the initial dataset are very promising ($\overline{ACC}=0.980$). The system is general and can easily be used for other groups of invertebrates as well. As such, our results pave the way for generating more data on spatial and temporal variation in invertebrate abundance, diversity and biomass. △ Less

Submitted 5 February, 2020; originally announced February 2020.

Cross-product Penalized Component Analysis (XCAN)

Authors: José Camacho, Evrim Acar, Morten A. Rasmussen, Rasmus Bro

Abstract: Matrix factorization methods are extensively employed to understand complex data. In this paper, we introduce the cross-product penalized component analysis (XCAN), a sparse matrix factorization based on the optimization of a loss function that allows a trade-off between variance maximization and structural preservation. The approach is based on previous developments, notably (i) the Sparse Princi… ▽ More Matrix factorization methods are extensively employed to understand complex data. In this paper, we introduce the cross-product penalized component analysis (XCAN), a sparse matrix factorization based on the optimization of a loss function that allows a trade-off between variance maximization and structural preservation. The approach is based on previous developments, notably (i) the Sparse Principal Component Analysis (SPCA) framework based on the LASSO, (ii) extensions of SPCA to constrain both modes of the factorization, like co-clustering or the Penalized Matrix Decomposition (PMD), and (iii) the Group-wise Principal Component Analysis (GPCA) method. The result is a flexible modeling approach that can be used for data exploration in a large variety of problems. We demonstrate its use with applications from different disciplines. △ Less

Submitted 28 June, 2019; originally announced July 2019.

Journal ref: Chemometrics and Intelligent Laboratory Systems, 2020, 203: 104038-

Fast hashing with Strong Concentration Bounds

Authors: Anders Aamand, Jakob B. T. Knudsen, Mathias B. T. Knudsen, Peter M. R. Rasmussen, Mikkel Thorup

Abstract: Previous work on tabulation hashing by Patrascu and Thorup from STOC'11 on simple tabulation and from SODA'13 on twisted tabulation offered Chernoff-style concentration bounds on hash based sums, e.g., the number of balls/keys hashing to a given bin, but under some quite severe restrictions on the expected values of these sums. The basic idea in tabulation hashing is to view a key as consisting of… ▽ More Previous work on tabulation hashing by Patrascu and Thorup from STOC'11 on simple tabulation and from SODA'13 on twisted tabulation offered Chernoff-style concentration bounds on hash based sums, e.g., the number of balls/keys hashing to a given bin, but under some quite severe restrictions on the expected values of these sums. The basic idea in tabulation hashing is to view a key as consisting of $c=O(1)$ characters, e.g., a 64-bit key as $c=8$ characters of 8-bits. The character domain $Σ$ should be small enough that character tables of size $|Σ|$ fit in fast cache. The schemes then use $O(1)$ tables of this size, so the space of tabulation hashing is $O(|Σ|)$. However, the concentration bounds by Patrascu and Thorup only apply if the expected sums are $\ll |Σ|$. To see the problem, consider the very simple case where we use tabulation hashing to throw $n$ balls into $m$ bins and want to analyse the number of balls in a given bin. With their concentration bounds, we are fine if $n=m$, for then the expected value is $1$. However, if $m=2$, as when tossing $n$ unbiased coins, the expected value $n/2$ is $\gg |Σ|$ for large data sets, e.g., data sets that do not fit in fast cache. To handle expectations that go beyond the limits of our small space, we need a much more advanced analysis of simple tabulation, plus a new tabulation technique that we call \emph{tabulation-permutation} hashing which is at most twice as slow as simple tabulation. No other hashing scheme of comparable speed offers similar Chernoff-style concentration bounds. △ Less

Submitted 10 August, 2020; v1 submitted 1 May, 2019; originally announced May 2019.

Comments: 54 pages, 3 figures. An extended abstract appeared at the 52nd Annual ACM Symposium on Theory of Computing (STOC20)

A Girsanov approach to slow parameterizing manifolds in the presence of noise

Authors: Mickaël D. Chekroun, Jeroen S. W. Lamb, Christian J. Pangerl, Martin Rasmussen

Abstract: We consider a three-dimensional slow-fast system arising in fluid dynamics with quadratic nonlinearity and additive noise. The associated deterministic system of this stochastic differential equation (SDE) exhibits a periodic orbit and a slow manifold. We show that in presence of noise, the deterministic slow manifold can be viewed as an approximate parameterization of the fast variable of the SDE… ▽ More We consider a three-dimensional slow-fast system arising in fluid dynamics with quadratic nonlinearity and additive noise. The associated deterministic system of this stochastic differential equation (SDE) exhibits a periodic orbit and a slow manifold. We show that in presence of noise, the deterministic slow manifold can be viewed as an approximate parameterization of the fast variable of the SDE in terms of the slow variables, for certain parameter regimes. We exploit this fact to obtain a two dimensional reduced model from the original stochastic system, which results into a Hopf normal form with additive noise. Both, the original as well as the reduced system admit ergodic invariant measures describing their respective long-time behaviour. It is then shown that for a suitable Wasserstein metric on a subset of the space of probability measures on the phase space, the discrepancy between the marginals along the radial component of each invariant measure is controlled by a parameterization defect measuring the quality of the parameterization. An important technical tool to arrive at this result is Girsanov's theorem that allows us to derive such error estimates in presence of an oscillatory instability. This approach is finally extended to parameter regimes for which the variable to parameterize is no longer evolving on a faster timescale than that of the resolved variables. There also, error estimates involving the Wasserstein metric are derived but this time for reduced systems obtained from stochastic parameterizing manifolds involving path-dependent coefficients to cope with such challenging regimes. △ Less

Submitted 29 November, 2022; v1 submitted 20 March, 2019; originally announced March 2019.

Comments: 49 pages, 6 figures

Classifying Convex Bodies by their Contact and Intersection Graphs

Authors: Anders Aamand, Mikkel Abrahamsen, Jakob Bæk Tejs Knudsen, Peter Michael Reichstein Rasmussen

Abstract: Suppose that $A$ is a convex body in the plane and that $A_1,\dots,A_n$ are translates of $A$. Such translates give rise to an intersection graph of $A$, $G=(V,E)$, with vertices $V=\{1,\dots,n\}$ and edges $E=\{uv\mid A_u\cap A_v\neq \emptyset\}$. The subgraph $G'=(V, E')$ satisfying that $E'\subset E$ is the set of edges $uv$ for which the interiors of $A_u$ and $A_v$ are disjoint is a unit dist… ▽ More Suppose that $A$ is a convex body in the plane and that $A_1,\dots,A_n$ are translates of $A$. Such translates give rise to an intersection graph of $A$, $G=(V,E)$, with vertices $V=\{1,\dots,n\}$ and edges $E=\{uv\mid A_u\cap A_v\neq \emptyset\}$. The subgraph $G'=(V, E')$ satisfying that $E'\subset E$ is the set of edges $uv$ for which the interiors of $A_u$ and $A_v$ are disjoint is a unit distance graph of $A$. If furthermore $G'=G$, i.e., if the interiors of $A_u$ and $A_v$ are disjoint whenever $u\neq v$, then $G$ is a contact graph of $A$. In this paper we study which pairs of convex bodies have the same contact, unit distance, or intersection graphs. We say that two convex bodies $A$ and $B$ are equivalent if there exists a linear transformation $B'$ of $B$ such that for any slope, the longest line segments with that slope contained in $A$ and $B'$, respectively, are equally long. For a broad class of convex bodies, including all strictly convex bodies and linear transformations of regular polygons, we show that the contact graphs of $A$ and $B$ are the same if and only if $A$ and $B$ are equivalent. We prove the same statement for unit distance and intersection graphs. △ Less

Submitted 5 February, 2019; originally announced February 2019.

Comments: 19 pages, 7 figures

MSC Class: 52C05

Dynamical characterization of stochastic bifurcations in a random logistic map

Authors: Yuzuru Sato, Thai Son Doan, Jeroen S. W. Lamb, Martin Rasmussen

Abstract: The emergence of noise-induced chaos in a random logistic map with bounded noise is understood as a two-step process consisting of a topological bifurcation flagged by a zero-crossing point of the supremum of the dichotomy spectrum and a subsequent dynamical bifurcation to a random strange attractor flagged by a zero crossing point of the Lyapunov exponent. The associated three consecutive dynamic… ▽ More The emergence of noise-induced chaos in a random logistic map with bounded noise is understood as a two-step process consisting of a topological bifurcation flagged by a zero-crossing point of the supremum of the dichotomy spectrum and a subsequent dynamical bifurcation to a random strange attractor flagged by a zero crossing point of the Lyapunov exponent. The associated three consecutive dynamical phases are characterized as a random periodic attractor, a random point attractor, and a random strange attractor, respectively. The first phase has a negative dichotomy spectrum reflecting uniform attraction to the random periodic attractor. The second phase no longer has a negative dichotomy spectrum - and the random point attractor is not uniformly attractive - but it retains a negative Lyapunov exponent reflecting the aggregate asymptotic contractive behaviour. For practical purposes, the extrema of the dichotomy spectrum equal that of the support of the spectrum of the finite-time Lyapunov exponents. We present detailed numerical results from various dynamical viewpoints, illustrating the dynamical characterisation of the three different phases. △ Less

Submitted 9 November, 2018; originally announced November 2018.

Expander Graphs are Non-Malleable Codes

Authors: Peter M. R. Rasmussen, Amit Sahai

Abstract: Any $d$-regular graph on $n$ vertices with spectral expansion $λ$ satisfying $n = Ω(d^3\log(d)/λ)$ yields a $O\left(\frac{λ^{3/2}}{d}\right)$-non-malleable code for single-bit messages in the split-state model. Any $d$-regular graph on $n$ vertices with spectral expansion $λ$ satisfying $n = Ω(d^3\log(d)/λ)$ yields a $O\left(\frac{λ^{3/2}}{d}\right)$-non-malleable code for single-bit messages in the split-state model. △ Less

Submitted 20 March, 2019; v1 submitted 28 September, 2018; originally announced October 2018.

Comments: 10 pages Resubmitted with revised introduction and acknowledgement

Ion Impact Induced Ultrafast Electron Dynamics in Correlated Materials and Finite Graphene Clusters

Authors: M. Bonitz, K. Balzer, N. Schlünzen, M. Rasmussen, J. -P. Joost

Abstract: Strongly correlated systems of fermions have an interesting phase diagram arising from the Hubbard gap. Excitation across the gap leads to the formation of doubly occupied lattice sites (doublons). This state offers interesting electronic and optical properties. Moreover, when the system is driven out of equilibrium interesting collective dynamics may arise that are related to the spatial propagat… ▽ More Strongly correlated systems of fermions have an interesting phase diagram arising from the Hubbard gap. Excitation across the gap leads to the formation of doubly occupied lattice sites (doublons). This state offers interesting electronic and optical properties. Moreover, when the system is driven out of equilibrium interesting collective dynamics may arise that are related to the spatial propagation of doublons. Here, a novel mechanism that was recently proposed by us [Balzer \textit{et al.}, submitted for publication] is verified by exact diagonalization and nonequilibrium Green functions (NEGF) simulations---fermionic doublon creation by the impact of energetic ions. We report the formation of a nonequilibrium steady state with homogeneous doublon distribution. A physically intuitive picture is given in terms of an analytical model for a two-site system where the doublon formation is explained in terms of a two-fold passage of an avoided crossing (Landau-Zener picture). The effect should be particularly important for strongly correlated finite systems, such as graphene nanoribbons, and directly observable with fermionic atoms in optical lattices. We demonstrate that doublon formation and propagation in correlated lattice systems can be accurately simulated with NEGF. In addition to two-time results we present single-time results within the generalized Kadanoff-Baym ansatz (GKBA) with Hartree-Fock propagators (HF-GKBA), and we present systematic improvements that use correlated propagators (correlated GKBA). △ Less

Submitted 19 September, 2018; v1 submitted 23 August, 2018; originally announced August 2018.

Journal ref: phys. stat. sol. (b) 257, issue 7, article 1800490 (2019)

Long-lasting field-free alignment of large molecules inside helium nanodroplets

Authors: Adam S. Chatterley, Constant Schouder, Lars Christiansen, Benjamin Shepperson, Mette H. Rasmussen, Henrik Stapelfeldt

Abstract: Molecules with their axes sharply confined in space, available through laser-induced alignment methods, are essential for many current experiments, including ultrafast molecular imaging. Most of these applications require both that the aligning laser field is turned-off, to avoid undesired perturbations, and that the molecules remain aligned sufficiently long that reactions and dynamics can be map… ▽ More Molecules with their axes sharply confined in space, available through laser-induced alignment methods, are essential for many current experiments, including ultrafast molecular imaging. Most of these applications require both that the aligning laser field is turned-off, to avoid undesired perturbations, and that the molecules remain aligned sufficiently long that reactions and dynamics can be mapped out. Presently, this is only possible for small, linear molecules and for times less than 1 picosecond. Here, we demonstrate strong, field-free alignment of large molecules inside helium nanodroplets, lasting tens of picoseconds. Molecular alignment in either one or three dimensions is created by a slowly switched-on laser pulse, made field-free through rapid pulse truncation, and retained thanks to the impeding effect of the helium environment on molecular rotation. We illustrate the opportunities that field-free aligned molecules open by measuring the alignment-dependent strong-field ionization yield of a thiophene oligomer. Our technique will enable molecular-frame experiments, including ultrafast excited state dynamics, on a variety of large molecules and complexes. △ Less

Submitted 19 July, 2018; originally announced July 2018.

Conditioned Lyapunov exponents for random dynamical systems

Authors: Maximilian Engel, Jeroen S. W. Lamb, Martin Rasmussen

Abstract: We introduce the notion of Lyapunov exponents for random dynamical systems, conditioned to trajectories that stay within a bounded domain for asymptotically long times. This is motivated by the desire to characterize local dynamical properties in the presence of unbounded noise (when almost all trajectories are unbounded). We illustrate its use in the analysis of local bifurcations in this context… ▽ More We introduce the notion of Lyapunov exponents for random dynamical systems, conditioned to trajectories that stay within a bounded domain for asymptotically long times. This is motivated by the desire to characterize local dynamical properties in the presence of unbounded noise (when almost all trajectories are unbounded). We illustrate its use in the analysis of local bifurcations in this context. The theory of conditioned Lyapunov exponents of stochastic differential equations builds on the stochastic analysis of quasi-stationary distributions for killed processes and associated quasi-ergodic distributions. We show that conditioned Lyapunov exponents describe the local stability behaviour of trajectories that remain within a bounded domain and - in particular - that negative conditioned Lyapunov exponents imply local synchronisation. Furthermore, a conditioned dichotomy spectrum is introduced and its main characteristics are established. △ Less

Submitted 18 May, 2018; originally announced May 2018.

MSC Class: 37A50; 37H10; 37H15; 60F99

Early-warning signals for bifurcations in random dynamical systems with bounded noise

Authors: Christian Kuehn, Giuseppe Malavolta, Martin Rasmussen

Abstract: We consider discrete-time one-dimensional random dynamical systems with bounded noise, which generate an associated set-valued dynamical system. We provide necessary and sufficient conditions for a discontinuous bifurcation of a minimal invariant set of the set-valued dynamical system in terms of the derivatives of the so-called extremal maps. We propose an algorithm for reconstructing the derivat… ▽ More We consider discrete-time one-dimensional random dynamical systems with bounded noise, which generate an associated set-valued dynamical system. We provide necessary and sufficient conditions for a discontinuous bifurcation of a minimal invariant set of the set-valued dynamical system in terms of the derivatives of the so-called extremal maps. We propose an algorithm for reconstructing the derivatives of the extremal maps from a time series that is generated by iterations of the original random dynamical system. We demonstrate that the derivative reconstructed for different parameters can be used as an early-warning signal to detect an upcoming bifurcation, and apply the algorithm to the bifurcation analysis of the stochastic return map of the Koper model, which is a three-dimensional multiple time scale ordinary differential equation used as prototypical model for the formation of mixed-mode oscillation patterns. We apply our algorithm to data generated by this map to detect an upcoming transition. △ Less

Submitted 6 April, 2018; v1 submitted 1 March, 2018; originally announced March 2018.

MSC Class: 37G35; 37H20; 37C70; 49K21; 70K70

Doublon formation by ions impacting a strongly correlated finite lattice system

Authors: Karsten Balzer, Maximilian Rodriguez Rasmussen, Niclas Schlünzen, Jan-Philip Joost, Michael Bonitz

Abstract: Strongly correlated systems of fermions have a number of exciting collective properties. Among them, the creation of a lattice that is occupied by doublons, i.e. two quantum particles with opposite spins, offers interesting electronic properties. In the past a variety of methods have been proposed to control doublon formation, both, spatially and temporally. Here, a novel mechanism is proposed and… ▽ More Strongly correlated systems of fermions have a number of exciting collective properties. Among them, the creation of a lattice that is occupied by doublons, i.e. two quantum particles with opposite spins, offers interesting electronic properties. In the past a variety of methods have been proposed to control doublon formation, both, spatially and temporally. Here, a novel mechanism is proposed and verified by exact diagonalization and nonequilibrium Green functions simulations---fermionic doublon creation by the impact of energetic ions. We report the formation of a nonequilibrium steady state with homogeneous doublon distribution. The effect should be observable in strongly correlated solids in contact with a high-pressure plasma and in fermionic atoms in optical lattices. △ Less

Submitted 16 January, 2018; originally announced January 2018.

Journal ref: Phys. Rev. Lett. 121, 267602 (2018)

Hopf bifurcation with additive noise

Authors: Thai Son Doan, Maximilian Engel, Jeroen S. W. Lamb, Martin Rasmussen

Abstract: We consider the dynamics of a two-dimensional ordinary differential equation exhibiting a Hopf bifurcation subject to additive white noise and identify three dynamical phases: (I) a random attractor with uniform synchronisation of trajectories, (II) a random attractor with non-uniform synchronisation of trajectories and (III) a random attractor without synchronisation of trajectories. The random a… ▽ More We consider the dynamics of a two-dimensional ordinary differential equation exhibiting a Hopf bifurcation subject to additive white noise and identify three dynamical phases: (I) a random attractor with uniform synchronisation of trajectories, (II) a random attractor with non-uniform synchronisation of trajectories and (III) a random attractor without synchronisation of trajectories. The random attractors in phases (I) and (II) are random equilibrium points with negative Lyapunov exponents while in phase (III) there is a so-called random strange attractor with positive Lyapunov exponent. We analyse the occurrence of the different dynamical phases as a function of the linear stability of the origin (deterministic Hopf bifurcation parameter) and shear (ampitude-phase coupling parameter). We show that small shear implies synchronisation and obtain that synchronisation cannot be uniform in the absence of linear stability at the origin or in the presence of sufficiently strong shear. We provide numerical results in support of a conjecture that irrespective of the linear stability of the origin, there is a critical strength of the shear at which the system dynamics loses synchronisation and enters phase (III). △ Less

Submitted 26 October, 2017; originally announced October 2017.

MSC Class: 37C75; 37D45; 37G35; 37H10; 37H15

HyperENTM: Evolving Scalable Neural Turing Machines through HyperNEAT

Authors: Jakob Merrild, Mikkel Angaju Rasmussen, Sebastian Risi

Abstract: Recent developments within memory-augmented neural networks have solved sequential problems requiring long-term memory, which are intractable for traditional neural networks. However, current approaches still struggle to scale to large memory sizes and sequence lengths. In this paper we show how access to memory can be encoded geometrically through a HyperNEAT-based Neural Turing Machine (HyperENT… ▽ More Recent developments within memory-augmented neural networks have solved sequential problems requiring long-term memory, which are intractable for traditional neural networks. However, current approaches still struggle to scale to large memory sizes and sequence lengths. In this paper we show how access to memory can be encoded geometrically through a HyperNEAT-based Neural Turing Machine (HyperENTM). We demonstrate that using the indirect HyperNEAT encoding allows for training on small memory vectors in a bit-vector copy task and then applying the knowledge gained from such training to speed up training on larger size memory vectors. Additionally, we demonstrate that in some instances, networks trained to copy bit-vectors of size 9 can be scaled to sizes of 1,000 without further training. While the task in this paper is simple, these results could open up the problems amendable to networks with external memories to problems with larger memory vectors and theoretically unbounded memory sizes. △ Less

Submitted 12 October, 2017; originally announced October 2017.

Classification of random circle homeomorphisms up to topological conjugacy

Authors: Thai Son Doan, Jeroen S. W. Lamb, Julian Newman, Martin Rasmussen

Abstract: We provide a classification of random orientation-preserving homeomorphisms of $\mathbb{S}^1$, up to topological conjugacy of the random dynamical systems generated by i.i.d. iterates of the random homeomorphism. This classification covers all random circle homeomorphisms for which the noise space is a connected Polish space and an additional extremely weak condition is satisfied. We provide a classification of random orientation-preserving homeomorphisms of $\mathbb{S}^1$, up to topological conjugacy of the random dynamical systems generated by i.i.d. iterates of the random homeomorphism. This classification covers all random circle homeomorphisms for which the noise space is a connected Polish space and an additional extremely weak condition is satisfied. △ Less

Submitted 17 July, 2017; originally announced July 2017.

MSC Class: 37H10

Combinatorial classification of quantum lens spaces

Authors: Peter Lunding Jensen, Frederik Ravn Klausen, Peter M. R. Rasmussen

Abstract: We answer the question of how large the dimension of a quantum lens space must be, compared to the primary parameter $r$, for the isomorphism class to depend on the secondary parameters. Since classification results in C*-algebra theory reduces this question to one concerning a certain kind of $SL$-equivalence of integer matrices of a special form, our approach is entirely combinatorial and based… ▽ More We answer the question of how large the dimension of a quantum lens space must be, compared to the primary parameter $r$, for the isomorphism class to depend on the secondary parameters. Since classification results in C*-algebra theory reduces this question to one concerning a certain kind of $SL$-equivalence of integer matrices of a special form, our approach is entirely combinatorial and based on the counting of certain paths in the graphs shown by Hong and Szymański to describe the quantum lens spaces. △ Less

Submitted 15 January, 2017; originally announced January 2017.

Comments: 27 pages, 2 figures

MSC Class: 46L35 (primary); 05C30 (secondary)

Journal ref: Pacific J. Math. 297 (2018) 339-365

Investigations of the effects of random sampling patterns on the stability of generalized sampling

Authors: Robert Dahl Jacobsen, Jesper Møller, Morten Nielsen, Morten Grud Rasmussen

Abstract: We investigate how the choice of spatial point process for generating random sampling patterns affects the numerical stability of non-uniform generalized sampling between Fourier bases and Daubechies scaling functions. Specifically, we consider binomial, Poisson and determinantal point processes and demonstrate that the more regular point patterns from the determinantal point process are superior. We investigate how the choice of spatial point process for generating random sampling patterns affects the numerical stability of non-uniform generalized sampling between Fourier bases and Daubechies scaling functions. Specifically, we consider binomial, Poisson and determinantal point processes and demonstrate that the more regular point patterns from the determinantal point process are superior. △ Less

Submitted 17 September, 2017; v1 submitted 15 July, 2016; originally announced July 2016.

Generalized Sampling in Julia

Authors: Robert Dahl Jacobsen, Morten Nielsen, Morten Grud Rasmussen

Abstract: Generalized sampling is a numerically stable framework for obtaining reconstructions of signals in different bases and frames from their samples. In this paper, we will introduce a carefully documented toolbox for performing generalized sampling in Julia. Julia is a new language for technical computing with focus on performance, which is ideally suited to handle the large size problems often encou… ▽ More Generalized sampling is a numerically stable framework for obtaining reconstructions of signals in different bases and frames from their samples. In this paper, we will introduce a carefully documented toolbox for performing generalized sampling in Julia. Julia is a new language for technical computing with focus on performance, which is ideally suited to handle the large size problems often encountered in generalized sampling. The toolbox provides specialized solutions for the setup of Fourier bases and wavelets. The performance of the toolbox is compared to existing implementations of generalized sampling in MATLAB. △ Less

Submitted 23 November, 2016; v1 submitted 14 July, 2016; originally announced July 2016.