Search

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., Comment: 8 pages, 5 figures, 1 table. Short summary of a full paper which can be found here: arXiv:2307.13293

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., Comment: 11 pages, 14 figures, 5 tables. Replaced with published version. Reproduced with permission from Astronomy & Astrophysics, $\copyright$ ESO

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., Comment: 12 pages, 8 figures

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., Comment: 13 pages, 3 figures, 2 tables

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., Comment: 10 pages, 27 figures, 3 tables

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, Comment: 27 pages, 4 figures

Twisted double- and mono-bilayer graphene are graphene-based moir\'e 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\'e 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., Comment: 35 pages, 12 figures

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., Comment: Main text: 6 pages, 4 figures. Supplementary information: 10 pages, 16 figures

We perform an extensive exact diagonalization study of interaction driven insulators in spin- and valley-polarized moir\'{e} 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\'{e} heterostructure and affirms the relevance of our results beyond the scope of graphene based materials., Comment: 16 pages, 18 figures

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., Comment: Submitted: October 24, 2018, Accepted: November 4, 2019

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., Comment: 9 pages, 10 figures

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 $\pi$-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 ${\epsilon = 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., Comment: 26 pages, 14 figures

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 $\pi$-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 ${\nu = 0.98(1)}$, ${\eta_{\phi} = 0.53(1)}$, ${\eta_{\psi} = 0.18(1)}$, ${\beta = 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., Comment: 6+6 pages, 12 figures