SuSeFLAV 1.2: Program for supersymmetric mass spectra with seesaw mechanism and rare lepton flavor violating decays

Accurate supersymmetric spectra are required to confront data from direct and indirect searches of supersymmetry. SuSeFLAV is a numerical tool capable of computing supersymmetric spectra precisely for various supersymmetric breaking scenarios applicable even in the presence of flavor violation. The program solves MSSM RGEs with complete 3 x 3 flavor mixing at 2-loop level and one loop finite threshold corrections to all MSSM parameters by incorporating radiative electroweak symmetry breaking conditions. The program also incorporates the Type-I seesaw mechanism with three massive right handed neutrinos at user defined mass scales and mixing. It also computes branching ratios of flavor violating processes such as l(j) -> l(i)gamma, l(j) -> 3 l(i), b -> s gamma and supersymmetric contributions to flavor conserving quantities such as (g(mu) - 2). A large choice of executables suitable for various operations of the program are provided. Program summary Program title: SuSeFLAV Catalogue identifier: AEOD_v1_0 Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AEOD_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: GNU General Public License No. of lines in distributed program, including test data, etc.: 76552 No. of bytes in distributed program, including test data, etc.: 582787 Distribution format: tar.gz Programming language: Fortran 95. Computer: Personal Computer, Work-Station. Operating system: Linux, Unix. Classification: 11.6. Nature of problem: Determination of masses and mixing of supersymmetric particles within the context of MSSM with conserved R-parity with and without the presence of Type-I seesaw. Inter-generational mixing is considered while calculating the mass spectrum. Supersymmetry breaking parameters are taken as inputs at a high scale specified by the mechanism of supersymmetry breaking. RG equations including full inter-generational mixing are then used to evolve these parameters up to the electroweak breaking scale. The low energy supersymmetric spectrum is calculated at the scale where successful radiative electroweak symmetry breaking occurs. At weak scale standard model fermion masses, gauge couplings are determined including the supersymmetric radiative corrections. Once the spectrum is computed, the program proceeds to various lepton flavor violating observables (e.g., BR(mu -> e gamma), BR(tau -> mu gamma) etc.) at the weak scale. Solution method: Two loop RGEs with full 3 x 3 flavor mixing for all supersymmetry breaking parameters are used to compute the low energy supersymmetric mass spectrum. An adaptive step size Runge-Kutta method is used to solve the RGEs numerically between the high scale and the electroweak breaking scale. Iterative procedure is employed to get the consistent radiative electroweak symmetry breaking condition. The masses of the supersymmetric particles are computed at 1-loop order. The third generation SM particles and the gauge couplings are evaluated at the 1-loop order including supersymmetric corrections. A further iteration of the full program is employed such that the SM masses and couplings are consistent with the supersymmetric particle spectrum. Additional comments: Several executables are presented for the user. Running time: 0.2 s on a Intel(R) Core(TM) i5 CPU 650 with 3.20 GHz. (c) 2012 Elsevier B.V. All rights reserved.

Some implications of lepton flavor violating processes in a supersymmetric Type II seesaw model at TeV scale

We have conceived a supersymmetric Type II seesaw model at TeV scale, which has some additional particles consisting of scalar and fermionic triplet Higgs states, whose masses are around a few hundred GeV. In this particular model, we have studied constraints on the masses of triplet states arising from the lepton flavor violating (LFV) processes, such as mu -> 3e and mu -> e gamma. We have analyzed the implications of these constraints on other observable quantities such as the muon anomalous magnetic moment and the decay patterns of scalar triplet Higgses. Scalar triplet Higgs states can decay into leptons and into supersymmetric fields. We have found that the constraints from LFV can affect these various decay modes.

Revisiting lepton flavor violation in a supersymmetric type II seesaw mechanism

In view of the recent measurement of the reactor mixing angle theta(13) and updated limit on BRd(mu -> e gamma) by the MEG experiment, we reexamine the charged lepton flavor violations in a framework of the supersymmetric type II seesaw mechanism. The supersymmetric type II seesaw predicts a strong correlation between BR(mu -> e gamma) and BR(tau -> mu gamma) mainly in terms of the neutrino mixing angles. We show that such a correlation can be determined accurately after the measurement of theta(13). We compute different factors that can affect this correlation and show that the minimal supergravity-like scenarios, in which slepton masses are taken to be universal at the high scale, predict 3.5 <= BR(tau -> mu gamma)/= BR(mu -> e gamma) <= 30 for normal hierarchical neutrino masses. Any experimental indication of deviation from this prediction would rule out the minimal models of the supersymmetric type II seesaw. We show that the current MEG limit puts severe constraints on the light sparticle spectrum in the minimal supergravity model if the seesaw scale lies within 10(13)-10(15) GeV. It is shown that these constraints can be relaxed and a relatively light sparticle spectrum can be obtained in a class of models in which the soft mass of a triplet scalar is taken to be nonuniversal at the high scale.

Implications of Higgs boson to diphoton decay rate in the bilinear R-parity violating supersymmetric model

The Large Hadron Collider has recently discovered a Higgs-like particle having a mass around 125 GeVand also indicated that there is an enhancement in the Higgs to diphoton decay rate as compared to that in the standard model. We have studied implications of these discoveries in the bilinear R-parity violating supersymmetric model, whose main motivation is to explain the nonzero masses for neutrinos. The R-parity violating parameters in this model are epsilon and b(epsilon), and these parameters determine the scale of neutrino masses. If the enhancement in the Higgs to diphoton decay rate is true, then we have found epsilon greater than or similar to 0.01 GeV and b epsilon similar to 1 GeV2 in order to be compatible with the neutrino oscillation data. Also, in the above mentioned analysis, we can determine the soft masses of sleptons (m(L)) and CP-odd Higgs boson mass (mA). We have estimated that m(L) greater than or similar to 300 GeV and m(A) greater than or similar to 700 GeV. We have also commented on the allowed values of epsilon and b(epsilon), in case there is no enhancement in the Higgs to diphoton decay rate. Finally, we present a model to explain the smallness of epsilon and b(epsilon).

Invisible decays of the lightest Higgs boson in supersymmetric models

We consider supersymmetric models in which the lightest Higgs scalar can decay invisibly consistent with the constraints on the 126 GeV state discovered at the CERN LHC. We consider the invisible decay in the minimal supersymmetric standard model (MSSM), as well its extension containing an additional chiral singlet superfield, the so-called next-to-minimal or nonminimal supersymmetric standard model (NMSSM). We consider the case of MSSM with both universal as well as nonuniversal gaugino masses at the grand unified scale, and find that only an E-6 grand unified model with unnaturally large representation can give rise to sufficiently light neutralinos which can possibly lead to the invisible decay h(0) -> (chi) over tilde (0)(1)(chi) over tilde (0)(1). Following this, we consider the case of NMSSM in detail, where we also find that it is not possible to have the invisible decay of the lightest Higgs scalar with universal gaugino masses at the grand unified scale. We delineate the regions of the NMSSM parameter space where it is possible for the lightest Higgs boson to have a mass of about 126 GeV, and then concentrate on the region where this Higgs can decay into light neutralinos, with the soft gaugino masses M-1 and M-2 as two independent parameters, unconstrained by grand unification. We also consider, simultaneously, the other important invisible Higgs decay channel in the NMSSM, namely the decay into the lightest CP-odd scalars, h(1) -> a(1)a(1), which is studied in detail. With the invisible Higgs branching ratio being constrained by the present LHC results, we find that mu(eff) < 170 GeV and M-1 < 80 GeV are disfavored in NMSSM for fixed values of the other input parameters. The dependence of our results on the parameters of NMSSM is discussed in detail.

Beam polarization effects in the radiative production of lightest neutralinos in e(+)e(-) collisions in supersymmetric grand unified models

We study the production of the lightest neutralinos in the process e(+)e(-) -> chi(0)(1)chi(0)(1)gamma in supersymmetric grand unified models for the International Linear Collider energies with longitudinally polarized beams. We consider cases where the standard model gauge group is unified into the grand unified gauge groups SU(5), or SO(10). We have carried out a comprehensive study of this process in the SU(5) and SO(10) grand unified theories which includes the QED radiative corrections. We compare and contrast the dependence of the signal cross section on the grand unified gauge group, and on the different representations of the grand unified gauge group, when the electron and positron beams are longitudinally polarized. To assess the feasibility of experimentally observing the radiative production process, we have also considered in detail the background to this process coming from the radiative neutrino production process e(+)e(-)-> nu(nu) over bar gamma with longitudinally polarized electron and positron beams. In addition we have also considered the supersymmetric background coming from the radiative production of scalar neutrinos in the process e(+)e(-) -> (nu) over tilde(nu) over tilde*gamma with longitudinally polarized beams. The process can be a major background to the radiative production of neutralinos when the scalar neutrinos decay invisibly.

Anatomy of Higgs mass in supersymmetric inverse seesaw models

We compute the one loop corrections to the CP-even Higgs mass matrix in the supersymmetric inverse seesaw model to single out the different cases where the radiative corrections from the neutrino sector could become important. It is found that there could be a significant enhancement in the Higgs mass even for Dirac neutrino masses of O(30) GeV if the left-handed sneutrino soft mass is comparable or larger than the right-handed neutrino mass. In the case where right-handed neutrino masses are significantly larger than the supersymmetry breaking scale, the corrections can utmost account to an upward shift of 3 GeV. For very heavy multi TeV sneutrinos, the corrections replicate the stop corrections at 1-loop. We further show that general gauge mediation with inverse seesaw model naturally accommodates a 125 GeV Higgs with TeV scale stops. (C) 2014 The Authors. Published by Elsevier B.V.

Role of supersymmetric heavy Higgs boson production in the self-coupling measurement of 125 GeV Higgs boson at the LHC

Measurement of the self-coupling of the 125 GeV Higgs boson is one of the most crucial tasks for a high luminosity run of the LHC, and it can only be measured in the di-Higgs final state. In the minimal supersymmetric standard model, heavy CP even Higgs (H) can decay into a lighter 125 GeV Higgs boson (h) and, therefore, can influence the rate of di-Higgs production. We investigate the role of single H production in the context of measuring the self-coupling of h. We have found that the H -> hh decay can change the value of Higgs (h) self-coupling substantially, in a low tan beta regime where the mass of the heavy Higgs boson lies between 250 and 600 GeV and, depending on the parameter space, it may be seen as an enhancement of the self-coupling of the 125 GeV Higgs boson.

Invisible decays of the heavier Higgs boson in the minimal supersymmetric standard model

We consider the possibility that the heavier CP-even Higgs boson (H-0) in the minimal supersymmetric standard model (MSSM) decays invisibly into neutralinos in the light of the recent discovery of the 126 GeV resonance at the CERN Large Hadron Collider (LHC). For this purpose we consider the minimal supersymmetric standard model with universal, nonuniversal and arbitrary boundary conditions on the supersymmetry breaking gaugino mass parameters at the grand unified scale. Typically, scenarios with universal and nonuniversal gaugino masses do not allow invisible decays of the lightest Higgs boson (h(0)), which is identified with the 126 GeV resonance, into the lightest neutralinos in the MSSM. With arbitrary gaugino masses at the grand unified scale, such an invisible decay is possible. The second lightest Higgs boson can decay into various invisible final states for a considerable region of the MSSM parameter space with arbitrary gaugino masses as well as with the gaugino masses restricted by universal and nonuniversal boundary conditions at the grand unified scale. The possibility of the second lightest Higgs boson of the MSSM decaying into invisible channels is more likely for arbitrary gaugino masses at the grand unified scale. The heavier Higgs boson decay into lighter particles leads to the intriguing possibility that the entire Higgs boson spectrum of the MSSM may be visible at the LHC even if it decays invisibly, during the searches for an extended Higgs boson sector at the LHC. In such a scenario the nonobservation of the extended Higgs sector of the MSSM may carefully be used to rule out regions of the MSSM parameter space at the LHC.

A method for coupling the phase-field model based on a grand-potential formalism to thermodynamic databases

In this paper we derive an approach for the effective utilization of thermodynamic data in phase-field simulations. While the most widely used methodology for multi-component alloys is following the work by Eiken et al. (2006), wherein, an extrapolative scheme is utilized in conjunction with the TQ interface for deriving the driving force for phase transformation, a corresponding simplistic method based on the formulation of a parabolic free-energy model incorporating all the thermodynamics has been laid out for binary alloys in the work by Folch and Plapp (2005). In the following, we extend this latter approach for multi-component alloys in the framework of the grand-potential formalism. The coupling is applied for the case of the binary eutectic solidification in the Cr-Ni alloy and two-phase solidification in the ternary eutectic alloy (Al-Cr-Ni). A thermodynamic justification entails the basis of the formulation and places it in context of the bigger picture of Integrated Computational Materials Engineering. (C) 2015 Elsevier Ltd. All rights reserved.

An evolutionary formalism for weak quantum measurements

Unitary evolution and projective measurement are fundamental axioms of quantum mechanics. Even though projective measurement yields one of the eigenstates of the measured operator as the outcome, there is no theory that predicts which eigenstate will be observed in which experimental run. There exists only an ensemble description, which predicts probabilities of various outcomes over many experimental runs. We propose a dynamical evolution equation for the projective collapse of the quantum state in individual experimental runs, which is consistent with the well-established framework of quantum mechanics. In case of gradual weak measurements, its predictions for ensemble evolution are different from those of the Born rule. It is an open question whether or not suitably designed experiments can observe this alternate evolution.

Structural and lithium ion transport studies in sol-gel-prepared lithium silicophosphate glasses

Lithium silicophosphate glasses have been prepared by a sol-gel route over a wide range of compositions. Their structural and electrical properties have been investigated. Infrared spectroscopic studies show the presence of hydroxyl groups attached to Si and P. MAS NMR investigations provide evidence for the presence of different phosphatic units in the structure. The variations of de conductivities at 423 K and activation energies have been studied as a function of composition, and both exhibit an increasing trend with the ratio of nonbridging oxygen to bridging oxygen in the structure. Ac conductivity behavior shows that the power law exponent, s, is temperature dependent and exhibits a minimum. Relaxation behavior has been examined in detail using an electrical modulus formalism, and modulus data were fitted to Kohlraush-William-Watts stretched exponential function. A structural model has been proposed and the unusual properties exhibited by this unique system of glasses have been rationalized using this model. Ion transport in these glasses appears to be confined to unidimensional conduits defined by modified phosphate chains and interspersed with unmodified silica units.

Superconductivity in the nonadiabatic limit

Isotope-shift exponent (cu) and the pressure coefficient of superconducting transition temperature (beta) have been studied in the nonadiabatic limit. We have considered the effect of nonadiabaticity in both within and beyond the Migdal-Eliashberg formalism. It reveals from our study that the pressure coefficient of superconducting transition is high for the low-T-c region and low for the high-T-c region and the minimum value of alpha is obtained where the transition temperature is maximum. Lowest value of isotope-shift exponent is obtained for small momentum exchange between the electrons and the bosonic field. Qualitative variation of beta with temperature is consistent with the experimental results of the hole doped superconductors for small momentum exchange.

Einstein Relation in n-Channel Inversion Layers of Nonlinear Optical, Optoelectronic and Related Materials: Simplified Theory, Relative Comparison and Suggestion for an Experimental Determination

In this paper, we study the Einstein relation for the diffusivity to mobility ratio (DMR) in n-channel inversion layers of non-linear optical materials on the basis of a newly formulated electron dispersion relation by considering their special properties within the frame work of k.p formalism. The results for the n-channel inversion layers of III-V, ternary and quaternary materials form a special case of our generalized analysis. The DMR for n-channel inversion layers of II-VI, IV-VI and stressed materials has been investigated by formulating the respective 2D electron dispersion laws. It has been found, taking n-channel inversion layers of CdGeAs2, Cd(3)AS(2), InAs, InSb, Hg1-xCdxTe, In1-xGaxAsyP1-y lattice matched to InP, CdS, PbTe, PbSnTe, Pb1-xSnxSe and stressed InSb as examples, that the DMR increases with the increasing surface electric field with different numerical values and the nature of the variations are totally band structure dependent. The well-known expression of the DMR for wide gap materials has been obtained as a special case under certain limiting conditions and this compatibility is an indirect test for our generalized formalism. Besides, an experimental method of determining the 2D DMR for n-channel inversion layers having arbitrary dispersion laws has been suggested.

Optimal barrier subdivision for Kramers' escape rate

We examine the effect of subdividing the potential barrier along the reaction coordinate on Kramers' escape rate for a model potential. Using the known supersymmetric potential approach, we show the existence of an optimal number of subdivisions that maximizes the rate.