On Blow-up Profile of Ground States of Boson Stars with External Potential

We study minimizers of the pseudo-relativistic Hartree functional \({\mathcal {E}}_{a}(u):=\Vert (-\varDelta +m^{2})^{1/4}u\Vert _{L^{2}}^{2}+\int _{{\mathbb {R}}^{3}}V(x)|u(x)|^{2}\mathrm{d}x-\frac{a}{2}\int _{{\mathbb {R}}^{3}}(\left| \cdot \right| ^{-1}\star |u|^{2})(x)|u(x)|^{2}\mathrm{d}x\) under the mass constraint \(\int _{{\mathbb {R}}^3}|u(x)|^2\mathrm{d}x=1\). Here \(m>0\) is the mass of particles and \(V\ge 0\) is an external potential. We prove that minimizers exist if and only if a satisfies \(0\le a<a^{*}\), and there is no minimizer if \(a\ge a^*\), where \(a^*\) is called the Chandrasekhar limit. When a approaches \(a^*\) from below, the blow-up behavior of minimizers is derived under some general external potentials V. Here we consider three cases of V: trapping potential, i.e. \(V\in L_{\mathrm{loc}}^{\infty }({\mathbb {R}}^3)\) satisfies \(\lim _{|x|\rightarrow \infty }V(x)=\infty \); periodic potential, i.e. \(V\in C({\mathbb {R}}^3)\) satisfies \(V(x+z)=V(x)\) for all \(z\in \mathbb {Z}^3\); and ring-shaped potential, e.g. \( V(x)=||x|-1|^p\) for some \(p>0\).     

Uniform in N global well-posedness of the time-dependent Hartree–Fock–Bogoliubov equations in $$\mathbb {R}^{1+1}$$

We prove the global well-posedness of the time-dependent Hartree–Fock–Bogoliubov (TDHFB) equations in \(\mathbb {R}^{1+1}\) with two-body interaction potential of the form \(N^{-1}v_N(x) = N^{\beta -1} v(N^\beta x)\) where \(v\ge 0\) is a sufficiently regular radial function, i.e., \(v \in L^1(\mathbb {R})\cap C^\infty (\mathbb {R})\). In particular, using methods of dispersive PDEs similar to the ones used in Grillakis and Machedon (Commun Partial Differ Equ 42:24–67, 2017), we are able to show for any scaling parameter \(\beta >0\) the TDHFB equations are globally well-posed in some Strichartz-type spaces independent of N, cf. (Bach et al. in The time-dependent Hartree–Fock–Bogoliubov equations for Bosons, 2016. arXiv:1602.05171).     

Reexploration of interacting holographic dark energy model: cases of interaction term excluding the Hubble parameter

In this paper, we make a deep analysis for the five typical interacting holographic dark energy models with the interaction terms \(Q=3\beta H_{0}\rho _\mathrm{{de}}\), \(Q=3\beta H_{0}\rho _\mathrm{{c}}\), \(Q=3\beta H_{0}(\rho _\mathrm{{de}}+\rho _\mathrm{c})\), \(Q=3\beta H_{0}\sqrt{\rho _\mathrm{{de}}\rho _\mathrm{c}}\), and \(Q=3\beta H_{0}\frac{\rho _\mathrm{{de}}\rho _{c}}{\rho _\mathrm{{de}}+\rho _\mathrm{c}}\), respectively. We obtain observational constraints on these models by using the type Ia supernova data (the Joint Light-Curve Analysis sample), the cosmic microwave background data (Planck 2015 distance priors), the baryon acoustic oscillations data, and the direct measurement of the Hubble constant. We find that the values of \(\chi _\mathrm{min}^2\) for all the five models are almost equal (around 699), indicating that the current observational data equally favor these IHDE models. In addition, a comparison with the cases of an interaction term involving the Hubble parameter H is also made.     

The decays of $$\bar{B}^0$$, $$\bar{B}^0_s$$B¯s0 and $$B^-$$B- into $$\eta _c$$ηc plus a scalar or vector meson

We investigate the decays of \(\bar{B}^0_s\), \(\bar{B}^0\) and \(B^-\) into \(\eta _c\) plus a scalar or vector meson in a theoretical framework by taking into account the dominant process for the weak decay of \(\bar{B}\) meson into \(\eta _c\) and a \(q\bar{q}\) pair. After hadronization of this \(q\bar{q}\) component into pairs of pseudoscalar mesons we obtain certain weights for the pseudoscalar meson-pseudoscalar meson components. In addition, the \(\bar{B}^0\) and \(\bar{B}^0_s\) decays into \(\eta _c\) and \(\rho ^0\), \(K^*\) are evaluated and compared to the \(\eta _c\) and \(\phi \) production. The calculation is based on the postulation that the scalar mesons \(f_0(500)\), \(f_0(980)\) and \(a_0(980)\) are dynamically generated states from the pseudoscalar meson-pseudoscalar meson interactions in S-wave. Up to a global normalization factor, the \(\pi \pi \), \(K \bar{K}\) and \(\pi \eta \) invariant mass distributions for the decays of \(\bar{B}^0_s \rightarrow \eta _c \pi ^+ \pi ^-\), \(\bar{B}^0_s \rightarrow \eta _c K^+ K^-\), \(\bar{B}^0 \rightarrow \eta _c \pi ^+ \pi ^-\), \(\bar{B}^0 \rightarrow \eta _c K^+ K^-\), \(\bar{B}^0 \rightarrow \eta _c \pi ^0 \eta \), \(B^- \rightarrow \eta _c K^0 K^-\) and \(B^- \rightarrow \eta _c \pi ^- \eta \) are predicted. Comparison is made with the limited experimental information available and other theoretical calcualtions. Further comparison of these results with coming LHCb measurements will be very valuable to make progress in our understanding of the nature of the low lying scalar mesons, \(f_0(500), f_0(980)\) and \(a_0(980)\).     

$$\varvec{B^0\rightarrow K^{*0}\mu ^+\mu ^-}$$ decay in the aligned two-Higgs-doublet model

In the aligned two-Higgs-doublet model, we perform a complete one-loop computation of the short-distance Wilson coefficients \(C_{7,9,10}^{(\prime )}\), which are the most relevant ones for \(b\rightarrow s\ell ^+\ell ^-\) transitions. It is found that, when the model parameter \(\left| \varsigma _{u}\right| \) is much smaller than \(\left| \varsigma _{d}\right| \), the charged scalar contributes mainly to chirality-flipped \(C_{9,10}^\prime \), with the corresponding effects being proportional to \(\left| \varsigma _{d}\right| ^2\). Numerically, the charged-scalar effects fit into two categories: (A) \(C_{7,9,10}^\mathrm {H^\pm }\) are sizable, but \(C_{9,10}^{\prime \mathrm {H^\pm }}\simeq 0\), corresponding to the (large \(\left| \varsigma _{u}\right| \), small \(\left| \varsigma _{d}\right| \)) region; (B) \(C_7^\mathrm {H^\pm }\) and \(C_{9,10}^{\prime \mathrm {H^\pm }}\) are sizable, but \(C_{9,10}^\mathrm {H^\pm }\simeq 0\), corresponding to the (small \(\left| \varsigma _{u}\right| \), large \(\left| \varsigma _{d}\right| \)) region. Taking into account phenomenological constraints from the inclusive radiative decay \(B\rightarrow X_{s}{\gamma }\), as well as the latest model-independent global analysis of \(b\rightarrow s\ell ^+\ell ^-\) data, we obtain the much restricted parameter space of the model. We then study the impact of the allowed model parameters on the angular observables \(P_2\) and \(P_5'\) of \(B^0\rightarrow K^{*0}\mu ^+\mu ^-\) decay, and we find that \(P_5'\) could be increased significantly to be consistent with the experimental data in case B.     

Universality for 1d Random Band Matrices: Sigma-Model Approximation

The paper continues the development of the rigorous supersymmetric transfer matrix approach to the random band matrices started in (J Stat Phys 164:1233–1260, 2016; Commun Math Phys 351:1009–1044, 2017). We consider random Hermitian block band matrices consisting of \(W\times W\) random Gaussian blocks (parametrized by \(j,k \in \Lambda =[1,n]^d\cap \mathbb {Z}^d\)) with a fixed entry’s variance \(J_{jk}=\delta _{j,k}W^{-1}+\beta \Delta _{j,k}W^{-2}\), \(\beta >0\) in each block. Taking the limit \(W\rightarrow \infty \) with fixed n and \(\beta \), we derive the sigma-model approximation of the second correlation function similar to Efetov’s one. Then, considering the limit \(\beta , n\rightarrow \infty \), we prove that in the dimension \(d=1\) the behaviour of the sigma-model approximation in the bulk of the spectrum, as \(\beta \gg n\), is determined by the classical Wigner–Dyson statistics.     

Convergence in High Probability of the Quantum Diffusion in a Random Band Matrix Model

We consider Hermitian random band matrices H in \(d \geqslant 1 \) dimensions. The matrix elements \(H_{xy},\) indexed by \(x, y \in \varLambda \subset \mathbb {Z}^d,\) are independent, uniformly distributed random variable if \(|x-y| \) is less than the band width W,  and zero otherwise. We update the previous results of the converge of quantum diffusion in a random band matrix model from convergence of the expectation to convergence in high probability. The result is uniformly in the size \(|\varLambda | \) of the matrix.     

Lifshits Tails for Randomly Twisted Quantum Waveguides

We consider the Dirichlet Laplacian \(H_\gamma \) on a 3D twisted waveguide with random Anderson-type twisting \(\gamma \). We introduce the integrated density of states \(N_\gamma \) for the operator \(H_\gamma \), and investigate the Lifshits tails of \(N_\gamma \), i.e. the asymptotic behavior of \(N_\gamma (E)\) as \(E \downarrow \inf \mathrm{supp}\, dN_\gamma \). In particular, we study the dependence of the Lifshits exponent on the decay rate of the single-site twisting at infinity.     

Theoretical diagnostic and prediction of physical properties of quaternary InGaAsP compound using artificial neural networks optimized by the Levenberg Maquardt algorithm

The quaternaries \(In_{1 - x} Ga_{x} As_{y} P_{1 - y}\) are the main promising elements for the fabrication of optoelectronic devices. The adjustment of their physical parameters is assumed by the change of the molar fraction \(x\) and \(y\). These parameters can be affected by the variation of temperature and pressure. To make the theoretical diagnosis of these materials, it is fundamental to know the energy gap ‘\(\varvec{E}_{\varvec{g}}\)’ and the lattice parameter ‘\(a\)’, over a wide range of chemical compositions \(0 \le x \le 0.47\) and \(0 \le y \le 1\), at different temperatures and pressures. We show that by using the Artificial Neural Network method optimized by the Levenberg Maquardt algorithm ANN-LM, it is possible to obtain results very close to the experiment. The scatter plot and error calculation show that the ANN-LM model provides more accurate values of the lattice parameter than those calculated by Vegard’s law. On the other hand, the energy gap values \(Eg (x, y, T)\) estimated, using the ANN-LM model, proved to be close to the experimental values that those calculated by the empirical equations. In addition, the ANN-LM method allowed us to estimate with great accuracy the values of the energy gap at different temperatures and pressures \(Eg (P, T)\). Our work provides crucial information on the physical properties of the quaternary without the use of approximations, and without taking into account the hypothesis of a perfect agreement between \(InGaAsP\) and \(InP\) substrate.     

Decay widths of the spin-2 partners of the <Emphasis Type="Italic">X</Emphasis>(3872)

We consider the X(3872) resonance as a \(J^\mathrm{{PC}}=1^{++}\) \(D\bar{D}^*\) hadronic molecule. According to heavy quark spin symmetry, there will exist a partner with quantum numbers \(2^{++}\), \(X_{2}\), which would be a \(D^*\bar{D}^*\) loosely bound state. The \(X_{2}\) is expected to decay dominantly into \(D\bar{D}\), \(D\bar{D}^*\) and \(\bar{D} D^*\) in d-wave. In this work, we calculate the decay widths of the \(X_{2}\) resonance into the above channels, as well as those of its bottom partner, \(X_{b2}\), the mass of which comes from assuming heavy flavor symmetry for the contact terms. We find partial widths of the \(X_{2}\) and \(X_{b2}\) of the order of a few MeV. Finally, we also study the radiative \(X_2\rightarrow D\bar{D}^{*}\gamma \) and \(X_{b2} \rightarrow \bar{B} B^{*}\gamma \) decays. These decay modes are more sensitive to the long-distance structure of the resonances and to the \(D\bar{D}^{*}\) or \(B\bar{B}^{*}\) final state interaction.     

Heuristic Relative Entropy Principles with Complex Measures: Large-Degree Asymptotics of a Family of Multi-variate Normal Random Polynomials

Let \(z\in \mathbb {C}\), let \(\sigma ^2>0\) be a variance, and for \(N\in \mathbb {N}\) define the integrals

$$\begin{aligned} E_N^{}(z;\sigma ) := \left\{ \begin{array}{ll} {\frac{1}{\sigma }} \!\!\!\displaystyle \int _{\mathbb {R}}\! (x^2+z^2) \frac{e^{-\frac{1}{2\sigma ^2} x^2}}{\sqrt{2\pi }}dx&{}\quad \text{ if }\, N=1,\\ {\frac{1}{\sigma }} \!\!\!\displaystyle \int _{\mathbb {R}^N}\! \prod \prod \limits _{1\le k<l\le N}\!\! e^{-\frac{1}{2N}(1-\sigma ^{-2}) (x_k-x_l)^2} \prod _{1\le n\le N}\!\!\!\!(x_n^2+z^2) \frac{e^{-\frac{1}{2\sigma ^2} x_n^2}}{\sqrt{2\pi }}dx_n &{}\quad \text{ if }\, N>1. \end{array}\right. \!\!\! \end{aligned}$$

These are expected values of the polynomials \(P_N^{}(z)=\prod _{1\le n\le N}(X_n^2+z^2)\) whose 2N zeros \(\{\pm i X_k\}^{}_{k=1,\ldots ,N}\) are generated by N identically distributed multi-variate mean-zero normal random variables \(\{X_k\}^{N}_{k=1}\) with co-variance \(\mathrm{{Cov}}_N^{}(X_k,X_l)=(1+\frac{\sigma ^2-1}{N})\delta _{k,l}+\frac{\sigma ^2-1}{N}(1-\delta _{k,l})\). The \(E_N^{}(z;\sigma )\) are polynomials in \(z^2\), explicitly computable for arbitrary N, yet a list of the first three \(E_N^{}(z;\sigma )\) shows that the expressions become unwieldy already for moderate N—unless \(\sigma = 1\), in which case \(E_N^{}(z;1) = (1+z^2)^N\) for all \(z\in \mathbb {C}\) and \(N\in \mathbb {N}\). (Incidentally, commonly available computer algebra evaluates the integrals \(E_N^{}(z;\sigma )\) only for N up to a dozen, due to memory constraints). Asymptotic evaluations are needed for the large-N regime. For general complex z these have traditionally been limited to analytic expansion techniques; several rigorous results are proved for complex z near 0. Yet if \(z\in \mathbb {R}\) one can also compute this “infinite-degree” limit with the help of the familiar relative entropy principle for probability measures; a rigorous proof of this fact is supplied. Computer algebra-generated evidence is presented in support of a conjecture that a generalization of the relative entropy principle to signed or complex measures governs the \(N\rightarrow \infty \) asymptotics of the regime \(iz\in \mathbb {R}\). Potential generalizations, in particular to point vortex ensembles and the prescribed Gauss curvature problem, and to random matrix ensembles, are emphasized.

     

Scaling exponents of forced polymer translocation through a nanopore

We investigate several properties of a translocating homopolymer through a thin pore driven by an external field present inside the pore only using Langevin Dynamics (LD) simulations in three dimensions (3D). Motivated by several recent theoretical and numerical studies that are apparently at odds with each other, we estimate the exponents describing the scaling with chain length (Nof the average translocation time \(\ensuremath \langle\tau\rangle\) , the average velocity of the center of mass \(\ensuremath \langle v_{{\rm CM}}\rangle\) , and the effective radius of gyration \(\ensuremath \langle {R}_g\rangle\) during the translocation process defined as \(\ensuremath \langle\tau\rangle \sim N^{\alpha}\) , \(\ensuremath \langle v_{{\rm CM}} \rangle \sim N^{-\delta}\) , and \(\ensuremath {R}_g \sim N^{\bar{\nu}}\) respectively, and the exponent of the translocation coordinate (s -coordinate) as a function of the translocation time \(\ensuremath \langle s^2(t)\rangle\sim t^{\beta}\) . We find \(\ensuremath \alpha=1.36 \pm 0.01\) , \(\ensuremath \beta=1.60 \pm 0.01\) for \(\ensuremath \langle s^2(t)\rangle\sim \tau^{\beta}\) and \(\ensuremath \bar{\beta}=1.44 \pm 0.02\) for \(\ensuremath \langle\Delta s^2(t)\rangle\sim\tau^{\bar{\beta}}\) , \(\ensuremath \delta=0.81 \pm 0.04\) , and \(\ensuremath \bar{\nu}\simeq\nu=0.59 \pm 0.01\) , where \( \nu\) is the equilibrium Flory exponent in 3D. Therefore, we find that \(\ensuremath \langle\tau\rangle\sim N^{1.36}\) is consistent with the estimate of \(\ensuremath \langle\tau\rangle\sim\langle R_g \rangle/\langle v_{{\rm CM}} \rangle\) . However, as observed previously in Monte Carlo (MC) calculations by Kantor and Kardar (Y. Kantor, M. Kardar, Phys. Rev. E 69, 021806 (2004)) we also find the exponent α = 1.36 ± 0.01 < 1 + ν. Further, we find that the parallel and perpendicular components of the gyration radii, where one considers the “cis” and “trans” parts of the chain separately, exhibit distinct out-of-equilibrium effects. We also discuss the dependence of the effective exponents on the pore geometry for the range of N studied here.     

Systematic study of rigid triaxiality in Ba–Pt nuclei and role of $$$$ subshell effect

A systematic variation of shape variables \(\beta \) and \(\gamma \) with N and \(N_\mathrm {p}N_\mathrm {n}\) is studied in the framework of an asymmetric rotor model of Davydov and Filippov for the \(Z=\) 50–82, \(N=\) 82–126 major shell space. The role of the \(Z=64\) subshell in producing smooth systematics has been discussed. The quadrupole moments are extracted after considering both axially symmetric and axially asymmetric nuclei. The correlation of \(\beta \) with \(\gamma \) together with the measured quadrupole moments indicates that \(\gamma \)-rigidity is better observed in nuclei with modest deformation.     

On quantum analogue of dynamical stabilisation of inverted harmonic oscillator by time periodical uniform field

Quantum analogue of stabilised forced oscillations around an unstable equilibrium position is explored by solving the non-stationary Schrödinger equation (NSE) of the inverted harmonic oscillator (IHO) driven periodically by spatial uniform field of frequency \(\Omega \), amplitude \(F_{0}\) and phase \(\phi \), i.e. the system with the Hamiltonian of \(\hat{{H}}=(\hat{{p}}^{2}/2m)-(m\omega ^{2}x^{2}/2)-F_0 x\sin \) \(\left( {\Omega t+\phi } \right) \). The NSE has been solved both analytically and numerically by Maple 15 in dimensionless variables \(\xi = x\sqrt{m\omega /\hbar }\hbox {, }f_0 =F_0 /\omega \sqrt{\hbar m\omega }\) and \(\tau =\omega t\). The initial condition (IC) has been specified by the wave function (w.f.) of a generalised Gaussian type which suits well the corresponding quantum IC operator. The solution obtained demonstrates the non-monotonous behaviour of the coordinate spreading \(\sigma \left( \tau \right) \hbox { =}\sqrt{\big ( {\overline{\Delta \xi ^{2}\big ( \tau \big )} } \big )}\) which decreases first from quite macroscopic values of \(\sigma _{0} =2^{12,\ldots ,25}\) to minimal one of \(\sim \!(1/\sqrt{2})\) at times \(\tau <\tau _0 =0.125\ln \!\left( {16\sigma _0^4 +1} \right) \) and then grows back unlimitedly. For certain phases \(\phi \) depending on the \(\Omega /\omega \) ratio and \(n=\log _2\!\sigma _0 \), the mass centre of the packet \(\xi _{\mathrm {av}}( \tau )= \overline{\hat{{x}}(\tau )} \cdot \sqrt{m\omega /\hbar }\) delays approximately two natural ‘periods’ \(\sim \!(4\pi /\omega )\) in the area of the stationary point and then escapes to ‘\(+\)’ or ‘?’ infinity in a bifurcating way.  For ‘resonant’ \(\Omega =\omega \), the bifurcation phases \(\phi \) fit well with the regression formula of Fermi–Dirac type of argument n with their asymptotic \(\phi ( {\Omega ,n\rightarrow \infty } )\) obeying the classical formula \(\phi _{\mathrm {cl}} ( \Omega )=-\hbox {arctg} \, \Omega \) for initial energy \(E = 0\) in the wide range of \(\Omega =2^{-4},...,2^{7}\).     

Likelihood analysis of the pMSSM11 in light of LHC 13-TeV data

We use MasterCode to perform a frequentist analysis of the constraints on a phenomenological MSSM model with 11 parameters, the pMSSM11, including constraints from \(\sim 36\)/fb of LHC data at 13 TeV and PICO, XENON1T and PandaX-II searches for dark matter scattering, as well as previous accelerator and astrophysical measurements, presenting fits both with and without the \((g-2)_\mu \) constraint. The pMSSM11 is specified by the following parameters: 3 gaugino masses \(M_{1,2,3}\), a common mass for the first-and second-generation squarks \(m_{\tilde{q}}\) and a distinct third-generation squark mass \(m_{\tilde{q}_3}\), a common mass for the first-and second-generation sleptons \(m_{\tilde{\ell }}\) and a distinct third-generation slepton mass \(m_{\tilde{\tau }}\), a common trilinear mixing parameter A, the Higgs mixing parameter \(\mu \), the pseudoscalar Higgs mass \(M_A\) and \(\tan \beta \). In the fit including \((g-2)_\mu \), a Bino-like \(\tilde{\chi }^0_{1}\) is preferred, whereas a Higgsino-like \(\tilde{\chi }^0_{1}\) is mildly favoured when the \((g-2)_\mu \) constraint is dropped. We identify the mechanisms that operate in different regions of the pMSSM11 parameter space to bring the relic density of the lightest neutralino, \(\tilde{\chi }^0_{1}\), into the range indicated by cosmological data. In the fit including \((g-2)_\mu \), coannihilations with \(\tilde{\chi }^0_{2}\) and the Wino-like \(\tilde{\chi }^\pm _{1}\) or with nearly-degenerate first- and second-generation sleptons are active, whereas coannihilations with the \(\tilde{\chi }^0_{2}\) and the Higgsino-like \(\tilde{\chi }^\pm _{1}\) or with first- and second-generation squarks may be important when the \((g-2)_\mu \) constraint is dropped. In the two cases, we present \(\chi ^2\) functions in two-dimensional mass planes as well as their one-dimensional profile projections and best-fit spectra. Prospects remain for discovering strongly-interacting sparticles at the LHC, in both the scenarios with and without the \((g-2)_\mu \) constraint, as well as for discovering electroweakly-interacting sparticles at a future linear \(e^+ e^-\) collider such as the ILC or CLIC.     

Isospin analysis of charmless <Emphasis Type="Italic">B</Emphasis>-meson decays

We discuss the determination of the CKM angle \(\alpha \) using the non-leptonic two-body decays \(B\rightarrow \pi \pi \), \(B\rightarrow \rho \rho \) and \(B\rightarrow \rho \pi \) using the latest data available. We illustrate the methods used in each case and extract the corresponding value of \(\alpha \). Combining all these elements, we obtain the determination \(\alpha _\mathrm{dir}={({86.2}_{-4.0}^{+4.4} \cup {178.4}_{-5.1}^{+3.9})}^{\circ }\). We assess the uncertainties associated to the breakdown of the isospin hypothesis and the choice of the statistical framework in detail. We also determine the hadronic amplitudes (tree and penguin) describing the QCD dynamics involved in these decays, briefly comparing our results with theoretical expectations. For each observable of interest in the \(B\rightarrow \pi \pi \), \(B\rightarrow \rho \rho \) and \(B\rightarrow \rho \pi \) systems, we perform an indirect determination based on the constraints from all the other observables available and we discuss the compatibility between indirect and direct determinations. Finally, we review the impact of future improved measurements on the determination of \(\alpha \).     

Scaling Limit of Symmetric Random Walk in High-Contrast Periodic Environment

It is shown that the deterministic infinite trigonometric products

$$\begin{aligned} \prod _{n\in \mathbb {N}}\left[ 1- p +p\cos \left( \textstyle n^{-s}_{_{}}t\right) \right] =: {\text{ Cl }_{p;s}^{}}(t) \end{aligned}$$

with parameters \( p\in (0,1]\ \& \ s>\frac{1}{2}\), and variable \(t\in \mathbb {R}\), are inverse Fourier transforms of the probability distributions for certain random series \(\Omega _{p}^\zeta (s)\) taking values in the real \(\omega \) line; i.e. the \({\text{ Cl }_{p;s}^{}}(t)\) are characteristic functions of the \(\Omega _{p}^\zeta (s)\). The special case \(p=1=s\) yields the familiar random harmonic series, while in general \(\Omega _{p}^\zeta (s)\) is a “random Riemann-\(\zeta \) function,” a notion which will be explained and illustrated—and connected to the Riemann hypothesis. It will be shown that \(\Omega _{p}^\zeta (s)\) is a very regular random variable, having a probability density function (PDF) on the \(\omega \) line which is a Schwartz function. More precisely, an elementary proof is given that there exists some \(K_{p;s}^{}>0\), and a function \(F_{p;s}^{}(|t|)\) bounded by \(|F_{p;s}^{}(|t|)|\!\le \! \exp \big (K_{p;s}^{} |t|^{1/(s+1)})\), and \(C_{p;s}^{}\!:=\!-\frac{1}{s}\int _0^\infty \ln |{1-p+p\cos \xi }|\frac{1}{\xi ^{1+1/s}}\mathrm{{d}}\xi \), such that

$$\begin{aligned} \forall \,t\in \mathbb {R}:\quad {\text{ Cl }_{p;s}^{}}(t) = \exp \bigl ({- C_{p;s}^{} \,|t|^{1/s}\bigr )F_{p;s}^{}(|t|)}; \end{aligned}$$

the regularity of \(\Omega _{p}^\zeta (s)\) follows. Incidentally, this theorem confirms a surmise by Benoit Cloitre, that \(\ln {\text{ Cl }_{{{1}/{3}};2}^{}}(t) \sim -C\sqrt{t}\; \left( t\rightarrow \infty \right) \) for some \(C>0\). Graphical evidence suggests that \({\text{ Cl }_{{{1}/{3}};2}^{}}(t)\) is an empirically unpredictable (chaotic) function of t. This is reflected in the rich structure of the pertinent PDF (the Fourier transform of \({\text{ Cl }_{{{1}/{3}};2}^{}}\)), and illustrated by random sampling of the Riemann-\(\zeta \) walks, whose branching rules allow the build-up of fractal-like structures.

     

Detection of relic gravitational waves in thermal case by using Adv.LIGO data of GW150914

The thermal spectrum of relic gravitational waves enhances the usual spectrum. Our analysis shows that there exist some chances for detection of the thermal spectrum in addition to the usual spectrum by comparison with sensitivity of Adv.LIGO of GW150914 and detector based on the maser light. The behavior of the inflation and reheating stages are often known as power law expansion like \(S(\eta )\propto \eta ^{1+\beta }\), \(S(\eta )\propto \eta ^{1+\beta _s}\), respectively, with constraints \(1+\beta <0, 1+\beta _{s}>0\). The \(\beta \) and \(\beta _s\) have an unique effect on the shape of the spectrum. We find some values of the \(\beta \) and \(\beta _s\) by considering the mentioned comparison. As obtained, the results give us more information as regards the evolution of inflation and reheating stages.