Dmitry Gorbunov
Current Organizations / Institutes
The organizations/institutes from the scholar's latest article
Previous Organizations / Institutes
The organizations / institutes from the scholar's other articles in the past ten years
Research Subjects
Research Interests
Dmitry Gorbunov has a total of 1000 co-authors. Here are the top co-authors by number of shared publications:
- Sergei V. Demidov
Institute for Nuclear Research Russian Academy of Sciences 117312 Moscow Russia
36 shared publications - Fedor Bezrukov
Department of Physics and Astronomy, University of Manchester, Manchester M13 9PL, United Kingdom
29 shared publications - Sergey Troitsky
Institute for Nuclear Research
29 shared publications - Igor' I. Tkachev
Institute for Nuclear Research of the Russian Academy of Sciences, Moscow 117312, Russia
18 shared publications - Anton Chudaykin
Université de Genève
17 shared publications - Grigorii I. Rubtsov
Institute for Nuclear Research of the Russian Academy of Sciences Moscow Russia
15 shared publications
Latest Publications
- Open Access
- 26 October 2026
- Published by Springer Nature in The European Physical Journal C
- Vol. 86 (10) , 1231
Measurement of di-muons from 400 GeV/c protons interacting in a thick molybdenum/tungsten target
Di-muon events emanating from a replica of the SHiP target at the CERN SPS (with a 400 GeV/c proton beam) contain a clear signal of $${J \hspace{-1.66656pt}/\hspace{-1.111pt}\psi \hspace{1.111pt}}$$ production. The target consisted of 13 interaction lengths of molybdenum and tungsten, followed by a 2.4m iron hadron absorber. The di-muon rate is in reasonable agreement with Monte Carlo simulation based on Pythia v8. However, the $${J \hspace{-1.66656pt}/\hspace{-1.111pt}\psi \hspace{1.111pt}} $$ rate observed in the data is much lower than predicted by the Monte Carlo simulation, which is advantageous for SHiP since muons with a large $$p_\textrm{T} $$ are more difficult to deflect by the active muon shield. For the center of mass rapidity ( $$y_{\textrm{cm}} $$ ) interval with the largest overlap with the NA50 measurement, $$0.3< y_{\textrm{cm}} < 0.6$$ , we obtain the production cross-section per nucleon $$\sigma ({J \hspace{-1.66656pt}/\hspace{-1.111pt}\psi \hspace{1.111pt}})/A$$ including the branching ratios of $${J \hspace{-1.66656pt}/\hspace{-1.111pt}\psi \hspace{1.111pt}} $$ decays into muons, $$B_{\mu ^+\mu ^-}$$ : $$B_{\mu ^+\mu ^-}\sigma ({J \hspace{-1.66656pt}/\hspace{-1.111pt}\psi \hspace{1.111pt}})/A=(1.18~\pm ~0.04~\pm ~0.10)~\textrm{nb}$$ /nucleon. This is to be compared to the NA50 extrapolated result (using a much thinner target): $$B_{\mu ^+\mu ^-}\sigma ({J \hspace{-1.66656pt}/\hspace{-1.111pt}\psi \hspace{1.111pt}})/A=(0.99~\pm ~0.04)~\textrm{nb}$$ /nucleon. Within the systematic errors, no significant enhancement due to secondary production of $${J \hspace{-1.66656pt}/\hspace{-1.111pt}\psi \hspace{1.111pt}}$$ inside the target is observed. An upper limit at 90% confidence level of possible contributions from secondary collisions of $$<32\%$$ is obtained.- 26 October 2026
- Published by Springer Nature in The European Physical Journal C
- Vol. 86 (10) , 1231
Open AccessMeasurement of di-muons from 400 GeV/c protons interacting in a thick molybdenum/tungsten target
Di-muon events emanating from a replica of the SHiP target at the CERN SPS (with a 400 GeV/c proton beam) contain a clear signal of $${J \hspace{-1.66656pt}/\hspace{-1.111pt}\psi \hspace{1.111pt}}$$ production. The target consisted of 13 interaction lengths of molybdenum and tungsten, followed by a 2.4m iron hadron absorber. The di-muon rate is in reasonable agreement with Monte Carlo simulation based on Pythia v8. However, the $${J \hspace{-1.66656pt}/\hspace{-1.111pt}\psi \hspace{1.111pt}} $$ rate observed in the data is much lower than predicted by the Monte Carlo simulation, which is advantageous for SHiP since muons with a large $$p_\textrm{T} $$ are more difficult to deflect by the active muon shield. For the center of mass rapidity ( $$y_{\textrm{cm}} $$ ) interval with the largest overlap with the NA50 measurement, $$0.3< y_{\textrm{cm}} < 0.6$$ , we obtain the production cross-section per nucleon $$\sigma ({J \hspace{-1.66656pt}/\hspace{-1.111pt}\psi \hspace{1.111pt}})/A$$ including the branching ratios of $${J \hspace{-1.66656pt}/\hspace{-1.111pt}\psi \hspace{1.111pt}} $$ decays into muons, $$B_{\mu ^+\mu ^-}$$ : $$B_{\mu ^+\mu ^-}\sigma ({J \hspace{-1.66656pt}/\hspace{-1.111pt}\psi \hspace{1.111pt}})/A=(1.18~\pm ~0.04~\pm ~0.10)~\textrm{nb}$$ /nucleon. This is to be compared to the NA50 extrapolated result (using a much thinner target): $$B_{\mu ^+\mu ^-}\sigma ({J \hspace{-1.66656pt}/\hspace{-1.111pt}\psi \hspace{1.111pt}})/A=(0.99~\pm ~0.04)~\textrm{nb}$$ /nucleon. Within the systematic errors, no significant enhancement due to secondary production of $${J \hspace{-1.66656pt}/\hspace{-1.111pt}\psi \hspace{1.111pt}}$$ inside the target is observed. An upper limit at 90% confidence level of possible contributions from secondary collisions of $$<32\%$$ is obtained. - Open Access
- 19 June 2026
- Published by Elsevier in Physics Letters B
- Vol. 879
Changing Constraints on ∑mν from SPT-3G 2018 to D1 in the Context of DESI DR2
- 19 June 2026
- Published by Elsevier in Physics Letters B
- Vol. 879
Open AccessChanging Constraints on ∑mν from SPT-3G 2018 to D1 in the Context of DESI DR2
- Open Access
- 19 March 2026
- Published by Springer Nature in Journal of High Energy Physics
- Vol. 2026 (3)
Pion bremsstrahlung in the splitting function formalism and the dark photon production
A bstract: We study the production of hypothetical vector portal mediators, dark photons γ ′ , with masses in the range 0.4–3.5 GeV in the negatively charged pion-proton collisions π − p → γ ′ X via inelastic pion bremsstrahlung and QCD Drell-Yan-like process. In both cases we estimate the value of total dark photon production cross section and obtain the energy distribution for dark photons that could be produced in the NA64h experiment. We also present the mean energies of dark photons produced in the same way by the secondary pions with momenta typical for T2K, DUNE and SHiP experiments.
- 19 March 2026
- Published by Springer Nature in Journal of High Energy Physics
- Vol. 2026 (3)
Open AccessPion bremsstrahlung in the splitting function formalism and the dark photon production
A bstract: We study the production of hypothetical vector portal mediators, dark photons γ ′ , with masses in the range 0.4–3.5 GeV in the negatively charged pion-proton collisions π − p → γ ′ X via inelastic pion bremsstrahlung and QCD Drell-Yan-like process. In both cases we estimate the value of total dark photon production cross section and obtain the energy distribution for dark photons that could be produced in the NA64h experiment. We also present the mean energies of dark photons produced in the same way by the secondary pions with momenta typical for T2K, DUNE and SHiP experiments.
- Open AccessRelaxing DESI DR2 BAO Constraints on $\sum m_ν$ with Planck and SPT-3G 2018 in the Context of SPT D1We present constraints on the sum of neutrino masses $\sum m_ν$ from a dataset incorporating the full SPT-3G 2018 TT/TE/EE+lensing spectra together with Planck PR4 lensing and low-$\ell$ parts of the Planck PR3 spectra. Using it as a baseline for the DESI DR2 BAO measurements, we arrive at a $95\%$ upper limit of $\sum m_ν< 0.11$ eV, relaxing the tension between $\rm Λ$CDM and lower bounds on $\sum m_ν$ based on neutrino oscillation experiments. When including DES Y1 weak lensing information and the Pantheon+ SNIa catalog, the limit is further loosened to $\sum m_ν0$. On contrast, replacing SPT-3G 2018 primary CMB and lensing spectra with ones from the SPT-3G 2019-2020 (D1) release tightens the overall constraint to $<0.082$ eV and pushes the $\sum m_ν$ posterior mode value to zero, indicating a preference for quasi-negative neutrino masses in line with the D1 analysis. This is a curious shift within SPT-3G measurements of the same field taken in 2018 and in 2019-2020 and processed with different analysis pipelines.Open AccessRelaxing DESI DR2 BAO Constraints on $\sum m_ν$ with Planck and SPT-3G 2018 in the Context of SPT D1We present constraints on the sum of neutrino masses $\sum m_ν$ from a dataset incorporating the full SPT-3G 2018 TT/TE/EE+lensing spectra together with Planck PR4 lensing and low-$\ell$ parts of the Planck PR3 spectra. Using it as a baseline for the DESI DR2 BAO measurements, we arrive at a $95\%$ upper limit of $\sum m_ν< 0.11$ eV, relaxing the tension between $\rm Λ$CDM and lower bounds on $\sum m_ν$ based on neutrino oscillation experiments. When including DES Y1 weak lensing information and the Pantheon+ SNIa catalog, the limit is further loosened to $\sum m_ν0$. On contrast, replacing SPT-3G 2018 primary CMB and lensing spectra with ones from the SPT-3G 2019-2020 (D1) release tightens the overall constraint to $<0.082$ eV and pushes the $\sum m_ν$ posterior mode value to zero, indicating a preference for quasi-negative neutrino masses in line with the D1 analysis. This is a curious shift within SPT-3G measurements of the same field taken in 2018 and in 2019-2020 and processed with different analysis pipelines.
- Open Access
Only slightly "Amplifying nonresonant production of dark sector particles in scattering dominance regime''
In recent article [PhysRevD.109.055041 (2024)] it has been argued that production of dark photons -- hypothetical massive vectors -- in nuclear reactors via mixing with the visible photons is considerably enhanced (by a factor of 10) due to Compton scattering of the latter. We revisit the production of dark photons in the reactor environment and find that although the scattering of photons indeed leads to a larger number of produced dark photons, the overall enhancement is considerably smaller than what has been obtained in Ref. [PhysRevD.109.055041 (2024)]. Our findings are validated using GEANT simulations, which take into account oscillations between ordinary and visible photons and interaction of the latter with matter. The correction to the limit on mixing parameter between dark and visible photons is expected to be below 30% for all masses. It's application must be accompanied with an update on the original spectrum of photons produced in nuclear reactions inside the nuclear reactor.
Open AccessOnly slightly "Amplifying nonresonant production of dark sector particles in scattering dominance regime''
In recent article [PhysRevD.109.055041 (2024)] it has been argued that production of dark photons -- hypothetical massive vectors -- in nuclear reactors via mixing with the visible photons is considerably enhanced (by a factor of 10) due to Compton scattering of the latter. We revisit the production of dark photons in the reactor environment and find that although the scattering of photons indeed leads to a larger number of produced dark photons, the overall enhancement is considerably smaller than what has been obtained in Ref. [PhysRevD.109.055041 (2024)]. Our findings are validated using GEANT simulations, which take into account oscillations between ordinary and visible photons and interaction of the latter with matter. The correction to the limit on mixing parameter between dark and visible photons is expected to be below 30% for all masses. It's application must be accompanied with an update on the original spectrum of photons produced in nuclear reactions inside the nuclear reactor.
- Open Access
- 10 December 2025
- Published by American Physical Society (APS) in Physical Review D
- Vol. 112 (12)
Cosmic domain walls on a lattice: Illusive effects of initial conditions
Evolution of cosmic domain walls (DWs) settles to the scaling solution, which is often assumed to be independent of initial conditions. However, lattice simulations performed in this work reveal a clear dependence of the scaling DW area on the initial configuration of the sourcing scalar field, specifically, its infrared (IR) properties. Namely, the DW area grows as one suppresses IR modes in the initial scalar field spectrum. This growth is saturated, when the area parameter ξ commonly used in the literature reaches the value ξ max ≈ 1.2 . The dependence of ξ on IR modes is argued to be of nonphysical origin: it is likely to be due to effects of the lattice boundary. Assuming that the memory of initial conditions is erased, one recognizes ξ ≈ 1.2 obtained in the situation with maximally suppressed IR modes as a genuine universal value of the area parameter in the scaling regime. We demonstrate that ignorance about initial conditions may affect predictions for the energy density of gravitational waves by the factor O ( 5 ) . The spectral shape of gravitational waves is also affected by the choice of initial conditions, most notably in the low-frequency part. Likewise, we revisit annihilation of DWs under the influence of a potential bias. It has been previously found in Babichev [.] that the annihilation happens significantly earlier compared to the estimate based on the simple balance between the potential bias and surface energy density. We further support this observation and show that the tendency toward an earlier annihilation gets even stronger upon removing IR modes in simulations.
- 10 December 2025
- Published by American Physical Society (APS) in Physical Review D
- Vol. 112 (12)
Open AccessCosmic domain walls on a lattice: Illusive effects of initial conditions
Evolution of cosmic domain walls (DWs) settles to the scaling solution, which is often assumed to be independent of initial conditions. However, lattice simulations performed in this work reveal a clear dependence of the scaling DW area on the initial configuration of the sourcing scalar field, specifically, its infrared (IR) properties. Namely, the DW area grows as one suppresses IR modes in the initial scalar field spectrum. This growth is saturated, when the area parameter ξ commonly used in the literature reaches the value ξ max ≈ 1.2 . The dependence of ξ on IR modes is argued to be of nonphysical origin: it is likely to be due to effects of the lattice boundary. Assuming that the memory of initial conditions is erased, one recognizes ξ ≈ 1.2 obtained in the situation with maximally suppressed IR modes as a genuine universal value of the area parameter in the scaling regime. We demonstrate that ignorance about initial conditions may affect predictions for the energy density of gravitational waves by the factor O ( 5 ) . The spectral shape of gravitational waves is also affected by the choice of initial conditions, most notably in the low-frequency part. Likewise, we revisit annihilation of DWs under the influence of a potential bias. It has been previously found in Babichev [.] that the annihilation happens significantly earlier compared to the estimate based on the simple balance between the potential bias and surface energy density. We further support this observation and show that the tendency toward an earlier annihilation gets even stronger upon removing IR modes in simulations.
- 19 November 2025
- Published by American Physical Society (APS) in Physical Review D
- Vol. 112 (10)
Search for a photon peak from keV-scale dark matter annihilation with NuSTAR: Constraints on ⟨ σ v ⟩ after 11 years of observations
- 19 November 2025
- Published by American Physical Society (APS) in Physical Review D
- Vol. 112 (10)
Search for a photon peak from keV-scale dark matter annihilation with NuSTAR: Constraints on ⟨ σ v ⟩ after 11 years of observations
- 1 October 2025
- Published by IOP Publishing in Journal of Cosmology and Astroparticle Physics
- Vol. 2025 (10)
Biased domain walls: faster annihilation, weaker gravitational waves
We study the evolution of domain wall networks and their phenomenological implications in a model of a real scalar χ , where a Z 2 -symmetry is slightly broken by a potential bias V bias . It is demonstrated that the latter triggers domain wall annihilation considerably earlier than previously thought. Namely, we observe that the scaling relation t ann ∝ 1/ V bias 2/3 for the annihilation time t ann fits to the simulation data better than a commonly assumed t ann ∝ 1/ V bias . As a result, the energy density of gravitational waves produced by the network of biased domain walls, for a given tiny V bias , is suppressed compared to naive expectations. The spectral shape of gravitational waves is similar to that resulting from unbiased domain walls, but with more power in the close-to-maximum ultraviolet part. In the far ultraviolet region, the spectrum of gravitational waves becomes nearly flat; such a plateau has been recognised earlier in the case of unbiased walls. In our investigation we mainly focus on the symmetry breaking potential V breaking ∝ χ 3 , and argue that no significant modifications of the domain walls evolution take place if one includes higher powers of χ .
- 1 October 2025
- Published by IOP Publishing in Journal of Cosmology and Astroparticle Physics
- Vol. 2025 (10)
Biased domain walls: faster annihilation, weaker gravitational waves
We study the evolution of domain wall networks and their phenomenological implications in a model of a real scalar χ , where a Z 2 -symmetry is slightly broken by a potential bias V bias . It is demonstrated that the latter triggers domain wall annihilation considerably earlier than previously thought. Namely, we observe that the scaling relation t ann ∝ 1/ V bias 2/3 for the annihilation time t ann fits to the simulation data better than a commonly assumed t ann ∝ 1/ V bias . As a result, the energy density of gravitational waves produced by the network of biased domain walls, for a given tiny V bias , is suppressed compared to naive expectations. The spectral shape of gravitational waves is similar to that resulting from unbiased domain walls, but with more power in the close-to-maximum ultraviolet part. In the far ultraviolet region, the spectrum of gravitational waves becomes nearly flat; such a plateau has been recognised earlier in the case of unbiased walls. In our investigation we mainly focus on the symmetry breaking potential V breaking ∝ χ 3 , and argue that no significant modifications of the domain walls evolution take place if one includes higher powers of χ .
- Open Access
Decays and annihilation of galactic dark matter: determine $D$-, $J_s$-, $J_p$- and $J_d$-factors with dark matter profiles inferred from GravSphere fit to stellar observations
Dark matter mass density profiles and velocity distributions for a set of dwarf spheroidal galaxies (dSphs) have recently been obtained (F.Bezrukov, D.Gorbunov, E.Koreshkova arXiv:2412.20585) by performing a multi-parametric fit to the stellar observations with the help of the GravSphere which solves the Jeans equation. We use these results to calculate the geometrical factors for estimation of the fluxes of cosmic rays expected from decay ($D$-factor) and annihilation ($J_s$-, $J_p$- and $J_d$-factors for $s$-, $p$- and $d$- wave processes) of dark matter particles in galaxies. The general novelty is the account for a possible anisotropy in velocities of dark matter particles. On the basis of this analysis we present empirical scaling approximations to these factors as functions of typical observables: distance to the galaxy $d$, it's half-radius $r_h$ and line-of-sight stellar velocity dispersion $σ_{LOS}$. They can be applied to any galaxy, and for $D$- and $J_s$-factors we refine the estimates of (Andrew B. Pace, Louis E. Strigari arXiv:1802.06811): the shifts in the central values remain within 1-2$σ$ error bars.Open AccessDecays and annihilation of galactic dark matter: determine $D$-, $J_s$-, $J_p$- and $J_d$-factors with dark matter profiles inferred from GravSphere fit to stellar observations
Dark matter mass density profiles and velocity distributions for a set of dwarf spheroidal galaxies (dSphs) have recently been obtained (F.Bezrukov, D.Gorbunov, E.Koreshkova arXiv:2412.20585) by performing a multi-parametric fit to the stellar observations with the help of the GravSphere which solves the Jeans equation. We use these results to calculate the geometrical factors for estimation of the fluxes of cosmic rays expected from decay ($D$-factor) and annihilation ($J_s$-, $J_p$- and $J_d$-factors for $s$-, $p$- and $d$- wave processes) of dark matter particles in galaxies. The general novelty is the account for a possible anisotropy in velocities of dark matter particles. On the basis of this analysis we present empirical scaling approximations to these factors as functions of typical observables: distance to the galaxy $d$, it's half-radius $r_h$ and line-of-sight stellar velocity dispersion $σ_{LOS}$. They can be applied to any galaxy, and for $D$- and $J_s$-factors we refine the estimates of (Andrew B. Pace, Louis E. Strigari arXiv:1802.06811): the shifts in the central values remain within 1-2$σ$ error bars. - Open Access
- 4 August 2025
- Published by American Physical Society (APS) in Physical Review D
- Vol. 112 (3)
Search for dark-matter axions beyond the quantum limit: The cosmological axion Sarov haloscope proposal
Firmly established in astrophysical observations, dark matter evades direct detection in experiments. Axions and axionlike particles are among the leading dark-matter candidates, and numerous attempts to detect them in laboratories have been performed. Here, we propose to advance these efforts substantially, extending the sensitivity for dark-matter axions in the mass range ( 38 – 54 ) μ eV down to the axion-photon couplings g a γ γ ≲ ( 10 − 14 − 10 − 15 ) GeV − 1 , motivated by generic models of quantum chromodynamics axion. Single-photon detectors operating at ultralow temperatures are key elements of the experiment. The projected sensitivity will be reached in one year of data taking with magnetic field of (1–10) T, making Cosmological Axion Sarov Haloscope (CASH) the most sensitive haloscope in this mass range.
- 4 August 2025
- Published by American Physical Society (APS) in Physical Review D
- Vol. 112 (3)
Open AccessSearch for dark-matter axions beyond the quantum limit: The cosmological axion Sarov haloscope proposal
Firmly established in astrophysical observations, dark matter evades direct detection in experiments. Axions and axionlike particles are among the leading dark-matter candidates, and numerous attempts to detect them in laboratories have been performed. Here, we propose to advance these efforts substantially, extending the sensitivity for dark-matter axions in the mass range ( 38 – 54 ) μ eV down to the axion-photon couplings g a γ γ ≲ ( 10 − 14 − 10 − 15 ) GeV − 1 , motivated by generic models of quantum chromodynamics axion. Single-photon detectors operating at ultralow temperatures are key elements of the experiment. The projected sensitivity will be reached in one year of data taking with magnetic field of (1–10) T, making Cosmological Axion Sarov Haloscope (CASH) the most sensitive haloscope in this mass range.