Showing 1–50 of 78 results for author: Lang, T

A detailed survey of the parallel mean free path of solar energetic particle protons and electrons

Authors: J. T. Lang, R. D. Strauss, N. E. Engelbrecht, J. P. van den Berg, N. Dresing, D. Ruffolo, R. Bandyopadhyay

Abstract: In this work, more than a dozen solar energetic particle (SEP) events are identified where the source region is magnetically well-connected to at least one spacecraft at 1~au. The observed intensity-time profiles, for all available proton and electron energy channels, are compared to results computed using a numerical 1D SEP transport model in order to derive the parallel mean free paths (pMFPs) a… ▽ More In this work, more than a dozen solar energetic particle (SEP) events are identified where the source region is magnetically well-connected to at least one spacecraft at 1~au. The observed intensity-time profiles, for all available proton and electron energy channels, are compared to results computed using a numerical 1D SEP transport model in order to derive the parallel mean free paths (pMFPs) as a function of energy (or rigidity) at 1~au. These inversion results are then compared to theoretical estimates of the pMFP, using observed turbulence quantities with observationally-motivated variations as input. For protons, a very good comparison between inversion and theoretical results is obtained. It is shown that the observed inter-event variations in the inversion pMFP values can be explained by natural variations in the background turbulence values. For electrons, there is relatively good agreement with pMFPs derived assuming the damping model of dynamical turbulence, although the theoretical values are extremely sensitive to the details of the turbulence dissipation range which themselves display a high level of variation. △ Less

Submitted 9 June, 2024; originally announced June 2024.

Comments: Accepted to ApJ

Minimalistic System Modelling: Behaviours, Interfaces, and Local Reasoning

Authors: Didier Galmiche, Timo Lang, David Pym

Abstract: The infrastructure upon which the functioning of society depends is composed of complex ecosystems of systems. Consequently, we must reason about the properties of such ecosystems, which requires that we construct models of them. There are very many approaches to systems modelling, typically building on complex structural and dynamic frameworks. Our purpose here is to explore a modelling framework… ▽ More The infrastructure upon which the functioning of society depends is composed of complex ecosystems of systems. Consequently, we must reason about the properties of such ecosystems, which requires that we construct models of them. There are very many approaches to systems modelling, typically building on complex structural and dynamic frameworks. Our purpose here is to explore a modelling framework based on minimal assumptions, starting from a primitive notion of behaviour, and to show that such an approach allows the recovery of the key ideas, including a generalized CAP theorem, required for effective modelling of and reasoning about ecosystems of systems. We establish a logic of behaviours and use it to express local reasoning principles for the compositional structure of systems. △ Less

Submitted 29 January, 2024; originally announced January 2024.

Comments: 23 pages

MSC Class: 68Q60 ACM Class: F.3.0

Predicting the activity of chemical compounds based on machine learning approaches

Authors: Do Hoang Tu, Tran Van Lang, Pham Cong Xuyen, Le Mau Long

Abstract: Exploring methods and techniques of machine learning (ML) to address specific challenges in various fields is essential. In this work, we tackle a problem in the domain of Cheminformatics; that is, providing a suitable solution to aid in predicting the activity of a chemical compound to the best extent possible. To address the problem at hand, this study conducts experiments on 100 different combi… ▽ More Exploring methods and techniques of machine learning (ML) to address specific challenges in various fields is essential. In this work, we tackle a problem in the domain of Cheminformatics; that is, providing a suitable solution to aid in predicting the activity of a chemical compound to the best extent possible. To address the problem at hand, this study conducts experiments on 100 different combinations of existing techniques. These solutions are then selected based on a set of criteria that includes the G-means, F1-score, and AUC metrics. The results have been tested on a dataset of about 10,000 chemical compounds from PubChem that have been classified according to their activity △ Less

Submitted 10 September, 2023; originally announced January 2024.

Gauge Symmetry Breaking Lattice Regularizations and their Continuum Limit

Authors: Thorsten Lang, Susanne Schander

Abstract: Lattice regularizations are pivotal in the non-perturbative quantization of gauge field theories. Wilson's proposal to employ group-valued link fields simplifies the regularization of gauge fields in principal fiber bundles, preserving gauge symmetry within the discretized lattice theory. Maintaining gauge symmetry is desirable as its violation can introduce unwanted degrees of freedom. However, n… ▽ More Lattice regularizations are pivotal in the non-perturbative quantization of gauge field theories. Wilson's proposal to employ group-valued link fields simplifies the regularization of gauge fields in principal fiber bundles, preserving gauge symmetry within the discretized lattice theory. Maintaining gauge symmetry is desirable as its violation can introduce unwanted degrees of freedom. However, not all theories with gauge symmetries admit gauge-invariant lattice regularizations, as observed in general relativity where the diffeomorphism group serves as the gauge symmetry. In such cases, gauge symmetry-breaking regularizations become necessary. In this paper, we argue that a broken lattice gauge symmetry is acceptable as long as gauge symmetry is restored in the continuum limit. We propose a method to construct the continuum limit for a class of lattice-regularized Hamiltonian field theories, where the regularization breaks the Lie algebra of first-class constraints. Additionally, we offer an approach to represent the exact gauge group on the Hilbert space of the continuum theory. The considered class of theories is limited to those with first-class constraints linear in momenta, excluding the entire gauge group of general relativity but encompassing its subgroup of spatial diffeomorphisms. We discuss potential techniques for extending this quantization to the full gauge group. △ Less

Submitted 31 October, 2023; originally announced November 2023.

The Influence of Satellite Trails on H.E.S.S. Gamma-Ray Astronomical Observations

Authors: Samuel T. Spencer, Thomas Lang, Alison M. W. Mitchell

Abstract: The number of satellites launched into low earth orbit has almost tripled (to over 4000) in the last three years due to the increasing commercialisation of space. Satellite constellations with a total of over 400,000 satellites are proposed to be launched in the near future. Many of these satellites are highly reflective, resulting in a high optical brightness that affects ground-based astronomica… ▽ More The number of satellites launched into low earth orbit has almost tripled (to over 4000) in the last three years due to the increasing commercialisation of space. Satellite constellations with a total of over 400,000 satellites are proposed to be launched in the near future. Many of these satellites are highly reflective, resulting in a high optical brightness that affects ground-based astronomical observations across the electromagnetic spectrum. Despite this, the potential effect of these satellites on Imaging Atmospheric Cherenkov Telescopes (IACTs) has so far been assumed to be negligible due to their nanosecond integration times. This has, however, never been verified. We aim to identify satellite trails in data taken by the High Energy Stereoscopic System (H.E.S.S.) IACT array in Namibia, using Night Sky Background (NSB) data from the CT5 camera installed in 2019. We determine which observation times and pointing directions are affected the most, and evaluate the impact on Hillas parameters used for classification and reconstruction of high-energy Extensive Air Shower events. Finally, we predict how future planned satellite launches will affect gamma-ray observations with IACTs. △ Less

Submitted 2 August, 2023; originally announced August 2023.

Comments: 8 pages, 5 figures, 1 table. Short summary of a full paper which can be found here: arXiv:2307.13293

Journal ref: Proceedings of the 38th International Cosmic Ray Conference (ICRC2023). PoS(ICRC2023)694

Impact of Satellite Trails on H.E.S.S. Astronomical Observations

Authors: Thomas Lang, Samuel T. Spencer, Alison M. W. Mitchell

Abstract: The number of satellites launched into Earth's orbit has almost tripled in the last three years due to the increasing commercialisation of space. Multiple satellite constellations, consisting of over 400,000 individual satellites, have either been partially launched or are proposed for launch in the near future. Many of these satellites are highly reflective, resulting in a high optical brightness… ▽ More The number of satellites launched into Earth's orbit has almost tripled in the last three years due to the increasing commercialisation of space. Multiple satellite constellations, consisting of over 400,000 individual satellites, have either been partially launched or are proposed for launch in the near future. Many of these satellites are highly reflective, resulting in a high optical brightness that affects ground-based astronomical observations. Despite this caveat, the potential effect of these satellites on gamma-ray-observing Imaging Atmospheric Cherenkov Telescopes (IACTs) has largely been assumed to be negligible due to their nanosecond-scale integration times. However, this assumption has not been verified to date. As IACTs are sensitive to optical wavelength light, we aim to identify satellite trails in data taken by the High Energy Stereoscopic System (H.E.S.S.) IACT array. In particular, this study is aimed at quantifying the potential effects on data quality and extensive air shower event classification and reconstruction. Using night sky background measurements from H.E.S.S., we determined which observation times and pointing directions are affected most by these satellite trails. We then evaluated their impact on the standard Hillas parameter variables used for event analysis. Due to the brightest trails, false trigger events can occur, however, for most modern analyses, the effect on astronomical results will be minimal. We observe a mild increase in the rate of trail detections over time, which is partially correlated with the number of satellite launches. Overall, the fraction of H.E.S.S. data affected is currently minimal. We note that these trails could still have a non-negligible effect on future Cherenkov Telescope Array observations if advanced analysis techniques designed to lower the energy threshold of the instrument are applied. △ Less

Submitted 21 September, 2023; v1 submitted 25 July, 2023; originally announced July 2023.

Comments: 11 pages, 14 figures, 5 tables. Replaced with published version. Reproduced with permission from Astronomy & Astrophysics, $©$ ESO

Journal ref: A&A 677, A141 (2023)

The rise of data-driven weather forecasting

Authors: Zied Ben-Bouallegue, Mariana C A Clare, Linus Magnusson, Estibaliz Gascon, Michael Maier-Gerber, Martin Janousek, Mark Rodwell, Florian Pinault, Jesper S Dramsch, Simon T K Lang, Baudouin Raoult, Florence Rabier, Matthieu Chevallier, Irina Sandu, Peter Dueben, Matthew Chantry, Florian Pappenberger

Abstract: Data-driven modeling based on machine learning (ML) is showing enormous potential for weather forecasting. Rapid progress has been made with impressive results for some applications. The uptake of ML methods could be a game-changer for the incremental progress in traditional numerical weather prediction (NWP) known as the 'quiet revolution' of weather forecasting. The computational cost of running… ▽ More Data-driven modeling based on machine learning (ML) is showing enormous potential for weather forecasting. Rapid progress has been made with impressive results for some applications. The uptake of ML methods could be a game-changer for the incremental progress in traditional numerical weather prediction (NWP) known as the 'quiet revolution' of weather forecasting. The computational cost of running a forecast with standard NWP systems greatly hinders the improvements that can be made from increasing model resolution and ensemble sizes. An emerging new generation of ML models, developed using high-quality reanalysis datasets like ERA5 for training, allow forecasts that require much lower computational costs and that are highly-competitive in terms of accuracy. Here, we compare for the first time ML-generated forecasts with standard NWP-based forecasts in an operational-like context, initialized from the same initial conditions. Focusing on deterministic forecasts, we apply common forecast verification tools to assess to what extent a data-driven forecast produced with one of the recently developed ML models (PanguWeather) matches the quality and attributes of a forecast from one of the leading global NWP systems (the ECMWF IFS). The results are very promising, with comparable skill for both global metrics and extreme events, when verified against both the operational analysis and synoptic observations. Increasing forecast smoothness and bias drift with forecast lead time are identified as current drawbacks of ML-based forecasts. A new NWP paradigm is emerging relying on inference from ML models and state-of-the-art analysis and reanalysis datasets for forecast initialization and model training. △ Less

Submitted 3 November, 2023; v1 submitted 19 July, 2023; originally announced July 2023.

Quantum Geometrodynamics Revived I. Classical Constraint Algebra

Authors: Thorsten Lang, Susanne Schander

Abstract: In this series of papers, we present a set of methods to revive quantum geometrodynamics which encountered numerous mathematical and conceptual challenges in its original form promoted by Wheeler and De Witt. In this paper, we introduce the regularization scheme on which we base the subsequent quantization and continuum limit of the theory. Specifically, we employ the set of piecewise constant fie… ▽ More In this series of papers, we present a set of methods to revive quantum geometrodynamics which encountered numerous mathematical and conceptual challenges in its original form promoted by Wheeler and De Witt. In this paper, we introduce the regularization scheme on which we base the subsequent quantization and continuum limit of the theory. Specifically, we employ the set of piecewise constant fields as the phase space of classical geometrodynamics, resulting in a theory with finitely many degrees of freedom of the spatial metric field. As this representation effectively corresponds to a lattice theory, we can utilize well-known techniques to depict the constraints and their algebra on the lattice. We are able to compute the lattice corrections to the constraint algebra. This model can now be quantized using the usual methods of finite--dimensional quantum mechanics, as we demonstrate in the following paper. The application of the continuum limit is the subject of a future publication. △ Less

Submitted 24 June, 2023; v1 submitted 17 May, 2023; originally announced May 2023.

Quantum Geometrodynamics Revived II. Hilbert Space of Positive Definite Metrics

Authors: Thorsten Lang, Susanne Schander

Abstract: This paper represents the second in a series of works aimed at reinvigorating the quantum geometrodynamics program. Our approach introduces a lattice regularization of the hypersurface deformation algebra, such that each lattice site carries a set of canonical variables given by the components of the spatial metric and the corresponding conjugate momenta. In order to quantize this theory, we descr… ▽ More This paper represents the second in a series of works aimed at reinvigorating the quantum geometrodynamics program. Our approach introduces a lattice regularization of the hypersurface deformation algebra, such that each lattice site carries a set of canonical variables given by the components of the spatial metric and the corresponding conjugate momenta. In order to quantize this theory, we describe a representation of the canonical commutation relations that enforces the positivity of the operators $\hat q_{ab} s^a s^b$ for all choices of $s$. Moreover, symmetry of $\hat q_{ab}$ and $\hat p^{ab}$ is ensured. This reflects the physical requirement that the spatial metric should be a positive definite, symmetric tensor. To achieve this end, we resort to the Cholesky decomposition of the spatial metric into upper triangular matrices with positive diagonal entries. Moreover, our Hilbert space also carries a representation of the vielbein fields and naturally separates the physical and gauge degrees of freedom. Finally, we introduce a generalization of the Weyl quantization for our representation. We want to emphasize that our proposed methodology is amenable to applications in other fields of physics, particularly in scenarios where the configuration space is restricted by complicated relationships among the degrees of freedom. △ Less

Submitted 25 June, 2023; v1 submitted 16 May, 2023; originally announced May 2023.

Cut-restriction: from cuts to analytic cuts

Authors: Agata Ciabattoni, Timo Lang, Revantha Ramanayake

Abstract: Cut-elimination is the bedrock of proof theory with a multitude of applications from computational interpretations to proof analysis. It is also the starting point for important meta-theoretical investigations including decidability, complexity, disjunction property, and interpolation. Unfortunately cut-elimination does not hold for the sequent calculi of most non-classical logics. It is well-know… ▽ More Cut-elimination is the bedrock of proof theory with a multitude of applications from computational interpretations to proof analysis. It is also the starting point for important meta-theoretical investigations including decidability, complexity, disjunction property, and interpolation. Unfortunately cut-elimination does not hold for the sequent calculi of most non-classical logics. It is well-known that the key to applications is the subformula property (a typical consequence of cut-elimination) rather than cut-elimination itself. With this in mind we introduce cut-restriction, a procedure to restrict arbitrary cuts to analytic cuts (when elimination is not possible). The algorithm applies to all sequent calculi satisfying language-independent and simple-to-check conditions, and it is obtained by adapting age-old cut-elimination. Our work encompasses existing results in a uniform way, and establishes novel analytic subformula properties. △ Less

Submitted 28 April, 2023; v1 submitted 26 April, 2023; originally announced April 2023.

Comments: 13 pages, conference preprint

Acousto-Optic Modulation in Ambient Air

Authors: Yannick Schrödel, Claas Hartmann, Tino Lang, Jiaan Zheng, Max Steudel, Matthias Rutsch, Sarper H. Salman, Martin Kellert, Mikhail Pergament, Thomas Hahn-Jose, Sven Suppelt, Jan Helge Dörsam, Anne Harth, Wim P. Leemans, Franz X. Kärtner, Ingmar Hartl, Mario Kupnik, Christoph M. Heyl

Abstract: Control over intensity, shape, direction, and phase of coherent light is essential in numerous fields, reaching from gravitational wave astronomy over quantum metrology and ultrafast sciences to semi-conductor fabrication. Modern laser optics, however, frequently demands parameter regimes where either the wavelength or the optical power restricts control due to linear absorption, light-induced dam… ▽ More Control over intensity, shape, direction, and phase of coherent light is essential in numerous fields, reaching from gravitational wave astronomy over quantum metrology and ultrafast sciences to semi-conductor fabrication. Modern laser optics, however, frequently demands parameter regimes where either the wavelength or the optical power restricts control due to linear absorption, light-induced damage or optical nonlinearity. The properties of solid media, upon which most photonic control schemes rely, impose these limitations. We propose to circumvent these constraints using gaseous media tailored by high-intensity ultrasound waves. We demonstrate a first implementation of this approach by deflecting ultrashort laser pulses using ultrasound waves in ambient air, entirely omitting transmissive solid media. At optical peak powers of 20 GW exceeding previous limits of solid-based acousto-optic modulation by about three orders of magnitude, we reach a deflection efficiency greater than 50% while preserving excellent beam quality. Our approach is not limited to laser pulse deflection via acousto-optic modulation: gas-phase photonic schemes controlled by sonic waves can prospectively be translated to various optical methods, e.g., lenses or waveguides, rendering them effectively invulnerable against damage and opening up new spectral regions. △ Less

Submitted 28 April, 2023; v1 submitted 13 April, 2023; originally announced April 2023.

Comments: 23 pages, including 9 pages of main text and 9 pages of Methods and Extended Data, 3 figures, 3 Extended Data figures, 1 Extended Data table

Report number: PUBDB-2023-01719

Clustering large 3D volumes: A sampling-based approach

Authors: Thomas Lang

Abstract: In many applications of X-ray computed tomography, an unsupervised segmentation of the reconstructed 3D volumes forms an important step in the image processing chain for further investigation of the digitized object. Therefore, the goal is to train a clustering algorithm on the volume, which produces a voxelwise classification by assigning a cluster index to each voxel. However, clustering methods… ▽ More In many applications of X-ray computed tomography, an unsupervised segmentation of the reconstructed 3D volumes forms an important step in the image processing chain for further investigation of the digitized object. Therefore, the goal is to train a clustering algorithm on the volume, which produces a voxelwise classification by assigning a cluster index to each voxel. However, clustering methods, e.g., K-Means, typically have an asymptotic polynomial runtime with respect to the dataset size, and thus, these techniques are rarely applicable to large volumes. In this work, we introduce a novel clustering technique based on random sampling, which allows for the voxelwise classification of arbitrarily large volumes. The presented method conducts efficient linear passes over the data to extract a representative random sample of a fixed size on which the classifier can be trained. Then, a final linear pass performs the segmentation and assigns a cluster index to each individual voxel. Quantitative and qualitative evaluations show that excellent results can be achieved even with a very small sample size. Consequently, the unsupervised segmentation by means of clustering becomes feasible for arbitrarily large volumes. △ Less

Submitted 7 March, 2023; originally announced March 2023.

Comments: 12 pages, 8 figures

MSC Class: 68W20 ACM Class: I.5.3; H.3.3

Feature-Adaptive Interactive Thresholding of Large 3D Volumes

Authors: Thomas Lang, Tomas Sauer

Abstract: Thresholding is the most widely used segmentation method in volumetric image processing, and its pointwise nature makes it attractive for the fast handling of large three-dimensional samples. However, global thresholds often do not properly extract components in the presence of artifacts, measurement noise or grayscale value fluctuations. This paper introduces Feature-Adaptive Interactive Threshol… ▽ More Thresholding is the most widely used segmentation method in volumetric image processing, and its pointwise nature makes it attractive for the fast handling of large three-dimensional samples. However, global thresholds often do not properly extract components in the presence of artifacts, measurement noise or grayscale value fluctuations. This paper introduces Feature-Adaptive Interactive Thresholding (FAITH), a thresholding technique that incorporates (geometric) features, local processing and interactive user input to overcome these limitations. Given a global threshold suitable for most regions, FAITH uses interactively selected seed voxels to identify critical regions in which that threshold will be adapted locally on the basis of features computed from local environments around these voxels. The combination of domain expert knowledge and a rigorous mathematical model thus enables a very exible way of local thresholding with intuitive user interaction. A qualitative analysis shows that the proposed model is able to overcome limitations typically occuring in plain thresholding while staying efficient enough to also allow the segmentation of big volumes. △ Less

Submitted 13 October, 2022; originally announced October 2022.

Comments: 13 pages, 3 figures, 2 tables

MSC Class: 68U10 ACM Class: I.4.6

Geometric Active Learning for Segmentation of Large 3D Volumes

Authors: Thomas Lang, Tomas Sauer

Abstract: Segmentation, i.e., the partitioning of volumetric data into components, is a crucial task in many image processing applications ever since such data could be generated. Most existing applications nowadays, specifically CNNs, make use of voxelwise classification systems which need to be trained on a large number of annotated training volumes. However, in many practical applications such data sets… ▽ More Segmentation, i.e., the partitioning of volumetric data into components, is a crucial task in many image processing applications ever since such data could be generated. Most existing applications nowadays, specifically CNNs, make use of voxelwise classification systems which need to be trained on a large number of annotated training volumes. However, in many practical applications such data sets are seldom available and the generation of annotations is time-consuming and cumbersome. In this paper, we introduce a novel voxelwise segmentation method based on active learning on geometric features. Our method uses interactively provided seed points to train a voxelwise classifier based entirely on local information. The combination of an ad hoc incorporation of domain knowledge and local processing results in a flexible yet efficient segmentation method that is applicable to three-dimensional volumes without size restrictions. We illustrate the potential and flexibility of our approach by applying it to selected computed tomography scans where we perform different segmentation tasks to scans from different domains and of different sizes. △ Less

Submitted 13 October, 2022; originally announced October 2022.

Comments: 10 pages, 27 figures, 3 tables

MSC Class: 68T10 ACM Class: I.5.2

A New Hip Fracture Risk Index Derived from FEA-Computed Proximal Femur Fracture Loads and Energies-to-Failure

Authors: Xuewei Cao, Joyce H Keyak, Sigurdur Sigurdsson, Chen Zhao, Weihua Zhou, Anqi Liu, Thomas Lang, Hong-Wen Deng, Vilmundur Gudnason, Qiuying Sha

Abstract: Hip fracture risk assessment is an important but challenging task. Quantitative CT-based patient specific finite element analysis (FEA) computes the force (fracture load) to break the proximal femur in a particular loading condition. It provides different structural information about the proximal femur that can influence a subject overall fracture risk. To obtain a more robust measure of fracture… ▽ More Hip fracture risk assessment is an important but challenging task. Quantitative CT-based patient specific finite element analysis (FEA) computes the force (fracture load) to break the proximal femur in a particular loading condition. It provides different structural information about the proximal femur that can influence a subject overall fracture risk. To obtain a more robust measure of fracture risk, we used principal component analysis (PCA) to develop a global FEA computed fracture risk index that incorporates the FEA-computed yield and ultimate failure loads and energies to failure in four loading conditions (single-limb stance and impact from a fall onto the posterior, posterolateral, and lateral aspects of the greater trochanter) of 110 hip fracture subjects and 235 age and sex matched control subjects from the AGES-Reykjavik study. We found that the first PC (PC1) of the FE parameters was the only significant predictor of hip fracture. Using a logistic regression model, we determined if prediction performance for hip fracture using PC1 differed from that using FE parameters combined by stratified random resampling with respect to hip fracture status. The results showed that the average of the area under the receive operating characteristic curve (AUC) using PC1 was always higher than that using all FE parameters combined in the male subjects. The AUC of PC1 and AUC of the FE parameters combined were not significantly different than that in the female subjects or in all subjects △ Less

Submitted 18 November, 2022; v1 submitted 3 October, 2022; originally announced October 2022.

Comments: 27 pages, 4 figures

Recommendations on test datasets for evaluating AI solutions in pathology

Authors: André Homeyer, Christian Geißler, Lars Ole Schwen, Falk Zakrzewski, Theodore Evans, Klaus Strohmenger, Max Westphal, Roman David Bülow, Michaela Kargl, Aray Karjauv, Isidre Munné-Bertran, Carl Orge Retzlaff, Adrià Romero-López, Tomasz Sołtysiński, Markus Plass, Rita Carvalho, Peter Steinbach, Yu-Chia Lan, Nassim Bouteldja, David Haber, Mateo Rojas-Carulla, Alireza Vafaei Sadr, Matthias Kraft, Daniel Krüger, Rutger Fick , et al. (5 additional authors not shown)

Abstract: Artificial intelligence (AI) solutions that automatically extract information from digital histology images have shown great promise for improving pathological diagnosis. Prior to routine use, it is important to evaluate their predictive performance and obtain regulatory approval. This assessment requires appropriate test datasets. However, compiling such datasets is challenging and specific recom… ▽ More Artificial intelligence (AI) solutions that automatically extract information from digital histology images have shown great promise for improving pathological diagnosis. Prior to routine use, it is important to evaluate their predictive performance and obtain regulatory approval. This assessment requires appropriate test datasets. However, compiling such datasets is challenging and specific recommendations are missing. A committee of various stakeholders, including commercial AI developers, pathologists, and researchers, discussed key aspects and conducted extensive literature reviews on test datasets in pathology. Here, we summarize the results and derive general recommendations for the collection of test datasets. We address several questions: Which and how many images are needed? How to deal with low-prevalence subsets? How can potential bias be detected? How should datasets be reported? What are the regulatory requirements in different countries? The recommendations are intended to help AI developers demonstrate the utility of their products and to help regulatory agencies and end users verify reported performance measures. Further research is needed to formulate criteria for sufficiently representative test datasets so that AI solutions can operate with less user intervention and better support diagnostic workflows in the future. △ Less

Submitted 21 April, 2022; originally announced April 2022.

Journal ref: Mod Pathol (2022)

Non-coplanar magnetism, topological density wave order and emergent symmetry at half-integer filling of moiré Chern bands

Authors: Patrick H. Wilhelm, Thomas C. Lang, Mathias S. Scheurer, Andreas M. Läuchli

Abstract: Twisted double- and mono-bilayer graphene are graphene-based moiré materials hosting strongly correlated fermions in a gate-tunable conduction band with a topologically non-trivial character. Using unbiased exact diagonalization complemented by unrestricted Hartree-Fock calculations, we find that the strong electron-electron interactions lead to a non-coplanar magnetic state, which has the same sy… ▽ More Twisted double- and mono-bilayer graphene are graphene-based moiré materials hosting strongly correlated fermions in a gate-tunable conduction band with a topologically non-trivial character. Using unbiased exact diagonalization complemented by unrestricted Hartree-Fock calculations, we find that the strong electron-electron interactions lead to a non-coplanar magnetic state, which has the same symmetries as the tetrahedral antiferromagnet on the triangular lattice and can be thought of as a skyrmion lattice commensurate with the moiré scale, competing with a set of ferromagnetic, topological charge density waves featuring an approximate emergent O(3) symmetry, "rotating" the different charge density wave states into each other. Direct comparison with exact diagonalization reveals that the ordered phases are accurately described within the unrestricted Hartree-Fock approximation. Exhibiting a finite charge gap and Chern number $|C|=1$, the formation of charge density wave order which is intimately connected to a skyrmion lattice phase is consistent with recent experiments on these systems. △ Less

Submitted 20 March, 2023; v1 submitted 11 April, 2022; originally announced April 2022.

Comments: 35 pages, 12 figures

Journal ref: SciPost Phys. 14, 040 (2023)

A theory of cut-restriction: first steps

Authors: Agata Ciabattoni, Timo Lang, Revantha Ramanayake

Abstract: Cut-elimination is the bedrock of proof theory. It is the algorithm that eliminates cuts from a sequent calculus proof that leads to cut-free calculi and applications. Cut-elimination applies to many logics irrespective of their semantics. Such is its influence that whenever cut-elimination is not provable in a sequent calculus the invariable response has been a move to a richer proof system to re… ▽ More Cut-elimination is the bedrock of proof theory. It is the algorithm that eliminates cuts from a sequent calculus proof that leads to cut-free calculi and applications. Cut-elimination applies to many logics irrespective of their semantics. Such is its influence that whenever cut-elimination is not provable in a sequent calculus the invariable response has been a move to a richer proof system to regain it. In this paper we investigate a radically different approach to the latter: adapting age-old cut-elimination to restrict the shape of the cut-formulas when elimination is not possible. We tackle the "first level" above cut-free: analytic cuts. Our methodology is applied to the sequent calculi for bi-intuitionistic logic and S5 where analytic cuts are already known to be required. This marks the first steps in a theory of cut-restriction. △ Less

Submitted 3 March, 2022; originally announced March 2022.

Comments: 18 pages

Real-time Street Human Motion Capture

Authors: Yanquan Chen, Fei Yang, Tianyu Lang, Guanfang Dong, Anup Basu

Abstract: In recent years, motion capture technology using computers has developed rapidly. Because of its high efficiency and excellent performance, it replaces many traditional methods and is being widely used in many fields. Our project is about street scene video human motion capturing and analysis. The primary goal of the project is to capture the human motion in a video and use the motion information… ▽ More In recent years, motion capture technology using computers has developed rapidly. Because of its high efficiency and excellent performance, it replaces many traditional methods and is being widely used in many fields. Our project is about street scene video human motion capturing and analysis. The primary goal of the project is to capture the human motion in a video and use the motion information for 3D animation (human) in real-time. We applied a neural network for motion capture and implement it in the unity under a street view scene. By analyzing the motion data, we will have a better estimation of the street condition, which is useful for other high-tech applications such as self-driving cars. △ Less

Submitted 21 December, 2021; originally announced December 2021.

Comments: 7 pages, 7 figures

Dynamic Lockstep Processors for Applications with Functional Safety Relevance

Authors: Hans Dermot Doran, Timo Lang

Abstract: Lockstep processing is a recognized technique for helping to secure functional-safety relevant processing against, for instance, single upset errors that might cause faulty execution of code. Lockstepping processors does however bind processing resources in a fashion not beneficial to architectures and applications that would benefit from multi-core/-processors. We propose a novel on-demand synchr… ▽ More Lockstep processing is a recognized technique for helping to secure functional-safety relevant processing against, for instance, single upset errors that might cause faulty execution of code. Lockstepping processors does however bind processing resources in a fashion not beneficial to architectures and applications that would benefit from multi-core/-processors. We propose a novel on-demand synchronizing of cores/processors for lock-step operation featuring post-processing resource release, a concept that facilitates the implementation of modularly redundant core/processor arrays. We discuss the fundamentals of the design and some implementation notes on work achieved to date. △ Less

Submitted 19 July, 2021; originally announced July 2021.

Comments: 4 pages, 8 figures

Cryogenic Penning-Trap Apparatus for Precision Experiments with Sympathetically Cooled (anti)protons

Authors: M. Niemann, T. Meiners, J. Mielke, N. Pulido, J. Schaper, M. J. Borchert, J. M. Cornejo, A. -G. Paschke, G. Zarantonello, H. Hahn, T. Lang, C. Manzoni, M. Marangoni, G. Cerullo, U. Morgner, J. -A. Fenske, A. Bautista-Salvador, R. Lehnert, S. Ulmer, C. Ospelkaus

Abstract: Current precision experiments with single (anti)protons to test CPT symmetry progress at a rapid pace, but are complicated by the need to cool particles to sub-thermal energies. We describe a cryogenic Penning-trap setup for $^9$Be$^+$ ions designed to allow coupling of single (anti)protons to laser-cooled atomic ions for sympathetic cooling and quantum logic spectroscopy. We report on trapping an… ▽ More Current precision experiments with single (anti)protons to test CPT symmetry progress at a rapid pace, but are complicated by the need to cool particles to sub-thermal energies. We describe a cryogenic Penning-trap setup for $^9$Be$^+$ ions designed to allow coupling of single (anti)protons to laser-cooled atomic ions for sympathetic cooling and quantum logic spectroscopy. We report on trapping and laser cooling of clouds and single $^9$Be$^+$ ions. We discuss prospects for a microfabricated trap to allow coupling of single (anti)protons to laser-cooled $^9$Be$^+$ ions for sympathetic laser cooling to sub-mK temperatures on ms time scales. △ Less

Submitted 18 July, 2021; originally announced July 2021.

Comments: Presented at the Eighth Meeting on CPT and Lorentz Symmetry, Bloomington, Indiana, May 12-16, 2019

Journal ref: World Scientific, Singapore, Proceedings of the Eighth Meeting on CPT and Lorentz Symmetry (2020)

Decidability and Complexity in Weakening and Contraction Hypersequent Substructural Logics

Authors: A. R. Balasubramanian, Timo Lang, Revantha Ramanayake

Abstract: We establish decidability for the infinitely many axiomatic extensions of the commutative Full Lambek logic with weakening FLew (i.e. IMALLW) that have a cut-free hypersequent proof calculus (specifically: every analytic structural rule extension). Decidability for the corresponding extensions of its contraction counterpart FLec was established recently but their computational complexity was left… ▽ More We establish decidability for the infinitely many axiomatic extensions of the commutative Full Lambek logic with weakening FLew (i.e. IMALLW) that have a cut-free hypersequent proof calculus (specifically: every analytic structural rule extension). Decidability for the corresponding extensions of its contraction counterpart FLec was established recently but their computational complexity was left unanswered. In the second part of this paper, we introduce just enough on length functions for well-quasi-orderings and the fast-growing complexity classes to obtain complexity upper bounds for both the weakening and contraction extensions. A specific instance of this result yields the first complexity bound for the prominent fuzzy logic MTL (monoidal t-norm based logic) providing an answer to a long-standing open problem. △ Less

Submitted 19 April, 2021; originally announced April 2021.

Comments: Accepted for publication in the proceedings of LICS 2021

Comparing spontaneous and pellet-triggered ELMs via non-linear extended MHD simulations

Authors: A. Cathey, M. Hoelzl, S. Futatani, P. T. Lang, K. Lackner, G. T. A. Huijsmans, S. J. P. Pamela, S. Günter, the JOREK team, the ASDEX Upgrade Team, the EUROfusion MST1 Team

Abstract: Injecting frozen deuterium pellets into an ELMy H-mode plasma is a well established scheme for triggering edge localized modes (ELMs) before they naturally occur. Based on an ASDEX Upgrade H-mode plasma, this article presents a comparison of extended MHD simulations of spontaneous type-I ELMs and pellet-triggered ELMs allowing to study their non-linear dynamics in detail. In particular, pellet-tri… ▽ More Injecting frozen deuterium pellets into an ELMy H-mode plasma is a well established scheme for triggering edge localized modes (ELMs) before they naturally occur. Based on an ASDEX Upgrade H-mode plasma, this article presents a comparison of extended MHD simulations of spontaneous type-I ELMs and pellet-triggered ELMs allowing to study their non-linear dynamics in detail. In particular, pellet-triggered ELMs are simulated by injecting deuterium pellets into different time points during the pedestal build-up described in [A. Cathey et al. Nuclear Fusion 60, 124007 (2020)]. Realistic ExB and diamagnetic background plasma flows as well as the time dependent bootstrap current evolution are included during the build-up to capture the balance between stabilising and destabilising terms for the edge instabilities accurately. Dependencies on the pellet size and injection times are studied. The spatio-temporal structures of the modes and the resulting divertor heat fluxes are compared in detail between spontaneous and triggered ELMs. We observe that the premature excitation of ELMs by means of pellet injection is caused by a helical perturbation described by a toroidal mode number of n = 1. In accordance with experimental observations, the pellet-triggered ELMs show reduced thermal energy losses and narrower divertor wetted area with respect to spontaneous ELMs. The peak divertor energy fluency is seen to decrease when ELMs are triggered by pellets injected earlier during the pedestal build-up. △ Less

Submitted 11 February, 2021; originally announced February 2021.

Programmable quantum simulation of 2D antiferromagnets with hundreds of Rydberg atoms

Authors: Pascal Scholl, Michael Schuler, Hannah J. Williams, Alexander A. Eberharter, Daniel Barredo, Kai-Niklas Schymik, Vincent Lienhard, Louis-Paul Henry, Thomas C. Lang, Thierry Lahaye, Andreas M. Läuchli, Antoine Browaeys

Abstract: Quantum simulation using synthetic systems is a promising route to solve outstanding quantum many-body problems in regimes where other approaches, including numerical ones, fail. Many platforms are being developed towards this goal, in particular based on trapped ions, superconducting circuits, neutral atoms or molecules. All of which face two key challenges: (i) scaling up the ensemble size, whil… ▽ More Quantum simulation using synthetic systems is a promising route to solve outstanding quantum many-body problems in regimes where other approaches, including numerical ones, fail. Many platforms are being developed towards this goal, in particular based on trapped ions, superconducting circuits, neutral atoms or molecules. All of which face two key challenges: (i) scaling up the ensemble size, whilst retaining high quality control over the parameters and (ii) certifying the outputs for these large systems. Here, we use programmable arrays of individual atoms trapped in optical tweezers, with interactions controlled by laser-excitation to Rydberg states to implement an iconic many-body problem, the antiferromagnetic 2D transverse field Ising model. We push this platform to an unprecedented regime with up to 196 atoms manipulated with high fidelity. We probe the antiferromagnetic order by dynamically tuning the parameters of the Hamiltonian. We illustrate the versatility of our platform by exploring various system sizes on two qualitatively different geometries, square and triangular arrays. We obtain good agreement with numerical calculations up to a computationally feasible size (around 100 particles). This work demonstrates that our platform can be readily used to address open questions in many-body physics. △ Less

Submitted 22 December, 2020; originally announced December 2020.

Comments: Main text: 6 pages, 4 figures. Supplementary information: 10 pages, 16 figures

Journal ref: Nature 595, 233 (2021)

Interplay of Fractional Chern Insulator and Charge-Density-Wave Phases in Twisted Bilayer Graphene

Authors: Patrick Wilhelm, Thomas C. Lang, Andreas M. Läuchli

Abstract: We perform an extensive exact diagonalization study of interaction driven insulators in spin- and valley-polarized moiré flat bands of twisted bilayer graphene aligned with its hexagonal boron nitride substrate. In addition to previously reported fractional Chern insulator phases, we provide compelling evidence for competing charge-density-wave phases at multiple fractional fillings of a realistic… ▽ More We perform an extensive exact diagonalization study of interaction driven insulators in spin- and valley-polarized moiré flat bands of twisted bilayer graphene aligned with its hexagonal boron nitride substrate. In addition to previously reported fractional Chern insulator phases, we provide compelling evidence for competing charge-density-wave phases at multiple fractional fillings of a realistic single-band model. A thorough analysis at different interlayer hopping parameters, motivated by experimental variability, and the role of kinetic energy at various Coulomb interaction strengths highlight the competition between these phases. The interplay of the single-particle and the interaction induced hole dispersion with the inherent Berry curvature of the Chern bands is intuitively understood to be the driving mechanism for the ground-state selection. The resulting phase diagram features remarkable agreement with experimental findings in a related moiré heterostructure and affirms the relevance of our results beyond the scope of graphene based materials. △ Less

Submitted 4 March, 2021; v1 submitted 17 December, 2020; originally announced December 2020.

Comments: 16 pages, 18 figures

Journal ref: Phys. Rev. B 103, 125406 (2021)

Transition from no-ELM response to pellet ELM triggering during pedestal build-up -- insights from extended MHD simulations

Authors: S Futatani, A Cathey, M Hoelzl, P T Lang, G T A Huijsmans, M. Dunne, JOREK Team, ASDEX Upgrade Team, EUROfusion MST1 Team

Abstract: Pellet ELM triggering is a well established scheme for decreasing the time between two successive ELM crashes below its natural value. Reliable ELM pacing has been demonstrated experimentally in several devices increasing the ELM frequency considerably. However, it was also shown that the frequency cannot be increased arbitrarily due to a so-called lag-time. During this time after a preceding natu… ▽ More Pellet ELM triggering is a well established scheme for decreasing the time between two successive ELM crashes below its natural value. Reliable ELM pacing has been demonstrated experimentally in several devices increasing the ELM frequency considerably. However, it was also shown that the frequency cannot be increased arbitrarily due to a so-called lag-time. During this time after a preceding natural or triggered ELM crash, neither a natural ELM crash occurs nor the triggering of an ELM crash by pellet injection is possible. For this article, pellet ELM triggering simulations are advanced beyond previous studies in two ways. Firstly, realistic ExB and diamagnetic background flows are included. And secondly, the pellet is injected at different stages of the pedestal build-up. This allows to recover the lag-time for the first time in simulations and investigate it in detail. A series of non-linear extended MHD simulations is performed to investigate the plasma dynamics resulting from an injection at different time points during the pedestal build-up. The experimentally observed lag-time is qualitatively reproduced well. In particular, a sharp transition is observed between the regime where no ELMs can be triggered and the regime where pellet injection causes an ELM crash. Via variations of pellet parameters and injection time, the two regimes are studied and compared in detail revealing pronounced differences in the non-linear dynamics. The toroidal mode spectrum is significantly broader when an ELM crash is triggered enhancing the stochasticity and therefore also the losses of thermal energy along magnetic field lines. In the heat fluxes to the divertor targets, pronounced toroidal asymmetries are observed. In case of high injection velocities leading to deep penetration, also the excitation of core modes like the $2/1$ neoclassical tearing mode is observed. △ Less

Submitted 24 December, 2020; v1 submitted 16 September, 2020; originally announced September 2020.

Comments: 45 pages

Post-compression of picosecond pulses into the few-cycle regime

Authors: Prannay Balla, Ammar Bin Wahid, Ivan Sytcevich, Chen Guo, Anne-Lise Viotti, Laura Silletti, Andrea Cartella, Skirmantas Alisauskas, Hamed Tavakol, Uwe Grosse-Wortmann, Arthur Schönberg, Marcus Seidel, Andrea Trabattoni, Bastian Manschwetus, Tino Lang, Francesca Calegari, Arnaud Couairon, Anne L'Huillier, Cord L. Arnold, Ingmar Hartl, Christoph M. Heyl

Abstract: In this work, we demonstrate post-compression of 1.2 picosecond laser pulses to 13 fs via gas-based multi-pass spectral broadening. Our results yield a single-stage compression factor of about 40 at 200 W in-burst average power and a total compression factor >90 at reduced power. The employed scheme represents a route towards compact few-cycle sources driven by industrial-grade Yb:YAG lasers at hi… ▽ More In this work, we demonstrate post-compression of 1.2 picosecond laser pulses to 13 fs via gas-based multi-pass spectral broadening. Our results yield a single-stage compression factor of about 40 at 200 W in-burst average power and a total compression factor >90 at reduced power. The employed scheme represents a route towards compact few-cycle sources driven by industrial-grade Yb:YAG lasers at high average power. △ Less

Submitted 24 March, 2020; originally announced March 2020.

Comment on "The role of electron-electron interactions in two-dimensional Dirac fermions''

Authors: Stephan Hesselmann, Thomas C. Lang, Michael Schuler, Stefan Wessel, Andreas M. Läuchli

Abstract: Tang et al. [Science 361, 570 (2018)] report on the properties of Dirac fermions with both on-site and Coulomb interactions. The substantial decrease up to ~40% of the Fermi velocity of Dirac fermions with on-site interaction is inconsistent with the numerical data near the Gross-Neveu quantum critical point. This results from an inappropriate finite-size extrapolation. Tang et al. [Science 361, 570 (2018)] report on the properties of Dirac fermions with both on-site and Coulomb interactions. The substantial decrease up to ~40% of the Fermi velocity of Dirac fermions with on-site interaction is inconsistent with the numerical data near the Gross-Neveu quantum critical point. This results from an inappropriate finite-size extrapolation. △ Less

Submitted 20 December, 2019; originally announced December 2019.

Comments: Submitted: October 24, 2018, Accepted: November 4, 2019

Journal ref: Science 366, eaav8877 (2019)

Quantifying the fragility of unprotected quadratic band crossing points

Authors: Stephan Hesselmann, Carsten Honerkamp, Stefan Wessel, Thomas C. Lang

Abstract: We examine a basic lattice model of interacting fermions that exhibits quadratic band crossing points (QBCPs) in the non-interacting limit. In particular, we consider spinless fermions on the honeycomb lattice with nearest neighbor hopping $t$ and third-nearest neighbor hopping $t''$, which exhibits fine-tuned QBCPs at the corners of the Brillouin zone for ${t'' = t/2}$. In this situation, the den… ▽ More We examine a basic lattice model of interacting fermions that exhibits quadratic band crossing points (QBCPs) in the non-interacting limit. In particular, we consider spinless fermions on the honeycomb lattice with nearest neighbor hopping $t$ and third-nearest neighbor hopping $t''$, which exhibits fine-tuned QBCPs at the corners of the Brillouin zone for ${t'' = t/2}$. In this situation, the density of states remains finite at the Fermi level of the half-filled band and repulsive nearest-neighbor interactions $V$ lead to a charge-density-wave (CDW) instability at infinitesimally small $V$ in the random-phase approximation or mean-field theory. We examine the fragility of the QBCPs against dispersion renormalizations in the ${t\mbox{-}t''\mbox{-}V}$ model using perturbation theory, and find that the $t''$-value needed for the QBCPs increases with $V$ due to the hopping renormalization. However, the instability toward CDW formation always requires a nonzero threshold interaction strength, i.e., one cannot fine-tune $t''$ to recover the QBCPs in the interacting system. These perturbative arguments are supported by quantum Monte Carlo simulations for which we carefully compare the corresponding threshold scales at and beyond the QBCP fine-tuning point. From this analysis, we thus gain a quantitative microscopic understanding of the fragility of the QBCPs in this basic interacting fermion system. △ Less

Submitted 20 February, 2020; v1 submitted 13 December, 2019; originally announced December 2019.

Comments: 9 pages, 10 figures

Journal ref: Phys. Rev. B 101, 075128 (2020)

Torus Spectroscopy of the Gross-Neveu-Yukawa Quantum Field Theory: Free Dirac versus Chiral Ising Fixed Point

Authors: Michael Schuler, Stephan Hesselmann, Seth Whitsitt, Thomas C. Lang, Stefan Wessel, Andreas M. Läuchli

Abstract: We establish the universal torus low-energy spectra at the free Dirac fixed point and at the strongly coupled chiral Ising fixed point and their subtle crossover behaviour in the Gross-Neuveu-Yukawa field theory with ${n_\text{D}=4}$ component Dirac spinors in $D=(2+1)$ dimensions. These fixed points and the field theories are directly relevant for the long-wavelength physics of certain interactin… ▽ More We establish the universal torus low-energy spectra at the free Dirac fixed point and at the strongly coupled chiral Ising fixed point and their subtle crossover behaviour in the Gross-Neuveu-Yukawa field theory with ${n_\text{D}=4}$ component Dirac spinors in $D=(2+1)$ dimensions. These fixed points and the field theories are directly relevant for the long-wavelength physics of certain interacting Dirac systems, such as repulsive spinless fermions on the honeycomb lattice or $π$-flux square lattice. The torus energy spectrum has been shown previously to serve as a characteristic fingerprint of relativistic fixed points and is a powerful tool to discriminate quantum critical behaviour in numerical simulations. Here, we use a combination of exact diagonalization and quantum Monte Carlo simulations of strongly interacting fermionic lattice models, to compute the critical torus energy spectrum on finite-size clusters with periodic boundaries and extrapolate them to the thermodynamic limit. Additionally, we compute the torus energy spectrum analytically using the perturbative expansion in ${ε= 4 - D}$, which is in good agreement with the numerical results, thereby validating the presence of the chiral Ising fixed point in the lattice models at hand. We show that the strong interaction between the spinor field and the scalar order-parameter field strongly influences the critical torus energy spectrum and we observe prominent multiplicity features related to an emergent symmetry predicted from the quantum field theory. Building on these results we are able to address the subtle crossover physics of the low-energy spectrum flowing from the chiral Ising fixed point to the Dirac fixed point, and analyze earlier flawed attempts to extract Fermi velocity renormalizations from the low-energy spectrum. △ Less

Submitted 15 March, 2021; v1 submitted 11 July, 2019; originally announced July 2019.

Comments: 26 pages, 14 figures

Journal ref: Phys. Rev. B 103, 125128 (2021)

A Game Model for Proofs with Costs

Authors: Timo Lang, Carlos Olarte, Elaine Pimentel, Christian Fermuller

Abstract: We look at substructural calculi from a game semantic point of view, guided by certain intuitions about resource conscious and, more specifically, cost conscious reasoning. To this aim, we start with a game, where player I defends a claim corresponding to a (single-conclusion) sequent, while player II tries to refute that claim. Branching rules for additive connectives are modeled by choices of II… ▽ More We look at substructural calculi from a game semantic point of view, guided by certain intuitions about resource conscious and, more specifically, cost conscious reasoning. To this aim, we start with a game, where player I defends a claim corresponding to a (single-conclusion) sequent, while player II tries to refute that claim. Branching rules for additive connectives are modeled by choices of II, while branching for multiplicative connectives leads to splitting the game into parallel subgames, all of which have to be won by player I to succeed. The game comes into full swing by adding cost labels to assumptions, and a corresponding budget. Different proofs of the same end-sequent are interpreted as more or less expensive strategies for I to defend the corresponding claim. This leads to a new kind of labelled calculus, which can be seen as a fragment of SELL (subexponential linear logic). Finally, we generalize the concept of costs in proofs by using a semiring structure, illustrate our interpretation by examples and investigate some proof-theoretical properties. △ Less

Submitted 27 June, 2019; originally announced June 2019.

Comments: To appear in TABLEAUX'19

Versatile control of $^9$Be$^+$ ions using a spectrally tailored UV frequency comb

Authors: A. -G. Paschke, G. Zarantonello, H. Hahn, T. Lang, C. Manzoni, M. Marangoni, G. Cerullo, U. Morgner, C. Ospelkaus

Abstract: We demonstrate quantum control of $^9$Be$^+$ ions directly implemented by an optical frequency comb. Based on numerical simulations of the relevant processes in $^9$Be$^+$ for different magnetic field regimes, we demonstrate a wide applicability when controlling the comb's spectral properties. We introduce a novel technique for the selective and efficient generation of a spectrally tailored narrow… ▽ More We demonstrate quantum control of $^9$Be$^+$ ions directly implemented by an optical frequency comb. Based on numerical simulations of the relevant processes in $^9$Be$^+$ for different magnetic field regimes, we demonstrate a wide applicability when controlling the comb's spectral properties. We introduce a novel technique for the selective and efficient generation of a spectrally tailored narrow-bandwidth optical frequency comb near 313 nm. We experimentally demonstrate internal state control and internal-motional state coupling of $^9$Be$^+$ ions implemented by stimulated-Raman manipulation using a spectrally optimized optical frequency comb. Our pulsed laser approach is a key enabling step for the implementation of quantum logic and quantum information experiments in Penning traps. △ Less

Submitted 7 March, 2019; originally announced March 2019.

Comments: Includes supplementary material

Multi-microjoule GaSe-based mid-infrared optical parametric amplifier with an ultra-broad idler spectrum covering 4.2-16 μm

Authors: Kun Liu, Houkun Liang, Lifeng Wang, Shizhen Qu, Tino Lang, Hao Li, Qi Jie Wang, Ying Zhang

Abstract: We report a multi-microjoule, ultra-broadband mid-infrared optical parametric amplifier based on a GaSe nonlinear crystal pumped at ~2 μm. The generated idler pulse has a flat spectrum spanning from 4.5 to 13.3 μm at -3 dB and 4.2 to 16 μm in the full spectral range, with a central wavelength of 8.8 μm. The proposed scheme supports a sub-cycle Fourier-transform-limited pulse width. A (2+1)-dimensi… ▽ More We report a multi-microjoule, ultra-broadband mid-infrared optical parametric amplifier based on a GaSe nonlinear crystal pumped at ~2 μm. The generated idler pulse has a flat spectrum spanning from 4.5 to 13.3 μm at -3 dB and 4.2 to 16 μm in the full spectral range, with a central wavelength of 8.8 μm. The proposed scheme supports a sub-cycle Fourier-transform-limited pulse width. A (2+1)-dimensional numerical simulation is employed to reproduce the obtained idler spectrum. To our best knowledge, this is the broadest -3 dB spectrum ever obtained by optical parametric amplifiers in this spectral region. The idler pulse energy is ~3.4 μJ with a conversion efficiency of ~2% from the ~2 μm pump to the idler pulse. △ Less

Submitted 28 July, 2019; v1 submitted 7 November, 2018; originally announced November 2018.

Comments: 5 pages, 5 figures

Journal ref: Opt. Lett. 44, 1003 (2019)

Quantum Monte Carlo simulation of the chiral Heisenberg Gross-Neveu-Yukawa phase transition with a single Dirac cone

Authors: Thomas C. Lang, Andreas M. Läuchli

Abstract: We present quantum Monte Carlo simulations for the chiral Heisenberg Gross-Neveu-Yukawa quantum phase transition of relativistic fermions with $N=4$ Dirac spinor components subject to a repulsive, local four fermion interaction in 2+1$d$. Here we employ a two dimensional lattice Hamiltonian with a single, spin-degenerate Dirac cone, which exactly reproduces a linear energy-momentum relation for al… ▽ More We present quantum Monte Carlo simulations for the chiral Heisenberg Gross-Neveu-Yukawa quantum phase transition of relativistic fermions with $N=4$ Dirac spinor components subject to a repulsive, local four fermion interaction in 2+1$d$. Here we employ a two dimensional lattice Hamiltonian with a single, spin-degenerate Dirac cone, which exactly reproduces a linear energy-momentum relation for all finite size lattice momenta in the absence of interactions. This allows us to significantly reduce finite size corrections compared to the widely studied honeycomb and $π$-flux lattices. A Hubbard term dynamically generates a mass beyond a critical coupling of ${U_c = 6.76(1)}$ as the system acquires antiferromagnetic order and SU(2) spin rotational symmetry is spontaneously broken. At the quantum phase transition we extract a self-consistent set of critical exponents ${ν= 0.98(1)}$, ${η_φ = 0.53(1)}$, ${η_ψ = 0.18(1)}$, ${β= 0.75(1)}$. We provide evidence for the continuous degradation of the quasi-particle weight of the fermionic excitations as the critical point is approached from the semimetallic phase. Finally we study the effective "speed of light" of the low-energy relativistic description, which depends on the interaction $U$, but is expected to be regular across the quantum phase transition. We illustrate that the strongly coupled bosonic and fermionic excitations share a common velocity at the critical point. △ Less

Submitted 15 August, 2018; v1 submitted 3 August, 2018; originally announced August 2018.

Comments: 6+6 pages, 12 figures

Journal ref: Phys. Rev. Lett. 123, 137602 (2019)

Hamiltonian Renormalisation II. Renormalisation Flow of 1+1 dimensional free scalar fields: Derivation

Authors: Thorsten Lang, Klaus Liegener, Thomas Thiemann

Abstract: In the companion paper we motivated a renormalisation flow on Osterwalder-Schrader data (OS-data) consisting of 1. a Hilbert space, 2. a cyclic vacuum and 3. a Hamiltonian annihilating that vacuum. As the name suggests, the motivation was via the OS reconstruction theorem which allows to reconstruct the OS data from an OS measure satisfying (a subset of) the OS axioms, in particular reflection pos… ▽ More In the companion paper we motivated a renormalisation flow on Osterwalder-Schrader data (OS-data) consisting of 1. a Hilbert space, 2. a cyclic vacuum and 3. a Hamiltonian annihilating that vacuum. As the name suggests, the motivation was via the OS reconstruction theorem which allows to reconstruct the OS data from an OS measure satisfying (a subset of) the OS axioms, in particular reflection positivity. The guiding principle was to map the usual Wilsonian path integral renormalisation flow onto a flow of the corresponding OS data. We showed that this induced flow on the OS data has an unwanted feature which disqualifies the associated coarse grained Hamiltonians from being the projections of a continuum Hamiltonian onto vectors in the coarse grained Hilbert space. This motivated the definition of a direct Hamiltonian renormalisation flow which follows the guiding principle but does not suffer from the afore mentioned caveat. In order to test our proposal, we apply it to the only known completely solvable model, namely the case of free scalar quantum fields. In this paper we focus on the Klein Gordon field in two spacetime dimensions and illustrate the difference between the path integral induced and direct Hamiltonian flow. Generalisations to more general models in higher dimensions will be discussed in our companion papers. △ Less

Submitted 6 July, 2019; v1 submitted 16 November, 2017; originally announced November 2017.

Comments: 30 pages

Journal ref: Class.Quant.Grav. 35 (2018) no.24, 245012

Hamiltonian Renormalisation IV. Renormalisation Flow of D+1 dimensional free scalar fields and Rotation Invariance

Authors: Thorsten Lang, Klaus Liegener, Thomas Thiemann

Abstract: In this article we extend the test of Hamiltonian Renormalisation proposed in this series of articles to the D-dimensional case using a massive free scalar field. The concepts we introduce are explicitly computed for the D=2 case but transfer immediately to higher dimensions. In this article we define and verify a criterion that monitors, at finite resolution defined by a cubic lattice, whether th… ▽ More In this article we extend the test of Hamiltonian Renormalisation proposed in this series of articles to the D-dimensional case using a massive free scalar field. The concepts we introduce are explicitly computed for the D=2 case but transfer immediately to higher dimensions. In this article we define and verify a criterion that monitors, at finite resolution defined by a cubic lattice, whether the flow approaches a rotationally invariant fixed point. △ Less

Submitted 6 July, 2019; v1 submitted 15 November, 2017; originally announced November 2017.

Comments: 28 pages

Journal ref: Class.Quant.Grav. 35 (2018) no.24, 245014

Hamiltonian Renormalization III. Renormalisation Flow of 1+1 dimensional free scalar fields: Properties

Authors: Thorsten Lang, Klaus Liegener, Thomas Thiemann

Abstract: This is the third paper in a series of four in which a renormalisation flow is introduced which acts directly on the Osterwalder-Schrader data (OS data) without recourse to a path integral. Here the OS data consist of a Hilbert space, a cyclic vacuum vector therein and a Hamiltonian annihilating the vacuum which can be obtained from an OS measure, that is a measure respecting (a subset of) the OS… ▽ More This is the third paper in a series of four in which a renormalisation flow is introduced which acts directly on the Osterwalder-Schrader data (OS data) without recourse to a path integral. Here the OS data consist of a Hilbert space, a cyclic vacuum vector therein and a Hamiltonian annihilating the vacuum which can be obtained from an OS measure, that is a measure respecting (a subset of) the OS axioms. In the previous paper we successfully tested our proposal for the two-dimensional massive Klein-Gordon model, that is, we could confirm that our framework finds the correct fixed point starting from a natural initial naive discretisation of the finite resolution Hamiltonians, in particular the underlying Laplacian on a lattice, and a natural coarse graining map that drives the renormalisation flow. However, several questions remained unanswered. How generic can the initial discretisation and the coarse graining map be in order that the fixed point is not changed or is at least not lost, in other words, how universal is the fixed point structure? Is the fix point in fact stable, that is, is the fixed point actually a limit of the renormalisation sequence? We will address these questions in the present paper. △ Less

Submitted 6 July, 2019; v1 submitted 15 November, 2017; originally announced November 2017.

Comments: 27 pages

Journal ref: Class.Quant.Grav. 35 (2018) no.24, 245013

Hamiltonian Renormalisation I: Derivation from Osterwalder-Schrader Reconstruction

Authors: Thorsten Lang, Klaus Liegener, Thomas Thiemann

Abstract: A possible avenue towards a non-perturbative Quantum Field Theory (QFT) on Minkowski space is the constructive approach which employs the Euclidian path integral formulation, in the presence of both ultraviolet (UV) and infrared (IR) regulators, as starting point. The UV regulator is to be taken away by renormalisation group techniques which in case of success leads to a measure on the space of ge… ▽ More A possible avenue towards a non-perturbative Quantum Field Theory (QFT) on Minkowski space is the constructive approach which employs the Euclidian path integral formulation, in the presence of both ultraviolet (UV) and infrared (IR) regulators, as starting point. The UV regulator is to be taken away by renormalisation group techniques which in case of success leads to a measure on the space of generalised Euclidian fields in finite volume. The IR regulator corresponds to the thermodynamic limit of the system in the statistical physics sense. If the resulting measure obeys the Osterwalder-Schrader axioms, the actual QFT on Minkowski space is then obtained by Osterwalder-Schrader reconstruction. In this work we study the question whether it is possible to reformulate the renormalisation group non-perturbatively directly at the operator (Hamiltonian) level. Hamiltonian renormalisation would be the natural route to follow if one had easier access to an interacting Hamiltonian operator rather than to a path integral, at least in the presence of UV and/or IR cut-off, which is generically the case in complicated gauge theories such as General Relativity. Our guiding principle for the definition of the direct Hamiltonian renormalisation group is that it results in the same continuum theory as the covariant (path integral) renormalisation group. In order to achieve this, we invert the Osterwalder-Schrader reconstruction, which may be called Osterwalder-Schrader construction of a Wiener measure from the underlying Hamiltonian theory. The resulting correspondence between reflection positive measures and Osterwalder-Schrader data consisting of a Hilbert space, a Hamiltonian and a ground state vector allows us to monitor the effect of the renormalisation flow of measures in terms of their Osterwalder-Schrader data which motivates a natural direct Hamiltonian renormalisation scheme. △ Less

Submitted 6 July, 2019; v1 submitted 15 November, 2017; originally announced November 2017.

Comments: 36 pages

Journal ref: Class.Quant.Grav. 35 (2018) no.24, 245011

Spontaneous particle-hole symmetry breaking of correlated fermions on the Lieb lattice

Authors: Martin Bercx, Johannes S. Hofmann, Fakher F. Assaad, Thomas C. Lang

Abstract: We study spinless fermions with nearest-neighbor repulsive interactions ($t$-$V$ model) on the two-dimensional three-band Lieb lattice. At half-filling, the free electronic band structure consists of a flat band at zero energy and a single cone with linear dispersion. The flat band is expected to be unstable upon inclusion of electronic correlations, and a natural channel is charge order. However,… ▽ More We study spinless fermions with nearest-neighbor repulsive interactions ($t$-$V$ model) on the two-dimensional three-band Lieb lattice. At half-filling, the free electronic band structure consists of a flat band at zero energy and a single cone with linear dispersion. The flat band is expected to be unstable upon inclusion of electronic correlations, and a natural channel is charge order. However, due to the three-orbital unit cell, commensurate charge order implies an imbalance of electron and hole densities and therefore doping away from half-filling. Our numerical results show that below a finite-temperature Ising transition a charge density wave with one electron and two holes per unit cell and its partner under particle-hole transformation are spontaneously generated. Our calculations are based on recent advances in auxiliary-field and continuous-time quantum Monte Carlo simulations that allow sign-free simulations of spinless fermions at half-filling. It is argued that particle-hole symmetry breaking provides a route to access levels of finite doping, without introducing a sign problem. △ Less

Submitted 9 January, 2017; v1 submitted 11 October, 2016; originally announced October 2016.

Comments: 9 pages, 6 figures, added data for strong Coulomb repulsion and classical Ising-limit

Journal ref: Phys. Rev. B 95, 035108 (2017)

Interaction induced Dirac fermions from quadratic band touching in bilayer graphene

Authors: Sumiran Pujari, Thomas C. Lang, Ganpathy Murthy, Ribhu K. Kaul

Abstract: We revisit the effect of local interactions on the quadratic band touching (QBT) of Bernal stacked bilayer graphene models using renormalization group (RG) arguments and quantum Monte Carlo simulations of the Hubbard model. We present an RG argument which predicts, contrary to previous studies, that weak interactions do not flow to strong coupling even if the free dispersion has a QBT. Instead the… ▽ More We revisit the effect of local interactions on the quadratic band touching (QBT) of Bernal stacked bilayer graphene models using renormalization group (RG) arguments and quantum Monte Carlo simulations of the Hubbard model. We present an RG argument which predicts, contrary to previous studies, that weak interactions do not flow to strong coupling even if the free dispersion has a QBT. Instead they generate a linear term in the dispersion, which causes the interactions to flow back to weak coupling. Consistent with this RG scenario, in unbiased quantum Monte Carlo simulations of the Hubbard model we find compelling evidence that antiferromagnetism turns on at a finite $U/t$, despite the $U=0$ hopping problem having a QBT. The onset of antiferromagnetism takes place at a continuous transition which is consistent with a dynamical critical exponent $z=1$ as expected for 2+1 d Gross-Neveu criticality. We conclude that generically in models of bilayer graphene, even if the free dispersion has a QBT, small local interactions generate a Dirac phase with no symmetry breaking and that there is a finite-coupling transition out of this phase to a symmetry-broken state. △ Less

Submitted 21 August, 2016; v1 submitted 13 April, 2016; originally announced April 2016.

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

Pellet refuelling of particle loss due to ELM mitigation with RMPs in the ASDEX Upgrade tokamak at low collisionality

Authors: M Valovič, P T Lang, A Kirk, W Suttrop, M Cavedon, L R Fischer, L Garzotti, L Guimarais, G Kocsis, G Cseh, B Plőckl, T Szepesi, A Thornton, A Mlynek, G Tardini, E Viezzer, R Scannell, E Wolfrum, the ASDEX Upgrade team, the EUROfusion MST1 team

Abstract: The complete refuelling of the plasma density loss (pump-out) caused by mitigation of Edge Localised Modes (ELMs) is demonstrated on the ASDEX Upgrade tokamak. The plasma is refuelled by injection of frozen deuterium pellets and ELMs are mitigated by external resonant magnetic perturbations (RMPs). In this experiment relevant dimensionless parameters, such as relative pellet size, relative RMP amp… ▽ More The complete refuelling of the plasma density loss (pump-out) caused by mitigation of Edge Localised Modes (ELMs) is demonstrated on the ASDEX Upgrade tokamak. The plasma is refuelled by injection of frozen deuterium pellets and ELMs are mitigated by external resonant magnetic perturbations (RMPs). In this experiment relevant dimensionless parameters, such as relative pellet size, relative RMP amplitude and pedestal collisionality are kept at the ITER like values. Refuelling of density pump out requires a factor of two increase of nominal fuelling rate. Energy confinement and pedestal temperatures are not restored to pre-RMP values by pellet refuelling. △ Less

Submitted 7 October, 2015; originally announced October 2015.

Comments: 10 pages, 5 figures. This is an author-created, un-copyedited version of an article submitted for publication in Nuclear Fusion. IoP Publishing Ltd is not responsible for any errors or omissions in this version of the manuscript or any version derived from it

Journal ref: 2016 Nucl. Fusion 56 066009

Explorations in Statistics Research: An Approach to Expose Undergraduates to Authentic Data Analysis

Authors: Deborah Nolan, Duncan Temple Lang

Abstract: The Explorations in Statistics Research workshop is a one-week NSF-funded summer program that introduces undergraduate students to current research problems in applied statistics. The goal of the workshop is to expose students to exciting, modern applied statistical research and practice, with the ultimate aim of interesting them in seeking more training in statistics at the undergraduate and grad… ▽ More The Explorations in Statistics Research workshop is a one-week NSF-funded summer program that introduces undergraduate students to current research problems in applied statistics. The goal of the workshop is to expose students to exciting, modern applied statistical research and practice, with the ultimate aim of interesting them in seeking more training in statistics at the undergraduate and graduate levels. The program is explicitly designed to engage students in the connections between authentic domain problems and the statistical ideas and approaches needed to address these problems, which is an important aspect of statistical thinking that is difficult to teach and sometimes lacking in our methodological courses and programs. Over the past nine years, we ran the workshop six times and a similar program in the sciences two times. We describe the program, summarize feedback from participants, and identify the key features to its success. We abstract these features and provide a set of recommendations for how faculty can incorporate important elements into their regular courses. △ Less

Submitted 22 August, 2015; originally announced August 2015.

Programming with models: writing statistical algorithms for general model structures with NIMBLE

Authors: Perry de Valpine, Daniel Turek, Christopher J. Paciorek, Clifford Anderson-Bergman, Duncan Temple Lang, Rastislav Bodik

Abstract: We describe NIMBLE, a system for programming statistical algorithms for general model structures within R. NIMBLE is designed to meet three challenges: flexible model specification, a language for programming algorithms that can use different models, and a balance between high-level programmability and execution efficiency. For model specification, NIMBLE extends the BUGS language and creates mode… ▽ More We describe NIMBLE, a system for programming statistical algorithms for general model structures within R. NIMBLE is designed to meet three challenges: flexible model specification, a language for programming algorithms that can use different models, and a balance between high-level programmability and execution efficiency. For model specification, NIMBLE extends the BUGS language and creates model objects, which can manipulate variables, calculate log probability values, generate simulations, and query the relationships among variables. For algorithm programming, NIMBLE provides functions that operate with model objects using two stages of evaluation. The first stage allows specialization of a function to a particular model and/or nodes, such as creating a Metropolis-Hastings sampler for a particular block of nodes. The second stage allows repeated execution of computations using the results of the first stage. To achieve efficient second-stage computation, NIMBLE compiles models and functions via C++, using the Eigen library for linear algebra, and provides the user with an interface to compiled objects. The NIMBLE language represents a compilable domain-specific language (DSL) embedded within R. This paper provides an overview of the design and rationale for NIMBLE along with illustrative examples including importance sampling, Markov chain Monte Carlo (MCMC) and Monte Carlo expectation maximization (MCEM). △ Less

Submitted 12 April, 2016; v1 submitted 19 May, 2015; originally announced May 2015.

Comments: 20 pages, 2 figures

Journal ref: Journal of Computational and Graphical Statistics (2017) 26: 403-413

Enhancing R with Advanced Compilation Tools and Methods

Authors: Duncan Temple Lang

Abstract: I describe an approach to compiling common idioms in R code directly to native machine code and illustrate it with several examples. Not only can this yield significant performance gains, but it allows us to use new approaches to computing in R. Importantly, the compilation requires no changes to R itself, but is done entirely via R packages. This allows others to experiment with different compila… ▽ More I describe an approach to compiling common idioms in R code directly to native machine code and illustrate it with several examples. Not only can this yield significant performance gains, but it allows us to use new approaches to computing in R. Importantly, the compilation requires no changes to R itself, but is done entirely via R packages. This allows others to experiment with different compilation strategies and even to define new domain-specific languages within R. We use the Low-Level Virtual Machine (LLVM) compiler toolkit to create the native code and perform sophisticated optimizations on the code. By adopting this widely used software within R, we leverage its ability to generate code for different platforms such as CPUs and GPUs, and will continue to benefit from its ongoing development. This approach potentially allows us to develop high-level R code that is also fast, that can be compiled to work with different data representations and sources, and that could even be run outside of R. The approach aims to both provide a compiler for a limited subset of the R language and also to enable R programmers to write other compilers. This is another approach to help us write high-level descriptions of what we want to compute, not how. △ Less

Submitted 9 September, 2014; originally announced September 2014.

Comments: Published in at http://dx.doi.org/10.1214/13-STS462 the Statistical Science (http://www.imstat.org/sts/) by the Institute of Mathematical Statistics (http://www.imstat.org)

Report number: IMS-STS-STS462

Journal ref: Statistical Science 2014, Vol. 29, No. 2, 181-200

Strangeness in Quark Matter 2013: Opening Talk

Authors: J. Steinheimer, T. Lang, H. van Hees, A. S. Botvina, K. K. Gudima, I. N. Mishustin, H. Stöcker, M. Bleicher

Abstract: We discuss several new developments in the field of strange and heavy flavor physics in high energy heavy ion collisions. As shown by many recent theoretical works, heavy flavored particles give us a unique opportunity to study the properties of systems created in these collisions. Two in particular important aspects, the production of (multi) strange hypernuclei and the properties of heavy flavor… ▽ More We discuss several new developments in the field of strange and heavy flavor physics in high energy heavy ion collisions. As shown by many recent theoretical works, heavy flavored particles give us a unique opportunity to study the properties of systems created in these collisions. Two in particular important aspects, the production of (multi) strange hypernuclei and the properties of heavy flavor mesons, are at the core of several future facilities and will be discussed in detail. △ Less

Submitted 14 March, 2014; v1 submitted 13 March, 2014; originally announced March 2014.

Comments: 6 pages 8 figures, Opening talk of the Strangeness in Quark Matter conference 2013

Planning with Noisy Probabilistic Relational Rules

Authors: Tobias Lang, Marc Toussaint

Abstract: Noisy probabilistic relational rules are a promising world model representation for several reasons. They are compact and generalize over world instantiations. They are usually interpretable and they can be learned effectively from the action experiences in complex worlds. We investigate reasoning with such rules in grounded relational domains. Our algorithms exploit the compactness of rules for e… ▽ More Noisy probabilistic relational rules are a promising world model representation for several reasons. They are compact and generalize over world instantiations. They are usually interpretable and they can be learned effectively from the action experiences in complex worlds. We investigate reasoning with such rules in grounded relational domains. Our algorithms exploit the compactness of rules for efficient and flexible decision-theoretic planning. As a first approach, we combine these rules with the Upper Confidence Bounds applied to Trees (UCT) algorithm based on look-ahead trees. Our second approach converts these rules into a structured dynamic Bayesian network representation and predicts the effects of action sequences using approximate inference and beliefs over world states. We evaluate the effectiveness of our approaches for planning in a simulated complex 3D robot manipulation scenario with an articulated manipulator and realistic physics and in domains of the probabilistic planning competition. Empirical results show that our methods can solve problems where existing methods fail. △ Less

Submitted 16 January, 2014; originally announced January 2014.

Journal ref: Journal Of Artificial Intelligence Research, Volume 39, pages 1-49, 2010

Entanglement Spectra of Interacting Fermions in Quantum Monte Carlo Simulations

Authors: Fakher F. Assaad, Thomas C. Lang, Francesco Parisen Toldin

Abstract: In a recent article T. Grover [Phys. Rev. Lett. 111, 130402 (2013)] introduced a simple method to compute Renyi entanglement entropies in the realm of the auxiliary field quantum Monte Carlo algorithm. Here, we further develop this approach and provide a stabilization scheme to compute higher order Renyi entropies and an extension to access the entanglement spectrum. The method is tested on system… ▽ More In a recent article T. Grover [Phys. Rev. Lett. 111, 130402 (2013)] introduced a simple method to compute Renyi entanglement entropies in the realm of the auxiliary field quantum Monte Carlo algorithm. Here, we further develop this approach and provide a stabilization scheme to compute higher order Renyi entropies and an extension to access the entanglement spectrum. The method is tested on systems of correlated topological insulators. △ Less

Submitted 25 March, 2014; v1 submitted 22 November, 2013; originally announced November 2013.

Comments: 7+ pages, 5 figures

Journal ref: Phys. Rev. B 89, 125121 (2014)

The characterization of topological properties in Quantum Monte Carlo simulations of the Kane-Mele-Hubbard model

Authors: Zi Yang Meng, Hsiang-Hsuan Hung, Thomas C. Lang

Abstract: Topological insulators present a bulk gap, but allow for dissipationless spin transport along the edges. These exotic states are characterized by the $Z_2$ topological invariant and are protected by time-reversal symmetry. The Kane-Mele model is one model to realize this topological class in two dimensions, also called the quantum spin Hall state. In this review, we provide a pedagogical introduct… ▽ More Topological insulators present a bulk gap, but allow for dissipationless spin transport along the edges. These exotic states are characterized by the $Z_2$ topological invariant and are protected by time-reversal symmetry. The Kane-Mele model is one model to realize this topological class in two dimensions, also called the quantum spin Hall state. In this review, we provide a pedagogical introduction to the influence of correlation effects in the quantum spin Hall states, with special focus on the half-filled Kane-Mele-Hubbard model, solved by means of unbiased determinant quantum Monte Carlo (QMC) simulations. We explain the idea of identifying the topological insulator via $π$-flux insertion, the $Z_2$ invariant and the associated behavior of the zero-frequency Green's function, as well as the spin Chern number in parameter-driven topological phase transitions. The examples considered are two descendants of the Kane-Mele-Hubbard model, the generalized and dimerized Kane-Mele-Hubbard model. From the $Z_2$ index, spin Chern numbers and the Green's function behavior, one can observe that correlation effects induce shifts of the topological phase boundaries. Although the implementation of these topological quantities has been successfully employed in QMC simulations to describe the topological phase transition, we also point out their limitations as well as suggest possible future directions in using numerical methods to characterize topological properties of strongly correlated condensed matter systems. △ Less

Submitted 7 December, 2013; v1 submitted 22 October, 2013; originally announced October 2013.

Comments: 35 pages, 16 figures, brief review for Mod. Phys. Lett B

Journal ref: Mod. Phys. Lett B, Vol 28, No 1 (2014) 143001

Dimerized Solids and Resonating Plaquette Order in SU(N)-Dirac Fermions

Authors: Thomas C. Lang, Zi Yang Meng, Alejandro Muramatsu, Stefan Wessel, Fakher F. Assaad

Abstract: We study the quantum phases of fermions with an explicit SU(N)-symmetric, Heisenberg-like nearest-neighbor flavor exchange interaction on the honeycomb lattice at half-filling. Employing projective (zero temperature) quantum Monte Carlo simulations for even values of N, we explore the evolution from a weak-coupling semimetal into the strong-coupling, insulating regime. Furthermore, we compare our… ▽ More We study the quantum phases of fermions with an explicit SU(N)-symmetric, Heisenberg-like nearest-neighbor flavor exchange interaction on the honeycomb lattice at half-filling. Employing projective (zero temperature) quantum Monte Carlo simulations for even values of N, we explore the evolution from a weak-coupling semimetal into the strong-coupling, insulating regime. Furthermore, we compare our numerical results to a saddle-point approximation in the large-N limit. From the large-N regime down to the SU(6) case, the insulating state is found to be a columnar valence bond crystal, with a direct transition to the semimetal at weak, finite coupling, in agreement with the mean-field result in the large-N limit. At SU(4) however, the insulator exhibits a subtly different valence bond crystal structure, stabilized by resonating valence bond plaquettes. In the SU(2) limit, our results support a direct transition between the semimetal and an antiferromagnetic insulator. △ Less

Submitted 12 August, 2013; v1 submitted 13 June, 2013; originally announced June 2013.

Comments: 5 pages, 6 figures

Journal ref: Phys. Rev. Lett. 111, 066401 (2013)

Correlated D-meson decays competing against thermal QGP dilepton radiation

Authors: Thomas Lang, Hendrik van Hees, Jan Steinheimer, Marcus Bleicher

Abstract: The QGP that might be created in ultrarelativistic heavy-ion collisions is expected to radiate thermal dilepton radiation. However, this thermal dilepton radiation interferes with dileptons originating from hadron decays. In the invariant mass region between the $φ$ and $J/Ψ$ peak ($1\,$GeV$\lesssim M_{\ell^+ \ell^-} \lesssim 3 \,$GeV) the most substantial background of hadron decays originates fr… ▽ More The QGP that might be created in ultrarelativistic heavy-ion collisions is expected to radiate thermal dilepton radiation. However, this thermal dilepton radiation interferes with dileptons originating from hadron decays. In the invariant mass region between the $φ$ and $J/Ψ$ peak ($1\,$GeV$\lesssim M_{\ell^+ \ell^-} \lesssim 3 \,$GeV) the most substantial background of hadron decays originates from correlated D$\bar{\mathrm{D}}$-meson decays. We evaluate this background using a Langevin simulation for charm quarks. As background medium we utilize the well-tested UrQMD-hybrid model. The required drag and diffusion coefficients are taken from a resonance approach. The decoupling of the charm quarks from the hot medium is performed at a temperature of $130\,$MeV and as hadronization mechanism a coalescence approach is chosen. This model for charm quark interactions with the medium has already been successfully applied to the study of the medium modification and the elliptic flow at FAIR, RHIC and LHC energies. In this proceeding we present our results for the dilepton radiation from correlated D$\bar{\mathrm{D}}$ decays at RHIC energy in comparison to PHENIX measurements in the invariant mass range between 1 and 3 GeV using different interaction scenarios. These results can be utilized to estimate the thermal QGP radiation. △ Less

Submitted 12 June, 2013; originally announced June 2013.

Comments: 4 pages, 1 figure