Bipartite temporal Bell inequalities are similar to the usual Bell inequalities except that, instead of changing the direction of the polariser at each measurement, one changes the time at which the measurement is being performed. By doing so, one is able to test for realism and locality, but relying on position measurements only. This is particularly useful in experimental setups where the momentum direction cannot be probed (such as in cosmology for instance). We study these bipartite temporal Bell inequalities for continuous systems placed in two-mode squeezed states, and find some regions in parameter space where they are indeed violated. We highlight the role played by the rotation angle, which is one of the three parameters characterising a two-mode squeezed state (the other two being the squeezing amplitude and the squeezing angle). In single-time measurements, it only determines the overall phase of the wavefunction and can therefore be discarded, but in multiple-time measurements, its time dynamics becomes relevant and crucially determines when bipartite temporal Bell inequalities can be violated. Our study opens up the possibility of new experimental designs for the observation of Bell inequality violations.

We performed the detailed analysis of the primordial black hole (PBH) formation mechanism in an axionlike curvaton model with a coupling to inflaton. The phase direction of the complex scalar works as a curvaton and produces enough PBHs to explain the black hole binaries ($\sim 30 M_{\odot}$) observed in the LIGO-Virgo Collaboration or PBHs as dark matter (DM) ($\sim 10^{-12} M_{\odot}$). We examine whether our model satisfies the current constraints on the PBH mass spectrum, the curvature perturbation and the secondarily produced gravitational waves. We also take into account ambiguity about the choice of the window functions and effect of the non-Gaussianity.

We revise primordial black hole (PBH) formation in the axion-like curvaton model and investigate whether PBHs formed in this model can be the origin of the gravtitational wave (GW) signals detected by the Advanced LIGO. In this model, small-scale curvature perturbations with large amplitude are generated, which is essential for PBH formation. On the other hand, large curvature perturbations also become a source of primordial GWs by their second-order effects. Severe constraints are imposed on such GWs by pulsar timing array (PTA) experiments. We also check the consistency of the model with these constraints. In this analysis, it is important to take into account the effect of non-Gaussianity, which is generated easily in the curvaton model. We see that, if there are non-Gaussianities, the fixed amount of PBHs can be produced with a smaller amplitude of the primordial power spectrum.

Primordial black holes (PBHs) can be produced by the perturbations that exit the horizon during inflationary phase. While inflation models predict the power spectrum of the perturbations in Fourier space, the PBH abundance depends on the probability distribution function (PDF) of density perturbations in real space. In order to estimate the PBH abundance in a given inflation model, we must relate the power spectrum in Fourier space to the PDF in real space by coarse-graining the perturbations with a window function. However, there are uncertainties on what window function should be used, which could change the relation between the PBH abundance and the power spectrum. This is particularly important in considering PBHs with mass $30 M_\odot$ that account for the LIGO events because the required power spectrum is severely constrained by the observations. In this paper, we investigate how large influence the uncertainties on the choice of a window function have over the power spectrum required for LIGO PBHs. As a result, it is found that the uncertainties significantly affect the prediction for the stochastic gravitational waves (GWs) induced by the second order effect of the perturbations. In particular, the pulsar timing array constraints on the produced GWs could disappear for the real-space top-hat window function.

Continuing the project described by Kato et al. (2009, arXiv:0905.1757), we collected times of superhump maxima for 128 SU UMa-type dwarf novae observed mainly during the 2015-2016 season and characterized these objects. The data have improved the distribution of orbital periods, the relation between the orbital period and the variation of superhumps, the relation between period variations and the rebrightening type in WZ Sge-type objects. Coupled with new measurements of mass ratios using growing stages of superhumps, we now have a clearer and statistically greatly improved evolutionary path near the terminal stage of evolution of cataclysmic variables. Three objects (V452 Cas, KK Tel, ASASSN-15cl) appear to have slowly growing superhumps, which is proposed to reflect the slow growth of the 3:1 resonance near the stability border. ASASSN-15sl, ASASSN-15ux, SDSS J074859.55+312512.6 and CRTS J200331.3-284941 are newly identified eclipsing SU UMa-type (or WZ Sge-type) dwarf novae. ASASSN-15cy has a short (~0.050 d) superhump period and appears to belong to EI Psc-type objects with compact secondaries having an evolved core. ASASSN-15gn, ASASSN-15hn, ASASSN-15kh and ASASSN-16bu are candidate period bouncers with superhump periods longer than 0.06 d. We have newly obtained superhump periods for 79 objects and 13 orbital periods, including periods from early superhumps. In order that the future observations will be more astrophysically beneficial and rewarding to observers, we propose guidelines how to organize observations of various superoutbursts.

A model for the asteroid resurfacing by regolith convection is built to estimate its timescale. In the model, regolith convection is driven by the impact-induced global seismic shaking. The model consists of three steps: (i) intermittent impact of meteoroids, (ii) impact-induced global vibration (seismic shaking), and (iii) vibration-induced regolith convection. In order to assess the feasibility of the resurfacing process driven by regolith convection, we estimate the resurfacing timescale as a function of the size of a target asteroid. In the estimate, a set of parameter values is assumed on the basis of previous works. However, some of them (e.g., seismic quality factor $Q$, seismic efficiency $\eta$, and seismic frequency $f$) are very uncertain. Although these parameter values might depend on asteroid size, we employ the standard values to estimate the representative behavior. To clarify the parameter dependences, we develop an approximated scaling form for the resurfacing timescale. According to the estimated result, we find that the regolith-convection-based resurfacing timescale is shorter than the mean collisional lifetime in most of the parameter uncertainty ranges. These parameter ranges are within those reported by previous works for small asteroids. This means that the regolith convection can be a possible mechanism for the asteroid resurfacing process.

We present analyses to determine the fundamental parameters of the Galaxy based on VLBI astrometry of 52 Galactic maser sources obtained with VERA, VLBA and EVN. We model the Galaxy's structure with a set of parameters including the Galaxy center distance R_0, the angular rotation velocity at the LSR Omega_0, mean peculiar motion of the sources with respect to Galactic rotation (U_src, V_src, W_src), rotation-curve shape index, and the V component of the Solar peculiar motions V_sun. Based on a Markov chain Monte Carlo method, we find that the Galaxy center distance is constrained at a 5% level to be R_0 = 8.05 +/- 0.45 kpc, where the error bar includes both statistical and systematic errors. We also find that the two components of the source peculiar motion U_src and W_src are fairly small compared to the Galactic rotation velocity, being U_src = 1.0 +/- 1.5 km/s and W_src = -1.4 +/- 1.2 km/s. Also, the rotation curve shape is found to be basically flat between Galacto-centric radii of 4 and 13 kpc. On the other hand, we find a linear relation between V_src and V_sun as V_src = V_sun -19 (+/- 2) km/s, suggesting that the value of V_src is fully dependent on the adopted value of V_sun. Regarding the rotation speed in the vicinity of the Sun, we also find a strong correlation between Omega_0 and V_sun. We find that the angular velocity of the Sun, Omega_sun, which is defined as Omega_sun = Omega_0 + V_sun/R_0, can be well constrained with the best estimate of Omega_sun = 31.09 +/- 0.78 km/s/kpc. This corresponds to Theta_0 = 238 +/- 14 km/s if one adopts the above value of R_0 and recent determination of V_sun ~ 12 km/s.

We have performed the astrometry of H2O masers in the Galactic star-forming region Onsala 2 North (ON2N) with the VLBI Exploration of Radio Astrometry.We obtained the trigonometric parallax of 0.261+/-0.009 mas, corresponding to the heliocentric distance of 3.83+/-0.13 kpc. ON2N is expected to be on the solar circle, because its radial velocity with respect to the Local Standard of Rest (LSR) is nearly zero. Using present parallax and proper motions of the masers, the Galactocentric distance of the Sun and the Galactic rotation velocity at the Sun are found to be R_0 = 7.80+/-0.26 kpc and Theta_0 = 213+/-5 km/s, respectively. The ratio of Galactic constants, namely the angular rotation velocity of the LSR can be determined more precisely, and is found to be Omega_0=Theta_0/R_0 = 27.3+/-0.8 km/s/kpc, which is consistent with the recent estimations but different from 25.9 km/s/kpc derived from the recommended values of Theta_0 and R_0 by the International Astronomical Union (1985).

A new computational scheme is developed to study gas accretion from a circumbinary disk. The scheme decomposes the gas velocity into two components one of which denotes the Keplerian rotation and the other of which does the deviation from it. This scheme enables us to solve the centrifugal balance of a gas disk against gravity with better accuracy, since the former inertia force cancels the gravity. It is applied to circumbinary disk rotating around binary of which primary and secondary has mass ratio, 1.4:0.95. The gravity is reduced artificially softened only in small circular regions around the primary and secondary. The radii are 7% of the binary separation and much smaller than those in the previous grid based simulations. 7 Models are constructed to study dependence on the gas temperature and the initial inner radius of the disk. The gas accretion shows both fast and slow time variations while the binary is assumed to have a circular orbit. The time variation is due to oscillation of spiral arms in the circumbinary disk. The masses of primary and secondary disks increase while oscillating appreciably. The mass accretion rate tends to be higher for the primary disk although the secondary disk has a higher accretion rate in certain periods. The primary disk is perturbed intensely by the impact of gas flow so that the outer part is removed. The secondary disk is quiet in most of time on the contrary. Both the primary and secondary disks have traveling spiral waves which transfer angular momentum within them.

We present a distance measurement for the semiregular variable S Crateris (S Crt) based on its annual parallax. With the unique dual beam system of the VLBI Exploration for Radio Astrometry (VERA) telescopes, we measured the absolute proper motion of a water maser spot associated with S Crt, referred to the quasar J1147-0724 located at an angular separation of 1.23$^{\circ}$. In observations spanning nearly two years, we have detected the maser spot at the LSR velocity of 34.7 km s$^{-1}$, for which we measured the annual parallax of 2.33$\pm$0.13 mas corresponding to a distance of 430$^{+25}_{-23}$ pc. This measurement has an accuracy one order of magnitude better than the parallax measurements of HIPPARCOS. The angular distribution and three-dimensional velocity field of maser spots indicate a bipolar outflow with the flow axis along northeast-southwest direction. Using the distance and photospheric temperature, we estimate the stellar radius of S Crt and compare it with those of Mira variables.

We present results of multi-epoch VLBI observations with VERA (VLBI Exploration of Radio Astrometry) of the 22 GHz H$_{2}$O masers associated with a young stellar object (YSO) IRAS 22198+6336 in a dark cloud L1204G. Based on the phase-referencing VLBI astrometry, we derive an annual parallax of IRAS 22198+6336 to be 1.309$\pm$0.047 mas, corresponding to the distance of 764$\pm$27 pc from the Sun. Although the most principal error source of our astrometry is attributed to the internal structure of the maser spots, we successfully reduce the errors in the derived annual parallax by employing the position measurements for all of the 26 detected maser spots. Based on this result, we reanalyze the spectral energy distribution (SED) of IRAS 22198+6336 and find that the bolometric luminosity and total mass of IRAS 22198+6336 are 450$L_{\odot}$ and 7$M_{\odot}$, respectively. These values are consistent with an intermediate-mass YSO deeply embedded in the dense dust core, which has been proposed to be an intermediate-mass counterpart of a low-mass Class 0 source. In addition, we obtain absolute proper motions of the H$_{2}$O masers for the most blue-shifted components. We propose that the collimated jets aligned along the east-west direction are the most plausible explanation for the origin of the detected maser features.

The air fluorescence technique is used to detect ultra-high energy cosmic rays (UHECR), and to estimate their energy. Of fundamental importance is the photon yield due to excitation by electrons, in air of various densities and temperatures. After our previous report, the experiment has been continued using a Sr90 $\beta$ source to study the pressure dependence of photon yields for radiation in nitrogen and dry air. The photon yields in 15 wave bands between 300 nm and 430 nm have been determined. The total photon yield between 300 nm and 406 nm (used in most experiments) in air excited by a 0.85 MeV electron is 3.81+-0.13 (+-13 % systematics) photons per meter at 1013 hPa and 20 $^{\circ}$C. The air density and temperature dependencies of 15 wave bands are given for application to UHECR observations.

In order to detect ultrahigh-energy cosmic rays (UHECR), atmospheric fluorescence light from the trajectory of the extensive air shower may be measured by mirror-photosensor systems. In this type of experiment the photon yield from electrons exciting air of various densities and temperatures is most fundamental information for estimating the primary energy of UHECR. An experiment has been undertaken using a Sr90 $\beta$ source to study the pressure dependence of photon yields, and the life times of the excited states, for radiation in nitrogen and dry air. The photon yield between 300 nm and 406 nm in air excited by 0.85 MeV electrons is 3.73+-0.15 (+-14% systematic) photons per meter at 1000 hPa and 20 $^{\circ}$C. The air density and temperature dependence is given for application to UHECR observations.