Sagittarius A* (Sgr A*) is the central supermassive black hole with the mass $\sim 4\times 10^6 M_{\odot}$ in the Milky Way and stars are orbiting around it. In May 2018, one of the nearest stars to Sgr A* named S0-2/S2 experienced the pericenter passage. The redshift of photons from S0-2 had varied from $4000\>{\rm km\>s^{-1}}$ to $-2000\>{\rm km\>s^{-1}}$ during the pericenter passage, which is within $0.5\>{\rm yr}$. In this paper, we show that this steep variation of the redshift gives a strong constraint on a dark mass distribution inside the orbit of S0-2. By applying a simple $\chi^2$ analysis to the observed redshift, we can easily distinguish between the two models, the point mass model and the point mass plus an extended mass model without the best-fitting parameter search. Our redshift data during the pericenter passage in 2018 with Subaru/IRCS bound the amount of the extended mass inside the orbit of S0-2 less than $0.5\,\%$ ($\sim 2 \times 10^4 M_\odot$) of the mass of Sgr A*. This constraint obtained by our simple analysis is comparable to previous works with the best-fitting parameter search to the motion of S0-2 including the effect of the extended mass. We consider both the power-law and the Plummer models for the dark mass distribution model, but the significant difference between these results is not found.

Searching for space-time variations of the constants of Nature is a promising way to search for new physics beyond General Relativity and the standard model motivated by unification theories and models of dark matter and dark energy. We propose a new way to search for a variation of the fine-structure constant using measurements of late-type evolved giant stars from the S-star cluster orbiting the supermassive black hole in our Galactic Center. A measurement of the difference between distinct absorption lines (with different sensitivity to the fine structure constant) from a star leads to a direct estimate of a variation of the fine structure constant between the star's location and Earth. Using spectroscopic measurements of 5 stars, we obtain a constraint on the relative variation of the fine structure constant below $10^{-5}$. This is the first time a varying constant of Nature is searched for around a black hole and in a high gravitational potential. This analysis shows new ways the monitoring of stars in the Galactic Center can be used to probe fundamental physics.

The star S0-2, orbiting the Galactic central massive black hole candidate Sgr A*, passed its pericenter in May 2018. This event is the first chance to detect the general relativistic (GR) effect of a massive black hole, free from non-gravitational physics. The observable GR evidence in the event is the difference between the GR redshift and the Newtonian redshift of photons coming from S0-2. Within the present observational precision, the 1st post-Newtonian (1PN) GR evidence is detectable. In this paper, we give a theoretical analysis of the time evolution of the 1PN GR evidence, under a presupposition that is different from used in previous papers. Our presupposition is that the GR/Newtonian redshift is always calculated with the parameter values (the mass of Sgr A*, the initial conditions of S0-2, and so on) determined by fitting the GR/Newtonian motion of S0-2 with the observational data. It is then revealed that the difference of the GR redshift and the Newtonian one shows two peaks before and after the pericenter passage. This double-peak-appearance is due to our presupposition, and reduces to a single peak if the same parameter values are used in both GR and Newtonian redshifts as considered in previous papers. In addition to this theoretical discussion, we report our observational data obtained with the Subaru telescope by 2018. The quality and the number of Subaru data in 2018 are not sufficient to confirm the detection of the double-peak-appearance.

General Relativity predicts that a star passing close to a supermassive black hole should exhibit a relativistic redshift. We test this using observations of the Galactic center star S0-2. We combine existing spectroscopic and astrometric measurements from 1995-2017, which cover S0-2's 16-year orbit, with measurements in 2018 March to September which cover three events during its closest approach to the black hole. We detect the combination of special relativistic- and gravitational-redshift, quantified using a redshift parameter, $\Upsilon$. Our result, $\Upsilon=0.88 \pm 0.17$, is consistent with General Relativity ($\Upsilon=1$) and excludes a Newtonian model ($\Upsilon=0$ ) with a statistical significance of 5 $\sigma$.

During the next closest approach of the orbiting star S2/S0-2 to the Galactic supermassive black hole (SMBH), it is estimated that RV uncertainties of ~ 10 km/s allow us to detect post-Newtonian effects throughout 2018. To evaluate an achievable uncertainty in RV and its stability, we have carried out near-infrared, high resolution (R ~ 20,000) spectroscopic monitoring observations of S2 using the Subaru telescope and the near-infrared spectrograph IRCS from 2014 to 2016. The Br-gamma absorption lines are used to determine the RVs of S2. The RVs we obtained are 497 km/s, 877 km/s, and 1108 km/s in 2014, 2015, and 2016, respectively. The statistical uncertainties are derived using the jackknife analysis. The wavelength calibrations in our three-year monitoring are stable: short-term (hours to days) uncertainties in RVs are < 0.5 km/s, and a long-term (three years) uncertainty is 1.2 km/s. The uncertainties from different smoothing parameter, and from the partial exclusion of the spectra, are found to be a few km/s. The final results using the Br-gamma line are 497 +- 17 (stat.) +- 3 (sys.) km/s in 2014, 877 +- 15 (stat.) +- 4 (sys.) km/s in 2015, and 1108 +- 12 (stat.) +- 4 (sys.) km/s in 2016. When we use two He I lines at 2.113\mum in addition to Br-gamma, the mean RVs are 513 km/s and 1114 km/s for 2014 and 2016, respectively. The standard errors of the mean are 16.2 km/s (2014) and 5.4 km/s (2016), confirming the reliability of our measurements. The difference between the RVs estimated by Newtonian mechanics and general relativity will reach about 200 km/s near the next pericenter passage in 2018. Therefore our RV uncertainties of 13 - 17 km/s with Subaru enable us to detect the general relativistic effects in the RV measurements with more than 10 sigma in 2018.

Stationary pulsar magnetospheres in the force-free system are governed by the pulsar equation. In 1999, Contopoulos, Kazanas, and Fendt (hereafter CKF) numerically solved the pulsar equation and obtained a pulsar magnetosphere model called the CKF solution that has both closed and open magnetic field lines. The CKF solution is a successful solution, but it contains a poloidal current sheet that flows along the last open field line. This current sheet is artificially added to make the current system closed. In this paper, we suggest an alternative method to solve the pulsar equation and construct pulsar magnetosphere models without a current sheet. In our method, the pulsar equation is decomposed into Ampere's law and the force-free condition. We numerically solve these equations simultaneously with a fixed poloidal current. As a result, we obtain a pulsar magnetosphere model without a current sheet, which is similar to the CKF solution near the neutron star and has a jet-like structure at a distance along the pole. In addition, we discuss physical properties of the model and find that the force-free condition breaks down in a vicinity of the light cylinder due to dissipation that is included implicitly in the numerical method.

We numerically construct an one-parameter family of initial data of an expanding inhomogeneous universe model which is composed of regularly aligned black holes with an identical mass. They are initial data for vacuum solutions of the Einstein equations. We call this universe model the "black hole universe" and analyze the structure of these initial data. We study the relation between the mean expansion rate of the 3-space, which corresponds to the Hubble parameter, and the mass density of black holes. The result implies that the same relation as that of the Einstein-de Sitter universe is realized in the limit of the large separation between neighboring black holes. The applicability of the cosmological Newtonian $N$-body simulation to the dark matter composed of black holes is also discussed. The deviation of the spatial metric of the cosmological Newtonian $N$-body system from that of the black hole universe is found to be smaller than about 1% in a region distant from the particles, if the separation length between neighboring particles is 20 times larger than their gravitational radius. By contrast, the deviation of the square of the Hubble parameter of the cosmological Newtonian $N$-body system from that of the black hole universe is about 20% for the same separation length.

It is known that the Meissner-like effect is seen in a magnetosphere without an electric current in black hole spacetime: no non-monopole component of magnetic flux penetrates the event horizon if the black hole is extreme. In this paper, in order to see how an electric current affects the Meissner-like effect, we study a force-free electromagnetic system in a static and spherically symmetric extreme black hole spacetime. By assuming that the rotational angular velocity of the magnetic field is very small, we construct a perturbative solution for the Grad-Shafranov equation, which is the basic equation to determine a stationary, axisymmetric electromagnetic field with a force-free electric current. Our perturbation analysis reveals that, if an electric current exists, higher multipole components may be superposed upon the monopole component on the event horizon, even if the black hole is extreme.