A fair knowledge of the atmospheric muon distributions at Earth is a prerequisite for the simulations of cosmic ray setups and rare event search detectors. A modified power law is proposed for atmospheric muon energy distribution which gives a good description of the cosmic muon data in low as well as high energy regime. Using this distribution, analytical forms for zenith angle ($\theta$) distribution are obtained. Assuming a flat Earth, it leads to the $\cos^{n-1}\theta$ form where it is shown that the parameter $n$ is nothing but the power of the energy distribution. Exact analytical function is obtained for inclined trajectory of muon. A new closed form for zenith angle distribution is obtained without assuming a flat Earth and which gives an improved description of the data at all angles even above $70^o$. These distributions are tested with the available atmospheric muon data of energy and angular distributions. The parameters of these distributions can be used to characterize the cosmic muon data as a function of energy, angle and altitude.

We construct explicit examples of fibre inflation models which are globally embedded in type IIB orientifolds with chiral matter on D7-branes and full closed string moduli stabilisation. The minimal setup involves a Calabi-Yau threefold with h^{1,1}=4 Kaehler moduli which features multiple K3 fibrations and a del Pezzo divisor supporting non-perturbative effects. We perform a consistent choice of orientifold involution, brane setup and gauge fluxes which leads to chiral matter and a moduli-dependent Fayet-Iliopoulos term. After D-term stabilisation, the number of Kaehler moduli is effectively reduced to 3 and the internal volume reduces to the one of fibre inflation models. The inflationary potential is generated by suitable string loop corrections in combination with higher derivative effects. We analyse the inflationary dynamics both in the single-field approximation and by numerically deriving the full multi-field evolution in detail. Interestingly, we find that the Kaehler cone conditions set strong constraints on the allowed inflaton field range.

We argue that in type IIB LVS string models, after including the leading order moduli stabilisation effects, the moduli space for the remaining flat directions is compact due the Calabi-Yau K\"ahler cone conditions. In cosmological applications, this gives an inflaton field range which is bounded from above, in analogy with recent results from the weak gravity and swampland conjectures. We support our claim by explicitly showing that it holds for all LVS vacua with $h^{1,1} = 3$ obtained from 4-dimensional reflexive polytopes. In particular, we first search for all Calabi-Yau threefolds from the Kreuzer-Skarke list with $h^{1,1}=2$, $3$ and $4$ which allow for LVS vacua, finding several new LVS geometries which were so far unknown. We then focus on the $h^{1,1} = 3$ cases and show that the K\"ahler cones of all toric hypersurface threefolds force the effective 1-dimensional LVS moduli space to be compact. We find that the moduli space size can generically be trans-Planckian only for K3 fibred examples.

We describe global embeddings of fractional D3 branes at orientifolded singularities in type IIB flux compactifications. We present an explicit Calabi-Yau example where the chiral visible sector lives on a local orientifolded quiver while non-perturbative effects, $\alpha'$ corrections and a T-brane hidden sector lead to full closed string moduli stabilisation in a de Sitter vacuum. The same model can also successfully give rise to inflation driven by a del Pezzo divisor. Our model represents the first explicit Calabi-Yau example featuring both an inflationary and a chiral visible sector.

We present concrete embeddings of fibre inflation models in globally consistent type IIB Calabi-Yau orientifolds with closed string moduli stabilisation. After performing a systematic search through the existing list of toric Calabi-Yau manifolds, we find several examples that reproduce the minimal setup to embed fibre inflation models. This involves Calabi-Yau manifolds with $h^{1,1}= 3$ which are K3 fibrations over a $\mathbb{P}^1$ base with an additional shrinkable rigid divisor. We then provide different consistent choices of the underlying brane set-up which generate a non-perturbative superpotential suitable for moduli stabilisation and string loop corrections with the correct form to drive inflation. For each Calabi-Yau orientifold setting, we also compute the effect of higher derivative contributions and study their influence on the inflationary dynamics.

By demanding the validity of an effective field theory description during inflation, in this note we derive some peculiar inequalities among the three interesting stringy and cosmological parameters, namely the tensor-to-scalar ratio ($r$), the string coupling ($g_s$) and the compactification volume (${\cal V}$). In deriving these inequalities, we explicitly demand that the inflationary scale and the Hubble parameter during inflation are well below the Kaluza-Klein (KK) mass scale, string scale, and the four dimensional Planck mass. For the inflationary models developed within the framework of type IIB orientifold comapctification, we investigate the regions of parameters space spanned by the three parameters $(r, g_s, {\cal V})$ by satisfying our inequalities, and we find that the same can reduce the size of available parameter space quite significantly. Moreover, we comment on obtaining further constraints on the parameters by comparing gravitino mass ($m_{3/2}$) with the Hubble scale ($H$), which also provides a lower bound on tensor-to-scalar ratio ($r$), for the cases when $m_{3/2} <H$. We also illustrate the outcome of our bounds in some specific class of string(-inspired) models.

In the lights of current BICEP2 observations accompanied with the PLANCK satellite results, it has been observed that the simple single field chaotic inflationary models provide a good agreement with their spectral index n_s and large tensor-to-scalar ratio r (0.15 <r <0.26). To explore the other simple models, we consider the fractional-chaotic inflationary potentials of the form V_0 phi^(a/b) where a and b are relatively prime. We show that such kind of inflaton potentials can be realized elegantly in the supergravity framework with generalized shift symmetry and a nature bound a/b < 4 for consistency. Especially, for the number of e-folding from 50 to 60 and some a/b from 2 to 3, our predictions are nicely within at least 1 $\sigma$ region in the r-n_s plane. We also present a systematic investigation of such chaotic inflationary models with fractional exponents to explore the possibilities for the enhancement in the magnitude of running of spectral index (\alpha_{n_s}) beyond the simplistic models.

We successfully embed the Kim-Nilles-Peloso (KNP) alignment mechanism for enhancing the axion decay constant in the context of large volume type IIB orientifolds. The flat direction is generated in the plane of ($C_0$-$C_2$) axions corresponding to the involutively even universal axion $C_0$ and odd axion $C_2$, respectively. The moduli stabilization with large volume scheme has been established as well.

Using $\delta N$ formalism, in the context of a generic multi-field inflation driven on a non-flat field space background, we revisit the analytic expressions of the various cosmological observables such as scalar/tensor power spectra, scalar/tensor spectral tilts, non-Gaussianity parameters, tensor-to-scalar ratio, and the various runnings of these observables. In our backward formalism approach, the subsequent expressions of observables automatically include terms beyond the leading order slow-roll expansion correcting many of the expression at subleading order. To connect our analysis properly with the earlier results, we rederive the (well) known (single field) expressions in the limiting cases of our generic formulae. Further, in the light of PLANCK results, we examine for the compatibility of the consistency relations within the slow-roll regime of a two-field roulette poly-instanton inflation realized in the context of large volume scenarios.

In the context of Type IIB LARGE volume orientifold setup equipped with poly-instanton corrections, the standard single-field poly-instanton inflation driven by a 'Wilson' divisor volume modulus is generalized by the inclusion of respective axion modulus. This two-field dynamics results in a "Roulette" type inflation with the presence of several inflationary trajectories which could produce 50 (or more) e-foldings. The evolution of various trajectories along with physical observables are studied. The possibility of generating primordial non-Gaussianities in the slow-roll as well as in the beyond slow-roll region is investigated. We find that although the non- linearity parameters are quite small during the slow-roll regime, the same are significantly enhanced in the beyond slow-roll regime investigated up to the end of inflation.

The generation of Rossby rogue waves (Rossby rogons), as well as the excitation of bright and dark Rossby envelpe solitons are demonstrated on the basis of the modulational instability (MI) of a coherent Rossby wave packet. The evolution of an amplitude modulated Rossby wave packet is governed by one-dimensional (1D) nonlinear Schr\"odinger equation (NLSE). The latter is used to study the amplitude modulation of Rossby wave packets for fluids in Earth's atmosphere and in the solar photosphere. It is found that an ampitude modulated Rossby wave packet becomes stable (unstable) against quasi-stationary, long wavelength (in comparision with the Rossby wave length) perturbations, when the carrier Rossby wave number satisfies $k^2 < 1/2$ or $\sqrt{2}+1<k^2<3$ ($k^2 >3$ or $1/2<k^2<\sqrt{2}+1$). It is also shown that a Rossby rogon or a bright Rossby envelope soliton may be excited in the shallow water approximation for the Rossby waves in solar photosphere. However, the excitation of small or large scale perturbations may be possible for magnetized plasmas in the ionosphereic $E-$layer.

We present a generalized nonlinear theory for large amplitude electrostatic (ES) ion shocks in collisional quantum plasmas composed of mildly coupled degenerate electron fluids of arbitrary degeneracy and non-degenerate strongly correlated ion fluids with arbitrary atomic number. For our purposes, we use the inertialess electron momentum equation including the electrostatic, pressure gradient and relevant quantum forces, as well as a generalized viscoelastic momentum (GVEM) equation for strongly correlated non-degenerate ions. The ion continuity equation, in the quasi-neutral approximation, then closes our nonlinear system of equations. When the electric field is eliminated from the GVEM equation by using the inertialess electron momentum equation, we then obtain a generalized GVEM and the ion continuity equations exhibiting nonlinear couplings between the ion number density and the ion fluid velocity. The pair of nonlinear equations is numerically solved to study the dynamics of arbitrary large amplitudes planar and non-planar ES shocks arising from the balance between harmonic generation nonlinearities and the ion fluid viscosity for a wide range of the plasma mass-density and the ion atomic-number that are relevant for the cores of giant planets (viz. the Jupiter) and compact stars (viz. white dwarfs). Our numerical results reveal that the ES shock density profiles strongly depend on the plasma number density and composition (the atomic-number) parameters. Furthermore, density perturbations propagate with Mach numbers which significantly dependent on the studied plasma fractional parameters. It is concluded that the dynamics of the ES shocks in the super-dense degenerate plasma is quite different in the core of a white dwarf star from that in the lower density crust region.

It is shown that a strongly magnetized isothermal pair plasma near the surface of a pulsar supports low-frequency (in comparison to electron cyclotron frequency) toroidal electrostatic plasma modes in the equatorial region. Physically, the thermal pressure coupled with the magnetic pressure creates the low frequency oscillations which may grow for particular case of inhomogeneities of the equilibrium magnetic field and the pair plasma density.

In the context of Type IIB compactified on a large volume Swiss-Cheese orientifold in the presence of a mobile space-time filling D3-brane and stacks of fluxed D7-branes wrapping the "big" divisor Sigma_B of a Swiss-Cheese Calabi Yau in WCP^4 [1,1,1,6,9], we explore various implications of moduli dynamics and discuss their couplings and decay into MSSM (-like) matter fields early in the history of universe to reach thermal equilibrium. Like finite temperature effects in O'KKLT, we observe that the local minimum of zero-temperature effective scalar potential is stable against any finite temperature corrections (up to two-loops) in large volume scenarios as well. Also, we find that moduli are heavy enough to avoid any cosmological moduli problem.

Developing on the ideas of (section 4 of) [1] and [2] and using the formalisms of [3] and [4], after inclusion of perturbative and non-perturbative alpha' corrections to the Kaehler potential and (D1- and D3-) instanton generated superpotential, we show the possibility of getting finite values for the non-linear parameter f_{NL} while looking for non-Gaussianities in type IIB compactifications on orientifolds of the Swiss Cheese Calabi-Yau WCP^4[1,1,1,6,9] in the L(arge) V(olume) S(cenarios) limit. First we show that in the context of multi-field slow-roll inflation, for the Calabi-Yau volume V~10^5 and D3-instanton number n^s~10 along with N_e~18, one can realize f_{NL}~0.03, and for Calabi-Yau volume V~10^6 with D3-instanton number n^s~10 resulting in number of e-foldings N_e~60, one can realize f_{NL}~0.01. Further we show that with the slow-roll conditions violated and for the number of the D3-instanton wrappings n^s ~O(1), one can realize f_{NL} ~O(1). Using general considerations and some algebraic geometric assumptions, we show that with requiring a "freezeout" of curvature perturbations at super horizon scales, it is possible to get tensor-scalar ratio r ~O(10^{-3}) with the loss of scale invariance |n_R-1|=0.01 and one can obtain f_{NL} ~O(10^{-2}) as well in the context of slow-roll inflation scenarios in the same Swiss-Cheese setup. For all our calculations the large values of the genus-zero Gopakumar-Vafa invariants for compact projective varieties, play in important role. We also make some observations pertaining to the possibility of the axionic inflaton also being a cold dark matter candidate as well as a quintessence field used for explaining dark energy.

Continuing with the ideas of (section 4 of) [1], after inclusion of perturbative and non-perturbative alpha' corrections to the Kaehler potential and (D1- and D3-) instanton generated superpotential, we show the possibility of slow-roll axionic inflation in the large volume limit of Swiss-Cheese Calabi-Yau orientifold compactifications of type IIB string theory. We also include one- and two-loop corrections to the Kaehler potential but find the same to be subdominant to the (perturbative and non-perturbative) alpha' corrections. The NS-NS axions provide a flat direction for slow-roll inflation to proceed from a saddle point to the nearest dS minimum.

The nonlinear propagation of a circularly polarized electromagnetic (CPEM) wave in a strongly magnetized electron-positron-ion plasma is investigated. Two coupled equations describing the interaction between a high-frequency CPEM wave and the low-frequency electrostatic wake field are derived. It is found that the generation of the wake fields depends on the presence of the ion species. The wake field generation in turn leads to de-acceleration and frequency down conversion of the electromagnetic pulse.

Supernova remnant (SNR) blast shells can reach the flow speed $v_s = 0.1 c$ and shocks form at its front. Instabilities driven by shock-reflected ion beams heat the plasma in the foreshock, which may inject particles into diffusive acceleration. The ion beams can have the speed $v_b \approx v_s$. For $v_b \ll v_s$ the Buneman or upper-hybrid instabilities dominate, while for $v_b \gg v_s$ the filamentation and mixed modes grow faster. Here the relevant waves for $v_b \approx v_s$ are examined and how they interact nonlinearly with the particles. The collision of two plasma clouds at the speed $v_s$ is modelled with particle-in-cell (PIC) simulations, which convect with them magnetic fields oriented perpendicular to their flow velocity vector. One simulation models equally dense clouds and the other one uses a density ratio of 2. Both simulations show upper-hybrid waves that are planar over large spatial intervals and that accelerate electrons to $\sim$ 10 keV. The symmetric collision yields only short oscillatory wave pulses, while the asymmetric collision also produces large-scale electric fields, probably through a magnetic pressure gradient. The large-scale fields destroy the electron phase space holes and they accelerate the ions, which facilitates the formation of a precursor shock.

The quantum effect on the Weibel instability in an unmagnetized plasma is presented. Our analysis shows that the quantum effect tends to stabilize the Weibel instability in the hydrodynamic regime, whereas it produces a new oscillatory instability in the kinetic regime. A novel effect the quantum damping, which is associated with the Landau damping, is disclosed. The new quantum Weibel instability may be responsible for the generation of non-stationary magnetic fields in compact astrophysical objects as well as in the forthcoming intense laser-solid density plasma experiments.

We present for the first time the nonlinear dynamics of quantum electrodynamic (QED) photon splitting in a strongly magnetized electron-positron (pair) plasma. By using a QED corrected Maxwell equation, we derive a set of equations that exhibit nonlinear couplings between electromagnetic (EM) waves due to nonlinear plasma currents and QED polarization and magnetization effects. Numerical analyses of our coupled nonlinear EM wave equations reveal the possibility of a more efficient decay channel, as well as new features of energy exchange among the three EM modes that are nonlinearly interacting in magnetized pair plasmas. Possible applications of our investigation to astrophysical settings, such as magnetars, are pointed out.

We model relativistically colliding plasma by PIC simulations in one and two spatial dimensions, taking an ion-to-electron mass ratio of 400. Energy dissipation by a wave precursor of mixed polarity and different densities of the colliding plasma slabs results in a relativistic forward shock forming on millisecond timescales. The forward shock accelerates electrons to ultrarelativistic energies and reflects upstream ions, which drag the electrons along to preserve the plasma quasi-neutrality. No reverse shock forms. The shock may be representative for internal gamma ray burst shocks.

Nonlinear interactions involving electrostatic upper-hybrid (UH), ion-cyclotron (IC), lower-hybrid (LH), and Alfven waves in quantum magnetoplasmas are considered. For this purpose, the quantum hydrodynamical equations are used to derive the governing equations for nonlinearly coupled UH, IC, LH, and Alfven waves. The equations are then Fourier analyzed to obtain nonlinear dispersion relations, which admit both decay and modulational instabilities of the UH waves at quantum scales. The growth rates of the instabilities are presented. They can be useful in applications of our work to diagnostics in laboratory and astrophysical settings.

This paper is withdrawn because it is a duplicate of astro-ph/0601311

Particle acceleration by means of non-linear plasma wave interactions is of great topical interest. Accordingly, in this paper we focus on the electron surfing process. Self-consistent kinetic simulations, using both relativistic Vlasov and PIC (Particle In Cell) approaches, show here that electrons can be accelerated to highly relativistic energies (up to 100 m_e c^2) if the phase speed of the electrostatic wave is mildly relativistic (0.6c to 0.9c for the magnetic field strengths considered). The acceleration is strong because of relativistic stabilisation of the nonlinearly saturated electrostatic wave, seen in both relativistic Vlasov and PIC simulations. An inverse power law momentum distribution can arise for the most strongly accelerated electrons. These results are of relevance to observed rapid changes in the radio synchrotron emission intensities from microquasars, gamma ray bursts and other astrophysical objects that require rapid acceleration mechanisms for electrons.

A dispersion relation for electromagnetic wave propagation in a strongly magnetized cold plasma is deduced, taking photon-photon scattering into account. It is shown that the combined plasma and quantum electrodynamic effect is important for understanding the mode-structures in magnetar and pulsar atmospheres. The implications of our results are discussed.

The nonlinear propagation of low-frequency circularly polarized waves in a magnetized dusty plasma is analyzed. It is found that wave steepening and shock formation can take place due to the presence of nonlinear quantum vacuum effects, thus giving rise to ultra-intense electromagnetic shocks. Moreover, it is shown that solitary wave structures are admitted even under moderate astrophysical conditions. The results may have applications to astrophysical plasmas, as well as next generation laser interactions with laboratory plasmas containing dust clusters.

A new nonlinear electromagnetic wave mode in a magnetized dusty plasma is predicted. Its existence depends on the interaction of an intense circularly polarized electromagnetic wave with a dusty plasma, where quantum electrodynamical photon-photon scattering is taken into account. Specifically, we consider a dusty electron-positron-ion plasma, and show that the propagation of the new mode is admitted. It could be of significance for the physics of supernova remnants and in neutron star formation.

A new nonlinear electromagnetic wave mode in a magnetized plasma is predicted. Its existence depends on the interaction of an intense circularly polarized electromagnetic wave with a plasma, where quantum electrodynamical photon--photon scattering is taken into account. This scattering gives rise to a new coupling between the matter and the radiation. Specifically, we consider an electron--positron plasma, and show that the propagation of the new mode is admitted. It could be of significance in pulsar magnetospheres, and result in energy transport between the pulsar poles.

The existence of the dust acoustic wave (DAW) in a strongly magnetized electron-positron (pair)-dust plasma is demonstrated. In the DAW, the restoring force comes from the pressure of inertialess electrons and positrons, and the dust mass provides the inertia. The waves could be of interest in astrophysical settings such as the supernovae and pulsars, as well as in cluster explosions by intense laser beams in laboratory plasmas.

It is well-known that a charged particle moving with constant velocity in vacuum does not radiate. In a medium the situation can be different. If the so called Cherenkov condition is satisfied, i.e. the particle velocity exceeds the phase speed in the medium, the particle will radiate. We show that a charge moving with a constant velocity in a gas of photons emits Cherenkov radiation, even in the gamma-ray regime, due to nonlinear quantum electrodynamic effects. Our result is evaluated with respect to the radiation background in the early universe, and it is argued that the effect can be significant.

We present a novel nonlinear mechanism for exciting a gravitational radiation pulse (or a gravitational wave) by dust magnetohydrodynamic (DMHD) waves in dusty astrophysical plasmas. We derive the relevant equations governing the dynamics of nonlinearly coupled DMHD waves and a gravitational wave (GW). The system of equations is used to investigate the generation of a GW by compressional Alfv\'{e}n waves in a type II supernova. The growth rate of our nonlinear process is estimated, and the results are discussed in the context of the gravitational radiation accompanying supernova explosions.

The nonlinear propagation of intense incoherent photons in a photon gas is considered. The photon-photon interactions are governed by a pair of equations comprising a wave-kinetic equation for the incoherent photons in the presence of the slowly varying energy density perturbations of sound-like waves, and an equation for the latter waves in a background where the photon coupling is caused by quantum electrodynamical effects. The coupled equations are used to derive a dispersion relation, which admits new classes of modulational instabilities of incoherent photons. The present instabilities can lead to fragmentation of broadband short photon pulses in astrophysical and laboratory settings.

We analyze the collective behavior of neutrinos and antineutrinos in a dense background. Using the Wigner transform technique, it is shown that the interaction can be modelled by a coupled system of nonlinear Vlasov-like equations. From these equations, we derive a dispersion relation for neutrino-antineutrino interactions on a general background. The dispersion relation admits a novel modulational instability. The results are examined, together with a numerical example, and we discuss the induced density inhomogeneities using parameters relevant to the early Universe.

It is shown that the collective nonlinear interactions between intense neutrino or anti-neutrino fluxes and a dense neutrino plasma are governed by a multi-dimensional coupled cubic Schr\"odinger equation in which the interaction potential is positive or negative depending on the neutrino type. The cubic Schr\"odinger equation describes the splitting and focusing of intense neutrino beams due to the nonlinear excitations associated with the modifications of the individual neutrino energies in a dense neutrino background.

The derivation of the ponderomotive force due to neutrinos by Hardy and Melrose, Phys.Rev.D v.54, p.6491 (1996) is based on a flawed analysis of the ponderomotive force concept. Their conclusions also contain an erroneous physical assumption related to the neutrino emission in supernovae. A correct analysis shows the importance of the ponderomotive force due to neutrinos in type II supernovae explosions.