Atmospheric Secondary Particles In Near Earth Space

The Alpha Magnetic Spectrometer detects a large amount of particles below rigidity cutoff. Those high energy particles create questions related to radiation belts and atmospheric neutrinos. To understand the origin of these particles, we use a trajectory tracing program to simulate particle trajectories in realistic geomagnetic field. The complex behaviors and large e^+/e^- are explained here.Comment: 6 pages, 5 figures; Submit to the 8th Asia Pacific Physics Conferenc

Hybrid moments of the Riemann zeta-function

The "hybrid" moments $$\int_T^{2T}|\zeta(1/2+it)|^k{(\int_{t-G}^{t+G}|\zeta(1/2+ix)|^\ell dx)}^m dt$$ of the Riemann zeta-function $\zeta(s)$ on the critical line $\Re s = 1/2$ are studied. The expected upper bound for the above expression is $O_\epsilon(T^{1+\epsilon}G^m)$. This is shown to be true for certain specific values of the natural numbers $k,\ell,m$, and the explicitly determined range of $G = G(T;k,\ell,m)$. The application to a mean square bound for the Mellin transform function of $|\zeta(1/2+ix)|^4$ is given.Comment: 27 page

Electroweak two-loop corrections to sin^2{\theta}(eff,bb) and R(b) using numerical Mellin-Barnes integrals

Multi-loop integrals can be evaluated numerically using Mellin-Barnes representations. Here this technique is applied to the calculation of electroweak two-loop correction with closed fermion loops for two observables: the effective weak mixing angle for bottom quarks, sin^2{\theta}(eff,bb), and the branching ratio of the Z boson into bottom quarks, R(b). Good agreement with a previous result for sin^2{\theta}(eff,bb) is found. The result for R(b) is new, and a simple parametrization formula is provided which approximates the full result within integration errors.Comment: 12 pages, 2 figures. V2: typos in eqs. (9)-(14) fixed; error in the treatment of the running bottom-quark mass corrected, resulting in small shifts of the results in Tabs. 3,4, Fig. 2, and eqs. (21),(23); conclusions unchanged. V3: bug in calculation of R(b) corrected, resulting in a sizable reduction of the size of the two-loop correctio

On the Distribution of Plasmoids In High-Lundquist-Number Magnetic Reconnection

The distribution function $f(\psi)$ of magnetic flux $\psi$ in plasmoids formed in high-Lundquist-number current sheets is studied by means of an analytic phenomenological model and direct numerical simulations. The distribution function is shown to follow a power law $f(\psi)\sim\psi^{-1}$, which differs from other recent theoretical predictions. Physical explanations are given for the discrepant predictions of other theoretical models.Comment: Accepted for publication in Phys. Rev. Let

Chromomagnetic instability in dense quark matter

The results for the Debye and Meissner screening masses of the gluons and the photon in the case of neutral and beta-equilibrated dense two-flavor quark matter are presented. In the limits of the normal phase and the ideal two-flavor color superconducting phase, the screening masses coincide with the known results. Most interestingly, we find that the Meissner screening masses squared can be negative, indicating a plasma-type chromomagnetic instability in dense quark matter.Comment: 5 pages and 2 figures. To appear in Phys. Rev. D; version slightly shortened on request of editor

Penalized Likelihood Methods for Estimation of Sparse High Dimensional Directed Acyclic Graphs

Directed acyclic graphs (DAGs) are commonly used to represent causal relationships among random variables in graphical models. Applications of these models arise in the study of physical, as well as biological systems, where directed edges between nodes represent the influence of components of the system on each other. The general problem of estimating DAGs from observed data is computationally NP-hard, Moreover two directed graphs may be observationally equivalent. When the nodes exhibit a natural ordering, the problem of estimating directed graphs reduces to the problem of estimating the structure of the network. In this paper, we propose a penalized likelihood approach that directly estimates the adjacency matrix of DAGs. Both lasso and adaptive lasso penalties are considered and an efficient algorithm is proposed for estimation of high dimensional DAGs. We study variable selection consistency of the two penalties when the number of variables grows to infinity with the sample size. We show that although lasso can only consistently estimate the true network under stringent assumptions, adaptive lasso achieves this task under mild regularity conditions. The performance of the proposed methods is compared to alternative methods in simulated, as well as real, data examples.Comment: 19 pages, 8 figure

Effects of the Shear Viscosity on the Character of Cosmological Evolution

Bianchi type I cosmological models are studied that contain a stiff fluid with a shear viscosity that is a power function of the energy density, such as $\zeta = \alpha \epsilon^n$. These models are analyzed by describing the cosmological evolutions as the trajectories in the phase plane of Hubble functions. The simple and exact equations that determine these flows are obtained when $n$ is an integer. In particular, it is proved that there is no Einstein initial singularity in the models of $0\leq n < 1$. Cosmologies are found to begin with zero energy density and in the course of evolution the gravitational field will create matter. At the final stage, cosmologies are driven to the isotropic Fnedmann universe. It is also pointed out that although the anisotropy will always be smoothed out asymptotically, there are solutions that simultaneously possess non-positive and non-negative Hubble functions for all time. This means that the cosmological dimensional reduction can work even if the matter fluid having shear viscosity. These characteristics can also be found in any-dimensional models