Closed MJKirk closed 2 years ago
Did you switch off 2-loop contributions in v2?
I assume you mean 1-loop (unless wilson is secretly much better than announced!).
I believe the default is tree matching, but I have explicitly called wilson.Wilson.set_default_option('smeft_matching_order', 0)
at the beginning to make sure and the results are the same.
@MJKirk
I believe @DavidMStraub meant 1-loop SMEFT-WET matching only.
What SMEFT basis did you use for matching CVLL_ucuc?
I think this matching (12 psi^2phi^2 D type SMEFT WCs in the Warsaw-down basis to sdsd WCs in WET) requires two insertions of SMEFT operators which I do not think is implemented in wilson. But the non-zero result that you are getting might be due to the running of SMEFT Wilson coefficients from 100 GeV (your input scale) to 91.1876 GeV (default matching scale in the wilson). To see this try setting scale= 91.1876 GeV instead of 100 GeV in your initial condition.
Still, the origin of the difference between the two versions should be investigated.
What SMEFT basis did you use for matching CVLL_ucuc?
I am matching from SMEFT Warsaw
to WET flavio
But the non-zero result that you are getting might be due to the running of SMEFT Wilson coefficients from 100 GeV (your input scale) to 91.1876 GeV (default matching scale in the wilson). To see this try setting scale= 91.1876 GeV instead of 100 GeV in your initial condition.
Ah good point, for some reason (that reason being I didn't check carefully), I assumed just setting the scales equal would be enough to prevent running.
Setting the scale to MZ to avoid this, the ucuc
coefficients are exactly 0 in v2
, while phiq3_12
gives a non-zero but tiny (10^-23) value for sdsd
.
I think this matching (12 psi^2phi^2 D type SMEFT WCs in the Warsaw-down basis to sdsd WCs in WET) requires two insertions of SMEFT operators which I do not think is implemented in wilson. Still, the origin of the difference between the two versions should be investigated.
Reading both papers, at the end of Sec 4 of 1709.04486, they comment about these sorts of double insertions and say that "the last term, from the product of two dimension-six corrections to the gauge coupling should formally be dropped as it is of higher order in the power counting". However, looking at the old code:
https://github.com/wilson-eft/wilson/blob/1686d08487881fe692f9d3c69b90dbcc16ba44c8/wilson/match/_smeft_old.py#L92
it seems these were included in v1
which would explain the result.
Whereas in the new code
https://github.com/wilson-eft/wilson/blob/1686d08487881fe692f9d3c69b90dbcc16ba44c8/wilson/match/smeft_tree.py#L98
it seems these quadratic terms are correctly dropped.
So I think that clears up my narrow problem at least - these terms should not have been included originally, and as of v2
are not.
And I see now that in the test David wrote when he updated the matching: https://github.com/wilson-eft/wilson/blob/1686d08487881fe692f9d3c69b90dbcc16ba44c8/wilson/match/test_smeft_tree_new.py#L46-L54 the different treatment of quadratic terms was allowed for, so I think it is safe to say this is resolved.
@MJKirk thanks for such an in-depth check! I would argue that, in general, the quadratic terms must be dropped. The power counting is not the right argument, though. Much more importantly, double insertions of dim-6 SMEFT operators are in general not finite, and you will need dim-8 counter terms. While in specific cases this might be avoidable, e.g., when you cannot write down a dim-8 term for some reason, in general it is not. Hence, you cannot have a consistent treatment of double insertions of dim-6 operators without also taking dim-8 operators into account, regardless of the power-counting argument.
I've found a change in the matching results wilson gives between
v1.8
andv2.0
(I think the v1 results are correct, and v2 is wrong). Inwilson v1.8
: Setting the initial conditionsphiq3_12 = 1e-8
@ 100 GeV, then matching to the WET (in the flavio basis, at the same scale) givesCVLL_ucuc = -2.47e-12
. Inwilson v2.0
: The same thing gives a result ofCVLL_ucuc = -1.77e-19
, so 10^7 smaller! The phiq3 operator (after symmetry breaking) affects the W and Z couplings, so phiq3_12 gives a LH Z-u-c vertex, and then the matching is relatively straightforward. Doing it by hand I agree with the larger result fromv1.8
.I see similar differences in the matching of
phiu_12
(which gives a RH Z-u-c, 10^11 smaller in v2), andphiq3_11
(gives LH Z-u-c after CKM rotation, smaller by 10^7). The flavour diagonal matching (phiq3_11 -> CVLL_uuuu) is the basically same between versions though.Similar effects in the down sector:
I haven't had a change to check the Mathematica notebook from the new source paper to see what happens there, or whether it is just a bug in the conversion to Python.