{{\psi(3S)}} and {{\Upsilon(5S)}} Originating in Heavy Meson Molecules: A Coupled Channel Analysis Based on an Effective Vector Quark–Quark Interaction

In our previous coupled channel analysis based on the Cornell effective quark–quark interaction, it was indicated that the ${\psi(3S)}$ solution corresponding to ${\psi(4040)}$ originates from a ${{\rm D}^{^{*}}\overline{{\rm D}}^{*}}$ channel state. In this article, we report on a simultaneous analysis of the ${\psi}$ - and ${\Upsilon}$ -family states. The most conspicuous outcome is a finding that the ${\Upsilon(5S)}$ solution corresponding to ${\Upsilon(10860)}$ originates from a ${{\rm B}^{*}\overline{{\rm B}}^{*}}$ channel state, very much like ${\psi(3S)}$ . Some other characteristics of the result, including the induced very large S–D mixing and relation of some of the solutions with newly observed heavy quarkonia-like states are discussed.     

Quantum billiards in multidimensional models with branes

A gravitational $D$ -dimensional model with $l$ scalar fields and several forms is considered. When a cosmological-type diagonal metric is chosen, an electromagnetic composite brane ansatz is adopted and certain restrictions on the branes are imposed; the conformally covariant Wheeler–DeWitt (WDW) equation for the model is studied. Under certain restrictions asymptotic solutions to WDW equation are found in the limit of the formation of the billiard walls which reduce the problem to the so-called quantum billiard on the $(D+ l -2)$ -dimensional Lobachevsky space. Two examples of quantum billiards are considered. The first one deals with $9$ -dimensional quantum billiard for $D = 11$ model with $330$ four-forms which mimic space-like $M2$ - and $M5$ -branes of $D=11$ supergravity. The second one deals with the $9$ -dimensional quantum billiard for $D =10$ gravitational model with one scalar field, $210$ four-forms and $120$ three-forms which mimic space-like $D2$ -, $D4$ -, $FS1$ - and $NS5$ -branes in $D = 10$ $II A$ supergravity. It is shown that in both examples wave functions vanish in the limit of the formation of the billiard walls (i.e. we get a quantum resolution of the singularity for $11D$ model) but magnetic branes could not be neglected in calculations of quantum asymptotic solutions while they are irrelevant for classical oscillating behavior when all $120$ electric branes are present.     

An Algebraic Construction of Boundary Quantum Field Theory

We build up local, time translation covariant Boundary Quantum Field Theory nets of von Neumann algebras ${\mathcal A_V}$ on the Minkowski half-plane M ₊ starting with a local conformal net ${\mathcal A}$ of von Neumann algebras on ${\mathbb R}$ and an element V of a unitary semigroup ${\mathcal E(\mathcal A)}$ associated with ${\mathcal A}$ . The case V?=?1 reduces to the net ${\mathcal A_+}$ considered by Rehren and one of the authors; if the vacuum character of ${\mathcal A}$ is summable, ${\mathcal A_V}$ is locally isomorphic to ${\mathcal A_+}$ . We discuss the structure of the semigroup ${\mathcal E(\mathcal A)}$ . By using a one-particle version of Borchers theorem and standard subspace analysis, we provide an abstract analog of the Beurling-Lax theorem that allows us to describe, in particular, all unitaries on the one-particle Hilbert space whose second quantization promotion belongs to ${\mathcal E(\mathcal A^{(0)})}$ with ${\mathcal A^{(0)}}$ the U(1)-current net. Each such unitary is attached to a scattering function or, more generally, to a symmetric inner function. We then obtain families of models via any Buchholz-Mack-Todorov extension of ${\mathcal A^{(0)}}$ . A further family of models comes from the Ising model.     

Charmonium rescattering effects in \psi(2S) \rightarrow \gamma \eta_{c}^{} decay

Charmonium rescattering effects in the M1 transition of $ \psi$ (2S) $ \rightarrow$ $ \gamma$ $ \eta_{c}^{}$ are investigated by modeling a $ \chi_{{cJ}}^{}$ or J/ $ \psi$ rescattering into a $ \eta_{c}^{}$ final state. The absorptive and dispersive part of the transition amplitudes for the rescattering loops of $ \eta$ $ \psi$ ( $ \gamma^{{\ast}}_{}$ ) and $ \gamma$ $ \chi$ ( $ \psi$ ) are separately evaluated. The numerical results show that the contribution from the $ \gamma$ $ \chi$ ( $ \psi$ ) rescattering process is negligible. Compared with the virtual D $ \bar{{D}}$ (D ^*) rescattering processes, the $ \eta$ $ \psi$ ( $ \gamma^{{\ast}}_{}$ ) process may be regarded as the next-leading order of the hadronic loop mechanism, which only offers the partial decay width of ～ 0.045 keV to the $ \psi$ (2S) $ \rightarrow$ $ \gamma$ $ \eta_{c}^{}$ .     

Large N approach to Kaon decays and mixing 28 years later: \Delta I=1/2 rule, \hat{B}_K,and \Delta M_K

We review and update our results for $K\rightarrow \pi \pi $ decays and $K^0$ – $\bar{K}^0$ mixing obtained by us in the 1980s within an analytic approximate approach based on the dual representation of QCD as a theory of weakly interacting mesons for large $N$ , where $N$ is the number of colors. In our analytic approach the Standard Model dynamics behind the enhancement of $\hbox {Re}A_0$ and suppression of $\hbox {Re}A_2$ , the so-called $\Delta I=1/2$ rule for $K\rightarrow \pi \pi $ decays, has a simple structure: the usual octet enhancement through the long but slow quark–gluon renormalization group evolution down to the scales $\mathcal{O}(1\, {\hbox { GeV}})$ is continued as a short but fast meson evolution down to zero momentum scales at which the factorization of hadronic matrix elements is at work. The inclusion of lowest-lying vector meson contributions in addition to the pseudoscalar ones and of Wilson coefficients in a momentum scheme improves significantly the matching between quark–gluon and meson evolutions. In particular, the anomalous dimension matrix governing the meson evolution exhibits the structure of the known anomalous dimension matrix in the quark–gluon evolution. While this physical picture did not yet emerge from lattice simulations, the recent results on $\hbox {Re}A_2$ and $\hbox {Re}A_0$ from the RBC-UKQCD collaboration give support for its correctness. In particular, the signs of the two main contractions found numerically by these authors follow uniquely from our analytic approach. Though the current–current operators dominate the $\Delta I=1/2$ rule, working with matching scales $\mathcal{O}(1 \, {\hbox { GeV}})$ we find that the presence of QCD-penguin operator $Q_6$ is required to obtain satisfactory result for $\hbox {Re}A_0$ . At NLO in $1/N$ we obtain $R=\hbox {Re}A_0/\hbox {Re}A_2= 16.0\pm 1.5$ which amounts to an order of magnitude enhancement over the strict large $N$ limit value $\sqrt{2}$ . We also update our results for the parameter $\hat{B}_K$ , finding $\hat{B}_K=0.73\pm 0.02$ . The smallness of $1/N$ corrections to the large $N$ value $\hat{B}_K=3/4$ results within our approach from an approximate cancelation between pseudoscalar and vector meson one-loop contributions. We also summarize the status of $\Delta M_K$ in this approach.     

Formation and crystallization kinetics of Nd–Fe–B-based bulk amorphous alloy

In order to improve the glass-forming ability (GFA) of Nd–Fe–B ternary alloys to obtain fully amorphous bulk Nd–Fe–B-based alloy, the effects of Mo and Y doping on GFA of the alloys were investigated. It was found that the substitution of Mo for Fe and Y for Nd enhanced the GFA of the Nd–Y–Fe–Mo–B alloys. It was also revealed that the GFA of the samples was optimized by 4 at.% Mo doping and increased with the Y content. The fully amorphous structures were all formed in the Nd $_{6-{x}}$ Y $_{{x}}$ Fe $_{68}$ Mo $_{4}$ B $_{22}$ (x $=$ 1–5) alloy rods with 1.5 mm-diameter. After subsequent crystallization, the devitrified Nd $_{3}$ Y $_{3}$ Fe $_{68}$ Mo $_{4}$ B $_{22}$ alloy rod exhibited a uniform distribution of grains with a coercivity of 364.1 kA/m. The crystallization behavior of Nd $_{3}$ Y $_{3}$ Fe $_{68}$ Mo $_{4}$ B $_{22}$ BMG was investigated in isothermal situation. The Avrami exponent n determined by JAM plot is lower than 2.5, implying that the crystallization is mainly governed by a growth of particles with decreasing nucleation rate.     

Influence of temperature on the electric,dielectric and AC conductivity properties of nano-crystalline zinc substituted cobalt ferrite synthesized by solution combustion technique

Cobalt–zinc nanoferrites with formulae Co $_{1-x}$ Zn $_{x}$ Fe $_{2}$ O $_{4}$ , where x = 0.0, 0.1, 0.2 and 0.3, have been synthesized by solution combustion technique. The variation of DC resistivity with temperature shows the semiconducting behavior of all nanoferrites. The dielectric properties such as dielectric constant ( $\varepsilon $ ’) and dielectric loss tangent (tan $\delta )$ are investigated as a function of temperature and frequency. Dielectric constant and loss tangent are found to be increasing with an increase in temperature while with an increase in frequency both, $\varepsilon $ ’ and tan $\delta $ , are found to be decreasing. The dielectric properties have been explained on the basis of space charge polarization according to Maxwell–Wagner’s two-layer model and the hopping of charge between Fe $^{2+}$ and Fe $^{3+}$ . Further, a very high value of dielectric constant and a low value of tan $\delta $ are the prime achievements of the present work. The AC electrical conductivity ( $\sigma _\mathrm{AC})$ is studied as a function of temperature as well as frequency and $\sigma _\mathrm{AC}$ is observed to be increasing with the increase in temperature and frequency.     

Mössbauer spectroscopy in the investigation of new mineral–related materials

New materials based on the composition of the mineral schafarzikite, FeSb $_{2}\textit {O}_{4}$ , have been synthesised. $^{57}$ Fe- and $^{121}$ Sb- Mössbauer spectroscopy shows that iron is present as Fe $^{2+}$ and that antimony is present as Sb $^{3+}$ . The presence of Pb $^{2+}$ on the antimony sites in materials of composition FeSb $_{1.5}$ Pb $_{0.5}\textit {O}_{4}$ induces partial oxidation of Fe $^{2+}_{}$ to Fe $^{3+}$ . The quasi-one-dimensional magnetic structure of schafarzikite is retained in FeSb $_{1.5}$ Pb $_{0.5}\textit {O}_{4}$ and gives rise to weakly coupled non-magnetic Fe $^{2+}$ ions coexisting with Fe $^{3+}$ ions in a magnetically ordered state. A similar model can be applied to account for the spectra recorded from the compound Co $_{0.5}$ Fe $_{0.5}$ Sb $_{1.5}$ Pb $_{0.5}\textit {O}_{4}$ .     

Fomenko–Mischenko Theory, Hessenberg Varieties, and Polarizations

The symmetric algebra ${S(\mathfrak{g})}$ over a Lie algebra ${\mathfrak{g}}$ has the structure of a Poisson algebra. Assume ${\mathfrak{g}}$ is complex semisimple. Then results of Fomenko–Mischenko (translation of invariants) and Tarasov construct a polynomial subalgebra ${{\mathcal {H}} = {\mathbb C}[q_1,\ldots,q_b]}$ of ${S(\mathfrak{g})}$ which is maximally Poisson commutative. Here b is the dimension of a Borel subalgebra of ${\mathfrak{g}}$ . Let G be the adjoint group of ${\mathfrak{g}}$ and let ? = rank ${\mathfrak{g}}$ . Using the Killing form, identify ${\mathfrak{g}}$ with its dual so that any G-orbit O in ${\mathfrak{g}}$ has the structure (KKS) of a symplectic manifold and ${S(\mathfrak{g})}$ can be identified with the affine algebra of ${\mathfrak{g}}$ . An element ${x\in \mathfrak{g}}$ will be called strongly regular if ${\{({\rm d}q_i)_x\},\,i=1,\ldots,b}$ , are linearly independent. Then the set ${\mathfrak{g}^{\rm{sreg}}}$ of all strongly regular elements is Zariski open and dense in ${\mathfrak{g}}$ and also ${\mathfrak{g}^{\rm{sreg}}\subset \mathfrak{g}^{\rm{ reg}}}$ where ${\mathfrak{g}^{\rm{reg}}}$ is the set of all regular elements in ${\mathfrak{g}}$ . A Hessenberg variety is the b-dimensional affine plane in ${\mathfrak{g}}$ , obtained by translating a Borel subalgebra by a suitable principal nilpotent element. Such a variety was introduced in Kostant (Am J Math 85:327–404, 1963). Defining Hess to be a particular Hessenberg variety, Tarasov has shown that ${{\rm{Hess}}\subset \mathfrak{g}^{\rm{sreg}}}$ . Let R be the set of all regular G-orbits in ${\mathfrak{g}}$ . Thus if ${O\in R}$ , then O is a symplectic manifold of dimension 2n where n = b ? ?. For any ${O\in R}$ let ${O^{\rm{sreg}} = \mathfrak{g}^{\rm{sreg}} \cap O}$ . One shows that O ^sreg is Zariski open and dense in O so that O ^sreg is again a symplectic manifold of dimension 2n. For any ${O\in R}$ let ${{\rm{Hess}}(O) = {\rm{Hess}}\cap O}$ . One proves that Hess(O) is a Lagrangian submanifold of O ^sreg and that $${\rm{Hess}} = \sqcup_{O\in R}{\rm{Hess}}(O).$$ The main result of this paper is to show that there exists simultaneously over all ${O\in R}$ , an explicit polarization (i.e., a “fibration” by Lagrangian submanifolds) of O ^sreg which makes O ^sreg simulate, in some sense, the cotangent bundle of Hess(O).     

Finite W-Superalgebras and Truncated Super Yangians

Let ${Y_{m|n}^{\ell}}$ be the super Yangian of general linear Lie superalgebra for ${\mathfrak{gl}_{m|n}}$ . Let ${e \in \mathfrak{gl}_{m\ell|n\ell}}$ be a “rectangular” nilpotent element and ${\mathcal{W}_e}$ be the finite W-superalgebra associated to e. We show that ${Y_{m|n}^{\ell}}$ is isomorphic to ${\mathcal{W}_e}$ .     

{\Xi} Hypernucler States Predicted by the New Interaction Model ESC08c

The features of the new interaction model ESC08c in ${\Lambda N}$ , ${\Sigma N}$ and ${\Xi N}$ channels are demonstrated single hyperon potentials ${U_Y(Y=\Lambda, \Sigma, \Xi)}$ in nuclear matter on the basis of the G-matrix theory. (K ^?, K ⁺) productions of ${\Xi}$ hypernuclei are studied with ${\Xi}$ -nucleus folding potentials.     

Asymptotic Behavior for a Version of Directed Percolation on the Triangular Lattice

We consider a version of directed bond percolation on the triangular lattice such that vertical edges are directed upward with probability $y$ , diagonal edges are directed from lower-left to upper-right or lower-right to upper-left with probability $d$ , and horizontal edges are directed rightward with probabilities $x$ and one in alternate rows. Let $\tau (M,N)$ be the probability that there is at least one connected-directed path of occupied edges from $(0,0)$ to $(M,N)$ . For each $x \in [0,1]$ , $y \in [0,1)$ , $d \in [0,1)$ but $(1-y)(1-d) \ne 1$ and aspect ratio $\alpha =M/N$ fixed for the triangular lattice with diagonal edges from lower-left to upper-right, we show that there is an $\alpha _c = (d-y-dy)/[2(d+y-dy)] + [1-(1-d)^2(1-y)^2x]/[2(d+y-dy)^2]$ such that as $N \rightarrow \infty $ , $\tau (M,N)$ is $1$ , $0$ and $1/2$ for $\alpha > \alpha _c$ , $\alpha < \alpha _c$ and $\alpha =\alpha _c$ , respectively. A corresponding result is obtained for the triangular lattice with diagonal edges from lower-right to upper-left. We also investigate the rate of convergence of $\tau (M,N)$ and the asymptotic behavior of $\tau (M_N^-,N)$ and $\tau (M_N^+ ,N)$ where $M_N^-/N\uparrow \alpha _c$ and $M_N^+/N\downarrow \alpha _c$ as $N\uparrow \infty $ .     

Radiative open charm decay of the Y(3940) , Z(3930) , X(4160) resonances

We determine the radiative decay amplitudes for the decay into D^* and $ \bar{{D}}$ $ \gamma$ , or D ^* _s and $ \bar{{D}}_{s}^{}$ $ \gamma$ of some of the charmonium-like states classified as X , Y , Z resonances, plus some other hidden charm states which are dynamically generated from the interaction of vector mesons with charm. The mass distributions as a function of the $ \bar{{D}}$ $ \gamma$ or $ \bar{{D}}_{s}^{}$ $ \gamma$ invariant mass show a peculiar behavior as a consequence of the D ^* $ \bar{{D}}^{*}_{}$ nature of these states. The experimental search of these magnitudes can shed light on the nature of these states.     

Influence of assortativity and degree-preserving rewiring on the spectra of networks

Newman’s measure for (dis)assortativity, the linear degree correlation coefficient $\rho _{D}$ , is reformulated in terms of the total number N _k of walks in the graph with k hops. This reformulation allows us to derive a new formula from which a degree-preserving rewiring algorithm is deduced, that, in each rewiring step, either increases or decreases $\rho _{D}$ conform our desired objective. Spectral metrics (eigenvalues of graph-related matrices), especially, the largest eigenvalue $\lambda _{1}$ of the adjacency matrix and the algebraic connectivity $\mu _{N-1}$ (second-smallest eigenvalue of the Laplacian) are powerful characterizers of dynamic processes on networks such as virus spreading and synchronization processes. We present various lower bounds for the largest eigenvalue $\lambda _{1}$ of the adjacency matrix and we show, apart from some classes of graphs such as regular graphs or bipartite graphs, that the lower bounds for $\lambda _{1}$ increase with $\rho _{D}$ . A new upper bound for the algebraic connectivity $\mu _{N-1}$ decreases with $\rho _{D}$ . Applying the degree-preserving rewiring algorithm to various real-world networks illustrates that (a) assortative degree-preserving rewiring increases $\lambda _{1}$ , but decreases $\mu _{N-1}$ , even leading to disconnectivity of the networks in many disjoint clusters and that (b) disassortative degree-preserving rewiring decreases $\lambda _{1}$ , but increases the algebraic connectivity, at least in the initial rewirings.     

Analysis of the vertexes \Xi_{Q}^{*}′QV , \Sigma_{Q}^{*} \Sigma_{Q}^{}V and radiative decays \Xi_{Q}^{*} \rightarrow ′Q \gamma , \Sigma_{Q}^{*} \rightarrow \Sigma_{Q}^{} \gamma

In this article, we study the vertexes $ \Xi_{Q}^{*}$ ′_Q V and $ \Sigma_{Q}^{*}$ $ \Sigma_{Q}^{}$ V with the light-cone QCD sum rules, then assume the vector meson dominance of the intermediate $ \phi$ (1020) , $ \rho$ (770) and $ \omega$ (782) , and calculate the radiative decays $ \Xi_{Q}^{*}$ $ \rightarrow$ ′_Q $ \gamma$ and $ \Sigma_{Q}^{*}$ $ \rightarrow$ $ \Sigma_{Q}^{}$ $ \gamma$ .     

Reaction pp \rightarrow pp \pi \pi \pi as a background for hadronic decays of the \eta^{{\prime}}_{} meson

Isospin violating hadronic decays of the $ \eta$ and $ \eta{^\prime}$ mesons into 3 $ \pi$ mesons are driven by a term in the QCD Lagrangian proportional to the mass difference of the d and u quarks. The source giving large yield of the mesons for such decay studies are pp interactions close to the respective kinematical thresholds. The most important physics background for $ \eta$ , $ \eta{^\prime}$ $ \rightarrow$ $ \pi$ $ \pi$ $ \pi$ is coming from direct three-pion production reactions. In case of the $ \eta$ meson the background for the decays is relatively low ( $ \approx$ 10% . The purpose of this article is to provide an estimate of the direct pion production background for the $ \eta{^\prime}$ $ \rightarrow$ 3 $ \pi$ decays. Using the inclusive data from the COSY-11 experiment we have extracted the differential cross-section for the pp $ \rightarrow$ pp -multipion production reactions with the invariant mass of the pions equal to the $ \eta{^\prime}$ meson mass and estimated an upper limit for the signal to background ratio for studies of the $ \eta{^\prime}$ $ \rightarrow$ $ \pi^{+}_{}$ $ \pi^{-}_{}$ $ \pi^{0}_{}$ decay.     

The Lie–Rinehart Universal Poisson Algebra of Classical and Quantum Mechanics

The Lie–Rinehart algebra of a (connected) manifold ${\mathcal {M}}$ , defined by the Lie structure of the vector fields, their action and their module structure over ${C^\infty({\mathcal {M}})}$ , is a common, diffeomorphism invariant, algebra for both classical and quantum mechanics. Its (noncommutative) Poisson universal enveloping algebra ${\Lambda_{R}({\mathcal {M}})}$ , with the Lie–Rinehart product identified with the symmetric product, contains a central variable (a central sequence for non-compact ${{\mathcal {M}}}$ ) ${Z}$ which relates the commutators to the Lie products. Classical and quantum mechanics are its only factorial realizations, corresponding to Z  =  i z, z  =  0 and ${z = \hbar}$ , respectively; canonical quantization uniquely follows from such a general geometrical structure. For ${z =\hbar \neq 0}$ , the regular factorial Hilbert space representations of ${\Lambda_{R}({\mathcal{M}})}$ describe quantum mechanics on ${{\mathcal {M}}}$ . For z  =  0, if Diff( ${{\mathcal {M}}}$ ) is unitarily implemented, they are unitarily equivalent, up to multiplicity, to the representation defined by classical mechanics on ${{\mathcal {M}}}$ .     

On the production of \pi^{+}_{} \pi^{+}_{} pairs in pp collisions at 0.8 GeV

Data accumulated recently for the exclusive measurement of the pp $ \rightarrow$ pp $ \pi^{+}_{}$ $ \pi^{-}_{}$ reaction at a beam energy of 0.793GeV using the COSY-TOF spectrometer have been analyzed with respect to possible events from the pp $ \rightarrow$ nn $ \pi^{+}_{}$ $ \pi^{+}_{}$ reaction channel. The latter is expected to be the only $ \pi$ $ \pi$ production channel, which contains no major contributions from resonance excitation close to threshold and hence should be a good testing ground for chiral dynamics in the $ \pi$ $ \pi$ production process. No single event has been found, which meets all conditions for being a candidate for the pp $ \rightarrow$ nn $ \pi^{+}_{}$ $ \pi^{+}_{}$ reaction. This gives an upper limit for the cross-section of 0.16μb (90% C.L.), which is more than an order of magnitude smaller than the cross-sections of the other two-pion production channels at the same incident energy.     

Analysis of the {3\over 2}^{+} heavy and doubly heavy baryon states with QCD sum rules

In this article, we study the ${3\over 2}^{+}$ heavy and doubly heavy baryon states $\varXi^{*}_{cc}$ , $\varOmega^{*}_{cc}$ , $\varXi^{*}_{bb}$ , $\varOmega^{*}_{bb}$ , $\varSigma_{c}^{*}$ , $\varXi_{c}^{*}$ , $\varOmega_{c}^{*}$ , $\varSigma_{b}^{*}$ , $\varXi_{b}^{*}$ and $\varOmega_{b}^{*}$ by subtracting the contributions from the corresponding ${3\over 2}^{-}$ heavy and doubly heavy baryon states with the QCD sum rules, and we make reasonable predictions for their masses.     

Pyrid-2-yl and 2-CyanoPhenyl fused heterocyclic compounds as human P2X_{3} inhibitors: a combined approach based on homology modelling,docking and QSAR analysis

P2X receptors are hetero-oligomeric proteins that function as membrane ion channels and are gated by extracellular ATP. The hP2X $_{3}$ subunit is a constituent of the channels on a subset of sensory neurons involved in pain signaling, where ATP released by damaged and inflamed tissue can initiate action potentials. Hence, the inhibition of ATP-activated P2X $_{3}$ receptor is an exciting approach for the treatment of inflammatory and neuropathic pain. Recently, the crystal structures of zebrafish P2X $_{4}$ (zP2X $_{4})$ were obtained in closed, apo state (PDB ID: 3I5D) and ATP-bound, open state (PDB ID: 4DW1). These structures were used to develop a homology model of human P2X $_{3}$ (hP2X $_{3})$ in order to identify through docking studies, the binding modes of known P2X $_{3}$ inhibitors and their key active site interactions, along with a pharmacophore-based 3D-QSAR model for a series of 136 Pyrid-2-yl and 2-CyanoPhenyl fused heterocyclic compounds. These 3D-QSAR models have been developed with different combinations of training and test set divisions obtained by random separation, Jarvis–Patrick clustering, K-means clustering and sphere exclusion methods. The best predictive 3D-QSAR model resulted in training set R $^{2 }$ of 0.75, internal test set Q $^{2}$ of 0.74, Pearson-R value of 0.87 and root mean square error of 0.37. The information generated by the pharmacophore model and docking analyses using the homology model provides valuable clues to design novel potent hP2X $_{3}$ inhibitors.