We discuss, without assuming asymptotic flatness, a gravitational lens for an observer and source that are within a finite distance from a lens object. The proposed lens equation is consistent with the deflection angle of light that is defined for nonasymptotic observer and source by Takizawa et al. [Phys. Rev. D 101, 104032 (2020)] based on the Gauss-Bonnet theorem with using the optical metric. This lens equation, though it is shown to be equivalent to the lens equation by Bozza [Phys. Rev. D 78, 103005 (2008)], is linear in the deflection angle. Therefore, the proposed equation is more convenient for the purpose of doing an iterative analysis. As an explicit example of an asymptotically nonflat spacetime, we consider a static and spherically symmetric solution in Weyl conformal gravity, especially a case that $\gamma$ parameter in the Weyl gravity model is of the order of the inverse of the present Hubble radius. For this case, we examine iterative solutions for the finite-distance lens equation up to the third order. The effect of the Weyl gravity on the lensed image position begins at the third order and it is linear in the impact parameter of light. The deviation of the lensed image position from the general relativistic one is $\sim 10^{-2}$ microarcsecond for the lens and source with a separation angle of $\sim 1$ arcminute, where we consider a cluster of galaxies with $10^{14} M_{\odot}$ at $\sim 1$ Gpc for instance. The deviation becomes $\sim 10^{-1}$ microarcseconds, even if the separation angle is $\sim 10$ arcminutes. Therefore, effects of the Weyl gravity model are negligible in current and near-future observations of gravitational lensing. On the other hand, the general relativistic corrections at the third order $\sim 0.1$ milliarcseconds can be relevant with VLBI observations.

The gravitational deflection angle of light for an observer and source at finite distance from a lens object has been studied by Ishihara et al. [Phys. Rev. D, 94, 084015 (2016)], based on the Gauss-Bonnet theorem with using the optical metric. Their approach to finite-distance cases is limited within an asymptotically flat spacetime. By making several assumptions, we give an interpretation of their definition from the observer's viewpoint: The observer assumes the direction of a hypothetical light emission at the observer position and makes a comparison between the fiducial emission direction and the direction along the real light ray. The angle between the two directions at the observer location can be interpreted as the deflection angle by Ishihara et al. The present interpretation does not require the asymptotic flatness. Motivated by this, we avoid such asymptotic regions to discuss another integral form of the deflection angle of light. This form makes it clear that the proposed deflection angle can be used not only for asymptotically flat spacetimes but also for asymptotically nonflat ones. We examine the proposed deflection angle in two models for the latter case; Kottler (Schwarzschild-de Sitter) solution in general relativity and a spherical solution in Weyl conformal gravity. Effects of finite distance on the light deflection in Weyl conformal gravity result in an extra term in the deflection angle, which may be marginally observable in a certain parameter region. On the other hand, those in Kottler spacetime are beyond reach of the current technology.

GW Ori is a hierarchical triple system which has a rare circumtriple disk. We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of 1.3 mm dust continuum and 12CO J=2-1 molecular gas emission of the disk. For the first time, we identify three dust rings in the disk at ~46, 188, and 338 AU, with estimated dust mass of ~70-250 Earth masses, respectively. To our knowledge, the outer ring in GW Ori is the largest dust ring ever found in protoplanetary disks. We use visibility modelling of dust continuum to show that the disk has misaligned parts and the innermost dust ring is eccentric. The disk misalignment is also suggested by the CO kinematics modelling. We interpret these substructures as evidence of ongoing dynamical interactions between the triple stars and the circumtriple disk.

We present optical observations of type Ia supernova (SN) 2019ein, starting at 2 days after the estimated explosion date. The spectra and the light curves show that SN 2019ein belongs to the High-Velocity (HV) and Bload Line groups with relatively rapid decline in the light curves (Delta m15(B) = 1.36 +- 0.02 mag) and the short rise time (15.37 +- 0.55 days). The Si II 6355 velocity, associated with a photospheric component but not with a detached high-velocity feature, reached ~ 20,000 km s-1 at 12 days before the B-band maximum. The line velocity however decreased very rapidly and smoothly toward the maximum light, where it was ~ 13,000 km s-1 as relatively low among HV SNe. This indicates that the speed of the spectral evolution of HV SNe Ia is correlated not only to the velocity at the maximum light, but also to the light curve decline rate like the case for Normal-Velocity (NV) SNe Ia. Spectral synthesis modeling shows that the outermost layer at > 17,000 km s-1 is well described by the O-Ne-C burning layer extending to at least 25,000 km s-1, and there is no unburnt carbon below 30,000 km s-1; these properties are largely consistent with the delayed detonation scenario, and are shared with the prototypical HV SN 2002bo despite the large difference in Delta m15(B). This structure is strikingly different from that derived for the well-studied NV SN 2011fe. We suggest that the relation between the mass of 56Ni (or Delta m15) and the extent of the O-Ne-C burning layer provides an important constraint on the explosion mechanism(s) of HV and NV SNe.

We report on optical polarimetric observations of two Apollo type near-Earth asteroids, (422699) 2000 PD3 and (3200) Phaethon, and BVRI photometric observations of 2000 PD3 using the 1.6m Pirka telescope in 2017. We derived the geometric albedo of pv = 0.22 +- 0.06 and the color indices (B-V = 0.282 +- 0.072, V-R = 0.198 +- 0.035 and V-I = 0.203 +- 0.022) for 2000 PD3 which are consistent with those of S-type asteroids (including Q-types). The effective diameter of 2000 PD3 was derived as 0.69 +- 0.15 km using our derived geometric albedo. We found that our polarimetric data of Phaethon in 2017 is deviated from the polarimetric profile taken at different epoch of 2016 using the identical instrument setting (Ito et al., 2018). This result suggests that Phaethon would have a regional heterogeneity in grain size and/or albedo on its surface.

By using a method improved with a generalized optical metric, the deflection of light for an observer and source at finite distance from a lens object in a stationary, axisymmetric and asymptotically flat spacetime has been recently discussed [Ono, Ishihara, Asada, Phys. Rev. D 96, 104037 (2017)]. By using this method, in the weak field approximation, we study the deflection angle of light for an observer and source at finite distance from a rotating Teo wormhole, especially by taking account of the contribution from the geodesic curvature of the light ray in a space associated with the generalized optical metric. Our result of the deflection angle of light is compared with a recent work on the same wormhole but limited within the asymptotic source and observer [Jusufi, Ovgun, Phys. Rev. D 97, 024042, (2018)], in which they employ another approach proposed by Werner with using the Nazim's osculating Riemannian construction method via the Randers-Finsler metric. We show that the two different methods give the same result in the asymptotic limit. We obtain also the corrections to the deflection angle due to the finite distance from the rotating wormhole.

Vortices in protoplanetary disks have attracted attention since the discovery of lopsided structures. One of the possible mechanisms for producing vortices is the Rossby Wave Instability (RWI). In our previous work, we have performed detailed linear stability analyses of the RWI with various initial conditions. In this paper, we perform numerical simulations of the vortex formation by the RWI in 2D barotropic disks using the Athena++ code. As initial conditions, we consider axisymmetric disks with a Gaussian surface density bump of various contrasts and half-widths. Perturbations grow as expected from the linear stability analyses in the linear and weakly non-linear regimes. After the saturation, multiple vortices are formed in accordance with the most unstable azimuthal mode and coalesce one after another. In the end, only one quasi-stationary vortex (the RWI vortex) remains, which migrates inward. During the RWI evolution, the axisymmetric component approaches the stable configuration. We find that the axisymmetric component reaches the marginally stable state for the most unstable azimuthal mode at the saturation and the marginally stable state for the m = 1 mode at the final vortex merger. We investigate the structure and evolution of the RWI vortices. We obtain some empirical relations between the properties of the RWI vortices and the initial conditions. Using tracer particle analyses, we find that the RWI vortex can be considered as a physical entity like a large fluid particle. Our results provide a solid theoretical ground for quantitative interpretation of the observed lopsided structures in protoplanetary disks.

By using the Gauss-Bonnet theorem, the bending angle of light in a static, spherically symmetric and asymptotically flat spacetime has been recently discussed, especially by taking account of the finite distance from a lens object to a light source and a receiver [Ishihara, Suzuki, Ono, Asada, Phys. Rev. D 95, 044017 (2017)]. We discuss a possible extension of the method of calculating the bending angle of light to stationary, axisymmetric and asymptotically flat spacetimes. For this purpose, we consider the light rays on the equatorial plane in the axisymmetric spacetime. We introduce a spatial metric to define the bending angle of light in the finite-distance situation. We show that the proposed bending angle of light is coordinate-invariant by using the Gauss-Bonnet theorem. The non-vanishing geodesic curvature of the photon orbit with the spatial metric is caused in gravitomagnetism, even though the light ray in the four-dimensional spacetime follows the null geodesic. Finally, we consider Kerr spacetime as an example in order to examine how the bending angle of light is computed by the present method. The finite-distance correction to the gravitomagnetic deflection angle due to the Sun's spin is around a pico-arcsecond level. The finite-distance corrections for Sgr A$^{\ast}$ also are estimated to be very small. Therefore, the gravitomagnetic finite-distance corrections for these objects are unlikely to be observed with present technology.

We present the early-phase spectra and the light curves of the broad-lined supernova (SN) 2016coi from $t=7$ to $67$ days after the estimated explosion date. This SN was initially reported as a broad-lined Type SN Ic (SN Ic-BL). However, we found that spectra up to $t=12$ days exhibited the He~{\sc i} $\lambda$5876, $\lambda$6678, and $\lambda$7065 absorption lines. We show that the smoothed and blueshifted spectra of normal SNe Ib are remarkably similar to the observed spectrum of SN 2016coi. The line velocities of SN 2016coi were similar to those of SNe Ic-BL and significantly faster than those of SNe Ib. Analyses of the line velocity and light curve suggest that the kinetic energy and the total ejecta mass of SN 2016coi are similar to those of SNe Ic-BL. Together with broad-lined SNe 2009bb and 2012ap for which the detection of He~{\sc i} were also reported, these SNe could be transitional objects between SNe Ic-BL and SNe Ib, and be classified as broad-lined Type `Ib' SNe (SNe `Ib'-BL). Our work demonstrates the diversity of the outermost layer in broad-lined SNe, which should be related to the variety of the evolutionary paths.

Continuing work initiated in an earlier publication [Ishihara, Suzuki, Ono, Kitamura, Asada, Phys. Rev. D {\bf 94}, 084015 (2016) ], we discuss a method of calculating the bending angle of light in a static, spherically symmetric and asymptotically flat spacetime, especially by taking account of the finite distance from a lens object to a light source and a receiver. For this purpose, we use the Gauss-Bonnet theorem to define the bending angle of light, such that the definition can be valid also in the strong deflection limit. Finally, this method is applied to Schwarzschild spacetime in order to discuss also possible observational implications. The proposed corrections for Sgr A$^{\ast}$ for instance are able to amount to $\sim 10^{-5}$ arcseconds for some parameter range, which may be within the capability of near-future astronomy, while also the correction for the Sun in the weak field limit is $\sim 10^{-5}$ arcseconds.

We discuss a possible extension of calculations of the bending angle of light in a static, spherically symmetric and asymptotically flat spacetime to a non-asymptotically flat case. We examine a relation between the bending angle of light and the Gauss-Bonnet theorem by using the optical metric. A correspondence between the deflection angle of light and the surface integral of the Gaussian curvature may allow us to take account of the finite distance from a lens object to a light source and a receiver. Using this relation, we propose a method for calculating the bending angle of light for such cases. Finally, this method is applied to two examples of the non-asymptotically flat spacetimes to suggest finite-distance corrections: Kottler (Schwarzschild-de Sitter) solution to the Einstein equation and an exact solution in Weyl conformal gravity.

We study linear nonradial perturbations and stability of a marginal stable circular orbit (MSCO) such as the innermost stable circular orbit (ISCO) of a test particle in stationary axisymmetric spacetimes which possess a reflection symmetry with respect to the equatorial plane. The proposed approach is applied to Kerr solution and Majumdar-Papapetrou solution to Einstein equation. Finally, we reexamine MSCOs for a modified metric of a rapidly spinning black hole that has been recently proposed by Johannsen and Psaltis [PRD, 83, 124015 (2011)]. We show that, for the Johannsen and Psaltis's model, circular orbits that are stable against radial perturbations for some parameter region become unstable against vertical perturbations. This suggests that the last circular orbit for this model may be larger than the ISCO.

Protoplanetary disks with non-axisymmetric structures have been observed. The Rossby wave instability (RWI) is considered as one of the origins of the non-axisymmetric structures. We perform linear stability analyses of the RWI in barotropic flow using four representative types of the background flow on a wide parameter space. We find that the co-rotation radius is located at the background vortensity minimum with large concavity if the system is marginally stable to the RWI, and this allows us to check the stability against the RWI easily. We newly derive the necessary and sufficient condition for the onset of the RWI in semi-analytic form. We discuss the applicability of the new condition in realistic systems and the physical nature of the RWI.

We study a marginally stable circular orbit (MSCO) such as the innermost stable circular orbit (ISCO) of a timelike geodesic in any spherically symmetric and static spacetime. It turns out that the metric components are separable from the constants of motion along geodesics. We show also that a metric component $g_{rr}$ with a radial coordinate $r$ does not affect MSCOs. This suggests that, as a test of gravity, any ISCO measurement may be put into the same category as gravitational redshift experiments. MSCOs for exact solutions to the Einstein's equation are also mentioned.

In terms of Sturm's theorem, we reexamine a marginal stable circular orbit (MSCO) such as the innermost stable circular orbit (ISCO) of a timelike geodesic in any spherically symmetric and static spacetime. MSCOs for some of exact solutions to the Einstein's equation are discussed. Strum's theorem is explicitly applied to the Kottler (often called Schwarzschild-de Sitter) spacetime. Moreover, we analyze MSCOs for a spherically symmetric, static and vacuum solution in Weyl conformal gravity.

We analytically calculate the marginally stable surface density profile for rotational instability of protoplanetary disks. The derived profile can be utilized for considering the region in a rotating disk where radial pressure gradient force is comparable to the gravitational force, such as an inner edge, steep gaps or bumps and an outer region of the disk. In this paper we especially focus on the rotational instability in the outer region of disks. We find an protoplanetary disk with a surface density profile of similarity solution becomes rotationally unstable at a certain radius, depending on its temperature profile and a mass of the central star. If the temperature is relatively low and the mass of the central star is high, disks have rotationally stable similarity profiles. Otherwise, deviation from the similarity profiles of surface density could be observable, using facilities with high sensitivity, such as ALMA.

We study the three-dimensional (3D) compact U(1) lattice gauge theory coupled with $N$-flavor Higgs fields by means of the Monte Carlo simulations. This model is relevant to multi-component superconductors, antiferromagnetic spin systems in easy plane, inflational cosmology, etc. It is known that there is no phase transition in the N=1 model. For N=2, we found that the system has a second-order phase transition line $\tilde{c}_1(c_2)$ in the $c_2$(gauge coupling)$-c_1$(Higgs coupling) plane, which separates the confinement phase and the Higgs phase. Numerical results suggest that the phase transition belongs to the universality class of the 3D XY model as the previous works by Babaev et al. and Smiseth et al. suggested. For N=3, we found that there exists a critical line similar to that in the N=2 model, but the critical line is separated into two parts; one for $c_2 < c_{2{\rm tc}}=2.4\pm 0.1$ with first-order transitions, and the other for $ c_{2{\rm tc}} < c_2$ with second-order transitions, indicating the existence of a tricritical point. We verified that similar phase diagram appears for the N=4 and N=5 systems. We also studied the case of anistropic Higgs coupling in the N=3 model and found that there appear two second-order phase transitions or a single second-order transition and a crossover depending on the values of the anisotropic Higgs couplings. This result indicates that an "enhancement" of phase transition occurs when multiple phase transitions coincide at a certain point in the parameter space.

We report on a global CCD time-series photometric campaign to decode the pulsations of the nucleus of the planetary nebula NGC1501. The star is hot and hydrogen-deficient, similar to the pre-white-dwarf PG 1159 stars. NGC1501 shows pulsational brightness variations of a few percent with periods ranging from 19 to 87 minutes. The variations are very complex, suggesting a pulsation spectrum that requires a long unbroken time series to resolve. Our CCD photometry of the star covers a two-week period in 1991 November, and used a global network of observatories. We obtained nearly continuous coverage over an interval of one week in the middle of the run. We have identified 10 pulsation periods, ranging from 5235 s down to 1154 s. We find strong evidence that the modes are indeed nonradial g-modes. The ratios of the frequencies of the largest-amplitude modes agree with those expected for modes that are trapped by a density discontinuity in the outer layers. We offer a model for the pulsation spectrum that includes a common period spacing of 22.3 s and a rotation period of 1.17 days; the period spacing allows us to assign a seismological mass of 0.55+/-0.03 Msun.