Source author record

Kyohei Mukaida

Kyohei Mukaida appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

33works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

33 published item(s)

preprint2024arXiv

Quantum Corrections to Higgs Inflation in Einstein-Cartan Gravity

This paper studies the quantum corrections to the Higgs inflation model in the context of the Einstein-Cartan (E-C) gravity in the large-$ N $ limit with $N$ being the number of real scalar components in Higgs. Recently, it is realized that the Higgs inflation in the E-C formalism smoothly connects those in the metric and the Palatini formalisms in the presence of a non-minimal coupling between the Higgs fields and the Nieh-Yan term. This motivates us to investigate the quantum corrections to the E-C Higgs inflation and to clarify how the Ricci curvature squared $ R^2 $ induced by the quantum corrections succeeds in Ultraviolet (UV)-extending the Higgs inflation in metric formalism while it fails in the Palatini case. We show that a generalized $ R^2 $-term required for the renormalization in the E-C formalism induces a new scalar degree of freedom (DoF), the scalaron, which gradually decouples with the system due to its increasing mass as approaching the Palatini limit. The presence of the scalaron extends the UV cutoff at vacuum of the original model except for the parameter space close to the Palatini limit. This UV-extension is expected to solve the strong coupling problem that may exist during (p)reheating in the absence of the scalaron.

preprint2022arXiv

Charge transfer between rotating complex scalar fields

We consider the transfer of a U(1) charge density between Bose-Einstein condensates of complex scalar fields coupled to a thermal bath, focusing on the case of a homogeneous Affleck-Dine field transmitting the charge stored in its angular motion to an axion field. We demonstrate that in the absence of additional symmetries this charge transfer, aided by cosmic expansion as well as the thermal effective potential of the Affleck-Dine field, can be very efficient. The charge redistribution between the scalar fields becomes possible if the interactions with the thermal bath break the original U(1) x U(1) symmetry down to a single U(1) symmetry; the charge distribution between the two fields is then determined by minimizing the free energy. We discuss implications for cosmological setups involving complex scalars, with applications to axion dark matter, baryogenesis, kination domination, and gravitational wave production.

preprint2022arXiv

Effective treatment of U(1) gauge field and charged particles in axion inflation

The axionic inflaton with the Chern-Simons coupling may generate U(1) gauge fields and charged particles simultaneously. In order to incorporate the backreaction from the charged particles on the gauge fields, we develop a procedure to obtain an equilibrium solution for the gauge fields by treating the induced current as effective electric and magnetic conductivities. Introducing mean field approximation, and numerically solving self-consistency equations, we find that the gauge field amplitudes are drastically suppressed. Interestingly, as the production becomes more efficient, the charged particles gain a larger part of the transferred energy from the inflaton and eventually dominate it. Our formalism offers a basis to connect this class of inflationary models to a rich phenomenology such as baryogenesis and magnetogenesis.

preprint2022arXiv

Leptoflavorgenesis: baryon asymmetry of the Universe from lepton flavor violation

Charged-lepton flavor violation (CLFV) is a smoking-gun signature of physics beyond the Standard Model. The discovery of CLFV in upcoming experiments would indicate that CLFV processes must have been efficient in the early Universe at relatively low temperatures. In this letter, we point out that such efficient CLFV interactions open up new ways of creating the baryon asymmetry of the Universe. First, we quote the two-loop corrections from charged-lepton Yukawa interactions to the chemical transport in the Standard Model plasma, which imply that nonzero lepton flavor asymmetries summing up to $B-L = 0$ are enough to generate the baryon asymmetry. Then, we describe two scenarios of what we call {\it leptoflavorgenesis}, where efficient CLFV processes are responsible for the generation of primordial lepton flavor asymmetries that are subsequently converted to a baryon asymmetry by weak sphaleron processes. Here, the conversion factor from lepton flavor asymmetry to baryon asymmetry is suppressed by charged-lepton Yukawa couplings squared, which provides a natural explanation for the smallness of the observed baryon-to-photon ratio.

preprint2022arXiv

Scalar field couplings to quadratic curvature and decay into gravitons

Any field, even if it lives in a completely hidden sector, interacts with the visible sector at least via the gravitational interaction. In this paper, we show that a scalar field in such a hidden sector generically couples to the quadratic curvature via dimension five operators, i.e., $ϕR^2$, $ϕR^{μν} R_{μν}$, $ϕR_{μνρσ} R^{μνρσ}$, $ϕR_{μνρσ} \tilde R^{μνρσ}$, because they can emerge if there exist particles heavier than it. We derive these scalar couplings to the quadratic curvature by integrating out heavy particles in a systematic manner. Such couplings are of phenomenological interest since some of them induce the scalar decay into the graviton pair. We point out that the decay of a scalar field can produce a substantial amount of the stochastic cosmic graviton background at high frequency since the suppression scale of these operators is given by the mass of heavy particles not by the Planck scale. The resultant graviton spectrum is computed in some concrete models.

preprint2022arXiv

SU($N$)-natural inflation

We study SU($N$) gauge fields that couple to the inflaton through the Chern-Simons term. We provide a general procedure to construct homogeneous, isotropic, and attractor solutions of the gauge fields during inflation. The gauge fields develop various VEVs corresponding to different spontaneous symmetry breaking patterns of SU($N$) where embedded SU($2$) subgroups are broken with the spatial rotation SO($3$) symmetry. As specific examples, we construct the stable solutions for $N = 3$ and $4$. We also numerically solve the gauge field dynamics and confirm that our analytic solutions are complete and attractor. It is straightforward to apply our procedure to the other simple Lie groups.

preprint2021arXiv

Dark Matter Production during the Thermalization Era

We revisit the non-thermal dark matter (DM) production during the thermalization and reheating era after inflation. The decay of inflaton produces high-energy particles that are thermalized to complete the reheating of the Universe. Before the thermalization is completed, DM can be produced from a collision between the high-energy particles and/or the ambient plasma. We calculate the DM abundance produced from these processes for the case where the cross section of the DM production is proportional to the $n$-th power of the center of mass energy. We find that the collision between the high-energy particles is almost always dominant for $n \gtrsim 4$ while it is subdominant for $n \lesssim 2$. The production from the ambient plasma is dominant when $n \lesssim 3$ and the reheating temperature is of the order of or larger than the DM mass. The production from a collision between the high-energy particle and the ambient plasma is important for $n \lesssim 2$ and the reheating temperature is much lower than the DM mass.

preprint2020arXiv

Chiral Anomaly, Schwinger Effect, Euler-Heisenberg Lagrangian, and application to axion inflation

Particle production in strong electromagnetic fields is a recurring theme in solid state physics, heavy ion collisions, early universe cosmology and formal quantum field theory. In this paper we discuss the Dirac equation in a background of parallel electric and magnetic fields. We review the Schwinger particle production rate, clarify the emergence of the chiral anomaly equation and compute the induced current of charged fermions. We distinguish the contributions from non-perturbative particle production, from the running of the gauge coupling constant and from non-linearities in the effective QED Lagrangian, and clarify how these contributions arise within a single framework. We apply these results to axion inflation. A Chern-Simons coupling between the pseudoscalar particle driving cosmic inflaton and an abelian gauge group induces a dual production of gauge fields and charged fermions. We show that the resulting scalar and gravitational wave power spectra strongly depend on the fermion mass.

preprint2020arXiv

Gravitational Wave Production right after a Primordial Black Hole Evaporation

We discuss the footprint of evaporation of primordial black holes (PBHs) on stochastic gravitational waves (GWs) induced by scalar perturbations. We consider the case where PBHs once dominated the Universe but eventually evaporated before the big bang nucleosynthesis. The reheating through the PBH evaporation could end with a sudden change in the equation of state of the Universe compared to the conventional reheating caused by particle decay. We show that this "sudden reheating" by the PBH evaporation enhances the induced GWs, whose amount depends on the length of the PBH-dominated era and the width of the PBH mass function. We explore the possibility to constrain the primordial abundance of the evaporating PBHs by observing the induced GWs. We find that the abundance parameter $β\gtrsim 10^{-5} \text{ - }10^{-8}$ for $\mathcal{O}(10^3 \text{ - } 10^5) \, \text{g}$ PBHs can be constrained by future GW observations if the width of the mass function is smaller than about a hundredth of the mass.

preprint2019arXiv

Rapid bound-state formation of Dark Matter in the Early Universe

The thermal decoupling description of dark matter (DM) and co-annihilating partners is reconsidered. If DM is realized at around the TeV-mass region or above, even the heaviest electroweak force carriers could act as long-range forces, leading to the existence of meta-stable DM bound states. The formation and subsequent decay of the latter further deplete the relic density during the freeze-out process on top of the Sommerfeld enhancement, allowing for larger DM masses. While so far the bound-state formation was described via the emission of an on-shell mediator ($W^{\pm}$, $Z$, $H$, $g$, photon or exotic), we point out that this particular process does not have to be the dominant scattering-bound state conversion channel in general. If the mediator is coupled in a direct way to any relativistic species present in the Early Universe, the bound-state formation can efficiently occur through particle scattering, where a mediator is exchanged virtually. To demonstrate that such a virtually stimulated conversion process can dominate the on-shell emission even for all temperatures, we analyze a simplified model where DM is coupled to only one relativistic species in the primordial plasma through an electroweak-scale mediator. We find that the bound-state formation cross section via particle scattering can exceed the on-shell emission by up to several orders of magnitude.

preprint2016arXiv

Electroweak Vacuum Stabilized by Moduli during/after Inflation

It is known that the present electroweak vacuum is likely to be metastable and it may lead to a serious instability during/after inflation. We propose a simple solution to the problem of vacuum instability during/after inflation. If there is a moduli field which has Planck-suppressed interactions with the standard model fields, the Higgs quartic coupling in the early universe naturally takes a different value from the present one. A slight change of the quartic coupling in the early universe makes the Higgs potential absolutely stable and hence we are free from the vacuum instability during/after inflation.

preprint2016arXiv

Fate of Electroweak Vacuum during Preheating

Our electroweak vacuum may be metastable in light of the current experimental data of the Higgs/top quark mass. If this is really the case, high-scale inflation models require a stabilization mechanism of our vacuum during inflation. A possible candidate is the Higgs-inflaton/-curvature coupling because it induces an additional mass term to the Higgs during the slow roll regime. However, after the inflation, the additional mass term oscillates, and it can potentially destabilize our electroweak vacuum via production of large Higgs fluctuations during the inflaton oscillation era. In this paper, we study whether or not the Higgs-inflaton/-curvature coupling can save our vacuum by properly taking account of Higgs production during the preheating stage. We put upper bounds on the Higgs-inflaton/-curvature coupling, and discuss possible dynamics that might relax them.

preprint2016arXiv

Nonthermal Gravitino Production after Large Field Inflation

We revisit the nonthermal gravitino production at the (p)reheating stage after inflation. Particular attention is paid to large field inflation models with a $\mathbb{Z}_2$ symmetry, for which the previous perturbative analysis is inapplicable; and inflation models with a stabilizer superfield, which have not been studied non-perturbatively. It is found that in single-superfield inflation models (without the stabilizer field), nonthermal production of the transverse gravitino can be cosmologically problematic while the abundance of the longitudinal gravitino is small enough. In multi-superfield inflation models (with the stabilizer field), production of the transverse and longitudinal gravitinos is significantly suppressed, and they are cosmologically harmless. We also clarify the relation between the background field method used in the preheating context and the standard perturbative decay method to estimate the gravitino abundance.

preprint2016arXiv

Simple Cosmological Solution to the Higgs Instability Problem in the Chaotic Inflation and Formation of Primordial Black Holes

We revisit the compatibility between the chaotic inflation, which provides a natural solution to the initial condition problem, and the metastable electroweak vacuum, which is suggested by the results of LHC and the current mass measurements of top quark and Higgs boson. It is known that the chaotic inflation poses a threat to the stability of the electroweak vacuum because it easily generates large Higgs fluctuations during inflation or preheating and triggers the catastrophic vacuum decay. In this paper, we propose a simple cosmological solution in which the vacuum is stabilized during chaotic inflation, preheating and after that. This simple solution naturally predicts the formation of primordial black holes. We find interesting parameter regions where the present dark matter density is provided by them. Also, the thermal leptogenesis can be accommodated in our scenario.

preprint2016arXiv

Thermalization Process after Inflation and Effective Potential of Scalar Field

We investigate the thermalization process of the Universe after inflation to determine the evolution of the effective temperature. The time scale of thermalization is found to be so long that it delays the evolution of the effective temperature, and the resulting maximal temperature of the Universe can be significantly lower than the one obtained in the literature. Our results clarify the finite density corrections to the effective potential of a scalar field and also processes of heavy particle production. In particular, we find that the maximum temperature of the Universe may be at most electroweak scale if the reheating temperature is as low as ${\cal O} (1)$ MeV, which implies that the electroweak symmetry may be marginally restored. In addition, it is noticeable that the dark matter may not be produced from thermal plasma in such a low reheating scenario, since the maximum temperature can be smaller than the conventional estimation by five orders of magnitude. We also give implications to the Peccei-Quinn mechanism and the Affleck-Dine baryogenesis.

preprint2016arXiv

WIMP Dark Matter in a Well-Tempered Regime: A case study on Singlet-Doublets Fermionic WIMP

Serious searches for the weakly interacting massive particle (WIMP) have now begun. In this context, the most important questions that need to be addressed are: "To what extent can we constrain the WIMP models in the future?" and "What will then be the remaining unexplored regions in the WIMP parameter space for each of these models?" In our quest to answer these questions, we classify WIMP in terms of quantum number and study each case adopting minimality as a guiding principle. As a first step, we study one of the simple cases of the minimal composition in the well-tempered fermionic WIMP regime, namely the singlet-doublets WIMP model. We consider all available constraints from direct and indirect searches and also the predicted constraints coming from the near future and the future experiments. We thus obtain the current status, the near future prospects and the future prospects of this model in all its generality. We find that in the future, this model will be constrained almost solely by the future direct dark matter detection experiments (as compared to the weaker indirect and collider constraints) and the cosmological (relic density) constraints and will hence be gradually pushed to the corner of the coannihilation region, if no WIMP signal is detected. Future lepton colliders will then be useful in exploring this region not constrained by any other experiments.

preprint2015arXiv

Gravitational Effects on Inflaton Decay

We point out that the inflaton inevitably couples to all non-conformally coupled matters gravitationally through an oscillation in the Hubble parameter or the cosmic scale factor. It leads to particle production during the inflaton oscillation regime, which is most efficient just after inflation. Moreover, the analysis is extended to the model with non-minimal inflaton couplings to gravity, in which the Hubble parameter oscillates more violently. We apply our results to the graviton production by the inflaton: gravitons are also produced just after inflation, but the non-minimal coupling does not induce inflaton decay into the graviton pair.

preprint2015arXiv

On adiabatic invariant in generalized Galileon theories

We consider background dynamics of generalized Galileon theories in the context of inflation, where gravity and inflaton are non-minimally coupled to each other. In the inflaton oscillation regime, the Hubble parameter and energy density oscillate violently in many cases, in contrast to the Einstein gravity with minimally coupled inflaton. However, we find that there is an adiabatic invariant in the inflaton oscillation regime in any generalized Galileon theory. This adiabatic invariant is useful in estimating the expansion law of the universe and also the particle production rate due to the oscillation of the Hubble parameter.

preprint2015arXiv

Particle Production after Inflation with Non-minimal Derivative Coupling to Gravity

We study cosmological evolution after inflation in models with non-minimal derivative coupling to gravity. The background dynamics is solved and particle production associated with rapidly oscillating Hubble parameter is studied in detail. In addition, production of gravitons through the non-minimal derivative coupling with the inflaton is studied. We also find that the sound speed squared of the scalar perturbation oscillates between positive and negative values when the non-minimal derivative coupling dominates over the minimal kinetic term. This may lead to an instability of this model. We point out that the particle production rates are the same as those in the Einstein gravity with the minimal kinetic term, if we require the sound speed squared is positive definite.

preprint2014arXiv

Affleck-Dine Baryogenesis and Dark Matter Production after High-scale Inflation

The discovery of the primordial B-mode polarisation by the BICEP2 experiment indicates inflation with a relatively high energy scale. Taking this indication into account, we propose consistent scenarios to account for the observed baryon and dark matter densities in gravity and gauge mediated supersymmetry breaking models. The baryon asymmetry is explained by the Afflck-Dine mechanism, which requires relatively low reheating temperature to avoid a sizable baryonic isocurvature perturbation. The low reheating temperature then requires non-thermal production of dark matter to account for the correct relic density of dark matter. Our scenarios can account for the observations of baryon and dark matter density in gravity and gauge mediation and predict some parameters, including the mass of dark matter.

preprint2014arXiv

Axion Models with High Scale Inflation

We revisit the cosmological aspects of axion models. In the high-scale inflation scenario, the Peccei-Quinn (PQ) symmetry is likely to be restored during/after inflation. If the curvature of the PQ scalar potential at the origin is smaller than its vacuum expectation value; for instance in a class of SUSY axion models, thermal inflation happens before the radial component of the PQ scalar (saxion) relaxes into the global minimum of the potential and the decay of saxion coherent oscillation would produce too much axion dark radiation. In this paper, we study how to avoid the overproduction of axion dark radiation with some concrete examples. We show that, by taking account of the finite-temperature dissipation effect appropriately, the overproduction constraint can be relaxed since the PQ scalar can take part in the thermal plasma again even after the PQ phase transition. We also show that it can be further relaxed owing to the late time decay of another heavy CP-odd scalar, if it is present.

preprint2014arXiv

Curvaton Dynamics Revisited

We revisit the dynamics of the curvaton in detail taking account of effects from thermal environment, effective potential and decay/dissipation rate for general field values and couplings. We also consider the curvature perturbation generated through combinations of various effects: large scale modulation of the oscillation epoch, the effective dissipation rate and the timing at which the equation of state changes. In particular, we find that it tends to be difficult to explain the observed curvature perturbation by the curvaton mechanism without producing too large non-Gaussianity if the curvaton energy density is dissipated through thermal effects. Therefore, the interaction between the curvaton and light elements in the thermal plasma should be suppressed in order for the curvaton to survive the thermal dissipation.

preprint2014arXiv

Dark Matter Chaotic Inflation in Light of BICEP2

We propose an economical model in which a singlet Z$_2$-odd scalar field accounts for the primordial inflation and the present dark matter abundance simultaneously in the light of recent BICEP2 result. Interestingly, the reheating temperature and the thermal dark matter abundance are closely connected by the same interaction between the singlet scalar and the standard model Higgs. In addition, the reheating temperature turns out to be quite high, $T_\text{R} \gtrsim 10^{12}$ GeV, and hence the thermal leptogenesis is compatible with this model. Therefore, it can be one of the simplest cosmological scenarios.

preprint2014arXiv

Dark Matter Production in Late Time Reheating

We estimate dark matter density for the Universe with a reheating temperature smaller than the mass of dark matter, assuming dark matter to be a weakly interacting massive particle. During the reheating process, an inflaton decays and releases high-energy particles, which are scattered inelastically by the thermal plasma and emit many particles. Dark matters are produced through these inelastic scattering processes and pair creation processes by high-energy particles. We properly take account of the Landau-Pomeranchuk-Migdal effect on inelastic processes and show that the resultant energy density of dark matter is much larger than that estimated in the literature and can be consistent with that observed when the mass of dark matter is larger than $O(100)$ GeV.

preprint2014arXiv

Thermalization after/during Reheating

If reheating of the Universe takes place via Planck-suppressed decay, it seems that the thermalization of produced particles might be delayed, since they have large energy/small number densities and number violating large angle scatterings which decrease the momentum of particles by large amount are inefficient correspondingly. In this paper, we study the thermalization of such "under occupied" decay products in detail, following recent developments in understanding the thermalization of non-abelian plasma. Contrary to the above naive expectation, it is shown that in most cases thermalization after/during reheating occurs instantaneously by properly taking account of scatterings with small angles and of particles with small momenta. In particular, the condition for instantaneous thermalization before the completion of reheating is found to be $α^{8/5} \gg (m_ϕ/ M_{\rm pl}) (M_{\rm pl}^2 Γ_ϕ/ m_ϕ^3)^{1/5}$, which is much milder than that obtained in previous works with small angle scatterings taken into account.

preprint2013arXiv

Axino dark matter with R-parity violation and 130 GeV gamma-ray line

We show that decaying axino dark matter with R-parity violation can explain the observed excess of the 130GeV gamma-ray line from the Galactic center in the Fermi data. The branching fraction of the axino decay into monochromatic photons can be O(1), and constraints from continuum gamma-rays and the anti-proton flux are ameliorated. The Peccei-Quinn scale of $O(10^{13}-10^{14})$\,GeV and the R-parity violation parameter of $O(10^{-12}-10^{-11})$ are cosmologically favored.

preprint2013arXiv

Dissipative Effects on Reheating after Inflation

The inflaton must convert its energy into radiation after inflation, which, in a conventional scenario, is caused by the perturbative inflaton decay. This reheating process would be much more complicated in some cases: the decay products obtain masses from an oscillating inflaton and thermal environment, and hence the conventional reheating scenario can be modified. We study in detail processes of particle production from the inflaton, their subsequent thermalization and evolution of inflaton/plasma system by taking dissipation of the inflaton in a hot plasma into account. It is shown that the reheating temperature is significantly affected by these effects.

preprint2013arXiv

Dynamics of oscillating scalar field in thermal environment

There often appear coherently oscillating scalar fields in particle physics motivated cosmological scenarios, which may have rich phenomenological consequences. Scalar fields should somehow interact with background thermal bath in order to decay into radiation at an appropriate epoch, but introducing some couplings to the scalar field makes the dynamics complicated. We investigate in detail the dynamics of a coherently oscillating scalar field, which has renormalizable couplings to another field interacting with thermal background. The scalar field dynamics and its resultant abundance are significantly modified by taking account of following effects : (1) thermal correction to the effective potential, (2) dissipation effect on the scalar field in thermal bath, (3) non-perturbative particle production events and (4) formation of non-topological solitons. There appear many time scales depending on the scalar mass, amplitude, couplings and the background temperature, which make the efficiencies of these effects non-trivial.

preprint2013arXiv

Fate of $Z_2$ Symmetric Scalar Field

The evolution of a coherently oscillating scalar field with $Z_2$ symmetry is studied in detail. We calculate the dissipation rate of the scalar field based on the closed time path formalism. Consequently, it is shown that the energy density of the coherent oscillation can be efficiently dissipated if the coupling constant is larger than the critical value, even though the scalar particle is stable due to the $Z_2$ symmetry.

preprint2013arXiv

Scalar Trapping and Saxion Cosmology

We study in detail the dynamics of a scalar field in thermal bath with symmetry breaking potential. In particular, we focus on the process of trapping of a scalar field at an enhanced symmetry point through the thermal/non-thermal particle production, taking into account the interactions of produced particles with the standard model particles. As an explicit example, we revisit the saxion dynamics with an initial amplitude much larger than the Peccei-Quinn scale and show that the saxion trapping phenomenon happens for the most cases and it often leads to thermal inflation. We also study the saxion dynamics after thermal inflation, and it is shown that thermal dissipation effect on the saxion can relax the axion overproduction problem from the saxion decay.

preprint2013arXiv

The universe dominated by oscillating scalar with non-minimal derivative coupling to gravity

We study the expansion law of the universe dominated by the oscillating scalar field with non-minimal derivative coupling to gravity as G^{μν} \partial_μ ϕ\partial_ν ϕ. In this system the Hubble parameter oscillates with a frequency of the effective mass of the scalar field, which formerly caused a difficulty in analyzing how the universe expands. We find an analytical solution for power law potentials and interpret the solution in an intuitive way by using a new invariant of the system. As a result, we find marginally accelerated expansion for the quadratic potential and no accelerated expansion for the potential with higher power.

preprint2012arXiv

Boltzmann equation for non-equilibrium particles and its application to non-thermal dark matter production

We consider a scalar field (called $ϕ$) which is very weakly coupled to thermal bath, and study the evolution of its number density. We use the Boltzmann equation derived from the Kadanoff-Baym equations, assuming that the degrees of freedom in the thermal bath are well described as "quasi-particles." When the widths of quasi-particles are negligible, the evolution of the number density of $ϕ$ is well governed by a simple Boltzmann equation, which contains production rates and distribution functions both evaluated with dispersion relations of quasi-particles with thermal masses. We pay particular attention to the case that dark matter is non-thermally produced by the decay of particles in thermal bath, to which the above mentioned formalism is applicable. When the effects of thermal bath are properly included, the relic abundance of dark matter may change by $O(10-100%)$ compared to the result without taking account of thermal effects.