Effects of Hoyle state de-excitation on νp–process nucleosynthesis and Galactic chemical evolution
Not Available
Engineering topics
Publications and source records attributed to Mathews, Grant J..
Not Available
As neutron stars merge they can approach very high nuclear density. Here, we summarized recent results for the evolution and gravitational wave emission from binary-neutron star mergers using a a variety of nuclear equations of state with and without a crossover transition to quark matter. We discuss how the late time gravitational wave emission from binary neutron star mergers may possibly reveal the existence of a crossover transition to quark matter.
It is anticipated that the gravitational radiation detected in future gravitational wave (GW) detectors from binary neutron star (NS) mergers can probe the high-density equation of state (EOS). In this work, we perform the first simulations of binary NS mergers which adopt various parametrizations of the quark-hadron crossover (QHC) EOS. These are constructed from combinations of a hadronic EOS ($n_b$<2 $n_0$) and a quark-matter EOS ($n_b$>5 $n_0$), where $n_b$ and $n_0$ are the baryon number density and the nuclear saturation density, respectively. At the crossover densities (2 $n_0$ < $n_b$ < 5 $n_0$) the QHC EOSs continuously soften, while remaining stiffer than hadronic and first-order phase transition EOSs, achieving the stiffness of strongly correlated quark matter. This enhanced stiffness leads to significantly longer lifetimes of the postmerger NS than that for a pure hadronic EOS. We find a dual nature of these EOSs such that their maximum chirp GW frequencies $f_{\text{max}}$ fall into the category of a soft EOS while the dominant peak frequencies ($f_{\text{peak}}$) of the postmerger stage fall in between that of a soft and stiff hadronic EOS. An observation of this kind of dual nature in the characteristic GW frequencies will provide crucial evidence for the existence of strongly interacting quark matter at the crossover densities for QCD.
We report on a study of the natural warm inflationary paradigm (WNI). We show two important new results arise in this model. One is that the observational constraints on the primordial power spectrum from the cosmic microwave background (CMB) can be satisfied without going beyond the Planck scale of the effective field theory. The second is that WNI can inevitably provide perfect conditions for the production of primordial black holes (PBHs) in the golden window of black-hole mass range (10 -16 - 10 -11 M ⊙ ) where it can account for all of the dark matter content of the universe while satisfying observational constraints.
We investigate the neutrino flavor change effects due to neutrino self-interaction and shock wave propagation, as well as the matter effects on the neutrino process in core-collapsing supernovae (CCSNe). For the hydrodynamics, we use two models: a simple thermal bomb model and a specified hydrodynamics model for SN1987A. For the presupernova model, we take an updated model, adjusted to explain SN1987A, which employs recent developments in the (n, γ) reaction rates for nuclei near the stability line (A ~ 100). As for the neutrino luminosity, we adopt two different models: equivalent neutrino luminosity and nonequivalent luminosity models. The latter is taken from a synthetic analysis of CCSN simulation data, which quantitatively presented the results obtained by various neutrino transport models. Relevant neutrino-induced reaction rates are calculated using a shell model for light nuclei and a quasiparticle random phase approximation model for heavy nuclei. For each model, we present abundances of the light nuclei ( 7 Li, 7 Be, 11 B, and 11 C) and the heavy nuclei ( 92 Nb, 98 Tc, 138 La, and 180 Ta) produced by the neutrino process. The light nuclei abundances turn out to be sensitive to the Mikheyev–Smirnov–Wolfenstein (MSW) region around O-Ne-Mg layer while the heavy nuclei are mainly produced prior to the MSW region. Through detailed analyses, we find that neutrino self-interaction becomes a key ingredient, in addition to the MSW effect, for understanding the neutrino process and the relevant nuclear abundances. The normal mass hierarchy is shown to be more compatible with the meteorite data. The main nuclear reactions for each nucleus are also investigated in detail.
We study the evolution of rapid neutron-capture process (r-process) isotopes in the galaxy. We analyze relative contributions from core-collapse supernovae (CCSNe), neutron star mergers, and collapsars under a range of astrophysical conditions and nuclear input data. Here we show that, although the r-process in each of these sites can lead to a similar (universal) elemental distribution, the detailed isotopic abundances can differ from one site to another. These differences may allow for the identification of which sources contributed to the early evolution of r-process material in the galaxy. Our simulations suggest that the early evolution was dominated by CCSNe and collapsar r-process nucleosynthesis. This conclusion may be testable if the next generation of observatories can deduce isotopic r-process abundances.
We calculate the Galactic Chemical Evolution of Mo and Ru by taking into account the contribution from νp-process nucleosynthesis. We estimate yields of p-nuclei such as 92,94 Mo and 96,98 Ru through the νp-process in various supernova progenitors based upon recent models. In particular, the νp-process in energetic hypernovae produces a large amount of p-nuclei compared to the yield in ordinary core-collapse SNe. Because of this, the abundances of 92,94 Mo and 96,98 Ru in the Galaxy are significantly enhanced at [Fe/H] = 0 by the νp-process. We find that the νp-process in hypernovae is the main contributor to the elemental abundance of 92 Mo at low metallicity [Fe/H] < −2. Our theoretical prediction of the elemental abundances in metal-poor stars becomes more consistent with observational data when the νp-process in hypernovae is taken into account.
Not Available
The amplitude of the primordial magnetic field (PMF) is constrained from observational limits on primordial nuclear abundances. Within this constraint, it is possible that nuclear motion is regulated by Coulomb scattering with electrons and positrons (e ± s), while e ± s are affected by a PMF rather than collisions. For example, at a temperature of 10 9 K, thermal nuclei typically experience ~ 10 21 scatterings per second that are dominated by very small angle scattering leading to minuscule changes in the nuclear kinetic energy of order $\mathcal{O}$(1) eV. In this paper the upper limit on the effects of a possible discretization of the e ± momenta by the PMF on the nuclear momentum distribution is estimated under the extreme assumptions that the momentum of the e ± is relaxed before and after Coulomb scattering to Landau levels, and that during Coulomb scattering the PMF is neglected. This assumption explicitly breaks the time reversal invariance of Coulomb scattering, and the Maxwell-Boltzmann distribution is not a trivial steady state solution of the Boltzmann equation under these assumptions. Here we numerically evaluate the collision terms in the Boltzmann equation, and show that the introduction of a special direction in the e ± distribution by the PMF generates no directional dependence of the collisional destruction term of nuclei. Large anisotropies in the nuclear distribution function are then constrained from big bang nucleosynthesis. Ultimately, we conclude that a PMF does not significantly affect the isotropy or big bang nucleosynthesis.
Multicomponent relativistic fluids have been studied for decades. However, simulating the dynamics of the particles and fluids in such a mixture has been a challenge due to the fact that such simulations are computationally expensive in three spatial dimensions. Here, we report on the development and application of a multidimensional relativistic Monte Carlo code to explore the thermalization process in a relativistic multicomponent environment in a computationally inexpensive way. As an illustration we simulate the fully relativistic three-dimensional Brownian-motion-like solution to the thermalization of a high-mass particle (proton) in a bath of relativistic low-mass particles (electrons). Here, we follow the thermalization and ultimate equilibrium distribution of the Brownian-like particle as can happen in the cosmic plasma during big-bang nucleosynthesis. We also simulate the thermalization of energetic particles injected into the plasma as can occur, for example, by the decay of massive unstable particles during the big bang.
The core–cusp problem remains a challenging discrepancy between observations and simulations in the standard Λ CDM model for the formation of galaxies. The problem is that Λ CDM simulations predict a steep power-law mass density profile at the center of galactic dark matter halos. However, observations of dwarf galaxies in the Local Group reveal a density profile consistent with a nearly flat distribution of dark matter near the center. A number of solutions to this dilemma have been proposed. In this paper we summarize investigations into the possibility that the dark matter particles themselves self interact and scatter. Such self-interacting dark matter (SIDM) particles can smooth out the dark-matter profile in high-density regions. We also review the theoretical proposal that self-interacting dark matter may arise as an additional Higgs scalar in the 3–3–1 extension of the Standard Model (SM). We present new simulations of galaxy formation and evolution for this formulation of self-interacting dark matter. Current constraints on this self-interacting dark matter are then summarized.
Big-bang nucleosynthesis (BBN) is valuable as a means to constrain the physics of the early universe and it is the only probe of the radiation-dominated epoch. A fundamental assumption in BBN is that the nuclear velocity distributions obey Maxwell-Boltzmann (MB) statistics as they do in stars. Specifically, the BBN epoch is characterized by a dilute baryon plasma for which the velocity distribution of nuclei is mainly determined by the dominant Coulomb elastic scattering with mildly relativistic electrons. One must therefore deduce the momentum distribution for reacting nuclei from the multicomponent relativistic Boltzmann equation. However, the full multicomponent relativistic Boltzmann equation has only recently been analyzed and its solution has only been worked out in special cases. Moreover, a variety of schemes have been proposed that introduce nonthermal components into the BBN environment which can alter the thermal distribution of reacting nuclei. Here, we construct the relativistic Boltzmann equation in the context of BBN. We also derive a Langevin model and perform relativistic Monte-Carlo simulations which clarify the baryon distribution during BBN and can be used to analyze any relaxation from a nonthermal injection. Further, we show by these analyses that the thermalization process leads to a nuclear distribution function that remains very close to MB statistics even during the most relativistic environment relevant to BBN. Hence, the predictions of standard BBN remain unchanged.
Not Available
In this proposal we performed the first relativistic hydrodynamic simulations of the binary neutron stars near coalescence. Calculations were performed in three spatial dimensions and were supplemented with magnetohydrodynamic simulations of the magnetic reconnection as the neutron star fluid responds to the relativistic forces of the binary. The ultimate goal of the proposed work was to develop a theoretical model for the temporal and spectral evolution of cosmological gamma-ray bursts produced during the magnetic reconnection. This is the first time such hydrodynamic calculations were performed and some surprising results were obtained. The neutrons stars were observed to compress and heat in the simulations. They could even collapse to black holes many orbits before coalescence. The question then arose as to whether this strange phenomenon could power a gamma-ray burst. The bulk of computing effort in this proposal was therefore devoted to an analysis of the implications of this possible new gamma-ray burst paradigm. The initial focus of the work performed under this proposal was therefore to utilize (3+1) dimensional and spherical numerical general relativistic hydrodynamics to study the origin, evolution, and parametric sensitivity of the collapse instability. We derived modified conditions of hydrostatic equilibrium for stars in the curved space of quasi-static orbits. We followed the magnetic field evolution for approximately one orbit (10 msec). The magnetic field energy was observed to exponentiate with an e-folding time of about 1 msec. Thus, the field should build up very quickly to a magnitude such that reconnection and back reaction of the fluid inhibits further growth.
We investigate the possibility that inhomogeneous nucleosynthesis may eventually be used to explain the abundances of Li-6, Be-9, and B in Population II stars. The present work differs from previous studies in that we have used a more extensive reaction network. It is demonstrated that in the simplest scenario the abundances of the light elements with A less than or = 7 constrain the separation of inhomogeneities to sufficently small scales that the model is indistinguishable from homogeneous nucleosynthesis and that the abundnace of Li-6, Be-9, and B are then below observations by several orders of magnitude. This conclusion does not depend on the Li-7 constraint. We also examine alternative scenarios which involve a post-nucleosynthesis reprocessing of the light elements to reproduce the observed abundances of Li and B, while allowing for a somewhat higher baryon density (still well below the cosmological critical density). Future B/H measurements may be able to exclude even this exotic scenario and further restrict primirdial nucleosynthesis to approach the homogeneous model conclusions.
We present analytical models for inhomogeneous chemical evolution (ICE) of systems in which the star formation history resembles a series of bursts, localized in space and/or time, with intermittent periods of remixing. The additional parameter of this model is the metallicity increment of bursting subsystems, but this parameter is constrained by the spread in the age-metallicity relation. We apply this model to the solar annulus in the Galactic disk and show that ICE models yield an improved fit to the observed shape of the stellar abundance distribution function (ADF). The G-dwarf problem can be alleviated with ICE models, but infall of metal poor gas and/or some preenrichment of the disk during the epoch of protogalactic evolution is still required to explain the paucity of low-metallicity dwarfs. ICE models also suggest an explanation of the reduced frequency of metal-rich G-dwarfs relative to the predictions of the simple model. It does not seem likely that chemical evolution of the solar annulus proceeded in a medium that was well-mixed at all times.