This is the third paper of a series of our works on the self-similar orbit-averaged Fokker-Planck (OAFP) equation. The first paper provided an accurate spectral solution of the equation for isotropic pre-collapse star clusters and the second detailed the physical feature of the model. Based on the works, the present work applies the solution to the observed structural profiles of Galactic globular clusters. For fitting to the profiles, the most fundamental (quasi-)stationary model, the King model, and the variants have shown successful results while they can not apply to core-collapsing and core-collapsed clusters at the late stage of the relaxation evolution. We propose an energy-truncated self-similar OAFP model that can apply to clusters at both the early and late stages of the evolution. This new model fits the structural profiles of at least half of Galactic globular clusters while it also applies to core-collapsed stars with resolved cores. As a main result, we provide the completion rate of core collapse against concentration for the clusters. Also, we show our new model can apply to the globular clusters even in a broad range of radii (0.01$\sim$10 arcminutes). However, since our model includes polytrope (elongated outer halo), the tidal radius of the model becomes unrealistically large for some clusters. To avoid the issue, we also propose an approximated form of the new model. Lastly, we report that Milky Way globular clusters with low concentrations have the same spatial structures as stellar polytropes and discuss whether such polytropic cluster is a reasonable concept.

This is the second paper of a series of our works on the self-similar orbit-averaged Fokker-Planck (OAFP) equation and details physical properties of isotropic pre-collapse solution. The fundamental core collapse process at the late stage of relaxation evolution of spherical star clusters can be described by the self-similar OAFP equation. The accurate spectral solution was found recently in the first paper. The present work details the thermodynamic aspects of the model based on the stellar DF obtained from the solution. Our calculation shows the following local properties (i) Equation of state is $p=1.0\rho/\chi_\text{esc}$ in the core where $p$ is the pressure, $\rho$ the density and $\chi_\text{esc}$ the scaled escape energy, while it is $p=0.5\rho^{1.1}/\chi_\text{esc}$ at large radii. (ii) If we consider the center of the core is polytropic, the polytropic index is 177. Also, as a global property we construct caloric curves of the model to discuss the heat capacity together with Virial. Special focus is the cause of negative heat capacity of the core; the well-relaxed core can be directly compared to the isothermal sphere with positive heat capacity. Comparing our results to the previous works, we conclude, in the self-similar evolution, the negative heat capacity in the core holds due to collisionless and high-temperature stars that experience a rapid change in mean field potential through stellar- and heat- flows, rather than due to the (quasi-)isolation of the core from surroundings.

This is the first paper of a series of our works on the self-similar orbit-averaged Fokker-Planck (OAFP) equation for distribution function of stars in dense isotropic star clusters. At the late stage of relaxation evolution of the clusters, standard stellar dynamics predicts that the clusters evolve in a self-similar fashion forming collapsing cores. However, the corresponding mathematical model, the self-similar OAFP equation, has never been solved on the whole energy domain $(-1< E < 0)$. The existing works based on kinds of finite difference methods provide solutions only on the truncated domain $-1< E<-0.2$. To broaden the range of the truncated domain, the present work resorts to a (highly accurate and efficient) Gauss-Chebyshev pseudo-spectral method. We provide a spectral solution, whose number of significant figures is four, on the whole domain. Also, the solution can reduce to a semi-analytical form whose degree of polynomials is only eighteen holding three significant figures. We also provide the new eigenvalues; $c_{1}=9.0925\times10^{-4}$, $c_{2}=1.1118\times10^{-4}$, $c_{3}=7.1975\times10^{-2}$ and $c_{4}=3.303\times10^{-2}$, corresponding to the core collapse rate $\xi=3.64\times10^{-3}$, scaled escape energy $\chi_\text{esc}=13.881$ and power-law exponent $\alpha=2.2305$. Since the solution on the whole domain is unstable against degree of Chebyshev polynomials, we also provide spectral solutions on truncated domains ( $-1< E<E_\text{max}$, where $-0.35<E_\text{max}<-0.03$) to explain how to handle the instability. By reformulating the OAFP equation in several ways, we improve the accuracy of the spectral solution and reproduce an existing self-similar solution, which infers that existing solutions have only one significant figure at best.

In this paper we report on the characterization of the Hamamatsu VUV4 (S/N: S13370-6152) Vacuum Ultra-Violet (VUV) sensitive Silicon Photo-Multipliers (SiPMs) as part of the development of a solution for the detection of liquid xenon scintillation light for the nEXO experiment. Various SiPM features, such as: dark noise, gain, correlated avalanches, direct crosstalk and Photon Detection Efficiency (PDE) were measured in a dedicated setup at TRIUMF. SiPMs were characterized in the range $163 \text{ } \text{K} \leq \text{T}\leq 233 \text{ } \text{K}$. At an over voltage of $3.1\pm0.2$ V and at $\text{T}=163 \text{ }\text{K}$ we report a number of Correlated Avalanches (CAs) per pulse in the $1 \upmu\text{s}$ interval following the trigger pulse of $0.161\pm0.005$. At the same settings the Dark-Noise (DN) rate is $0.137\pm0.002 \text{ Hz/mm}^{2}$. Both the number of CAs and the DN rate are within nEXO specifications. The PDE of the Hamamatsu VUV4 was measured for two different devices at $\text{T}=233 \text{ }\text{K}$ for a mean wavelength of $189\pm7\text{ nm}$. At $3.6\pm0.2$ V and $3.5\pm0.2$ V of over voltage we report a PDE of $13.4\pm2.6\text{ }\%$ and $11\pm2\%$, corresponding to a saturation PDE of $14.8\pm2.8\text{ }\%$ and $12.2\pm2.3\%$, respectively. Both values are well below the $24\text{ }\%$ saturation PDE advertised by Hamamatsu. More generally, the second device tested at $3.5\pm0.2$ V of over voltage is below the nEXO PDE requirement. The first one instead yields a PDE that is marginally close to meeting the nEXO specifications. This suggests that with modest improvements the Hamamatsu VUV4 MPPCs could be considered as an alternative to the FBK-LF SiPMs for the final design of the nEXO detector.

We present a 190 ks observation of the Galactic supernova remnant (SNR) G306.3-0.9 with Suzaku. To study ejecta properties of this possible Type Ia SNR, the absolute energy scale at the Fe-K band was calibrated to a level of uncertainty less than 10 eV by a cross-calibration with the Hitomi microcalorimeter using the Perseus cluster spectra. This enabled us for the first time to accurately determine the ionization state of the Fe K$\alpha$ line of this SNR. The ionization timescale ($\tau$) of the Fe ejecta was measured to be $\log_{10} \tau$ (cm$^{-3}$ s) $=10.24\pm0.03$, significantly smaller than previous measurements. Marginally detected K$\alpha$ lines of Cr and Mn have consistent ionization timescales with Fe. The global spectrum was well fitted with shocked interstellar matter (ISM) and at least two ejecta components with different ionization timescales for Fe and intermediate mass elements (IME) such as S and Ar. One plausible interpretation of the one-order-of-magnitude shorter timescale of Fe than that of IME ($\log_{10} \tau = 11.17\pm0.07$) is a chemically stratified structure of ejecta. By comparing the X-ray absorption column to the HI distribution decomposed along the line of sight, we refined the distance to be $\sim$20 kpc. The large ISM-to-ejecta shocked mass ratio of $\sim$100 and dynamical timescale of $\sim$6 kyr place the SNR in the late Sedov phase. These properties are consistent with a stratified ejecta structure that has survived the mixing processes expected in an evolved supernova remnant.

The secular evolution of a collisional star cluster of $N$-'point' stars have been conventionally discussed based on cumulative two-body relaxation process. The relaxation process requires a cut-off on the range of two-body encounter between stars in physical space and the relaxation time is characterized by Coulomb logarithm $\ln[N]$; the conventional cut-off on the encounter distance in the literature gives "dominant" effect. In addition, incorrect cut-offs exposed a mathematical "infinite-density" problem in the late stage of core-collapse. The present paper shows these are merely the results due to incorrect cut-off process. If one correctly constrains the cut-off on interaction range between stars based on truncated BBGKY hierarchy, one must introduce a self-consistent 'truncated' Newtonian mean-field (m.f.) acceleration of star. The present paper shows the effect of total number on the structure of finite star cluster (in which stars undergo only two-body encounters) for the first time after establishing a mathematical formulation of a generalized Landau kinetic equation that includes the cut-off effect on both of collision term and m.f. potential by employing a BBGKY hierarchy for truncated DF stars. The cut-off effect increases the typical relaxation time by a few of percentage, which means the effect of cut-off itself is "non-dominant" on the relaxation time. On the other hand, the cut-off on m.f. acceleration is necessary to avoid "infinite-density" problem at the center of the system; the effect of total number $N$ on density profile and m.f. acceleration are shown by applying the truncated-DF BBGKY hierarchy to a toy model (a quasi-static modified Hubble density profile) for a core-halo structure of a star cluster at the late stage of evolution.

In secular evolution of finite dense star clusters, the statistical acceleration of stars is the most essential non-collective relaxation process whose effect can be mathematically modeled by Kandrup's generalized-Landau (g-Landau) kinetic equation for distribution function of stars. Understanding of the g-Landau equation is of significance in finite system since only the equation can correctly define the total- energy and - number of stars in phase spaces in case the effect of gravitational polarization is neglected. The present paper shows a kinetic formulation of an orbit-averaged generalized Landau equation in action-angle spaces beginning with BBGKY hierarchy and shows the conservation laws, $H$-theorem and anti-normalization condition. Furthermore, the orbit-averaged g-Landau equation is rewritten for anisotropic spherical system. It is shown that the statistical acceleration can be replaced by typical acceleration of star for the relaxation process in secular evolution of any anisotropic spherical systems. The derivation of the equations is made by generalizing the formulation for the inhomogeneous Landau equation done in (Polyachenko 1982) to the g-Landau equation.

Fundamental relaxation processes in the secular evolution of a collisional star cluster of $N$-'point' stars have been conventionally discussed based on either of collision kinetic theory (for strong two-body encounters) and wave one (for statistical acceleration and gravitational polarization). If combining the both theories together, one must introduce a self-consistent 'truncated' Newtonian mean-field (m.f.) acceleration of star at position $r$ and time $t$ due to a phase-space distribution function $f\left(r', p',t\right)$ for stars $A^{\triangle}(r,t)=-Gm\left(1-\frac{1}{N}\right)\int_{\mid r-r' \mid > \triangle}\frac{r-r'}{\mid r-r' \mid^{3}} f\left(r',p',t\right)\text{d}^{3}{r'}\text{d}^{3}{p'},$ where $G$ is the gravitational constant and $m$ the mass of stars. The lower limit $\triangle$ of the distance between two stars is order of the Landau distance. The truncated m.f. acceleration is a necessary consequence due to the strong encounters and m.f. acceleration being not able to 'coexist' at specific distance between stars. The present paper aims at initiating a star-cluster convergent kinetic theory to self-consistently derive kinetic equations of star clusters, mathematically non-divergent in distance- and wavenumber- spaces based on the truncated m.f. acceleration, correct at time scales of the secular evolution. This will be achieved by focusing on mathematical formulations of the Kandrup's generalised-Landau equation including the effect of the strong encounters and by extending the Grad's truncated distribution function and Klimontovich's theory of non-ideal systems.

Motivated by recent detection of transiting high-density super-Earths, we explore the detectability of hot rocky super-Earths orbiting very close to their host stars. In the environment hot enough for their rocky surfaces to be molten, they would have the atmosphere composed of gas species from the magma oceans. In this study, we investigate the radiative properties of the atmosphere that is in the gas/melt equilibrium with the underlying magma ocean. Our equilibrium calculations yield Na, K, Fe, Si, SiO, O, and O$_2$ as the major atmospheric species. We compile the radiative-absorption line data of those species available in literature, and calculate their absorption opacities in the wavelength region of 0.1--100~$\mathrm{\mu m}$. Using them, we integrate the thermal structure of the atmosphere. Then, we find that thermal inversion occurs in the atmosphere because of the UV absorption by SiO. In addition, we calculate the ratio of the planetary to stellar emission fluxes during secondary eclipse, and find prominent emission features induced by SiO at 4~$\mathrm{\mu m}$ detectable by Spitzer, and those at 10 and 100~$\mathrm{\mu m}$ detectable by near-future space telescopes.

Continuing the project described by Kato et al. (2009, PASJ, 61, S395, arXiv:0905.1757), we collected times of superhump maxima for 56 SU UMa-type dwarf novae mainly observed during the 2013-2014 season and characterized these objects. We detected negative superhumps in VW Hyi and indicated that the low number of normal outbursts in some supercycle can be interpreted as a result of the disk tilt. This finding, combined with the Kepler observation of V1504 Cyg and V344 Lyr, suggests that the disk tilt is responsible for modulating the outburst pattern in SU UMa-type dwarf novae. We also studied the deeply eclipsing WZ Sge-type dwarf nova MASTER OT J005740.99+443101.5 and found evidence of a sharp eclipse during the phase of early superhumps. The profile can be reproduced by a combination of the eclipse of the axisymmetric disk and the uneclipsed light source of early superhumps. This finding confirms the lack of evince of a greatly enhanced hot spot during the early stage of WZ Sge-type outburst. We detected growing (stage A) superhumps in MN Dra and give a suggestion that some of SU UMa-type dwarf novae situated near the critical condition of tidal instability may show long-lasting stage A superhumps. The large negative period derivatives reported in such systems can be understood a result of the combination of stage A and B superhumps. The WZ Sge-type dwarf novae AL Com and ASASSN-13ck showed a long-lasting (plateau-type) rebrightening. In the early phase of the rebrightening, both objects showed a precursor-like outburst, suggesting that the long-lasting rebrightening is triggered by a precursor outburst.

Violent variation of transit depths and an ingress-egress asymmetry of the transit light curve discovered in KIC 12557548 have been interpreted as evidences of a catastrophic evaporation of atmosphere with dust (M_p gtrsim 1 M_oplus/Gyr) from a close-in small planet. To explore what drives the anomalous atmospheric escape, we perform time-series analysis of the transit depth variation of Kepler archival data for ~ 3.5 yr. We find a ~ 30% periodic variation of the transit depth with P1 = 22.83 pm 0.21 days, which is within the error of the rotation period of the host star estimated using the light curve modulation, Prot = 22.91 pm 0.24 days. We interpret the results as evidence that the atmospheric escape of KIC 12557548b correlates with stellar activity. We consider possible scenarios that account for both the mass loss rate and the correlation with stellar activity. X-ray and ultraviolet (XUV)-driven evaporation is possible if one accepts a relatively high XUV flux and a high efficiency for converting the input energy to the kinetic energy of the atmosphere. Star-planet magnetic interaction is another possible scenario though huge uncertainty remains for the mass loss rate.

We develop the general reconstruction scheme in two scalar model. The quintom-like theory which may describe (different) non-singular Little Rip or de Sitter cosmology is reconstructed. (In)stability of such dark energy cosmologies as well as the flow to fixed points is studied. The stability of Little Rip universe which leads to dissolution of bound objects sometime in future indicates that no classical transition to de Sitter space occurs.

We have undertaken a search for the infrared emission from the intracluster dust in the Coma cluster of galaxies by the Multiband Imaging Photometer for Spitzer. Our observations yield the deepest mid and far-infrared images of a galaxy cluster ever achieved. In each of the three bands, we have not detected a signature of the central excess component in contrast to the previous report on the detection by Infrared Space Observatory (ISO). We still find that the brightness ratio between 70 and 160 microns shows a marginal sign of the central excess, in qualitative agreement with the ISO result. Our analysis suggests that the excess ratio is more likely due to faint infrared sources lying on fluctuating cirrus foreground. Our observations yield the 2 sigma upper limits on the excess emission within 100 kpc of the cluster center as 5 x 10^-3 MJy/sr, 6 x 10^-2 MJy/sr, and 7 x 10^-2 MJy/sr, at 24, 70, and 160 microns, respectively. These values are in agreement with those found in other galaxy clusters and suggest that dust is deficient near the cluster center by more than 3 orders of magnitude compared to the interstellar medium.