Noncommutative products of Euclidean spaces

We present natural families of coordinate algebras on noncommutative products of Euclidean spaces \({\mathbb {R}}^{N_1} \times _{\mathcal {R}} {\mathbb {R}}^{N_2}\). These coordinate algebras are quadratic ones associated with an \(\mathcal {R}\)-matrix which is involutive and satisfies the Yang–Baxter equations. As a consequence, they enjoy a list of nice properties, being regular of finite global dimension. Notably, we have eight-dimensional noncommutative euclidean spaces \({\mathbb {R}}^{4} \times _{\mathcal {R}} {\mathbb {R}}^{4}\). Among these, particularly well behaved ones have deformation parameter \(\mathbf{u} \in {\mathbb {S}}^2\). Quotients include seven spheres \({\mathbb {S}}^{7}_\mathbf{u}\) as well as noncommutative quaternionic tori \({\mathbb {T}}^{{\mathbb {H}}}_\mathbf{u} = {\mathbb {S}}^3 \times _\mathbf{u} {\mathbb {S}}^3\). There is invariance for an action of \({{\mathrm{SU}}}(2) \times {{\mathrm{SU}}}(2)\) on the torus \({\mathbb {T}}^{{\mathbb {H}}}_\mathbf{u}\) in parallel with the action of \(\mathrm{U}(1) \times \mathrm{U}(1)\) on a ‘complex’ noncommutative torus \({\mathbb {T}}^2_\theta \) which allows one to construct quaternionic toric noncommutative manifolds. Additional classes of solutions are disjoint from the classical case.     

A Heisenberg Double Addition to the Logarithmic Kazhdan–Lusztig Duality

For a Hopf algebra B, we endow the Heisenberg double \({\mathcal{H}(B^*)}\) with the structure of a module algebra over the Drinfeld double \({\mathcal{D}(B)}\). Based on this property, we propose that \({\mathcal{H}(B^*)}\) is to be the counterpart of the algebra of fields on the quantum-group side of the Kazhdan–Lusztig duality between logarithmic conformal field theories and quantum groups. As an example, we work out the case where B is the Taft Hopf algebra related to the \({\overline{\mathcal{U}}_{\mathfrak{q}} s\ell(2)}\) quantum group that is Kazhdan–Lusztig-dual to (p,1) logarithmic conformal models. The corresponding pair \({(\mathcal{D}(B),\mathcal{H}(B^*))}\) is “truncated” to \({(\overline{\mathcal{U}}_{\mathfrak{q}} s\ell2,\overline{\mathcal{H}}_{\mathfrak{q}} s\ell(2))}\), where \({\overline{\mathcal{H}}_{\mathfrak{q}} s\ell(2)}\) is a \({\overline{\mathcal{U}}_{\mathfrak{q}} s\ell(2)}\) module algebra that turns out to have the form \({\overline{\mathcal{H}}_{\mathfrak{q}} s\ell(2)=\mathbb{C}_{\mathfrak{q}}[z,\partial]\otimes\mathbb{C}[\lambda]/(\lambda^{2p}-1)}\), where \({\mathbb{C}_{\mathfrak{q}}[z,\partial]}\) is the \({\overline{\mathcal{U}}_{\mathfrak{q}} s\ell(2)}\)-module algebra with the relations z ^p  = 0, ? ^p  = 0, and \({\partial z = \mathfrak{q}-\mathfrak{q}^{-1} + \mathfrak{q}^{-2} z\partial}\).     

The Structure of the Kac–Wang–Yan Algebra

The Lie algebra \({\mathcal{D}}\) of regular differential operators on the circle has a universal central extension \({\hat{\mathcal{D}}}\). The invariant subalgebra \({\hat{\mathcal{D}}^+}\) under an involution preserving the principal gradation was introduced by Kac, Wang, and Yan. The vacuum \({\hat{\mathcal{D}}^+}\)-module with central charge \({c \in \mathbb{C}}\), and its irreducible quotient \({\mathcal{V}_c}\), possess vertex algebra structures, and \({\mathcal{V}_c}\) has a nontrivial structure if and only if \({c \in \frac{1}{2}\mathbb{Z}}\). We show that for each integer \({n > 0}\), \({\mathcal{V}_{n/2}}\) and \({\mathcal{V}_{-n}}\) are \({\mathcal{W}}\)-algebras of types \({\mathcal{W}(2, 4,\dots,2n)}\) and \({\mathcal{W}(2, 4,\dots, 2n^2 + 4n)}\), respectively. These results are formal consequences of Weyl’s first and second fundamental theorems of invariant theory for the orthogonal group \({{\rm O}(n)}\) and the symplectic group \({{\rm Sp}(2n)}\), respectively. Based on Sergeev’s theorems on the invariant theory of \({{\rm Osp}(1, 2n)}\) we conjecture that \({\mathcal{V}_{-n+1/2}}\) is of type \({\mathcal{W}(2, 4,\dots, 4n^2 + 8n + 2)}\), and we prove this for \({n = 1}\). As an application, we show that invariant subalgebras of \({\beta\gamma}\)-systems and free fermion algebras under arbitrary reductive group actions are strongly finitely generated.     

Holography of massive M2-brane theory: non-linear extension

We investigate the gauge/gravity duality between the \(\mathcal{N} = 6\) mass-deformed ABJM theory with \(\hbox {U}_k(N)\times \hbox {U}_{-k}(N)\) gauge symmetry and the 11-dimensional supergravity on LLM geometries with SO(2,1)\(\times \)SO(4)/\({\mathbb {Z}}_k\) \(\times \)SO(4)/\({\mathbb {Z}}_k\) isometry, in terms of a KK holography, which involves quadratic order field redefinitions. We establish the quadratic order KK mappings for various gauge invariant fields in order to obtain the canonical 4-dimensional gravity equations of motion and to reduce the LLM solutions to an asymptotically AdS\(_4\) gravity solutions. The non-linearity of the KK maps indicates that we can observe the true purpose of the non-linear KK holography of the LLM solutions. We read the vacuum expectation value of conformal dimension two operator from the asymptotically AdS\(_4\) gravity solutions. For the LLM solutions which are represented by square-shaped Young diagrams, we compare the vacuum expectation value obtained from the holographic procedure with the result obtained from the field theory, which is given by \(\langle \mathcal{O}^{(\Delta =2)}\rangle =\sqrt{k}N^{\frac{3}{2}}f_{(\Delta =2)}+\mathcal{O}(N)\), where \(f_{\Delta }\) is independent of N. Based on this result, we examine the gauge/gravity duality in the large N limit and finite k. We also show that the vacuum expectation values of the massive KK graviton modes are vanishing as expected by the supersymmetry.     

On non-consensus motions of dynamical linear multiagent systems

The first-principle density functional theory (DFT) calculations were employed to investigate the electronic structures, magnetic properties and half-metallicity of \(\text {Ti}_{2}\text {IrZ}\) (Z \(=\) B, Al, Ga, and In) Heusler alloys with \(\text {AlCu}_{2}\text {Mn}\)- and \(\text {CuHg}_{2}\text {Ti}\)-type structures within local density approximation and generalised gradient approximation for the exchange correlation potential. It was found that \(\text {CuHg}_{2}\text {Ti}\)-type structure in ferromagnetic state was energetically more favourable than \(\text {AlCu}_{2}\text {Mn}\)-type structure in all compounds except \(\text {Ti}_{2}\text {IrB}\) which was stable in \(\text {AlCu}_{2}\text {Mn}\)-type structure in non-magnetic state. \(\text {Ti}_{2}\text {IrZ}\) (Z \(=\) B, Al, Ga, and In) alloys in \(\text {CuHg}_{2}\text {Ti}\)-type structure were half-metallic ferromagnets at their equilibrium lattice constants. Half-metallic band gaps were respectively equal to 0.87, 0.79, 0.75, and 0.73 eV for \(\text {Ti}_{2}\text {IrB}\), \(\text {Ti}_{2}\text {IrAl}\), \(\text {Ti}_{2}\text {IrGa}\), and \(\text {Ti}_{2}\text {IrIn}\). The origin of half-metallicity was discussed for \(\text {Ti}_{2}\text {IrGa}\) using the energy band structure. The total magnetic moments of \(\text {Ti}_{2}\text {IrZ}\) (Z \(=\) B, Al, Ga, and In) compounds in \(\text {CuHg}_{2}\text {Ti}\)-type structure were obtained as \(2\mu _{\mathrm{B}}\) per formula unit, which were in agreement with Slater–Pauling rule (\(M_{\mathrm{tot}} =Z_{\mathrm{tot}}-\)18). All the four compounds were half-metals in a wide range of lattice constants indicating that they may be suitable and promising materials for future spintronic applications.     

Dynamic Equations for a Free Particle Motion in a Discrete Space-Time Geometry

One considers the discrete space-time geometry \(\mathcal {G}_{\mathrm {d}}\), which is given on the set of points (events), where the geometry of Minkowski is given. This discrete geometry is not a geometry on lattice. Motion of a free particle is considered in \(\mathcal {G}_{\mathrm {d}}\). Free motion in \(\mathcal {G}_{\mathrm {d}}\) can be reduced to a motion in geometry of Minkowski \(\mathcal {G}_{\mathrm {M}}\) in some force field. Primordial free motion in \(\mathcal {G}_{\mathrm {d}}\) appears to be stochastic. In \(\mathcal {G }_{\mathrm {M}}\) it is difficult to describe the force field responsible for stochastic motion of a particle. The nature of this force field appears to be geometrical.     

Further studies on the exclusive productions of $$J/\psi +\chi _{cJ}$$ ($$J=0,1,2$$J=0,1,2) via $$e^+e^-$$e+e- annihilation at the B factories

By including the interference effect between the QCD and the QED diagrams, we carry out a complete analysis on the exclusive productions of \(e^+e^- \rightarrow J/\psi +\chi _{cJ}\) (\(J=0,1,2\)) at the B factories with \(\sqrt{s}=10.6\) GeV at the next-to-leading-order (NLO) level in \(\alpha _s\), within the nonrelativistic QCD framework. It is found that the \({\mathcal {O}} (\alpha ^3\alpha _s)\)-order terms that represent the tree-level interference are comparable with the usual NLO QCD corrections, especially for the \(\chi _{c1}\) and \(\chi _{c2}\) cases. To explore the effect of the higher-order terms, namely \({\mathcal {O}} (\alpha ^3\alpha _s^2)\), we perform the QCD corrections to these \({\mathcal {O}} (\alpha ^3\alpha _s)\)-order terms for the first time, which are found to be able to significantly influence the \({\mathcal {O}} (\alpha ^3\alpha _s)\)-order results. In particular, in the case of \(\chi _{c1}\) and \(\chi _{c2}\), the newly calculated \({\mathcal {O}} (\alpha ^3\alpha _s^2)\)-order terms can to a large extent counteract the \({\mathcal {O}} (\alpha ^3\alpha _s)\) contributions, evidently indicating the indispensability of the corrections. In addition, we find that, as the collision energy rises, the percentage of the interference effect in the total cross section will increase rapidly, especially for the \(\chi _{c1}\) case.     

Antihydrogen formation via antiproton scattering on positronium between the $e^{-} +\bar {H}(n = 2)$ and e−+H̄(n=3)$e^{-}+\bar {H}(n = 3)$ thresholds

The charge exchange reaction \(\bar {\mathrm {p}} + \text {Ps} \rightarrow \mathrm {e}^{-} + \bar {\mathrm {H}} \), of interest for the future experiments (GBAR, AEGIS, ATRAP, ...) aiming to produce antihydrogen atoms, is investigated in the energy range between the \(\mathrm {e}^{-}+\bar {\mathrm {H}}(n = 2)\) and \(\mathrm {e}^{-}+\bar {\mathrm {H}}(n = 3)\) thresholds. An ab-initio method based on the solution of the Faddeev-Merkuriev equations is used. Special focus is put on the impact of the Feshbach resonances and the Gailitis-Damburg oscillations, appearing in the vicinity of the \(\bar {\mathrm {p}} +\text {Ps}(n = 2)\) threshold, on the \(\bar {\mathrm {H}}\) production cross section.     

Regularity of the Speed of Biased Random Walk in a One-Dimensional Percolation Model

We consider biased random walks on the infinite cluster of a conditional bond percolation model on the infinite ladder graph. Axelson-Fisk and Häggström established for this model a phase transition for the asymptotic linear speed \(\overline{\hbox {v}}\) of the walk. Namely, there exists some critical value \(\lambda _{\hbox {c}}>0\) such that \(\overline{\hbox {v}}>0\) if \(\lambda \in (0,\lambda _{\hbox {c}})\) and \(\overline{\hbox {v}}=0\) if \(\lambda \ge \lambda _{\hbox {c}}\). We show that the speed \(\overline{\hbox {v}}\) is continuous in \(\lambda \) on \((0,\infty )\) and differentiable on \((0,\lambda _{\hbox {c}}/2)\). Moreover, we characterize the derivative as a covariance. For the proof of the differentiability of \(\overline{\hbox {v}}\) on \((0,\lambda _{\hbox {c}}/2)\), we require and prove a central limit theorem for the biased random walk. Additionally, we prove that the central limit theorem fails to hold for \(\lambda \ge \lambda _{\hbox {c}}/2\).     

Graph-Counting Polynomials for Oriented Graphs

If \(\mathcal{F}\) is a set of subgraphs F of a finite graph E we define a graph-counting polynomial \(p_\mathcal{F}(z)=\sum _{F\in \mathcal{F}}z^{|F|}\) In the present note we consider oriented graphs and discuss some cases where \(\mathcal{F}\) consists of unbranched subgraphs E. We find several situations where something can be said about the location of the zeros of \(p_\mathcal{F}\).     

Fully constrained Majorana neutrino mass matrices using $$\varvec{\varSigma (72\times 3)}$$

In 2002, two neutrino mixing ansatze having trimaximally mixed middle (\(\nu _2\)) columns, namely tri-chi-maximal mixing (\(\text {T}\chi \text {M}\)) and tri-phi-maximal mixing (\(\text {T}\phi \text {M}\)), were proposed. In 2012, it was shown that \(\text {T}\chi \text {M}\) with \(\chi =\pm \,\frac{\pi }{16}\) as well as \(\text {T}\phi \text {M}\) with \(\phi = \pm \,\frac{\pi }{16}\) leads to the solution, \(\sin ^2 \theta _{13} = \frac{2}{3} \sin ^2 \frac{\pi }{16}\), consistent with the latest measurements of the reactor mixing angle, \(\theta _{13}\). To obtain \(\text {T}\chi \text {M}_{(\chi =\pm \,\frac{\pi }{16})}\) and \(\text {T}\phi \text {M}_{(\phi =\pm \,\frac{\pi }{16})}\), the type I see-saw framework with fully constrained Majorana neutrino mass matrices was utilised. These mass matrices also resulted in the neutrino mass ratios, \(m_1:m_2:m_3=\frac{\left( 2+\sqrt{2}\right) }{1+\sqrt{2(2+\sqrt{2})}}:1:\frac{\left( 2+\sqrt{2}\right) }{-1+\sqrt{2(2+\sqrt{2})}}\). In this paper we construct a flavour model based on the discrete group \(\varSigma (72\times 3)\) and obtain the aforementioned results. A Majorana neutrino mass matrix (a symmetric \(3\times 3\) matrix with six complex degrees of freedom) is conveniently mapped into a flavon field transforming as the complex six-dimensional representation of \(\varSigma (72\times 3)\). Specific vacuum alignments of the flavons are used to arrive at the desired mass matrices.     

Stochastic Averaging Principle for Spatial Birth-and-Death Evolutions in the Continuum

We study a spatial birth-and-death process on the phase space of locally finite configurations \({\varGamma }^+ \times {\varGamma }^-\) over \({\mathbb {R}}^d\). Dynamics is described by an non-equilibrium evolution of states obtained from the Fokker-Planck equation and associated with the Markov operator \(L^+(\gamma ^-) + \frac{1}{\varepsilon }L^-\), \(\varepsilon > 0\). Here \(L^-\) describes the environment process on \({\varGamma }^-\) and \(L^+(\gamma ^-)\) describes the system process on \({\varGamma }^+\), where \(\gamma ^-\) indicates that the corresponding birth-and-death rates depend on another locally finite configuration \(\gamma ^- \in {\varGamma }^-\). We prove that, for a certain class of birth-and-death rates, the corresponding Fokker-Planck equation is well-posed, i.e. there exists a unique evolution of states \(\mu _t^{\varepsilon }\) on \({\varGamma }^+ \times {\varGamma }^-\). Moreover, we give a sufficient condition such that the environment is ergodic with exponential rate. Let \(\mu _{\mathrm {inv}}\) be the invariant measure for the environment process on \({\varGamma }^-\). In the main part of this work we establish the stochastic averaging principle, i.e. we prove that the marginal of \(\mu _t^{\varepsilon }\) onto \({\varGamma }^+\) converges weakly to an evolution of states on \({\varGamma }^+\) associated with the averaged Markov birth-and-death operator \({\overline{L}} = \int _{{\varGamma }^-}L^+(\gamma ^-)d \mu _{\mathrm {inv}}(\gamma ^-)\).     

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.

     

BCFT and OSFT moduli: an exact perturbative comparison

Starting from the pseudo-\({\mathcal {B}}_0\) gauge solution for marginal deformations in OSFT, we analytically compute the relation between the perturbative deformation parameter \(\tilde{\lambda }\) in the solution and the BCFT marginal parameter \(\lambda \), up to fifth order, by evaluating the Ellwood invariants. We observe that the microscopic reason why \(\tilde{\lambda }\) and \(\lambda \) are different is that the OSFT propagator renormalizes contact-term divergences differently from the contour deformation used in BCFT.     

Phase Transition and Uniqueness of Levelset Percolation

The main purpose of this paper is to introduce and establish basic results of a natural extension of the classical Boolean percolation model (also known as the Gilbert disc model). We replace the balls of that model by a positive non-increasing attenuation function \(l:(0,\infty ) \rightarrow [0,\infty )\) to create the random field \(\Psi (y)=\sum _{x\in \eta }l(|x-y|),\) where \(\eta \) is a homogeneous Poisson process in \({\mathbb {R}}^d.\) The field \(\Psi \) is then a random potential field with infinite range dependencies whenever the support of the function l is unbounded. In particular, we study the level sets \(\Psi _{\ge h}(y)\) containing the points \(y\in {\mathbb {R}}^d\) such that \(\Psi (y)\ge h.\) In the case where l has unbounded support, we give, for any \(d\ge 2,\) a necessary and sufficient condition on l for \(\Psi _{\ge h}(y)\) to have a percolative phase transition as a function of h. We also prove that when l is continuous then so is \(\Psi \) almost surely. Moreover, in this case and for \(d=2,\) we prove uniqueness of the infinite component of \(\Psi _{\ge h}\) when such exists, and we also show that the so-called percolation function is continuous below the critical value \(h_c\).     

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\).     

Modularity of logarithmic parafermion vertex algebras

The parafermionic cosets \(\mathsf {C}_{k} = {\text {Com}} ( \mathsf {H} , \mathsf {L}_{k}(\mathfrak {sl}_{2}) )\) are studied for negative admissible levels k, as are certain infinite-order simple current extensions \(\mathsf {B}_{k}\) of \(\mathsf {C}_{k}\). Under the assumption that the tensor theory considerations of Huang, Lepowsky and Zhang apply to \(\mathsf {C}_{k}\), irreducible \(\mathsf {C}_{k}\)- and \(\mathsf {B}_{k}\)-modules are obtained from those of \(\mathsf {L}_{k}(\mathfrak {sl}_{2})\). Assuming the validity of a certain Verlinde-type formula likewise gives the Grothendieck fusion rules of these irreducible modules. Notably, there are only finitely many irreducible \(\mathsf {B}_{k}\)-modules. The irreducible \(\mathsf {C}_{k}\)- and \(\mathsf {B}_{k}\)-characters are computed and the latter are shown, when supplemented by pseudotraces, to carry a finite-dimensional representation of the modular group. The natural conjecture then is that the \(\mathsf {B}_{k}\) are \(C_2\)-cofinite vertex operator algebras.     

Two Higgs doublet models and $$b\rightarrow s$$ exclusive decays

We computed the leading order Wilson coefficients relevant to all the exclusive \(b\rightarrow s\ell ^+\ell ^-\) decays in the framework of the two Higgs doublet model (2HDM) with a softly broken \(\mathbb {Z}_2\) symmetry by including the \(\mathcal {O}(m_b)\) corrections. We elucidate the issue of appropriate matching between the full and the effective theory when dealing with the (pseudo-)scalar operators for which keeping the external momenta different from zero is necessary. We then make a phenomenological analysis by using the measured \({\mathcal {B}}(B_s\rightarrow \mu ^+\mu ^-)\) and \({\mathcal {B}}(B\rightarrow K \mu ^+\mu ^-)_{\mathrm {high}-q^2}\), for which the hadronic uncertainties are well controlled, and we discuss their impact on various types of 2HDM. A brief discussion of the decays with \(\tau \)-leptons in the final state is provided too.