Showing 1–50 of 70 results for author: Larsen, P

Exploring the Core-galaxy Connection

Authors: Isabele Souza Vitório, Michael Buehlmann, Eve Kovacs, Patricia Larsen, Nicholas Frontiere, Katrin Heitmann

Abstract: Halo core tracking is a novel concept designed to efficiently follow halo substructure in large simulations. We have recently developed this concept in gravity-only simulations to investigate the galaxy-halo connection in the context of empirical and semi-analytic models. Here, we incorporate information from hydrodynamics simulations, with an emphasis on establishing a connection between cores an… ▽ More Halo core tracking is a novel concept designed to efficiently follow halo substructure in large simulations. We have recently developed this concept in gravity-only simulations to investigate the galaxy-halo connection in the context of empirical and semi-analytic models. Here, we incorporate information from hydrodynamics simulations, with an emphasis on establishing a connection between cores and galaxies. We compare cores across gravity-only, adiabatic hydrodynamics, and subgrid hydrodynamics simulations with the same initial phases. We demonstrate that cores are stable entities whose halo-centric radial profiles match across the simulations. We further develop a methodology that uses merging and infall mass cuts to group cores in the hydrodynamics simulation, creating on average, a one-to-one match to corresponding galaxies. We apply this methodology to cores from the gravity-only simulation, thus creating a proxy for galaxies which approximate the populations from the hydrodynamics simulation. Our results pave the way to incorporate inputs from smaller-scale hydrodynamics simulations directly into large-scale gravity-only runs in a principled manner. △ Less

Submitted 3 July, 2024; v1 submitted 28 June, 2024; originally announced July 2024.

Comments: 14 pages, 13 figures

Pathological Regularization Regimes in Classification Tasks

Authors: Maximilian Wiesmann, Paul Larsen

Abstract: In this paper we demonstrate the possibility of a trend reversal in binary classification tasks between the dataset and a classification score obtained from a trained model. This trend reversal occurs for certain choices of the regularization parameter for model training, namely, if the parameter is contained in what we call the pathological regularization regime. For ridge regression, we give nec… ▽ More In this paper we demonstrate the possibility of a trend reversal in binary classification tasks between the dataset and a classification score obtained from a trained model. This trend reversal occurs for certain choices of the regularization parameter for model training, namely, if the parameter is contained in what we call the pathological regularization regime. For ridge regression, we give necessary and sufficient algebraic conditions on the dataset for the existence of a pathological regularization regime. Moreover, our results provide a data science practitioner with a hands-on tool to avoid hyperparameter choices suffering from trend reversal. We furthermore present numerical results on pathological regularization regimes for logistic regression. Finally, we draw connections to datasets exhibiting Simpson's paradox, providing a natural source of pathological datasets. △ Less

Submitted 20 June, 2024; originally announced June 2024.

MSC Class: 62J07 (Primary); 62J12; 62R07 (Secondary)

Phase Transitions in Quantum Many-Body Scars

Authors: Peter Græns Larsen, Anne E. B. Nielsen

Abstract: We propose a type of phase transition in quantum many-body systems, which occurs in highly excited quantum many-body scar states, while the rest of the spectrum is largely unaffected. Such scar state phase transitions can be realized by embedding a matrix product state, known to undergo a phase transition, as a scar state into the thermal spectrum of a parent Hamiltonian. We find numerically that… ▽ More We propose a type of phase transition in quantum many-body systems, which occurs in highly excited quantum many-body scar states, while the rest of the spectrum is largely unaffected. Such scar state phase transitions can be realized by embedding a matrix product state, known to undergo a phase transition, as a scar state into the thermal spectrum of a parent Hamiltonian. We find numerically that the mechanism for the scar state phase transition involves the formation or presence of low-entropy states at energies similar to the scar state in the vicinity of the phase transition point. △ Less

Submitted 30 May, 2024; originally announced May 2024.

Comments: 5 pages, 2 figures

State of the Art Report: Verified Computation

Authors: Jim Woodcock, Mikkel Schmidt Andersen, Diego F. Aranha, Stefan Hallerstede, Simon Thrane Hansen, Nikolaj Kuhne Jakobsen, Tomas Kulik, Peter Gorm Larsen, Hugo Daniel Macedo, Carlos Ignacio Isasa Martin, Victor Alexander Mtsimbe Norrild

Abstract: This report describes the state of the art in verifiable computation. The problem being solved is the following: The Verifiable Computation Problem (Verifiable Computing Problem) Suppose we have two computing agents. The first agent is the verifier, and the second agent is the prover. The verifier wants the prover to perform a computation. The verifier sends a description of the computation to t… ▽ More This report describes the state of the art in verifiable computation. The problem being solved is the following: The Verifiable Computation Problem (Verifiable Computing Problem) Suppose we have two computing agents. The first agent is the verifier, and the second agent is the prover. The verifier wants the prover to perform a computation. The verifier sends a description of the computation to the prover. Once the prover has completed the task, the prover returns the output to the verifier. The output will contain proof. The verifier can use this proof to check if the prover computed the output correctly. The check is not required to verify the algorithm used in the computation. Instead, it is a check that the prover computed the output using the computation specified by the verifier. The effort required for the check should be much less than that required to perform the computation. This state-of-the-art report surveys 128 papers from the literature comprising more than 4,000 pages. Other papers and books were surveyed but were omitted. The papers surveyed were overwhelmingly mathematical. We have summarised the major concepts that form the foundations for verifiable computation. The report contains two main sections. The first, larger section covers the theoretical foundations for probabilistically checkable and zero-knowledge proofs. The second section contains a description of the current practice in verifiable computation. Two further reports will cover (i) military applications of verifiable computation and (ii) a collection of technical demonstrators. The first of these is intended to be read by those who want to know what applications are enabled by the current state of the art in verifiable computation. The second is for those who want to see practical tools and conduct experiments themselves. △ Less

Submitted 16 February, 2024; v1 submitted 29 August, 2023; originally announced August 2023.

Comments: 54 pages

Optimization and Quality Assessment of Baryon Pasting for Intracluster Gas using the Borg Cube Simulation

Authors: F. Kéruzoré, L. E. Bleem, M. Buehlmann, J. D. Emberson, N. Frontiere, S. Habib, K. Heitmann, P. Larsen

Abstract: Synthetic datasets generated from large-volume gravity-only simulations are an important tool in the calibration of cosmological analyses. Their creation often requires accurate inference of baryonic observables from the dark matter field. We explore the effectiveness of a baryon pasting algorithm in providing precise estimations of three-dimensional gas thermodynamic properties based on gravity-o… ▽ More Synthetic datasets generated from large-volume gravity-only simulations are an important tool in the calibration of cosmological analyses. Their creation often requires accurate inference of baryonic observables from the dark matter field. We explore the effectiveness of a baryon pasting algorithm in providing precise estimations of three-dimensional gas thermodynamic properties based on gravity-only simulations. We use the Borg Cube, a pair of simulations originating from identical initial conditions, with one run evolved as a gravity-only simulation, and the other incorporating non-radiative hydrodynamics. Matching halos in both simulations enables comparisons of gas properties on an individual halo basis. This comparative analysis allows us to fit for the model parameters that yield the closest agreement between the gas properties in both runs. To capture the redshift evolution of these parameters, we perform the analysis at five distinct redshift steps, spanning from $z=0$ to $2$. We find that the investigated algorithm, utilizing information solely from the gravity-only simulation, achieves few-percent accuracy in reproducing the median intracluster gas pressure and density, albeit with a scatter of approximately 20%, for cluster-scale objects up to $z=2$. We measure the scaling relation between integrated Compton parameter and cluster mass ($Y_{500c} | M_{500c}$), and find that the imprecision of baryon pasting adds less than 5% to the intrinsic scatter measured in the hydrodynamic simulation. We provide best-fitting values and their redshift evolution, and discuss future investigations that will be undertaken to extend this work. △ Less

Submitted 27 November, 2023; v1 submitted 23 June, 2023; originally announced June 2023.

Comments: 14 pages, 8 figures, 3 tables; accepted in the Open Journal of Astrophysics

Digital Twin as a Service (DTaaS): A Platform for Digital Twin Developers and Users

Authors: Prasad Talasila, Cláudio Gomes, Peter Høgh Mikkelsen, Santiago Gil Arboleda, Eduard Kamburjan, Peter Gorm Larsen

Abstract: Establishing digital twins is a non-trivial endeavour especially when users face significant challenges in creating them from scratch. Ready availability of reusable models, data and tool assets, can help with creation and use of digital twins. A number of digital twin frameworks exist to facilitate creation and use of digital twins. In this paper we propose a digital twin framework to author digi… ▽ More Establishing digital twins is a non-trivial endeavour especially when users face significant challenges in creating them from scratch. Ready availability of reusable models, data and tool assets, can help with creation and use of digital twins. A number of digital twin frameworks exist to facilitate creation and use of digital twins. In this paper we propose a digital twin framework to author digital twin assets, create digital twins from reusable assets and make the digital twins available as a service to other users. The proposed framework automates the management of reusable assets, storage, provision of compute infrastructure, communication and monitoring tasks. The users operate at the level of digital twins and delegate rest of the work to the digital twin as a service framework. △ Less

Submitted 13 June, 2023; v1 submitted 12 May, 2023; originally announced May 2023.

Comments: 8 pages, 6 figures. Accepted at Digital Twin 2023

ACM Class: D.2.11

Model-Based Monitoring and State Estimation for Digital Twins: The Kalman Filter

Authors: Hao Feng, Cláudio Gomes, Peter Gorm Larsen

Abstract: A digital twin (DT) monitors states of the physical twin (PT) counterpart and provides a number of benefits such as advanced visualizations, fault detection capabilities, and reduced maintenance cost. It is the ability to be able to detect the states inside the DT that enable such benefits. In order to estimate the desired states of a PT, we propose the use of a Kalman Filter (KF). In this tutoria… ▽ More A digital twin (DT) monitors states of the physical twin (PT) counterpart and provides a number of benefits such as advanced visualizations, fault detection capabilities, and reduced maintenance cost. It is the ability to be able to detect the states inside the DT that enable such benefits. In order to estimate the desired states of a PT, we propose the use of a Kalman Filter (KF). In this tutorial, we provide an introduction and detailed derivation of the KF. We demonstrate the use of KF to monitor and anomaly detection through an incubator system. Our experimental result shows that KF successfully can detect the anomaly during monitoring. △ Less

Submitted 29 April, 2023; originally announced May 2023.

Comments: 14 pages

Bidirectional UML Visualisation of VDM Models

Authors: Jonas Lund, Lucas Bjarke Jensen, Nick Battle, Peter Gorm Larsen, Hugo Daniel Macedo

Abstract: The VDM-PlantUML Plugin enables translations between the text based UML tool PlantUML and VDM++ and has been released as a part of the VDM VSCode extension. This enhances already extensive feature-set of VDM VSCode with support for UML. The link between VDM and UML is thoroughly described with a set of translation rules that serve as the base of the implementation of the translation plugin. This i… ▽ More The VDM-PlantUML Plugin enables translations between the text based UML tool PlantUML and VDM++ and has been released as a part of the VDM VSCode extension. This enhances already extensive feature-set of VDM VSCode with support for UML. The link between VDM and UML is thoroughly described with a set of translation rules that serve as the base of the implementation of the translation plugin. This is however still an early rendition of the plugin with limited usability due to the loss of information between translations and a lack of workflow optimisations, which we plan to solve in the future. △ Less

Submitted 13 April, 2023; originally announced April 2023.

VDM recursive functions in Isabelle/HOL

Authors: Leo Freitas, Peter Gorm Larsen

Abstract: For recursive functions general principles of induction needs to be applied. Instead of verifying them directly using the Vienna Development Method Specification Language (VDM-SL), we suggest a translation to Isabelle/HOL. In this paper, the challenges of such a translation for recursive functions are presented. This is an extension of an existing translation and a VDM mathematical toolbox in Isab… ▽ More For recursive functions general principles of induction needs to be applied. Instead of verifying them directly using the Vienna Development Method Specification Language (VDM-SL), we suggest a translation to Isabelle/HOL. In this paper, the challenges of such a translation for recursive functions are presented. This is an extension of an existing translation and a VDM mathematical toolbox in Isabelle/HOL enabling support for recursive functions. △ Less

Submitted 30 March, 2023; originally announced March 2023.

Modelling Chess in VDM++

Authors: Morten Haahr Kristensen, Peter Gorm Larsen

Abstract: The game of chess is well-known and widely played all over the world. However, the rules for playing it are rather complex since there are different types of pieces and the ways they are allowed to move depend upon the type of the piece. In this paper we discuss alternative paradigms that can be used for modelling the rule of the chess game using VDM++ and show what we believe is the best model. I… ▽ More The game of chess is well-known and widely played all over the world. However, the rules for playing it are rather complex since there are different types of pieces and the ways they are allowed to move depend upon the type of the piece. In this paper we discuss alternative paradigms that can be used for modelling the rule of the chess game using VDM++ and show what we believe is the best model. It is also illustrated how this model can be connected to a standard textual notation for the moves in a chess game. This can be used to combine the formal model to a more convenient interface. △ Less

Submitted 18 March, 2023; originally announced March 2023.

The catalog-to-cosmology framework for weak lensing and galaxy clustering for LSST

Authors: J. Prat, J. Zuntz, Y. Omori, C. Chang, T. Tröster, E. Pedersen, C. García-García, E. Phillips-Longley, J. Sanchez, D. Alonso, X. Fang, E. Gawiser, K. Heitmann, M. Ishak, M. Jarvis, E. Kovacs, P. Larsen, Y. -Y. Mao, L. Medina Varela, M. Paterno, S. D. Vitenti, Z. Zhang, The LSST Dark Energy Science Collaboration

Abstract: We present TXPipe, a modular, automated and reproducible pipeline for ingesting catalog data and performing all the calculations required to obtain quality-assured two-point measurements of lensing and clustering, and their covariances, with the metadata necessary for parameter estimation. The pipeline is developed within the Rubin Observatory Legacy Survey of Space and Time (LSST) Dark Energy Sci… ▽ More We present TXPipe, a modular, automated and reproducible pipeline for ingesting catalog data and performing all the calculations required to obtain quality-assured two-point measurements of lensing and clustering, and their covariances, with the metadata necessary for parameter estimation. The pipeline is developed within the Rubin Observatory Legacy Survey of Space and Time (LSST) Dark Energy Science Collaboration (DESC), and designed for cosmology analyses using LSST data. In this paper, we present the pipeline for the so-called 3x2pt analysis -- a combination of three two-point functions that measure the auto- and cross-correlation between galaxy density and shapes. We perform the analysis both in real and harmonic space using TXPipe and other LSST-DESC tools. We validate the pipeline using Gaussian simulations and show that it accurately measures data vectors and recovers the input cosmology to the accuracy level required for the first year of LSST data under this simplified scenario. We also apply the pipeline to a realistic mock galaxy sample extracted from the CosmoDC2 simulation suite (Korytov et al. 2019). TXPipe establishes a baseline framework that can be built upon as the LSST survey proceeds. Furthermore, the pipeline is designed to be easily extended to science probes beyond the 3x2pt analysis. △ Less

Submitted 21 April, 2023; v1 submitted 19 December, 2022; originally announced December 2022.

Comments: 20+11 pages, 10+11 figures. Version accepted in The Open Journal of Astrophysics

Comparison between the HUBCAP and DIGITBrain Platforms for Model-Based Design and Evaluation of Digital Twins

Authors: Prasad Talasila, Daniel-Cristian Crăciunean, Pirvu Bogdan-Constantin, Peter Gorm Larsen, Constantin Zamfirescu, Alea Scovill

Abstract: Digital twin technology is an essential approach to managing the lifecycle of industrial products. Among the many approaches used to manage digital twins, co-simulation has proven to be a reliable one. There have been multiple attempts to create collaborative and sustainable platforms for management of digital twins. This paper compares two such platforms, namely the HUBCAP and the DIGITbrain. Bot… ▽ More Digital twin technology is an essential approach to managing the lifecycle of industrial products. Among the many approaches used to manage digital twins, co-simulation has proven to be a reliable one. There have been multiple attempts to create collaborative and sustainable platforms for management of digital twins. This paper compares two such platforms, namely the HUBCAP and the DIGITbrain. Both these platforms have been and continue to be used among a stable group of researchers and industrial product manufacturers of digital twin technologies. This comparison of the HUBCAP and the DIGITbrain platforms is illustrated with an example use case of industrial factory to be used for manufacturing of agricultural robots. △ Less

Submitted 15 December, 2022; originally announced December 2022.

Computational exfoliation of atomically thin 1D materials with application to Majorana bound states

Authors: Hadeel Moustafa, Peter Mahler Larsen, Morten N. Gjerding, Jens Jørgen Mortensen, Kristian S. Thygesen, Karsten W. Jacobsen

Abstract: We introduce a computational database with calculated structural, thermodynamic, electronic, magnetic, and optical properties of 820 one-dimensional materials. The materials are systematically selected and exfoliated from experimental databases of crystal structures based on a dimensionality scoring parameter. The database is furthermore expanded by chemical element substitution in the materials.… ▽ More We introduce a computational database with calculated structural, thermodynamic, electronic, magnetic, and optical properties of 820 one-dimensional materials. The materials are systematically selected and exfoliated from experimental databases of crystal structures based on a dimensionality scoring parameter. The database is furthermore expanded by chemical element substitution in the materials. The materials are investigated in both their bulk form and as isolated one-dimensional components. We discuss the methodology behind the database, give an overview of some of the calculated properties, and look at patterns and correlations in the data. The database is furthermore applied in computational screening to identify materials, which could exhibit Majorana bound states. △ Less

Submitted 1 April, 2022; originally announced April 2022.

Comments: 17 pages, 20 figures

Transitioning from Stage-III to Stage-IV: Cosmology from galaxy$\times$CMB lensing and shear$\times$CMB lensing

Authors: Zhuoqi Zhang, Chihway Chang, Patricia Larsen, Lucas F. Secco, Joe Zuntz, the LSST Dark Energy Science Collaboration

Abstract: We examine the cosmological constraining power from two cross-correlation probes between galaxy and CMB surveys: the cross-correlation of lens galaxy density with CMB lensing convergence $\langleδκ\rangle$, and source galaxy weak lensing shear with CMB lensing convergence $\langleγκ\rangle$. These two cross-correlation probes provide an independent cross-check of other large-scale structure constr… ▽ More We examine the cosmological constraining power from two cross-correlation probes between galaxy and CMB surveys: the cross-correlation of lens galaxy density with CMB lensing convergence $\langleδκ\rangle$, and source galaxy weak lensing shear with CMB lensing convergence $\langleγκ\rangle$. These two cross-correlation probes provide an independent cross-check of other large-scale structure constraints and are insensitive to galaxy-only or CMB-only systematic effects. In addition, when combined with other large-scale structure probes, the cross-correlations can break degeneracies in cosmological and nuisance parameters, improving both the precision and robustness of the analysis. In this work, we study how the constraining power of $\langleδκ\rangle+\langleγκ\rangle$ changes from Stage-III (ongoing) to Stage-IV (future) surveys. Given the flexibility in selecting the lens galaxy sample, we also explore systematically the impact on cosmological constraints when we vary the redshift range and magnitude limit of the lens galaxies using mock galaxy catalogs. We find that in our setup, the contribution to cosmological constraints from $\langleδκ\rangle$ and $\langleγκ\rangle$ are comparable in the Stage-III datasets; but in Stage-IV surveys, the noise in $\langleδκ\rangle$ becomes subdominant to cosmic variance, preventing $\langleδκ\rangle$ to further improve the constraints. This implies that to maximize the cosmological constraints from future $\langleδκ\rangle+\langleγκ\rangle$ analyses, we should focus more on the requirements on $\langleγκ\rangle$ instead of $\langleδκ\rangle$. Furthermore, the selection of the lens sample should be optimized in terms of our ability to characterize its redshift or galaxy bias instead of its number density. △ Less

Submitted 3 September, 2023; v1 submitted 8 November, 2021; originally announced November 2021.

Comments: 18 pages, 13 figures

Constructing Neural Network-Based Models for Simulating Dynamical Systems

Authors: Christian Møldrup Legaard, Thomas Schranz, Gerald Schweiger, Ján Drgoňa, Basak Falay, Cláudio Gomes, Alexandros Iosifidis, Mahdi Abkar, Peter Gorm Larsen

Abstract: Dynamical systems see widespread use in natural sciences like physics, biology, chemistry, as well as engineering disciplines such as circuit analysis, computational fluid dynamics, and control. For simple systems, the differential equations governing the dynamics can be derived by applying fundamental physical laws. However, for more complex systems, this approach becomes exceedingly difficult. D… ▽ More Dynamical systems see widespread use in natural sciences like physics, biology, chemistry, as well as engineering disciplines such as circuit analysis, computational fluid dynamics, and control. For simple systems, the differential equations governing the dynamics can be derived by applying fundamental physical laws. However, for more complex systems, this approach becomes exceedingly difficult. Data-driven modeling is an alternative paradigm that seeks to learn an approximation of the dynamics of a system using observations of the true system. In recent years, there has been an increased interest in data-driven modeling techniques, in particular neural networks have proven to provide an effective framework for solving a wide range of tasks. This paper provides a survey of the different ways to construct models of dynamical systems using neural networks. In addition to the basic overview, we review the related literature and outline the most significant challenges from numerical simulations that this modeling paradigm must overcome. Based on the reviewed literature and identified challenges, we provide a discussion on promising research areas. △ Less

Submitted 22 July, 2022; v1 submitted 2 November, 2021; originally announced November 2021.

Validating Synthetic Galaxy Catalogs for Dark Energy Science in the LSST Era

Authors: Eve Kovacs, Yao-Yuan Mao, Michel Aguena, Anita Bahmanyar, Adam Broussard, James Butler, Duncan Campbell, Chihway Chang, Shenming Fu, Katrin Heitmann, Danila Korytov, François Lanusse, Patricia Larsen, Rachel Mandelbaum, Christopher B. Morrison, Constantin Payerne, Marina Ricci, Eli Rykoff, F. Javier Sánchez, Ignacio Sevilla-Noarbe, Melanie Simet, Chun-Hao To, Vinu Vikraman, Rongpu Zhou, Camille Avestruz , et al. (14 additional authors not shown)

Abstract: Large simulation efforts are required to provide synthetic galaxy catalogs for ongoing and upcoming cosmology surveys. These extragalactic catalogs are being used for many diverse purposes covering a wide range of scientific topics. In order to be useful, they must offer realistically complex information about the galaxies they contain. Hence, it is critical to implement a rigorous validation proc… ▽ More Large simulation efforts are required to provide synthetic galaxy catalogs for ongoing and upcoming cosmology surveys. These extragalactic catalogs are being used for many diverse purposes covering a wide range of scientific topics. In order to be useful, they must offer realistically complex information about the galaxies they contain. Hence, it is critical to implement a rigorous validation procedure that ensures that the simulated galaxy properties faithfully capture observations and delivers an assessment of the level of realism attained by the catalog. We present here a suite of validation tests that have been developed by the Rubin Observatory Legacy Survey of Space and Time (LSST) Dark Energy Science Collaboration (DESC). We discuss how the inclusion of each test is driven by the scientific targets for static ground-based dark energy science and by the availability of suitable validation data. The validation criteria that are used to assess the performance of a catalog are flexible and depend on the science goals. We illustrate the utility of this suite by showing examples for the validation of cosmoDC2, the extragalactic catalog recently released for the LSST DESC second Data Challenge. △ Less

Submitted 13 January, 2022; v1 submitted 7 October, 2021; originally announced October 2021.

Comments: 46 pages, 33 figures

Farpoint: A High-Resolution Cosmology Simulation at the Gigaparsec Scale

Authors: Nicholas Frontiere, Katrin Heitmann, Esteban Rangel, Patricia Larsen, Adrian Pope, Imran Sultan, Thomas Uram, Salman Habib, Silvio Rizzi, Joe Insley

Abstract: In this paper we introduce the Farpoint simulation, the latest member of the Hardware/Hybrid Accelerated Cosmology Code (HACC) gravity-only simulation family. The domain covers a volume of (1000$h^{-1}$Mpc)$^3$ and evolves close to two trillion particles, corresponding to a mass resolution of $m_p\sim 4.6\cdot 10^7 h^{-1}$M$_\odot$. These specifications enable comprehensive investigations of the g… ▽ More In this paper we introduce the Farpoint simulation, the latest member of the Hardware/Hybrid Accelerated Cosmology Code (HACC) gravity-only simulation family. The domain covers a volume of (1000$h^{-1}$Mpc)$^3$ and evolves close to two trillion particles, corresponding to a mass resolution of $m_p\sim 4.6\cdot 10^7 h^{-1}$M$_\odot$. These specifications enable comprehensive investigations of the galaxy-halo connection, capturing halos down to small masses. Further, the large volume resolves scales typical of modern surveys with good statistical coverage of high mass halos. The simulation was carried out on the GPU-accelerated system Summit, one of the fastest supercomputers currently available. We provide specifics about the Farpoint run and present an initial set of results. The high mass resolution facilitates precise measurements of important global statistics, such as the halo concentration-mass relation and the correlation function down to small scales. Selected subsets of the simulation data products are publicly available via the HACC Simulation Data Portal. △ Less

Submitted 28 February, 2022; v1 submitted 4 September, 2021; originally announced September 2021.

Journal ref: The Astrophysical Journal Supplement Series, 259(1):15, 2022

A Survey of Practical Formal Methods for Security

Authors: Tomas Kulik, Brijesh Dongol, Peter Gorm Larsen, Hugo Daniel Macedo, Steve Schneider, Peter Würtz Vinther Tran-Jørgensen, Jim Woodcock

Abstract: In today's world, critical infrastructure is often controlled by computing systems. This introduces new risks for cyber attacks, which can compromise the security and disrupt the functionality of these systems. It is therefore necessary to build such systems with strong guarantees of resiliency against cyber attacks. One way to achieve this level of assurance is using formal verification, which pr… ▽ More In today's world, critical infrastructure is often controlled by computing systems. This introduces new risks for cyber attacks, which can compromise the security and disrupt the functionality of these systems. It is therefore necessary to build such systems with strong guarantees of resiliency against cyber attacks. One way to achieve this level of assurance is using formal verification, which provides proofs of system compliance with desired cyber security properties. The use of Formal Methods (FM) in aspects of cyber security and safety-critical systems are reviewed in this article. We split FM into the three main classes: theorem proving, model checking and lightweight FM. To allow the different uses of FM to be compared, we define a common set of terms. We further develop categories based on the type of computing system FM are applied in. Solutions in each class and category are presented, discussed, compared and summarised. We describe historical highlights and developments and present a state-of-the-art review in the area of FM in cyber security. This review is presented from the point of view of FM practitioners and researchers, commenting on the trends in each of the classes and categories. This is achieved by considering all types of FM, several types of security and safety critical systems and by structuring the taxonomy accordingly. The article hence provides a comprehensive overview of FM and techniques available to system designers of security-critical systems, simplifying the process of choosing the right tool for the task. The article concludes by summarising the discussion of the review, focusing on best practices, challenges, general future trends and directions of research within this field. △ Less

Submitted 3 September, 2021; originally announced September 2021.

Comments: Technical Report, Long survey version

The Specification Language Server Protocol: A Proposal for Standardised LSP Extensions

Authors: Jonas Kjær Rask, Frederik Palludan Madsen, Nick Battle, Hugo Daniel Macedo, Peter Gorm Larsen

Abstract: The Language Server Protocol (LSP) changed the field of Integrated Development Environments(IDEs), as it decouples core (programming) language features functionality from editor smarts, thus lowering the effort required to extend an IDE to support a language. The concept is a success and has been adopted by several programming languages and beyond. This is shown by the emergence of several LSP imp… ▽ More The Language Server Protocol (LSP) changed the field of Integrated Development Environments(IDEs), as it decouples core (programming) language features functionality from editor smarts, thus lowering the effort required to extend an IDE to support a language. The concept is a success and has been adopted by several programming languages and beyond. This is shown by the emergence of several LSP implementations for the many programming and specification languages (languages with a focus on modelling, reasoning, or proofs). However, for such languages LSP has been ad-hocly extended with the additional functionalities that are typically not found for programming languages and thus not supported in LSP. This foils the original LSP decoupling goal, because the move towards a new IDE requires yet another re-implementation of the ad-hoc LSP extension. In this paper we contribute with a conservative extension of LSP providing a first proposal towards a standard protocol decoupling the support of specification languages from the IDE. We hope our research attracts the larger community and motivates the need of a joint task force leading to a standardised LSP extension serving the particular needs of specification languages. △ Less

Submitted 6 August, 2021; originally announced August 2021.

Comments: In Proceedings F-IDE 2021, arXiv:2108.02369

Journal ref: EPTCS 338, 2021, pp. 3-18

The Incubator Case Study for Digital Twin Engineering

Authors: Hao Feng, Cláudio Gomes, Casper Thule, Kenneth Lausdahl, Michael Sandberg, Peter Gorm Larsen

Abstract: To demystify the Digital Twin concept, we built a simple yet representative thermal incubator system. The incubator is an insulated box fitted with a heatbed, and complete with a software system for communication, a controller, and simulation models. We developed two simulation models to predict the temperature inside the incubator, one with two free parameters and one with four free parameters. O… ▽ More To demystify the Digital Twin concept, we built a simple yet representative thermal incubator system. The incubator is an insulated box fitted with a heatbed, and complete with a software system for communication, a controller, and simulation models. We developed two simulation models to predict the temperature inside the incubator, one with two free parameters and one with four free parameters. Our experiments showed that the latter model was better at predicting the thermal inertia of the heatbed itself, which makes it more appropriate for further development of the digital twin. The hardware and software used in this case study are available open source, providing an accessible platform for those who want to develop and verify their own techniques for digital twins. △ Less

Submitted 20 February, 2021; originally announced February 2021.

Comments: 18 pages, 13 figures

DESC DC2 Data Release Note

Authors: LSST Dark Energy Science Collaboration, Bela Abolfathi, Robert Armstrong, Humna Awan, Yadu N. Babuji, Franz Erik Bauer, George Beckett, Rahul Biswas, Joanne R. Bogart, Dominique Boutigny, Kyle Chard, James Chiang, Johann Cohen-Tanugi, Andrew J. Connolly, Scott F. Daniel, Seth W. Digel, Alex Drlica-Wagner, Richard Dubois, Eric Gawiser, Thomas Glanzman, Salman Habib, Andrew P. Hearin, Katrin Heitmann, Fabio Hernandez, Renée Hložek , et al. (32 additional authors not shown)

Abstract: In preparation for cosmological analyses of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST), the LSST Dark Energy Science Collaboration (LSST DESC) has created a 300 deg$^2$ simulated survey as part of an effort called Data Challenge 2 (DC2). The DC2 simulated sky survey, in six optical bands with observations following a reference LSST observing cadence, was processed with th… ▽ More In preparation for cosmological analyses of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST), the LSST Dark Energy Science Collaboration (LSST DESC) has created a 300 deg$^2$ simulated survey as part of an effort called Data Challenge 2 (DC2). The DC2 simulated sky survey, in six optical bands with observations following a reference LSST observing cadence, was processed with the LSST Science Pipelines (19.0.0). In this Note, we describe the public data release of the resulting object catalogs for the coadded images of five years of simulated observations along with associated truth catalogs. We include a brief description of the major features of the available data sets. To enable convenient access to the data products, we have developed a web portal connected to Globus data services. We describe how to access the data and provide example Jupyter Notebooks in Python to aid first interactions with the data. We welcome feedback and questions about the data release via a GitHub repository. △ Less

Submitted 13 June, 2022; v1 submitted 12 January, 2021; originally announced January 2021.

Comments: 25 pages, 3 figures; 9 tables. A detailed changelog can be found in Appendix A. To obtain data, visit the DESC Data Portal at https://data.lsstdesc.org/

The Last Journey. II. SMACC -- Subhalo Mass-loss Analysis using Core Catalogs

Authors: Imran Sultan, Nicholas Frontiere, Salman Habib, Katrin Heitmann, Eve Kovacs, Patricia Larsen, Esteban Rangel

Abstract: This paper introduces SMACC -- Subhalo Mass-loss Analysis using Core Catalogs. SMACC adds a mass model to substructure merger trees based on halo "core tracking." Our approach avoids the need for running expensive subhalo finding algorithms and instead uses subhalo mass-loss modeling to assign masses to halo cores. We present details of the SMACC methodology and demonstrate its excellent performan… ▽ More This paper introduces SMACC -- Subhalo Mass-loss Analysis using Core Catalogs. SMACC adds a mass model to substructure merger trees based on halo "core tracking." Our approach avoids the need for running expensive subhalo finding algorithms and instead uses subhalo mass-loss modeling to assign masses to halo cores. We present details of the SMACC methodology and demonstrate its excellent performance in describing halo substructure and its evolution. Validation of the approach is carried out using cosmological simulations at significantly different resolutions. We apply SMACC to the 1.24 trillion-particle Last Journey simulation and construct core catalogs with the additional mass information. These catalogs can be readily used as input to semi-analytic models or subhalo abundance matching approaches to determine approximate galaxy distributions, as well as for in-depth studies of small-scale structure evolution. △ Less

Submitted 16 December, 2020; originally announced December 2020.

Comments: 15 pages, 10 figures. Submitted to ApJ

dMVX: Secure and Efficient Multi-Variant Execution in a Distributed Setting

Authors: Alexios Voulimeneas, Dokyung Song, Per Larsen, Michael Franz, Stijn Volckaert

Abstract: Multi-variant execution (MVX) systems amplify the effectiveness of software diversity techniques. The key idea is to run multiple diversified program variants in lockstep while providing them with the same input and monitoring their run-time behavior for divergences. Thus, adversaries have to compromise all program variants simultaneously to mount an attack successfully. Recent work proposed distr… ▽ More Multi-variant execution (MVX) systems amplify the effectiveness of software diversity techniques. The key idea is to run multiple diversified program variants in lockstep while providing them with the same input and monitoring their run-time behavior for divergences. Thus, adversaries have to compromise all program variants simultaneously to mount an attack successfully. Recent work proposed distributed, heterogeneous MVX systems that leverage different ABIs and ISAs to increase the diversity between program variants further. However, existing distributed MVX system designs suffer from high performance overhead due to time-consuming network transactions for the MVX system's operations. This paper presents dMVX, a novel hybrid distributed MVX design, which incorporates new techniques that significantly reduce the overhead of MVX systems in a distributed setting. Our key insight is that we can intelligently reduce the MVX operations that use expensive network transfers. First, we can limit the monitoring of system calls that are not security-critical. Second, we observe that, in many circumstances, we can also safely cache or avoid replication operations needed for I/O related system calls. Our evaluation shows that dMVX reduces the performance degradation from over 50% to 3.1% for realistic server benchmarks. △ Less

Submitted 3 November, 2020; originally announced November 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.

Comments: 39 pages, 19 figures, version accepted for publication in ApJS

The Last Journey. I. An Extreme-Scale Simulation on the Mira Supercomputer

Authors: Katrin Heitmann, Nicholas Frontiere, Esteban Rangel, Patricia Larsen, Adrian Pope, Imran Sultan, Thomas Uram, Salman Habib, Hal Finkel, Danila Korytov, Eve Kovacs, Silvio Rizzi, Joe Insley

Abstract: The Last Journey is a large-volume, gravity-only, cosmological N-body simulation evolving more than 1.24 trillion particles in a periodic box with a side-length of 5.025Gpc. It was implemented using the HACC simulation and analysis framework on the BG/Q system, Mira. The cosmological parameters are chosen to be consistent with the results from the Planck satellite. A range of analysis tools have b… ▽ More The Last Journey is a large-volume, gravity-only, cosmological N-body simulation evolving more than 1.24 trillion particles in a periodic box with a side-length of 5.025Gpc. It was implemented using the HACC simulation and analysis framework on the BG/Q system, Mira. The cosmological parameters are chosen to be consistent with the results from the Planck satellite. A range of analysis tools have been run in situ to enable a diverse set of science projects, and at the same time, to keep the resulting data amount manageable. Analysis outputs have been generated starting at redshift z~10 to allow for construction of synthetic galaxy catalogs using a semi-analytic modeling approach in post-processing. As part of our in situ analysis pipeline we employ a new method for tracking halo sub-structures, introducing the concept of subhalo cores. The production of multi-wavelength synthetic sky maps is facilitated by generating particle lightcones in situ, also beginning at z~10. We provide an overview of the simulation set-up and the generated data products; a first set of analysis results is presented. A subset of the data is publicly available. △ Less

Submitted 8 January, 2021; v1 submitted 2 June, 2020; originally announced June 2020.

Comments: 14 pages, 9 figures. Accepted for publication in ApJS. New visualization and new results for the matter correlation function added, minor edits. The Last Journey data products can be accessed here: https://cosmology.alcf.anl.gov/

A Cloud-Based Collaboration Platform for Model-Based Design of Cyber-Physical Systems

Authors: Peter Gorm Larsen, Hugo Daniel Macedo, John Fitzgerald, Holger Pfeifer, Martin Benedikt, Stefano Tonetta, Angelo Marguglio, Sergio Gusmeroli, George Suciu Jr

Abstract: Businesses, particularly small and medium-sized enterprises, aiming to start up in Model-Based Design (MBD) face difficult choices from a wide range of methods, notations and tools before making the significant investments in planning, procurement and training necessary to deploy new approaches successfully. In the development of Cyber-Physical Systems (CPSs) this is exacerbated by the diversity o… ▽ More Businesses, particularly small and medium-sized enterprises, aiming to start up in Model-Based Design (MBD) face difficult choices from a wide range of methods, notations and tools before making the significant investments in planning, procurement and training necessary to deploy new approaches successfully. In the development of Cyber-Physical Systems (CPSs) this is exacerbated by the diversity of formalisms covering computation, physical and human processes. In this paper, we propose the use of a cloud-enabled and open collaboration platform that allows businesses to offer models, tools and other assets, and permits others to access these on a pay-per-use basis as a means of lowering barriers to the adoption of MBD technology, and to promote experimentation in a sandbox environment. △ Less

Submitted 5 May, 2020; originally announced May 2020.

Revisiting the Common Neighbour Analysis and the Centrosymmetry Parameter

Authors: Peter M Larsen

Abstract: We review two standard methods for structural classification in simulations of crystalline phases, the Common Neighbour Analysis and the Centrosymmetry Parameter. We explore the definitions and implementations of each of their common variants, and investigate their respective failure modes and classification biases. Simple modifications to both methods are proposed, which improve their robustness,… ▽ More We review two standard methods for structural classification in simulations of crystalline phases, the Common Neighbour Analysis and the Centrosymmetry Parameter. We explore the definitions and implementations of each of their common variants, and investigate their respective failure modes and classification biases. Simple modifications to both methods are proposed, which improve their robustness, interpretability, and applicability. We denote these variants the Interval Common Neighbour Analysis, and the Minimum-Weight Matching Centrosymmetry Parameter. △ Less

Submitted 19 March, 2020; originally announced March 2020.

Comments: 11 pages, 8 figures

Resolving pseudosymmetry in tetragonal ZrO2 using EBSD with a modified dictionary indexing approach

Authors: Edward L. Pang, Peter M. Larsen, Christopher A. Schuh

Abstract: Resolving pseudosymmetry has long presented a challenge for electron backscatter diffraction (EBSD) and has been notoriously challenging in the case of tetragonal ZrO2 in particular. In this work, a method is proposed to resolve pseudosymmetry by building upon the dictionary indexing method and augmenting it with the application of global optimization to fit accurate pattern centers, clustering of… ▽ More Resolving pseudosymmetry has long presented a challenge for electron backscatter diffraction (EBSD) and has been notoriously challenging in the case of tetragonal ZrO2 in particular. In this work, a method is proposed to resolve pseudosymmetry by building upon the dictionary indexing method and augmenting it with the application of global optimization to fit accurate pattern centers, clustering of the Hough-indexed orientations to focus the dictionary in orientation space, and interpolation to improve the accuracy of the indexed solution. The proposed method is demonstrated to resolve pseudosymmetry with 100% accuracy in simulated patterns of tetragonal ZrO2, even with high degrees of binning and noise. The method is then used to index an experimental dataset, which confirms its ability to efficiently and accurately resolve pseudosymmetry in these materials. The present method can be applied to resolve pseudosymmetry in a wide range of materials, possibly even some more challenging than tetragonal ZrO2. Source code for this implementation is available online. △ Less

Submitted 9 March, 2020; originally announced March 2020.

Minimum-Strain Symmetrization of Bravais Lattices

Authors: Peter M. Larsen, Edward L. Pang, Pablo A. Parrilo, Karsten W. Jacobsen

Abstract: Bravais lattices are the most fundamental building blocks of crystallography. They are classified into groups according to their translational, rotational, and inversion symmetries. In computational analysis of Bravais lattices, fulfilment of symmetry conditions is usually determined by analysis of the metric tensor, using either a numerical tolerance to produce a binary (i.e. yes or no) classific… ▽ More Bravais lattices are the most fundamental building blocks of crystallography. They are classified into groups according to their translational, rotational, and inversion symmetries. In computational analysis of Bravais lattices, fulfilment of symmetry conditions is usually determined by analysis of the metric tensor, using either a numerical tolerance to produce a binary (i.e. yes or no) classification, or a distance function which quantifies the deviation from an ideal lattice type. The metric tensor, though, is not scale-invariant, which complicates the choice of threshold and the interpretation of the distance function. Here, we quantify the distance of a lattice from a target Bravais class using strain. For an arbitrary lattice, we find the minimum-strain transformation needed to fulfil the symmetry conditions of a desired Bravais lattice type; the norm of the strain tensor is used to quantify the degree of symmetry breaking. The resulting distance is invariant to scale and rotation, and is a physically intuitive quantity. By symmetrizing to all Bravais classes, each lattice can be placed in a 14 dimensional space, which we use to create a map of the space of Bravais lattices and the transformation paths between them. A software implementation is available online under a permissive license. △ Less

Submitted 8 October, 2019; originally announced October 2019.

Comments: 12 pages, 9 figures

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

Global optimization for accurate determination of EBSD pattern centers

Authors: Edward L. Pang, Peter M. Larsen, Christopher A. Schuh

Abstract: Accurate pattern center determination has long been a challenge for the electron backscatter diffraction (EBSD) community and is becoming critically accuracy-limiting for more recent advanced EBSD techniques. Here, we study the parameter landscape over which a pattern center must be fitted in quantitative detail and reveal that it is both sloppy and noisy, which limits the accuracy to which patter… ▽ More Accurate pattern center determination has long been a challenge for the electron backscatter diffraction (EBSD) community and is becoming critically accuracy-limiting for more recent advanced EBSD techniques. Here, we study the parameter landscape over which a pattern center must be fitted in quantitative detail and reveal that it is both sloppy and noisy, which limits the accuracy to which pattern centers can be determined. To locate the global optimum in this challenging landscape, we propose a combination of two approaches: the use of a global search algorithm and averaging the results from multiple patterns. We demonstrate the ability to accurately determine pattern centers of simulated patterns, inclusive of effects of binning and noise on the error of the fitted pattern center. We also demonstrate the ability of this method to accurately detect changes in pattern center in an experimental dataset with noisy and highly binned patterns. Source code for our pattern center fitting algorithm is available online. △ Less

Submitted 28 August, 2019; originally announced August 2019.

CosmoDC2: A Synthetic Sky Catalog for Dark Energy Science with LSST

Authors: Danila Korytov, Andrew Hearin, Eve Kovacs, Patricia Larsen, Esteban Rangel, Joseph Hollowed, Andrew J. Benson, Katrin Heitmann, Yao-Yuan Mao, Anita Bahmanyar, Chihway Chang, Duncan Campbell, Joseph Derose, Hal Finkel, Nicholas Frontiere, Eric Gawiser, Salman Habib, Benjamin Joachimi, François Lanusse, Nan Li, Rachel Mandelbaum, Christopher Morrison, Jeffrey A. Newman, Adrian Pope, Eli Rykoff , et al. (5 additional authors not shown)

Abstract: This paper introduces cosmoDC2, a large synthetic galaxy catalog designed to support precision dark energy science with the Large Synoptic Survey Telescope (LSST). CosmoDC2 is the starting point for the second data challenge (DC2) carried out by the LSST Dark Energy Science Collaboration (LSST DESC). The catalog is based on a trillion-particle, 4.225 Gpc^3 box cosmological N-body simulation, the `… ▽ More This paper introduces cosmoDC2, a large synthetic galaxy catalog designed to support precision dark energy science with the Large Synoptic Survey Telescope (LSST). CosmoDC2 is the starting point for the second data challenge (DC2) carried out by the LSST Dark Energy Science Collaboration (LSST DESC). The catalog is based on a trillion-particle, 4.225 Gpc^3 box cosmological N-body simulation, the `Outer Rim' run. It covers 440 deg^2 of sky area to a redshift of z=3 and is complete to a magnitude depth of 28 in the r-band. Each galaxy is characterized by a multitude of properties including stellar mass, morphology, spectral energy distributions, broadband filter magnitudes, host halo information and weak lensing shear. The size and complexity of cosmoDC2 requires an efficient catalog generation methodology; our approach is based on a new hybrid technique that combines data-driven empirical approaches with semi-analytic galaxy modeling. A wide range of observation-based validation tests has been implemented to ensure that cosmoDC2 enables the science goals of the planned LSST DESC DC2 analyses. This paper also represents the official release of the cosmoDC2 data set, including an efficient reader that facilitates interaction with the data. △ Less

Submitted 27 July, 2019; v1 submitted 15 July, 2019; originally announced July 2019.

Comments: 27 pages, 17 figures, submitted to APJS

HACC Cosmological Simulations: First Data Release

Authors: Katrin Heitmann, Thomas D. Uram, Hal Finkel, Nicholas Frontiere, Salman Habib, Adrian Pope, Esteban Rangel, Joseph Hollowed, Danila Korytov, Patricia Larsen, Benjamin S. Allen, Kyle Chard, Ian Foster

Abstract: We describe the first major public data release from cosmological simulations carried out with Argonne's HACC code. This initial release covers a range of datasets from large gravity-only simulations. The data products include halo information for multiple redshifts, down-sampled particles, and lightcone outputs. We provide data from two very large LCDM simulations as well as beyond-LCDM simulatio… ▽ More We describe the first major public data release from cosmological simulations carried out with Argonne's HACC code. This initial release covers a range of datasets from large gravity-only simulations. The data products include halo information for multiple redshifts, down-sampled particles, and lightcone outputs. We provide data from two very large LCDM simulations as well as beyond-LCDM simulations spanning eleven w0-wa cosmologies. Our release platform uses Petrel, a research data service, located at the Argonne Leadership Computing Facility. Petrel offers fast data transfer mechanisms and authentication via Globus, enabling simple and efficient access to stored datasets. Easy browsing of the available data products is provided via a web portal that allows the user to navigate simulation products efficiently. The data hub will be extended by adding more types of data products and by enabling computational capabilities to allow direct interactions with simulation results. △ Less

Submitted 3 October, 2019; v1 submitted 26 April, 2019; originally announced April 2019.

Comments: 8 pages, 5 figures. Final version published in ApJS. The HACC Simulation Data Portal can be accessed here: https://cosmology.alcf.anl.gov/

Journal ref: Astrophys.J.Suppl. 244 (2019) no.1, 17

DMON: A Distributed Heterogeneous N-Variant System

Authors: Alexios Voulimeneas, Dokyung Song, Fabian Parzefall, Yeoul Na, Per Larsen, Michael Franz, Stijn Volckaert

Abstract: N-Variant Execution (NVX) systems utilize software diversity techniques for enhancing software security. The general idea is to run multiple different variants of the same program alongside each other while monitoring their run-time behavior. If the internal disparity between the running variants causes observable differences in response to malicious inputs, the monitor can detect such divergences… ▽ More N-Variant Execution (NVX) systems utilize software diversity techniques for enhancing software security. The general idea is to run multiple different variants of the same program alongside each other while monitoring their run-time behavior. If the internal disparity between the running variants causes observable differences in response to malicious inputs, the monitor can detect such divergences in execution and then raise an alert and/or terminate execution. Existing NVX systems execute multiple, artificially diversified program variants on a single host. This paper presents a novel, distributed NVX design that executes program variants across multiple heterogeneous host computers; our prototype implementation combines an x86-64 host with an ARMv8 host. Our approach greatly increases the level of "internal different-ness" between the simultaneously running variants that can be supported, encompassing different instruction sets, endianness, calling conventions, system call interfaces, and potentially also differences in hardware security features. A major challenge to building such a heterogeneous distributed NVX system is performance. We present solutions to some of the main performance challenges. We evaluate our prototype system implementing these ideas to show that it can provide reasonable performance on a wide range of realistic workloads. △ Less

Submitted 8 March, 2019; originally announced March 2019.

Dark Energy Survey Year 1 Results: Constraints on Intrinsic Alignments and their Colour Dependence from Galaxy Clustering and Weak Lensing

Authors: S. Samuroff, J. Blazek, M. A. Troxel, N. MacCrann, E. Krause, C. D. Leonard, J. Prat, D. Gruen, S. Dodelson, T. F. Eifler, M. Gatti, W. G. Hartley, B. Hoyle, P. Larsen, J. Zuntz, T. M. C. Abbott, S. Allam, J. Annis, G. M. Bernstein, E. Bertin, S. L. Bridle, D. Brooks, A. Carnero Rosell, M. Carrasco Kind, J. Carretero , et al. (48 additional authors not shown)

Abstract: We perform a joint analysis of intrinsic alignments and cosmology using tomographic weak lensing, galaxy clustering and galaxy-galaxy lensing measurements from Year 1 (Y1) of the Dark Energy Survey. We define early- and late-type subsamples, which are found to pass a series of systematics tests, including for spurious photometric redshift error and point spread function correlations. We analyse th… ▽ More We perform a joint analysis of intrinsic alignments and cosmology using tomographic weak lensing, galaxy clustering and galaxy-galaxy lensing measurements from Year 1 (Y1) of the Dark Energy Survey. We define early- and late-type subsamples, which are found to pass a series of systematics tests, including for spurious photometric redshift error and point spread function correlations. We analyse these split data alongside the fiducial mixed Y1 sample using a range of intrinsic alignment models. In a fiducial Nonlinear Alignment Model (NLA) analysis, assuming a flat \lcdm~cosmology, we find a significant difference in intrinsic alignment amplitude, with early-type galaxies favouring $A_\mathrm{IA} = 2.38^{+0.32}_{-0.31}$ and late-type galaxies consistent with no intrinsic alignments at $0.05^{+0.10}_{-0.09}$. We find weak evidence of a diminishing alignment amplitude at higher redshifts in the early-type sample. The analysis is repeated using a number of extended model spaces, including a physically motivated model that includes both tidal torquing and tidal alignment mechanisms. In multiprobe likelihood chains in which cosmology, intrinsic alignments in both galaxy samples and all other relevant systematics are varied simultaneously, we find the tidal alignment and tidal torquing parts of the intrinsic alignment signal have amplitudes $A_1 = 2.66 ^{+0.67}_{-0.66}$, $A_2=-2.94^{+1.94}_{-1.83}$, respectively, for early-type galaxies and $A_1 = 0.62 ^{+0.41}_{-0.41}$, $A_2 = -2.26^{+1.30}_{-1.16}$ for late-type galaxies. In the full (mixed) Y1 sample the best constraints are $A_1 = 0.70 ^{+0.41}_{-0.38}$, $A_2 = -1.36 ^{+1.08}_{-1.41}$. For all galaxy splits and IA models considered, we report cosmological parameter constraints that are consistent with the results of Troxel et al. (2017) and Dark Energy Survey Collaboration (2017). △ Less

Submitted 6 August, 2019; v1 submitted 16 November, 2018; originally announced November 2018.

Comments: 31 pages, 23 figures; accepted by MNRAS

Report number: FERMILAB-PUB-18-622-AE

The Borg Cube Simulation: Cosmological Hydrodynamics with CRK-SPH

Authors: J. D. Emberson, Nicholas Frontiere, Salman Habib, Katrin Heitmann, Patricia Larsen, Hal Finkel, Adrian Pope

Abstract: A challenging requirement posed by next-generation observations is a firm theoretical grasp of the impact of baryons on structure formation. Cosmological hydrodynamic simulations modeling gas physics are vital in this regard. A high degree of modeling flexibility exists in this space making it important to explore a range of methods in order to gauge the accuracy of simulation predictions. We pres… ▽ More A challenging requirement posed by next-generation observations is a firm theoretical grasp of the impact of baryons on structure formation. Cosmological hydrodynamic simulations modeling gas physics are vital in this regard. A high degree of modeling flexibility exists in this space making it important to explore a range of methods in order to gauge the accuracy of simulation predictions. We present results from the first cosmological simulation using Conservative Reproducing Kernel Smoothed Particle Hydrodynamics (CRK-SPH). We employ two simulations: one evolved purely under gravity and the other with non-radiative hydrodynamics. Each contains 2x2304^3 cold dark matter plus baryon particles in an 800 Mpc/h box. We compare statistics to previous non-radiative simulations including power spectra, mass functions, baryon fractions, and concentration. We find self-similar radial profiles of gas temperature, entropy, and pressure and show that a simple analytic model recovers these results to better than 40% over two orders of magnitude in mass. We quantify the level of non-thermal pressure support in halos and demonstrate that hydrostatic mass estimates are biased low by 24% (10%) for halos of mass 10^15 (10^13) Msun/h. We compute angular power spectra for the thermal and kinematic Sunyaev-Zel'dovich effects and find good agreement with the low-l Planck measurements. Finally, artificial scattering between particles of unequal mass is shown to have a large impact on the gravity-only run and we highlight the importance of better understanding this issue in hydrodynamic applications. This is the first in a simulation campaign using CRK-SPH with future work including subresolution gas treatments. △ Less

Submitted 3 June, 2019; v1 submitted 8 November, 2018; originally announced November 2018.

Comments: Accepted for publication in the Astrophysical Journal

Journal ref: ApJ, 877, 85, 2019

Co-simulation of Continuous Systems: A Tutorial

Authors: Cláudio Gomes, Casper Thule, Peter Gorm Larsen, Joachim Denil, Hans Vangheluwe

Abstract: Co-simulation consists of the theory and techniques to enable global simulation of a coupled system via the composition of simulators. Despite the large number of applications and growing interest in the challenges, the field remains fragmented into multiple application domains, with limited sharing of knowledge. This tutorial aims at introducing co-simulation of continuous systems, targeted at… ▽ More Co-simulation consists of the theory and techniques to enable global simulation of a coupled system via the composition of simulators. Despite the large number of applications and growing interest in the challenges, the field remains fragmented into multiple application domains, with limited sharing of knowledge. This tutorial aims at introducing co-simulation of continuous systems, targeted at researchers new to the field. △ Less

Submitted 22 September, 2018; originally announced September 2018.

MSC Class: 65Y10 ACM Class: I.6.1; I.6.7

Definition of a scoring parameter to identify low-dimensional materials components

Authors: Peter Mahler Larsen, Mohnish Pandey, Mikkel Strange, Karsten Wedel Jacobsen

Abstract: The last decade has seen intense research in materials with reduced dimensionality. The low dimensionality leads to interesting electronic behavior due to electronic confinement and reduced screening. The investigations have to a large extent focused on 2D materials both in their bulk form, as individual layers a few atoms thick, and through stacking of 2D layers into heterostructures. The identif… ▽ More The last decade has seen intense research in materials with reduced dimensionality. The low dimensionality leads to interesting electronic behavior due to electronic confinement and reduced screening. The investigations have to a large extent focused on 2D materials both in their bulk form, as individual layers a few atoms thick, and through stacking of 2D layers into heterostructures. The identification of low-dimensional compounds is therefore of key interest. Here, we perform a geometric analysis of material structures, demonstrating a strong clustering of materials depending on their dimensionalities. Based on the geometric analysis, we propose a simple scoring parameter to identify materials of a particular dimension or of mixed dimensionality. The method identifies spatially connected components of the materials and gives a measure of the degree of "1D-ness," "2D-ness," etc., for each component. The scoring parameter is applied to the Inorganic Crystal Structure Database and the Crystallography Open Database ranking the materials according to their degree of dimensionality. In the case of 2D materials the scoring parameter is seen to clearly separate 2D from non-2D materials and the parameter correlates well with the bonding strength in the layered materials. About 3000 materials are identified as one-dimensional, while more than 9000 are mixed-dimensionality materials containing a molecular (0D) component. The charge states of the components in selected highly ranked materials are investigated using density functional theory and Bader analysis showing that the spatially separated components have either zero charge, corresponding to weak interactions, or integer charge, indicating ionic bonding. △ Less

Submitted 19 March, 2019; v1 submitted 6 August, 2018; originally announced August 2018.

Journal ref: Phys. Rev. Materials 3, 034003 (2019)

SoK: Sanitizing for Security

Authors: Dokyung Song, Julian Lettner, Prabhu Rajasekaran, Yeoul Na, Stijn Volckaert, Per Larsen, Michael Franz

Abstract: The C and C++ programming languages are notoriously insecure yet remain indispensable. Developers therefore resort to a multi-pronged approach to find security issues before adversaries. These include manual, static, and dynamic program analysis. Dynamic bug finding tools --- henceforth "sanitizers" --- can find bugs that elude other types of analysis because they observe the actual execution of a… ▽ More The C and C++ programming languages are notoriously insecure yet remain indispensable. Developers therefore resort to a multi-pronged approach to find security issues before adversaries. These include manual, static, and dynamic program analysis. Dynamic bug finding tools --- henceforth "sanitizers" --- can find bugs that elude other types of analysis because they observe the actual execution of a program, and can therefore directly observe incorrect program behavior as it happens. A vast number of sanitizers have been prototyped by academics and refined by practitioners. We provide a systematic overview of sanitizers with an emphasis on their role in finding security issues. Specifically, we taxonomize the available tools and the security vulnerabilities they cover, describe their performance and compatibility properties, and highlight various trade-offs. △ Less

Submitted 12 June, 2018; originally announced June 2018.

The Computational 2D Materials Database: High-Throughput Modeling and Discovery of Atomically Thin Crystals

Authors: Sten Haastrup, Mikkel Strange, Mohnish Pandey, Thorsten Deilmann, Per S. Schmidt, Nicki F. Hinsche, Morten N. Gjerding, Daniele Torelli, Peter M. Larsen, Anders C. Riis-Jensen, Jakob Gath, Karsten W. Jacobsen, Jens Jørgen Mortensen, Thomas Olsen, Kristian S. Thygesen

Abstract: We introduce the Computational 2D Materials Database (C2DB), which organises a variety of structural, thermodynamic, elastic, electronic, magnetic, and optical properties of around 1500 two-dimensional materials distributed over more than 30 different crystal structures. Material properties are systematically calculated by state-of-the art density functional theory and many-body perturbation theor… ▽ More We introduce the Computational 2D Materials Database (C2DB), which organises a variety of structural, thermodynamic, elastic, electronic, magnetic, and optical properties of around 1500 two-dimensional materials distributed over more than 30 different crystal structures. Material properties are systematically calculated by state-of-the art density functional theory and many-body perturbation theory (G$_0\!$W$\!_0$ and the Bethe-Salpeter Equation for $\sim$200 materials) following a semi-automated workflow for maximal consistency and transparency. The C2DB is fully open and can be browsed online or downloaded in its entirety. In this paper, we describe the workflow behind the database, present an overview of the properties and materials currently available, and explore trends and correlations in the data. Moreover, we identify a large number of new potentially synthesisable 2D materials with interesting properties targeting applications within spintronics, (opto-)electronics, and plasmonics. The C2DB offers a comprehensive and easily accessible overview of the rapidly expanding family of 2D materials and forms an ideal platform for computational modeling and design of new 2D materials and van der Waals heterostructures. △ Less

Submitted 9 October, 2018; v1 submitted 8 June, 2018; originally announced June 2018.

Comments: Add journal reference and DOI; Minor updates to figures and wording

Journal ref: 2d Materials 2018, Volume 5, Issue 4, pp. 042002

Robotic design choice overview using co-simulation

Authors: Martin Peter Christiansen, Peter Gorm Larsen, Rasmus Nyholm Jørgensen

Abstract: Rapid robotic system development sets a demand for multi-disciplinary methods and tools to explore and compare design alternatives. In this paper, we present collaborative modeling that combines discrete-event models of controller software with continuous-time models of physical robot components. The presented co-modeling method utilized VDM for discrete-event and 20-sim for continuous-time modeli… ▽ More Rapid robotic system development sets a demand for multi-disciplinary methods and tools to explore and compare design alternatives. In this paper, we present collaborative modeling that combines discrete-event models of controller software with continuous-time models of physical robot components. The presented co-modeling method utilized VDM for discrete-event and 20-sim for continuous-time modeling. The collaborative modeling method is illustrated with a concrete example of collaborative model development of a mobile robot animal feeding system. Simulations are used to evaluate the robot model output response in relation to operational demands. The result of the simulations provides the developers with an overview of the impacts of each solution instance in the chosen design space. Based on the solution overview the developers can select candidates that are deemed viable to be deployed and tested on an actual physical robot. △ Less

Submitted 17 February, 2018; originally announced February 2018.

Comments: 5 pages, 4 figures, conference

Report number: Agromek and NJF joint seminar (2014), 41-45

Rich Ground State Chemical Ordering in Nanoparticles: Exact Solution of a Model for Ag-Au Clusters

Authors: Peter Mahler Larsen, Karsten Wedel Jacobsen, Jakob Schiøtz

Abstract: We show that nanoparticles can have very rich ground state chemical order. This is illustrated by determining the chemical ordering of Ag-Au 309-atom Mackay icosahedral nanoparticles. The energy of the nanoparticles is described using a cluster expansion model, and a Mixed Integer Programming (MIP) approach is used to find the exact ground state configurations for all stoichiometries. The chemical… ▽ More We show that nanoparticles can have very rich ground state chemical order. This is illustrated by determining the chemical ordering of Ag-Au 309-atom Mackay icosahedral nanoparticles. The energy of the nanoparticles is described using a cluster expansion model, and a Mixed Integer Programming (MIP) approach is used to find the exact ground state configurations for all stoichiometries. The chemical ordering varies widely between the different stoichiometries, and display a rich zoo of structures with non-trivial ordering. △ Less

Submitted 13 June, 2018; v1 submitted 20 December, 2017; originally announced December 2017.

Comments: Revised version. New figure added, discussion expanded, some material moved into supplementary file

Journal ref: Phys. Rev. Lett. 120, 256101 (2018)

PartiSan: Fast and Flexible Sanitization via Run-time Partitioning

Authors: Julian Lettner, Dokyung Song, Taemin Park, Stijn Volckaert, Per Larsen, Michael Franz

Abstract: Sanitizers can detect security vulnerabilities in C/C++ code that elude static analysis. Current practice is to continuously fuzz and sanitize internal pre-release builds. Sanitization-enabled builds are rarely released publicly. This is in large part due to the high memory and processing requirements of sanitizers. We present PartiSan, a run-time partitioning technique that speeds up sanitizers… ▽ More Sanitizers can detect security vulnerabilities in C/C++ code that elude static analysis. Current practice is to continuously fuzz and sanitize internal pre-release builds. Sanitization-enabled builds are rarely released publicly. This is in large part due to the high memory and processing requirements of sanitizers. We present PartiSan, a run-time partitioning technique that speeds up sanitizers and allows them to be used in a more flexible manner. Our core idea is to partition the execution into sanitized slices that incur a run-time overhead, and unsanitized slices running at full speed. With PartiSan, sanitization is no longer an all-or-nothing proposition. A single build can be distributed to every user regardless of their willingness to enable sanitization and the capabilities of their host system. PartiSan can automatically adjust the amount of sanitization to fit within a performance budget or disable sanitization if the host lacks sufficient resources. The flexibility afforded by run-time partitioning also means that we can alternate between different types of sanitizers dynamically; today, developers have to pick a single type of sanitizer ahead of time. Finally, we show that run-time partitioning can speed up fuzzing by running the sanitized partition only when the fuzzer discovers an input that causes a crash or uncovers new execution paths. △ Less

Submitted 14 May, 2018; v1 submitted 21 November, 2017; originally announced November 2017.

Improved Orientation Sampling for Indexing Diffraction Patterns of Polycrystalline Materials

Authors: Peter Mahler Larsen, Søren Schmidt

Abstract: Orientation mapping is a widely used technique for revealing the microstructure of a polycrystalline sample. The crystalline orientation at each point in the sample is determined by analysis of the diffraction pattern, a process known as pattern indexing. A recent development in pattern indexing is the use of a brute-force approach, whereby diffraction patterns are simulated for a large number of… ▽ More Orientation mapping is a widely used technique for revealing the microstructure of a polycrystalline sample. The crystalline orientation at each point in the sample is determined by analysis of the diffraction pattern, a process known as pattern indexing. A recent development in pattern indexing is the use of a brute-force approach, whereby diffraction patterns are simulated for a large number of crystalline orientations, and compared against the experimentally observed diffraction pattern in order to determine the most likely orientation. Whilst this method can robust identify orientations in the presence of noise, it has very high computational requirements. In this article, the computational burden is reduced by developing a method for nearly-optimal sampling of orientations. By using the quaternion representation of orientations, it is shown that the optimal sampling problem is equivalent to that of optimally distributing points on a four-dimensional sphere. In doing so, the number of orientation samples needed to achieve a indexing desired accuracy is significantly reduced. Orientation sets at a range of sizes are generated in this way for all Laue groups, and are made available online for easy use. △ Less

Submitted 10 October, 2017; v1 submitted 25 July, 2017; originally announced July 2017.

Comments: 11 pages, 7 figures

Time-of-Flight Three Dimensional Neutron Diffraction in Transmission Mode for Mapping Crystal Grain Structures

Authors: Alberto Cereser, Markus Strobl, Stephen Hall, Axel Steuwer, Ryoji Kiyanagi, Anton Tremsin, Erik Bergbäck Knudsen, Takenao Shinohara, Peter Willendrup, Alice Bastos da Silva Fanta, Srinivasan Iyengar, Peter Mahler Larsen, Takayasu Hanashima, Taketo Moyoshi, Peter M. Kadletz, Philip Krooß, Thomas Niendorf, Morten Sales, Wolfgang W. Schmahl, Søren Schmidt

Abstract: The physical properties of polycrystalline materials depend on their microstructure, which is the nano-to-centimeter-scale arrangement of phases and defects in their interior. Such microstructure depends on the shape, crystallographic phase and orientation, and interfacing of the grains constituting the material. This article presents a new non-destructive 3D technique to study bulk samples with s… ▽ More The physical properties of polycrystalline materials depend on their microstructure, which is the nano-to-centimeter-scale arrangement of phases and defects in their interior. Such microstructure depends on the shape, crystallographic phase and orientation, and interfacing of the grains constituting the material. This article presents a new non-destructive 3D technique to study bulk samples with sizes in the cm range with a resolution of hundred micrometers: time-of-flight three-dimensional neutron diffraction (ToF 3DND). Compared to existing analogous X-ray diffraction techniques, ToF 3DND enables studies of samples that can be both larger in size and made of heavier elements. Moreover, ToF 3DND facilitates the use of complicated sample environments. The basic ToF 3DND setup, utilizing an imaging detector with high spatial and temporal resolution, can easily be implemented at a time-of-flight neutron beamline. The technique was developed and tested with data collected at the Materials and Life Science Experimental Facility of the Japan Proton Accelerator Complex (J-PARC) for an iron sample. We successfully reconstructed the shape of 108 grains and developed an indexing procedure. The reconstruction algorithms have been validated by reconstructing two stacked Co-Ni-Ga single crystals and by comparison with a grain map obtained by post-mortem electron backscatter diffraction (EBSD). △ Less

Submitted 20 April, 2017; originally announced April 2017.

Comments: Submitted

Journal ref: Nature Scientific Reports 2017

Co-simulation: State of the art

Authors: Cláudio Gomes, Casper Thule, David Broman, Peter Gorm Larsen, Hans Vangheluwe

Abstract: It is essential to find new ways of enabling experts in different disciplines to collaborate more efficient in the development of ever more complex systems, under increasing market pressures. One possible solution for this challenge is to use a heterogeneous model-based approach where different teams can produce their conventional models and carry out their usual mono-disciplinary analysis, but in… ▽ More It is essential to find new ways of enabling experts in different disciplines to collaborate more efficient in the development of ever more complex systems, under increasing market pressures. One possible solution for this challenge is to use a heterogeneous model-based approach where different teams can produce their conventional models and carry out their usual mono-disciplinary analysis, but in addition, the different models can be coupled for simulation (co-simulation), allowing the study of the global behavior of the system. Due to its potential, co-simulation is being studied in many different disciplines but with limited sharing of findings. Our aim with this work is to summarize, bridge, and enhance future research in this multidisciplinary area. We provide an overview of co-simulation approaches, research challenges, and research opportunities, together with a detailed taxonomy with different aspects of the state of the art of co-simulation and classification for the past five years. The main research needs identified are: finding generic approaches for modular, stable and accurate coupling of simulation units; and expressing the adaptations required to ensure that the coupling is correct. △ Less

Submitted 1 February, 2017; originally announced February 2017.

Comments: 157 pages, about 30 figures

MSC Class: 65Y10 ACM Class: I.6.1; I.6.7

Predicting the Plant Root-Associated Ecological Niche of 21 Pseudomonas Species Using Machine Learning and Metabolic Modeling

Authors: Jennifer Chien, Peter Larsen

Abstract: Plants rarely occur in isolated systems. Bacteria can inhabit either the endosphere, the region inside the plant root, or the rhizosphere, the soil region just outside the plant root. Our goal is to understand if using genomic data and media dependent metabolic model information is better for training machine learning of predicting bacterial ecological niche than media independent models or pure g… ▽ More Plants rarely occur in isolated systems. Bacteria can inhabit either the endosphere, the region inside the plant root, or the rhizosphere, the soil region just outside the plant root. Our goal is to understand if using genomic data and media dependent metabolic model information is better for training machine learning of predicting bacterial ecological niche than media independent models or pure genome based species trees. We considered three machine learning techniques: support vector machine, non-negative matrix factorization, and artificial neural networks. In all three machine-learning approaches, the media-based metabolic models and flux balance analyses were more effective at predicting bacterial niche than the genome or PRMT models. Support Vector Machine trained on a minimal media base with Mannose, Proline and Valine was most predictive of all models and media types with an f-score of 0.8 for rhizosphere and 0.97 for endosphere. Thus we can conclude that media-based metabolic modeling provides a holistic view of the metabolome, allowing machine learning algorithms to highlight the differences between and categorize endosphere and rhizosphere bacteria. There was no single media type that best highlighted differences between endosphere and rhizosphere bacteria metabolism and therefore no single enzyme, reaction, or compound that defined whether a bacteria's origin was of the endosphere or rhizosphere. △ Less

Submitted 11 January, 2017; originally announced January 2017.

Comments: 15 pages, keywords: Pseudomonas, SVM, ANN, NMF, FBA, endosphere, rhizosphere, metabolic model, machine learning, KBase Comments: (e.g.: 10 pages, 5 figures, conference or other essential info)

Multi-Variant Execution of Parallel Programs

Authors: Stijn Volckaert, Bjorn De Sutter, Koen De Bosschere, Per Larsen

Abstract: Multi-Variant Execution Environments (MVEEs) are a promising technique to protect software against memory corruption attacks. They transparently execute multiple, diversified variants (often referred to as replicae) of the software receiving the same inputs. By enforcing and monitoring the lock-step execution of the replicae's system calls, and by deploying diversity techniques that prevent an att… ▽ More Multi-Variant Execution Environments (MVEEs) are a promising technique to protect software against memory corruption attacks. They transparently execute multiple, diversified variants (often referred to as replicae) of the software receiving the same inputs. By enforcing and monitoring the lock-step execution of the replicae's system calls, and by deploying diversity techniques that prevent an attacker from simultaneously compromising multiple replicae, MVEEs can block attacks before they succeed. Existing MVEEs cannot handle non-trivial multi-threaded programs because their undeterministic behavior introduces benign system call inconsistencies in the replicae, which trigger false positive detections and deadlocks in the MVEEs. This paper for the first time extends the generality of MVEEs to protect multi-threaded software by means of secure and efficient synchronization replication agents. On the PARSEC 2.1 parallel benchmarks running with four worker threads, our prototype MVEE incurs a run-time overhead of only 1.32x. △ Less

Submitted 26 July, 2016; originally announced July 2016.

Demonstration of cosmic microwave background delensing using the cosmic infrared background

Authors: Patricia Larsen, Anthony Challinor, Blake D. Sherwin, Daisy Mak

Abstract: Delensing is an increasingly important technique to reverse the gravitational lensing of the cosmic microwave background (CMB) and thus reveal primordial signals the lensing may obscure. We present a first demonstration of delensing on Planck temperature maps using the cosmic infrared background (CIB). Reversing the lensing deflections in Planck CMB temperature maps using a linear combination of t… ▽ More Delensing is an increasingly important technique to reverse the gravitational lensing of the cosmic microwave background (CMB) and thus reveal primordial signals the lensing may obscure. We present a first demonstration of delensing on Planck temperature maps using the cosmic infrared background (CIB). Reversing the lensing deflections in Planck CMB temperature maps using a linear combination of the 545 and 857GHz maps as a lensing tracer, we find that the lensing effects in the temperature power spectrum are reduced in a manner consistent with theoretical expectations. In particular, the characteristic sharpening of the acoustic peaks of the temperature power spectrum resulting from successful delensing is detected at a significance of 16$\rmσ$, with an amplitude of $A_{\rm{delens}} = 1.12 \pm 0.07$ relative to the expected value of unity. This first demonstration on data of CIB delensing, and of delensing techniques in general, is significant because lensing removal will soon be essential for achieving high-precision constraints on inflationary B-mode polarization. △ Less

Submitted 6 October, 2016; v1 submitted 19 July, 2016; originally announced July 2016.

Comments: 6 pages, 5 figures, accepted for publication in PRL

Journal ref: Phys. Rev. Lett. 117, 151102 (2016)

Robust Structural Identification via Polyhedral Template Matching

Authors: Peter Mahler Larsen, Søren Schmidt, Jakob Schiøtz

Abstract: Successful scientific applications of large-scale molecular dynamics often rely on automated methods for identifying the local crystalline structure of condensed phases. Many existing methods for structural identification, such as Common Neighbour Analysis, rely on interatomic distances (or thresholds thereof) to classify atomic structure. As a consequence they are sensitive to strain and thermal… ▽ More Successful scientific applications of large-scale molecular dynamics often rely on automated methods for identifying the local crystalline structure of condensed phases. Many existing methods for structural identification, such as Common Neighbour Analysis, rely on interatomic distances (or thresholds thereof) to classify atomic structure. As a consequence they are sensitive to strain and thermal displacements, and preprocessing such as quenching or temporal averaging of the atomic positions is necessary to provide reliable identifications. We propose a new method, Polyhedral Template Matching (PTM), which classifies structures according to the topology of the local atomic environment, without any ambiguity in the classification, and with greater reliability than e.g. Common Neighbour Analysis in the presence of thermal fluctuations. We demonstrate that the method can reliably be used to identify structures even in simulations near the melting point, and that it can identify the most common ordered alloy structures as well. In addition, the method makes it easy to identify the local lattice orientation in polycrystalline samples, and to calculate the local strain tensor. An implementation is made available under a Free and Open Source Software license. △ Less

Submitted 18 April, 2016; v1 submitted 16 March, 2016; originally announced March 2016.

Comments: 20 pages, 14 figures. Revised version: algorithm improved slightly

Journal ref: Modell. SImul. Mater. Sci. Eng. 24, 055007 (2016)

Control-Flow Integrity: Precision, Security, and Performance

Authors: Nathan Burow, Scott A. Carr, Joseph Nash, Per Larsen, Michael Franz, Stefan Brunthaler, Mathias Payer

Abstract: Memory corruption errors in C/C++ programs remain the most common source of security vulnerabilities in today's systems. Control-flow hijacking attacks exploit memory corruption vulnerabilities to divert program execution away from the intended control flow. Researchers have spent more than a decade studying and refining defenses based on Control-Flow Integrity (CFI), and this technique is now int… ▽ More Memory corruption errors in C/C++ programs remain the most common source of security vulnerabilities in today's systems. Control-flow hijacking attacks exploit memory corruption vulnerabilities to divert program execution away from the intended control flow. Researchers have spent more than a decade studying and refining defenses based on Control-Flow Integrity (CFI), and this technique is now integrated into several production compilers. However, so far no study has systematically compared the various proposed CFI mechanisms, nor is there any protocol on how to compare such mechanisms. We compare a broad range of CFI mechanisms using a unified nomenclature based on (i) a qualitative discussion of the conceptual security guarantees, (ii) a quantitative security evaluation, and (iii) an empirical evaluation of their performance in the same test environment. For each mechanism, we evaluate (i) protected types of control-flow transfers, (ii) the precision of the protection for forward and backward edges. For open-source compiler-based implementations, we additionally evaluate (iii) the generated equivalence classes and target sets, and (iv) the runtime performance. △ Less

Submitted 27 January, 2017; v1 submitted 12 February, 2016; originally announced February 2016.

Comments: Version submitted to ACM CSUR 01/27/17