Search

We study the conformality loss of theories with long-range interactions. We consider the $O(2)\times O(N)$ multiscalar model with coupling $r^{-d-\delta}$ in $d=4-\epsilon$ dimension. We compute the critical exponents of the long-range fixed points (LRFPs) to three loops. The phase diagram of the model is dominated by two processes: the short-range crossover and merger-annihilation of LRFPs. The two processes intersect at the lower edge of the conformal window, below which the LRFPs disappear into the complex plane. We propose a novel scenario for the short-range crossover of complex LRFPs, in which the short-range crossover occurs on a vertical line in the complex plane of $\delta$. The complex marginal operator generates renormalization group flow on the transition line and significantly enriches the short-range crossover of complex LRFPs., Comment: 4+6 pages, 2 figures

Satellite altimetry has been widely utilized to monitor global sea surface dynamics, enabling investigation of upper ocean variability from basin-scale to localized eddy ranges. However, the sparse spatial resolution of observational altimetry limits our understanding of oceanic submesoscale variability, prevalent at horizontal scales below 0.25o resolution. Here, we introduce a state-of-the-art generative diffusion model to train high-resolution sea surface height (SSH) reanalysis data and demonstrate its advantage in observational SSH downscaling over the eddy-rich Kuroshio Extension region. The diffusion-based model effectively downscales raw satellite-interpolated data from 0.25o resolution to 1/16o, corresponding to approximately 12-km wavelength. This model outperforms other high-resolution reanalysis datasets and neural network-based methods. Also, it successfully reproduces the spatial patterns and power spectra of satellite along-track observations. Our diffusion-based results indicate that eddy kinetic energy at horizontal scales less than 250 km has intensified significantly since 2004 in the Kuroshio Extension region. These findings underscore the great potential of deep learning in reconstructing satellite altimetry and enhancing our understanding of ocean dynamics at eddy scales., Comment: 28 pages, 7 figures, and 1 table

We study the $\mathbb{Z}_2$ and $U(1)$ Abelian lattice gauge theories using a bootstrap method, in which the loop equations and positivity conditions are employed for Wilson loops with lengths $L\leqslant L_{\textrm{max}}$ to derive two-sided bounds on the Wilson loop averages. We address a fundamental question that whether the constraints from loop equations and positivity are strong enough to solve lattice gauge theories. We answer this question by bootstrapping the 2D $U(1)$ lattice gauge theory. We show that with sufficiently large $L_{\textrm{max}}=60$, the two-sided bounds provide estimates for the plaquette averages with precision near $10^{-8}$ or even higher, suggesting the bootstrap constraints are sufficient to numerically pin down this theory. We compute the bootstrap bounds on the plaquette averages in the 3D $\mathbb{Z}_2$ and $U(1)$ lattice gauge theories with $L_{\textrm{max}}=16$. In the regions with weak or strong coupling, the two-sided bootstrap bounds converge quickly and coincide with the perturbative results to high precision. The bootstrap bounds are well consistent with the Monte Carlo results in the nonperturbative region. We observe interesting connections between the bounds generated by the bootstrap computations and the Griffiths' inequalities. We present results towards bootstrapping the string tension and glueball mass in Abelian lattice gauge theories., Comment: 40 pages, 6 figures; v2: typos corrected, refs added, more related data

We initiate the analytical functional bootstrap study of conformal field theories with large $N$ limits. In this first paper we particularly focus on the 1D $O(N)$ vector bootstrap. We obtain a remarkably simple bootstrap equation from the $O(N)$ vector crossing equations in the large $N$ limit. The bootstrap bound is saturated by the generalized free field theory. We study the analytical extremal functionals of this crossing equation, for which the total positivity of the $SL(2,\mathbb{R})$ conformal block plays a critical role. We prove the $SL(2,\mathbb{R})$ conformal block is totally positive for large scaling dimension $\Delta$ and show that the total positivity is violated below a critical value $\Delta_{\textrm{TP}}^*\approx 0.32315626$. The conformal block forms a surprisingly sophisticated mathematical structure, which for instance can violate total positivity at the order $10^{-5654}$ for a normal value $\Delta=0.1627$! We construct a series of analytical functionals $\{\alpha_M\}$ which satisfy the bootstrap positive conditions up to a range $\Delta\leqslant \Lambda_M$. The functionals $\{\alpha_M\}$ have a trivial large $M$ limit. Surprisingly, due to total positivity, they can approach the large $M$ limit in a way consistent with the bootstrap positive conditions for arbitrarily high $\Lambda_M$, therefore proving the bootstrap bound analytically. Our result provides a concrete example to illustrate how the analytical properties of the conformal block lead to nontrivial bootstrap bounds. We expect this work paves the way for large $N$ analytical functional bootstrap in higher dimensions., Comment: v1: 48 pages, 9 figures; v2: refs added, typos corrected, added an estimation of the critical value for total positivity

Background \& purpose: The recent emergence of neural networks models for the analysis of breast images has been a breakthrough in computer aided diagnostic. This approach was not yet developed in Contrast Enhanced Spectral Mammography (CESM) where access to large databases is complex. This work proposes a deep-learning-based Computer Aided Diagnostic development for CESM recombined images able to detect lesions and classify cases. Material \& methods: A large CESM diagnostic dataset with biopsy-proven lesions was collected from various hospitals and different acquisition systems. The annotated data were split on a patient level for the training (55%), validation (15%) and test (30%) of a deep neural network with a state-of-the-art detection architecture. Free Receiver Operating Characteristic (FROC) was used to evaluate the model for the detection of 1) all lesions, 2) biopsied lesions and 3) malignant lesions. ROC curve was used to evaluate breast cancer classification. The metrics were finally compared to clinical results. Results: For the evaluation of the malignant lesion detection, at high sensitivity (Se>0.95), the false positive rate was at 0.61 per image. For the classification of malignant cases, the model reached an Area Under the Curve (AUC) in the range of clinical CESM diagnostic results. Conclusion: This CAD is the first development of a lesion detection and classification model for CESM images. Trained on a large dataset, it has the potential to be used for helping the management of biopsy decision and for helping the radiologist detecting complex lesions that could modify the clinical treatment.

We present the results of a conformal bootstrap study of the presumed unitary IR fixed point of quantum electrodynamics in three dimensions (QED$_3$) coupled to $N_f=4$ two-component Dirac fermions. Specifically, we study the four-point correlators of the $SU(4)$ adjoint fermion bilinear $r$ and the monopole of lowest topological charge $\mathcal{M}_{1/2}$. Most notably, the scaling dimensions of the fermion bilinear $r$ and the monopole $\mathcal{M}_{1/2}$ are found to be constrained into a closed island with a combination of spectrum assumptions inspired by the $1/N_f$ perturbative results as well as a novel interval positivity constraint on the next-lowest-charge monopole $\mathcal{M}_1$. Bounds in this island on the $SU(4)$ and topological $U(1)_t$ conserved current central charges $c_J$, $c_J^t$, as well as on the stress tensor central charge $c_T$, are comfortably consistent with the perturbative results. Together with the scaling dimensions, this suggests that a part of estimates from the $1/N_f$ expansion -- even at $N_f=4$ -- provide a self-consistent solution to the bootstrap crossing relations, despite some of our assumptions not being strictly justified., Comment: v1: 66 pages, 14 figures, 7 tables; v2: revtex-2 version, typos corrected

We study the IR phase of three dimensional quantum electrodynamics (QED$_3$) coupled to $N_f=2$ flavors of two-component Dirac fermions, which has been controversial for decades. This theory has been proposed to be self-dual with symmetry enhancement $(SU(2)_f\times U(1)_t )/\mathbb{Z}_2\rightarrow O(4)$ at the IR fixed point. We focus on the four-point correlator of monopole operators with unit topological charge of $U(1)_t$. We illustrate the $O(4)\rightarrow SU(2)_f\times U(1)_t $ branching rules based on an $O(4)$ symmetric positive structure in the monopole four-point crossing equations. We use conformal bootstrap method to derive nonperturbative constraints on the CFT data and test the conformality and self-duality of $N_f=2$ QED$_3$. In particular we find the CFT data obtained from previous lattice simulations can be ruled out by introducing irrelevant assumptions in the spectrum, indicating the IR phase of $N_f=2$ QED$_3$ is not conformal., Comment: 4+5 pages, 2+1 figures; v2: references added, typos corrected; v3: 3+1 figures, explanation on the parity charges added, matches version to appear in PLB

A program of wide interest in modern conformal bootstrap studies is to numerically solve general conformal field theories, based on a critical assumption that the dynamics is encoded in the conformal four-point crossing equations and positivity condition. In this letter we propose and verify a novel algebraic property of the crossing equations which provides strong restriction for this program. We show for various types of symmetries $\cal G$, the crossing equations can be linearly converted into the $SO(N)$ vector crossing equations associated with the $SO(N)\rightarrow \cal G$ branching rules and the transformations satisfy positivity condition. The dynamics constrained by the $\cal G$-symmetric crossing equations combined with positivity condition degenerates to the $SO(N)$ symmetric cases, while the non-$SO(N)$ symmetric theories are not directly solvable without introducing the $SO(N)$ symmetry breaking assumptions on the spectrum., Comment: 16 pages, no figures; v2: 12 pages, no figures, SO(N)-ization of WZW model added, discussion improved, to match the published version

Infrared fixed points of gauge theories provide intriguing targets for the modern conformal bootstrap program. In this work we provide some preliminary evidence that a family of gauged fermionic CFTs saturate bootstrap bounds and can potentially be solved with the conformal bootstrap. We start by considering the bootstrap for $SO(N)$ vector 4-point functions in general dimension $D$. In the large $N$ limit, upper bounds on the scaling dimensions of the lowest $SO(N)$ singlet and traceless symmetric scalars interpolate between two solutions at $\Delta =D/2-1$ and $\Delta =D-1$ via generalized free field theory. In 3D the critical $O(N)$ vector models are known to saturate the bootstrap bounds and correspond to the kinks approaching $\Delta =1/2$ at large $N$. We show that the bootstrap bounds also admit another infinite family of kinks ${\cal T}_D$, which at large $N$ approach solutions containing free fermion bilinears at $\Delta=D-1$ from below. The kinks ${\cal T}_D$ appear in general dimensions with a $D$-dependent critical $N^*$ below which the kink disappears. We also study relations between the bounds obtained from the bootstrap with $SO(N)$ vectors, $SU(N)$ fundamentals, and $SU(N)\times SU(N)$ bi-fundamentals. We provide a proof for the coincidence between bootstrap bounds with different global symmetries. We show evidence that the proper symmetries of the underlying theories of ${\cal T}_D$ are subgroups of $SO(N)$, and we speculate that the kinks ${\cal T}_D$ relate to the fixed points of gauge theories coupled to fermions., Comment: v1, 59 pages, 10 figures; v2, references added; v3, discussion on the critical flavor number improved, references added, version to appear in JHEP

Three-dimensional conformal field theories (CFTs) with slightly broken higher spin symmetry provide an interesting laboratory to study general properties of CFTs and their roles in the AdS/CFT correspondence. In this work we compute the four-point functions at arbitrary 't Hooft coupling $\lambda$ in the CFTs with slightly broken higher spin symmetry. We use a bootstrap approach based on the approximate higher spin Ward identity. We show that the bootstrap equation is separated into two parts with opposite parity charges, and it leads to a recursion relation for the $\lambda$ expansions of the correlation functions. The $\lambda$ expansions terminate at order $\lambda^2$ and the solutions are exact in $\lambda$. Our work generalizes the approach proposed by Maldacena and Zhiboedov to four-point correlators, and it amounts to an on-shell study for the $3D$ Chern-Simons vector models and the Vasiliev theory in $AdS_4$. Besides, we show that the same results can also be obtained rather simply from bosonization duality of $3D$ Chern-Simons vector models. The odd term at order $O(\lambda)$ in the spinning four-point function relates to the free boson correlator through a Legendre transformation. This provides new evidence on the $3D$ bosonization duality at the spinning four-point function level. We expect this work can be generalized to a complete classification of general four-point functions of single trace currents., Comment: 40 pages, 3 figures; v2: references added

We bootstrap the deconfined quantum critical point (DQCP) and 3D Quantum Electrodynamics (QED$_3$) coupled to $N_f$ flavors of two-component Dirac fermions. We show the lattice and perturbative results on the $SO(5)$ symmetric DQCP are excluded by the bootstrap bounds combined with an irrelevant condition of the lowest singlet scalar. Remarkably, we discover a new family of kinks in the 3D $SO(N)$ vector bootstrap bounds with $N\geqslant6$. We demonstrate bound coincidences between $SU(N_f)$ adjoint and $SO(N_f^2-1)$ vector bootstrap which result from a novel algebraic relation between their crossing equations. By introducing gap assumptions breaking the $SO(N_f^2-1)$ symmetry, the $SU(N_f)$ adjoint bootstrap bounds with large $N_f$ converge to the $1/N_f$ perturbative results of QED$_3$. Our results provide strong evidence that the $SO(5)$ DQCP is not continuous and the critical flavor number of QED$_3$ is slightly above $2$: $N_f^*\in(2,4)$. Bootstrap results near $N_f^*$ are well consistent with the merger and annihilation mechanism for the loss of conformality in QED$_3$., Comment: v1: 6 pages, 3 figures; v2: 4+9 pages, 4+3 figures, quality of figures improved, lattice and perturbative results for DQCP added, more results for QED$_3$ from bootstrap with $SO(N)$-symmetry breaking gaps, discussions improved according to new results

Mammography is currently the primary imaging modality for breast cancer screening and plays an important role in cancer diagnostics. A standard mammographic image acquisition always includes the compression of the breast prior x-ray exposure. The breast is compressed between two plates (the image receptor and the compression paddle) until a nearly uniform breast thickness is obtained. The breast flattening improves diagnostic image quality 1 and reduces the absorbed dose 2. However, this technique can also be a source of discomfort and might deter some women from attending breast screening by mammography 3,4. Therefore, the characterization of the pain perceived during breast compression is of potential interest to compare different compression approaches. The aim of this work is to develop simulation tools enabling the characterization of existing breast compression techniques in terms of patient comfort, dose delivered to the patient and resulting image quality. A 3D biomechanical model of the breast was developed providing physics-based predictions of tissue motion and internal stress and strain intensity. The internal stress and strain intensity are assumed to be directly correlated with the patient discomfort. The resulting compressed breast model is integrated in an image simulation framework to assess both image quality and average glandular dose. We present the results of compression simulations on two breast geometries, under different compression paddles (flex or rigid).

We compute the superconformal partial waves of the four-point correlator $\langle JJJJ\rangle$, in which the external operator $J$ is the superconformal primary of the $4D$ $\mathcal{N}=2$ stress-tensor multiplet $\mathcal{J}$. We develop the superembedding formalism for the superconformal field theories (SCFTs) with extended supersymmetry. In $\mathcal{N}=2$ SCFTs, the three-point functions $\langle \mathcal{J}\mathcal{J}\mathcal{O}\rangle$ with general multiplet $\mathcal{O}$ contain two independent nilpotent superconformal invariants and new superconformal tensor structures, which can be nicely constructed from variables in superembedding space, and the three-point functions can be solved in compact forms. We compute the superconformal partial waves corresponding to the exchange of long multiplets using supershadow approach. The results are consistent with the non-trivial constraints by decomposing the $\mathcal{N}=2$ superconformal blocks into $\mathcal{N}=1$ superconformal blocks. Our results provide the necessary ingredient to study the fascinating $4D$ $\mathcal{N}=2$ SCFTs using conformal bootstrap., Comment: 50 Pages, 0 figure; v2: references added

We show that the scaling dimensions of lowest operators in conformal field theories (CFTs) can be isolated in small and closed regions from single correlator bootstrap. We find the conserved currents play crucial roles in bootstrapping the crossing equation. By imposing a mild gap between the scaling dimensions of the conserved current and its next operator, the scaling dimensions of lowest operators are forced to lie in small isolated regions, i.e., these CFTs can be almost fixed by few lowest operators in certain channels. For CFTs with extended supersymmetry, the single correlator crossing equation involves several conserved or shorted operators and by imposing gaps in these sectors it is possible to isolate different CFTs. Specifically, we bootstrap the isolated regions corresponding to the 3D Ising model, $O(N)$ vector model, $N=1,2$ supersymmetric Ising models by introducing mild gaps in certain sectors with conserved or shorted operators., Comment: 4 figures, 5 pages

We use the conformal bootstrap approach to explore $5D$ CFTs with $O(N)$ global symmetry, which contain $N$ scalars $\phi_i$ transforming as $O(N)$ vector. Specifically, we study multiple four-point correlators of the leading $O(N)$ vector $\phi_i$ and the $O(N)$ singlet $\sigma$. The crossing symmetry of the four-point functions and the unitarity condition provide nontrivial constraints on the scaling dimensions ($\Delta_\phi$, $\Delta_\sigma$) of $\phi_i$ and $\sigma$. With reasonable assumptions on the gaps between scaling dimensions of $\phi_i$ ($\sigma$) and the next $O(N)$ vector (singlet) scalar, we are able to isolate the scaling dimensions $(\Delta_\phi$, $\Delta_\sigma)$ in small islands. In particular, for large $N=500$, the isolated region is highly consistent with the result obtained from large $N$ expansion. We also study the interacting $O(N)$ CFTs for $1\leqslant N\leqslant100$. Isolated regions on $(\Delta_\phi,\Delta_\sigma)$ plane are obtained using conformal bootstrap program with lower order of derivatives $\Lambda$; however, they disappear after increasing $\Lambda$. We think these islands are corresponding to interacting but nonunitary $O(N)$ CFTs. Our results provide a lower bound on the critical value $N_c>100$, below which the interacting $O(N)$ CFTs turn into nonunitary. The critical value is unexpectedly large comparing with previous estimations., Comment: 28 pages, 4 figures

We compute the most general superconformal blocks for scalar operators in $4\mathcal{D}$ $\mathcal{N}=1$ superconformal field theories. Specifically we employ the supershadow formalism to study the four-point correlator $\langle\Phi_1 \Phi_2 \Phi_3 \Phi_4\rangle$, in which the four scalars $\Phi_i$ have arbitrary scaling dimensions and R-charges with the only constraint from R-symmetry invariance of the four-point function. The exchanged operators can have arbitrary R-charges. Our results extend previous studies on $4\mathcal{D}$ $\mathcal{N}=1$ superconmformal blocks to the most general case, which are the essential ingredient for superconformal bootstrap, especially for bootstrapping mixed correlators of scalars with independent scaling dimensions and R-charges., Comment: 33 pages, no figure

We show that a new class of helical phase inflation models can be simply realized in minimal supergravity, wherein the inflaton is the phase component of a complex field and its potential admits a deformed helicoid structure. We find a new unique complex-valued index $\chi$ that characterizes almost the entire region of the $n_s-r$ plane favored by new Planck observations. Continuously varying the index $\chi$, predictions interpolate from quadratic/natural inflation parameterized by a phase/axion decay constant to Starobinsky-like inflation parameterized by the $\alpha$-parameter. We demonstrate that the simple supergravity construction realizing Starobinsky-like inflation can be obtained from a more microscopic model by integrating out heavy fields, and that the flat phase direction for slow-roll inflation is protected by a mildly broken global $U(1)$ symmetry. %, which is mildly broken at the inflation energy scale. We study the geometrical origin of the index $\chi$, and find that it corresponds to a linear constraint relating \kah moduli. We argue that such a linear constraint is a natural result of moduli stabilization in Type \MyRoman{2} orientifold compactifications on Calabi-Yau threefolds with geometric and non-geometric fluxes. Possible choices for the index $\chi$ are discrete points on the complex plane that relate to the distribution of supersymmetric Minkowski vacua on moduli space. More precise observations of the inflationary epoch in the future may provide a better estimation of the index $\chi$. Since $\chi$ is determined by the fluxes and vacuum expectation values of complex structure moduli, such observations would characterize the geometry of the internal space as well., Comment: 26 pages, 4 figures; 4+1 figure, discussion on several energy scales added, references added, to appear in JHEP

The supergravity (SUGRA) theories with exact global $U(1)$ symmetry or shift symmetry in K\"ahler potential provide the natural frameworks for inflation. However, the quadratic inflation is disfavoured by the new results on primordial tensor fluctuations from the Planck Collaboration. To be consistent with the new Planck data, we point out that the explicit symmetry breaking is needed, and study these two SUGRA inflation in details. For the SUGRA inflation with global $U(1)$ symmetry, the symmetry breaking term leads to a trigonometric modulation on inflaton potential. The coefficient of the $U(1)$ symmetry breaking term is of the order $10^{-2}$, which is sufficient large to improve the inflationary predictions while its higher order corrections are negligible. Such models predict sizeable tensor fluctuations and highly agree with the Planck results. In particular, the model with a linear $U(1)$ symmetry breaking term predicts the tensor-to-scalar ratio around $\textbf{r}\sim0.01$ and running spectral index $\alpha_s\sim-0.004$, which comfortably fit with the Planck observations. For the SUGRA inflation with breaking shift symmetry, the inflaton potential is modulated by an exponential factor. The modulated linear and quadratic models are consistent with the Planck observations. In both kinds of models the tensor-to-scalar ratio can be of the order $10^{-2}$, which will be tested by the near future observations., Comment: 20 pages, 10 figures, references added

We study helical phase inflation in supergravity theory in details. The inflation is driven by the phase component of a complex field along helical trajectory. The helicoid structure originates from the monodromy of superpotential with an singularity at origin. We show that such monodromy can be formed by integrating out heavy fields in supersymmetric field theory. The supergravity corrections to the potential provide strong field stabilizations for the scalars except inflaton, therefore the helical phase inflation accomplishes the "monodromy inflation" within supersymmetric field theory. The phase monodromy can be easily generalized for natural inflation, in which the super-Planckian phase decay constant is realized with consistent field stabilization. The phase-axion alignment is fulfilled indirectly in the process of integrating out the heavy fields. Besides, we show that the helical phase inflation can be naturally realized in no-scale supergravity with $SU(2,1)/SU(2)\times U(1)$ symmetry since the no-scale K\"ahler potential provides symmetry factors of phase monodromy directly. We also demonstrate that the helical phase inflation can reduce to the shift symmetry realization of supergravity inflation. The super-Planckian field excursion is accomplished by the phase component, which admits no dangerous polynomial higher order corrections. The helical phase inflation process is free from the UV-sensitivity problem, and it suggests that inflation can be effectively studied in supersymmetric field theory close to the unification scale in Grand Unified Theory and a UV-completed frame is not prerequisite., Comment: 26 pages, 3 figures

We show that the quadratic inflation can be realized by the phase of a complex field with helicoid potential. Remarkably, this helicoid potential can be simply realized in minimal supergravity. The global $U(1)$ symmetry of the K\"ahler potential introduces a flat direction and evades the $\eta$ problem automatically. So such inflation is technically natural. The phase excursion is super-Planckian as required by the Lyth bound, while the norm of the complex field can be suppressed in the sub-Planckian region. This model resolves the ultraviolet sensitive problem of the large field inflation, besides, it also provides a new type of monodromy inflation in supersymmetric field theory with consistent field stabilization., Comment: 4 pages+references, 2 figures, phase monodromy in supersymmetric field theory is added

To obtain natural inflation with large tensor-to-scalar ratio in string framework, we need a special moduli stabilization mechanism which can separate the masses of real and imaginary components of K\"ahler moduli at different scales, and achieve a trans-Planckian axion decay constant from sub-Planckian axion decay constants. In this work, we stabilize the matter fields by F-terms and the real components of K\"ahler moduli by D-terms of two anomalous $U(1)_X\times U(1)_A$ symmetries strongly at high scales, while the corresponding axions remain light due to their independence on the Fayet-Iliopoulos (FI) term in moduli stabilization. The racetrack-type axion superpotential is obtained from gaugino condensations of the hidden gauge symmetries $SU(n)\times SU(m)$ with massive matter fields in the bi-fundamental respresentations. The axion alignment via Kim-Nilles-Pelroso (KNP) mechanism corresponds to an approximate $S_2$ exchange symmetry of two K\"ahler moduli in our model, and a slightly $S_2$ symmetry breaking leads to the natural inflation with super-Planckian decay constant., Comment: 17 pages, no figure

We propose a natural inflation model driven by an imaginary or axionic component of a K\"ahler modulus in string-inspired supergravity. The shift symmetry of the axion is gauged under an anomalous $U(1)_X$ symmetry, which leads to a modulus-dependent Fayet-Iliopoulos (FI) term. The matter fields are stabilized by F-terms, and the real component of the modulus is stabilized by the $U(1)_X$ D-term, while its axion remains light. Therefore, the masses of real and imaginary components of the modulus are separated at different scales. The scalar potential for natural inflation is realized by the superpotential from the non-perturbative effects. The trans-Planckian axion decay constant, which is needed to fit with BICEP2 observations, can be obtained naturally in this model., Comment: 14 pages, no figure, references added, version published in JHEP