CPsuperH2.3: an Updated Tool for Phenomenology in the MSSM with Explicit CP Violation

We describe the Fortran code CPsuperH2.3, which incorporates the following updates compared with its predecessor CPsuperH2.0. It implements improved calculations of the Higgs-boson masses and mixing including stau contributions and finite threshold effects on the tau-lepton Yukawa coupling. It incorporates the LEP limits on the processes e^+ e^- to H_i Z, H_i H_j and the CMS limits on H_i to tau^+ tau^- obtained from 4.6/fb of data at a centre-of-mass energy of 7 TeV. It also includes the decay mode H_i to Z gamma and the Schiff-moment contributions to the electric dipole moments of Mercury and Radium225, with several calculational options for the case of Mercury. These additions make CPsuperH2.3 a suitable tool for analyzing possible CP-violating effects in the MSSM in the era of the LHC and a new generation of EDM experimentsComment: 31 pages, 10 eps figures, 7 tables; H to Z gamma and SM BRs included; To appear in CPC; Typos in Eq.(A.2) corrected;The program may be obtained from http://www.hep.man.ac.uk/u/jslee/CPsuperH.html, or by contacting the first author at [email protected]; A comment added after Eq.(15) and a typo in Eq.(A.4) correcte

2D Full-Wave Simulation of Waves in Space and Tokamak Plasmas

Simulation results using a 2D full-wave code (FW2D) for space and NSTX fusion plasmas are presented. The FW2D code solves the cold plasma wave equations using the finite element method. The wave code has been successfully applied to describe low frequency waves in planetary magnetospheres (i.e., dipole geometry) and the results include generation and propagation of externally driven ultra-low frequency waves via mode conversion at Mercury and mode coupling, refraction and reflection of internally driven field-aligned propagating left-handed electromagnetic ion cyclotron (EMIC) waves at Earth. In this paper, global structure of linearly polarized EMIC waves is examined and the result shows such resonant wave modes can be localized near the equatorial plane. We also adopt the FW2D code to tokamak geometry and examine radio frequency (RF) waves in the scape-off layer (SOL) of tokamaks. By adopting the rectangular and limiter boundary, we compare the results with existing AORSA simulations. The FW2D code results for the high harmonic fast wave heating case on NSTX with a rectangular vessel boundary shows excellent agreement with the AORSA code

Region-based memory management for Mercury programs

Region-based memory management (RBMM) is a form of compile time memory management, well-known from the functional programming world. In this paper we describe our work on implementing RBMM for the logic programming language Mercury. One interesting point about Mercury is that it is designed with strong type, mode, and determinism systems. These systems not only provide Mercury programmers with several direct software engineering benefits, such as self-documenting code and clear program logic, but also give language implementors a large amount of information that is useful for program analyses. In this work, we make use of this information to develop program analyses that determine the distribution of data into regions and transform Mercury programs by inserting into them the necessary region operations. We prove the correctness of our program analyses and transformation. To execute the annotated programs, we have implemented runtime support that tackles the two main challenges posed by backtracking. First, backtracking can require regions removed during forward execution to be "resurrected"; and second, any memory allocated during a computation that has been backtracked over must be recovered promptly and without waiting for the regions involved to come to the end of their life. We describe in detail our solution of both these problems. We study in detail how our RBMM system performs on a selection of benchmark programs, including some well-known difficult cases for RBMM. Even with these difficult cases, our RBMM-enabled Mercury system obtains clearly faster runtimes for 15 out of 18 benchmarks compared to the base Mercury system with its Boehm runtime garbage collector, with an average runtime speedup of 24%, and an average reduction in memory requirements of 95%. In fact, our system achieves optimal memory consumption in some programs.Comment: 74 pages, 23 figures, 11 tables. A shorter version of this paper, without proofs, is to appear in the journal Theory and Practice of Logic Programming (TPLP

A Comprehensive Analysis of Electric Dipole Moment Constraints on CP-violating Phases in the MSSM

We analyze the constraints placed on individual, flavor diagonal CP-violating phases in the minimal supersymmetric extension of the Standard Model (MSSM) by current experimental bounds on the electric dipole moments (EDMs) of the neutron, Thallium, and Mercury atoms. We identify the four CP-violating phases that are individually highly constrained by current EDM bounds, and we explore how these phases and correlations among them are constrained by current EDM limits. We also analyze the prospective implications of the next generation of EDM experiments. We point out that all other CP-violating phases in the MSSM are not nearly as tightly constrained by limits on the size of EDMs. We emphasize that a rich set of phenomenological consequences is potentially associated with these generically large EDM-allowed phases, ranging from B physics, electroweak baryogenesis, and signals of CP-violation at the CERN Large Hadron Collider and at future linear colliders. Our numerical study takes into account the complete set of contributions from one- and two-loop EDMs of the electron and quarks, one- and two-loop Chromo-EDMs of quarks, the Weinberg 3-gluon operator, and dominant 4-fermion CP-odd operator contributions, including contributions which are both included and not included yet in the CPsuperH2.0 package. We also introduce an open-source numerical package, 2LEDM, which provides the complete set of two-loop electroweak diagrams contributing to the electric dipole moments of leptons and quarks.Comment: 23 pages, 11 figures; v2: references added, minor change

First-order integer programming for MAP problems

Finding the most probable (MAP) model in SRL frameworks such as Markov logic and Problog can, in principle, be solved by encoding the problem as a `grounded-out' mixed integer program (MIP). However, useful first-order structure disappears in this process motivating the development of first-order MIP approaches. Here we present mfoilp, one such approach. Since the syntax and semantics of mfoilp is essentially the same as existing approaches we focus here mainly on implementation and algorithmic issues. We start with the (conceptually) simple problem of using a logic program to generate a MIP instance before considering more ambitious exploitation of first-order representations.Comment: corrected typo

Supersymmetric Electroweak Baryogenesis Via Resonant Sfermion Sources

We calculate the baryon asymmetry produced at the electroweak phase transition by quasi-degenerate third generation sfermions in the minimal supersymmetric extension of the Standard Model. We evaluate constraints from Higgs searches, from collider searches for supersymmetric particles, and from null searches for the permanent electric dipole moment (EDM) of the electron, of the neutron and of atoms. We find that resonant sfermion sources can in principle provide a large enough baryon asymmetry in various corners of the sfermion parameter space, and we focus, in particular, on the case of large $\tan\beta$, where third-generation down-type (s)fermions become relevant. We show that in the case of stop and sbottom sources, the viable parameter space is ruled out by constraints from the non-observation of the Mercury EDM. We introduce a new class of CP violating sources, quasi-degenerate staus, that escapes current EDM constraints while providing large enough net chiral currents to achieve successful "slepton-mediated" electroweak baryogenesis.Comment: 35 pages, 9 figures; v2: several revisions, but conclusions unchanged. Matches version published in PR