Abstract
This work investigates the bulk viscosity of warm, dense, neutrino-transparent, color-superconducting quark matter, where damping of density oscillations in the kHz frequency range arises from weak-interaction-driven direct Urca processes involving quarks. We study the two-flavor red-green paired color-superconducting (2SC) phase, while allowing for the presence of unpaired strange quarks and blue color light quarks of all flavors. Our calculations are based on the SU(3) Nambu–Jona–Lasinio model, extended to include both vector interactions and the ‘t Hooft determinant term. The primary focus is on how variations in the NJL Lagrangian parameters—specifically, the diquark and vector coupling strengths—affect both the static properties of quark matter, such as its equation of state and composition, and its dynamical behavior, including bulk viscosity and associated damping timescales. We find that the bulk viscosity and corresponding damping timescale can change by more than an order of magnitude upon varying the vector coupling by a factor of two at high densities and by a lesser degree at lower densities. This sensitivity primarily arises from the susceptibility of 2SC matter, with a smaller contribution from modifications to the weak interaction rates. In comparison, changes in the diquark coupling have a more limited impact. The damping of density oscillations in 2SC matter is similar quantitatively to nucleonic matter and can be a leading mechanism of dissipation in merging hybrid stars containing color superconducting cores.
1 Introduction
The exploration of matter within the extreme environments of neutron star mergers presents a fascinating way to understand the fundamental properties of nuclear and quark matter at densities substantially exceeding the nuclear saturation density. Although recent gravitational wave observations, such as the GW170817 event (The LIGO Scientific Collaboration et al., 2017), which was accompanied by electromagnetic counterparts, have only captured the inspiral phase of the merger, it is anticipated that the post-merger phase could become observable with the next-generation of gravitational wave detectors. Nevertheless, numerical simulations of mergers offer clues into the dynamics of the merger process and the spectrum of gravitational waves that are emitted (Faber and Rasio, 2012; Baiotti and Rezzolla, 2017; Baiotti, 2019; Bauswein et al., 2019; Radice and Bernuzzi, 2023; Radice and Hawke, 2024). The gravitational waves, emitted within the first tens to hundreds of milliseconds after the merger, carry unique information about the state of matter under extreme conditions. In particular, the post-merger matter is extremely hot, produces a large number of neutrinos, and, unlike a supernova, is highly neutron-rich.
The accuracy of the predictions of simulations of binary neutron star (BNS) mergers depends on a number of factors, among which we would like to address the need for including dissipation in the typically nondissipative general-relativistic hydrodynamics. Dissipation may control the oscillations that source the gravitational waves. The possible influence of various transport phenomena has been addressed in recent years, including the electrical and thermal conductivities (Alford et al., 2018; Harutyunyan and Sedrakian, 2016; Schmitt and Shternin, 2018; Harutyunyan et al., 2018; Harutyunyan and Sedrakian, 2024; Harutyunyan et al., 2024), bulk and shear viscosities (Alford et al., 2018), etc. In recent years, much of the attention has been focused on the bulk viscosity-driven weak -decay processes on nucleons (Alford and Harris, 2019; Alford et al., 2020; Alford et al., 2021a; Alford et al., 2021b; Alford et al., 2023), hyperons (Alford and Haber, 2021) and quarks (Cruz Rojas et al., 2024; Hernández et al., 2024; Alford et al., 2024). The microscopic computation of rates allows one to estimate the damping timescales using local snapshots of matter conditions provided by simulations (Alford et al., 2019; Alford et al., 2020; Alford et al., 2021b; Alford et al., 2023).
The explicit inclusion of bulk viscosity in numerical simulations is still an area under active development (Most et al., 2022; Celora et al., 2022; Chabanov and Rezzolla, 2025a; Chabanov and Rezzolla, 2025b). See also Hammond et al. (2021), Radice et al. (2022), Camelio et al. (2023a), Camelio et al. (2023b) for related studies that provide a more qualitative assessment of bulk viscous effects based on simulations. Some works, such as (Most et al., 2022; Celora et al., 2022), account for the bulk viscosity within frameworks that evolve the system using ideal hydrodynamics, thus neglecting the back-reaction of the bulk viscosity on the fluid motion. Other studies incorporate bulk viscosity dynamically, but assume it to be constant throughout the evolution (Chabanov and Rezzolla, 2025a; Chabanov and Rezzolla, 2025b). Although many details, such as the dependence of bulk viscosity on temperature, density, and the composition of matter at supranuclear densities, remain uncertain, it is generally expected that bulk viscosity can significantly affect the gravitational wave emission during inspiral (Ripley et al., 2023; Ghosh et al., 2025) and merger. In particular, it may damp the oscillations of the remnant stellar core more rapidly, potentially leading to a faster decay of the gravitational wave signal.
Within this context, the potential presence and behavior of quark matter remains largely unexplored. Studies of cold quark matter in hybrid and strange stars, particularly concerning the damping of -mode oscillations in cold stellar configurations have spanned several decades (Madsen, 1992; Drago et al., 2005; Alford and Schmitt, 2007; Blaschke et al., 2007; Sa’d et al., 2007a; Sa’d et al., 2007b; Huang et al., 2010; Wang et al., 2010; Wang and Shovkovy, 2010). Notably, Alford and Schmitt (2007) examined 2SC quark matter at finite temperatures using a model parameterized by the 2SC gap and chemical potentials, though without enforcing electric charge neutrality. That study focused primarily on non-leptonic processes. Importantly, Alford and Schmitt (2007) also derived the bulk viscosity resulting from coupled non-leptonic and leptonic reactions—an approach similar in spirit to the one used us recently (Alford et al., 2024) and expanded further here.
The investigation of bulk viscosity in quark matter within the context of BNS mergers remains in its early stages (Cruz Rojas et al., 2024; Hernández et al., 2024; Alford et al., 2024). Cruz Rojas et al. (2024), Hernández et al. (2024) focused on bulk viscosity in unpaired quark matter arising from non-leptonic and semi-leptonic weak processes. Cruz Rojas et al. (2024) employed both perturbative QCD and holographic methods to obtain improved weak- and strong-coupling estimates of the bulk viscosity. Hernández et al. (2024) computed the bulk viscosity using the MIT bag model and perturbative QCD. Our previous work (Alford et al., 2024) presented the first calculation of bulk viscosity that incorporates the effects of color superconductivity at intermediate densities and finite temperatures in the BNS context. In doing so, the conditions specific to BNS mergers, such as the charge and color neutrality, the density and temperature dependence of pairing gaps and chemical potentials, were treated self-consistently within the vector-interaction-enhanced NJL model. We assume that, in the temperature range MeV and for densities relevant to the quark cores of neutron stars, the matter remains transparent to neutrinos. This assumption does not contradict the existing studies of the neutrino mean free path in quark matter, although the precise condition for trapping depends on the specific model employed—such as the quark matter equation of state (EoS), composition, and possible pairing patterns (see Carter and Reddy (2000), Steiner et al. (2001), Colvero and Lugones (2014)). We also expect that the bulk viscosity will be significantly suppressed once neutrinos become trapped on microscopic scales and reach equilibrium with the surrounding quark matter, consistent with our findings in the neutrino-trapped regime of nucleonic matter (Alford et al., 2020; Alford et al., 2021b; Alford et al., 2019).
Cruz Rojas et al. (2024), Hernández et al. (2024) found a peak in the bulk viscosity at temperatures MeV, which is due to non-leptonic processes, in agreement with earlier studies of cold quark matter (Alford and Schmitt, 2007). Furthermore, Hernández et al. (2024) identified a potential second peak at higher temperatures, around MeV, which emerges for sufficiently large strange quark masses and is driven by semi-leptonic processes. In contrast, Alford et al. (2024) focuses on the high-temperature regime, MeV, and investigates semi-leptonic processes after verifying that non-leptonic ones equilibrate too rapidly to contribute significantly at these temperatures. Their analysis shows a peak within this range, driven by semi-leptonic Urca processes. Collectively, these studies indicate the presence of two distinct peaks in the bulk viscosity: the first at low temperatures ( MeV) driven by non-leptonic processes, and the second at higher temperatures (on the order of a few MeV) associated with semi-leptonic interactions.
This work aims to build upon our previous study (Alford et al., 2024), preserving the essential modeling framework while focusing on variations of the Lagrangian parameters within the vector-interaction-enhanced NJL model (Bonanno and Sedrakian, 2012). First, we investigate how the static composition and EoS of 2SC quark matter with strange quarks depend on the strengths of the vector and diquark couplings. Second, we examine dynamical properties, such as the bulk viscosity and corresponding damping time scales. In this context, we also analyze how Urca process rates are modified by variations in pairing strength and the repulsive vector interaction.
Before proceeding, we point out that the 2SC phase studied in this work competes with alternatives, which we briefly mention here. One example is the quarkyonic phase which features elements of the baryon spectrum for large momenta and quark spectrum at small momenta (McLerran and Reddy, 2019; Han et al., 2019; Kovensky and Schmitt, 2020; Kojo, 2024; Fujimoto et al., 2024; Bluhm et al., 2025). The quarkyonic phase is expected to occupy the finite-temperature, moderate density segment of the QCD phase diagram and its competition with the 2SC phase depends on several factors that are hard to pin down. Furthermore, it has been realized that the chiral phase transition at low temperatures and moderate densities may proceed via an inhomogeneous phase, instead of the first-order homogeneous transition line (Buballa and Carignano, 2016). This phase may have various realizations such as plane-wave sinusoidal modulation or amplitude-modulated “kink crystal” condensate, with periodically alternating positive and negative values of (real-valued) solitonic profile (Karasawa et al., 2016; Abuki, 2018; Ferrer and de la Incera, 2021; Tabatabaee, 2023; Motta et al., 2025).
Although lattice QCD at finite density is hindered by the sign problem, existing results at zero density (Bazavov et al., 2012; Aoki et al., 2024) allow one to speculate that if the anomaly remains strong, it could suppress the emergence of certain exotic phases through modifications of meson masses (Kono et al., 2021). Conversely, if the anomaly becomes weaker at finite chemical potential, this may favor the appearance or extension of inhomogeneous chiral phases mentioned above. In particular, the axial anomaly contributes to the effective potential governing chiral condensate modulations. For spatially inhomogeneous phases that rely on chiral spirals or solitonic structures, any change in the strength of breaking could alter the relative stability of competing condensate configurations (Carignano and Buballa, 2020; Gao et al., 2022). Below we include the ’t Hooft interaction which breaks the symmetry and effectively resembles the axial anomaly.
The present paper is organized as follows. In Section 2 we discuss the equilibrium state of the 2SC phase, including the key thermodynamic parameters that are required for computing the bulk viscosity of quark matter. In Section 3 we discuss the computations of the semi-leptonic Urca rates in the 2SC phase. Section 4 presents the general formalism for computing the bulk viscosity in quark matter based on Urca processes. Our numerical results are presented in Section 5. We conclude and summarize in Section 6. Our calculations use natural units where .
2 Finite temperature 2SC phase
To describe the properties of 2SC quark matter, we adopt a local vector-interaction-enhanced NJL Lagrangian, which is given by (ignoring electromagnetism)where the quark spinor fields carry color and flavor () indices, the matrix of quark current masses is given by , where are the Gell-Mann matrices in the color space, and , where is the unit matrix in the flavor space. The charge conjugate spinors are defined as and , where is the charge conjugation matrix. The form of Equation 1 differs from the original NJL Lagrangian (which contains only the kinetic energy terms and terms proportional to ). It has been extended by three additional terms: (i) the vector interaction with coupling , which is assumed to be repulsive, (ii) the pairing interaction with coupling , and (iii) the ’t Hooft interaction with coupling , which breaks the symmetry. The numerical values of the parameters of the Lagrangian are MeV, MeV, MeV, , (Rehberg et al., 1996). These parameters of the SU(3) NJL model were determined by fitting to empirical data, ensuring that the model accurately reproduces known properties of mesons and quarks. The specific physical observables that are reproduced by the model are the vacuum masses and decay constants of the pion, kaon, eta-meson, and the mass of the eta prime meson, which is particularly sensitive to the anomaly modeled by the ’t Hooft determinant term. We note that the temperature range where the matter is neutrino transparent is limited to MeV, which implies that the condition is always fulfilled for chemical potentials of quarks of all flavors. We have checked that at the low temperatures considered in this work, there are no cutoff artifacts, as the quantities of interest are dominated by unpaired quarks near their thermally smeared Fermi surfaces. In the temperature regime of interest, the corresponding Fermi energies—approximately 300 MeV for light quarks and 500 MeV for strange quarks—remain well below the ultraviolet cutoff, placing our analysis in a regime where such effects are safely suppressed. Nevertheless, should the cutoff approach the Fermi surface, a renormalization group–based framework could be employed to systematically account for and remove any emerging artifacts (Gholami et al., 2025).
At intermediate densities, the 2SC phase is expected to be the dominant pairing channel. More intricate pairing patterns, such as crystalline and gapless color superconducting states, emerge at the low temperatures typical of cold neutron stars. However, the 2SC phase remains a robust feature at finite temperatures. In the 2SC phase, pairing occurs in a color- and flavor-antisymmetric manner between up () and down () quarks, while strange () quarks remain unpaired. The pairing gap in the quasiparticle spectrum is
The overall quark pair wave function must be antisymmetric under exchange of quarks, according to the Pauli exclusion principle. It is seen that the color part of the wave function is antisymmetric, specifically in the anti-triplet configuration of the SU(3) color; the flavor part is anti-symmetric as well and involves light up and down quarks; finally, the spin part is anti-symmetric, implying spin-0 pairing, that is, the Cooper pairs form a spin-singlet state.
The quark-antiquark condensates are defined as
and the constituent mass of each quark flavor is given by
The quark dispersion relations include energy shifts arising from the quark interaction with the vector mean fields and which are defined as and and are the mean field expectation values of the vector mesons and in quark matter. The quark (thermodynamic) chemical potentials are given bywith and being the baryon and charge chemical potential, andbeing the quark charge matrix in flavor space and
the diagonal generators of the SU(3) color gauge group related to the third and eighth gluons. The values of , , and are determined by the requirement of electrical and color neutrality. For some purposes, such as calculating Fermi-Dirac factors, these can be absorbed into “effective” quark chemical potentials
Starting from the Lagrangian (Equation 1) the partition function and the thermodynamic potential of the 2SC phase can be computed in the mean-field approximation; see, for example, Alford et al. (2008), Rüster et al. (2005), Blaschke et al. (2005), Gómez Dumm et al. (2006), Bonanno and Sedrakian (2012). We find the values of the chemical potentials , , by requiring neutrality, i.e., setting . The mean fields, including the pairing gaps, are also obtained by stationarizing with respect to their values. For the computation of the incompressibility of quark matter, we will need the expression for the pressure which is given bywhere are the quasiparticle energies of quarks in quark matter, is lepton pressure (which is approximated as corresponding to an ideal relativistic gas), is the vacuum pressure, and is an effective bag constant, which controls to deconfinement phase transition density.
In the following, we will explore how the Urca rates and resultant bulk viscosity depend on key parameters in the Lagrangian. In doing so, we will vary the temperature in the range MeV, in which we assume the neutrinos are free-streaming. As a main point of this work, we will vary the two additional couplings in the Lagrangian that enhance the ordinary NJL model. Specifically, to understand the role of pairing interaction strength, we consider two values of diquark coupling and 1.25. Similarly, we vary the vector interaction in the range to highlight its role in the physics of the bulk viscosity of 2SC matter.
Figure 1 shows the density dependence of the gap for fixed temperature MeV and several values of and in the range indicated above. The gaps predicted by the model are much larger than the characteristic temperatures relevant for untrapped neutrinos, i.e., we are working essentially in the limit . Therefore, the excitations are exponentially suppressed, and the variations of the gap with temperature are insignificant. The enhancement of the gap with increasing density may be associated with the increase of the density of states at the Fermi surfaces of quarks, while coupling being fixed. Further, it is seen that the increase of the attractive pairing strength from 1 to 1.25 increases the gap by . Finally, the repulsive vector interaction acts oppositely by reducing the pairing gap by at most a few percent when going from to . The energy gap indicates how effectively the phase space of green and red light quarks is suppressed in weak reactions, thereby limiting their contribution to bulk viscosity. The gap variations across model parameters seen in Figure 1 consistently maintain the effective shutdown of these quark color-flavor channels.
FIGURE 1
Figure 2 shows the pressure of the NJL model for fixed temperature MeV and several values of and . As expected, the pressure increases as both the diquark and vector couplings are increased. The significance of the pressure in the present context is that its derivative with respect to baryon density enters the computation of compressibility of matter, which enters the expression for damping timescale via the bulk viscosity, see Equation 76 below. Note that we use in Equation 9 the value as we are interested only in the pressure derivative.
FIGURE 2
Figure 3 displays the composition of the 2SC phase in -equilibrium for varying diquark and vector couplings. It focuses on the populations of unpaired (blue-colored) up and down quarks, along with blue strange quarks, which are relevant for the computation of the bulk viscosity. Note that blue strange quarks have slightly different densities compared to their red-green counterparts, with chemical potential differences of approximately (roughly 10 MeV). The substantial strange quark population rapidly equilibrates with the light flavor sector through non-leptonic processes. Additionally, strange quarks are shown to reduce the lepton fraction in matter, a phenomenon that was identified decades ago (Duncan et al., 1983; Duncan et al., 1984). Examining the variations in coupling constants, we observe that increasing the attractive pairing strength from 1 to 1.25 slightly raises the -quark population while reducing the -quark population, accompanied by increased electron and muon populations. Additionally, strengthening the repulsive vector interaction increases the strange and -quark populations while reducing the -blue quark populations. Through -equilibrium, electron and muon populations are suppressed as the vector coupling increases.
FIGURE 3
Before discussing the reaction rates, we first examine the (effective) chemical potentials and masses of quarks and leptons shown in Figures 4, 5. For quarks, the masses are modified by medium effects through the chiral condensate, which generally depends on both density and temperature. This comparison reveals to what degree the different particles are relativistic. Nevertheless, the formalism applied in the following sections is fully relativistic. Figure 4 shows the effective chemical potentials defined by (Equation 8) as functions of density for fixed values of temperatures as well as coupling constants and . We show only the chemical potentials of blue quarks, which are relevant for the computation of the bulk viscosity. The temperature dependence of the chemical potentials is weak in the shown density range. In addition, the lower panels show the chemical potentials of electrons and muons, which are connected through -equilibrium conditions to the chemical potentials of quarks. The general trend is that the quark chemical potentials rise with density, as one would expect, and the lepton chemical potentials drop, as the rising fraction of strange quarks means that fewer charged leptons are needed to establish electrical neutrality. Rapid non-leptonic processes enforce down-strange flavor equilibrium, ensuring . Therefore, the observed differences between the effective chemical potentials of and quarks in Figure 4 arise from their respective couplings to and mesons. Increasing the pairing interaction increases the chemical potentials of strange and down quarks while reducing that of up quarks. Simultaneously, the chemical potentials of electrons and muons (which are equal) also increase. Finally, increasing the repulsive vector interaction increases -quark and -quark chemical potentials and decreases the -quark chemical potential. In parallel, through -equilibrium, the electron and muon chemical potentials decrease.
FIGURE 4
FIGURE 5
Figure 5 shows the masses of blue quarks as functions of density. An increase in the diquark coupling leads to a rise in the mass of the and a decrease in the dynamically generated masses of - and -quarks. Additionally, an increase in the vector coupling reduces the -quark mass and raises the - and -quark masses. The temperature dependence of the quark masses is observed to be very weak. Comparing effective quark chemical potentials shown in Figure 4 to the respective masses of quarks, we conclude that the light quarks are ultrarelativistic in the entire density range considered. The strange quarks are mildly relativistic in the low-density limit , where is the saturation density, but become strongly relativistic as the density increases.
3 Urca reaction rates for ude and udse compositions
We examine neutrino-transparent and matter, i.e., matter consisting of or quarks and electrons. Muons are also included in the composition whenever they become energetically favorable. Muonic contribution to the reaction rates is subdominant and will be neglected below; the case of nuclear matter has been studied in Alford et al. (2023). The fundamental (semi-leptonic) -equilibration processes in this system include and -quark decay along with the electron capture processes of the direct Urca typewhere and are the electron neutrino and antineutrino, respectively. The processes (Equations 10–13) proceed exclusively from left to right because in neutrino-transparent matter, neutrinos/antineutrinos can only appear in final states. The rates of the Urca processes (Equations 10–13) can be written aswhere etc., are the Fermi distribution functions of fermions (quarks and leptons) with being the single-particle spectrum of particles with mass , and . The effective chemical potentials for quarks are given by Equation 8, and for leptons we have simply . The mapping between the particle labeling and their momenta is as follows: , , , and .
The spin-averaged relativistic matrix element of the Urca processes readswhere is the Fermi coupling constant and is the Cabibbo angle with . For the matrix element of the Urca processes including -quark, one simply needs to replace with . The twelve-dimensional phase-space integrals in Equations 14, 15 can then be reduced to four-dimensional integrals. We follow the method of Alford et al. (2021b).
Before doing so, we note that, as has previously been noted (Alford and Schmitt, 2007; Cruz Rojas et al., 2024; Hernández et al., 2024), the non-leptonic processes
equilibrate much faster than the other Urca processes. As a result, in the regime where the semi-leptonic Urca process contributes to the bulk viscosity significantly, the non-leptonic equilibration is already completed and one can assume at all times [see (Alford and Schmitt, 2007)]. Under this condition, we have only a single equilibrating quantity , which is the relevant measure of how much the system is driven out of -equilibrium state by a cycle of compression and rarefaction.
Returning to Equations 14, 15, we carry out part of the integrations, as described in Alford et al. (2021b), and obtainwhere , is the electron mass, , with , and ; with , and is the Fermi distribution function of the dimensionless variable . The -functions in Equations 18, 19 imply
The integration variables and are normalized-by-temperature transferred energy and momentum, respectively; the variable is the normalized-by-temperature electron energy, computed from its chemical potential, is the normalized lepton momentum, and is the normalized neutrino/antineutrino energy.
Alternative forms of Equations 18, 19 can be obtained by exploiting the relation after which the inner two integrals can be done analytically. We find then
The energy integration limits and are determined by solving Equations 20–22 and then comparing the results to the initial energy integration bounds specified in Equations 18, 19.
As discussed above, the light-flavor quarks are ultrarelativistic under the considered conditions, therefore, we take this limit in Equations 20–22, after which they simplify to (note that ).which, together with the limits of integration in Equations 18, 19, imply
Further simplifications of Equations 23, 24 can be achieved in the low temperature limit, which corresponds to . In this case , , but because . Using the limiting formulawe can approximate the logarithms in Equations 23, 24 asafter which we obtain
It is interesting to note that these expressions can be rewritten in terms of the shift between the chemical potentials of - and -quarks (this shift vanishes in the case discussed here as and quarks are coupled to the -meson due to Equation 8, but it might be nonzero if isovector channel of interaction is included in analogy to the -meson in the case of nucleonic matter discussed in Alford et al. (2021b)). Indeed, using the relationwe obtainwhich demonstrates that for the -quark decay rate vanishes in the ultrarelativistic limit according to (Equation 38). For the -capture rate, we findwhere we made a variable change and made sure that in order to compute the derivative with respect to . If we now assume that -equilibrium is established with , and also we obtain from Equation 40where we dropped the smaller than terms, and replaced the lower boundary of integration by (as the logarithm is suppressed at small values of ).
To compute the bulk viscosity, consider next small departures from -equilibrium . Then, the net -quark production rate can be approximated as , where is the equilibration coefficient defined as
Note that when computing , we also take the derivative of the terms containing .
The integrals in Equations 41, 42 can be computed numerically, after which we obtain
We observe that, as has been noted in other treatments of Urca processes in quark matter (Burrows, 1980; Duncan et al., 1984), the rate contains an additional power of as compared to the low- Urca reaction rates of nonrelativistic baryons. This is because the , , and are all ultrarelativistic, so their Fermi momenta are on the borderline between the Urca process being allowed and being forbidden: the and momenta have to be collinear in order to create a . It is only the thermal blurring of the Fermi surfaces that creates phase space for the process to occur. Consequently, the product in Equation 16 contributes an additional power of beyond the scaling that emerges when the phase space includes a range of angles even in the low temperature limit. It is worth noting that the product does not introduce an additional power of because the angle between the -quark and neutrino momenta can be arbitrary, since the neutrino is thermal.
Thus, the low-temperature direct Urca rates for light quarks (when the isospin chemical shift is absent) differ qualitatively from those of massive particles. For example, Urca reactions involving the quark instead of the quark will have low-temperature rates similar to the nucleonic Urca rates (Alford et al., 2023)
The relevant coefficient in the low- limit is given by
In this limit we have also the relations , , , therefore for the square brackets in Equation 45 we can write , where we used the chemical equilibrium condition . Numerically, we find that , therefore, the square brackets in Equation 45 can be approximated as , which leads to the simple expressions
To estimate the rate of non-leptonic processes (Equation 17) we will use below the low-temperature expression (Madsen, 1993).
4 Bulk viscosity of matter
In this section, we will derive the bulk viscosity that arises from processes (Equations 10–13) in matter, i.e., matter consisting of quarks and electrons, with paired red-green light quarks; the unpaired excitations are blue light quarks, strange quarks of all colors, and leptons (electrons, and muons if energetically favored). Consider small-amplitude density oscillations with a frequency . Separating the oscillating parts from the static equilibrium values of particle densities, we can write , where , and labels the particles.
Oscillations drive the system out of chemical equilibrium, leading to nonzero chemical imbalance in the case of matter. To include strange quarks, note that the non-leptonic reaction proceeds much faster than the Urca processes; therefore, the relation always holds, and the shift in chemical potentials is given by , which can be written aswhere the susceptibilities and are given aswhere , and the derivatives are computed in the beta-equilibrium state. Note that off-diagonal elements do not vanish because of strong interactions between quarks.
If the weak processes were switched off, then the number of all particle species would be conserved separately, which implieswhere is the fluid velocity divergence. Once the weak reactions are switched on, there is a net production of particles that should be the neutrino production rates by quarks, given bywhich define the equilibration coefficients and . Therefore, the rate equations for each fermion species can be written in this case aswhere denotes the rate of the non-leptonic reaction , which is driven by a nearly vanishing chemical potential difference, . Despite its small magnitude, this shift cannot be neglected because the corresponding -coefficient may be very large; see Jones (2001) for a discussion of this point.
Among the resulting balance equations, only one is independent due to the presence of three constraints: charge neutrality (both color and electric) and baryon number conservation. These constraints can be expressed in the formwhere we denote with the densities of only blue quarks, and with – the summed densities of red and green quarks, and is the baryon density. Then we findwhere .
Substituting Equation 61 into Equation 55, assuming and using Equation 50 for , , and (the paired quarks and muons do not participate in reactions), we obtainwhere is the summed rate of the -quark production by Equations 10–13. To eliminate from this equation, we use the condition of chemical equilibrium with respect to non-leptonic reactionwhich can be solved for , and where we introduced shorthand notations
Substituting from Equation 63 back into Equation 62 and using Equation 50 gives
Subtracting from this expression we obtain the nonequilibrium shift.
Solving Equation 63 for (recall that the non-equilibrium shifts of , and the paired quarks are zero) we find
Then the bulk viscous pressure will be given by (using short-hand notation )where we used Equations 59, 60. Assuming isothermal perturbations to compute the pressure derivatives, and further using Equation 66 and the symmetry relation we find
The bulk viscosity is the real part of and is thus defined aswhich has the classic resonant form depending on two quantities: the susceptibility prefactor , which depends only on the EoS, and the relaxation rate , which depends on the EoS and the microscopic interaction rates. In the case where the non-diagonal susceptibilities can be neglected, the quantities and are given bywhere we redefined which are computed in chemical equilibrium. If we neglect the contribution from -quarks, then , , and we find the appropriate quantities for the quark matter
We note that can be interpreted as the beta-disequilibrium–-quark-fraction susceptibility: it quantifies how the out-of-beta-equilibrium chemical potential responds to a change in the -quark fraction. Similarly, is the beta-disequilibrium–baryon-density susceptibility: it characterizes the response of the out-of-beta-equilibrium chemical potential to a change in the baryon density, while keeping the particle fractions fixed.
5 Numerical results
5.1 Equilibration coefficients and Urca rates
As expected, the non-leptonic processes described by Equation 17 occur at significantly higher rates than the Urca processes. This is evident in Figure 6, which shows the equilibration coefficients for both the Urca processes (Equations 10–13) and the non-leptonic processes (Equation 17). Notably, in the temperature range MeV the rates of the -Urca and -Urca processes are comparable for the matter composition predicted by the vector-enhanced NJL model. As a result, the bulk viscosity associated with the nonleptonic channels is expected to peak at much lower temperatures, well below the relevant range for BNS mergers, whose oscillation frequencies lie in the kHz regime. This observation supports our assumption that, in the temperature range MeV, the bulk viscosity can be reliably calculated from the Urca processes alone, under the additional condition . The effects of diquark and vector couplings enter the equilibration coefficients through the composition (chemical potentials, pairing gap, etc.) of the participating particles. It is seen that the combined effect on the equilibration coefficients is to decrease both non-leptonic and Urca processes with increasing vector coupling. The changes associated with variations of the diquark coupling are insignificant.
FIGURE 6
Analytical expressions for and coefficients were derived in the low-temperature limit, see Equations 43, 46 under the assumptions that , and . They are accurate to within a few percent when compared to numerical results that do not rely on these approximations and provide insight into the scaling of these coefficients with various parameters. In particular, it is seen that is suppressed by a factor of compared to which has ; it also has a smaller numerical prefactor from the phase space integration. Due to the influence of interactions on the chemical potentials of light quarks mediated by the isoscalar -field, the -equilibrium for massless particles requires that . This condition implies that the direct Urca process is only thermally allowed. In contrast, the -Urca channels are kinematically open, with a substantial available energy range given by , primarily due to the large mass of the strange quark (Duncan et al., 1983; Duncan et al., 1984). We previously discussed the additional power of arising from the matrix element involving the four-product of the -quark and electron momenta, which are massless. In contrast, the phase-space contributes the standard scaling, which would be the only temperature-dependent component if all particles were massive (implicit weak dependence of masses and other thermodynamic parameters is understood). As a consequence, the difference in temperature scaling leads to the intersection of the and functions at a temperature in the MeV range, as seen in Figure 6.
We show in Figure 7 the temperature dependence of the total -relaxation rate due to the Urca processes (Equations 10–13) as a function of temperature for different fixed values of density, diquark, and vector couplings. It follows a power law scaling in the temperature range MeV. A larger vector coupling enhances the rate slightly, this enhancement being more pronounced at higher density. There are no significant changes with the diquark coupling in the limit . Due to the Lorentzian structure of bulk viscosity in the frequency domain, see Equation 72, its maximum is located at the “resonance” temperature determined by equating to the characteristic angular frequency . For illustration, we consider two representative frequencies. At kHz, the intersection occurs near MeV, while for kHz, it shifts to approximately MeV.
FIGURE 7
5.2 Bulk viscosity and damping timescales
Next, we turn to the main quantity of interest of this study - the bulk viscosity of neutrino-transparent quark matter in the 2SC phase. Figure 8 shows the temperature dependence of bulk viscosity for density oscillations at two frequencies kHz and kHz and at two baryon densities and . As a key feature, the figure also shows the influence of the variation of the vector coupling . The diquark coupling is set to ; results for were also calculated but are not shown because they are almost identical to those for . This is understandable, since for either of these values of , the gap in the spectrum of the paired quarks is much greater than the temperature, so they are frozen out.
FIGURE 8
The curves of bulk viscosity as a function of temperature have the Lorentzian shape expected from (72). For density oscillations of angular frequency , the maximum is reached when , so its location along the axis is determined by the temperature dependence of the relaxation rate .
As density increases, the maximum slowly moves to lower temperatures because the relaxation rate rises with density (as the increase in phase space near the Fermi surfaces leads to faster Urca rates), so is achieved at lower temperatures. The leftward shift in the curve is small because rises quickly, roughly as , so it only takes a small reduction in to compensate for the effect of rising density. This temperature dependence differs from the scaling seen in Hernández et al. (2024) because it is a combination of scaling for the channel (Equation 44) and for the Urca channel (Equation 43). The position of the maximum is even less affected by changes in the vector coupling , since it has only a small effect on the Urca rates, see Figure 6. We also find that the effect of changing is quite small, therefore, we do not show this case explicitly.
The overall scale of the bulk viscosity curve (e.g., the value attained at the maximum) is controlled by the prefactor in Equation 72. The susceptibilities and are affected by the density and the couplings. As density increases, the overall scale decreases considerably, mainly due to a decrease in the susceptibility of 2SC matter. We also see a noticeable effect of the couplings. At lower density , doubling the vector coupling reduces the bulk viscosity at all temperatures by a factor of 3. At higher density , doubling the vector coupling reduces the bulk viscosity by more than an order of magnitude. Without explicitly showing in the figure, let us point out that increasing the magnitude of the diquark coupling leads to a further decrease in the bulk viscosity, which is more pronounced for a larger density . The magnitude of the decrease is about a factor of two. The temperature dependence of the bulk viscosity is seen to be self-similar for the curves shown. From Equation 72 the rate goes as well below the maximum and as well above the maximum, so from the scaling given in the previous paragraph one would expect and respectively. However, there is also a small temperature dependence from the susceptibilities, so the scaling is closer to on the descending side.
We also reiterate that, as noticed in Alford et al. (2024), including/excluding the strange quarks increases/decreases the total relaxation rate by a factor of 1–2, resulting in only a slight shift of the viscosity peak toward lower/higher temperatures.
For comparison, we also show the bulk viscosity of neutrino-transparent nucleonic matter, as calculated in Alford et al. (2023), using the DDME2 (dotted line) and NL3 (dash-dotted line) models at a baryon density of . The key difference between these two models lies in the behavior of the direct Urca process at low temperatures: it is allowed for the NL3 model but remains blocked for DDME2 at this density. Consequently, the bulk viscosities predicted by the two models differ markedly at temperatures MeV, but are rather similar at MeV. In this temperature range, the bulk viscosity of 2SC quark matter is similar to that of nucleonic matter, and variations in the vector or diquark couplings do not significantly affect this outcome.
At lower temperatures, the bulk viscosity of 2SC quark matter falls between the values predicted by the two nucleonic models, regardless of the specific values of the vector and diquark couplings. This behavior can be attributed to the -quark fraction being slightly below, but very close to the threshold for the direct Urca processes described in Equations 10, 11. Specifically, the difference in Fermi momenta between initial-state and final-state particles, , is less than 1 MeV for densities in the range . As a result, thermal smearing at temperatures as low as MeV is sufficient to provide the phase space required for light-quark Urca processes to occur.
Since the bulk viscosity arises from semi-leptonic weak interactions involving the unpaired blue quarks, one would expect similar behavior in unpaired quark matter, specifically a peak in the bulk viscosity at MeV-scale temperatures. This feature is indeed observed in Alford and Schmitt (2007) (Figure 8), Sa’d et al. (2007b) (Figure 3) and in Hernández et al. (2024) (Figure 4), although the maximum viscosity reported in those studies is smaller due to their assumption of a lighter strange quark.
In the aftermath of a BNS merger, the resulting remnant undergoes rapid and large-amplitude density oscillations, driven by strong differential rotation, internal thermal gradients, and potentially turbulent motion. The remnant’s temperature is highly variable over its density range, and some regions may pass through or remain within the temperature window where the bulk viscous damping is especially efficient. If so, the bulk viscosity can significantly damp oscillation modes, converting mechanical energy into heat (or radiated neutrinos) and altering the thermal evolution of the remnant. The characteristic time scale of damping of local oscillations (on the hydrodynamic scales characterized by fixed density, temperature, entropy, etc.) is given by (Alford et al., 2018; Alford and Harris, 2019; Alford et al., 2020)where is the incompressibility and is the pressure. Using our results for the bulk viscosity, we have estimated the damping timescale from Equation 76. The results are shown in Figure 9. As , we see an inversion in the temperature dependence in the damping time scale, which implies that it is shortest at the resonant maximum of the bulk viscosity. We have verified that the compressibility in (Equation 76) is insensitive to the temperature within the relevant range MeV. Turning to the dependence of the damping timescale on the diquark and vector couplings, we find that its behavior is determined by that of the bulk viscosity . At lower density , varying the vector coupling by a factor of two changes the maximum value of by roughly a factor of 3, with the shortest damping times corresponding to the smallest values of the vector and diquark couplings. At higher density (), the damping timescale becomes even more sensitive to the vector coupling: varying the coupling by a factor of two changes by more than an order of magnitude. While a smaller diquark coupling also leads to shorter damping times for a fixed vector coupling, the variation it introduces is, in comparison, not large.
FIGURE 9
Overall, we find that the damping timescales range from a few milliseconds up to several hundred milliseconds, which are comparable to the characteristic timescales of the short-term evolution of the post-merger remnant. For longer-lived remnants, which can survive for several seconds or more, bulk viscous dissipation may play an even more significant role in their dynamical evolution.
Furthermore, a comparison with the results obtained for nucleonic matter reveals that the damping timescales in 2SC quark matter are remarkably similar. This suggests that, within the framework of our NJL model, it would be challenging to distinguish 2SC quark matter from nuclear matter based solely on their bulk viscous dissipation properties in the temperature range of 1–10 MeV.
6 Conclusion
Expanding on our previous work (Alford et al., 2024), we have studied the effects of variations in diquark and vector couplings–specifically and – on the input parameters used to compute the weak interaction rates and the bulk viscosity of the 2SC phase of finite-temperature quark matter within the vector-enhanced NJL model. We focus on the low temperature range and the baryon density range , within which neutrinos are expected to be free-streaming.
The primary effect of varying the vector coupling is a shift in the quark chemical potentials, which affects material properties such as bulk viscosity in two ways. Firstly, it changes the relevant static susceptibilities in the EoS; secondly, it changes flavor-changing Urca rates by altering the available phase space. In contrast, changing the diquark coupling mainly affects the degree of suppression of red-green quark contributions by changing the pairing gap. However, because red-green quarks are indirectly coupled to the unpaired blue-quark sector through -equilibrium and charge neutrality conditions, their influence on the overall thermodynamics is nontrivial and cannot simply be neglected.
Quantitatively, we find that varying the vector coupling by a factor of two changes the bulk viscosity and corresponding damping timescale by a factor of 3–20 at densities from to This sensitivity primarily arises from the susceptibility of 2SC matter, see Equation 74, with a smaller contribution from modifications to the weak interaction rates. In comparison, changes in the diquark coupling have a more limited impact. The bulk viscosity and damping time for 2SC quark matter closely resemble those of nucleonic matter (Alford et al., 2023), making it difficult to distinguish these states via their bulk viscous behavior.
It is worth noting that at temperatures below the range studied here, , one can no longer rely on our assumption that the nonleptonic flavor equilibration process can be treated as instantaneous. When its finite rate is taken into account, the bulk viscosity shows a second peak at MeV (Alford and Schmitt, 2007; Cruz Rojas et al., 2024; Hernández et al., 2024). This peak is absent in nuclear matter. The peak we observe at temperatures in the MeV range, arising from semi-leptonic processes, is also seen in Hernández et al. (2024) (Figure 4) for their heaviest strange quark, MeV. In our model the strange quark is even heavier, MeV so the semileptonic peak can be treated as a separate feature. Note that the picture may be quantitatively different in other partly-paired phases. For example, Urca rates in the gapless 2SC phase, as computed in Jaikumar et al. (2006), may allow the red-green quarks to contribute to the Urca rates at levels comparable to blue quarks. Furthermore, as in our setup, blue quarks are unpaired, one could anticipate that semi-leptonic rates for unpaired quark matter may be similar to those for 2SC matter studied here. By extrapolation, this suggests that the unpaired quark matter at to 10 MeV may also be difficult to distinguish from nuclear matter, as was found for the 2SC phase.
Statements
Data availability statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.
Author contributions
MA: Conceptualization, Investigation, Software, Writing – original draft, Writing – review and editing. AH: Conceptualization, Investigation, Software, Writing – original draft, Writing – review and editing. AS: Conceptualization, Investigation, Software, Writing – original draft, Writing – review and editing. ST: Conceptualization, Investigation, Software, Writing – original draft, Writing – review and editing.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This work was partly supported by the U.S. Department of Energy, Office of Science, Office of Nuclear Physics under Award No. DE-FG02-05ER41375 (MA), the Higher Education and Science Committee (HESC) of the Republic of Armenia through “Remote Laboratory” program, Grant No. 24RL-1C010 (AH and AS), the Polish National Science Centre (NCN) Grant 2023/51/B/ST9/02798 (AS and ST) and the Deutsche For\-schungs\-gemeinschaft (DFG) Grant No. SE 1836/6-1 (AS). ST is a member of the IMPRS for “Quantum Dynamics and Control” at the Max Planck Institute for the Physics of Complex Systems, Dresden, Germany, and acknowledges its partial support.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1
AbukiH. (2018). Chiral crystallization in an external magnetic background: chiral spiral versus real kink crystal. Phys. Rev. D.98, 054006. 10.1103/PhysRevD.98.054006
2
AlfordM.HarutyunyanA.SedrakianA. (2019). Bulk viscosity of baryonic matter with trapped neutrinos. Phys. Rev. D.100, 103021. 10.1103/PhysRevD.100.103021
3
AlfordM.HarutyunyanA.SedrakianA. (2020). Bulk viscous damping of density oscillations in neutron star mergers. Particles3, 500–517. 10.3390/particles3020034
4
AlfordM.HarutyunyanA.SedrakianA. (2021b). Bulk viscosity from urca processes: npeμ matter in the neutrino-trapped regime. Phys. Rev. D.104, 103027. 10.1103/PhysRevD.104.103027
5
AlfordM.HarutyunyanA.SedrakianA. (2023). Bulk viscosity from Urca processes: n p e μ matter in the neutrino-transparent regime. Phys. Rev. D.108, 083019. 10.1103/PhysRevD.108.083019
6
AlfordM.HarutyunyanA.SedrakianA.TsiopelasS. (2024). Bulk viscosity of two-color superconducting quark matter in neutron star mergers. Phys. Rev. D.110, L061303. 10.1103/PhysRevD.110.L061303
7
AlfordM. G.BovardL.HanauskeM.RezzollaL.SchwenzerK. (2018). Viscous dissipation and heat conduction in binary neutron-star mergers. Phys. Rev. Lett.120, 041101. 10.1103/PhysRevLett.120.041101
8
AlfordM. G.HaberA. (2021). Strangeness-changing rates and hyperonic bulk viscosity in neutron star mergers. Phys. Rev. C103, 045810. 10.1103/PhysRevC.103.045810
9
AlfordM. G.HaberA.HarrisS. P.ZhangZ. (2021a). Beta equilibrium under neutron star merger conditions. Universe7, 399. 10.3390/universe7110399
10
AlfordM. G.HarrisS. P. (2019). Damping of density oscillations in neutrino-transparent nuclear matter. Phys. Rev. C100, 035803. 10.1103/PhysRevC.100.035803
11
AlfordM. G.SchmittA. (2007). Bulk viscosity in 2SC quark matter. J. Phys.G34, 67–101. 10.1088/0954-3899/34/1/005
12
AlfordM. G.SchmittA.RajagopalK.SchäferT. (2008). Color superconductivity in dense quark matter. Rev. Mod. Phys.80, 1455–1515. 10.1103/RevModPhys.80.1455
13
AokiS.AokiY.FukayaH.HashimotoS.KanamoriI.KanekoT.et al (2024). Axial U(1) symmetry near the pseudocritical temperature in Nf = 2 + 1 lattice QCD with chiral fermions. PoS, 185. 10.22323/1.453.0185
14
BaiottiL. (2019). Gravitational waves from neutron star mergers and their relation to the nuclear equation of state. Prog. Part. Nucl. Phys.109, 103714. 10.1016/j.ppnp.2019.103714
15
BaiottiL.RezzollaL. (2017). Binary neutron-star mergers: a review of Einstein’s richest laboratory. Rept. Prog. Phys.80, 096901. 10.1088/1361-6633/aa67bb
16
BausweinA.BastianN. U. F.BlaschkeD. B.ChatziioannouK.ClarkJ. A.FischerT.et al (2019). Identifying a first-order phase transition in neutron-star mergers through gravitational waves. Phys. Rev. Lett.122, 061102. 10.1103/PhysRevLett.122.061102
17
BazavovA.BhattacharyaT.BuchoffM. I.ChengM.ChristN. H.DingH. T.et al (2012). Chiral transition and U(1)A symmetry restoration from lattice QCD using domain wall fermions. Phys. Rev. D.86, 094503. 10.1103/PhysRevD.86.094503
18
BlaschkeD.FredrikssonS.GrigorianH.ÖztaşA. M.SandinF. (2005). Phase diagram of three-flavor quark matter under compact star constraints. Phys. Rev. D.72, 065020. 10.1103/PhysRevD.72.065020
19
BlaschkeD. B.BerdermannJ.ColangeloP.CreanzaD.De FazioF.FiniR. A.et al (2007). Neutrino emissivity and bulk viscosity of iso-CSL quark matter in neutron stars. Am. Inst. Phys. Conf. Ser.964, 290–295. 10.1063/1.2823866
20
BluhmM.FujimotoY.McLerranL.NahrgangM. (2025). Quark saturation in the qcd phase diagram. Phys. Rev. C111, 044914. 10.1103/PhysRevC.111.044914
21
BonannoL.SedrakianA. (2012). Composition and stability of hybrid stars with hyperons and quark color-superconductivity. A&A539, A16. 10.1051/0004-6361/201117832
22
BuballaM.CarignanoS. (2016). Inhomogeneous chiral symmetry breaking in dense neutron-star matter. Eur. Phys. J. A52, 57. 10.1140/epja/i2016-16057-6
23
BurrowsA. (1980). Beta decay in quark stars. Phys. Rev. Lett.44, 1640–1643. 10.1103/PhysRevLett.44.1640
24
CamelioG.GavassinoL.AntonelliM.BernuzziS.HaskellB. (2023a). Simulating bulk viscosity in neutron stars. I. Formalism. Phys. Rev. D.107, 103031. 10.1103/PhysRevD.107.103031
25
CamelioG.GavassinoL.AntonelliM.BernuzziS.HaskellB. (2023b). Simulating bulk viscosity in neutron stars. II. Evolution in spherical symmetry. Phys. Rev. D.107, 103032. 10.1103/PhysRevD.107.103032
26
CarignanoS.BuballaM. (2020). Inhomogeneous chiral condensates in three-flavor quark matter. Phys. Rev. D.101, 014026. 10.1103/PhysRevD.101.014026
27
CarterG. W.ReddyS. (2000). Neutrino propagation in color superconducting quark matter. Phys. Rev. D.62, 103002. 10.1103/PhysRevD.62.103002
28
CeloraT.HawkeI.HammondP. C.AnderssonN.ComerG. L. (2022). Formulating bulk viscosity for neutron star simulations. Phys. Rev. D.105, 103016. 10.1103/PhysRevD.105.103016
29
ChabanovM.RezzollaL. (2025a). Impact of bulk viscosity on the postmerger gravitational-wave signal from merging neutron stars. Phys. Rev. Lett.134, 071402. 10.1103/PhysRevLett.134.071402
30
ChabanovM.RezzollaL. (2025b). Numerical modeling of bulk viscosity in neutron stars. Phys. Rev. D.111, 044074. 10.1103/PhysRevD.111.044074
31
ColveroG. C.LugonesG. (2014). Neutrino diffusive transport in hot quark matter: a detailed analysis. Phys. Rev. C89, 055803. 10.1103/PhysRevC.89.055803
32
Cruz RojasJ.GordaT.HoyosC.JokelaN.JärvinenM.KurkelaA.et al (2024). Estimate for the bulk viscosity of strongly coupled quark matter using perturbative QCD and holography. Phys. Rev. Lett.133, 071901. 10.1103/PhysRevLett.133.071901
33
DragoA.LavagnoA.PagliaraG. (2005). Bulk viscosity in hybrid stars. Phys. Rev. D.71, 103004. 10.1103/PhysRevD.71.103004
34
DuncanR. C.ShapiroS. L.WassermanI. (1983). Equilibrium composition and neutrino emissivity of interacting quark matter in neutron stars. ApJ267, 358–370. 10.1086/160875
35
DuncanR. C.WassermanI.ShapiroS. L. (1984). Neutrino emissivity of interacting quark matter in neutron stars. II - finite neutrino momentum effects. ApJ278, 806–812. 10.1086/161850
36
FaberJ. A.RasioF. A. (2012). Binary neutron star mergers. Living Rev. rel.15, 8. 10.12942/lrr-2012-8
37
FerrerE. J.de la InceraV. (2021). Magnetic dual chiral density wave: a candidate quark matter phase for the interior of neutron stars. Universe7, 458. 10.3390/universe7120458
38
FujimotoY.KojoT.McLerranL. D. (2024). Momentum shell in quarkyonic matter from explicit duality: a dual model for cold, dense qcd. Phys. Rev. Lett.132, 112701. 10.1103/PhysRevLett.132.112701
39
GaoB.MinamikawaT.KojoT.HaradaM. (2022). Impacts of the U (1)A anomaly on nuclear and neutron star equation of state based on a parity doublet model. Phys. Rev. C106, 065205. 10.1103/PhysRevC.106.065205
40
GholamiH.HofmannM.BuballaM. (2025). Renormalization-group consistent treatment of color superconductivity in the NJL model. Phys. Rev. D.111, 014006. 10.1103/PhysRevD.111.014006
41
GhoshS.HernándezJ. L.Keshari PradhanB.ManuelC.ChatterjeeD.TolosL. (2025). Tidal heating in binary inspiral of strange quark stars. arXiv e-prints arXiv:2504.07659. 10.48550/arXiv.2504.07659
42
Gómez DummD.BlaschkeD. B.GrunfeldA. G.ScoccolaN. N. (2006). Phase diagram of neutral quark matter in nonlocal chiral quark models. Phys. Rev. D.73, 114019. 10.1103/PhysRevD.73.114019
43
HammondP.HawkeI.AnderssonN. (2021). Thermal aspects of neutron star mergers. Phys. Rev. D.104, 103006. 10.1103/PhysRevD.104.103006
44
HanS.MamunM. A. A.LalitS.ConstantinouC.PrakashM. (2019). Treating quarks within neutron stars. Phys. Rev. D.100, 103022. 10.1103/PhysRevD.100.103022
45
HarutyunyanA.NathanailA.RezzollaL.SedrakianA. (2018). Electrical resistivity and hall effect in binary neutron star mergers. Eur. Phys. J. A54, 191. 10.1140/epja/i2018-12624-1
46
HarutyunyanA.SedrakianA. (2016). Electrical conductivity of a warm neutron star crust in magnetic fields. Phys. Rev. C94, 025805. 10.1103/PhysRevC.94.025805
47
HarutyunyanA.SedrakianA. (2024). Thermal conductivity and thermal hall effect in dense electron-ion plasma. Particles7, 967–983. 10.3390/particles7040059
48
HarutyunyanA.SedrakianA.GevorgyanN. T.HayrapetyanM. V. (2024). Electrical conductivity of a warm neutron star crust in magnetic fields: neutron-Drip regime. Phys. Rev. C109, 055804. 10.1103/PhysRevC.109.055804
49
HernándezJ. L.ManuelC.TolosL. (2024). Damping of density oscillations from bulk viscosity in quark matter. Phys. Rev. D.109, 123022. 10.1103/PhysRevD.109.123022
50
HuangX. G.HuangM.RischkeD. H.SedrakianA. (2010). Anisotropic hydrodynamics, bulk viscosities, and r-modes of strange quark stars with strong magnetic fields. Phys. Rev. D.81, 045015. 10.1103/PhysRevD.81.045015
51
JaikumarP.RobertsC. D.SedrakianA. (2006). Direct urca neutrino rate in color superconducting quark matter. Phys. Rev. C73, 042801. 10.1103/PhysRevC.73.042801
52
JonesP. B. (2001). Bulk viscosity of neutron-star matter. Phys. Rev. D.64, 084003. 10.1103/PhysRevD.64.084003
53
KarasawaS.LeeT. G.TatsumiT. (2016). Brazovskii–dyugaev effect on the inhomogeneous chiral transition in quark matter. Prog. Theor. Exp. Phys.2016, 043D02. 10.1093/ptep/ptw025
54
KojoT. (2024). Stiffening of matter in quark-hadron continuity: a mini-review. arXiv e-prints arXiv:2412.20442. 10.48550/arXiv.2412.20442
55
KonoS.JidoD.KurodaY.HaradaM. (2021). The role of the ua(1) breaking term in dynamical chiral symmetry breaking of chiral effective theories. Prog. Theor. Exp. Phys.2021, 093D02. 10.1093/ptep/ptab084
56
KovenskyN.SchmittA. (2020). Holographic quarkyonic matter. J. High Energy Phys.2020, 112. 10.1007/JHEP09(2020)112
57
MadsenJ. (1992). Bulk viscosity of strange quark matter, damping of quark star vibration, and the maximum rotation rate of pulsars. Phys. Rev. D.46, 3290–3295. 10.1103/PhysRevD.46.3290
58
MadsenJ. (1993). Rate of the weak reaction s + u → u + d in quark matter. Phys. Rev. D.47, 325–330. 10.1103/PhysRevD.47.325
59
McLerranL.ReddyS. (2019). Quarkyonic matter and neutron stars. Phys. Rev. Lett.122, 122701. 10.1103/PhysRevLett.122.122701
60
MostE. R.HarrisS. P.PlumbergC.AlfordM. G.NoronhaJ.Noronha-HostlerJ.et al (2022). Projecting the likely importance of weak-interaction-driven bulk viscosity in neutron s tar mergers. Mon. Not. Ras.509, 1096–1108. 10.1093/mnras/stab2793
61
MottaT. F.BernhardtJ.BuballaM.FischerC. S. (2025). New tool to detect inhomogeneous chiral-symmetry breaking. Phys. Rev. D.111, 074030. 10.1103/PhysRevD.111.074030
62
RadiceD.BernuzziS. (2023). Ab-initio general-relativistic neutrino-radiation hydrodynamics simulations of long-lived neutron star merger remnants to neutrino cooling timescales. Astrophys. J.959, 46. 10.3847/1538-4357/ad0235
63
RadiceD.BernuzziS.PeregoA.HaasR. (2022). A new moment-based general-relativistic neutrino-radiation transport code: methods and first applications to neutron star mergers. Mon. Not. Ras.512, 1499–1521. 10.1093/mnras/stac589
64
RadiceD.HawkeI. (2024). Turbulence modelling in neutron star merger simulations. Liv. Rev. Comput. Astrophys.10, 1. 10.1007/s41115-023-00019-9
65
RehbergP.KlevanskyS. P.HüfnerJ. (1996). Hadronization in the su(3) nambu–jona-lasinio model. Phys. Rev. C53, 410–429. 10.1103/PhysRevC.53.410
66
RipleyJ. L.HegadeK. R. A.YunesN. (2023). Probing internal dissipative processes of neutron stars with gravitational waves during the inspiral of neutron star binaries. Phys. Rev. D.108, 103037. 10.1103/PhysRevD.108.103037
67
RüsterS. B.WerthV.BuballaM.ShovkovyI. A.RischkeD. H. (2005). Phase diagram of neutral quark matter: self-Consistent treatment of quark masses. Phys. Rev. D.72, 034004. 10.1103/PhysRevD.72.034004
68
Sa’dB. A.ShovkovyI. A.RischkeD. H. (2007a). Bulk viscosity of spin-one color superconductors with two quark flavors. Phys. Rev. D.75, 065016. 10.1103/PhysRevD.75.065016
69
Sa’dB. A.ShovkovyI. A.RischkeD. H. (2007b). Bulk viscosity of strange quark matter: urca versus nonleptonic processes. Phys. Rev. D.75, 125004. 10.1103/PhysRevD.75.125004
70
SchmittA.ShterninP. (2018). Reaction rates and transport in neutron stars. Astrophys. Space Sci. Libr.457, 455–574. 10.1007/978-3-319-97616-7_9
71
SteinerA. W.PrakashM.LattimerJ. M. (2001). Diffusion of neutrinos in proto-neutron star matter with quarks. Phys. Lett. B509, 10–18. 10.1016/S0370-2693(01)00434-8
72
TabatabaeeM. (2023). Chiral symmetry breaking and phase diagram of dual chiral density wave in a rotating quark matter. Phys. Rev. D.108, 094042. 10.1103/PhysRevD.108.094042
73
The LIGO Scientific Collaboration, The Virgo CollaborationAbbottR.AbbottT.AcerneseF.AckleyK.AdamsC.et al (2017). Gw170817: observation of gravitational waves from a binary neutron star inspiral. Phys. Rev. Lett.119, 161101. 10.1103/PhysRevLett.119.161101
74
WangX.MalekzadehH.ShovkovyI. A. (2010). Nonleptonic weak processes in spin-one color superconducting quark matter. Phys. Rev. D.81, 045021. 10.1103/PhysRevD.81.045021
75
WangX.ShovkovyI. A. (2010). Bulk viscosity of spin-one color superconducting strange quark matter. Phys. Rev. D.82, 085007. 10.1103/PhysRevD.82.085007
Summary
Keywords
neutron stars, neutrino interactions, quark matter, gravitational waves, transport coefficients
Citation
Alford M, Harutyunyan A, Sedrakian A and Tsiopelas S (2025) Bulk viscosity of two-flavor color superconducting quark matter in neutron star mergers. Front. Astron. Space Sci. 12:1648066. doi: 10.3389/fspas.2025.1648066
Received
16 June 2025
Accepted
21 July 2025
Published
10 September 2025
Volume
12 - 2025
Edited by
Masayuki Matsuzaki, Fukuoka University of Education, Japan
Reviewed by
Kouji Kashiwa, Fukuoka Institute of Technology, Japan
Jianmin Dong, Chinese Academy of Sciences, China
Updates
Copyright
© 2025 Alford, Harutyunyan, Sedrakian and Tsiopelas.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Armen Sedrakian, armen.sedrakian@uwr.edu.pl
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.