We summarize the derivation of the finite temperature, finite chemical potential thermodynamic potential in the bag-model approximation to quantum chromodynamics (QCD) that includes a finite $s$-quark mass in the Feynman diagram contributions for both zero-order and two-loop corrections to the quark interaction. The thermodynamic potential for quarks in QCD is a desired ingredient for computations of the equation of state in the early universe, supernovae, neutron stars, and heavy-ion collisions. The 2-loop contributions are normally divergent and become even more difficult in the limit of finite quark masses and finite chemical potential. We introduce various means to interpolate between the low and high chemical potential limits. Although physically well motivated, we show that the infinite series Pad\'e rational polynomial interpolation scheme introduces spurious poles. Nevertheless, we show that lower order interpolation schemes such as polynomial interpolation reproduce the Pad\'e result without the presence of spurious poles. We propose that in this way one can determine the equation of state for the two-loop corrections for arbitrary chemical potential, temperature and quark mass. This provides a new realistic bag-model treatment of the QCD equation of state. We compute the QCD phase diagram with up to the two-loop corrections. We show that the two-loop corrections decrease the pressure of the quark-gluon plasma and therefore increase the critical temperature and chemical potential of the phase transition. We also show, however, that the correction for finite $s$-quark mass in the two-loop correction serves to decrease the critical temperature for the quark-hadron phase transition in the early universe.

The General Antiparticle Spectrometer (GAPS) will carry out a sensitive dark matter search by measuring low-energy ($\mathrm{E} < 0.25 \mathrm{GeV/nucleon}$) cosmic ray antinuclei. The primary targets are low-energy antideuterons produced in the annihilation or decay of dark matter. At these energies antideuterons from secondary/tertiary interactions are expected to have very low fluxes, significantly below those predicted by well-motivated, beyond the standard model theories. GAPS will also conduct low-energy antiproton and antihelium searches. Combined, these observations will provide a powerful search for dark matter and provide the best observations to date on primordial black hole evaporation on Galactic length scales. The GAPS instrument detects antinuclei using the novel exotic atom technique. It consists of a central tracker with a surrounding time-of-flight (TOF) system. The tracker is a one cubic meter volume containing 10 cm-diameter lithium-drifted silicon (Si(Li)) detectors. The TOF is a plastic scintillator system that will both trigger the Si(Li) tracker and enable better reconstruction of particle tracks. After coming to rest in the tracker, antinuclei will form an excited exotic atom. This will then de-excite via characteristic X-ray transitions before producing a pion/proton star when the antiparticle annihilates with the nucleus. This unique event topology will give GAPS the nearly background-free detection capability required for a rare-event search. Here we present the scientific motivation for the GAPS experiment, its design and its current status as it prepares for flight in the austral summer of 2021-22.

We present an updated general purpose nuclear equation of state (EoS) for use in simulations of core-collapse supernovae, neutron star mergers and black hole collapse. This EoS is formulated in the context of Density Functional Theory (DFT) and is generalized to include all DFT EoSs consistent with known nuclear and astrophysical constraints. This EoS also allows for the possibility of the formation of material with a net proton excess ($Y_p > 0.5$) and has an improved treatment of the nuclear statistical equilibrium and the transition to heavy nuclei as the density approaches nuclear matter density. We include the effects of pions in the regime above nuclear matter density and incorporate all of the known mesonic and baryonic states at high temperature. We analyze how a 3-body nuclear force term in the DFT at high densities stiffens the EoS to satisfy the maximum neutron star constraint, however the density dependence of the symmetry anergy and the formation of pions at high temperatures allows for a softening of the central core in supernova collapse calculations leading to a robust explosion. We also add the possibility of a transition to a QCD chiral-symmetry-restoration and deconfinement phase at densities above nuclear matter density. This paper details the physics, and constraints on, this new EoS and presents an illustration of its implementation in both neutron stars and core-collapse supernova simulations. We present the first results from core-collapse supernova simulations with this EoS.

The Gamma-Ray Imager/Polarimeter for Solar flares (GRIPS) is a balloon-borne telescope designed to study solar-flare particle acceleration and transport. We describe GRIPS's first Antarctic long-duration flight in Jan 2016 and report preliminary calibration and science results. Electron and ion dynamics, particle abundances and the ambient plasma conditions in solar flares can be understood by examining hard X-ray (HXR) and gamma-ray emission (20 keV to 10 MeV) with enhanced imaging, spectroscopy and polarimetry. GRIPS is specifically designed to answer questions including: What causes the spatial separation between energetic electrons producing HXRs and energetic ions producing gamma-ray lines? How anisotropic are the relativistic electrons, and why can they dominate in the corona? How do the compositions of accelerated and ambient material vary with space and time, and why? GRIPS's key technological improvements over the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) include 3D position-sensitive germanium detectors (3D-GeDs) and a single-grid, multi-pitch rotating modulator (MPRM) collimator. The 3D-GeDs have spectral FWHM resolution of a few hundred keV and spatial resolution $<$1 mm$^3$. For photons that Compton scatter, usually $\gtrsim$150 keV, the energy deposition sites can be tracked, providing polarization measurements as well as enhanced background reduction. The MPRM single-grid design provides twice the throughput of a bi-grid imaging system like RHESSI. The grid is composed of 2.5 cm thick W/Cu slats with 1-13 mm variable slit pitch, achieving quasi-continuous FWHM angular coverage over 12.5-162 arcsecs. This resolution is capable of imaging the separate magnetic loop footpoint emissions in a variety of flare sizes. (Abstract edited down from source.)

The Wisconsin Plasma Astrophysics Laboratory (WiPAL) is a flexible user facility designed to study a range of astrophysically relevant plasma processes as well as novel geometries that mimic astrophysical systems. A multi-cusp magnetic bucket constructed from strong samarium cobalt permanent magnets now confines a 10 m$^3$, fully ionized, magnetic-field free plasma in a spherical geometry. Plasma parameters of $ T_{e}\approx5$ to $20$ eV and $n_{e}\approx10^{11}$ to $5\times10^{12}$ cm$^{-3}$ provide an ideal testbed for a range of astrophysical experiments including self-exciting dynamos, collisionless magnetic reconnection, jet stability, stellar winds, and more. This article describes the capabilities of WiPAL along with several experiments, in both operating and planning stages, that illustrate the range of possibilities for future users.

A non-minimally coupled scalar field can have, in principle, a negative effective Planck mass squared which depends on the scalar field. Surprisingly, an isotropic and homogeneous cosmological universe with a non-minimally coupled scalar field is perfectly smooth as the rolling scalar field causes the effective Planck mass to change sign and pass through zero. However, we show that any small deviations from homogeneity diverge as the effective Planck mass vanishes, with catastrophic consequences for the cosmology. The physical origin of the divergence is due to the presence of non-zero scalar anisotropic stress from the non-minimally coupled scalar field. Thus, while the homogeneous and isotropic cosmology appears surprisingly sensible when the effective Planck mass vanishes, inhomogeneities tell a different story.

We present an updated and improved equation of state (which we call the NDL EoS) for use in neutron-star structure and supernova simulations. This EoS is based upon a framework originally developed by Bowers & Wilson, but there are numerous changes. Among them are: (1) a reformulation in the context of density functional theory; (2) the possibility of the formation of material with a net proton excess (Ye > 0.5); (3) an improved treatment of the nuclear statistical equilibrium and the transition to heavy nuclei as the density approaches nuclear matter density; (4) an improved treatment of the effects of pions in the regime above nuclear matter density including the incorporation of all the known mesonic and baryonic states at high temperature; (5) the effects of 3-body nuclear forces at high densities; and (6) the possibility of a first-order or crossover transition to a QCD chiral symmetry restoration and deconfinement phase at densities above nuclear matter density. This paper details the physics of, and constraints on, this new EoS and describes its implementation in numerical simulations. We show comparisons of this EoS with other equations of state commonly used in supernova collapse simulations.

We summarize several new developments in the nuclear equation of state for supernova simulations and neutron stars. We discuss an updated and improved Notre-Dame-Livermore Equation of State (NDL EoS) for use in supernovae simulations. This Eos contains many updates. Among them are the effects of 3- body nuclear forces at high densities and the possible transition to a QCD chiral and/or super-conducting color phase at densities. We also consider the neutron star equation of state and neutrino transport in the presence of strong magnetic fields. We study a new quantum hadrodynamic (QHD) equation of state for neutron stars (with and without hyperons) in the presence of strong magnetic fields. The parameters are constrained by deduced masses and radii. The calculated adiabatic index for these magnetized neutron stars exhibit rapid changes with density. This may provide a mechanism for star-quakes and flares in magnetars. We also investigate the strong magnetic field effects on the moments of inertia and spin down of neutron stars. The change of the moment of inertia associated with emitted magnetic flares is shown to match well with observed glitches in some magnetars. We also discuss a perturbative calculation of neutrino scattering and absorption in hot and dense hyperonic neutron-star matter in the presence of a strong magnetic field. The absorption cross-sections show a remarkable angular dependence in that the neutrino absorption strength is reduced in a direction parallel to the magnetic field and enhanced in the opposite direction. The pulsar kick velocities associated with this asymmetry comparable to observed pulsar velocities and may affect the early spin down rate of proto-neutron star magnetars with a toroidal field configuration.