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.

Monte Carlo modeling of polarized light propagation: Stokes vs Jones - Part I

This bipartite comparative study aims at inspecting the similarities and differences between the Jones and Stokes–Mueller formalisms when modeling polarized light propagation with numerical simulations of the Monte Carlo type. In this first part, we review the theoretical concepts that concern light propagation and detection with both pure and partially/totally unpolarized states. The latter case involving fluctuations, or “depolarizing effects,” is of special interest here: Jones and Stokes–Mueller are equally apt to model such effects and are expected to yield identical results. In a second, ensuing paper, empirical evidence is provided by means of numerical experiments, using both formalisms.

The Utilization of Classical Spin Monte Carlo Methods to Simulate the Magnetic Behavior of Extended Three-Dimensional Cubic Networks Incorporating M(II) Ions with an S = 5/2 Ground State Spin

The numerical simulations of the magnetic properties of extended three-dimensional networks containing M(II) ions with an S = 5/2 ground-state spin have been carried out within the framework of the isotropic Heisenberg model. Analytical expressions fitting the numerical simulations for the primitive cubic, diamond, together with (10−3) cubic networks have all been derived. With these empirical formulas in hands, we can now extract the interaction between the magnetic ions from the experimental data for these networks. In the case of the primitive cubic network, these expressions are directly compared with those from the high-temperature expansions of the partition function. A fit of the experimental data for three complexes, namely [(N(CH3)4][Mn(N3)] 1, [Mn(CN4)]n 2, and [FeII(bipy)3][MnII2(ox)3] 3, has been carried out. The best fits were those obtained using the following parameters, J = −3.5 cm-1, g = 2.01 (1); J = −8.3 cm-1, g = 1.95 (2); and J = −2.0 cm-1, g = 1.95 (3).

Stabilization of subcritical bypass transition in the leading-edge boundary layer by suction

This work investigates the subcritical spatial transition in the swept Hiemenz boundary layer by means of direct numerical simulations (DNS). A pair of steady co-rotating vortices located at the attachment line is enforced as a primary disturbance leading to streaks which are stable. A small secondary, time-dependent disturbance interacts with these streaks such that instability and breakdown to turbulence may occur. The instability only occurs for a certain band of secondary disturbance frequencies. Positive secondary instability growth rates could be observed for Reynolds numbers as low as , whereas the linear critical Reynolds number is. Uniform wall suction is shown to stabilise this transition mechanism, analogously to results from linear stability theory. The effects of suction on the formation of primary streaks and on the secondary growth rate are decoupled. For streaks of different suction whose amplitude is held constant by adjusting the Reynolds number, the suction is shown to increase the growth rate of the secondary instability. The stabilising influence of wall suction consists in decreasing the streak amplitude only. Depending on the Reynolds number and the suction strength, breakdown may either occur locally and may be convected along the far-field streamlines, or occur globally and cover broad regions in the downstream direction.

Reconstructing 4000 years of mass movement and tsunami history in a deep peri-Alpine lake (Lake Geneva, France-Switzerland)

The study of mass movements in lake sediments provides insights into past natural hazards at historic and prehistoric timescales. Sediments from the deep basin of Lake Geneva reveal a succession of six large-scale (volumes of 22 × 106 to 250 × 106 m3) mass-transport deposits, associated with five mass-movement events within 2600 years (4000 cal bp to 563 ad). The mass-transport deposits result from: (i) lateral slope failures (mass-transport deposit B at 3895 ± 225 cal bp and mass-transport deposits A and C at 3683 ± 128 cal bp); and (ii) Rhône delta collapses (mass-transport deposits D to G dated at 2650 ± 150 cal bp, 2185 ± 85 cal bp, 1920 ± 120 cal bp and 563 ad, respectively). Mass-transport deposits A and C were most probably triggered by an earthquake, whereas the Rhône delta collapses were likely to be due to sediment overload with a rockfall as the external trigger (mass-transport deposit G, the Tauredunum event in 563 ad known from historical records), an earthquake (mass-transport deposit E) or unknown external triggers (mass-transport deposits D and F). Independent of their origin and trigger mechanisms, numerical simulations show that all of these recorded mass-transport deposits are large enough to have generated at least metre-scale tsunamis during mass movement initiation. Since the Tauredunum event in 563 ad, two small-scale (volumes of 1 to 2 × 106 m3) mass-transport deposits (H and I) are present in the seismic record, both of which are associated with small lateral slope failures. Mass-transport deposits H and I might be related to earthquakes in Lausanne/Geneva (possibly) 1322 ad and Aigle 1584 ad, respectively. The sedimentary record of the deep basin of Lake Geneva, in combination with the historical record, show that during the past 3695 years, at least six tsunamis were generated by mass movements, indicating that the tsunami hazard in the Lake Geneva region should not be neglected, although such events are not frequent with a recurrence time of 0·0016 yr−1.

Reconstruction of in-situ porosity and porewater compositions of low-permeability crystalline rocks: Magnitude of artefacts induced by drilling and sample recovery

Geological site characterisation programmes typically rely on drill cores for direct information on subsurface rocks. However, porosity, transport properties and porewater composition measured on drill cores can deviate from in-situ values due to two main artefacts caused by drilling and sample recovery: (1) mechanical disruption that increases porosity and (2) contamination of the porewater by drilling fluid. We investigated the effect and magnitude of these perturbations on large drill core samples (12–20 cm long, 5 cmdiameter) of high-grade, granitic gneisses obtained from 350 to 600 m depth in a borehole on Olkiluoto Island (SW Finland). The drilling fluid was traced with sodium–iodide. By combining out-diffusion experiments, gravimetry, UV-microscopy and iodide mass balance calculations, we successfully quantified the magnitudes of the artefacts: 2–6% increase in porosity relative to the bulk connected porosity and 0.9 to 8.9 vol.% contamination by drilling fluid. The spatial distribution of the drilling-induced perturbations was revealed by numerical simulations of 2D diffusion matched to the experimental data. This showed that the rims of the samples have a mechanically disrupted zone 0.04 to 0.22 cm wide, characterised by faster transport properties compared to the undisturbed centre (1.8 to 7.7 times higher pore diffusion coefficient). Chemical contamination was shown to affect an even wider zone in all samples, ranging from 0.15 to 0.60 cm, inwhich iodide enrichmentwas up to 180 mg/kgwater, compared to 0.5 mg/kgwater in the uncontaminated centre. For all samples in the present case study, it turned out that the magnitude of the artefacts caused by drilling and sample recovery is so small that no correction is required for their effects. Therefore, the standard laboratory measurements of porosity, transport properties and porewater composition can be taken as valid in-situ estimates. However, it is clear that the magnitudes strongly depend on site- and drilling-specific factors and therefore our results cannot be transferred simply to other locations. We recommend the approach presented in this study as a route to obtain reliable values in future drilling campaigns aimed at characterising in-situ bedrock properties.

Comparison of 3D kinetic and hydrodynamic models to ROSINA-COPS measurements of the neutral coma of 67P/Churyumov-Gerasimenko

67P/Churyumov-Gerasimenko (67P) is a Jupiter-family comet and the object of investigation of the European Space Agency mission Rosetta. This report presents the first full 3D simulation results of 67P’s neutral gas coma. In this study we include results from a direct simulation Monte Carlo method, a hydrodynamic code, and a purely geometric calculation which computes the total illuminated surface area on the nucleus. All models include the triangulated 3D shape model of 67P as well as realistic illumination and shadowing conditions. The basic concept is the assumption that these illumination conditions on the nucleus are the main driver for the gas activity of the comet. As a consequence, the total production rate of 67P varies as a function of solar insolation. The best agreement between the model and the data is achieved when gas fluxes on the night side are in the range of 7% to 10% of the maximum flux, accounting for contributions from the most volatile components. To validate the output of our numerical simulations we compare the results of all three models to in situ gas number density measurements from the ROSINA COPS instrument. We are able to reproduce the overall features of these local neutral number density measurements of ROSINA COPS for the time period between early August 2014 and January 1 2015 with all three models. Some details in the measurements are not reproduced and warrant further investigation and refinement of the models. However, the overall assumption that illumination conditions on the nucleus are at least an important driver of the gas activity is validated by the models. According to our simulation results we find the total production rate of 67P to be constant between August and November 2014 with a value of about 1 × 10²⁶ molecules s⁻¹.

Enabling the Next Generation of Spaceborne Quadrupole Mass Spectrometers

The quadrupole mass spectrometer (QMS) has over 30 years of spaceflight heritage in making important neutral gas and low energy ion observations. Given their geometrical constraints, these instruments are currently operated at the extreme limit of their capabilities. However, a technique called higher order auxiliary excitation provides a set of novel, robust, electronics-based solutions for improving the performance of these sensors. By driving the quadrupole rods with an additional frequency nearly twice that of the normal RF operating frequency, substantially increased abundance sensitivity, maximum attainable mass resolution, and peak stability can be achieved through operation of voltage scan lines through the center of formed upper stability islands. Such improvements are modeled using numerical simulations of ion trajectories in a quadrupole field with and without applied higher order auxiliary excitation. When compared to a traditional QMS with a mass range up to 500Da, sensors can be designed with the same precision electronics to have expected mass ranges beyond 1500Da with a power increase of less than twice that of its heritage implementations.

Collisional Formation and Modeling of Asteroid Families

In the last decade, thanks to the development of sophisticated numerical codes, major breakthroughs have been achieved in our understanding of the formation of asteroid families by catastrophic disruption of large parent bodies. In this review, we describe numerical simulations of asteroid collisions that reproduced the main properties of families, accounting for both the fragmentation of an asteroid at the time of impact and the subsequent gravitational interactions of the generated fragments. The simulations demonstrate that the catastrophic disruption of bodies larger than a few hundred meters in diameter leads to the formation of large aggregates due to gravitational reaccumulation of smaller fragments, which helps explain the presence of large members within asteroid families. Thus, for the first time, numerical simulations successfully reproduced the sizes and ejection velocities of members of representative families. Moreover, the simulations provide constraints on the family dynamical histories and on the possible internal structure of family members and their parent bodies.