-
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.
-
Uniqueness up to Inner Automorphism of Regular Exact Borel Subalgebras
Authors:
Anna Rodriguez Rasmussen
Abstract:
Külshammer, König and Ovsienko proved that for any quasi-hereditary algebra $(A,\leq_A)$ there exists a Morita equivalent quasi-hereditary algebra $(R, \leq_R)$ containing a basic exact Borel subalgebra $B$. The obtained Borel subalgebra is in fact a regular exact Borel subalgebra. Later, Conde showed that given a quasi-hereditary algebra $(R,\leq_R)$ with a basic regular exact Borel subalgebra…
▽ More
Külshammer, König and Ovsienko proved that for any quasi-hereditary algebra $(A,\leq_A)$ there exists a Morita equivalent quasi-hereditary algebra $(R, \leq_R)$ containing a basic exact Borel subalgebra $B$. The obtained Borel subalgebra is in fact a regular exact Borel subalgebra. Later, Conde showed that given a quasi-hereditary algebra $(R,\leq_R)$ with a basic regular exact Borel subalgebra $B$ and a Morita equivalent quasi-hereditary algebra $(R',\leq_{R'})$ with a basic regular exact Borel subalgebra $B'$, the algebras $R$ and $R'$ are isomorphic, and Külshammer and Miemietz showed that there is even an isomorphism $\varphi:R\rightarrow R'$ such that $\varphi(B)=B'$.
In this article, we show that if $R=R'$, then $\varphi$ can be chosen to be an inner automorphism. Moreover, instead of just proving this for regular exact Borel subalgebras of quasi-hereditary algebras, we generalize this to an appropriate class of subalgebras of arbitrary finite-dimensional algebras. As an application, we show that if $(A, \leq_A)$ is a finite-dimensional algebra and $G$ is a finite group acting on $A$ via automorphisms, then under some natural compatibility conditions, there is a Morita equivalent quasi-hereditary algebra $(R, \leq_R)$ with a basic regular exact Borel subalgebra $B$ such that $g(B)=B$ for every $g\in G$.
△ Less
Submitted 26 March, 2024; v1 submitted 22 March, 2024;
originally announced March 2024.
-
Irrational rotations and 2-filling rays
Authors:
Lvzhou Chen,
Alexander J. Rasmussen
Abstract:
We study a skew product transformation associated to an irrational rotation of the circle [0,1]/~. This skew product keeps track of the number of times an orbit of the rotation lands in the two complementary intervals of {0,1/2} in the circle. We show that under certain conditions on the continued fraction expansion of the irrational number defining the rotation, the skew product transformation ha…
▽ More
We study a skew product transformation associated to an irrational rotation of the circle [0,1]/~. This skew product keeps track of the number of times an orbit of the rotation lands in the two complementary intervals of {0,1/2} in the circle. We show that under certain conditions on the continued fraction expansion of the irrational number defining the rotation, the skew product transformation has certain dense orbits. This is in spite of the presence of numerous non-dense orbits. We use this to construct laminations on infinite type surfaces with exotic properties. In particular, we show that for every infinite type surface with an isolated planar end, there is an infinite clique of 2-filling rays based at that end. These 2-filling rays are relevant to Bavard--Walker's loop graphs.
△ Less
Submitted 24 May, 2024; v1 submitted 29 February, 2024;
originally announced March 2024.
-
Mitigation of the Brighter-Fatter Effect in the LSST Camera
Authors:
Alex Broughton,
Yousuke Utsumi,
Andrés Plazas Malagón,
Christopher Waters,
Craig Lage,
Adam Snyder,
Andrew Rasmussen,
Stuart Marshall,
Jim Chiang,
Simona Murgia,
Aaron Roodman
Abstract:
Thick, fully depleted charge-coupled devices (CCDs) are known to exhibit non-linear behavior at high signal levels due to the dynamic behavior of charges collecting in the potential wells of pixels, called the brighter-fatter effect (BFE). This particularly impacts bright calibration stars, which appear larger than their intrinsic shape, creating a flux-dependent point-spread function (PSF) that i…
▽ More
Thick, fully depleted charge-coupled devices (CCDs) are known to exhibit non-linear behavior at high signal levels due to the dynamic behavior of charges collecting in the potential wells of pixels, called the brighter-fatter effect (BFE). This particularly impacts bright calibration stars, which appear larger than their intrinsic shape, creating a flux-dependent point-spread function (PSF) that if left unmitigated, could make up a large fraction of the error budget in Stage IV weak-lensing (WL) surveys such as the Legacy Survey of Space and Time (LSST). In this paper, we analyze image measurements of flat fields and artificial stars taken at different illumination levels with the LSST Camera (LSSTCam) at SLAC National Accelerator Laboratory in order to quantify this effect in the LSST Camera before and after a previously introduced correction technique. We observe that the BFE evolves anisotropically as a function of flux due to higher-order BFEs, which violates the fundamental assumption of this correction method. We then introduce a new sampling method based on a physically motivated model to account these higher-order terms in the correction, and then we test the modified correction on both datasets. We find that the new method corrects the effect in flat fields better than it corrects the effect in artificial stars which we conclude is the result of a unmodeled curl component of the deflection field by the correction. We use these results to define a new metric for the full-well capacity of our sensors and advise image processing strategies to further limit the impact of the effect on LSST WL science pathways.
△ Less
Submitted 5 December, 2023;
originally announced December 2023.
-
The Bayesian approach to inverse Robin problems
Authors:
Aksel Kaastrup Rasmussen,
Fanny Seizilles,
Mark Girolami,
Ieva Kazlauskaite
Abstract:
In this paper we investigate the Bayesian approach to inverse Robin problems. These are problems for certain elliptic boundary value problems of determining a Robin coefficient on a hidden part of the boundary from Cauchy data on the observable part. Such a nonlinear inverse problem arises naturally in the initialisation of large-scale ice sheet models that are crucial in climate and sea-level pre…
▽ More
In this paper we investigate the Bayesian approach to inverse Robin problems. These are problems for certain elliptic boundary value problems of determining a Robin coefficient on a hidden part of the boundary from Cauchy data on the observable part. Such a nonlinear inverse problem arises naturally in the initialisation of large-scale ice sheet models that are crucial in climate and sea-level predictions. We motivate the Bayesian approach for a prototypical Robin inverse problem by showing that the posterior mean converges in probability to the data-generating ground truth as the number of observations increases. Related to the stability theory for inverse Robin problems, we establish a logarithmic convergence rate for Sobolev-regular Robin coefficients, whereas for analytic coefficients we can attain an algebraic rate. The use of rescaled analytic Gaussian priors in posterior consistency for nonlinear inverse problems is new and may be of separate interest in other inverse problems. Our numerical results illustrate the convergence property in two observation settings.
△ Less
Submitted 29 November, 2023;
originally announced November 2023.
-
A Bayesian approach for consistent reconstruction of inclusions
Authors:
Babak Maboudi Afkham,
Kim Knudsen,
Aksel Kaastrup Rasmussen,
Tanja Tarvainen
Abstract:
This paper considers a Bayesian approach for inclusion detection in nonlinear inverse problems using two known and popular push-forward prior distributions: the star-shaped and level set prior distributions. We analyze the convergence of the corresponding posterior distributions in a small measurement noise limit. The methodology is general; it works for priors arising from any Hölder continuous t…
▽ More
This paper considers a Bayesian approach for inclusion detection in nonlinear inverse problems using two known and popular push-forward prior distributions: the star-shaped and level set prior distributions. We analyze the convergence of the corresponding posterior distributions in a small measurement noise limit. The methodology is general; it works for priors arising from any Hölder continuous transformation of Gaussian random fields and is applicable to a range of inverse problems. The level set and star-shaped prior distributions are examples of push-forward priors under Hölder continuous transformations that take advantage of the structure of inclusion detection problems. We show that the corresponding posterior mean converges to the ground truth in a proper probabilistic sense. Numerical tests on a two-dimensional quantitative photoacoustic tomography problem showcase the approach. The results highlight the convergence properties of the posterior distributions and the ability of the methodology to detect inclusions with sufficiently regular boundaries.
△ Less
Submitted 25 August, 2023;
originally announced August 2023.
-
Photometry, Centroid and Point-Spread Function Measurements in the LSST Camera Focal Plane Using Artificial Stars
Authors:
Johnny H. Esteves,
Yousuke Utsumi,
Adam Snyder,
Theo Schutt,
Alex Broughton,
Bahrudin Trbalic,
Sidney Mau,
Andrew Rasmussen,
Andrés A. Plazas Malagón,
Andrew Bradshaw,
Stuart Marshall,
Seth Digel,
James Chiang,
Marcelle Soares-Santos,
Aaron Roodman
Abstract:
The Vera C. Rubin Observatory's LSST Camera (LSSTCam) pixel response has been characterized using laboratory measurements with a grid of artificial stars. We quantify the contributions to photometry, centroid, point-spread function size, and shape measurement errors due to small anomalies in the LSSTCam CCDs. The main sources of those anomalies are quantum efficiency variations and pixel area vari…
▽ More
The Vera C. Rubin Observatory's LSST Camera (LSSTCam) pixel response has been characterized using laboratory measurements with a grid of artificial stars. We quantify the contributions to photometry, centroid, point-spread function size, and shape measurement errors due to small anomalies in the LSSTCam CCDs. The main sources of those anomalies are quantum efficiency variations and pixel area variations induced by the amplifier segmentation boundaries and "tree-rings" - circular variations in silicon doping concentration. This laboratory study using artificial stars projected on the sensors shows overall small effects. The residual effects on point-spread function (PSF) size and shape are below $0.1\%$, meeting the ten-year LSST survey science requirements. However, the CCD mid-line presents distortions that can have a moderate impact on PSF measurements. This feature can be avoided by masking the affected regions. Effects of tree-rings are observed on centroids and PSFs of the artificial stars and the nature of the effect is confirmed by a study of the flat-field response. Nevertheless, further studies of the full-focal plane with stellar data should more completely probe variations and might reveal new features, e.g. wavelength-dependent effects. The results of this study can be used as a guide for the on-sky operation of LSSTCam.
△ Less
Submitted 3 November, 2023; v1 submitted 1 August, 2023;
originally announced August 2023.
-
Largest hyperbolic actions of 3--manifold groups
Authors:
Carolyn Abbott,
Hoang Thanh Nguyen,
Alexander J. Rasmussen
Abstract:
The set of equivalence classes of cobounded actions of a group G on different hyperbolic metric spaces carries a natural partial order. Following Abbott--Balasubramanya--Osin, the group G is H--accessible if the resulting poset has a largest element. In this paper, we prove that every non-geometric 3--manifold has a finite cover with H--inaccessible fundamental group and give conditions under whic…
▽ More
The set of equivalence classes of cobounded actions of a group G on different hyperbolic metric spaces carries a natural partial order. Following Abbott--Balasubramanya--Osin, the group G is H--accessible if the resulting poset has a largest element. In this paper, we prove that every non-geometric 3--manifold has a finite cover with H--inaccessible fundamental group and give conditions under which the fundamental group of the original manifold is H--inaccessible. We also prove that every Croke--Kleiner admissible group (a class of graphs of groups that generalizes fundamental groups of 3--dimensional graph manifolds) has a finite index subgroup that is H--inaccessible.
△ Less
Submitted 12 May, 2023;
originally announced May 2023.
-
Quasi-Hereditary Skew Group Algebras
Authors:
Anna Rodriguez Rasmussen
Abstract:
Given an algebra and a finite group acting on it via automorphisms, a natural object of study is the associated skew group algebra. In this article, we study the relationship between quasi-hereditary structures on the original algebra and on the corresponding skew group algebra. Assuming a natural compatibility condition on the partial order, we show that the skew group algebra is quasi-hereditary…
▽ More
Given an algebra and a finite group acting on it via automorphisms, a natural object of study is the associated skew group algebra. In this article, we study the relationship between quasi-hereditary structures on the original algebra and on the corresponding skew group algebra. Assuming a natural compatibility condition on the partial order, we show that the skew group algebra is quasi-hereditary if and only if the original algebra is. Moreover, we show that in this setting an exact Borel subalgebra of the original algebra which is invariant as a set under the group action gives rise to an exact Borel subalgebra of the skew group algebra, and that under this construction, properties such as normality and regularity of the exact Borel subalgebra are preserved.
△ Less
Submitted 11 May, 2023;
originally announced May 2023.
-
Emergent properties of van der Waals bilayers revealed by computational stacking
Authors:
Sahar Pakdel,
Asbjørn Rasmussen,
Alireza Taghizadeh,
Mads Kruse,
Thomas Olsen,
Kristian S. Thygesen
Abstract:
Stacking of two-dimensional (2D) materials has emerged as a facile strategy for realising exotic quantum states of matter and engineering electronic properties. Yet, developments beyond the proof-of-principle level are impeded by the vast size of the configuration space defined by layer combinations and stacking orders. Here we employ a density functional theory (DFT) workflow to calculate interla…
▽ More
Stacking of two-dimensional (2D) materials has emerged as a facile strategy for realising exotic quantum states of matter and engineering electronic properties. Yet, developments beyond the proof-of-principle level are impeded by the vast size of the configuration space defined by layer combinations and stacking orders. Here we employ a density functional theory (DFT) workflow to calculate interlayer binding energies of 8451 homobilayers created by stacking 1052 different monolayers in various configurations. Analysis of the stacking orders in 247 experimentally known van der Waals crystals is used to validate the workflow and determine the criteria for realizable bilayers. For the 2586 most stable bilayer systems, we calculate a range of electronic, magnetic, and vibrational properties, and explore general trends and anomalies. We identify an abundance of bistable bilayers with stacking order-dependent magnetic or electrical polarisation states making them candidates for slidetronics applications.
△ Less
Submitted 3 April, 2023;
originally announced April 2023.
-
Transverse measures to infinite type laminations
Authors:
Mladen Bestvina,
Alexander J. Rasmussen
Abstract:
We study the cone of transverse measures to a fixed geodesic lamination on an infinite type hyperbolic surface. Under simple hypotheses on the metric, we give an explicit description of this cone as an inverse limit of finite-dimensional cones. We study the problem of when the cone of transverse measures admits a base and show that such a base exists for many laminations. Moreover, the base is a (…
▽ More
We study the cone of transverse measures to a fixed geodesic lamination on an infinite type hyperbolic surface. Under simple hypotheses on the metric, we give an explicit description of this cone as an inverse limit of finite-dimensional cones. We study the problem of when the cone of transverse measures admits a base and show that such a base exists for many laminations. Moreover, the base is a (typically infinite-dimensional) simplex (called a Choquet simplex) and can be described explicitly as an inverse limit of finite-dimensional simplices. We show that on any fixed infinite type hyperbolic surface, every Choquet simplex arises as a base for some lamination. We use our inverse limit description and a new construction of geodesic laminations to give other explicit examples of cones with exotic properties.
△ Less
Submitted 16 August, 2023; v1 submitted 31 August, 2022;
originally announced September 2022.
-
Valuations, completions, and hyperbolic actions of metabelian groups
Authors:
Carolyn R. Abbott,
Sahana Balasubramanya,
Sam Payne,
Alexander J. Rasmussen
Abstract:
Actions on hyperbolic metric spaces are an important tool for studying groups, and so it is natural, but difficult, to attempt to classify all such actions of a fixed group. In this paper, we build strong connections between hyperbolic geometry and commutative algebra in order to classify the cobounded hyperbolic actions of numerous metabelian groups up to a coarse equivalence. In particular, we t…
▽ More
Actions on hyperbolic metric spaces are an important tool for studying groups, and so it is natural, but difficult, to attempt to classify all such actions of a fixed group. In this paper, we build strong connections between hyperbolic geometry and commutative algebra in order to classify the cobounded hyperbolic actions of numerous metabelian groups up to a coarse equivalence. In particular, we turn this classification problem into the problems of classifying ideals in the completions of certain rings and calculating invariant subspaces of matrices. We use this framework to classify the cobounded hyperbolic actions of many abelian-by-cyclic groups associated to expanding integer matrices. Each such action is equivalent to an action on a tree or on a Heintze group (a classically studied class of negatively curved Lie groups). Our investigations incorporate number systems, factorization in formal power series rings, completions, and valuations.
△ Less
Submitted 26 July, 2022;
originally announced July 2022.
-
(In-)Approximability Results for Interval, Resource Restricted, and Low Rank Scheduling
Authors:
Marten Maack,
Simon Pukrop,
Anna Rodriguez Rasmussen
Abstract:
We consider variants of the restricted assignment problem where a set of jobs has to be assigned to a set of machines, for each job a size and a set of eligible machines is given, and the jobs may only be assigned to eligible machines with the goal of makespan minimization. For the variant with interval restrictions, where the machines can be arranged on a path such that each job is eligible on a…
▽ More
We consider variants of the restricted assignment problem where a set of jobs has to be assigned to a set of machines, for each job a size and a set of eligible machines is given, and the jobs may only be assigned to eligible machines with the goal of makespan minimization. For the variant with interval restrictions, where the machines can be arranged on a path such that each job is eligible on a subpath, we present the first better than $2$-approximation and an improved inapproximability result. In particular, we give a $(2-\frac{1}{24})$-approximation and show that no better than $9/8$-approximation is possible, unless P=NP. Furthermore, we consider restricted assignment with $R$ resource restrictions and rank $D$ unrelated scheduling. In the former problem, a machine may process a job if it can meet its resource requirements regarding $R$ (renewable) resources. In the latter, the size of a job is dependent on the machine it is assigned to and the corresponding processing time matrix has rank at most $D$. The problem with interval restrictions includes the 1 resource variant, is encompassed by the 2 resource variant, and regarding approximation the $R$ resource variant is essentially a special case of the rank $R+1$ problem. We show that no better than $3/2$, $8/7$, and $3/2$-approximation is possible (unless P=NP) for the 3 resource, 2 resource, and rank 3 variant, respectively. Both the approximation result for the interval case and the inapproximability result for the rank 3 variant are solutions to open challenges stated in previous works. Lastly, we also consider the reverse objective, that is, maximizing the minimal load any machine receives, and achieve similar results.
△ Less
Submitted 11 March, 2022;
originally announced March 2022.
-
Fixation Maximization in the Positional Moran Process
Authors:
Joachim Brendborg,
Panagiotis Karras,
Andreas Pavlogiannis,
Asger Ullersted Rasmussen,
Josef Tkadlec
Abstract:
The Moran process is a classic stochastic process that models invasion dynamics on graphs. A single "mutant" (e.g., a new opinion, strain, social trait etc.) invades a population of residents spread over the nodes of a graph. The mutant fitness advantage $δ\geq 0$ determines how aggressively mutants propagate to their neighbors. The quantity of interest is the fixation probability, i.e., the proba…
▽ More
The Moran process is a classic stochastic process that models invasion dynamics on graphs. A single "mutant" (e.g., a new opinion, strain, social trait etc.) invades a population of residents spread over the nodes of a graph. The mutant fitness advantage $δ\geq 0$ determines how aggressively mutants propagate to their neighbors. The quantity of interest is the fixation probability, i.e., the probability that the initial mutant eventually takes over the whole population. However, in realistic settings, the invading mutant has an advantage only in certain locations. E.g., a bacterial mutation allowing for lactose metabolism only confers an advantage on places where dairy products are present. In this paper we introduce the positional Moran process, a natural generalization in which the mutant fitness advantage is only realized on specific nodes called active nodes. The associated optimization problem is fixation maximization: given a budget $k$, choose a set of $k$ active nodes that maximize the fixation probability of the invading mutant. We show that the problem is NP-hard, while the optimization function is not submodular, thus indicating strong computational hardness. Then we focus on two natural limits. In the limit of $δ\to\infty$ (strong selection), although the problem remains NP-hard, the optimization function becomes submodular and thus admits a constant-factor approximation using a simple greedy algorithm. In the limit of $δ\to 0$ (weak selection), we show that in $O(m^ω)$ time we can obtain a tight approximation, where $m$ is the number of edges and $ω$ is the matrix-multiplication exponent. Finally, we present an experimental evaluation of the new algorithms together with some proposed heuristics.
△ Less
Submitted 25 April, 2022; v1 submitted 6 January, 2022;
originally announced January 2022.
-
Higher rank confining subsets and hyperbolic actions of solvable groups
Authors:
Carolyn R. Abbott,
Sahana Balasubramanya,
Alexander J. Rasmussen
Abstract:
Recent papers of the authors have completely described the hyperbolic actions of several families of classically studied solvable groups. A key tool for these investigations is the machinery of confining subsets of Caprace, Cornulier, Monod, and Tessera, which applies, in particular, to solvable groups with virtually cyclic abelianizations. In this paper, we extend this machinery and give a corres…
▽ More
Recent papers of the authors have completely described the hyperbolic actions of several families of classically studied solvable groups. A key tool for these investigations is the machinery of confining subsets of Caprace, Cornulier, Monod, and Tessera, which applies, in particular, to solvable groups with virtually cyclic abelianizations. In this paper, we extend this machinery and give a correspondence between the hyperbolic actions of certain solvable groups with higher rank abelianizations and confining subsets of these more general groups. We then apply this extension to give a complete description of the hyperbolic actions of generalized solvable Baumslag-Solitar groups and to reprove a result of Sgobbi-Wong computing their Bieri-Neumann-Strebel invariants.
△ Less
Submitted 17 May, 2023; v1 submitted 18 August, 2021;
originally announced August 2021.
-
Stable commutator length on big mapping class groups
Authors:
Elizabeth Field,
Priyam Patel,
Alexander J. Rasmussen
Abstract:
We study stable commutator length on mapping class groups of certain infinite-type surfaces. In particular, we show that stable commutator length defines a continuous function on the commutator subgroups of such infinite-type mapping class groups. We furthermore show that the commutator subgroups are open and closed subgroups and that the abelianizations are finitely generated in many cases. Our r…
▽ More
We study stable commutator length on mapping class groups of certain infinite-type surfaces. In particular, we show that stable commutator length defines a continuous function on the commutator subgroups of such infinite-type mapping class groups. We furthermore show that the commutator subgroups are open and closed subgroups and that the abelianizations are finitely generated in many cases. Our results apply to many popular infinite-type surfaces with locally coarsely bounded mapping class groups.
△ Less
Submitted 9 June, 2022; v1 submitted 4 August, 2021;
originally announced August 2021.
-
Direct regularized reconstruction for the three-dimensional Calderón problem
Authors:
Kim Knudsen,
Aksel K. Rasmussen
Abstract:
Electrical Impedance Tomography gives rise to the severely ill-posed Calderón problem of determining the electrical conductivity distribution in a bounded domain from knowledge of the associated Dirichlet-to-Neumann map for the governing equation. The uniqueness and stability questions for the three-dimensional problem were largely answered in the affirmative in the 1980's using complex geometrica…
▽ More
Electrical Impedance Tomography gives rise to the severely ill-posed Calderón problem of determining the electrical conductivity distribution in a bounded domain from knowledge of the associated Dirichlet-to-Neumann map for the governing equation. The uniqueness and stability questions for the three-dimensional problem were largely answered in the affirmative in the 1980's using complex geometrical optics solutions, and this led further to a direct reconstruction method relying on a non-physical scattering transform. In this paper, the reconstruction problem is taken one step further towards practical applications by considering data contaminated by noise. Indeed, a regularization strategy for the three-dimensional Calderón problem is presented based on a suitable and explicit truncation of the scattering transform. This gives a certified, stable and direct reconstruction method that is robust to small perturbations of the data. Numerical tests on simulated noisy data illustrate the feasibility and regularizing effect of the method, and suggest that the numerical implementation performs better than predicted by theory.
△ Less
Submitted 25 January, 2022; v1 submitted 15 June, 2021;
originally announced June 2021.
-
Atomic Simulation Recipes -- a Python framework and library for automated workflows
Authors:
Morten Gjerding,
Thorbjørn Skovhus,
Asbjørn Rasmussen,
Fabian Bertoldo,
Ask Hjorth Larsen,
Jens Jørgen Mortensen,
Kristian Sommer Thygesen
Abstract:
The Atomic Simulation Recipes (ASR) is an open source Python framework for working with atomistic materials simulations in an efficient and sustainable way that is ideally suited for high-throughput projects. Central to ASR is the concept of a Recipe: a high-level Python script that performs a well defined simulation task robustly and accurately while keeping track of the data provenance. The ASR…
▽ More
The Atomic Simulation Recipes (ASR) is an open source Python framework for working with atomistic materials simulations in an efficient and sustainable way that is ideally suited for high-throughput projects. Central to ASR is the concept of a Recipe: a high-level Python script that performs a well defined simulation task robustly and accurately while keeping track of the data provenance. The ASR leverages the functionality of the Atomic Simulation Environment (ASE) to interface with external simulation codes and attain a high abstraction level. We provide a library of Recipes for common simulation tasks employing density functional theory and many-body perturbation schemes. These Recipes utilize the GPAW electronic structure code, but may be adapted to other simulation codes with an ASE interface. Being independent objects with automatic data provenance control, Recipes can be freely combined through Python scripting giving maximal freedom for users to build advanced workflows. ASR also implements a command line interface that can be used to run Recipes and inspect results. The ASR Migration module helps users maintain their data while the Database and App modules makes it possible to create local databases and present them as customized web pages.
△ Less
Submitted 27 April, 2021;
originally announced April 2021.
-
Improvement of Noise-Robust Single-Channel Voice Activity Detection with Spatial Pre-processing
Authors:
Max Væhrens,
Andreas Jonas Fuglsig,
Anders Post Jacobsen,
Nicolai Almskou Rasmussen,
Victor Mølbach Nissen,
Joachim Roland Hejslet,
Zheng-Hua Tan
Abstract:
Voice activity detection (VAD) remains a challenge in noisy environments. With access to multiple microphones, prior studies have attempted to improve the noise robustness of VAD by creating multi-channel VAD (MVAD) methods. However, MVAD is relatively new compared to single-channel VAD (SVAD), which has been thoroughly developed in the past. It might therefore be advantageous to improve SVAD meth…
▽ More
Voice activity detection (VAD) remains a challenge in noisy environments. With access to multiple microphones, prior studies have attempted to improve the noise robustness of VAD by creating multi-channel VAD (MVAD) methods. However, MVAD is relatively new compared to single-channel VAD (SVAD), which has been thoroughly developed in the past. It might therefore be advantageous to improve SVAD methods with pre-processing to obtain superior VAD, which is under-explored. This paper improves SVAD through two pre-processing methods, a beamformer and a spatial target speaker detector. The spatial detector sets signal frames to zero when no potential speaker is present within a target direction. The detector may be implemented as a filter, meaning the input signal for the SVAD is filtered according to the detector's output; or it may be implemented as a spatial VAD to be combined with the SVAD output. The evaluation is made on a noisy reverberant speech database, with clean speech from the Aurora 2 database and with white and babble noise. The results show that SVAD algorithms are significantly improved by the presented pre-processing methods, especially the spatial detector, across all signal-to-noise ratios. The SVAD algorithms with pre-processing significantly outperform a baseline MVAD in challenging noise conditions.
△ Less
Submitted 12 April, 2021;
originally announced April 2021.
-
Recent Progress of the Computational 2D Materials Database (C2DB)
Authors:
M. N. Gjerding,
A. Taghizadeh,
A. Rasmussen,
S. Ali,
F. Bertoldo,
T. Deilmann,
U. P. Holguin,
N. R. Knøsgaard,
M. Kruse,
S. Manti,
T. G. Pedersen,
T. Skovhus,
M. K. Svendsen,
J. J. Mortensen,
T. Olsen,
K. S. Thygesen
Abstract:
The C2DB is a highly curated open database organizing a wealth of computed properties for more than 4000 atomically thin two-dimensional (2D) materials. Here we report on new materials and properties that were added to the database since its first release in 2018. The set of new materials comprise several hundred monolayers exfoliated from experimentally known layered bulk materials, (homo)bilayer…
▽ More
The C2DB is a highly curated open database organizing a wealth of computed properties for more than 4000 atomically thin two-dimensional (2D) materials. Here we report on new materials and properties that were added to the database since its first release in 2018. The set of new materials comprise several hundred monolayers exfoliated from experimentally known layered bulk materials, (homo)bilayers in various stacking configurations, native point defects in semiconducting monolayers, and chalcogen/halogen Janus monolayers. The new properties include exfoliation energies, Bader charges, spontaneous polarisations, Born charges, infrared polarisabilities, piezoelectric tensors, band topology invariants, exchange couplings, Raman- and second harmonic generation spectra. We also describe refinements of the employed material classification schemes, upgrades of the computational methodologies used for property evaluations, as well as significant enhancements of the data documentation and provenance. Finally, we explore the performance of Gaussian process-based regression for efficient prediction of mechanical and electronic materials properties. The combination of open access, detailed documentation, and extremely rich materials property data sets make the C2DB a unique resource that will advance the science of atomically thin materials.
△ Less
Submitted 5 February, 2021;
originally announced February 2021.
-
Laboratory Measurements of Instrumental Signatures of the LSST Camera Focal Plane
Authors:
Adam Snyder,
Aurelien Barrau,
Andrew Bradshaw,
Boyd Bowdish,
James Chiang,
Celine Combet,
Seth Digel,
Richard Dubois,
Ludovic Eraud,
Claire Juramy,
Craig Lage,
Travis Lange,
Myriam Migliore,
Andrei Nomerotski,
Paul O'Connor,
HyeYun Park,
Andrew Rasmussen,
Kevin Reil,
Aaron Roodman,
Adrian Shestakov,
Yousuke Utsumi,
Duncan Wood
Abstract:
Electro-optical testing and characterization of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) Camera focal plane, consisting of 205 charge-coupled devices (CCDs) arranged into 21 stand-alone Raft Tower Modules (RTMs) and 4 Corner Raft Tower Modules (CRTMs), is currently being performed at the SLAC National Accelerator Laboratory. Testing of the camera sensors is performed us…
▽ More
Electro-optical testing and characterization of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) Camera focal plane, consisting of 205 charge-coupled devices (CCDs) arranged into 21 stand-alone Raft Tower Modules (RTMs) and 4 Corner Raft Tower Modules (CRTMs), is currently being performed at the SLAC National Accelerator Laboratory. Testing of the camera sensors is performed using a set of custom-built optical projectors, designed to illuminate the full focal plane or specific regions of the focal plane with a series of light illumination patterns: the crosstalk projector, the flat illuminator projector, and the spot grid projector. In addition to measurements of crosstalk, linearity and full well, the ability to project realistically-sized sources, using the spot grid projector, makes possible unique measurements of instrumental signatures such as deferred charge distortions, astrometric shifts due to sensor effects, and the brighter-fatter effect, prior to camera first light. Here we present the optical projector designs and usage, the electro-optical measurements and how these results have been used in testing and improving the LSST Camera instrumental signature removal algorithms.
△ Less
Submitted 4 January, 2021;
originally announced January 2021.
-
Laminations and 2-filling rays on infinite type surfaces
Authors:
Lvzhou Chen,
Alexander J. Rasmussen
Abstract:
The loop graph of an infinite type surface is an infinite diameter hyperbolic graph first studied in detail by Juliette Bavard. An important open problem in the study of infinite type surfaces is to describe the boundary of the loop graph as a space of geodesic laminations. We approach this problem by constructing the first examples of 2-filling rays on infinite type surfaces. Such rays accumulate…
▽ More
The loop graph of an infinite type surface is an infinite diameter hyperbolic graph first studied in detail by Juliette Bavard. An important open problem in the study of infinite type surfaces is to describe the boundary of the loop graph as a space of geodesic laminations. We approach this problem by constructing the first examples of 2-filling rays on infinite type surfaces. Such rays accumulate onto geodesic laminations which are in some sense filling, but without strong enough properties to correspond to points in the boundary of the loop graph. We give multiple constructions using both a hands-on combinatorial approach and an approach using train tracks and automorphisms of flat surfaces. In addition, our approaches are sufficiently robust to describe all 2-filling rays with certain other basic properties as well as to produce uncountably many distinct mapping class group orbits.
△ Less
Submitted 1 June, 2022; v1 submitted 12 October, 2020;
originally announced October 2020.
-
The LSST DESC DC2 Simulated Sky Survey
Authors:
LSST Dark Energy Science Collaboration,
Bela Abolfathi,
David Alonso,
Robert Armstrong,
Éric Aubourg,
Humna Awan,
Yadu N. Babuji,
Franz Erik Bauer,
Rachel Bean,
George Beckett,
Rahul Biswas,
Joanne R. Bogart,
Dominique Boutigny,
Kyle Chard,
James Chiang,
Chuck F. Claver,
Johann Cohen-Tanugi,
Céline Combet,
Andrew J. Connolly,
Scott F. Daniel,
Seth W. Digel,
Alex Drlica-Wagner,
Richard Dubois,
Emmanuel Gangler,
Eric Gawiser
, et al. (55 additional authors not shown)
Abstract:
We describe the simulated sky survey underlying the second data challenge (DC2) carried out in preparation for analysis of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) by the LSST Dark Energy Science Collaboration (LSST DESC). Significant connections across multiple science domains will be a hallmark of LSST; the DC2 program represents a unique modeling effort that stresses…
▽ More
We describe the simulated sky survey underlying the second data challenge (DC2) carried out in preparation for analysis of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) by the LSST Dark Energy Science Collaboration (LSST DESC). Significant connections across multiple science domains will be a hallmark of LSST; the DC2 program represents a unique modeling effort that stresses this interconnectivity in a way that has not been attempted before. This effort encompasses a full end-to-end approach: starting from a large N-body simulation, through setting up LSST-like observations including realistic cadences, through image simulations, and finally processing with Rubin's LSST Science Pipelines. This last step ensures that we generate data products resembling those to be delivered by the Rubin Observatory as closely as is currently possible. The simulated DC2 sky survey covers six optical bands in a wide-fast-deep (WFD) area of approximately 300 deg^2 as well as a deep drilling field (DDF) of approximately 1 deg^2. We simulate 5 years of the planned 10-year survey. The DC2 sky survey has multiple purposes. First, the LSST DESC working groups can use the dataset to develop a range of DESC analysis pipelines to prepare for the advent of actual data. Second, it serves as a realistic testbed for the image processing software under development for LSST by the Rubin Observatory. In particular, simulated data provide a controlled way to investigate certain image-level systematic effects. Finally, the DC2 sky survey enables the exploration of new scientific ideas in both static and time-domain cosmology.
△ Less
Submitted 26 January, 2021; v1 submitted 12 October, 2020;
originally announced October 2020.
-
Towards fully automatized GW band structure calculations: What we can learn from 60.000 self-energy evaluations
Authors:
Asbjørn Rasmussen,
Thorsten Deilmann,
Kristian S. Thygesen
Abstract:
We analyze a data set comprising 370 GW band structures composed of 61716 quasiparticle (QP) energies of two-dimensional (2D) materials spanning 14 crystal structures and 52 elements. The data results from PAW plane wave based one-shot G$_0$W$_0$@PBE calculations with full frequency integration. We investigate the distribution of key quantities like the QP self-energy corrections and renormalizati…
▽ More
We analyze a data set comprising 370 GW band structures composed of 61716 quasiparticle (QP) energies of two-dimensional (2D) materials spanning 14 crystal structures and 52 elements. The data results from PAW plane wave based one-shot G$_0$W$_0$@PBE calculations with full frequency integration. We investigate the distribution of key quantities like the QP self-energy corrections and renormalization factor $Z$ and explore their dependence on chemical composition and magnetic state. The linear QP approximation is identified as a significant error source and propose schemes for controlling and drastically reducing this error at low computational cost. We analyze the reliability of the $1/N_\text{PW}$ basis set extrapolation and find that is well-founded with narrow distributions of $r^2$ peaked very close to 1. Finally, we explore the validity of the scissors operator approximation concluding that it is generally not valid for reasonable error tolerances. Our work represents a step towards the development of automatized workflows for high-throughput G$_0$W$_0$ band structure calculations for solids.
△ Less
Submitted 4 January, 2021; v1 submitted 1 September, 2020;
originally announced September 2020.
-
Higher-form Gauge Symmetries in Multipole Topological Phases
Authors:
Oleg Dubinkin,
Alex Rasmussen,
Taylor L. Hughes
Abstract:
In this article we study field-theoretical aspects of multipolar topological insulators. Previous research has shown that such systems naturally couple to higher-rank tensor gauge fields that arise as a result of gauging dipole or subsystem $U(1)$ symmetries. Here we propose a complementary framework using electric higher-form symmetries. We utilize the fact that gauging 1-form electric symmetries…
▽ More
In this article we study field-theoretical aspects of multipolar topological insulators. Previous research has shown that such systems naturally couple to higher-rank tensor gauge fields that arise as a result of gauging dipole or subsystem $U(1)$ symmetries. Here we propose a complementary framework using electric higher-form symmetries. We utilize the fact that gauging 1-form electric symmetries results in a 2-form gauge field which couples naturally to extended line-like objects: Wilson lines. In our context the Wilson lines are electric flux lines associated to the electric polarization of the system. This allows us to define a generalized 2-form Peierls' substitution for dipoles that shows that the off-diagonal components of a rank-2 tensor gauge field $A_{ij}$ can arise as a lattice Peierls factor generated by the background antisymmetric 2-form gauge field. This framework has immediate applications: (i) it allows us to construct a manifestly topological quadrupolar response action given by a Dixmier-Douady invariant -- a generalization of a Chern number for 2-form gauge fields -- which makes plain the quantization of the quadrupole moment in the presence of certain crystal symmetries; (ii) it allows for a clearer interpretation of the rank-2 Berry phase calculation of the quadrupole moment; (iii) it allows for a proof of a generic Lieb-Schultz-Mattis theorem for dipole-conserving systems.
△ Less
Submitted 10 July, 2020;
originally announced July 2020.
-
Anisotropic properties of monolayer 2D materials: an overview from the C2DB database
Authors:
Luca Vannucci,
Urko Petralanda,
Asbjørn Rasmussen,
Thomas Olsen,
Kristian S. Thygesen
Abstract:
We analyze the occurrence of in-plane anisotropy in the electronic, magnetic, elastic and transport properties of more than one thousand 2D materials from the C2DB database. We identify hundreds of anisotropic materials and classify them according to their point group symmetry and degree of anisotropy. A statistical analysis reveals that a lower point group symmetry and a larger amount of differen…
▽ More
We analyze the occurrence of in-plane anisotropy in the electronic, magnetic, elastic and transport properties of more than one thousand 2D materials from the C2DB database. We identify hundreds of anisotropic materials and classify them according to their point group symmetry and degree of anisotropy. A statistical analysis reveals that a lower point group symmetry and a larger amount of different elements in the structure favour all types of anisotropies, which could be relevant for future materials design approaches. Besides, we identify novel compounds, predicted to be easily exfoliable from a parent bulk compound, with anisotropies that largely outscore those of already known 2D materials. Our findings provide a comprehensive reference for future studies of anisotropic response in atomically-thin crystals and point to new previously unexplored materials for the next generation of anisotropic 2D devices.
△ Less
Submitted 19 August, 2020; v1 submitted 8 July, 2020;
originally announced July 2020.
-
Sensor Distortion Effects in Photon Monte Carlo Simulations
Authors:
J. R. Peterson,
P. O'Connor,
A. Nomerotski,
E. Magnier,
J. G. Jernigan,
J. Cheng,
W. Cui,
E. Peng,
A. Rasmussen,
G. Sembroski
Abstract:
We present a detailed method to simulating sensor distortions using a photon and electron Monte Carlo method. We use three dimensional electrostatic simulations to parameterize the perturbed electric field profile for non-ideal sensor details. We follow the conversion of simulated photons, and the subsequent response of the converted electrons to the electric field pattern. These non-ideal sensor…
▽ More
We present a detailed method to simulating sensor distortions using a photon and electron Monte Carlo method. We use three dimensional electrostatic simulations to parameterize the perturbed electric field profile for non-ideal sensor details. We follow the conversion of simulated photons, and the subsequent response of the converted electrons to the electric field pattern. These non-ideal sensor details can be implemented efficiently in a Monte Carlo approach. We demonstrate that the non-ideal sensor distortions have a variety of observable consequence including the modification of the astrometric pattern, the distortion of the electron diffusion size and shape, and the distortion of flats. We show analytic validation of the diffusion physics, reproduce two kinds of edge distortion, and show qualitative validation of field-free regions, lithography errors, and fringing. We also demonstrate that there are two related effects of doping variation having different observable consequences. We show that field distortions from accumulated electrons lead to intensity-dependent point-spread-functions and the sub-linear variance in flats. The method is implemented in the Photon Simulator (PhoSim) and the code is publically available.
△ Less
Submitted 21 January, 2020; v1 submitted 9 January, 2020;
originally announced January 2020.
-
Efficient Reordered Nonlinear Gauss-Seidel Solvers With Higher Order For Black-Oil Models
Authors:
Øystein S. Klemetsdal,
Atgeirr F. Rasmussen,
Olav Møyner,
Knut-Andreas Lie
Abstract:
The fully implicit method is the most commonly used approach to solve black-oil problems in reservoir simulation. The method requires repeated linearization of large nonlinear systems and produces ill-condi\-tioned linear systems. We present a strategy to reduce computational time that relies on two key ideas: (\textit{i}) a sequential formulation that decouples flow and transport into separate su…
▽ More
The fully implicit method is the most commonly used approach to solve black-oil problems in reservoir simulation. The method requires repeated linearization of large nonlinear systems and produces ill-condi\-tioned linear systems. We present a strategy to reduce computational time that relies on two key ideas: (\textit{i}) a sequential formulation that decouples flow and transport into separate subproblems, and (\textit{ii}) a highly efficient Gauss--Seidel solver for the transport problems. This solver uses intercell fluxes to reorder the grid cells according to their upstream neighbors, and groups cells that are mutually dependent because of counter-current flow into local clusters. The cells and local clusters can then be solved in sequence, starting from the inflow and moving gradually downstream, since each new cell or local cluster will only depend on upstream neighbors that have already been computed. Altogether, this gives optimal localization and control of the nonlinear solution process.
This method has been successfully applied to real-field problems using the standard first-order finite volume discretization. Here, we extend the idea to first-order dG methods on fully unstructured grids. We also demonstrate proof of concept for the reordering idea by applying it to the full simulation model of the Norne oil field, using a prototype variant of the open-source OPM Flow simulator.
△ Less
Submitted 6 January, 2020;
originally announced January 2020.
-
Largest hyperbolic actions and quasi-parabolic actions in groups
Authors:
Carolyn R. Abbott,
Alexander J. Rasmussen
Abstract:
The set of equivalence classes of cobounded actions of a group on different hyperbolic metric spaces carries a natural partial order. The resulting poset thus gives rise to a notion of the "best" hyperbolic action of a group as the largest element of this poset, if such an element exists. We call such an action a largest hyperbolic action. While hyperbolic groups admit largest hyperbolic actions,…
▽ More
The set of equivalence classes of cobounded actions of a group on different hyperbolic metric spaces carries a natural partial order. The resulting poset thus gives rise to a notion of the "best" hyperbolic action of a group as the largest element of this poset, if such an element exists. We call such an action a largest hyperbolic action. While hyperbolic groups admit largest hyperbolic actions, we give evidence in this paper that this phenomenon is rare for non-hyperbolic groups. In particular, we prove that many families of groups of geometric origin do not have largest hyperbolic actions, including for instance many 3-manifold groups and most mapping class groups. Our proofs use the quasi-trees of metric spaces of Bestvina--Bromberg--Fujiwara, among other tools. In addition, we give a complete characterization of the poset of hyperbolic actions of Anosov mapping torus groups, and we show that mapping class groups of closed surfaces of genus at least two have hyperbolic actions which are comparable only to the trivial action.
△ Less
Submitted 12 March, 2022; v1 submitted 30 October, 2019;
originally announced October 2019.
-
Geometry of the graphs of nonseparating curves: covers and boundaries
Authors:
Alexander J. Rasmussen
Abstract:
We investigate the geometry of the graphs of nonseparating curves for surfaces of finite positive genus with potentially infinitely many punctures. This graph has infinite diameter and is known to be Gromov hyperbolic by work of the author. We study finite covers between such surfaces and show that lifts of nonseparating curves to the nonseparating curve graph of the cover span quasiconvex subgrap…
▽ More
We investigate the geometry of the graphs of nonseparating curves for surfaces of finite positive genus with potentially infinitely many punctures. This graph has infinite diameter and is known to be Gromov hyperbolic by work of the author. We study finite covers between such surfaces and show that lifts of nonseparating curves to the nonseparating curve graph of the cover span quasiconvex subgraphs which are infinite diameter and not coarsely equal to the nonseparating curve graph of the cover. In the finite type case, we also reprove a theorem of Hamenstädt identifying the Gromov boundary with the space of ending laminations on full genus subsurfaces. We introduce several tools based around the analysis of bicorn curves and laminations which may be of independent interest for studying the geometry of nonseparating curve graphs of infinite type surfaces and their boundaries.
△ Less
Submitted 6 August, 2020; v1 submitted 29 October, 2019;
originally announced October 2019.
-
The Open Porous Media Flow Reservoir Simulator
Authors:
Atgeirr Flø Rasmussen,
Tor Harald Sandve,
Kai Bao,
Andreas Lauser,
Joakim Hove,
Bård Skaflestad,
Robert Klöfkorn,
Markus Blatt,
Alf Birger Rustad,
Ove Sævareid,
Knut-Andreas Lie,
Andreas Thune
Abstract:
The Open Porous Media (OPM) initiative is a community effort that encourages open innovation and reproducible research for simulation of porous media processes. OPM coordinates collaborative software development, maintains and distributes open-source software and open data sets, and seeks to ensure that these are available under a free license in a long-term perspective.
In this paper, we presen…
▽ More
The Open Porous Media (OPM) initiative is a community effort that encourages open innovation and reproducible research for simulation of porous media processes. OPM coordinates collaborative software development, maintains and distributes open-source software and open data sets, and seeks to ensure that these are available under a free license in a long-term perspective.
In this paper, we present OPM Flow, which is a reservoir simulator developed for industrial use, as well as some of the individual components used to make OPM Flow. The descriptions apply to the 2019.10 release of OPM.
△ Less
Submitted 4 October, 2019;
originally announced October 2019.
-
WWPD elements of big mapping class groups
Authors:
Alexander J. Rasmussen
Abstract:
We study mapping class groups of infinite type surfaces with isolated punctures and their actions on the loop graphs introduced by Bavard-Walker. We classify all of the mapping classes in these actions which are loxodromic with a WWPD action on the corresponding loop graph. The WWPD property is a weakening of Bestvina-Fujiwara's weak proper discontinuity and is useful for constructing non-trivial…
▽ More
We study mapping class groups of infinite type surfaces with isolated punctures and their actions on the loop graphs introduced by Bavard-Walker. We classify all of the mapping classes in these actions which are loxodromic with a WWPD action on the corresponding loop graph. The WWPD property is a weakening of Bestvina-Fujiwara's weak proper discontinuity and is useful for constructing non-trivial quasimorphisms. We use this classification to give a sufficient criterion for subgroups of big mapping class groups to have infinite-dimensional second bounded cohomology and use this criterion to give simple proofs that certain natural subgroups of big mapping class groups have infinite-dimensional second bounded cohomology.
△ Less
Submitted 6 August, 2020; v1 submitted 14 September, 2019;
originally announced September 2019.
-
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.
-
Actions of solvable Baumslag-Solitar groups on hyperbolic metric spaces
Authors:
Carolyn R. Abbott,
Alexander J. Rasmussen
Abstract:
We give a complete list of the cobounded actions of solvable Baumslag-Solitar groups on hyperbolic metric spaces up to a natural equivalence relation. The set of equivalence classes carries a natural partial order first introduced by Abbott-Balasubramanya-Osin, and we describe the resulting poset completely. There are finitely many equivalence classes of actions, and each equivalence class contain…
▽ More
We give a complete list of the cobounded actions of solvable Baumslag-Solitar groups on hyperbolic metric spaces up to a natural equivalence relation. The set of equivalence classes carries a natural partial order first introduced by Abbott-Balasubramanya-Osin, and we describe the resulting poset completely. There are finitely many equivalence classes of actions, and each equivalence class contains the action on a point, a tree, or the hyperbolic plane.
△ Less
Submitted 2 November, 2021; v1 submitted 10 June, 2019;
originally announced June 2019.
-
Intrinsically interacting topological crystalline insulators and superconductors
Authors:
Alex Rasmussen,
Yuan-Ming Lu
Abstract:
Motivated by recent progress in crystalline symmetry protected topological (SPT) phases of interacting bosons, we study topological crystalline insulators/superconductors (TCIs) of strongly interacting fermions. We construct a class of intrinsically interacting fermionic TCIs, and show that they are beyond both free-fermion TCIs and bosonic crystalline SPT phases. We also show how these phases can…
▽ More
Motivated by recent progress in crystalline symmetry protected topological (SPT) phases of interacting bosons, we study topological crystalline insulators/superconductors (TCIs) of strongly interacting fermions. We construct a class of intrinsically interacting fermionic TCIs, and show that they are beyond both free-fermion TCIs and bosonic crystalline SPT phases. We also show how these phases can be characterized by symmetry protected gapless fermion modes on the corners/hinges of an open system.
△ Less
Submitted 29 October, 2018;
originally announced October 2018.
-
Classification and construction of higher-order symmetry protected topological phases of interacting bosons
Authors:
Alex Rasmussen,
Yuan-Ming Lu
Abstract:
Motivated by the recent discovery of higher-order topological insulators, we study their counterparts in strongly interacting bosons: `higher-order symmetry protected topological (HOSPT) phases'. While the usual (1st-order) SPT phases in d spatial dimensions support anomalous (d-1)-dimensional surface states, HOSPT phases in d dimensions are characterized by topological boundary states of dimensio…
▽ More
Motivated by the recent discovery of higher-order topological insulators, we study their counterparts in strongly interacting bosons: `higher-order symmetry protected topological (HOSPT) phases'. While the usual (1st-order) SPT phases in d spatial dimensions support anomalous (d-1)-dimensional surface states, HOSPT phases in d dimensions are characterized by topological boundary states of dimension (d-2) or smaller, protected by certain global symmetries and robust against disorders. Based on a dimensional reduction analysis, we show that HOSPT phases can be built from lower-dimensional SPT phases in a way that preserves the associated crystalline symmetries. When the total symmetry is a direct product of global and crystalline symmetry groups, we are able to classify the HOSPT phases using the Künneth formula of group cohomology. Based on a decorated domain wall picture of the Künneth formula, we show how to systematically construct the HOSPT phases, and demonstrate our construction with many examples in two and three dimensions.
△ Less
Submitted 28 October, 2018; v1 submitted 19 September, 2018;
originally announced September 2018.
-
Dyonic zero-energy modes
Authors:
Morten I. K. Munk,
Asbjørn Rasmussen,
Michele Burrello
Abstract:
One-dimensional systems with topological order are intimately related to the appearance of zero-energy modes localized on their boundaries. The most common example is the Kitaev chain, which displays Majorana zero-energy modes and it is characterized by a two-fold ground state degeneracy related to the global $\mathbb{Z}_2$ symmetry associated with fermionic parity. By extending the symmetry to th…
▽ More
One-dimensional systems with topological order are intimately related to the appearance of zero-energy modes localized on their boundaries. The most common example is the Kitaev chain, which displays Majorana zero-energy modes and it is characterized by a two-fold ground state degeneracy related to the global $\mathbb{Z}_2$ symmetry associated with fermionic parity. By extending the symmetry to the $\mathbb{Z}_N$ group, it is possible to engineer systems hosting topological parafermionic modes. In this work, we address one-dimensional systems with a generic discrete symmetry group $G$. We define a ladder model of gauge fluxes that generalizes the Ising and Potts models and displays a symmetry broken phase. Through a non-Abelian Jordan-Wigner transformation, we map this flux ladder into a model of dyonic operators, defined by the group elements and irreducible representations of $G$. We show that the so-obtained dyonic model has topological order, with zero-energy modes localized at its boundary. These dyonic zero-energy modes are in general weak topological modes, but strong dyonic zero modes appear when suitable position-dependent couplings are considered.
△ Less
Submitted 4 January, 2019; v1 submitted 24 July, 2018;
originally announced July 2018.
-
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.
-
Deconfined Quantum Critical Point on the Triangular Lattice
Authors:
Chao-Ming Jian,
Alex Thomson,
Alex Rasmussen,
Zhen Bi,
Cenke Xu
Abstract:
We first propose a topological term that captures the "intertwinement" between the standard "$\sqrt{3} \times \sqrt{3}$" antiferromagnetic order (or the so-called 120$^\circ$ state) and the "$\sqrt{12}\times \sqrt{12}$" valence solid bond (VBS) order for spin-1/2 systems on a triangular lattice. Then using a controlled renormalization group calculation, we demonstrate that there exists an unfine-t…
▽ More
We first propose a topological term that captures the "intertwinement" between the standard "$\sqrt{3} \times \sqrt{3}$" antiferromagnetic order (or the so-called 120$^\circ$ state) and the "$\sqrt{12}\times \sqrt{12}$" valence solid bond (VBS) order for spin-1/2 systems on a triangular lattice. Then using a controlled renormalization group calculation, we demonstrate that there exists an unfine-tuned direct continuous deconfined quantum critical point (dQCP) between the two ordered phases mentioned above. This dQCP is described by the $N_f = 4$ quantum electrodynamics (QED) with an emergent PSU(4)=SU(4)/$Z_4$ symmetry only at the critical point. The topological term aforementioned is also naturally derived from the $N_f = 4 $ QED. We also point out that physics around this dQCP is analogous to the boundary of a $3d$ bosonic symmetry protected topological state with on-site symmetries only.
△ Less
Submitted 12 October, 2017;
originally announced October 2017.
-
Emergent Symmetry and Tricritical Points near the deconfined Quantum Critical Point
Authors:
Chao-Ming Jian,
Alex Rasmussen,
Yi-Zhuang You,
Cenke Xu
Abstract:
Recent proposal of the duality between the $N=2$ noncompact QED$_3$ and the easy-plane noncompact CP$^1$ (NCCP$^1$) model suggests that the deconfined quantum critical point (dQCP) between the easy-plane antiferromagnet and the VBS order on the square lattice may have an emergent O(4) symmetry, due to the self-duality of the $N=2$ noncompact QED$_3$. Recent numerical progresses suggest that this e…
▽ More
Recent proposal of the duality between the $N=2$ noncompact QED$_3$ and the easy-plane noncompact CP$^1$ (NCCP$^1$) model suggests that the deconfined quantum critical point (dQCP) between the easy-plane antiferromagnet and the VBS order on the square lattice may have an emergent O(4) symmetry, due to the self-duality of the $N=2$ noncompact QED$_3$. Recent numerical progresses suggest that this easy-plane dQCP does exist and it has an emergent O(4) symmetry. But for the O(4) symmetry to really emerge at the dQCP, certain O(4) symmetry breaking perturbations need to be irrelevant at the putative O(4) fixed point. It is more convenient to study these symmetry breaking perturbations in the $N=2$ noncompact QED$_3$. We demonstrate that a natural large-$N$ generalization and a controlled $1/N$ expansion supports the stability of the O(4) fixed point against the symmetry breaking perturbations. We also develop the theory for two tricritical points close to the easy-plane dQCP. One tricritical point is between the dQCP and a {\it self-dual} $Z_2$ topological order; the other is the tricritical point that connects the continuous dQCP and a first order Néel-VBS transition, motivated by recent numerical results.
△ Less
Submitted 9 August, 2017;
originally announced August 2017.
-
Uniform hyperbolicity of the graphs of nonseparating curves via bicorn curves
Authors:
Alexander J. Rasmussen
Abstract:
We show that the graphs of nonseparating curves for oriented finite type surfaces are uniformly hyperbolic. Our proof follows the proof of uniform hyperbolicity of the graphs of curves for closed surfaces due to Przytycki-Sisto, while introducing new arguments using homology to certify that certain curves are nonseparating. As demonstrated by Aramayona-Valdez, this proves also that the graph of no…
▽ More
We show that the graphs of nonseparating curves for oriented finite type surfaces are uniformly hyperbolic. Our proof follows the proof of uniform hyperbolicity of the graphs of curves for closed surfaces due to Przytycki-Sisto, while introducing new arguments using homology to certify that certain curves are nonseparating. As demonstrated by Aramayona-Valdez, this proves also that the graph of nonseparating curves for any oriented infinite type surface with finite positive genus is hyperbolic.
△ Less
Submitted 2 April, 2020; v1 submitted 26 July, 2017;
originally announced July 2017.
-
Linear feature detection algorithm for astronomical surveys - II. Defocusing effects on meteor tracks
Authors:
Dino Bektešević,
Dejan Vinković,
Andrew Rasmussen,
Željko Ivezić
Abstract:
Given the current limited knowledge of meteor plasma micro-physics and its interaction with the surrounding atmosphere and ionosphere, meteors are a highly interesting observational target for high-resolution wide-field astronomical surveys. Such surveys are capable of resolving the physical size of meteor plasma heads, but they produce large volumes of images that need to be automatically inspect…
▽ More
Given the current limited knowledge of meteor plasma micro-physics and its interaction with the surrounding atmosphere and ionosphere, meteors are a highly interesting observational target for high-resolution wide-field astronomical surveys. Such surveys are capable of resolving the physical size of meteor plasma heads, but they produce large volumes of images that need to be automatically inspected for possible existence of long linear features produced by meteors. Here we show how big aperture sky survey telescopes detect meteors as defocused tracks with a central brightness depression. We derive an analytic expression for a defocused point source meteor track and use it to calculate brightness profiles of meteors modeled as uniform brightness disks. We apply our modeling to meteor images as seen by the SDSS and LSST telescopes. The expression is validated by Monte Carlo ray-tracing simulations of photons traveling through the atmosphere and the LSST telescope optics. We show that estimates of the meteor distance and size can be extracted from the measured FWHM and the strength of the central dip in the observed brightness profile. However, this extraction becomes difficult when the defocused meteor track is distorted by the atmospheric seeing or contaminated by a long lasting glowing meteor trail. The FWHM of satellites tracks is distinctly narrower than meteor values, which enables removal of a possible confusion between satellites and meteors.
△ Less
Submitted 22 July, 2017;
originally announced July 2017.
-
Gapless Topological Order, Gravity, and Black Holes
Authors:
Alex Rasmussen,
Adam Jermyn
Abstract:
In this work we demonstrate that linearized gravity exhibits gapless topological order with an extensive ground state degeneracy. This phenomenon is closely related both to the topological order of the pyrochlore U(1) spin liquid and to recent work by Hawking et. al. who used the soft photon and graviton theorems to demonstrate that the vacuum in linearized gravity is not unique. We first consider…
▽ More
In this work we demonstrate that linearized gravity exhibits gapless topological order with an extensive ground state degeneracy. This phenomenon is closely related both to the topological order of the pyrochlore U(1) spin liquid and to recent work by Hawking et. al. who used the soft photon and graviton theorems to demonstrate that the vacuum in linearized gravity is not unique. We first consider lattice models whose low-energy behavior are described by electromagnetism and linearized gravity, and then argue that the topological nature of these models carries over into the continuum. We demonstrate that these models can have many ground states without making assumptions about the topology of spacetime or about the high-energy nature of the theory, and show that the infinite family of symmetries described by Hawking et. al. are simply the difierent topological sectors. We argue that in this context black holes appear as topological defects in the IR theory, and that this suggests a potential approach to understanding both the firewall paradox and information encoding in gravitational theories. Finally, we use insights from the soft boson theorems to make connections between deconfined gauge theories with continuous gauge groups and gapless topological order.
△ Less
Submitted 28 August, 2018; v1 submitted 14 March, 2017;
originally announced March 2017.
-
Estimation of Received Signal Strength Distribution for Smart Meters with Biased Measurement Data Set
Authors:
Mathias Rønholt Kielgast,
Anders Charly Rasmussen,
Mathias Hjorth Laursen,
Jimmy Jessen Nielsen,
Petar Popovski,
Rasmus Krigslund
Abstract:
This letter presents an experimental study and a novel modelling approach of the wireless channel of smart utility meters placed in basements or sculleries. The experimental data consist of signal strength measurements of consumption report packets. Since such packets are only registered if they can be decoded by the receiver, the part of the signal strength distribution that falls below the recei…
▽ More
This letter presents an experimental study and a novel modelling approach of the wireless channel of smart utility meters placed in basements or sculleries. The experimental data consist of signal strength measurements of consumption report packets. Since such packets are only registered if they can be decoded by the receiver, the part of the signal strength distribution that falls below the receiver sensitivity threshold is not observable. We combine a Rician fading model with a bias function that captures the cut-off in the observed signal strength measurements. Two sets of experimental data are analysed. It is shown that the proposed method offers an approximation of the distribution of the signal strength measurements that is better than a naïve Rician fitting.
△ Less
Submitted 17 October, 2016;
originally announced October 2016.
-
Majorana box qubits
Authors:
Stephan Plugge,
Asbjørn Rasmussen,
Reinhold Egger,
Karsten Flensberg
Abstract:
Quantum information protected by the topology of the storage medium is expected to exhibit long coherence times. Another feature are topologically protected gates generated through braiding of Majorana bound states. However, braiding requires structures with branched topological segments which have inherent difficulties in the semiconductor-superconductor heterostructures now believed to host Majo…
▽ More
Quantum information protected by the topology of the storage medium is expected to exhibit long coherence times. Another feature are topologically protected gates generated through braiding of Majorana bound states. However, braiding requires structures with branched topological segments which have inherent difficulties in the semiconductor-superconductor heterostructures now believed to host Majorana bound states. In this paper, we construct quantum bits taking advantage of the topological protection and non-local properties of Majorana bound states in a network of parallel wires, but without relying on braiding for quantum gates. The elementary unit is made from three topological wires, two wires coupled by a trivial superconductor and the third acting as an interference arm. Coulomb blockade of the combined wires spawns a fractionalized spin, non-locally addressable by quantum dots used for single-qubit readout, initialization, and manipulation. We describe how the same tools allow for measurement-based implementation of the Clifford gates, in total making the architecture universal. Proof-of-principle demonstration of topologically protected qubits using existing techniques is therefore within reach.
△ Less
Submitted 21 December, 2016; v1 submitted 6 September, 2016;
originally announced September 2016.
-
High fidelity point-spread function retrieval in the presence of electrostatic, hysteretic pixel response
Authors:
Andrew Rasmussen,
Augustin Guyonnet,
Craig Lage,
Pierre Antilogus,
Pierre Astier,
Peter Doherty,
Kirk Gilmore,
Ivan Kotov,
Robert Lupton,
Andrei Nomerotski,
Paul O'Connor,
Christopher Stubbs,
Anthony Tyson,
Christopher Walter
Abstract:
We employ electrostatic conversion drift calculations to match CCD pixel signal covariances observed in flat field exposures acquired using candidate sensor devices for the LSST Camera. We thus constrain pixel geometry distortions present at the end of integration, based on signal images recorded. We use available data from several operational voltage parameter settings to validate our understandi…
▽ More
We employ electrostatic conversion drift calculations to match CCD pixel signal covariances observed in flat field exposures acquired using candidate sensor devices for the LSST Camera. We thus constrain pixel geometry distortions present at the end of integration, based on signal images recorded. We use available data from several operational voltage parameter settings to validate our understanding. Our primary goal is to optimize flux point-spread function (FPSF) estimation quantitatively, and thereby minimize sensor-induced errors which may limit performance in precision astronomy applications. We consider alternative compensation scenarios that will take maximum advantage of our understanding of this underlying mechanism in data processing pipelines currently under development.
To quantitatively capture the pixel response in high-contrast/high dynamic range operational extrema, we propose herein some straightforward laboratory tests that involve altering the time order of source illumination on sensors, within individual test exposures. Hence the word {\it hysteretic} in the title of this paper.
△ Less
Submitted 12 September, 2016; v1 submitted 5 August, 2016;
originally announced August 2016.
-
Stable Interacting (2 + 1)d Conformal Field Theories at the Boundary of a class of (3 + 1)d Symmetry Protected Topological Phases
Authors:
Zhen Bi,
Alex Rasmussen,
Yoni BenTov,
Cenke Xu
Abstract:
Motivated by recent studies of symmetry protected topological (SPT) phases, we explore the possible gapless quantum disordered phases in the $(2+1)d$ nonlinear sigma model defined on the Grassmannian manifold $\frac{U(N)}{U(n)\times U(N - n)}$ with a Wess-Zumino-Witten (WZW) term at level $k$, which is the effective low energy field theory of the boundary of certain $(3+1)d$ SPT states. With…
▽ More
Motivated by recent studies of symmetry protected topological (SPT) phases, we explore the possible gapless quantum disordered phases in the $(2+1)d$ nonlinear sigma model defined on the Grassmannian manifold $\frac{U(N)}{U(n)\times U(N - n)}$ with a Wess-Zumino-Witten (WZW) term at level $k$, which is the effective low energy field theory of the boundary of certain $(3+1)d$ SPT states. With $k = 0$, this model has a well-controlled large-$N$ limit, $i.e.$ its renormalization group equations can be computed exactly with large-$N$. However, with the WZW term, the large-$N$ and large-$k$ limit alone is not sufficient for a reliable study of the nature of the quantum disordered phase. We demonstrate that through a combined large-$N$, large-$k$ and $ε-$generalization, a stable fixed point in the quantum disordered phase can be reliably located in the large$-N$ limit and leading order $ε-$expansion, which corresponds to a $(2+1)d$ strongly interacting conformal field theory.
△ Less
Submitted 30 May, 2016; v1 submitted 17 May, 2016;
originally announced May 2016.
-
Defect Tolerant Monolayer Transition Metal Dichalcogenides
Authors:
Mohnish Pandey,
Filip A. Rasmussen,
Korina Kuhar,
Thomas Olsen,
Karsten W. Jacobsen,
Kristian S. Thygesen
Abstract:
Localized electronic states formed inside the band gap of a semiconductor due to crystal defects can be detrimental to the material's optoelectronic properties. Semiconductors with lower tendency to form defect induced deep gap states are termed defect tolerant. Here we provide a systematic first principles investigation of defect tolerance in 29 monolayer transition metal dichalcogenides (TMDs) o…
▽ More
Localized electronic states formed inside the band gap of a semiconductor due to crystal defects can be detrimental to the material's optoelectronic properties. Semiconductors with lower tendency to form defect induced deep gap states are termed defect tolerant. Here we provide a systematic first principles investigation of defect tolerance in 29 monolayer transition metal dichalcogenides (TMDs) of interest for nanoscale optoelectronics. We find that the TMDs based on group VI and X metals form deep gap states upon creation of a chalcogen (S, Se, Te) vacancy while the TMDs based on group IV metals form only shallow defect levels and are thus predicted to be defect tolerant. Interestingly, all the defect sensitive TMDs have valence and conduction bands with very similar orbital composition. This indicates a bonding/anti-bonding nature of the gap which in turn suggests that dangling bonds will fall inside the gap. These ideas are made quantitative by introducing a descriptor that measures the degree of similarity of the conduction and valence band manifolds. Finally, the study is generalized to non-polar nanoribbons of the TMDs where we find that only the defect sensitive materials form edge states within the band gap.
△ Less
Submitted 11 April, 2016;
originally announced April 2016.
-
Stable Gapless Bose Liquid Phases without any Symmetry
Authors:
Alex Rasmussen,
Yi-Zhuang You,
Cenke Xu
Abstract:
It is well-known that a stable algebraic spin liquid state (or equivalently an algebraic Bose liquid (ABL) state) with emergent gapless photon excitations can exist in quantum spin ice systems, or in a quantum dimer model on a bipartite $3d$ lattice. This photon phase is stable against any weak perturbation without assuming any symmetry. Further works concluded that certain lattice models give ris…
▽ More
It is well-known that a stable algebraic spin liquid state (or equivalently an algebraic Bose liquid (ABL) state) with emergent gapless photon excitations can exist in quantum spin ice systems, or in a quantum dimer model on a bipartite $3d$ lattice. This photon phase is stable against any weak perturbation without assuming any symmetry. Further works concluded that certain lattice models give rise to more exotic stable algebraic Bose liquid phases with graviton-like excitations. In this paper we will show how these algebraic Bose liquid states can be generalized to stable phases with even more exotic types of gapless excitations and then argue that these new phases are stable against weak perturbations. We also explicitly show that these theories have an (algebraic) topological ground state degeneracy on a torus, and construct the corresponding topological invariants.
△ Less
Submitted 29 January, 2016;
originally announced January 2016.
-
Efficient many-body calculations of 2D materials using exact limits for the screened potential: Band gaps of MoS$_2$, hBN, and phosphorene
Authors:
Filip A. Rasmussen,
Per S. Schmidt,
Kirsten T. Winther,
Kristian S. Thygesen
Abstract:
Calculating the quasiparticle (QP) band structure of two-dimensional (2D) materials within the GW self-energy approximation has proven to be a rather demanding computational task. The main reason is the strong $\mathbf{q}$-dependence of the 2D dielectric function around $\mathbf{q} = \mathbf{0}$ that calls for a much denser sampling of the Brillouin zone than is necessary for similar 3D solids. He…
▽ More
Calculating the quasiparticle (QP) band structure of two-dimensional (2D) materials within the GW self-energy approximation has proven to be a rather demanding computational task. The main reason is the strong $\mathbf{q}$-dependence of the 2D dielectric function around $\mathbf{q} = \mathbf{0}$ that calls for a much denser sampling of the Brillouin zone than is necessary for similar 3D solids. Here we use an analytical expression for the small $\mathbf{q}$-limit of the 2D response function to perform the BZ integral over the critical region around $\mathbf{q} = \mathbf{0}$. This drastically reduces the requirements on the $\mathbf{q}$-point mesh and implies a significant computational speed-up. For example, in the case of monolayer MoS$_2$, convergence of the $G_0W_0$ band gap to within $\sim 0.1\,\mathrm{eV}$ is achieved with $12\times 12$ $\mathbf{q}$-points rather than the $36\times 36$ mesh required with discrete BZ sampling techniques. We perform a critical assessment of the band gap of the three prototypical 2D semiconductors MoS$_2$, hBN, and phosphorene including the effect of self-consistency at the GW$_0$ and GW level. The method is implemented in the open source GPAW code.
△ Less
Submitted 19 June, 2016; v1 submitted 31 October, 2015;
originally announced November 2015.