Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a comparison with observational data. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a mm VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects. Based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed case studies to assess the attainable image quality with the current and future EHT array for different weather conditions. The results show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of the input models can be recovered robustly after performing calibration steps. With the planned addition of new stations to the EHT array, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images.

We investigate the second-order gravitational scalar perturbations for a barotropic fluid. We derive the effective energy-momentum tensor described by the quadratic terms of the gravitational and the matter perturbations. We show that the second-order effective energy-momentum tensor is gauge dependent. We impose three gauge conditions (longitudinal, spatially-flat, and comoving gauges) for dust and radiation. The resulting energy-momentum tensor is described only by a gauge invariant variable, but the functional form depends on the gauge choice. In the matter-dominated epoch with dust-like fluid background, the second-order effective energy density and pressure of the perturbations evolve as 1/a^2 in all three gauge choices, like the curvature density of the Universe, but they do not provide the correct equation of state. The value of this parameter depends also on the gauge choice. In the radiation-dominated epoch, the perturbations in the short-wave limit behave in the same way as the radiation-like fluid in the longitudinal and the spatially-flat gauges. However, they behave in a different way in the comoving gauge. As a whole, we conclude that the second-order effective energy-momentum tensor of the scalar perturbation is strictly gauge dependent.

We study the kinematics of the M87 jet using the first year data of the KVN and VERA Array (KaVA) large program, which has densely monitored the jet at 22 and 43 GHz since 2016. We find that the apparent jet speeds generally increase from $\approx0.3c$ at $\approx0.5$ mas from the jet base to $\approx2.7c$ at $\approx20$ mas, indicating that the jet is accelerated from subluminal to superluminal speeds on these scales. We perform a complementary jet kinematic analysis by using archival Very Long Baseline Array monitoring data observed in $2005-2009$ at 1.7 GHz and find that the jet is moving at relativistic speeds up to $\approx5.8c$ at distances of $200-410$ mas. We combine the two kinematic results and find that the jet is gradually accelerated over a broad distance range that coincides with the jet collimation zone, implying that conversion of Poynting flux to kinetic energy flux takes place. If the jet emission consists of a single streamline, the observed trend of jet acceleration ($\Gamma\propto z^{0.16\pm0.01}$) is relatively slow compared to models of a highly magnetized jet. This indicates that Poynting flux conversion through the differential collimation of poloidal magnetic fields may be less efficient than expected. However, we find a non-negligible dispersion in the observed speeds for a given jet distance, making it difficult to describe the jet velocity field with a single power-law acceleration function. We discuss the possibility that the jet emission consists of multiple streamlines following different acceleration profiles, resulting in jet velocity stratification.

A new type of chromospheric jets in a polar limb of a coronal hole is discovered in the Ca II filtergram of the Solar Optical Telescope on board the \textit{Hinode}. We identify 30 jets in the Ca II movie of duration of 53 min. The average speed at their maximum heights is found to be 132$\pm$44 km s$^{-1}$ ranging from 57 km s$^{-1}$ to 264 km s$^{-1}$ along the propagation direction. The average lifetime is 20$\pm$6 ranging from 11 seconds to 36 seconds. The speed and lifetime of the jets are located at end-tails of those parameters determined for type II spicules, hence implying a new type of jets. To confirm whether these jets are different from conventional spicules, we construct a time-height image averaged over horizontal region of 1$\arcsec$, and calculate lagged cross-correlations of intensity profiles at each height with the intensity at 2 Mm. From this, we obtain a cross-correlation map as a function of lag and height. We find that the correlation curve as a function of lag time is well fitted into three different Gaussian functions whose standard deviations of the lag time are 193 seconds, 42 seconds, and 17 seconds. The corresponding propagation speeds are calculated to be 9 km s$^{-1}$, 67 km s$^{-1}$, and 121 km s$^{-1}$, respectively. The kinematic properties of the former two components seem to correspond to the 3 minutes oscillations and type II spicules, while the latter component to the jets addressed in this study.

We present the jet kinematics of the flat spectrum radio quasar (FSRQ) 4C +21.35 using time-resolved KaVA very long baseline interferometry array radio maps obtained from September 2014 to July 2016. During two out of three observing campaigns, observations were performed bi-weekly at 22 and 43 GHz quasi-simultaneously. At 22 GHz, we identified three jet components near the core with apparent speeds up to (14.4+/-2.1)c. The timing of the ejection of a new component detected in 2016 is consistent with a gamma-ray flare in November 2014. At 43 GHz, we found four inner jet (<3 mas) components with speeds from (3.5+/-1.4)c to (6.8+/-1.5)c. Jet component speeds tend to be higher with increasing distances from the core. We compared our data with archival Very Long Baseline Array (VLBA) data from the Boston University (BU) 43 GHz and the Monitoring Of Jets in Active galactic nuclei with VLBA Experiments (MOJAVE) 15.4 GHz monitoring programs. Whereas MOJAVE data and our data are in good agreement, jet speeds obtained from the BU Program data in the same time period are about twice as high as the ones we obtain from the KaVA data. The discrepancy at 43 GHz indicates that radio arrays with different angular resolution identify and trace different jet features even when the data are obtained at the same frequency and at the same time. The flux densities of jet components decay exponentially, in agreement with a synchrotron cooling time scale of about 1 year. Using known electron Lorentz factor values (about 9,000), we estimate the magnetic field strength to be around 1-3 micro-Tesla. When adopting a jet viewing angle of 5 degrees, the intrinsic jet speed is of order 0.99c.

The KVN(Korean VLBI Network)-style simultaneous multi-frequency receiving mode is demonstrated to be promising for mm-VLBI observations. Recently, other Very long baseline interferometry (VLBI) facilities all over the globe start to implement compatible optics systems. Simultaneous dual/multi-frequency VLBI observations at mm wavelengths with international baselines are thus possible. In this paper, we present the results from the first successful simultaneous 22/43 GHz dual-frequency observation with KaVA(KVN and VERA array), including images and astrometric results. Our analysis shows that the newly implemented simultaneous receiving system has brought a significant extension of the coherence time of the 43 GHz visibility phases along the international baselines. The astrometric results obtained with KaVA are consistent with those obtained with the independent analysis of the KVN data. Our results thus confirm the good performance of the simultaneous receiving systems for the non-KVN stations. Future simultaneous observations with more global stations bring even higher sensitivity and micro-arcsecond level astrometric measurements of the targets.

The Galactic Center supermassive black hole Sagittarius A* (Sgr A*) is one of the most promising targets to study the dynamics of black hole accretion and outflow via direct imaging with very long baseline interferometry (VLBI). At 3.5 mm (86 GHz), the emission from Sgr A* is resolvable with the Global Millimeter VLBI Array (GMVA). We present the first observations of Sgr A* with the phased Atacama Large Millimeter/submillimeter Array (ALMA) joining the GMVA. Our observations achieve an angular resolution of ~87{\mu}as, improving upon previous experiments by a factor of two. We reconstruct a first image of the unscattered source structure of Sgr A* at 3.5 mm, mitigating effects of interstellar scattering. The unscattered source has a major axis size of 120 $\pm$ 34{\mu}as (12 $\pm$ 3.4 Schwarzschild radii), and a symmetrical morphology (axial ratio of 1.2$^{+0.3}_{-0.2}$), which is further supported by closure phases consistent with zero within 3{\sigma}. We show that multiple disk-dominated models of Sgr A* match our observational constraints, while the two jet-dominated models considered are constrained to small viewing angles. Our long-baseline detections to ALMA also provide new constraints on the scattering of Sgr A*, and we show that refractive scattering effects are likely to be weak for images of Sgr A* at 1.3 mm with the Event Horizon Telescope. Our results provide the most stringent constraints to date for the intrinsic morphology and refractive scattering of Sgr A*, demonstrating the exceptional contribution of ALMA to millimeter VLBI.

The Alfv\'{e}n speed and plasma beta in photospheric bright points observed by the Broadband Filter Imager (BFI) of the Solar Optical Telescope (SOT) on board the \textit{Hinode} satellite, are estimated seismologically. The diagnostics is based on the theory of slow magnetoacoustic waves in a non-isothermally stratified photosphere with a uniform vertical magnetic field. We identify and track bright points in a G-band movie by using the 3D region growing method, and align them with blue continuum images to derive their brightness temperatures. From the Fourier power spectra of 118 continuum light curves made in the bright points, we find that light curves of 91 bright points have oscillations with properties which are significantly different from oscillation in quiet regions, with the periods ranging 2.2--16.2~min. We find that the model gives a moderate value of the plasma beta when $\gamma$ lies at around 5/3. The calculated Alfv\'{e}n speed is 9.68$\pm$2.02~km~s$^{-1}$, ranging in 6.3--17.4~km~s$^{-1}$. The plasma beta is estimated to be of 0.93$\pm$0.36, ranging in 0.2--1.9.

Measurement of the solar wind speed near the Sun is important for understanding the acceleration mechanism of the solar wind. In this study, we determine 2D solar wind speeds from 6 to 26 solar radii by applying Fourier motion filters to \textit{SOHO}/LASCO C3 movies observed from 1999 to 2010. Our method successfully reproduces the original flow speeds in the artificially generated data as well as streamer blobs. We measure 2D solar wind speeds from 1-day to 1-year timescales and their variation in solar cycle 24. We find that the solar wind speeds at timescales longer than a month in the solar maximum period are relatively uniform in the azimuthal direction, while they are clearly bimodal in the minimum period, as expected from the \textit{Ulysses} observations and IPS reconstruction. The bimodal structure appears at around 2006, becomes most distinctive in 2009, and abruptly disappears in 2010. The radial evolution of the solar wind speeds resembles the Parker's solar wind solution.

We present the results of comparative study of amplitude calibrations for East-Asia VLBI Network (EAVN) at 22 and 43 GHz using two different methods of an "a-priori" and a "template spectrum", particularly on lower declination sources. Using observational data sets of early EAVN observations, we investigated the elevation-dependence of the gain values at seven stations of the KaVA (KVN and VERA Array) and three additional telescopes in Japan (Takahagi 32m, Yamaguchi 32m and Nobeyama 45m). By comparing the independently obtained gain values based on these two methods, we found that the gain values from each method were consistent within 10% at elevations higher than 10 degree. We also found that the total flux densities of two images produced from the different amplitude calibrations were in agreement within 10% at both 22 and 43 GHz. By using the template spectrum method, furthermore, the additional radio telescopes can participate in the KaVA (i.e. EAVN) so that it can give a notable sensitivity increase. Therefore, our results will constrain the detailed conditions to reliably measure the VLBI amplitude using EAVN and give a potential to extend possible telescopes comprising EAVN.

We report the initial results of our high-cadence monitoring program on the radio jet in the active galaxy M87, obtained by the KVN and VERA Array (KaVA) at 22 GHz. This is a pilot study that preceded a larger KaVA-M87 monitoring program, which is currently ongoing. The pilot monitoring was mostly performed every two to three weeks from December 2013 to June 2014, at a recording rate of 1 Gbps, obtaining the data for a total of 10 epochs. We successfully obtained a sequence of good quality radio maps that revealed the rich structure of this jet from <~1 mas to 20 mas, corresponding to physical scales (projected) of ~0.1-2 pc (or ~140-2800 Schwarzschild radii). We detected superluminal motions at these scales, together with a trend of gradual acceleration. The first evidence for such fast motions and acceleration near the jet base were obtained from recent VLBA studies at 43 GHz, and the fact that very similar kinematics are seen at a different frequency and time with a different instrument suggests these properties are fundamental characteristics of this jet. This pilot program demonstrates that KaVA is a powerful VLBI array for studying the detailed structural evolution of the M87 jet and also other relativistic jets.

We investigate the scalar and tensor spectral indices of the quadratic inflation model in Eddington-inspired Born-Infeld (EiBI) gravity. We find that the EiBI corrections to the spectral indices are of second and first order in the slow-roll approximation for the scalar and tensor perturbations respectively. This is very promising since the quadratic inflation model in general relativity provides a very nice fit for the spectral indices. Together with the suppression of the tensor-to-scalar ratio EiBI inflation agrees well with the observational data.

We investigate the scalar and the tensor perturbations of the $\varphi^2$ inflation model in the strong-gravity limit of Eddington-inspired Born-Infeld (EiBI) theory. In order to consider the strong EiBI-gravity effect, we take the value of $\kappa$ large, where $\kappa$ is the EiBI theory parameter. The energy density of the Universe at the early stage is very high, and the Universe is in a strong-gravity regime. Therefore, the perturbation feature is not altered from what was investigated earlier. At the attractor inflationary stage, however, the feature is changed in the strong EiBI-gravity limit. The correction to the scalar perturbation in this limit comes mainly via the background matter field, while that to the tensor perturbation comes directly from the gravity ($\kappa$) effect. The change in the value of the scalar spectrum is little compared with that in the weak EiBI-gravity limit, or in GR. The form of the tensor spectrum is the same with that in the weak limit, but the value of the spectrum can be suppressed down to zero in the strong limit. Therefore, the resulting tensor-to-scalar ratio can also be suppressed in the same way, which makes $\varphi^2$ model in EiBI theory viable.

We investigate the tensor perturbation in the inflation model driven by a massive-scalar field in Eddington-inspired Born-Infeld gravity. For short wave-length modes, the perturbation feature is very similar to that of the usual chaotic inflation. For long wave-length modes, the perturbation exhibits a peculiar rise in the power spectrum which may leave a signature in the cosmic microwave background radiation.

We investigate a nonsingular initial state of the Universe which leads to inflation naturally. The model is described by a scalar field with a quadratic potential in Eddington-inspired Born-Infeld gravity. The curvature of this initial state is given by the mass scale of the scalar field which is much smaller than the Planck scale. Therefore, in this model, quantum gravity is not necessary in understanding this pre-inflationary stage, no matter how large the energy density becomes. The initial state in this model evolves eventually to a long inflationary period which is similar to the usual chaotic inflation.

We investigate the tensor and the scalar perturbations in the symmetric bouncing universe driven by one ordinary field and its Lee-Wick partner field which is a ghost. We obtain the even- and the odd-mode functions of the tensor perturbation in the matter-dominated regime. The tensor perturbation grows in time during the contracting phase of the Universe, and decays during the expanding phase. The power spectrum for the tensor perturbation is evaluated and the spectral index is given by $n_{\rm T} =6$. We add the analysis on the scalar perturbation by inspecting the even- and the odd-mode functions in the matter-dominated regime, which was studied numerically in our previous work. We conclude that the comoving curvature by the scalar perturbation is constant in the super-horizon scale and starts to decay in the far sub-horizon scale while the Universe expands.

We investigate the scalar perturbation in the Lee-Wick bouncing universe driven by an ordinary scalar field plus a ghost field. We consider only a symmetric evolution of the universe and the scalar fields about the bouncing point. The gauge invariant Sasaki-Mukhanov variable is numerically solved in the spatially flat gauge. We find a new form of the initial perturbation growing during the contracting phase. After the bouncing, this growing mode stabilizes to a constant mode which is responsible for the late-time power spectrum.

We explore a possibility that the daily sea-level pressure (SLP) over South Korea responds to the high-speed solar wind event. This is of interest in two aspects: First, if there is a statistical association this can be another piece of evidence showing that various meteorological observables indeed respond to variations in the interplanetary environment. Second, this can be a very crucial observational constraint since most models proposed so far are expected to preferentially work in higher latitude regions than the low latitude region studied here. We have examined daily solar wind speed ${\rm V}$, daily SLP difference ${\rm \Delta SLP}$, and daily ${\rm \log(BV^{2})}$ using the superposed epoch analysis in which the key date is set such that the daily solar wind speed exceeds 800 ${\rm kms^{-1}}$. We find that the daily ${\rm \Delta SLP}$ averaged out of 12 events reaches its peak at day +1 and gradually decreases back to its normal level. The amount of positive deviation of ${\rm \Delta SLP}$ is +2.5 hPa. The duration of deviation is a few days. We also find that ${\rm \Delta SLP}$ is well correlated with both the speed of solar wind and ${\rm \log(BV^{2})}$. The obtained linear correlation coefficients and chance probabilities with one-day lag for two cases are $r \simeq 0.81$ with $P> 99.9%$, and $r \simeq 0.84$ with $P> 99.9%$, respectively. We conclude by briefly discussing future direction to pursue.

We investigate the dependence of the amount of the observed galactic cosmic ray (GCR) influx on the solar North-South asymmetry using the neutron count rates obtained from four stations and sunspot data in archives spanning six solar cycles from 1953 to 2008. We find that the observed GCR influxes at Moscow, Kiel, Climax and Huancayo stations are more suppressed when the solar activity in the southern hemisphere is dominant compared with when the solar activity in the northern hemisphere is dominant. Its reduction rates at four stations are all larger than those of the suppression due to other factors including the solar polarity effect on the GCR influx. We perform the student's t-test to see how significant these suppressions are. It is found that suppressions due to the solar North-South asymmetry as well as the solar polarity are significant and yet the suppressions associated with the former are larger and more significant.

We investigate the BPS configuration of the multi D-vortices produced from the D2${\bar {\rm D}}$2 system. Based on the DBI-type action with a Gaussian-type runaway potential for a complex tachyon field, the BPS limit is achieved when the tachyon profile is thin. The solution states randomly-distributed $n$ static D-vortices with zero interaction. With the obtained BPS configuration, we derive the relativistic Lagrangian which describes the dynamics of free massive D-vortices. We also discuss the 90${}^{\circ}$ and 180${}^{\circ}$ scattering of two identical D-vortices, and present its implications on the reconnection in the dynamics of cosmic superstrings.

Moduli space dynamics of multi-D-vortices from D2${\bar {\rm D}}$ (equivalently, parallel straight D-strings from D3${\bar {\rm D}}$3) is systematically studied. For the BPS D-vortices, we show through exact calculations that the classical motion of randomly-distributed $n$ D-vortices is governed by a relativistic Lagrangian of free massive point-particles. When the head-on collision of two identical BPS D-vortices of zero radius is considered, it predicts either 90${}^{\circ}$ scattering or 0${}^{\circ}$ scattering equivalent to 180${}^{\circ}$ scattering. Since the former leads to a reconnection of two identical D-strings and the latter does to a case of their passing through each other, two possibilities are consistent with the prediction of string theory. It is also shown that the force between two non-BPS vortices is repulsive. Although the obtained moduli space dynamics of multi-BPS-D-vortices is exact in classical regime, the quantum effect of an F-string pair production should be included in determining the probabilities of the reconnection and the passing through for fast-moving cosmic superstrings.

We investigate the decay of an inhomogeneous D1-brane wrapped on a $S^1$ with an electric field. The model that we consider consists of an array of tachyon kink and anti-kink with a constant electric flux. Beginning with an initially static configuration, we numerically evolve the tachyon field with some perturbations under a fixed boundary condition at diametrically opposite points on the circle $S^1$. When the electric flux is smaller than the critical value, the tachyon kink becomes unstable; the tachyon field rolls down the potential, and the lower dimensional D0- and $\bar {\rm D}0$-brane become thin, which resembles the caustic formation known for this type of the system in the literature. For the supercritical values of the electric flux, the tachyon kink remains stable.

We analyze homogeneous anisotropic cosmology driven by the dilaton and the self-interacting ``massive'' antisymmetric tensor field which are indispensable bosonic degrees with the graviton in the NS-NS sector of string theories with D-branes. We found the attractor solutions for this system, which show the overall features of general solutions, and confirmed it through numerical analysis. The dilaton possesses the potential due to the presence of the D-brane and the curvature of extra dimensions. In the presence of the non-vanishing antisymmetric tensor field, the homogeneous universe expands anisotropically while the D-brane term dominates. The isotropy is recovered as the dilaton rolls down and the curvature term dominates. With the stabilizing potential for the dilaton, the isotropy can also be recovered.

We discuss the cosmological implications of an antisymmetric tensor field when it experiences a kind of ``Higgs mechanism,'' including nonlocal interaction terms, on D-branes. Even when a huge magnetic-condensation-inducing anisotropy is assumed in the early universe, the expansion of the universe leads to an isotropic $B$-matter-dominated universe.

We study five dimensional brane physics induced by an O(2) texture formed in one extra dimension. The model contains two 3-branes of nonzero tension, and the extra dimension is compact. The symmetry-breaking scale of the texture controls the particle hierarchy between the two branes. The TeV-scale particles are confined to the negative-tension brane where the observer sees gravity as essentially four dimensional. The effect of massive Kaluza-Klein gravitons is suppressed.

We investigate O(4) textures in a background with a positive cosmological constant. We find static solutions which co-move with the expanding background. There exists a solution in which the scalar field is regular at the horizon. This solution has a noninteger winding number smaller than one. There also exist solutions in which scalar-field derivatives are singular at the horizon. Such solutions can complete one winding within the horizon. If the winding number is larger than some critical value, static solutions including the regular one are unstable under perturbations.

A charged superconductiong cosmic string produces an extremely large electric field in its vicinity. This leads to vacuum instability and to the formation of a charged vacuum condensate which screens the electric charge of the string. We analyze the structure of this condensate using the Thomas-Fermi method.

Topological defects can arise in symmetry breaking models where the scalar field potential $V(\phi)$ has no minima and is a monotonically decreasing function of $|\phi|$. The properties of such vacuumless defects are quite different from those of the ``usual'' strings and monopoles. In some models such defects can serve as seeds for structure formation, or produce an appreciable density of mini-black holes.

It has been recently shown that topological defects can arise in symmetry breaking models where the scalar field potential $V(\phi)$ has no minima and is a monotonically decreasing function of $|\phi|$. Here we study the gravitational fields produced by such vacuumless defects in the cases of both global and gauge symmetry breaking. We find that a global monopole has a strongly repulsive gravitational field, and its spacetime has an event horizon similar to that in de Sitter space. A gauge monopole spacetime is essentially that of a magnetically charged black hole. The gravitational field of a global string is repulsive and that of a gauge string is attractive at small distances and repulsive at large distances. Both gauge and global string spacetimes have singularities at a finite distance from the string core.

Inflation of cosmic gauge and global strings is investigated by numerically solving the combined Einstein and field equations. Above some critical symmetry-breaking scales, the strings undergo inflation along the radial direction as well as the axial direction at the core. The nonsingular nature of the spacetimes around supercritical gauge and global strings is discussed and contrasted to the singular static solutions that have been discussed in the literature.

A thin wall approximation is exploited to describe a global monopole coupled to gravity. The core is modelled by de Sitter space; its boundary by a thin wall with a constant energy density; its exterior by the asymptotic Schwarzschild solution with negative gravitational mass $M$ and solid angle deficit, $\Delta\Omega/4\pi = 8\pi G\eta^2$, where $\eta$ is the symmetry breaking scale. The deficit angle equals $4\pi$ when $\eta=1/\sqrt{8\pi G} \equiv M_p$. We find that: (1) if $\eta <M_p$, there exists a unique globally static non-singular solution with a well defined mass, $M_0<0$. $M_0$ provides a lower bound on $M$. If $M_0<M<0$, the solution oscillates. There are no inflating solutions in this symmetry breaking regime. (2) if $\eta \ge M_p$, non-singular solutions with an inflating core and an asymptotically cosmological exterior will exist for all $M<0$. (3) if $\eta$ is not too large, there exists a finite range of values of $M$ where a non-inflating monopole will also exist. These solutions appear to be metastable towards inflation. If $M$ is positive all solutions are singular. We provide a detailed description of the configuration space of the model for each point in the space of parameters, $(\eta, M)$ and trace the wall trajectories on both the interior and the exterior spacetimes. Our results support the proposal that topological defects can undergo inflation.

The evolution of a global monopole with an inflating core is investigated. An analytic expression for the exterior metric at large distances from the core is obtained. The overall spacetime structure is studied numerically, both in vacuum and in a radiation background.