A Deep Dive into Baryon Asymmetry -- the C2HDM

This paper presents a new implementation of baryon asymmetry calculations in the BSMPT code based on a generalized WKB ansatz for transport equations, which is validated and applied to the CP-violating 2-Higgs-Doublet Model to analyze key dependencies, uncertainties, and the interplay with gravitational wave signals.

The Gist

The Big Mystery: Why is there more stuff than "anti-stuff"?

Imagine the universe as a giant party. According to the laws of physics, when the party started (the Big Bang), it should have created equal amounts of "matter" (the good guys) and "antimatter" (the bad guys). If they met, they would have annihilated each other, leaving behind only pure energy and no one to tell the story.

But here we are. We exist. There is a tiny, tiny excess of matter over antimatter. Scientists call this the Baryon Asymmetry. The paper asks: How did this tiny imbalance happen?

The Scenario: A Cosmic Bubble Party

The authors propose a scenario called Electroweak Baryogenesis. Imagine the early universe as a pot of boiling water. As it cools down, bubbles of a new state of matter start forming inside the water (like steam bubbles in boiling water).

1. The Bubble Wall: As these bubbles expand, they have a "wall" moving through the hot plasma.
2. The Reflection: When particles hit this moving wall, they bounce off. Because of a subtle rule-breaking in physics called CP Violation (think of it as a slight bias in how the universe treats left-handed vs. right-handed particles), the wall reflects "good guys" and "bad guys" differently.
3. The Result: This creates a pile-up of particles just outside the bubble wall.
4. The Capture: Inside the bubble, a "cleaning crew" (called sphalerons) usually wipes out any imbalance. But if the bubble forms fast enough and the wall is strong enough, this cleaning crew gets suppressed inside the bubble, trapping the imbalance. The universe ends up with a little extra matter.

What This Paper Actually Did

The authors didn't discover a new particle; they built a better calculator to figure out exactly how much extra matter is produced in this scenario. They updated a software tool called BSMPT (which stands for "Beyond Standard Model Phase Transitions").

Think of their work as upgrading a weather simulation. Previous versions might have guessed the wind speed or the shape of the storm. This new version tries to calculate those things with much higher precision.

The Two Main Upgrades

The paper highlights two major improvements to their calculator:

1. The "Moment" Expansion (Counting the Details)
To predict how particles move, the authors use a mathematical trick called a "moment expansion."

• The Analogy: Imagine trying to describe the traffic on a highway.
  • Low precision: You just say, "There are 1,000 cars."
  • Medium precision: You say, "There are 1,000 cars, and 60% are going 60mph."
  • High precision: You track the speed, direction, and acceleration of every single car in every lane.
• The Paper's Claim: They upgraded their code to track up to 50 different "moments" (layers of detail) instead of just a few. They found that while adding more details makes the math harder, it changes the answer. Surprisingly, the answer keeps changing even after 50 layers, suggesting we might need even more detail to get the "true" answer.

2. The Shape of the Bubble Wall (The Kink vs. The Real Thing)
The bubble wall isn't a sharp line; it's a transition zone.

• The Old Way (Kink Profile): Scientists used to assume the wall looked like a perfect, smooth "S" curve (a mathematical kink). It's a nice, simple shape to draw.
• The New Way (Field Profile): The authors now solve the actual equations of motion to see what the wall really looks like.
• The Discovery: The real wall is often "fatter" and more complex than the simple "S" curve. This shape matters because it changes how particles bounce off it. They found that using the simple "S" curve often overestimates how much matter is created.

The "C2HDM" Model

They tested their new calculator using a specific theory called the CP-violating 2-Higgs-Doublet Model (C2HDM).

• The Analogy: The Standard Model of physics is like a car with one engine. The C2HDM is like a car with two engines (two Higgs fields).
• The Goal: They wanted to see if having two engines creates enough "CP violation" (bias) to explain why we have matter.

Key Findings & Warnings

The paper is very honest about the uncertainties in their calculation. Here is what they found:

• The "Goldilocks" Problem: To get a stable, reliable answer, the bubble wall needs to be very wide and the universe needs to be expanding at a specific speed. If the wall is too thin or the expansion is too slow, the math gets messy and the answer jumps around wildly.
• The Trade-off: The conditions that make the math stable (wide walls, fast expansion) actually result in less matter being created. The conditions that create more matter (thin walls, slow expansion) make the math unstable and unreliable.
• The CP Violation: They confirmed that the more "bias" (CP violation) you put into the model, the more matter is created. This is a crucial guide for future model builders: if you want to explain our universe, your theory needs a lot of this specific type of bias.
• Gravitational Waves: They checked if these bubble collisions would create ripples in space-time (gravitational waves) that the LISA telescope could detect.
  • Type I Model: Some scenarios produce detectable waves, but they don't produce enough matter to explain our universe.
  • Type II Model: The rules are too strict; they produce neither enough matter nor detectable waves.

The Bottom Line

The authors have built a more powerful, more consistent engine to simulate the birth of matter in the universe. They found that:

1. We need to look at the problem with extreme mathematical detail (many "moments") to get a reliable answer.
2. The shape of the bubble wall is more complex than we thought, and using simple shapes gives the wrong answer.
3. There is a tension: the scenarios that are mathematically safe to calculate often predict too little matter, while the scenarios that predict enough matter are mathematically risky to calculate.

They conclude that while their tool is a big step forward, we still need to refine our math to be sure exactly how the universe got its extra matter.

Technical Summary

Technical Summary: A Deep Dive into Baryon Asymmetry – the C2HDM

Problem Statement
The origin of the baryon asymmetry of the Universe (BAU), quantified as $\eta_{obs} \approx 8.69 \times 10^{-11}$, remains a central open question in particle physics. Electroweak baryogenesis (EWBG) offers a dynamical explanation satisfying the Sakharov conditions, specifically requiring a strong first-order phase transition, baryon number violation, and CP violation. However, calculating the generated BAU involves complex non-equilibrium transport processes across a bubble wall. Previous implementations often relied on simplified approximations, such as limited moment expansions of the transport equations, fixed "kink" profiles for vacuum expectation values (VEVs), and arbitrary truncation schemes. This paper addresses the need for a more rigorous, self-consistent, and flexible computational framework to evaluate the BAU in extended Higgs sectors, specifically the CP-violating Two-Higgs-Doublet Model (C2HDM).

Methodology
The authors present a new implementation within the code BSMPT (specifically version 4), which performs the entire chain from a zero-temperature particle physics model to the computation of the BAU at the electroweak phase transition. The core methodological advancements include:

1. WKB Ansatz and Moment Expansion: The transport equations are derived using the WKB ansatz, generalizing the system to an arbitrary number of moments ($n$). The authors implement up to $n=50$ moment equations to solve for the chemical potentials and out-of-equilibrium perturbations of particle species (top quarks, bottom quarks, Higgs bosons).
2. Truncation Schemes: Two distinct schemes for closing the system of moment equations are implemented:
   • Constant Truncation: Relates the highest moment derivative to the previous one via a constant factor $R$ (e.g., $R = -v_w$).
   • Variance Truncation: A more sophisticated scheme where the $n$-th variance is set to zero, relating the highest moment to all previous moments.
3. VEV Profiles: Beyond the standard "kink" profile ansatz, the code now derives VEV profiles by solving the equations of motion (EOMs) for the scalar fields, including a friction term dependent on wall velocity. This allows for a more realistic modeling of the bubble wall structure and the associated complex mass phases.
4. Collision Terms: The implementation includes a detailed collision network involving top and bottom quarks and the Higgs boson. A key update is the authors' own recalculation of the top-Yukawa rate ($\hat{\Gamma}_y$) at non-zero temperature using finite temperature field theory, yielding $\hat{\Gamma}_y \approx 6 \times 10^{-3} T$.
5. Model Framework: The analysis focuses on the C2HDM, which introduces a second Higgs doublet with a softly broken $Z_2$ symmetry, allowing for spontaneous CP violation. The study considers both Type-I and Type-II Yukawa couplings.

Key Results
The paper validates the new implementation using a benchmark model and then performs a comprehensive numerical analysis within the C2HDM:

• Convergence and Uncertainties: The study reveals that the computed BAU is highly sensitive to the number of moment equations ($n$) and the truncation scheme. While increasing $n$ generally reduces the magnitude of the BAU, convergence to a stable value is not guaranteed for $n \le 50$ in all parameter regions.
• WKB Validity: The validity of the WKB approximation, which assumes a large wall width relative to the thermal scale ($L_w T \gg 1$), is scrutinized. The authors find that stable, convergent results typically require $L_w T_n \gtrsim 50$, a significantly stricter condition than the $L_w T_n \gtrsim 2$ often cited in literature. In the C2HDM, most viable points yield $L_w T_p < 7$, indicating a large inherent uncertainty in the BAU calculation for these models.
• Profile Dependence: Comparing the "kink" profile with the "field" profile (derived from EOMs) shows that the field solution often results in a larger effective wall width, leading to a smaller calculated BAU. Furthermore, the field solution can exhibit complex behaviors in the CP-violating phase ($\omega_{CP}$) that differ significantly from the simple kink ansatz.
• CP Violation Impact: A strong positive correlation is observed between the magnitude of the CP-violating parameters (specifically $\text{Im}(m_{12}^2)$ and the induced VEV phase) and the generated BAU.
• Gravitational Waves (GW): The study investigates the correlation between the BAU and the gravitational wave signal detectable by LISA. For C2HDM Type-I, parameter points exist that satisfy the observed BAU but produce GW signals below LISA's detection threshold, or vice versa (detectable GWs but insufficient BAU). For Type-II, experimental constraints at zero temperature severely limit the parameter space, resulting in neither sufficient BAU nor detectable GW signals.
• Benchmark Points: Four specific benchmark points (BP1–BP4) are analyzed to illustrate the interplay between wall velocity, transition temperature, and profile choices. Notably, the choice of transition temperature (critical vs. percolation) and the behavior of the Euclidean action $S_3/T$ significantly impact the final BAU values.

Significance and Claims
The authors claim that this work provides the most comprehensive and self-consistent code (BSMPTv4) currently available for computing the BAU from first principles within an extended Higgs sector. By integrating the full chain from model definition to cosmological observables, the code ensures algorithmic consistency and allows for robust uncertainty estimates.

The paper emphasizes that while the implementation is applicable to any extended Higgs sector, the current results highlight significant theoretical uncertainties, particularly regarding the convergence of the moment expansion and the validity of the WKB approximation in regions of interest. The authors conclude that future improvements in the computation of the BAU are crucial for reliably deducing model parameters from cosmological observations. The study serves as a guide for model building, suggesting that models with large CP violation are favored, but also cautioning that current approximations may overestimate the BAU in regions where the phase transition is strong but the wall width is small. The uncertainty analysis presented is posited as the necessary basis for any future attempt to constrain model parameters using cosmological data.