The Colored Hofstadter Butterfly for the Honeycomb Lattice

We rely on a recent method for determining edge spectra and we use it to compute the Chern numbers for Hofstadter models on the honeycomb lattice having rational magnetic flux per unit cell. Based on the bulk-edge correspondence, the Chern number \(\sigma _\mathrm{H}\) is given as the winding number of an eigenvector of a \(2 \times 2\) transfer matrix, as a function of the quasi-momentum \(k\in (0,2\pi )\) . This method is computationally efficient (of order \(\mathcal {O}(n^4)\) in the resolution of the desired image). It also shows that for the honeycomb lattice the solution for \(\sigma _\mathrm{H}\) for flux \(p/q\) in the \(r\) -th gap conforms with the Diophantine equation \(r=\sigma _\mathrm{H}\cdot p+ s\cdot q\) , which determines \(\sigma _\mathrm{H}\mod q\) . A window such as \(\sigma _\mathrm{H}\in (-q/2,q/2)\) , or possibly shifted, provides a natural further condition for \(\sigma _\mathrm{H}\) , which however turns out not to be met. Based on extensive numerical calculations, we conjecture that the solution conforms with the relaxed condition \(\sigma _\mathrm{H}\in (-q,q)\) .     

Optimal N-Point Configurations on the Sphere: “Magic” Numbers and Smale’s 7th Problem

This paper inquires into the concavity of the map \(N\mapsto v_s(N)\) from the integers \(N\ge 2\) into the minimal average standardized Riesz pair-energies \(v_s(N)\) of \(N\) -point configurations on the sphere \(\mathbb {S}^2\) for various \(s\in \mathbb {R}\) . The standardized Riesz pair-energy of a pair of points on \(\mathbb {S}^2\) a chordal distance \(r\) apart is \(V_s(r)= s^{-1}\left( r^{-s}-1 \right) \) , \(s \ne 0\) , which becomes \(V_0(r) = \ln \frac{1}{r}\) in the limit \(s\rightarrow 0\) . Averaging it over the \(\left( \begin{array}{c} N\\ 2\end{array}\right) \) distinct pairs in a configuration and minimizing over all possible \(N\) -point configurations defines \(v_s(N)\) . It is known that \(N\mapsto v_s(N)\) is strictly increasing for each \(s\in \mathbb {R}\) , and for \(s<2\) also bounded above, thus “overall concave.” It is (easily) proved that \(N\mapsto v_{-2}^{}(N)\) is even locally strictly concave, and that so is the map \(2n\mapsto v_s(2n)\) for \(s<-2\) . By analyzing computer-experimental data of putatively minimal average Riesz pair-energies \(v_s^x(N)\) for \(s\in \{-1,0,1,2,3\}\) and \(N\in \{2,\ldots ,200\}\) , it is found that the map \(N\mapsto {v}_{-1}^x(N)\) is locally strictly concave, while \(N\mapsto {v}_s^x(N)\) is not always locally strictly concave for \(s\in \{0,1,2,3\}\) : concavity defects occur whenever \(N\in {\mathcal {C}}^{x}_+(s)\) (an \(s\) -specific empirical set of integers). It is found that the empirical map \(s\mapsto {\mathcal {C}}^{x}_+(s),\ s\in \{-2,-1,0,1,2,3\}\) , is set-theoretically increasing; moreover, the percentage of odd numbers in \({\mathcal {C}}^{x}_+(s),\ s\in \{0,1,2,3\}\) is found to increase with \(s\) . The integers in \({\mathcal {C}}^{x}_+(0)\) are few and far between, forming a curious sequence of numbers, reminiscent of the “magic numbers” in nuclear physics. It is conjectured that these new “magic numbers” are associated with optimally symmetric optimal-log-energy \(N\) -point configurations on \(\mathbb {S}^2\) . A list of interesting open problems is extracted from the empirical findings, and some rigorous first steps toward their solutions are presented. It is emphasized how concavity can assist in the solution to Smale’s \(7\) th Problem, which asks for an efficient algorithm to find near-optimal \(N\) -point configurations on \(\mathbb {S}^2\) and higher-dimensional spheres.     

Dispersion dependence of two-photon absorption transition on frequency in Si PIN photodetector

The variation of two-photon absorption (TPA) coefficient \(\beta _{\mathrm{TPA}} (\omega )\) of Si excited at difference photon energy was investigated. The TPA coefficient was measured by using a picosecond pulsed laser with the wavelength could be tuned in a wide photon-energy range. An equivalent RC circuit model was adapted to derive the TPA coefficient \(\beta _{\mathrm{TPA}} (\omega )\) . The results showed that \(\beta _{\mathrm{TPA}} (\omega )\) varied from \(4.2 \times 10^{-4}\) to \(1.17 \times 10^{-3 }\)  cm/GW in the transparent wavelength region \(1.80<\lambda <1.36\,\upmu \) m of Si. The increasing tendency of \(\beta _{\mathrm{TPA}} (\omega )\) with the incident photon energy can be qualitatively interpreted as the photon energy increases from \(E_{\mathrm{ig}}/2\) to nearly \(E_{\mathrm{ig}}\) , the electrons excited from the valance band find an increasing availability of conduction band states. Comparing with the high-energy side transitions, the TPA coefficient in low-energy side is about 10 times too small. This can be attributed that the TPA transition in low-energy side is the process of photon-assisted electron transitions from valence to conduction band occurring between different points in k-space, while is direct transition in high-energy side.     

Pesin’s Entropy Formula for Systems Between C^1 and C^{1+\alpha }

In this article we give a new observation of Pesin’s entropy formula, motivated from Mañé’s proof of (Ergod Theory Dyn Sys 1:95–102, 1981). Let \(M\) be a compact Riemann manifold and \(f:\,M\rightarrow M\) be a \(C^1\) diffeomorphism on \(M\) . If \(\mu \) is an \(f\) -invariant probability measure which is absolutely continuous relative to Lebesgue measure and nonuniformly-H \(\ddot{\text {o}}\) lder-continuous(see Definition 1.1), then we have Pesin’s entropy formula, i.e., the metric entropy \(h_\mu (f)\) satisfies $$\begin{aligned} h_{\mu }(f)=\int \sum _{\lambda _i(x)> 0}\lambda _i(x)d\mu , \end{aligned}$$ where \(\lambda _1(x)\ge \lambda _2(x)\ge \cdots \ge \lambda _{dim\,M}(x)\) are the Lyapunov exponents at \(x\) with respect to \(\mu .\) Nonuniformly-H \(\ddot{\text {o}}\) lder-continuous is a new notion from probabilistic perspective weaker than \(C^{1+\alpha }.\)      

Split supersymmetry under GUT and current dark matter constraints

There are four types of two-Higgs doublet models under a discrete \(Z_2\) symmetry imposed to avoid tree-level flavor-changing neutral current, i.e. type-I, type-II, type-X, and type-Y models. We investigate the possibility to discriminate the four models in the light of the flavor physics data, including \(B_s\) – \(\bar{B}_s\) mixing, \(B_{s,d} \rightarrow \mu ^+ \mu ^-\) , \(B\rightarrow \tau \nu \) and \(\bar{B} \rightarrow X_s \gamma \) decays, the recent LHC Higgs data, the direct search for charged Higgs at LEP, and the constraints from perturbative unitarity and vacuum stability. After deriving the combined constraints on the Yukawa interaction parameters, we have shown that the correlation between the mass eigenstate rate asymmetry \(A_{\Delta \Gamma }\) of \(B_{s} \rightarrow \mu ^+ \mu ^-\) and the ratio \(R=\mathcal{B}(B_{s} \rightarrow \mu ^+ \mu ^-)_\mathrm{exp}/ \mathcal{B}(B_{s} \rightarrow \mu ^+ \mu ^-)_\mathrm{SM}\) could be a sensitive probe to discriminate the four models with future precise measurements of the observables in the \(B_{s} \rightarrow \mu ^+ \mu ^-\) decay at LHCb.     

Is the CNMSSM more credible than the CMSSM?

“Post-sphaleron baryogenesis”, a fresh and profound mechanism of baryogenesis accounts for the matter–antimatter asymmetry of our present universe in a framework of Pati–Salam symmetry. We attempt here to embed this mechanism in a non-SUSY SO(10) grand unified theory by reviving a novel symmetry breaking chain with Pati–Salam symmetry as an intermediate symmetry breaking step and as well to address post-sphaleron baryogenesis and neutron–antineutron oscillation in a rational manner. The Pati–Salam symmetry based on the gauge group \(\mathrm{SU}(2)_L \times \mathrm{SU}(2)_{R} \times \mathrm{SU}(4)_C\) is realized in our model at \(10^{5}\) – \(10^{6}\)  GeV and the mixing time for the neutron–antineutron oscillation process having \(\Delta B=2\) is found to be \(\tau _{n-\bar{n}} \simeq 10^{8}\) – \(10^{10}\)  s with the model parameters, which is within the reach of forthcoming experiments. Other novel features of the model include low scale right-handed \(W^{\pm }_R\) , \(Z_R\) gauge bosons, explanation for neutrino oscillation data via the gauged inverse (or extended) seesaw mechanism and most importantly TeV scale color sextet scalar particles responsible for an observable \(n\) – \(\bar{n}\) oscillation which may be accessible to LHC. We also look after gauge coupling unification and an estimation of the proton lifetime with and without the addition of color sextet scalars.     

Toward a Mathematical Holographic Principle

In work started in [17] and continued in this paper our objective is to study selectors of multivalued functions which have interesting dynamical properties, such as possessing absolutely continuous invariant measures. We specify the graph of a multivalued function by means of lower and upper boundary maps \(\tau _{1}\) and \(\tau _{2}.\) On these boundary maps we define a position dependent random map \(R_{p}=\{\tau _{1},\tau _{2};p,1-p\},\) which, at each time step, moves the point \(x\) to \(\tau _{1}(x)\) with probability \(p(x)\) and to \(\tau _{2}(x)\) with probability \(1-p(x)\) . Under general conditions, for each choice of \(p\) , \(R_{p}\) possesses an absolutely continuous invariant measure with invariant density \(f_{p}.\) Let \(\varvec{\tau }\) be a selector which has invariant density function \(f.\) One of our objectives is to study conditions under which \(p(x)\) exists such that \(R_{p}\) has \(f\) as its invariant density function. When this is the case, the long term statistical dynamical behavior of a selector can be represented by the long term statistical behavior of a random map on the boundaries of \(G.\) We refer to such a result as a mathematical holographic principle. We present examples and study the relationship between the invariant densities attainable by classes of selectors and the random maps based on the boundaries and show that, under certain conditions, the extreme points of the invariant densities for selectors are achieved by bang-bang random maps, that is, random maps for which \(p(x)\in \{0,1\}.\)      

Comprehensive Bayesian analysis of rare (semi)leptonic and radiative \varvec{B} decays

The available data on \(|\Delta B| = |\Delta S| = 1\) decays are in good agreement with the Standard Model when permitting subleading power corrections of about \(15\,\%\) at large hadronic recoil. Constraining new-physics effects in \(\mathcal {C}_{7}^{\mathrm {}}\) , \(\mathcal {C}_{9}^{\mathrm {}}\) , \(\mathcal {C}_{10}^{\mathrm {}}\) , the data still demand the same size of power corrections as in the Standard Model. In the presence of chirality-flipped operators, all but one of the power corrections reduce substantially. The Bayes factors are in favor of the Standard Model. Using new lattice inputs for \(B\rightarrow K^*\) form factors and under our minimal prior assumption for the power corrections, the favor shifts toward models with chirality-flipped operators. We use the data to further constrain the hadronic form factors in \(B\rightarrow K\) and \(B\rightarrow K^*\) transitions.     

QCD analysis of forward neutron production in DIS

Observing light-by-light scattering at the large hadron collider (LHC) has received quite some attention and it is believed to be a clean and sensitive channel to possible new physics. In this paper, we study the diphoton production at the LHC via the process \({{pp}}\rightarrow {{p}}\gamma \gamma {{p}}\rightarrow {{p}}\gamma \gamma {{p}}\) through graviton exchange in the large extra dimension (LED) model. Typically, when we do the background analysis, we also study the double Pomeron exchange of \(\gamma \gamma \) production. We compare its production in the quark–quark collision mode to the gluon–gluon collision mode and find that contributions from the gluon–gluon collision mode are comparable to the quark–quark one. Our result shows, for extra dimension \(\delta =4\) , with an integrated luminosity \(\mathcal{L} = 200\,\mathrm{fb}^{-1}\) at the 14 TeV LHC, that diphoton production through graviton exchange can probe the LED effects up to the scale \({M}_{S}=5.06 (4.51, 5.11)\,\mathrm{TeV}\) for the forward detector acceptance \(\xi _1 (\xi _2, \xi _3)\) , respectively, where \(0.0015<\xi _1<0.5\) , \(0.1<\xi _2<0.5\) , and \(0.0015<\xi _3<0.15\) .     

Deformed spinor networks for loop gravity: towards hyperbolic twisted geometries

In the context of a canonical quantization of general relativity, one can deform the loop gravity phase space on a graph by replacing the \(T^*\mathrm{SU}(2)\) phase space attached to each edge by \(\mathrm{SL}(2,{\mathbb C})\) seen as a phase space. This deformation is supposed to encode the presence of a non-zero cosmological constant.Here we show how to parametrize this phase space in terms of spinor variables, thus obtaining deformed spinor networks for loop gravity, with a deformed action of the gauge group \(\mathrm{SU}(2)\) at the vertices. These are to be formally interpreted as the generalization of loop gravity twisted geometries to a hyperbolic curvature.     

The CMSSM and NUHM1 after LHC Run 1

We analyze the impact of data from the full Run 1 of the LHC at 7 and 8 TeV on the CMSSM with \(\mu > 0\) and \(<0\) and the NUHM1 with \(\mu > 0\) , incorporating the constraints imposed by other experiments such as precision electroweak measurements, flavour measurements, the cosmological density of cold dark matter and the direct search for the scattering of dark matter particles in the LUX experiment. We use the following results from the LHC experiments: ATLAS searches for events with \({E\!\!/}_{T}\) accompanied by jets with the full 7 and 8 TeV data, the ATLAS and CMS measurements of the mass of the Higgs boson, the CMS searches for heavy neutral Higgs bosons and a combination of the LHCb and CMS measurements of \(\mathrm{BR}(B_s \rightarrow \mu ^+\mu ^-)\) and \(\mathrm{BR}(B_d \rightarrow \mu ^+\mu ^-)\) . Our results are based on samplings of the parameter spaces of the CMSSM for both \(\mu >0\) and \(\mu <0\) and of the NUHM1 for \(\mu > 0\) with 6.8 \(\times 10^6\) , 6.2 \(\times 10^6\) and 1.6 \(\times 10^7\) points, respectively, obtained using the MultiNest tool. The impact of the Higgs-mass constraint is assessed using FeynHiggs 2.10.0, which provides an improved prediction for the masses of the MSSM Higgs bosons in the region of heavy squark masses. It yields in general larger values of \(M_h\) than previous versions of FeynHiggs, reducing the pressure on the CMSSM and NUHM1. We find that the global \(\chi ^2\) functions for the supersymmetric models vary slowly over most of the parameter spaces allowed by the Higgs-mass and the \({E\!\!/}_{T}\) searches, with best-fit values that are comparable to the \(\chi ^2/\mathrm{dof}\) for the best Standard Model fit. We provide 95 % CL lower limits on the masses of various sparticles and assess the prospects for observing them during Run 2 of the LHC.     

Standard Model Gauge Couplings from Gauge-Dilatation Symmetry Breaking

It is well known that the self-energy of the gauge bosons is quadratically divergent in the Standard Model when a simple cutoff is imposed. We demonstrate phenomenologically that the quadratic divergences in fact unify. The unification occurs at a surprisingly low scale, \(\Lambda _\mathrm {u}\approx 4\times 10^7\)  GeV. Suppose now that there is a spontaneously broken rotational symmetry between the space-time coordinates and gauge theoretical phases. The symmetry-breaking pattern is such that the gauge bosons arise as the massless Goldstone bosons, whereas the dilatonic mode acts as the massive (Higgs) boson, whose vacuum expectation value determines the gauge couplings. In this case, the quadratic divergences or the tadpoles of the gauge boson self-energy should indeed unify because these divergences need to be cancelled by a universal dilatonic contribution, assuming dynamical symmetry breaking. If there is dynamical symmetry breaking, we are in principle able to calculate the value of the gauge couplings as well as the scale hierarchy \(\Lambda _\mathrm {cut}/\Lambda _\mathrm {u}\) . We perform this calculation by adopting a naive quartic symmetry-breaking potential which unfortunately violates local gauge invariance. Using tadpole-cancellation and dilatonic self-energy conditions, the value of \(\Lambda _\mathrm {cut}\) is then found to be approximately \(4\times 10^{18}\)  GeV in the Feynman gauge and \(5\times 10^{15}\)  GeV in the Landau gauge. The cancellation of an anomaly in the dilaton self-energy requires that the number of fermionic generations equals three. The symmetry-breaking needs to be driven by some other mass-generating mechanism such as electroweak symmetry breaking. Our estimation for \(\Lambda _\mathrm {u}\) is of the correct order if \(\Lambda _\mathrm {cut}\approx 5\times 10^{15}\)  GeV.     

Phase Transitions in Layered Systems

We consider the Ising model on \(\mathbb Z\times \mathbb Z\) where on each horizontal line \(\{(x,i), x\in \mathbb Z\}\) , called “layer”, the interaction is given by a ferromagnetic Kac potential with coupling strength \(J_{ \gamma }(x,y)={ \gamma }J({ \gamma }(x-y))\) , where \(J(\cdot )\) is smooth and has compact support; we then add a nearest neighbor ferromagnetic vertical interaction of strength \({ \gamma }^{A}\) , where \(A\ge 2\) is fixed, and prove that for any \(\beta \) larger than the mean field critical value there is a phase transition for all \({ \gamma }\) small enough.     

Reunion Probabilities of N One-Dimensional Random Walkers with Mixed Boundary Conditions

In this work we extend the results of the reunion probability of \(N\) one-dimensional random walkers to include mixed boundary conditions between their trajectories. The level of the mixture is controlled by a parameter \(c\) , which can be varied from \(c=0\) (independent walkers) to \(c\rightarrow \infty \) (vicious walkers). The expressions are derived by using Quantum Mechanics formalism (QMf) which allows us to map this problem into a Lieb-Liniger gas (LLg) of \(N\) one-dimensional particles. We use Bethe ansatz and Gaudin’s conjecture to obtain the normalized wave-functions and use this information to construct the propagator. As it is well-known, depending on the boundary conditions imposed at the endpoints of a line segment, the statistics of the maximum heights of the reunited trajectories have some connections with different ensembles in Random Matrix Theory. Here we seek to extend those results and consider four models: absorbing, periodic, reflecting, and mixed. In all four cases, the probability that the maximum height is less or equal than \(L\) takes the form \(F_N(L)=A_N\sum _{\varvec{k}\in \Omega _{\text {B}}} \mathrm{e}^{-\sum _{j=1}^Nk_j^2}\mathcal {V}_N(\varvec{k})\) , where \(A_N\) is a normalization constant, \(\mathcal {V}_N(\varvec{k})\) contains a deformed and weighted Vandermonde determinant, and \(\Omega _{\text {B}}\) is the solution set of quasi-momenta \(\varvec{k}\) obeying the Bethe equations for that particular boundary condition.     

Ga-vacancy-induced room-temperature ferromagnetic and adjusted-band-gap behaviors in GaN nanoparticles

Room-temperature ferromagnetism has been found in Ga-deficient GaN grown using the direct reaction of Ga \(_{2}\) O \(_{3}\) powder with NH \(_{3}\) gas. The observed magnetism in GaN induced by Ga vacancies is investigated both experimentally and theoretically. First-principles calculations reveal that the spontaneous spin polarization is created by the 3.0  \(\mu _\mathrm{B}\) local moment for GaN and magnetism originates from the polarization of the unpaired 2 \(p\) electrons of N surrounding the Ga vacancy. At the same time, the band gap can be also adjusted by changing the Ga-vacancy concentration.     

Dilaton gravity, (quasi-) black holes,and scalar charge

Electrically charged dust is considered in the framework of Einstein–Maxwell–dilaton gravity with a Lagrangian containing the interaction term \(P(\chi )F_{\mu \nu }F^{\mu \nu }\) , where \(P(\chi )\) is an arbitrary function of the dilaton scalar field \(\chi \) , which can be normal or phantom. Without assumption of spatial symmetry, we show that static configurations exist for arbitrary functions \(g_{00} = \exp (2\gamma (x^{i}))\) ( \(i=1,2,3\) ) and \(\chi =\chi (\gamma )\) . If \(\chi = \mathrm{const}\) , the classical Majumdar–Papapetrou (MP) system is restored. We discuss solutions that represent black holes (BHs) and quasi-black holes (QBHs), deduce some general results and confirm them by examples. In particular, we analyze configurations with spherical and cylindrical symmetries. It turns out that cylindrical BHs and QBHs cannot exist without negative energy density somewhere in space. However, in general, BHs and QBHs can be phantom-free, that is, can exist with everywhere nonnegative energy densities of matter, scalar and electromagnetic fields.     

Semiclassical strings in supergravity PFT

In this paper, we introduce the bulk viscosity in the formalism of modified gravity theory in which the gravitational action contains a general function \(f(R,T)\) , where \(R\) and \(T\) denote the curvature scalar and the trace of the energy–momentum tensor, respectively, within the framework of a flat Friedmann–Robertson–Walker model. As an equation of state for a prefect fluid, we take \(p=(\gamma -1)\rho \) , where \(0 \le \gamma \le 2\) and a viscous term as a bulk viscosity due to the isotropic model, of the form \(\zeta =\zeta _{0}+\zeta _{1}H\) , where \(\zeta _{0}\) and \(\zeta _{1}\) are constants, and \(H\) is the Hubble parameter. The exact non-singular solutions to the corresponding field equations are obtained with non-viscous and viscous fluids, respectively, by assuming a simplest particular model of the form of \(f(R,T) = R+2f(T)\) , where \(f(T)=\alpha T\) ( \(\alpha \) is a constant). A big-rip singularity is also observed for \(\gamma <0\) at a finite value of cosmic time under certain constraints. We study all possible scenarios with the possible positive and negative ranges of \(\alpha \) to analyze the expansion history of the universe. It is observed that the universe accelerates or exhibits a transition from a decelerated phase to an accelerated phase under certain constraints of \(\zeta _0\) and \(\zeta _1\) . We compare the viscous models with the non-viscous one through the graph plotted between the scale factor and cosmic time and find that the bulk viscosity plays a major role in the expansion of the universe. A similar graph is plotted for the deceleration parameter with non-viscous and viscous fluids and we find a transition from decelerated to accelerated phase with some form of bulk viscosity.     

Antiproton low-energy collisions with Ps-atoms and true muonium atoms (μ + μ −)

Three-charge-particle collisions with participation of ultra-slow antiprotons ( \(\overline {\rm {p}}\) ) is the subject of this work. Specifically we compute the total cross sections and corresponding thermal rates of the following three-body reactions: \(\overline {\rm p}+(e^+e^-) \rightarrow \overline {\rm {H}} + e^-\) and \(\overline {\rm p}+(\mu ^+\mu ^-) \rightarrow \overline {\rm {H}}_{\mu } + \mu ^-\) , where \(e^-(\mu ^-)\) is an electron (muon) and \(e^+(\mu ^+)\) is a positron (antimuon) respectively, \(\overline {\rm {H}}=(\overline {\rm p}e^+)\) is an antihydrogen atom and \(\overline {\rm {H}}_{\mu }=(\overline {\rm p}\mu ^+)\) is a muonic antihydrogen atom, i.e. a bound state of \(\overline {\rm {p}}\) and μ ⁺. A set of two-coupled few-body Faddeev-Hahn-type (FH-type) equations is numerically solved in the framework of a modified close-coupling expansion approach.     

A flexible scintillation light apparatus for rare event searches

Compelling experimental evidences of neutrino oscillations and their implication that neutrinos are massive particles have given neutrinoless double beta decay ( \(\beta \beta 0\nu \) ) a central role in astroparticle physics. In fact, the discovery of this elusive decay would be a major breakthrough, unveiling that neutrino and antineutrino are the same particle and that the lepton number is not conserved. It would also impact our efforts to establish the absolute neutrino mass scale and, ultimately, understand elementary particle interaction unification. All current experimental programs to search for \(\beta \beta 0\nu \) are facing with the technical and financial challenge of increasing the experimental mass while maintaining incredibly low levels of spurious background. The new concept described in this paper could be the answer which combines all the features of an ideal experiment: energy resolution, low cost mass scalability, isotope choice flexibility and many powerful handles to make the background negligible. The proposed technology is based on the use of arrays of silicon detectors cooled to 120 K to optimize the collection of the scintillation light emitted by ultra-pure crystals. It is shown that with a 54 kg array of natural CaMoO \(_4\) scintillation detectors of this type it is possible to yield a competitive sensitivity on the half-life of the \(\beta \beta 0\nu \) of \(^{100}\) Mo as high as \(\sim \) \(10^{24}\)  years in only 1 year of data taking. The same array made of \(^{40}\) Ca \(^{\mathrm {nat}}\) MoO \(_4\) scintillation detectors (to get rid of the continuous background coming from the two neutrino double beta decay of \(^{48}\) Ca) will instead be capable of achieving the remarkable sensitivity of \(\sim \) \(10^{25}\)  years on the half-life of \(^{100}\) Mo \(\beta \beta 0\nu \) in only 1 year of measurement.