We develop a relativistic multifluid dynamics appropriate for describing neutron star cores at finite temperatures based on Carter's convective variational procedure. The model includes seven fluids, accounting for both normal and superfluid/superconducting neutrons and protons, leptons (electrons and muons) and entropy. The formulation is compared to the non-variational relativistic multifluid hydrodynamics of Gusakov and collaborators and shown to be equivalent. Vortex lines and flux tubes, mutual friction, vortex pinning, heat conduction and viscosity are incorporated into the model in steps after the basic hydrodynamics is described. The multifluid system is then considered at the mesoscopic scale where the currents around individual vortex lines and flux tubes are important, and this mesoscopic theory is averaged to determine the detailed vortex line/flux tube contributions to the macroscopic ``effective'' theory. This matching procedure is partially successful, though obtaining full agreement between the averaged mesoscopic and macroscopic effective theory requires discarding subdominant terms. The matching procedure allow us to interpret the magnetic $H$-field inside a neutron star in a way that is consistent with condensed matter physics literature, and to clarify the difference between this interpretation and that in previous astrophysical works.

We investigate the compressional modes of cold neutron stars with cores consisting of superfluid neutrons, superconducting protons and normal fluid electrons and muons, and crusts that contain superfluid neutrons plus a normal fluid of (spherical) nuclei and electrons. We develop a two-fluid formalism for the core that accounts for leptonic buoyancy, and an analogous treatment for the crust. We adopt the Cowling approximation, neglecting gravitational perturbations, but include all effects of the background space-time. We introduce a phenomenological, easily-modified nuclear equation of state which contains all of the thermodynamic information required to compute the coupled fluid oscillations, with parameters that are constrained by nuclear physics and the requirement that the maximum mass of a neutron star is $\geq 2M_{\odot}$. Using four parametrizations of this equation of state with nuclear compressibilities $K=230$-$280$ MeV, we calculate the Brunt-V\"{a}is\"{a}l\"{a} frequency due to leptonic buoyancy, and find the corresponding $g$-mode frequencies and eigenfunctions. We find that the WKB approximation reproduces $g$-mode frequencies closely. We examine the dependence of $g$-mode frequencies on stellar mass, nuclear compressibility and strength of neutron-proton entrainment, and compare to previous calculations of $g$-mode frequencies due to leptonic buoyancy. We also compute the $p$-mode spectra, confirming previous findings that the two fluids behave as if uncoupled except the case of large entrainment, and show the existence of nearly resonant mode pairs which could lead to nonlinear $p$-$g$ instabilities even at zero temperature.

Many current and future dark matter and neutrino detectors are designed to measure scintillation light with a large array of photomultiplier tubes (PMTs). The energy resolution and particle identification capabilities of these detectors depend in part on the ability to accurately identify individual photoelectrons in PMT waveforms despite large variability in pulse amplitudes and pulse pileup. We describe a Bayesian technique that can identify the times of individual photoelectrons in a sampled PMT waveform without deconvolution, even when pileup is present. To demonstrate the technique, we apply it to the general problem of particle identification in single-phase liquid argon dark matter detectors. Using the output of the Bayesian photoelectron counting algorithm described in this paper, we construct several test statistics for rejection of backgrounds for dark matter searches in argon. Compared to simpler methods based on either observed charge or peak finding, the photoelectron counting technique improves both energy resolution and particle identification of low energy events in calibration data from the DEAP-1 detector and simulation of the larger MiniCLEAN dark matter detector.