Loop formulation of supersymmetric Yang-Mills quantum mechanics

We derive the fermion loop formulation of N=4 supersymmetric SU(N) Yang-Mills quantum mechanics on the lattice. The loop formulation naturally separates the contributions to the partition function into its bosonic and fermionic parts with fixed fermion number and provides a way to control potential fermion sign problems arising in numerical simulations of the theory. Furthermore, we present a reduced fermion matrix determinant which allows the projection into the canonical sectors of the theory and hence constitutes an alternative approach to simulate the theory on the lattice.

Supersymmetric quantum mechanics on the lattice: I. Loop formulation

Simulations of supersymmetric field theories on the lattice with (spontaneously) broken supersymmetry suffer from a fermion sign problem related to the vanishing of the Witten index. We propose a novel approach which solves this problem in low dimensions by formulating the path integral on the lattice in terms of fermion loops. For N=2 supersymmetric quantum mechanics the loop formulation becomes particularly simple and in this paper – the first in a series of three – we discuss in detail the reformulation of this model in terms of fermionic and bosonic bonds for various lattice discretisations including one which is Q-exact.

Development of a balanced experimental–computational approach to understanding the mechanics of proximal femur fractures

The majority of people who sustain hip fractures after a fall to the side would not have been identified using current screening techniques such as areal bone mineral density. Identifying them, however, is essential so that appropriate pharmacological or lifestyle interventions can be implemented. A protocol, demonstrated on a single specimen, is introduced, comprising the following components; in vitro biofidelic drop tower testing of a proximal femur; high-speed image analysis through digital image correlation; detailed accounting of the energy present during the drop tower test; organ level finite element simulations of the drop tower test; micro level finite element simulations of critical volumes of interest in the trabecular bone. Fracture in the femoral specimen initiated in the superior part of the neck. Measured fracture load was 3760 N, compared to 4871 N predicted based on the finite element analysis. Digital image correlation showed compressive surface strains as high as 7.1% prior to fracture. Voxel level results were consistent with high-speed video data and helped identify hidden local structural weaknesses. We found using a drop tower test protocol that a femoral neck fracture can be created with a fall velocity and energy representative of a sideways fall from standing. Additionally, we found that the nested explicit finite element method used allowed us to identify local structural weaknesses associated with femur fracture initiation.

Minor inheritance inhibits the calibration of the 10Be production rate from the AD 1717 Val Ferret rock avalanche, European Alps

To calibrate the in situ 10Be production rate, we collected surface samples from nine large granitic boulders within the deposits of a rock avalanche that occurred in AD 1717 in the upper Ferret Valley, Mont Blanc Massif, Italy. The 10Be concentrations were extremely low and successfully measured within 10% analytical uncertainty or less. The concentrations vary from 4829 ± 448 to 5917 ± 476 at g−1. Using the historical age exposure time, we calculated the local and sea level-high latitude (i.e. ≥60°) cosmogenic 10Be spallogenic production rates. Depending on the scaling schemes, these vary between 4.60 ± 0.38 and 5.26 ± 0.43 at g−1 a−1. Although they correlate well with global values, our production rates are clearly higher than those from more recent calibration sites. We conclude that our 10Be production rate is a mean and an upper bound for production rates in the Massif region over the past 300 years. This rate is probably influenced by inheritance and will yield inaccurate (e.g. too young) exposure ages when applied to surface-exposure studies in the area. Other independently dated rock-avalanche deposits in the region that are approximately 103 years old could be considered as possible calibration sites.

Benchmark problems for reactive transport modeling of the generation and attenuation of acid rock drainage

Acid rock drainage (ARD) is a problem of international relevance with substantial environmental and economic implications. Reactive transport modeling has proven a powerful tool for the process-based assessment of metal release and attenuation at ARD sites. Although a variety of models has been used to investigate ARD, a systematic model intercomparison has not been conducted to date. This contribution presents such a model intercomparison involving three synthetic benchmark problems designed to evaluate model results for the most relevant processes at ARD sites. The first benchmark (ARD-B1) focuses on the oxidation of sulfide minerals in an unsaturated tailing impoundment, affected by the ingress of atmospheric oxygen. ARD-B2 extends the first problem to include pH buffering by primary mineral dissolution and secondary mineral precipitation. The third problem (ARD-B3) in addition considers the kinetic and pH-dependent dissolution of silicate minerals under low pH conditions. The set of benchmarks was solved by four reactive transport codes, namely CrunchFlow, Flotran, HP1, and MIN3P. The results comparison focused on spatial profiles of dissolved concentrations, pH and pE, pore gas composition, and mineral assemblages. In addition, results of transient profiles for selected elements and cumulative mass loadings were considered in the intercomparison. Despite substantial differences in model formulations, very good agreement was obtained between the various codes. Residual deviations between the results are analyzed and discussed in terms of their implications for capturing system evolution and long-term mass loading predictions.

Supersymmetric quantum mechanics on the lattice: II. Exact results

Simulations of supersymmetric field theories with spontaneously broken supersymmetry require in addition to the ultraviolet regularisation also an infrared one, due to the emergence of the massless Goldstino. The intricate interplay between ultraviolet and infrared effects towards the continuum and infinite volume limit demands careful investigations to avoid potential problems. In this paper – the second in a series of three – we present such an investigation for N=2 supersymmetric quantum mechanics formulated on the lattice in terms of bosonic and fermionic bonds. In one dimension, the bond formulation allows to solve the system exactly, even at finite lattice spacing, through the construction and analysis of transfer matrices. In the present paper we elaborate on this approach and discuss a range of exact results for observables such as the Witten index, the mass spectra and Ward identities.

Supersymmetric quantum mechanics on the lattice: III. Simulations and algorithms

In the fermion loop formulation the contributions to the partition function naturally separate into topological equivalence classes with a definite sign. This separation forms the basis for an efficient fermion simulation algorithm using a fluctuating open fermion string. It guarantees sufficient tunnelling between the topological sectors, and hence provides a solution to the fermion sign problem affecting systems with broken supersymmetry. Moreover, the algorithm shows no critical slowing down even in the massless limit and can hence handle the massless Goldstino mode emerging in the supersymmetry broken phase. In this paper – the third in a series of three – we present the details of the simulation algorithm and demonstrate its efficiency by means of a few examples.

Local adaptations in populations of rock-dwelling haplochromines (Pisces:Cichlidae) from southern Lake Victoria

The radiation of cichlid fishes in the African great lakes is often described as adaptive, because, at a superficial level, cichlid fishes seem adapted to the ecological niches they occupy. However, adaptiveness has rarely been studied. We investigated to what extent island populations of three species of the rock-dwelling genus Neochromis, endemic to Lake Victoria, are adapted anatomically to exploit locally abundant resources. Specifically, we asked whether different resource environments were reflected in differences in the feeding apparatus, both within species and between species. In populations of two specialized biters, the algae scrapers N. rufocaudalis and N. omnicaeruleus, the biting force of the lower jaw increased with increasing amount of items that require biting in the diet. N. greenwoodi is a less specialized biter; we found differences between two populations in the hyoid position and in the premaxilla that enhance suction feeding. These adaptations were related to the amount of items requiring suction. Comparing across three sympatric pairs of species, in each case different diets were reflected in differences in anatomy.

GRAIL gravity field determination using the Celestial Mechanics Approach

The NASA mission GRAIL (Gravity Recovery and Interior Laboratory) inherited its concept from the GRACE (Gravity Recovery and Climate Experiment) mission to determine the gravity field of the Moon. We present lunar gravity fields based on the data of GRAIL’s primary mission phase. Gravity field recovery is realized in the framework of the Celestial Mechanics Approach, using a development version of the Bernese GNSS Software along with Ka-band range-rate data series as observations and the GNI1B positions provided by NASA JPL as pseudo-observations. By comparing our results with the official level-2 GRAIL gravity field models we show that the lunar gravity field can be recovered with a high quality by adapting the Celestial Mechanics Approach, even when using pre-GRAIL gravity field models as a priori fields and when replacing sophisticated models of non-gravitational accelerations by appropriately spaced pseudo-stochastic pulses (i.e., instantaneous velocity changes). We present and evaluate two lunar gravity field solutions up to degree and order 200 – AIUB-GRL200A and AIUB-GRL200B. While the first solution uses no gravity field information beyond degree 200, the second is obtained by using the official GRAIL field GRGM900C up to degree and order 660 as a priori information. This reduces the omission errors and demonstrates the potential quality of our solution if we resolved the gravity field to higher degree.