On the strength of a double Nobel prize winning experiment (Super)Kamiokande and an extremely successful long baseline neutrino programme, the third generation Water Cherenkov detector, Hyper-Kamiokande, is being developed by an international collaboration as a leading worldwide experiment based in Japan. The Hyper-Kamiokande detector will be hosted in the Tochibora mine, about 295 km away from the J-PARC proton accelerator research complex in Tokai, Japan. The currently existing accelerator will be steadily upgraded to reach a MW beam by the start of the experiment. A suite of near detectors will be vital to constrain the beam for neutrino oscillation measurements. A new cavern will be excavated at the Tochibora mine to host the detector. The experiment will be the largest underground water Cherenkov detector in the world and will be instrumented with new technology photosensors, faster and with higher quantum efficiency than the ones in Super-Kamiokande. The science that will be developed will be able to shape the future theoretical framework and generations of experiments. Hyper-Kamiokande will be able to measure with the highest precision the leptonic CP violation that could explain the baryon asymmetry in the Universe. The experiment also has a demonstrated excellent capability to search for proton decay, providing a significant improvement in discovery sensitivity over current searches for the proton lifetime. The atmospheric neutrinos will allow to determine the neutrino mass ordering and, together with the beam, able to precisely test the three-flavour neutrino oscillation paradigm and search for new phenomena. A strong astrophysical programme will be carried out at the experiment that will detect supernova neutrinos and will measure precisely solar neutrino oscillation.

Magnetically confined mountains on accreting neutron stars are promising sources of continuous-wave gravitational radiation and are currently the targets of directed searches with long-baseline detectors like the Laser Interferometer Gravitational Wave Observatory (LIGO). In this paper, previous ideal-magnetohydrodynamic models of isothermal mountains are generalized to a range of physically motivated, adiabatic equations of state. It is found that the mass ellipticity drops substantially, from \epsilon ~ 3e-4 (isothermal) to \epsilon ~ 9e-7 (non-relativistic degenerate neutrons), 6e-8 (relativistic degenerate electrons) and 1e-8 (non-relativistic degenerate electrons) (assuming a magnetic field of 3e12 G at birth). The characteristic mass M_{c} at which the magnetic dipole moment halves from its initial value is also modified, from M_{c}/M_{\sun} ~ 5e-4 (isothermal) to M_{c}/M_{\sun} ~ 2e-6, 1e-7, and 3e-8 for the above three equations of state, respectively. Similar results are obtained for a realistic, piecewise-polytropic nuclear equation of state. The adiabatic models are consistent with current LIGO upper limits, unlike the isothermal models. Updated estimates of gravitational-wave detectability are made. Monte Carlo simulations of the spin distribution of accreting millisecond pulsars including gravitational-wave stalling agree better with observations for certain adiabatic equations of state, implying that X-ray spin measurements can probe the equation of state when coupled with magnetic mountain models.

It is widely assumed that the observed reduction of the magnetic field of millisecond pulsars can be connected to the accretion phase during which the pulsar is spun up by mass accretion from a companion. A wide variety of reduction mechanisms have been proposed, including the burial of the field by a magnetic mountain, formed when the accreted matter is confined to the poles by the tension of the stellar magnetic field. A magnetic mountain effectively screens the magnetic dipole moment. On the other hand, observational data suggests that accreting neutron stars are sources of gravitational waves, and magnetic mountains are a natural source of a time-dependent quadrupole moment. We show that the emission is sufficiently strong to be detectable by current and next generation long-baseline interferometers. Preliminary results from fully three-dimensional magnetohydrodynamic (MHD) simulations are presented. We find that the initial axisymmetric state relaxes into a nearly axisymmetric configuration via toroidal magnetic modes. A substantial quadrupole moment is still present in the final state, which is stable (in ideal MHD) yet highly distorted.

The theory of polar magnetic burial in accreting neutron stars predicts that a mountain of accreted material accumulates at the magnetic poles of the star, and that, as the mountain spreads equatorward, it is confined by, and compresses, the equatorial magnetic field. Here, we extend previous, axisymmetric, Grad-Shafranov calculations of the hydromagnetic structure of a magnetic mountain up to accreted masses as high as $\Ma = 6\times 10^{-4}\Msun$, by importing the output from previous calculations (which were limited by numerical problems and the formation of closed bubbles to $\Ma < 10^{-4}\Msun$) into the time-dependent, ideal-magnetohydrodynamic code ZEUS-3D and loading additional mass onto the star dynamically. The rise of buoyant magnetic bubbles through the accreted layer is observed in these experiments. We also investigate the stability of the resulting hydromagnetic equilibria by perturbing them in ZEUS-3D. Surprisingly, it is observed that the equilibria are marginally stable for all $\Ma\leq 6\times 10^{-4}\Msun$; the mountain oscillates persistently when perturbed, in a combination of Alfv\'en and acoustic modes, without appreciable damping or growth, and is therefore not disrupted (apart from a transient Parker instability initially, which expels $< 1 %$ of the mass and magnetic flux).

Matter accreting onto the magnetic poles of a neutron star spreads under gravity towards the magnetic equator, burying the polar magnetic field and compressing it into a narrow equatorial belt. Steady-state, Grad-Shafranov calculations with a self-consistent mass-flux distribution (and a semi-quantitative treatment of Ohmic diffusion) show that, for $\Ma \gtrsim 10^{-5}\Msun$, the maximum field strength and latitudinal half-width of the equatorial magnetic belt are $B_{\rm max} = 5.6\times 10^{15} (\Ma/10^{-4}\Msun)^{0.32}$ G and $\Delta\theta = \max[3^{\circ} (\Ma/10^{-4}\Msun)^{-1.5},3^{\circ} (\Ma/10^{-4}\Msun)^{0.5}(\dot{M}_{\rm a}/10^{-8}\Msun {\rm yr}^{-1})^{-0.5}]$ respectively, where $\Ma$ is the total accreted mass and $\dot{M}_{\rm a}$ is the accretion rate. It is shown that the belt prevents north-south heat transport by conduction, convection, radiation, and ageostrophic shear. This may explain why millisecond oscillations observed in the tails of thermonuclear (type I) X-ray bursts in low-mass X-ray binaries are highly sinusoidal: the thermonuclear flame is sequestered in the magnetic hemisphere which ignites first. The model is also consistent with the occasional occurrence of closely spaced pairs of bursts. Time-dependent, ideal-magnetohydrodynamic simulations confirm that the equatorial belt is not disrupted by Parker and interchange instabilities.

Recent time-dependent, ideal-magnetohydrodynamic (ideal-MHD) simulations of polar magnetic burial in accreting neutron stars have demonstrated that stable, magnetically confined mountains form at the magnetic poles, emitting gravitational waves at $f_{*}$ (stellar spin frequency) and $2 f_{*}$. Global MHD oscillations of the mountain, whether natural or stochastically driven, act to modulate the gravitational wave signal, creating broad sidebands (full-width half-maximum $\sim 0.2f_*$) in the frequency spectrum around $f_{*}$ and $2 f_{*}$. The oscillations can enhance the signal-to-noise ratio achieved by a long-baseline interferometer with coherent matched filtering by up to 15 per cent, depending on where $f_*$ lies relative to the noise curve minimum. Coherent, multi-detector searches for continuous waves from nonaxisymmetric pulsars should be tailored accordingly.

We present a time-lapse sequence of 20 near-infrared (J- and K'-band) snapshots of the central 20x20 arcsec of the Crab pulsar wind nebula, taken at subarcsecond resolution with the Hokupa'a/QUIRC adaptive optics camera on the Gemini North Telescope, and sampled at intervals of 10 min and 24 hr. It is observed that the equatorial wisps and polar knots in the termination shock of the pulsar wind appear to fluctuate in brightness on kilosecond time-scales. Maximum flux variations of +/-24 +/-4 and +/-14 +/-4 per cent relative to the mean (in 1.2 ks) are measured for the wisps and knots respectively, with greatest statistical significance in J band where the nebula background is less prominent. The J and K' flux densities imply different near-infrared spectra for the nonthermal continuum emission from the wisps and outermost polar knot (sprite), giving F_nu ~ nu^{-0.56+/-0.12} and F_nu ~ nu^{-0.21+/-0.13) respectively. The data are compared with existing optical and UV photometry and applied to constrain theories of the variability of the wisps (relativistic ion-cyclotron instability) and knots (relativistic fire hose instability).

The amplitude of the gravitational radiation from an accreting neutron star undergoing polar magnetic burial is calculated. During accretion, the magnetic field of a neutron star is compressed into a narrow belt at the magnetic equator by material spreading equatorward from the polar cap. In turn, the compressed field confines the accreted material in a polar mountain which is misaligned with the rotation axis in general, producing gravitational waves. The equilibrium hydromagnetic structure of the polar mountain, and its associated mass quadrupole moment, are computed as functions of the accreted mass, M_a, by solving a Grad-Shafranov boundary value problem. The orientation- and polarization-averaged gravitational wave strain, h_c ~ 6e-24 for a 0.6 kHz source at 1 kpc with M_a > 10^-5 M_Sun, exceeds previous estimates that failed to treat equatorward spreading and flux freezing self-consistently. It is concluded that an accreting millisecond pulsar emits a persistent, sinusoidal gravitational wave signal at levels detectable, in principle, by long baseline interferometers after phase-coherent integration, provided that the polar mountain is hydromagnetically stable. Magnetic burial also reduces the magnetic dipole moment, mu, monotonically, implying a novel, observationally testable scaling h_c(mu). The implications for the rotational evolution of (accreting) X-ray and (isolated) radio millisecond pulsars are explored.

The hydromagnetic structure of a neutron star accreting symmetrically at both magnetic poles is calculated as a function of accreted mass, M_a, and polar cap radius,starting from a centered magnetic dipole and evolving through a quasistatic sequence of two-dimensional, Grad-Shafranov equilibria. The calculation is the first to track fully the growth of high-order magnetic multipoles, due to equatorward hydromagnetic spreading, while simultaneously preserving flux freezing and a self-consistent mass-flux distribution. Equilibria are constructed numerically by an iterative scheme and analytically by Green functions. Two key results are obtained, with implications for recycled pulsars. (i) The mass required to significantly reduce the magnetic dipole moment, 10^{-5} Msun, greatly exceeds previous estimates (~ 10^{-10} Msun), which ignored the confining stress exerted by the compressed equatorial magnetic field. (ii) Magnetic bubbles, disconnected from the stellar surface, form in the later stages of accretion (M_a > 10^{-4} Msun).