Pith. sign in

REVIEW 2 major objections 5 minor 30 references

The one-loop H±→W±Z amplitude in the two-Higgs-doublet model carries a CP-violating charge asymmetry that decomposes into bosonic, fermionic, and interference parts, and in the alignment limit only the fermionic part survives when the addit

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

T0 review · deepseek-v4-flash

2026-08-03 05:41 UTC pith:I72XEAVE

load-bearing objection Competent, genuinely additive one-loop 2HDM calculation; the alignment-limit vanishing claim has an addressable but real analyticity gap. the 2 major comments →

arxiv 2602.01332 v2 pith:I72XEAVE submitted 2026-02-01 hep-ph

CP violation in H^pm to W^pm Z: A physical approach for the 2HDM

classification hep-ph PACS 12.60.Fr11.30.Er
keywords CP violationtwo-Higgs-doublet modelcharged Higgs bosoncharge asymmetryH±→W±Z vertexalignment limitYukawa phasesloop-induced decay
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper tries to prove that the one-loop decay H± → W±Z in the two-Higgs-doublet model carries a CP-violating charge asymmetry that can be decomposed into three independent sources: interference inside the purely bosonic loops, inside the purely fermionic loops, and between bosonic and fermionic loops. It further shows that in the alignment limit, the purely bosonic contribution to the asymmetry disappears as long as the two additional neutral scalars satisfy m2 + m3 > mZ, leaving fermionic loops and their tadpoles as the only possible origin of a surviving asymmetry. This is worth knowing because the entire calculation is expressed in terms of physical masses and couplings, so a measured charge asymmetry would translate directly into constraints on the complex phases of the underlying scalar and Yukawa sector, without needing to fix the scalar potential. The authors argue the mechanism is generic to any CP-violating multi-Higgs-doublet model.

Core claim

The central claim is that the charge asymmetry δ = (Γ(H+→W+Z) − Γ(H−→W−Z)) / (Γ(H+→W+Z) + Γ(H−→W−Z)) arises at one loop from three distinct interference families. In addition to the previously identified interference between bosonic and fermionic amplitudes, the paper identifies CP violation generated entirely within the bosonic loops (through the complex phases of the couplings fi that link each neutral Higgs to the charged Higgs and W) and entirely within the fermionic loops (through relative phases among the Yukawa couplings ρ~). All amplitudes are written directly in terms of the physical couplings ei, fi, qi and ρ, with loop integrals identical for both charge states. In the exact align

What carries the argument

The argument is carried by a reparametrisation of the 2HDM in terms of physical couplings rather than potential parameters: the neutral-scalar couplings to vector boson pairs (ei, with e1²+e2²+e3² = v²), the charged-Higgs–neutral-scalar–W couplings (fi, with phases fixed by the ei), the neutral-scalar–charged-pair couplings (qi), and the generic Yukawa matrices ρ~ that carry independent phases for up-type, down-type and lepton sectors. The charge asymmetry is organised by writing the one-loop amplitude as a sum over diagrams with a common prefactor fi or ρ~, so that complex conjugation of the external charge reverses only these coupling phases while leaving the Passarino–Veltman integrals un

Load-bearing premise

The vanishing of the purely bosonic asymmetry in the alignment limit rests on the premise, stated in section 4.4, that for m2+m3>mZ the only potentially complex loop integral among the surviving bosonic diagrams is C00 of diagram C-1, with all others real; if any other diagram or subleading term acquires an imaginary part, the asymmetry need not vanish.

What would settle it

Compute the full one-loop bosonic amplitude numerically in the exact alignment limit for a scalar spectrum with m2+m3>mZ, without imposing the 'only C-1 is complex' premise; if the imaginary parts of any other triangle or bubble integrals are non-zero, the claim that the asymmetry vanishes is refuted. A more direct experimental falsifier would be a measured non-zero charge asymmetry in a parameter region where fermionic contributions are kinematically or parametrically suppressed (e.g., mH± < mt+mb and |ρ| → 0).

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • If the decomposition is correct, a measured charge asymmetry in H±→W±Z directly constrains combinations of physical phases (phases of fi, phases of ρ, and the relative boson–fermion phase), not just the parameters of the scalar potential.
  • In the alignment limit with m2+m3>mZ, any observed asymmetry is evidence of CP violation in the fermionic couplings or in the fermion-tadpole sector, since the purely bosonic contribution is predicted to vanish.
  • The three sources of asymmetry can in principle be separated by scanning the charged-Higgs mass: fermionic amplitudes become complex immediately above the tb threshold, while bosonic amplitudes require the heavier mHi + mW threshold.
  • The construction extends trivially to any multi-Higgs-doublet model, suggesting the charge asymmetry is a generic signature of extended scalar sectors, not a special feature of the 2HDM.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • A dedicated numerical scan that drops the analyticity premise of section 4.4 could test whether subleading bosonic diagrams develop imaginary parts below the m2+m3>mZ threshold; if they do, the vanishing of the bosonic asymmetry in the alignment limit would be an artifact of the assumption.
  • The same charge-asymmetry variable could be adapted to H±→W±γ, where the photon coupling forbids the scalar-mediated diagrams but shares the WZ gauge structure; a combined measurement of both asymmetries would help isolate the fermionic contribution.
  • One could use the relation between the asymmetry and the CP-odd invariants ImJi to design sum rules linking the asymmetry size to known collider constraints on e2, e3 and qi, turning a future observation into a measurement of the invariants themselves.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

2 major / 5 minor

Summary. The paper computes the one-loop amplitude and charge asymmetry δ for H^± → W^± Z in the general 2HDM, expressing all amplitudes in terms of physical couplings (e_i, f_i, q_i, and the Yukawa matrices ρ̃). The amplitude is split into bosonic and fermionic loop contributions, including fermion tadpoles in the bosonic category. The central claims are: (i) there are CP-violating charge asymmetries from pure-bosonic interference, pure-fermionic interference, and boson–fermion interference; (ii) this qualitatively confirms the Kanemura–Mura result for boson–fermion interference, while also providing new sources; (iii) in the exact alignment limit with no fermionic loops and with m2 + m3 > mZ, the purely bosonic asymmetry vanishes; and (iv) away from alignment, or with fermionic loops, a charge asymmetry can survive. Extensive one-loop amplitudes are tabulated, cross-checked with FormCalc and by hand, and code is made available on GitHub.

Significance. If the claims are correct, this work is a valuable addition to 2HDM phenomenology: it isolates contributions to CP violation in H^± → W^± Z in terms of physical couplings rather than potential parameters, and it identifies new sources of the charge asymmetry beyond those studied by Kanemura and Mura. The paper has notable strengths: the analytic tabulation of loop amplitudes is detailed and machine-checked, the physical-coupling parametrization makes the phase structure transparent, the CPT/Cutkosky discussion gives a useful consistency check, and the GitHub repository provides reproducible analytic results. However, the alignment-limit vanishing claim — which is highlighted in the abstract — rests on an analyticity assertion that is not proven, and one of the analyticity statements in the fermionic section is incorrect as written. These issues are local and fixable, but they directly affect the paper's central claims.

major comments (2)
  1. [Section 4.4, Eq. (4.13)] The alignment-limit vanishing claim is load-bearing for the abstract, but it rests on the unproved assertion that 'it turns out that only diagram C-1 can possibly contain loop integrals with an imaginary part' and that its C00 is real for m2+m3>mZ. The surviving amplitudes in Table 4.1 also involve the combinations C1, C11, C12, C22 for C-1 and C-2, as well as B0 and B1 for the bubbles. No analyticity argument is given for these coefficients, and the numerical illustrations in Fig. 8 include fermionic tadpoles; they do not plot the no-fermion asymmetry in the exact alignment limit. Please provide a complete proof, or a numerical scan of the imaginary parts of all surviving Passarino–Veltman coefficients, together with a no-fermion asymmetry plot in the alignment limit. Without this, Eq. (4.13) and the abstract conclusion remain conditional.
  2. [Section 5.1, just before Eqs. (5.3)] The statement 'Since the loop integrals involved in the H amplitudes are all real' is not correct. For the tb triangles in Table 5.4, the coefficients C0 and C1 become complex when m_H± > m_t + m_b, which is precisely the region where the fermionic F and G amplitudes develop their imaginary parts. The conclusion that the H amplitude does not contribute to the charge asymmetry can instead be justified by the fact that the H+ and H- amplitudes are proportional to ρ^* Z and ρ Z with a common loop factor Z for diagonal ρ, so |H+|^2 = |H-|^2. The erroneous analyticity claim should be removed and replaced by this argument.
minor comments (5)
  1. [Notation, Eq. (4.10) and Section 2.3] The symbol ρ is used both for the kinematic tadpole factor in Eq. (4.10) and (in the form ρ̃) for the Yukawa matrices. This is confusing in Tables 4.3 and 4.4. Please rename one of them.
  2. [Figure 8 caption] The phrase 'the thin curves below do not' is ambiguous. Please specify exactly which contributions are omitted in the thin curves.
  3. [Section 4.5] There is a typo: 'threshol' should be 'threshold'.
  4. [Abstract / Comparison with Ref. [4]] The paper says it 'qualitatively confirms' Kanemura and Mura, but no quantitative comparison or parameter-overlap plot with Ref. [4] is provided. Please add one, so that the confirmation claim can be checked and the differences identified in Sections 4.6 and 5.3 are put in context.
  5. [Table 4.1] The remark that 'i' in the K_d column denotes the imaginary unit should be made in the table caption, not only in a footnote in the text.

Circularity Check

0 steps flagged

No significant circularity: the one-loop H±→W±Z asymmetry is derived from explicit loop integrals; self-citations are notational, and the alignment-limit claim rests on an unproved analyticity assertion, which is a correctness gap, not circularity.

full rationale

The paper's central quantity, the charge asymmetry δ, is computed from one-loop amplitudes whose dependence on the couplings e_i, f_i, q_i and ρ~ is set by Feynman rules and evaluated with Passarino–Veltman integrals (Tables 4.1–5.8). No parameter appearing in the final asymmetry is fitted to the asymmetry itself: the physical couplings are either algebraic definitions from prior work (Refs. [12,13]) or illustrative benchmark values (e.g., |ρ~|=0.1, imported from Ref. [4] and used only for numerical examples). The claimed new results—bosonic self-interference and fermionic self-interference contributions—follow from the structure of Eqs. (4.7), (5.4) and (3.7), not from any input equivalent to the output. The self-citations to Refs. [12,13,31] provide notation, Feynman rules, and independent CP-conditions; they do not supply the asymmetry result, and the custodial-symmetry equivalence is credited to Ref. [5] and Haber–O'Neil [14], which are not the present authors. The only notable weakness is in Section 4.4, where the vanishing of the purely bosonic asymmetry in the alignment limit for m2+m3>mZ rests on the asserted statement that 'only diagram C-1 can possibly contain loop integrals with an imaginary part' and that the surviving C00, B0, B1 are then real. This is an analyticity claim that is not fully demonstrated, but it is an unproved premise about loop-integral reality, not an input that has been renamed as a prediction. A failure of that premise would weaken the claim, but would not make the derivation circular. The paper is self-contained against the independent Kanemura–Mura result [4] and explicitly identifies differences, so no fitted-input-called-prediction or self-definitional circularity is present.

Axiom & Free-Parameter Ledger

4 free parameters · 6 axioms · 0 invented entities

The calculation is built on standard QFT/2HDM assumptions. No new particles or forces are introduced. The numerical inputs are illustrative benchmarks, not fitted parameters. The most fragile premise is the analyticity/reality of loop integrals in the alignment limit.

free parameters (4)
  • Neutral scalar masses m_i = 125, 200, 250 GeV; 125, 200, 400 GeV
    Benchmark inputs for Figures 8–18; the central analytic results are not fitted to these masses.
  • Gauge/Higgs couplings e_i = e1 = 0.95v or 0.99v; e2,e3 split with r=0.5 or scanned
    Chosen by hand to represent near-alignment; the symmetry behavior depends on these but they are not fitted to data.
  • Trilinear couplings q_i = {100,150,10} GeV and {100,150,175} GeV
    Benchmark choices used for illustrations, selected to satisfy boundedness-from-below constraints.
  • Yukawa magnitudes and phases = |ρ̃|=0.1 or 0.01; ζ_t=0,π/4,π/2; ζ_b=ζ_τ=0; θ_BF=0,π/4
    Values adopted from Ref. [4] and flavor-fit literature; used only for numerical illustrations.
axioms (6)
  • domain assumption The 2HDM is the correct effective framework, with three neutral scalars, one charged pair, and the most general Yukawa couplings.
    The entire calculation is framed in the 2HDM (Section 1); the result is only relevant if nature realizes such an extended Higgs sector.
  • domain assumption The one-loop amplitude is the leading contribution, since H±W∓Z vanishes at tree level in SU(2)×U(1) with doublets.
    Reference [1] and Section 1 establish this known result.
  • standard math Passarino–Veltman reduction and dimensional reduction provide a consistent finite one-loop result.
    FormCalc and dimensional reduction are used in Section 4.3.1.
  • domain assumption For on-shell W/Z decaying to light fermions, terms proportional to p2^α or p3^β can be dropped.
    Section 3 restricts the Lorentz structure; this is what reduces the observable to a charge asymmetry.
  • domain assumption A loop integral develops an imaginary part only when √s_i > m_i + m_{i+1}; and for m2+m3>mZ the relevant C00/B0/B1 integrals in the alignment limit are real.
    Sections 2.5 and 4.4; this is the load-bearing premise for the alignment-limit vanishing claim.
  • standard math CPT requires total widths to be equal, so the asymmetry in H±→W±Z is balanced by other channels listed in Eq. (1.7).
    Section 1 and the summary use CPT and Cutkosky rules for consistency.

pith-pipeline@v1.3.0-alltime-deepseek · 30655 in / 12706 out tokens · 127229 ms · 2026-08-03T05:41:26.558379+00:00 · methodology

0 comments
read the original abstract

We investigate CP violation in the process H^\pm \to W^\pm Z within the framework of the two-Higgs-Doublet model (2HDM). Amplitudes are expressed transparently in terms of physical couplings. Our analysis qualitatively confirms recent results by Kanemura and Mura, for the interference between one-loop bosonic and fermionic amplitudes. Furthermore, we identify additional sources of CP violation arising from internal interference among the bosonic and also among the fermionic loop amplitudes. In the alignment limit, the asymmetry would vanish in the absence of fermionic loop contributions when the sum of the masses of the additional neutral scalars is higher than the mass of the Z. Our results allow for generic Yukawa couplings, and cover both explicit and spontaneous CP violation. For inclusive W^\pm and Z decays, the CP violation will only manifest itself as a charge asymmetry.

Figures

Figures reproduced from arXiv: 2602.01332 by Margarida Nesbitt Rebelo, Odd Magne Ogreid, Per Osland, Wafaa Khater.

Figure 1
Figure 1. Figure 1: Left: circles with constant e 2 2 + e 2 3 ; centre and right: charged-Higgs couplings f2 and f3. The phases of the fi are unphysical, the phases of f2 and f3 are shown with respect to that of f1, which is taken real and positive. Three degrees of alignment, κ, are considered (see left panel): outer, middle and inner e2-e3 circles refer to e1 = 0.90v, e1 = 0.95v, and e1 = 0.99v, respectively. The correspond… view at source ↗
Figure 2
Figure 2. Figure 2: Boson triangle diagrams involving the qi (or qj ) couplings. m1 m2 m3 q+k1 q q+k2 Z H+ Hi H+ W− Z p1 p2, α p3, β 1 fi ei (C-3) m1 m2 m3 q+k1 q q+k2 Hi W+, G+ H Z + W− Z p1 p2, α p3, β fi 1 ei (C-4) m1 m2 m3 q+k1 q q+k2 Hi W+, G+ Hj H+ W− Z p1 p2, α p3, β fi ej ek (C-5) m1 m2 m3 q+k1 q q+k2 Hi G+ H G0 + W− Z p1 p2, α p3, β fi 1 ei (C-6) m1 m2 m3 q+k1 q q+k2 W− Hi W− H , G− + W− Z p1 p2, α p3, β fi ei 1 (C-7… view at source ↗
Figure 3
Figure 3. Figure 3: Other boson triangle diagrams. Some represent two cases, with a [PITH_FULL_IMAGE:figures/full_fig_p013_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Boson loops—parts (B-1)–(B-4). 14 [PITH_FULL_IMAGE:figures/full_fig_p015_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Boson loops—parts (B-5)–(B-10). m1 m2 q+k1 q p1 p2, α p3, β fi H+ qi W− Z H+ H+ Hi (B-a) m1 m2 q+k1 q p1 p2, α p3, β ei H+ fi W− Z H+ G+ Hi (B-b) m1 m2 q+k1 q p1 p2, α p3, β ei H+ fi W− Z H+ W+ Hi (B-c) [PITH_FULL_IMAGE:figures/full_fig_p016_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Boson loops—parts (B-a) to (B-c). Bubble diagrams on the W line, like those shown in [PITH_FULL_IMAGE:figures/full_fig_p016_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Boson/fermion tadpole diagrams. Left: parts (T-1)–(T-4). Centre and right: [PITH_FULL_IMAGE:figures/full_fig_p018_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Decay rates in the alignment limit at which [PITH_FULL_IMAGE:figures/full_fig_p019_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Decay rate (left) and asymmetry (right) for some near-alignment cases. The [PITH_FULL_IMAGE:figures/full_fig_p021_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Similar to the lower panels of Fig. 9, where the parameter [PITH_FULL_IMAGE:figures/full_fig_p023_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Left and centre: Complex values of f2 and f3 for 40 random values of r ∈ {0, 1}. Right: Corresponding asymmetries vs r, sorted into 4 sets, “pp” (e2 > 0, e3 > 0), “pm” (e2 > 0, e3 < 0), “mp” (e2 < 0, e3 > 0), “mm” (e2 < 0, e3 < 0). The mass of H± is here 350 GeV. 5 Fermion contributions The boson loop integrations involved up to two momentum-of-integration factors (q α and/or q β ) in the numerator. This … view at source ↗
Figure 12
Figure 12. Figure 12: Quark loops. The arrows indicate momentum flow, not charge flow. [PITH_FULL_IMAGE:figures/full_fig_p024_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Decay rate and asymmetry resulting from fermion loops only ( [PITH_FULL_IMAGE:figures/full_fig_p029_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: Complex-plane vectors representing different contributions to the fermion am [PITH_FULL_IMAGE:figures/full_fig_p029_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: Decay rate and asymmetry resulting from bosonic and fermionic loops ( [PITH_FULL_IMAGE:figures/full_fig_p031_15.png] view at source ↗
Figure 16
Figure 16. Figure 16: Similar to Fig. 15, for [PITH_FULL_IMAGE:figures/full_fig_p032_16.png] view at source ↗
Figure 17
Figure 17. Figure 17: Similar to Fig. 15, for [PITH_FULL_IMAGE:figures/full_fig_p033_17.png] view at source ↗
Figure 18
Figure 18. Figure 18: Similar to Fig. 17, for [PITH_FULL_IMAGE:figures/full_fig_p034_18.png] view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

30 extracted references · 21 linked inside Pith

  1. [1]

    J. A. Grifols and A. Mendez,TheW ZH ± Coupling inSU(2)×U(1)Gauge Models,Phys. Rev. D22(1980) 1725

  2. [2]

    Mendez and A

    A. Mendez and A. Pomarol,One loop inducedH +W +Zvertex in the minimal supersymmetry model,Nucl. Phys. B349(1991) 369–380

  3. [3]

    Mendez and A

    A. Mendez and A. Pomarol,Signals of CP violation in the Higgs sector,Phys. Lett. B272(1991) 313–318

  4. [4]

    Kanemura and Y

    S. Kanemura and Y. Mura,Loop induced H ±W∓Z vertices in the general two Higgs doublet model with CP violation,JHEP10(2024) 041, [2408.06863]

  5. [5]

    Grzadkowski, M

    B. Grzadkowski, M. Maniatis and J. Wudka,The bilinear formalism and the custodial symmetry in the two-Higgs-doublet model,JHEP11(2011) 030, [1011.5228]

  6. [6]

    J. F. Gunion and H. E. Haber,Conditions for CP-violation in the general two-Higgs-doublet model,Phys. Rev. D72(2005) 095002, [hep-ph/0506227]

  7. [7]

    G. C. Branco, M. N. Rebelo and J. I. Silva-Marcos,CP-odd invariants in models with several Higgs doublets,Phys. Lett. B614(2005) 187–194, [hep-ph/0502118]

  8. [8]

    Lavoura and J

    L. Lavoura and J. P. Silva,Fundamental CP violating quantities in a SU(2)×U(1) model with many Higgs doublets,Phys. Rev. D50(1994) 4619–4624, [hep-ph/9404276]

  9. [9]

    F. J. Botella and J. P. Silva,Jarlskog - like invariants for theories with scalars and fermions,Phys. Rev. D51(1995) 3870–3875, [hep-ph/9411288]

  10. [10]

    R. F. Streater and A. S. Wightman,PCT, spin and statistics, and all that. Benjamin/Cummings, Reading, Mass., 1964

  11. [11]

    R. E. Cutkosky,Singularities and discontinuities of Feynman amplitudes,J. Math. Phys.1(1960) 429–433

  12. [12]

    Grzadkowski, O

    B. Grzadkowski, O. M. Ogreid and P. Osland,Measuring CP violation in Two-Higgs-Doublet models in light of the LHC Higgs data,JHEP11(2014) 084, [1409.7265]

  13. [13]

    Grzadkowski, H

    B. Grzadkowski, H. E. Haber, O. M. Ogreid and P. Osland,Heavy Higgs boson decays in the alignment limit of the 2HDM,JHEP12(2018) 056, [1808.01472]

  14. [14]

    H. E. Haber and D. O’Neil,Basis-independent methods for the two-Higgs-doublet model III: The CP-conserving limit, custodial symmetry, and the oblique parameters S, T, U,Phys. Rev. D83(2011) 055017, [1011.6188]

  15. [15]

    M. Aoki, S. Kanemura, K. Tsumura and K. Yagyu,Models of Yukawa interaction in the two Higgs doublet model, and their collider phenomenology,Phys. Rev. D80 (2009) 015017, [0902.4665]. 37

  16. [16]

    G. C. Branco, P. M. Ferreira, L. Lavoura, M. N. Rebelo, M. Sher and J. P. Silva, Theory and phenomenology of two-Higgs-doublet models,Phys. Rept.516(2012) 1–102, [1106.0034]

  17. [17]

    Hahn and M

    T. Hahn and M. Perez-Victoria,Automatized one loop calculations in four-dimensions and D-dimensions,Comput.Phys.Commun.118(1999) 153–165, [hep-ph/9807565]

  18. [18]

    F. J. Botella, G. C. Branco, M. Nebot and M. N. Rebelo,New physics and evidence for a complex CKM,Nucl. Phys. B725(2005) 155–172, [hep-ph/0502133]. [19]CKMfitter Groupcollaboration, J. Charles, A. Hocker, H. Lacker, S. Laplace, F. R. Le Diberder, J. Malcles et al.,CP violation and the CKM matrix: Assessing the impact of the asymmetricBfactories,Eur. Phys....

  19. [20]

    Tree-level metastability bounds of the 2HDM in terms of physical parameters

    O. M. Ogreid, “Tree-level metastability bounds of the 2HDM in terms of physical parameters.” Talk at MultiHiggs workshop, Lisboa September 2024, http://cftp.tecnico.ulisboa.pt/~2hdmwork/index.php?n=Main.HomePage, 2024

  20. [21]

    I. P. Ivanov and J. P. Silva,Tree-level metastability bounds for the most general two Higgs doublet model,Phys. Rev. D92(2015) 055017, [1507.05100]

  21. [22]

    Passarino and M

    G. Passarino and M. J. G. Veltman,One Loop Corrections for e+ e- Annihilation Into mu+ mu- in the Weinberg Model,Nucl. Phys. B160(1979) 151–207

  22. [23]

    G. J. van Oldenborgh and J. A. M. Vermaseren,New Algorithms for One Loop Integrals,Z. Phys. C46(1990) 425–438

  23. [24]

    Hahn,Generating Feynman diagrams and amplitudes with FeynArts 3,Comput

    T. Hahn,Generating Feynman diagrams and amplitudes with FeynArts 3,Comput. Phys. Commun.140(2001) 418–431, [hep-ph/0012260]

  24. [25]

    Kublbeck, M

    J. Kublbeck, M. Bohm and A. Denner,Feyn Arts: Computer Algebraic Generation of Feynman Graphs and Amplitudes,Comput. Phys. Commun.60(1990) 165–180. [26]Particle Data Groupcollaboration, S. Navas et al.,Review of particle physics, Phys. Rev. D110(2024) 030001

  25. [27]

    Athron, A

    P. Athron, A. Crivellin, T. E. Gonzalo, S. Iguro and C. Sierra,Global fit to the 2HDM with generic sources of flavour violation using GAMBIT,JHEP11(2024) 133, [2410.10493]

  26. [28]

    Crivellin, A

    A. Crivellin, A. Kokulu and C. Greub,Flavor-phenomenology of two-Higgs-doublet models with generic Yukawa structure,Phys. Rev. D87(2013) 094031, [1303.5877]

  27. [29]

    Pomarol and R

    A. Pomarol and R. Vega,Constraints on CP violation in the Higgs sector from the rho parameter,Nucl. Phys. B413(1994) 3–15, [hep-ph/9305272]

  28. [30]

    Pilaftsis,On the Classification of Accidental Symmetries of the Two Higgs Doublet Model Potential,Phys

    A. Pilaftsis,On the Classification of Accidental Symmetries of the Two Higgs Doublet Model Potential,Phys. Lett. B706(2012) 465–469, [1109.3787]. 38

  29. [31]

    Grzadkowski, O

    B. Grzadkowski, O. M. Ogreid and P. Osland,Spontaneous CP violation in the 2HDM: physical conditions and the alignment limit,Phys. Rev. D94(2016) 115002, [1609.04764]

  30. [32]

    P. M. Ferreira, B. Grzadkowski, O. M. Ogreid and P. Osland,Symmetries of the 2HDM: an invariant formulation and consequences,JHEP02(2021) 196, [2010.13698]. 39