Survey of vector-like fermion extensions of the Standard Model and their phenomenological implications

With the renewed interest in vector-like fermion extensions of the Standard Model, we present here a study of multiple vector-like theories and their phenomenological implications. Our focus is mostly on minimal flavor conserving theories that couple the vector-like fermions to the SM gauge fields and mix only weakly with SM fermions so as to avoid flavor problems. We present calculations for precision electroweak and vector-like state decays, which are needed to investigate compatibility with currently known data. We investigate the impact of vector-like fermions on Higgs boson production and decay, including loop contributions, in a wide variety of vector-like extensions and their parameter spaces.

Discrete-Rate Adaptation and Selection in Energy Harvesting Wireless Systems

In a system with energy harvesting (EH) nodes, the design focus shifts from minimizing energy consumption by infrequently transmitting less information to making the best use of available energy to efficiently deliver data while adhering to the fundamental energy neutrality constraint. We address the problem of maximizing the throughput of a system consisting of rate-adaptive EH nodes that transmit to a destination. Unlike related literature, we focus on the practically important discrete-rate adaptation model. First, for a single EH node, we propose a discrete-rate adaptation rule and prove its optimality for a general class of stationary and ergodic EH and fading processes. We then study a general system with multiple EH nodes in which one is opportunistically selected to transmit. We first derive a novel and throughput-optimal joint selection and rate adaptation rule (TOJSRA) when the nodes are subject to a weaker average power constraint. We then propose a novel rule for a multi-EH node system that is based on TOJSRA, and we prove its optimality for stationary and ergodic EH and fading processes. We also model the various energy overheads of the EH nodes and characterize their effect on the adaptation policy and the system throughput.

A matrix model for QCD: QCD color is mixed

We use general arguments to show that colored QCD states when restricted to gauge invariant local observables are mixed. This result has important implications for confinement: a pure colorless state can never evolve into two colored states by unitary evolution. Furthermore, the mean energy in such a mixed colored state is infinite. Our arguments are confirmed in a matrix model for QCD that we have developed using the work of Narasimhan and Ramadas(3) and Singer.(2) This model, a (0 + 1)-dimensional quantum mechanical model for gluons free of divergences and capturing important topological aspects of QCD, is adapted to analytical and numerical work. It is also suitable to work on large N QCD. As applications, we show that the gluon spectrum is gapped and also estimate some low-lying levels for N = 2 and 3 (colors). Incidentally the considerations here are generic and apply to any non-Abelian gauge theory.

Modelling ternary effects on antiphase boundary energy of Ni3Al

The shearing of ordered gamma' precipitates by matrix dislocations results in the formation of antiphase boundaries (APB) in Ni-base superalloys. The APB energy is an important source of order-strengthening in disk and blade alloys where Ti and Ta substitute for Al in gamma'. While the importance of APB energy is well-acknowledged, the effect of alloying on APB energy is not fully understood. In the present study, the effect of Ti and Ta additions on the {111} and {010} APB energies was probed via electronic structure calculations. Results suggest that at low levels of Ti/Ta, APB energies on either plane increases with alloying. However, at higher Ti/Ta levels, the APB energies decrease with alloying. These trends understood by accounting for nearest neighbour violations about the APB and additionally, invoking the effect of precipitate composition on the energy penalty of the violations. We propose an Environment Dependent Nearest Neighbour Bond (EDNNB) model that predicts APB energies that are in close agreement to calculated values.

A matrix model for QCD

Gribov's observation that global gauge fixing is impossible has led to suggestions that there may be a deep connection between gauge fixing and confinement. We find an unexpected relation between the topological nontriviality of the gauge bundle and colored states in SU(N) Yang-Mills theory, and show that such states are necessarily impure. We approximate QCD by a rectangular matrix model that captures the essential topological features of the gauge bundle, and demonstrate the impure nature of colored states explicitly. Our matrix model also allows the inclusion of the QCD theta-term, as well as to perform explicit computations of low-lying glueball masses. This mass spectrum is gapped. Since an impure state cannot evolve to a pure one by a unitary transformation, our result shows that the solution to the confinement problem in pure QCD is fundamentally quantum information-theoretic.

Maximum group velocity in a one-dimensional model with a sinusoidally varying staggered potential

We use Floquet theory to study the maximum value of the stroboscopic group velocity in a one-dimensional tight-binding model subjected to an on-site staggered potential varying sinusoidally in time. The results obtained by numerically diagonalizing the Floquet operator are analyzed using a variety of analytical schemes. In the low-frequency limit we use adiabatic theory, while in the high-frequency limit the Magnus expansion of the Floquet Hamiltonian turns out to be appropriate. When the magnitude of the staggered potential is much greater or much less than the hopping, we use degenerate Floquet perturbation theory; we find that dynamical localization occurs in the former case when the maximum group velocity vanishes. Finally, starting from an ``engineered'' initial state where the particles (taken to be hard-core bosons) are localized in one part of the chain, we demonstrate that the existence of a maximum stroboscopic group velocity manifests in a light-cone-like spreading of the particles in real space.

Phase diagram of the half-filled ionic Hubbard model

We study the phase diagram of the ionic Hubbard model (IHM) at half filling on a Bethe lattice of infinite connectivity using dynamical mean-field theory (DMFT), with two impurity solvers, namely, iterated perturbation theory (IPT) and continuous time quantum Monte Carlo (CTQMC). The physics of the IHM is governed by the competition between the staggered ionic potential Delta and the on-site Hubbard U. We find that for a finite Delta and at zero temperature, long-range antiferromagnetic (AFM) order sets in beyond a threshold U = U-AF via a first-order phase transition. For U smaller than U-AF the system is a correlated band insulator. Both methods show a clear evidence for a quantum transition to a half-metal (HM) phase just after the AFM order is turned on, followed by the formation of an AFM insulator on further increasing U. We show that the results obtained within both methods have good qualitative and quantitative consistency in the intermediate-to-strong-coupling regime at zero temperature as well as at finite temperature. On increasing the temperature, the AFM order is lost via a first-order phase transition at a transition temperature T-AF(U,Delta) or, equivalently, on decreasing U below U-AF(T,Delta)], within both methods, for weak to intermediate values of U/t. In the strongly correlated regime, where the effective low-energy Hamiltonian is the Heisenberg model, IPT is unable to capture the thermal (Neel) transition from the AFM phase to the paramagnetic phase, but the CTQMC does. At a finite temperature T, DMFT + CTQMC shows a second phase transition (not seen within DMFT + IPT) on increasing U beyond U-AF. At U-N > U-AF, when the Neel temperature T-N for the effective Heisenberg model becomes lower than T, the AFM order is lost via a second-order transition. For U >> Delta, T-N similar to t(2)/U(1 - x(2)), where x = 2 Delta/U and thus T-N increases with increase in Delta/U. In the three-dimensional parameter space of (U/t, T/t, and Delta/t), as T increases, the surface of first-order transition at U-AF(T,Delta) and that of the second-order transition at U-N(T,Delta) approach each other, shrinking the range over which the AFM order is stable. There is a line of tricritical points that separates the surfaces of first- and second-order phase transitions.

A Statistical Model of Current Loops and Magnetic Monopoles

We formulate a natural model of loops and isolated vertices for arbitrary planar graphs, which we call the monopole-dimer model. We show that the partition function of this model can be expressed as a determinant. We then extend the method of Kasteleyn and Temperley-Fisher to calculate the partition function exactly in the case of rectangular grids. This partition function turns out to be a square of a polynomial with positive integer coefficients when the grid lengths are even. Finally, we analyse this formula in the infinite volume limit and show that the local monopole density, free energy and entropy can be expressed in terms of well-known elliptic functions. Our technique is a novel determinantal formula for the partition function of a model of isolated vertices and loops for arbitrary graphs.

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.

Phase-field elasticity model based on mechanical jump conditions

Computational models based on the phase-field method typically operate on a mesoscopic length scale and resolve structural changes of the material and furthermore provide valuable information about microstructure and mechanical property relations. An accurate calculation of the stresses and mechanical energy at the transition region is therefore indispensable. We derive a quantitative phase-field elasticity model based on force balance and Hadamard jump conditions at the interface. Comparing the simulated stress profiles calculated with Voigt/Taylor (Annalen der Physik 274(12):573, 1889), Reuss/Sachs (Z Angew Math Mech 9:49, 1929) and the proposed model with the theoretically predicted stress fields in a plate with a round inclusion under hydrostatic tension, we show the quantitative characteristics of the model. In order to validate the elastic contribution to the driving force for phase transition, we demonstrate the absence of excess energy, calculated by Durga et al. (Model Simul Mater Sci Eng 21(5):055018, 2013), in a one-dimensional equilibrium condition of serial and parallel material chains. To validate the driving force for systems with curved transition regions, we relate simulations to the Gibbs-Thompson equilibrium condition

Homogeneous isotropic superfluid turbulence in two dimensions: Inverse and forward cascades in the Hall-Vinen-Bekharevich-Khalatnikov model

We present the first direct-numerical-simulation study of the statistical properties of two-dimensional superfluid turbulence in the simplified, Hall-Vinen-Bekharevich-Khalatnikov two-fluid model. We show that both normalfluid and superfluid energy spectra can exhibit two power-law regimes, the first associated with an inverse cascade of energy and the second with the forward cascade of enstrophy. We quantify the mutual-friction-induced alignment of normal and superfluid velocities by obtaining probability distribution functions of the angle between them and the ratio of their moduli.

A simplified model for dark matter interacting primarily with gluons

We consider a simple renormalizable model providing a UV completion for dark matter whose interactions with the Standard Model are primarily via the gluons. The model consists of scalar dark matter interacting with scalar colored mediator particles. A novel feature is the fact that (in contrast to more typical models containing dark matter whose interactions are mediated via colored scalars) the colored scalars typically decay into multi-quark final states, with no associated missing energy. We construct this class of models and examine associated phenomena related to dark matter annihilation, scattering with nuclei, and production at colliders.

Effect of Surface Energy Anisotropy on the Stability of Growth Fronts in Multiphase Alloys

Eutectic growth offers a variety of examples for pattern formation which are interesting both for theoreticians as well as experimentalists. One such example of patterns is ternary eutectic colonies which arise as a result of instabilities during growth of two solid phases. Here, in addition to the two major components being exchanged between the solid phases during eutectic growth, there is an impurity component which is rejected by both solid phases. During progress of solidification, there develops a boundary layer of the third impurity component ahead of the solidification front of the two solid phases. Similar to Mullins-Sekerka type instabilities, such a boundary layer tends to make the global solidification envelope unstable to morphological perturbations giving rise to two-phase cells. This phenomenon has been studied numerically in two dimensions for the conditions of directional solidification, by Plapp and Karma (Phys Rev E 66:061608, 2002) using phase-field simulations. While, in the work by Plapp and Karma (Phys Rev E 66:061608, 2002) all interfaces are isotropic, in our presentation, we extend the phase-field model by considering interfacial anisotropy in the solid-solid and solid-liquid interfaces and characterize the role of interfacial anisotropy on the stability of the growth front through phase-field simulations in two dimensions.

Estimation of activation energy for electroporation and pore growth rate in liquid crystalline and gel phases of lipid bilayers using molecular dynamics simulations

Molecular dynamics simulations of electroporation in POPC and DPPC lipid bilayers have been carried out at different temperatures ranging from 230 K to 350 K for varying electric fields. The dynamics of pore formation, including threshold field, pore initiation time, pore growth rate, and pore closure rate after the field is switched off, was studied in both the gel and liquid crystalline (L-alpha) phases of the bilayers. Using an Arrhenius model of pore initiation kinetics, the activation energy for pore opening was estimated to be 25.6 kJ mol(-1) and 32.6 kJ mol(-1) in the L-alpha phase of POPC and DPPC lipids respectively at a field strength of 0.32 V nm(-1). The activation energy decreases to 24.2 kJ mol(-1) and 23.7 kJ mol(-1) respectively at a higher field strength of 1.1 V nm(-1). At temperatures below the melting point, the activation energy in the gel phase of POPC and DPPC increases to 28.8 kJ mol(-1) and 34.4 kJ mol(-1) respectively at the same field of 1.1 V nm(-1). The pore closing time was found to be higher in the gel than in the L-alpha phase. The pore growth rate increases linearly with temperature and quadratically with field, consistent with viscosity limited growth models.

ASYMPTOTIC PRESERVING SCHEME FOR A KINETIC MODEL DESCRIBING IN COMPRESSIBLE FLUIDS

The kinetic theory of fluid turbulence modeling developed by Degond and Lemou in 7] is considered for further study, analysis and simulation. Starting with the Boltzmann like equation representation for turbulence modeling, a relaxation type collision term is introduced for isotropic turbulence. In order to describe some important turbulence phenomenology, the relaxation time incorporates a dependency on the turbulent microscopic energy and this makes difficult the construction of efficient numerical methods. To investigate this problem, we focus here on a multi-dimensional prototype model and first propose an appropriate change of frame that makes the numerical study simpler. Then, a numerical strategy to tackle the stiff relaxation source term is introduced in the spirit of Asymptotic Preserving Schemes. Numerical tests are performed in a one-dimensional framework on the basis of the developed strategy to confirm its efficiency.