We analyze the IceCube four-year neutrino data in search of a signal from the Fermi bubbles. No signal is found from the bubbles or from their dense shell, even when taking into account the softer background. This imposes a conservative $\xi_i<8\%$ upper limit on the cosmic-ray ion (CRI) acceleration efficiency, and an $\eta\equiv \xi_e/\xi_i \gtrsim0.006$ lower limit on the electron-to-ion ratio of acceleration efficiencies (at the $2\sigma$ confidence level). For typical $\xi_i$, a signal should surface once the number of IceCube neutrinos increases by $\sim$an order of magnitude, unless there is a $<$PeV cutoff on the CRI spectrum.

The nature of the bipolar, $\gamma$-ray Fermi bubbles (FB) is still unclear, in part because their faint, high-latitude X-ray counterpart has until now eluded a clear detection. We stack ROSAT data at varying distances from the FB edges, thus boosting the signal and identifying an expanding shell behind the southwest, southeast, and northwest edges, albeit not in the dusty northeast sector near Loop I. A Primakoff-like model for the underlying flow is invoked to show that the signals are consistent with halo gas heated by a strong, forward shock to $\sim$keV temperatures. Assuming ion--electron thermal equilibrium then implies a $\sim10^{56}$ erg event near the Galactic centre $\sim7$ Myr ago. However, the reported high absorption-line velocities suggest a preferential shock-heating of ions, and thus more energetic ($\sim 10^{57}$ erg), younger ($\lesssim 3$ Myr) FBs.

Current measurements of the $\gamma$-ray Fermi bubbles (FB) are based on model-dependent tracers, carry substantial systematic uncertainties, and are at some tension with each other. We show that gradient filters pick out the FB edges, which are found to smoothly connect to the bipolar X-ray structure emanating from the Galactic center, thus supporting the interpretation of the FBs as a Galactic-scale phenomenon. The sharp edges facilitate a direct, model-free measurement of the peripheral FB spectrum. The result is strikingly similar to the full FB-integrated spectrum, softened by a power law of index $\eta\simeq (0.2\mbox{--}0.3)$. This is naturally explained, in both hadronic and leptonic models, if cosmic rays are injected at the edge, and diffuse away preferentially at higher energies $E$. The inferred, averaged diffusion function in the (more plausible) leptonic model, $D(E)\simeq 10^{29.5}(E/10\mbox{ GeV})^{0.48\pm0.02}\mbox{ cm}^2\mbox{ s}^{-1}$, is consistent with estimates for Kraichnan-like turbulence. Our results, in particular the minute spatial variations in $\eta$, indicate that the FB edge is a strong, Mach $\gtrsim5$, forward shock.

A new class of neutrino dark energy models is presented. The new models are characterized by the lack of exotic particles or couplings that violate the standard model symmetry. It is shown that these models lead to several concrete predictions for the dark energy equation of state, as well as possible effects on the cosmic structure formation. These predictions, can be verified (or disproved) with future experiments. At this point, the strongest constraints on these models are obtained from big bang nucleosynthesis, and lead to new bounds on the mass of the lightest neutrino.