Source author record

Martin Bauer

Martin Bauer appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

53works
16topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

53 published item(s)

preprint2023arXiv

Limits on the cosmic neutrino background

We present the first comprehensive discussion of constraints on the cosmic neutrino background (C$ν$B) overdensity, including theoretical, experimental and cosmological limits for a wide range of neutrino masses and temperatures. Additionally, we calculate the sensitivities of future direct and indirect relic neutrino detection experiments and compare the results with the existing constraints, extending several previous analyses by taking into account that the C$ν$B reference frame may not be aligned with that of the Earth. The Pauli exclusion principle strongly disfavours overdensities $η_ν\gg 1$ at small neutrino masses, but allows for overdensities $η_ν\lesssim 125$ at the KATRIN mass bound $m_ν \simeq 0.8\,\mathrm{eV}$. On the other hand, cosmology strongly favours $0.2 \lesssim η_ν \lesssim 3.5$ in all scenarios. We find that direct detection proposals are capable of observing the C$ν$B without a significant overdensity for neutrino masses $m_ν \gtrsim 50\,\mathrm{meV}$, but require an overdensity $η_ν \gtrsim 3\times 10^5$ outside of this range. We also demonstrate that relic neutrino detection proposals are sensitive to the helicity composition of the C$ν$B, whilst some may be able to distinguish between Dirac and Majorana neutrinos.

preprint2022arXiv

Feebly-Interacting Particles:FIPs 2020 Workshop Report

With the establishment and maturation of the experimental programs searching for new physics with sizeable couplings at the LHC, there is an increasing interest in the broader particle and astrophysics community for exploring the physics of light and feebly-interacting particles as a paradigm complementary to a New Physics sector at the TeV scale and beyond. FIPs 2020 has been the first workshop fully dedicated to the physics of feebly-interacting particles and was held virtually from 31 August to 4 September 2020. The workshop has gathered together experts from collider, beam dump, fixed target experiments, as well as from astrophysics, axions/ALPs searches, current/future neutrino experiments, and dark matter direct detection communities to discuss progress in experimental searches and underlying theory models for FIPs physics, and to enhance the cross-fertilisation across different fields. FIPs 2020 has been complemented by the topical workshop "Physics Beyond Colliders meets theory", held at CERN from 7 June to 9 June 2020. This document presents the summary of the talks presented at the workshops and the outcome of the subsequent discussions held immediately after. It aims to provide a clear picture of this blooming field and proposes a few recommendations for the next round of experimental results.

preprint2022arXiv

The Square Root Normal Field Distance and Unbalanced Optimal Transport

This paper explores a novel connection between two areas: shape analysis of surfaces and unbalanced optimal transport. Specifically, we characterize the square root normal field (SRNF) shape distance as the pullback of the Wasserstein-Fisher-Rao (WFR) unbalanced optimal transport distance. In addition, we propose a new algorithm for computing the WFR distance and present numerical results that highlight the effectiveness of this algorithm. As a consequence of our results we obtain a precise method for computing the SRNF shape distance directly on piecewise linear surfaces and gain new insights about the degeneracy of this distance.

preprint2021arXiv

Can we run to infinity? The diameter of the diffeomorphism group with respect to right-invariant Sobolev metrics

The group $\text{Diff}(\mathcal{M})$ of diffeomorphisms of a closed manifold $\mathcal{M}$ is naturally equipped with various right-invariant Sobolev norms $W^{s,p}$. Recent work showed that for sufficiently weak norms, the geodesic distance collapses completely (namely, when $sp\le \text{dim}\mathcal{M}$ and $s<1$). But when there is no collapse, what kind of metric space is obtained? In particular, does it have a finite or infinite diameter? This is the question we study in this paper. We show that the diameter is infinite for strong enough norms, when $(s-1)p\ge \text{dim}\mathcal{M}$, and that for spheres the diameter is finite when $(s-1)p<1$. In particular, this gives a full characterization of the diameter of $\text{Diff}(S^1)$. In addition, we show that for $\text{Diff}_c(\mathbb{R}^n)$, if the diameter is not zero, it is infinite.

preprint2021arXiv

Light Dark Matter Annihilation and Scattering in LHC Detectors

We systematically study models with light scalar and pseudoscalar dark matter candidates and their potential signals at the LHC. First, we derive cosmological bounds on models with the Standard Model Higgs mediator and with a new weak-scale mediator. Next, we study two processes inspired by the indirect and direct detection process topologies, now happening inside the LHC detectors. We find that LHC can observe very light dark matter over a huge mass range if it is produced in mediator decays and then scatters with the detector material to generate jets in the nuclear recoil.

preprint2020arXiv

A diffeomorphism-invariant metric on the space of vector-valued one-forms

In this article we introduce a diffeomorphism-invariant Riemannian metric on the space of vector valued one-forms. The particular choice of metric is motivated by potential future applications in the field of functional data and shape analysis and by connections to the Ebin metric on the space of all Riemannian metrics. In the present work we calculate the geodesic equations and obtain an explicit formula for the solutions to the corresponding initial value problem. Using this we show that it is a geodesically and metrically incomplete space and study the existence of totally geodesic subspaces. Furthermore, we calculate the sectional curvature and observe that, depending on the dimension of the base manifold and the target space, it either has a semidefinite sign or admits both signs.

preprint2020arXiv

A Standard-based Open Source IoT Platform: FIWARE

The ever-increasing acceleration of technology evolution in all fields is rapidly changing the architectures of data-driven systems towards the Internet-of-Things concept. Many general and specific-purpose IoT platforms are already available. This article introduces the capabilities of the FIWARE framework that is transitioning from a research to a commercial level. We base our exposition on the analysis of three real-world use cases (global IoT market, analytics in smart cities, and IoT augmented autonomous driving) and their requirements that are addressed with the usage of FIWARE. We highlight the lessons learnt during the design, implementation and deployment phases for each of the use cases and their critical issues. Finally we give two examples showing that FIWARE still maintains openness to innovation: semantics and privacy.

preprint2020arXiv

Axion-Like Particles at Future Colliders

Axion-like particles (ALPs) are pseudo Nambu-Goldstone bosons of spontaneously broken global symmetries in high-energy extensions of the Standard Model (SM). This makes them a prime target for future experiments aiming to discover new physics which addresses some of the open questions of the SM. While future high-precision experiments can discover ALPs with masses well below the GeV scale, heavier ALPs can be searched for at future high-energy lepton and hadron colliders. We discuss the reach of the different proposed colliders, focusing on resonant ALP production, ALP production in the decay of heavy SM resonances, and associate ALP production with photons, Z bosons or Higgs bosons. We consider the leading effective operators mediating interactions between the ALP and SM particles and discuss search strategies for ALPs decaying promptly as well as ALPs with delayed decays. Projections for the high-luminosity run of the LHC and its high-energy upgrade, CLIC, the future $e^+e^-$ ring-colliders CEPC and FCC-ee, the future pp colliders SPPC and FCC-hh, and for the MATHUSLA surface array are presented. We further discuss the constraining power of future measurements of electroweak precision parameters on the relevant ALP couplings.

preprint2020arXiv

Geodesic completeness of the $H^{3/2}$ metric on $\mathrm{Diff}(S^{1})$

Of concern is the study of the long-time existence of solutions to the Euler--Arnold equation of the right-invariant $H^{3/2}$-metric on the diffeomorphism group of the circle. In previous work by Escher and Kolev it has been shown that this equation admits long-time solutions if the order $s$ of the metric is greater than $3/2$, the behaviour for the critical Sobolev index $s=3/2$ has been left open. In this article we fill this gap by proving the analogous result also for the boundary case. The behaviour of the $H^{3/2}$-metric is, however, still different from its higher order counter parts, as it does not induce a complete Riemannian metric on any group of Sobolev diffeomorphisms.

preprint2020arXiv

lbmpy: Automatic code generation for efficient parallel lattice Boltzmann methods

Lattice Boltzmann methods are a popular mesoscopic alternative to macroscopic computational fluid dynamics solvers. Many variants have been developed that vary in complexity, accuracy, and computational cost. Extensions are available to simulate multi-phase, multi-component, turbulent, or non-Newtonian flows. In this work we present lbmpy, a code generation package that supports a wide variety of different methods and provides a generic development environment for new schemes as well. A high-level domain-specific language allows the user to formulate, extend and test various lattice Boltzmann schemes. The method specification is represented in a symbolic intermediate representation. Transformations that operate on this intermediate representation optimize and parallelize the method, yielding highly efficient lattice Boltzmann compute kernels not only for single- and two-relaxation-time schemes but also for multi-relaxation-time, cumulant, and entropically stabilized methods. An integration into the HPC framework waLBerla makes massively parallel, distributed simulations possible, which is demonstrated through scaling experiments on the SuperMUC-NG supercomputing system

preprint2019arXiv

Axion-like particles, lepton-flavor violation and a new explanation of $a_μ$ and $a_e$

Axion-like particles (ALPs) with lepton-flavor violating couplings can be probed in exotic muon and tau decays. The sensitivity of different experiments depends strongly on the ALP mass and its couplings to leptons and photons. For ALPs that can be resonantly produced, the sensitivity of three-body decays such as $μ\to 3e$ and $τ\to 3μ$ exceeds by many orders of magnitude that of radiative decays like $μ\to eγ$ and $τ\toμγ$. Searches for these two types of processes are therefore highly complementary. We discuss experimental constraints on ALPs with a single dominant lepton-flavor violating coupling. Allowing for one or more such couplings offers qualitatively new ways to explain the anomalies related to the magnetic moments of the muon and electron.

preprint2019arXiv

Long-Lived Particles at the Energy Frontier: The MATHUSLA Physics Case

We examine the theoretical motivations for long-lived particle (LLP) signals at the LHC in a comprehensive survey of Standard Model (SM) extensions. LLPs are a common prediction of a wide range of theories that address unsolved fundamental mysteries such as naturalness, dark matter, baryogenesis and neutrino masses, and represent a natural and generic possibility for physics beyond the SM (BSM). In most cases the LLP lifetime can be treated as a free parameter from the $μ$m scale up to the Big Bang Nucleosynthesis limit of $\sim 10^7$m. Neutral LLPs with lifetimes above $\sim$ 100m are particularly difficult to probe, as the sensitivity of the LHC main detectors is limited by challenging backgrounds, triggers, and small acceptances. MATHUSLA is a proposal for a minimally instrumented, large-volume surface detector near ATLAS or CMS. It would search for neutral LLPs produced in HL-LHC collisions by reconstructing displaced vertices (DVs) in a low-background environment, extending the sensitivity of the main detectors by orders of magnitude in the long-lifetime regime. In this white paper we study the LLP physics opportunities afforded by a MATHUSLA-like detector at the HL-LHC. We develop a model-independent approach to describe the sensitivity of MATHUSLA to BSM LLP signals, and compare it to DV and missing energy searches at ATLAS or CMS. We then explore the BSM motivations for LLPs in considerable detail, presenting a large number of new sensitivity studies. While our discussion is especially oriented towards the long-lifetime regime at MATHUSLA, this survey underlines the importance of a varied LLP search program at the LHC in general. By synthesizing these results into a general discussion of the top-down and bottom-up motivations for LLP searches, it is our aim to demonstrate the exceptional strength and breadth of the physics case for the construction of the MATHUSLA detector.

preprint2019arXiv

waLBerla: A block-structured high-performance framework for multiphysics simulations

Programming current supercomputers efficiently is a challenging task. Multiple levels of parallelism on the core, on the compute node, and between nodes need to be exploited to make full use of the system. Heterogeneous hardware architectures with accelerators further complicate the development process. waLBerla addresses these challenges by providing the user with highly efficient building blocks for developing simulations on block-structured grids. The block-structured domain partitioning is flexible enough to handle complex geometries, while the structured grid within each block allows for highly efficient implementations of stencil-based algorithms. We present several example applications realized with waLBerla, ranging from lattice Boltzmann methods to rigid particle simulations. Most importantly, these methods can be coupled together, enabling multiphysics simulations. The framework uses meta-programming techniques to generate highly efficient code for CPUs and GPUs from a symbolic method formulation. To ensure software quality and performance portability, a continuous integration toolchain automatically runs an extensive test suite encompassing multiple compilers, hardware architectures, and software configurations.

preprint2016arXiv

A Numerical Framework for Sobolev Metrics on the Space of Curves

Statistical shape analysis can be done in a Riemannian framework by endowing the set of shapes with a Riemannian metric. Sobolev metrics of order two and higher on shape spaces of parametrized or unparametrized curves have several desirable properties not present in lower order metrics, but their discretization is still largely missing. In this paper, we present algorithms to numerically solve the geodesic initial and boundary value problems for these metrics. The combination of these algorithms enables one to compute Karcher means in a Riemannian gradient-based optimization scheme and perform principal component analysis and clustering. Our framework is sufficiently general to be applicable to a wide class of metrics. We demonstrate the effectiveness of our approach by analyzing a collection of shapes representing HeLa cell nuclei.

preprint2016arXiv

Analyzing the CP Nature of a New Scalar Particle via S->Zh Decay

Scalar particles $S$ which are singlets under the Standard Model gauge group are generic features of many models of fundamental physics, in particular as possible mediators to a hidden sector. We show that the decay $S\to Zh$ provides a powerful probe of the CP nature of the scalar, because it is allowed only if $S$ has CP-odd interactions. We perform a model-independent analysis of this decay using an effective Lagrangian and compute the relevant Wilson coefficients arising from integrating out heavy fermions to one-loop order.

preprint2016arXiv

Diphoton Resonance from a Warped Extra Dimension

We argue that extensions of the Standard Model (SM) with a warped extra dimension, which successfully address the hierarchy and flavor problems of elementary particle physics, can provide an elegant explanation of the 750 GeV diphoton excess recently reported by ATLAS and CMS. A gauge-singlet bulk scalar with ${\cal O}(1)$ couplings to fermions is identified as the new resonance $S$, and the vector-like Kaluza-Klein excitations of the SM quarks and leptons mediate its loop-induced couplings to photons and gluons. The electroweak gauge symmetry almost unambiguously dictates the bulk matter content and hence the hierarchies of the $S\to γγ$, $WW$, $ZZ$, $Zγ$, $t\bar t$ and dijet decay rates. We find that the $S\to Zγ$ decay mode is strongly suppressed, such that $\mbox{Br}(S\to Zγ)/\mbox{Br}(S\toγγ)<0.1$. The hierarchy problem for the new scalar boson is solved in analogy with the Higgs boson by localizing it near the infrared brane. The infinite sums over the Kaluza-Klein towers of fermion states are finite and can be calculated in closed form with a remarkably simple result. Reproducing the observed $pp\to S\toγγ$ signal requires Kaluza-Klein masses in the multi-TeV range, consistent with bounds from flavor physics and electroweak precision observables. Useful side products of our analysis, which can be adapted to almost any model for the diphoton resonance, are the calculation of the gluon-fusion production cross section $σ(pp\to S)$ at NNLO in QCD, an exact expression for the inclusive $S\to gg$ decay rate at N$^3$LO, a study of the $S\to t\bar t h$ three-body decay and a phenomenological analysis of portal couplings connecting $S$ with the Higgs field.

preprint2016arXiv

Extreme Scale-out SuperMUC Phase 2 - lessons learned

In spring 2015, the Leibniz Supercomputing Centre (Leibniz-Rechenzentrum, LRZ), installed their new Peta-Scale System SuperMUC Phase2. Selected users were invited for a 28 day extreme scale-out block operation during which they were allowed to use the full system for their applications. The following projects participated in the extreme scale-out workshop: BQCD (Quantum Physics), SeisSol (Geophysics, Seismics), GPI-2/GASPI (Toolkit for HPC), Seven-League Hydro (Astrophysics), ILBDC (Lattice Boltzmann CFD), Iphigenie (Molecular Dynamic), FLASH (Astrophysics), GADGET (Cosmological Dynamics), PSC (Plasma Physics), waLBerla (Lattice Boltzmann CFD), Musubi (Lattice Boltzmann CFD), Vertex3D (Stellar Astrophysics), CIAO (Combustion CFD), and LS1-Mardyn (Material Science). The projects were allowed to use the machine exclusively during the 28 day period, which corresponds to a total of 63.4 million core-hours, of which 43.8 million core-hours were used by the applications, resulting in a utilization of 69%. The top 3 users were using 15.2, 6.4, and 4.7 million core-hours, respectively.

preprint2016arXiv

Flavor Anomalies, the Diphoton Excess and a Dark Matter Candidate

We argue that the diphoton excess recently reported by ATLAS and CMS can be explained, along with several anomalies seen in the flavor sector, in models where a Standard-Model singlet scalar resonance with mass $M \approx 750$ GeV is produced in gluon fusion via loops containing a scalar color- triplet leptoquark $ϕ$. For a leptoquark mass $M_ϕ\lesssim$ 1 TeV, the production cross section is naturally in the 10 fb range. A large $S \to γγ$ branching ratio can be obtained by coupling the scalar S to new color-singlet fermions $χ$ with electroweak scale masses, which can be part of an $SU(2)_L$ multiplet, whose neutral component has the right mass and quantum numbers to be a dark matter candidate. Our model reveals a connection between flavor anomalies, the nature of dark matter and a new scalar, which acts as a mediator to the dark sector. The loop-mediated decay $S \to τ^+τ^-$ could be a striking signature of this model.

preprint2016arXiv

Hunting the Flavon

The next generation of experiments in particle physics will for the first time systematically test flavor physics models based on flavon fields. Starting from the current quark-flavor constrains on such models we show how the new generation of lepton flavor experiments will dominate indirect searches in the coming decades. A future 100 TeV hadron collider will then be the first experiment to probe flavons as propagating degrees of freedom. Our estimate of the collider reach relies on a proper treatment of backgrounds and detector effects. Complementary searches for indirect effects in lepton flavor experiments and propagating degrees of freedom at colliders are very limited at the LHC, but will be a new feature at a 100 TeV hadron collider.

preprint2016arXiv

One Leptoquark to Rule Them All: A Minimal Explanation for $R_{D^{(*)}}$, $R_K$ and $(g-2)_μ$

We show that by adding a single new scalar particle to the Standard Model, a TeV-scale leptoquark with the quantum numbers of a right-handed down quark, one can explain in a natural way three of the most striking anomalies of particle physics: the violation of lepton universality in $\bar B\to\bar K\ell^+\ell^-$ decays, the enhanced $\bar B\to D^{(*)}τ\barν$ decay rates, and the anomalous magnetic moment of the muon. Constraints from other precision measurements in the flavor sector can be satisfied without fine-tuning. Our model predicts enhanced $\bar B\to\bar K^{(*)}ν\barν$ decay rates and a new-physics contribution to $B_s\!-\!\bar B_s$ mixing close to the current central fit value.

preprint2016arXiv

The "forgotten" decay S -> Z+h as a CP analyzer

Scalar particles $S$ which are gauge singlets under the Standard Model are generic features of many models of fundamental physics, in particular as possible mediators to a hidden or dark sector. We show that the decay $S\to Zh$ provides a powerful probe of the CP nature of the scalar, because it is allowed only if $S$ has CP-odd interactions. We perform a model-independent analysis of this decay in the context of an effective Lagrangian and compute the relevant Wilson coefficients arising from integrating out heavy fermions to one-loop order. We illustrate our findings with the example of the 750 GeV diphoton resonance seen by ATLAS and CMS and show that the $S\to Zh$ decay rate could naturally be of similar magnitude or larger than the diphoton rate.

preprint2015arXiv

A Python Extension for the Massively Parallel Multiphysics Simulation Framework waLBerla

We present a Python extension to the massively parallel HPC simulation toolkit waLBerla. waLBerla is a framework for stencil based algorithms operating on block-structured grids, with the main application field being fluid simulations in complex geometries using the lattice Boltzmann method. Careful performance engineering results in excellent node performance and good scalability to over 400,000 cores. To increase the usability and flexibility of the framework, a Python interface was developed. Python extensions are used at all stages of the simulation pipeline: They simplify and automate scenario setup, evaluation, and plotting. We show how our Python interface outperforms the existing text-file-based configuration mechanism, providing features like automatic nondimensionalization of physical quantities and handling of complex parameter dependencies. Furthermore, Python is used to process and evaluate results while the simulation is running, leading to smaller output files and the possibility to adjust parameters dependent on the current simulation state. C++ data structures are exported such that a seamless interfacing to other numerical Python libraries is possible. The expressive power of Python and the performance of C++ make development of efficient code with low time effort possible.

preprint2015arXiv

Creating the Fermion Mass Hierarchies with Multiple Higgs Bosons

After the Higgs boson discovery, it is established that the Higgs mechanism explains electroweak symmetry breaking and generates the masses of all particles in the Standard Model, with the possible exception of neutrino masses. The hierarchies among fermion masses and mixing angles remain however unexplained. We propose a new class of two Higgs doublet models in which a flavor symmetry broken at the electroweak scale addresses this problem. The models are strongly constrained by electroweak precision tests and the fact that they produce modifications to Higgs couplings and flavor changing neutral currents; they are also constrained by collider searches for extra scalar bosons. The surviving models are very predictive, implying unavoidable new physics signals at the CERN Large Hadron Collider, e.g. extra Higgs Bosons with masses $M < 700$ GeV.

preprint2015arXiv

Curve Matching with Applications in Medical Imaging

In the recent years, Riemannian shape analysis of curves and surfaces has found several applications in medical image analysis. In this paper we present a numerical discretization of second order Sobolev metrics on the space of regular curves in Euclidean space. This class of metrics has several desirable mathematical properties. We propose numerical solutions for the initial and boundary value problems of finding geodesics. These two methods are combined in a Riemannian gradient-based optimization scheme to compute the Karcher mean. We apply this to a study of the shape variation in HeLa cell nuclei and cycles of cardiac deformations, by computing means and principal modes of variations.

preprint2015arXiv

Diffeomorphic density matching by optimal information transport

We address the following problem: given two smooth densities on a manifold, find an optimal diffeomorphism that transforms one density into the other. Our framework builds on connections between the Fisher-Rao information metric on the space of probability densities and right-invariant metrics on the infinite-dimensional manifold of diffeomorphisms. This optimal information transport, and modifications thereof, allows us to construct numerical algorithms for density matching. The algorithms are inherently more efficient than those based on optimal mass transport or diffeomorphic registration. Our methods have applications in medical image registration, texture mapping, image morphing, non-uniform random sampling, and mesh adaptivity. Some of these applications are illustrated in examples.

preprint2015arXiv

Flavor from the Electroweak Scale

We discuss the possibility that flavor hierarchies arise from the electroweak scale in a two Higgs doublet model, in which the two Higgs doublets jointly act as the flavon. Quark masses and mixing angles are explained by effective Yukawa couplings, generated by higher dimensional operators involving quarks and Higgs doublets. Modified Higgs couplings yield important effects on the production cross sections and decay rates of the light Standard Model like Higgs. In addition, flavor changing neutral currents arise at tree-level and lead to strong constraints from meson-antimeson mixing. Remarkably, flavor constraints turn out to prefer a region in parameter space that is in excellent agreement with the one preferred by recent Higgs precision measurements at the Large Hadron Collider (LHC). Direct searches for extra scalars at the LHC lead to further constraints. Precise predictions for the production and decay modes of the additional Higgs bosons are derived, and we present benchmark scenarios for searches at the LHC Run II. Flavor breaking at the electroweak scale as well as strong coupling effects demand a UV completion at the scale of a few TeV, possibly within the reach of the LHC.

preprint2015arXiv

Hörmander's condition for normal bundles on spaces of immersions

Several representations of geometric shapes involve quotients of mapping spaces. The projection onto the quotient space defines two sub-bundles of the tangent bundle, called the horizontal and vertical bundle. We investigate in these notes the sub-Riemannian geometries of these bundles. In particular, we show for a selection of bundles which naturally occur in applications that they are either bracket generating or integrable.

preprint2015arXiv

Landmark-Guided Elastic Shape Analysis of Human Character Motions

Motions of virtual characters in movies or video games are typically generated by recording actors using motion capturing methods. Animations generated this way often need postprocessing, such as improving the periodicity of cyclic animations or generating entirely new motions by interpolation of existing ones. Furthermore, search and classification of recorded motions becomes more and more important as the amount of recorded motion data grows. In this paper, we will apply methods from shape analysis to the processing of animations. More precisely, we will use the by now classical elastic metric model used in shape matching, and extend it by incorporating additional inexact feature point information, which leads to an improved temporal alignment of different animations.

preprint2015arXiv

Massively Parallel Phase-Field Simulations for Ternary Eutectic Directional Solidification

Microstructures forming during ternary eutectic directional solidification processes have significant influence on the macroscopic mechanical properties of metal alloys. For a realistic simulation, we use the well established thermodynamically consistent phase-field method and improve it with a new grand potential formulation to couple the concentration evolution. This extension is very compute intensive due to a temperature dependent diffusive concentration. We significantly extend previous simulations that have used simpler phase-field models or were performed on smaller domain sizes. The new method has been implemented within the massively parallel HPC framework waLBerla that is designed to exploit current supercomputers efficiently. We apply various optimization techniques, including buffering techniques, explicit SIMD kernel vectorization, and communication hiding. Simulations utilizing up to 262,144 cores have been run on three different supercomputing architectures and weak scalability results are shown. Additionally, a hierarchical, mesh-based data reduction strategy is developed to keep the I/O problem manageable at scale.

preprint2015arXiv

Metrics with prescribed horizontal bundle on spaces of curve

We study metrics on the shape space of curves that induce a prescribed splitting of the tangent bundle. More specifically, we consider reparametrization invariant metrics $G$ on the space $\operatorname{Imm}(S^1,\mathbb R^2)$ of parametrized regular curves. For many metrics the tangent space $T_c\operatorname{Imm}(S^1,\mathbb R^2)$ at each curve $c$ splits into vertical and horizontal components (with respect to the projection onto the shape space $B_i(S^1,\mathbb R^2)=\operatorname{Imm}(S^1,\mathbb R^2)/\operatorname{Diff}(S^1)$ of unparametrized curves and with respect to the metric $G$). In a previous article we characterized all metrics $G$ such that the induced splitting coincides with the natural splitting into normal and tangential parts. In these notes we extend this analysis to characterize all metrics that induce any prescribed splitting of the tangent bundle.

preprint2015arXiv

Second order elastic metrics on the shape space of curves

Second order Sobolev metrics on the space of regular unparametrized planar curves have several desirable completeness properties not present in lower order metrics, but numerics are still largely missing. In this paper, we present algorithms to numerically solve the initial and boundary value problems for geodesics. The combination of these algorithms allows to compute Karcher means in a Riemannian gradient-based optimization scheme. Our framework has the advantage that the constants determining the weights of the zero, first, and second order terms of the metric can be chosen freely. Moreover, due to its generality, it could be applied to more general spaces of mapping. We demonstrate the effectiveness of our approach by analyzing a collection of shapes representing physical objects.

preprint2014arXiv

A New Riemannian Setting for Surface Registration

We present a new approach for matching regular surfaces in a Riemannian setting. We use a Sobolev type metric on deformation vector fields which form the tangent bundle to the space of surfaces. In this article we compare our approach with the diffeomorphic matching framework. In the latter approach a deformation is prescribed on the ambient space, which then drags along an embedded surface. In contrast our metric is defined directly on the deformation vector field and can therefore be called an inner metric. We also show how to discretize the corresponding geodesic equation and compute the gradient of the cost functional using finite elements.

preprint2014arXiv

Constructing reparametrization invariant metrics on spaces of plane curves

Metrics on shape space are used to describe deformations that take one shape to another, and to determine a distance between them. We study a family of metrics on the space of curves, that includes several recently proposed metrics, for which the metrics are characterised by mappings into vector spaces where geodesics can be easily computed. This family consists of Sobolev-type Riemannian metrics of order one on the space $\text{Imm}(S^1,\mathbb R^2)$ of parametrized plane curves and the quotient space $\text{Imm}(S^1,\mathbb R^2)/\text{Diff}(S^1)$ of unparametrized curves. For the space of open parametrized curves we find an explicit formula for the geodesic distance and show that the sectional curvatures vanish on the space of parametrized and are non-negative on the space of unparametrized open curves. For the metric, which is induced by the "R-transform", we provide a numerical algorithm that computes geodesics between unparameterised, closed curves, making use of a constrained formulation that is implemented numerically using the RATTLE algorithm. We illustrate the algorithm with some numerical tests that demonstrate it's efficiency and robustness.

preprint2014arXiv

Light Dark Matter, Naturalness, and the Radiative Origin of the Electroweak Scale

We study classically scale invariant models in which the Standard Model Higgs mass term is replaced in the Lagrangian by a Higgs portal coupling to a complex scalar field of a dark sector. We focus on models that are weakly coupled with the quartic scalar couplings nearly vanishing at the Planck scale. The dark sector contains fermions and scalars charged under dark SU(2) x U(1) gauge interactions. Radiative breaking of the dark gauge group triggers electroweak symmetry breaking through the Higgs portal coupling. Requiring both a Higgs boson mass of 125.5 GeV and stability of the Higgs potential up to the Planck scale implies that the radiative breaking of the dark gauge group occurs at the TeV scale. We present a particular model which features a long-range abelian dark force. The dominant dark matter component is neutral dark fermions, with the correct thermal relic abundance, and in reach of future direct detection experiments. The model also has lighter stable dark fermions charged under the dark force, with observable effects on galactic-scale structure. Collider signatures include a dark sector scalar boson with mass < 250 GeV that decays through mixing with the Higgs boson, and can be detected at the LHC. The Higgs boson, as well as the new scalar, may have significant invisible decays into dark sector particles.

preprint2014arXiv

Metrics on Spaces of Immersions where Horizontality Equals Normality

We study metrics on shape space of immersions that have a particularly simple horizontal bundle. More specifically, we consider reparametrization invariant Sobolev metrics $G$ on the space $\operatorname{Imm}(M,N)$ of immersions of a compact manifold $M$ in a Riemannian manifold $(N,\overline{g})$. The tangent space $T_f\operatorname{Imm}(M,N)$ at each immersion $f$ has two natural splittings: one into components that are tangential/normal to the surface $f$ (with respect to $\overline{g}$) and another one into vertical/horizontal components (with respect to the projection onto the shape space $B_i(M,N)=\operatorname{Imm}(M,N)/\operatorname{Diff}(M)$ of unparametrized immersions and with respect to the metric $G$). The first splitting can be easily calculated numerically, while the second splitting is important because it mirrors the geometry of shape space and geodesics thereon. Motivated by facilitating the numerical calculation of geodesics on shape space, we characterize all metrics $G$ such that the two splittings coincide. In the special case of planar curves, we show that the regularity of curves in the metric completion can be controlled by choosing a strong enough metric within this class. We demonstrate in several examples that our approach allows us to efficiently calculate numerical solutions of the boundary value problem for geodesics on shape space.

preprint2014arXiv

Optical Flow on Moving Manifolds

Optical flow is a powerful tool for the study and analysis of motion in a sequence of images. In this article we study a Horn-Schunck type spatio-temporal regularization functional for image sequences that have a non-Euclidean, time varying image domain. To that end we construct a Riemannian metric that describes the deformation and structure of this evolving surface. The resulting functional can be seen as natural geometric generalization of previous work by Weickert and Schnörr (2001) and Lefèvre and Baillet (2008) for static image domains. In this work we show the existence and wellposedness of the corresponding optical flow problem and derive necessary and sufficient optimality conditions. We demonstrate the functionality of our approach in a series of experiments using both synthetic and real data.

preprint2014arXiv

Soft Collinear Effective Theory for Heavy WIMP Annihilation

In a large class of models for Weakly Interacting Massive Particles (WIMPs), the WIMP mass $M$ lies far above the weak scale $m_W$. This work identifies universal Sudakov-type logarithms $\sim α\log^2 (2\,M/m_W)$ that spoil the naive convergence of perturbation theory for annihilation processes. An effective field theory (EFT) framework is presented, allowing the systematic resummation of these logarithms. Another impact of the large separation of scales is that a long-distance wave-function distortion from electroweak boson exchange leads to observable modifications of the cross section. Careful accounting of momentum regions in the EFT allows the rigorously disentanglement of this so-called Sommerfeld enhancement from the short distance hard annihilation process. The WIMP is modeled as a heavy-particle field, while the light, energetic, final-state electroweak gauge bosons are treated as soft and collinear fields. Hard matching coefficients are computed at renormalization scale $μ\sim 2\,M$, then evolved down to $μ\sim m_W$, where electroweak symmetry breaking is incorporated and the matching onto the relevant quantum mechanical Hamiltonian is performed. The example of an $SU(2)_W$ triplet scalar dark matter candidate annihilating to line photons is used for concreteness, allowing the numerical exploration of the impact of next-to-leading order corrections and log resummation. For $M \simeq 3$ TeV, the resummed Sommerfeld enhanced cross section is reduced by a factor of $\sim 3$ with respect to the tree-level fixed order result.

preprint2014arXiv

The homogeneous Sobolev metric of order one on diffeomorphism groups on the real line

In this article we study Sobolev metrics of order one on diffeomorphism groups on the real line. We prove that the space $\operatorname{Diff}_{1}(\mathbb R)$ equipped with the homogenous Sobolev metric of order one is a flat space in the sense of Riemannian geometry, as it is isometric to an open subset of a mapping space equipped with the flat $L^2$-metric. Here $\operatorname{Diff}_{1}(\mathbb R)$ denotes the extension of the group of all either compactly supported, rapidly decreasing or $H^\infty$-diffeomorphisms, that allows for a shift towards infinity. In particular this result provides an analytic solution formula for the corresponding geodesic equation, the non-periodic Hunter-Saxton equation. In addition we show that one can obtain a similar result for the two-component Hunter-Saxton equation and discuss the case of the non-homogenous Sobolev one metric which is related to the Camassa-Holm equation.

preprint2013arXiv

$R$-transforms for Sobolev $H^2$-metrics on spaces of plane curves

We consider spaces of smooth immersed plane curves (modulo translations and/or rotations), equipped with reparameterization invariant weak Riemannian metrics involving second derivatives. This includes the full $H^2$-metric without zero order terms. We find isometries (called $R$-transforms) from some of these spaces into function spaces with simpler weak Riemannian metrics, and we use this to give explicit formulas for geodesics, geodesic distances, and sectional curvatures. We also show how to utilise the isometries to compute geodesics numerically.

preprint2013arXiv

Almost local metrics on shape space of hypersurfaces in n-space

This paper extends parts of the results from [P.W.Michor and D. Mumford, \emph{Appl. Comput. Harmon. Anal.,} 23 (2007), pp. 74--113] for plane curves to the case of hypersurfaces in $\mathbb R^n$. Let $M$ be a compact connected oriented $n-1$ dimensional manifold without boundary like the sphere or the torus. Then shape space is either the manifold of submanifolds of $\mathbb R^n$ of type $M$, or the orbifold of immersions from $M$ to $\mathbb R^n$ modulo the group of diffeomorphisms of $M$. We investigate almost local Riemannian metrics on shape space. These are induced by metrics of the following form on the space of immersions: $$ G_f(h,k) = \int_{M} Φ(\on{Vol}(f),\operatorname{Tr}(L))\g(h, k) \operatorname{vol}(f^*\g),$$ where $\g$ is the Euclidean metric on $\mathbb R^n$, $f^*\g$ is the induced metric on $M$, $h,k\in C^\infty(M,\mathbb R^n)$ are tangent vectors at $f$ to the space of embeddings or immersions, where $Φ:\mathbb R^2\to \mathbb R_{>0}$ is a suitable smooth function, $\operatorname{Vol}(f) = \int_M\operatorname{vol}(f^*\g)$ is the total hypersurface volume of $f(M)$, and the trace $\operatorname{Tr}(L)$ of the Weingarten mapping is the mean curvature. For these metrics we compute the geodesic equations both on the space of immersions and on shape space, the conserved momenta arising from the obvious symmetries, and the sectional curvature. For special choices of $Φ$ we give complete formulas for the sectional curvature. Numerical experiments illustrate the behavior of these metrics.

preprint2013arXiv

Constraining RS Models by Future Flavor and Collider Measurements: A Snowmass Whitepaper

Randall-Sundrum models are models of quark flavor, because they explain the hierarchies in the quark masses and mixings in terms of order one localization parameters of extra dimensional wavefunctions. The same small numbers which generate the light quark masses suppress contributions to flavor violating tree level amplitudes. In this note we update universal constraints from electroweak precision parameters and demonstrate how future measurements of flavor violation in ultra rare decay channels of Kaons and B mesons will constrain the parameter space of this type of models. We show how collider signatures are correlated with these flavor measurements and compute projected limits for direct searches at the 14 TeV LHC run, a 14 TeV LHC luminosity upgrade, a 33 TeV LHC energy upgrade, and a potential 100 TeV machine. We further discuss the effects of a warped model of leptons in future measurements of lepton flavor violation.

preprint2013arXiv

Overview of the Geometries of Shape Spaces and Diffeomorphism Groups

This article provides an overview of various notions of shape spaces, including the space of parametrized and unparametrized curves, the space of immersions, the diffeomorphism group and the space of Riemannian metrics. We discuss the Riemannian metrics that can be defined thereon, and what is known about the properties of these metrics. We put particular emphasis on the induced geodesic distance, the geodesic equation and its well-posedness, geodesic and metric completeness and properties of the curvature.

preprint2012arXiv

A Solution to the Flavor Problem of Warped Extra-Dimension Models

A minimal solution to the flavor problem of warped extra-dimension models, i.e. the excessive mixed-chirality contribution to CP violation in K-Kbar mixing arising from Kaluza-Klein (KK) gluon exchange, is proposed. Extending the strong-interaction gauge group in the bulk by an additional SU(3), and breaking this symmetry to QCD via boundary conditions, the constraints arising from the epsilon_K parameter are significantly relaxed. As a result, KK scales M_KK ~ 2 TeV are consistent with all flavor observables without significant fine-tuning. The model predicts the existence of an extended Higgs sector featuring massive color-octet scalars and a tower of KK pseudo-axial gluon resonances, whose existence is not in conflict with recent LHC dijet bounds.

preprint2012arXiv

Geodesic distance for right invariant Sobolev metrics of fractional order on the diffeomorphism group

We study Sobolev-type metrics of fractional order $s\geq0$ on the group $\Diff_c(M)$ of compactly supported diffeomorphisms of a manifold $M$. We show that for the important special case $M=S^1$ the geodesic distance on $\Diff_c(S^1)$ vanishes if and only if $s\leq\frac12$. For other manifolds we obtain a partial characterization: the geodesic distance on $\Diff_c(M)$ vanishes for $M=\R\times N, s<\frac12$ and for $M=S^1\times N, s\leq\frac12$, with $N$ being a compact Riemannian manifold. On the other hand the geodesic distance on $\Diff_c(M)$ is positive for $\dim(M)=1, s>\frac12$ and $\dim(M)\geq2, s\geq1$. For $M=\R^n$ we discuss the geodesic equations for these metrics. For $n=1$ we obtain some well known PDEs of hydrodynamics: Burgers' equation for $s=0$, the modified Constantin-Lax-Majda equation for $s=\frac 12$ and the Camassa-Holm equation for $s=1$.

preprint2012arXiv

Local Uniqueness of the Circular Integral Invariant

This article is concerned with the representation of curves by means of integral invariants. In contrast to the classical differential invariants they have the advantage of being less sensitive with respect to noise. The integral invariant most common in use is the circular integral invariant. A major drawback of this curve descriptor, however, is the absence of any uniqueness result for this representation. This article serves as a contribution towards closing this gap by showing that the circular integral invariant is injective in a neighbourhood of the circle. In addition, we provide a stability estimate valid on this neighbourhood. The proof is an application of Riesz-Schauder theory and the implicit function theorem in a Banach space setting.

preprint2012arXiv

Sobolev Metrics on Shape Space, II: Weighted Sobolev Metrics and Almost Local Metrics

In continuation of [3] we discuss metrics of the form $$ G^P_f(h,k)=\int_M \sum_{i=0}^pΦ_i(\Vol(f)) \g((P_i)_fh,k) \vol(f^*\g) $$ on the space of immersions $\Imm(M,N)$ and on shape space $B_i(M,N)=\Imm(M,N)/\on{Diff}(M)$. Here $(N,\g)$ is a complete Riemannian manifold, $M$ is a compact manifold, $f:M\to N$ is an immersion, $h$ and $k$ are tangent vectors to $f$ in the space of immersions, $f^*\g$ is the induced Riemannian metric on $M$, $\vol(f^*\g)$ is the induced volume density on $M$, $\Vol(f)=\int_M\vol(f^*\g)$, $Φ_i$ are positive real-valued functions, and $(P_i)_f$ are operators like some power of the Laplacian $Δ^{f^*\g}$. We derive the geodesic equations for these metrics and show that they are sometimes well-posed with the geodesic exponential mapping a local diffeomorphism. The new aspect here are the weights $\Ph_i(\Vol(f))$ which we use to construct scale invariant metrics and order 0 metrics with positive geodesic distance. We treat several concrete special cases in detail.

preprint2012arXiv

Sobolev metrics on the manifold of all Riemannian metrics

On the manifold $\Met(M)$ of all Riemannian metrics on a compact manifold $M$ one can consider the natural $L^2$-metric as described first by \cite{Ebin70}. In this paper we consider variants of this metric which in general are of higher order. We derive the geodesic equations, we show that they are well-posed under some conditions and induce a locally diffeomorphic geodesic exponential mapping. We give a condition when Ricci flow is a gradient flow for one of these metrics.

preprint2011arXiv

Curvature weighted metrics on shape space of hypersurfaces in $n$-space

Let $M$ be a compact connected oriented $n-1$ dimensional manifold without boundary. In this work, shape space is the orbifold of unparametrized immersions from $M$ to $\mathbb R^n$. The results of \cite{Michor118}, where mean curvature weighted metrics were studied, suggest incorporating Gauß curvature weights in the definition of the metric. This leads us to study metrics on shape space that are induced by metrics on the space of immersions of the form $$ G_f(h,k) = \int_{M} Φ. \bar g(h, k) \vol(f^*\bar{g}).$$ Here $f \in \Imm(M,\R^n)$ is an immersion of $M$ into $\R^n$ and $h,k\in C^\infty(M,\mathbb R^n)$ are tangent vectors at $f$. $\bar g$ is the standard metric on $\mathbb R^n$, $f^*\bar g$ is the induced metric on $M$, $\vol(f^*\bar g)$ is the induced volume density and $Φ$ is a suitable smooth function depending on the mean curvature and Gauß curvature. For these metrics we compute the geodesic equations both on the space of immersions and on shape space and the conserved momenta arising from the obvious symmetries. Numerical experiments illustrate the behavior of these metrics.

preprint2011arXiv

Sobolev metrics on shape space of surfaces

Let $M$ and $N$ be connected manifolds without boundary with $\dim(M) < \dim(N)$, and let $M$ compact. Then shape space in this work is either the manifold of submanifolds of $N$ that are diffeomorphic to $M$, or the orbifold of unparametrized immersions of $M$ in $N$. We investigate the Sobolev Riemannian metrics on shape space: These are induced by metrics of the following form on the space of immersions: $$ G^P_f(h,k) = \int_{M} \g(P^f h, k)\, \vol(f^*\g)$$ where $\g$ is some fixed metric on $N$, $f^*\g$ is the induced metric on $M$, $h,k \in Γ(f^*TN)$ are tangent vectors at $f$ to the space of embeddings or immersions, and $P^f$ is a positive, selfadjoint, bijective scalar pseudo differential operator of order $2p$ depending smoothly on $f$. We consider later specifically the operator $P^f=1 + AΔ^p$, where $Δ$ is the Bochner-Laplacian on $M$ induced by the metric $f^*\bar g$. For these metrics we compute the geodesic equations both on the space of immersions and on shape space, and also the conserved momenta arising from the obvious symmetries. We also show that the geodesic equation is well-posed on spaces of immersions, and also on diffeomorphism groups. We give examples of numerical solutions.

preprint2010arXiv

Top-Quark Forward-Backward Asymmetry in Randall-Sundrum Models Beyond the Leading Order

We calculate the ttbar forward-backward asymmetry, A_FB^t, in Randall-Sundrum (RS) models taking into account the dominant next-to-leading order (NLO) corrections in QCD. At Born level we include the exchange of Kaluza-Klein (KK) gluons and photons, the Z boson and its KK excitations, as well as the Higgs boson, whereas beyond the leading order (LO) we consider the interference of tree-level KK-gluon exchange with one-loop QCD box diagrams and the corresponding bremsstrahlungs corrections. We find that the strong suppression of LO effects, that arises due to the elementary nature and the mostly vector-like couplings of light quarks, is lifted at NLO after paying the price of an additional factor of alpha_s/(4 pi). In spite of this enhancement, the resulting RS corrections in A_FB^t remain marginal, leaving the predicted asymmetry SM-like. As our arguments are solely based on the smallness of the axial-vector couplings of light quarks to the strong sector, our findings are model-independent and apply to many scenarios of new physics that address the flavor problem via geometrical sequestering.

preprint2009arXiv

Randall-Sundrum Models and Precision Observables

I present a review of phenomenological implications of the Randall-Sundrum (RS) model with bulk fermions and brane-localized Higgs boson. Modifications to the W-boson mass, corrections to the Peskin-Takeuchi parameters and to the Z -> bb couplings will be discussed. From these observables severe bounds on the mass scale of Kaluza-Klein (KK) modes arise. Constraints from all three observables are very sensitive to the exact value of the Higgs boson mass and the bounds can be significantly lowered by allowing for a heavy Higgs boson (m_h ~ 1 TeV). Consequences thereof, as well as other approaches like "little RS" models and models with custodial symmetry will also be briefly discussed.