We report on the core-collapse supernova simulation we conducted for a 11.2M progenitor model in three-dimensional space up to 20 ms after bounce, using a radiation hydrodynamics code with full Boltzmann neutrino transport. We solve the six-dimensional Boltzmann equations for three neutrino species and the three-dimensional compressible Euler equations with Furusawa and Togashi's nuclear equation of state. We focus on the prompt convection at ~10 ms after bounce and investigate how neutrinos are transported in the convective matter. We apply a new analysis based on the eigenvalues and eigenvectors of the Eddington tensor and make a comparison between the Boltzmann transport results and the M1 closure approximation in the transition regime between the optically thick and thin limits. We visualize the eigenvalues and eigenvectors using an ellipsoid, in which each principal axis is parallel to one of the eigenvectors and has a length proportional to the corresponding eigenvalue. This approach enables us to understand the difference between the Eddington tensor derived directly from the Boltzmann simulation and the one given by the M1 prescription from a new perspective. We find that the longest principal axis of the ellipsoid is almost always nearly parallel to the energy flux in the M1 closure approximation whereas in the Boltzmann simulation it becomes perpendicular in some transition regions, where the mean free path is ~0.1 times the radius. In three spatial dimensions, the convective motions make it difficult to predict where this happens and possibly improves the closure relation there.

Using the Boltzmann-radiation-hydrodynamics code which solves the Boltzmann equation for the neutrino transport, we present the results of the simulations with the nuclear equations of state (EOSs) of Lattimer-Swesty (LS) and Furusawa-Shen (FS). We extend the simulation time of the LS model and conducted thorough investigations, though the previous paper reported some of the results briefly. Only the LS model shows the shock revival. This seems to originate from the nuclear composition: the different nuclear composition results in the different energy loss by the photodissociation and hence the different strength of the prompt convection and the later neutrino-driven convection. The proto-neutron star seen in the FS model is more compact than that in the LS model since the existence of multi-nuclear species softens the EOS. For the behavior of neutrinos, we examined the flux and the Eddington tensor of neutrinos. In the optically thick region, the diffusion of neutrinos and the dragging by the matter motion determine the flux. In the optically thin region, the free-streaming determines it. The Eddington tensor is compared with that obtained from the M1-closure relation. The M1-closure scheme overestimates the contribution from the velocity-dependent terms in the semi-transparent region.

Neutrinos are believed to have a key role in the explosion mechanism of core-collapse supernovae as they carry most of the energy released by the gravitational collapse of a massive star. If their flavor is converted fast inside the neutrino sphere, the supernova explosion may be influenced. This paper is reporting the results of the extended work of our previous paper. We perform a thorough survey of the ELN crossing in one of our self-consistent, realistic Boltzmann simulations in two spatial dimensions under axisymmetry for the existence of the crossings between $\nu_e$ and $\bar\nu_e$ angular distributions, or the electron lepton number (ELN) crossing. We report for the first time the positive detections deep inside the core of the massive star in the vicinity of neutrino sphere at $r$ $\approx$ 16 - 21 km. We find that low values of the electron fraction $Y_e$ produced by convective motions together with the appearance of light elements are critically important to give rise to the ELN crossing by enhancing the chemical potential difference between proton and neutron, and hence by mitigating the Fermi-degeneracy of $\nu_e$. Since the region of positive detection are sustained and, in fact, expanding with time, it may have an impact on the explosion of core-collapse supernovae, observational neutrino astronomy and nucleosynthesis of heavy nuclei.

With the Boltzmann-radiation-hydrodynamics code, which we have developed to solve numerically the Boltzmann equations for neutrino transfer, the Newtonian hydrodynamics equations, and the Newtonian self-gravity simultaneously and consistently, we simulate the collapse of a rotating core of the progenitor with a zero-age-main-sequence mass of $11.2\,M_\odot$ and a shelluler rotation of $1\,{\rm rad\,s^{-1}}$ at the center. We pay particular attention in this paper to the neutrino distribution in phase space, which is affected by the rotation. By solving the Boltzmann equations directly, we can assess the rotation-induced distortion of the angular distribution in momentum space, which gives rise to the rotational component of the neutrino flux. We compare the Eddington tensors calculated both from the raw data and from the M1-closure approximation. We demonstrate that the Eddington tensor is determined by complicated interplays of the fluid velocity and the neutrino interactions and that the M1-closure, which assumes that the Eddington factor is determined by the flux factor, fails to fully capture this aspect, especially in the vicinity of the shock. We find that the error in the Eddington factor reaches $\sim 20\%$ in our simulation. This is due not to the resolution but to the different dependence of the Eddington and flux factors on the angular profile of the neutrino distribution function, and hence modification to the closure relation is needed.

We present the first results of our spatially axisymmetric core-collapse supernova simulations with full Boltzmann neutrino transport, which amount to a time-dependent 5-dimensional (2 in space and 3 in momentum space) problem in fact. Special relativistic effects are fully taken into account with a two-energy-grid technique. We performed two simulations for a progenitor of 11.2M, employing different nuclear equations-of-state (EOS's): Lattimer and Swesty's EOS with the incompressibility of K = 220MeV (LS EOS) and Furusawa's EOS based on the relativistic mean field theory with the TM1 parameter set (FS EOS). In the LS EOS the shock wave reaches ~700km at 300ms after bounce and is still expanding whereas in the FS EOS it stalled at ~200km and has started to recede by the same time. This seems to be due to more vigorous turbulent motions in the former during the entire post-bounce phase, which leads to higher neutrino-heating efficiency in the neutrino-driven convection. We also look into the neutrino distributions in momentum space, which is the advantage of the Boltzmann transport over other approximate methods. We find non-axisymmetric angular distributions with respect to the local radial direction, which also generate off-diagonal components of the Eddington tensor. We find that the r {\theta}-component reaches ~10% of the dominant rr-component and, more importantly, it dictates the evolution of lateral neutrino fluxes, dominating over the {\theta}{\theta}-component, in the semi-transparent region. These data will be useful to further test and possibly improve the prescriptions used in the approximate methods.

We investigate the criterion for the acoustic mechanism to work successfully in core-collapse supernovae. The acoustic mechanism is an alternative to the neutrino-heating mechanism. It was proposed by Burrows et al., who claimed that acoustic waves emitted by $g$-mode oscillations in proto-neutron stars (PNS) energize a stalled shock wave and eventually induce an explosion. Previous works mainly studied to which extent the $g$-modes are excited in the PNS. In this paper, on the other hand, we investigate how strong the acoustic wave needs to be if it were to revive a stalled shock wave. By adding the acoustic power as a new axis, we draw a critical surface, an extension of the critical curve commonly employed in the context of neutrino heating. We perform both 1D and 2D parametrized simulations, in which we inject acoustic waves from the inner boundary. In order to quantify the power of acoustic waves, we use the extended Myers theory to take neutrino reactions into proper account. We find for the 1D simulations that rather large acoustic powers are required to relaunch the shock wave, since the additional heating provided by the secondary shocks developed from acoustic waves is partially canceled by the neutrino cooling that is also enhanced. In 2D, the required acoustic powers are consistent with those of Burrows et al. Our results seem to imply, however, that it is the sum of neutrino heating and acoustic powers that matters for shock revival.

We perform a set of cosmological simulations of structure formation in a mixed dark matter (MDM) model. Our model is motivated by the recently identified $3.5\,{\rm keV}$ X-ray line, which can be explained by the decay of non-resonantly produced sterile neutrinos accounting for $20$-$60\%$ of the dark matter in the Universe. These non-resonantly produced sterile neutrinos have a sizable free-streaming length and hence behave effectively as warm dark matter (WDM). Assuming the rest of dark matter is composed of some cold dark matter (CDM) particles, we follow the coevolution of a mixed WDM plus CDM cosmology. Specifically, we consider the models with the warm component fraction of $r_{\rm warm}=0.25$ and $0.50$. Our MDM models predict that the comoving Jeans length at the matter-radiation equality is close to that of the thermally produced warm dark matter model with particle mass $m_{\rm WDM}=2.4\,{\rm keV}$, but the suppression in the fluctuation power spectrum is weaker. We perform large $N$-body simulations to study the structure of non-linear dark halos in the MDM models. The abundance of substructure is significantly reduced in the MDM models, and hence the so-called small-scale crisis is mitigated. The cumulative maximum circular velocity function (CVF) of at least one halo in the MDM models is in good agreement with the CVFs of the observed satellites in the Milky Way and the Andromeda galaxy. We argue that the MDM models open an interesting possibility to reconcile the reported $3.5\,{\rm keV}$ line and the internal structure of galaxies.

We have observed the lightcurves of 13 V-type asteroids ((1933) Tinchen, (2011) Veteraniya, (2508) Alupka, (3657) Ermolova, (3900) Knezevic, (4005) Dyagilev, (4383) Suruga, (4434) Nikulin, (4796) Lewis, (6331) 1992 $\mathrm{FZ_{1}}$, (8645) 1998 TN, (10285) Renemichelsen, and (10320) Reiland). Using these observations we determined the rotational rates of the asteroids, with the exception of Nikulin and Renemichelsen. The distribution of rotational rates of 59 V-type asteroids in the inner main belt, including 29 members of the Vesta family that are regarded as ejecta from the asteroid (4) Vesta, is inconsistent with the best-fit Maxwellian distribution. This inconsistency may be due to the effect of thermal radiation Yarkovsky--O'Keefe--Radzievskii--Paddack (YORP) torques, and implies that the collision event that formed V-type asteroids is sub-billion to several billion years in age.

We present new Ryle Telescope (RT) observations of the Sunyaev Zel'dovich (SZ) decrement from the cluster Abell 773. The field contains a number of faint radio sources that required careful subtraction. We use ASCA observations to measure the gas temperature and a ROSAT HRI image to model the gas distribution. Normalising the gas distribution to fit the RT visibilities returns a value of H_0 of 77 (+19,-15) km/s/Mpc (1-sigma errors) for an Einstein-de-Sitter universe, or 85 (+20,-17) km/s/Mpc for a flat model with Omega_Lambda = 0.7. The errors quoted include estimates of the effects of the principal errors: noise in the SZ measurement, gas temperature uncertainty, and line-of sight depth uncertainty.

We have observed the Sunyaev-Zel'dovich effect in a sample of five moderate-redshift clusters with the Ryle Telescope, and used them in conjunction with X-ray imaging and spectral data from ROSAT and ASCA to measure the Hubble constant. This sample was chosen with a strict X-ray flux limit using both the BCS and NORAS cluster catalogues to be well above the surface-brightness limit of the ROSAT All-Sky Survey, and hence to be unbiased with respect to the orientation of the cluster. This controls the major potential systematic effect in the SZ/X-ray method of measureing H0. Taking the weighted geometric mean of the results and including the main sources of random error, namely the noise in the SZ measurement, the uncertainty in the X-ray temperatures and the unknown ellipticity of the clusters, we find H0 = 59 +8/-7 km/s/Mpc assuming a standard CDM model with Omega_M = 1.0, Omega_Lambda = 0.0, or H0 = 65 +8/-7 km/s/Mpc if Omega_M = 0.3, Omega_Lambda = 0.7.