Net Baryon Fluctuations from a Crossover Equation of State

We have constructed an equation of state which smoothly interpolates between an excluded volume hadron resonance gas at low energy density to a plasma of quarks and gluons at high energy density. This crossover equation of state agrees very well with lattice calculations at both zero and nonzero baryon chemical potential. We use it to compute the variance, skewness, and kurtosis of fluctuations of baryon number, and compare to measurements of proton number fluctuations in central Au-Au collisions as measured by the STAR collaboration in a beam energy scan at the Relativistic Heavy Ion Collider. The crossover equation of state can reproduce the data if the fluctuations are frozen out at temperatures well below than the average chemical freeze-out.Comment: 5 pages, 7 figures. arXiv admin note: substantial text overlap with arXiv:1506.0340

Matching Excluded Volume Hadron Resonance Gas Models and Perturbative QCD to Lattice Calculations

We match three hadronic equations of state at low energy densities to a perturbatively computed equation of state of quarks and gluons at high energy densities. One of them includes all known hadrons treated as point particles, which approximates attractive interactions among hadrons. The other two include, in addition, repulsive interactions in the form of excluded volumes occupied by the hadrons. A switching function is employed to make the crossover transition from one phase to another without introducing a thermodynamic phase transition. A chi-square fit to accurate lattice calculations with temperature $100 < T < 1000$ MeV determines the parameters. These parameters quantify the behavior of the QCD running gauge coupling and the hard core radius of protons and neutrons, which turns out to be $0.62 \pm 0.04$ fm. The most physically reasonable models include the excluded volume effect. Not only do they include the effects of attractive and repulsive interactions among hadrons, but they also achieve better agreement with lattice QCD calculations of the equation of state. The equations of state constructed in this paper do not result in a phase transition, at least not for the temperatures and baryon chemical potentials investigated. It remains to be seen how well these equations of state will represent experimental data on high energy heavy ion collisions when implemented in hydrodynamic simulations.Comment: 19 pages, 9 figure

Realization of the Large Mixing Angle Solar Neutrino Solution in an SO(10) Supersymmetric Grand Unified Model

An SO(10) supersymmetric grand unified model proposed earlier leading to the solar solution involving ``just-so'' vacuum oscillations is reexamined to study its ability to obtain the other possible solar solutions. It is found that all four viable solar neutrino oscillation solutions can be achieved in the model simply by modification of the right-handed Majorana neutrino mass matrix, M_R. Whereas the small mixing and vacuum solutions are easily obtained with several texture zeros in M_R, the currently-favored large mixing angle solution requires a nearly geometric hierarchical form for M_R that leads by the seesaw formula to a light neutrino mass matrix which has two or three texture zeros. The form of the matrix which provides the ``fine-tuning'' necessary to achieve the large mixing angle solution can be understood in terms of Froggatt-Nielsen diagrams for the Dirac and right-handed Majorana neutrino mass matrices. The solution fulfils several leptogenesis requirements which in turn can be responsible for the baryon asymmetry in the universe.Comment: 14 pages including 2 figure

Lifting a Realistic SO(10) Grand Unified Model to Five Dimensions

It has been shown recently that the problem of rapid proton decay induced by dimension five operators arising from the exchange of colored Higgsinos can be simply avoided in grand unified models where a fifth spatial dimension is compactified on an orbifold. Here we demonstrate that this idea can be used to solve the Higgsino-mediated proton decay problem in any realistic SO(10) model by lifting that model to five dimensions. A particular SO(10) model that has been proposed to explain the pattern of quark and lepton masses and mixings is used as an example. The idea is to break the SO(10) down to the Pati-Salam symmetry by the orbifold boundary conditions. The entire four-dimensional SO(10) model is placed on the physical SO(10) brane except for the gauge fields, the 45 and a single 10 of Higgs fields, which are placed in the five-dimensional bulk. The structure of the Higgs superpotential can be somewhat simplified in doing so, while the Yukawa superpotential and mass matrices derived from it remain essentially unaltered.Comment: 17 pages, version to be published in Phys. Rev. D with expanded discussion of the suppression of dim-5 proton decay operator

Resonant leptogenesis in a predictive SO(10) grand unified model

An SO(10) grand unified model considered previously by the authors featuring lopsided down quark and charged lepton mass matrices is successfully predictive and requires that the lightest two right-handed Majorana neutrinons be nearly degenerate in order to obtain the LMA solar neutrino solution. Here we use this model to test its predictions for baryogenesis through resonant-enhanced leptogenesis. With the conventional type I seesaw mechanism, the best predictions for baryogenesis appear to fall a factor of three short of the observed value. However, with a proposed type III seesaw mechanism leading to three pairs of massive pseudo-Dirac neutrinos, resonant leptogenesis is decoupled from the neutrino mass and mixing issues with successful baryogenesis easily obtained.Comment: 22 pages including 1 figure; published version with reference adde

Radiological Risks of Neutron Interrogation of Food

In recent years there has been growing interest in the use of neutron scanning techniques for security. Neutron techniques with a range of energy spectra including thermal, white and fast neutrons have been shown to work in different scenarios. As international interest in neutron scanning increases the risk of activating cargo, especially foodstuffs must be considered. There has been a limited amount of research into the activation of foods by neutron beams and we have sought to improve the amount of information available. In this paper we show that for three important metrics; Activity, Ingestion Dose and Time to Background there is a strong dependence on the food being irradiated and a weak dependence on the energy of irradiation. Previous studies into activation used results based on irradiation of pharmaceuticals as the basis for research into activation of food. The earlier work reports that 24Na production is the dominant threat which motivated the search for 23(n;\gamma)24Na in highly salted foods. We show that 42K can be more significant than 24Na in low salt foods such as Bananas and Potatoes

On the Numerical Dispersion of Electromagnetic Particle-In-Cell Code : Finite Grid Instability

The Particle-In-Cell (PIC) method is widely used in relativistic particle beam and laser plasma modeling. However, the PIC method exhibits numerical instabilities that can render unphysical simulation results or even destroy the simulation. For electromagnetic relativistic beam and plasma modeling, the most relevant numerical instabilities are the finite grid instability and the numerical Cherenkov instability. We review the numerical dispersion relation of the electromagnetic PIC algorithm to analyze the origin of these instabilities. We rigorously derive the faithful 3D numerical dispersion of the PIC algorithm, and then specialize to the Yee FDTD scheme. In particular, we account for the manner in which the PIC algorithm updates and samples the fields and distribution function. Temporal and spatial phase factors from solving Maxwell's equations on the Yee grid with the leapfrog scheme are also explicitly accounted for. Numerical solutions to the electrostatic-like modes in the 1D dispersion relation for a cold drifting plasma are obtained for parameters of interest. In the succeeding analysis, we investigate how the finite grid instability arises from the interaction of the numerical 1D modes admitted in the system and their aliases. The most significant interaction is due critically to the correct represenation of the operators in the dispersion relation. We obtain a simple analytic expression for the peak growth rate due to this interaction.Comment: 25 pages, 6 figure