The intergalactic medium (IGM) at $z\sim$ 5 to 6 is largely ionized, and yet the main source for the IGM ionization in the early universe is uncertain. Of the possible contributors are faint quasars with $-26 \lesssim M_{\rm 1450} \lesssim -23$, but their number density is poorly constrained at $z\sim5$. In this paper, we present our survey of faint quasars at $z\sim5$ in the European Large-Area {\it ISO} Survey-North 1 (ELAIS-N1) field over a survey area of 6.51 deg$^2$ and examine if such quasars can be the dominant source of the IGM ionization. We use the deep optical/near-infrared data of the ELAIS-N1 field as well as the additional medium-band observations to find $z \sim 5$ quasars through a two-step approach using the broadband color selection, and SED fitting with the medium-band information included. Adopting Bayesian information criterion, we identify ten promising quasar candidates. Spectra of three of the candidates are obtained, confirming all of them to be quasars at $z\sim5$ and supporting the reliability of the quasar selection. Using the promising candidates, we derive the $z\sim5$ quasar luminosity function at $-26 \lesssim M_{\rm 1450} \lesssim -23$. The number density of faint $z\sim5$ quasars in the ELAIS-N1 field is consistent with several previous results that quasars are not the main contributors to the IGM-ionizing photons at $z\sim5$

Active galactic nuclei (AGNs) show a correlation between the size of the broad line region (BLR) and the monochromatic continuum luminosity at 5100 \AA, allowing black hole mass estimation based on single-epoch spectra. However, the validity of the correlation is yet to be clearly tested for high-luminosity AGNs. We present the first reverberation-mapping results of the Seoul National University AGN monitoring program (SAMP), which is designed to focus on luminous AGNs for probing the high end of the size-luminosity relation. We report time lag measurements of two AGNs, namely, 2MASS J10261389+5237510 and SDSS J161911.24+501109.2, using the light curves obtained over a $\sim$1000 day period with an average cadence of $\sim$10 and $\sim$20 days, respectively for photometry and spectroscopy monitoring. Based on a cross-correlation analysis and H$\beta$ line width measurements, we determine the H$\beta$ lag as $41.8^{+4.9}_{-6.0}$ and $52.6^{+17.6}_{-14.7}$ days in the observed-frame, and black hole mass as $3.65^{+0.49}_{-0.57} \times 10^7 M_{\odot}$ and $23.02^{+7.81}_{-6.56} \times 10^7 M_{\odot}$, respectively for 2MASS J1026 and SDSS J1619.

While the reverberation mapping technique is the best available method for measuring black hole mass in active galactic nuclei (AGN) beyond the local volume, this method has been mainly applied to relatively low-to-moderate luminosity AGNs at low redshift. We present the strategy of the Seoul National University AGN Monitoring Project, which aims at measuring the time delay of the \Hb\ line emission with respect to AGN continuum, using a sample of relatively high luminosity AGNs out to z $\sim$0.5. We present simulated cross correlation results based on a number of mock light curves, in order to optimally determine monitoring duration and cadence. We describe our campaign strategy based on the simulation results and the availability of observing facilities. We present the sample selection, and the properties of the selected 100 AGNs, including the optical luminosity, expected time lag, black hole mass, and Eddington ratio.

Intensive Monitoring Survey of Nearby Galaxies (IMSNG) is a high cadence observation program monitoring nearby galaxies with high probabilities of hosting supernovae (SNe). IMSNG aims to constrain the SN explosion mechanism by inferring sizes of SN progenitor systems through the detection of the shock-heated emission that lasts less than a few days after the SN explosion. To catch the signal, IMSNG utilizes a network of 0.5-m to 1-m class telescopes around the world and monitors the images of 60 nearby galaxies at distances D < 50 Mpc to a cadence as short as a few hours. The target galaxies are bright in near-ultraviolet (NUV) with M_NUV < -18.4 AB mag and have high probabilities of hosting SNe (0.06 SN/yr per galaxy). With this strategy, we expect to detect the early light curves of 3.4 SNe per year to a depth of R ~ 19.5 mag, enabling us to detect the shock-heated emission from a progenitor star with a radius as small as 0.1 R_sun. The accumulated data will be also useful for studying faint features around the target galaxies and other science projects. So far, 18 SNe have occurred in our target fields (16 in IMSNG galaxies) over 5 years, confirming our SN rate estimate of 0.06 SN/yr per galaxy.

The faint quasars with $M_{1450}>-24$ mag are known to hold the key to the determination of the ultraviolet emissivity for the cosmic re-ionization. But only a few have been identified so far because of the limitations on the survey data. Here, we present the first results of the $z\sim5$ faint quasar survey with the Infrared Medium-deep Survey (IMS), which covers $\sim100$ deg$^{2}$ areas in $J$-band to the depths of $J_{\rm AB}\sim23$ mag. To improve selection methods, the medium-band follow-up imaging has been carried out using the SED camera for QUasars in Early uNiverse (SQUEAN) on the Otto Struve 2.1 m Telescope. The optical spectra of the candidates were obtained with 8-m class telescopes. We newly discovered 10 quasars with $-25<M_{1450}<-23$ at $z\sim5$, among which three have been missed in a previous survey using the same optical data over the same area, implying the necessity for improvements in high redshift faint quasars selection. We derived photometric redshifts from the medium-band data, and find that they have high accuracies of $\langle|\Delta z|/(1+z)\rangle=0.016$. The medium-band-based approach allows us to rule out many of the interlopers that contaminate $\gtrsim20~\%$ of the broad-band-selected quasar candidates. These results suggest that the medium-band-based approach is a powerful way to identify $z\sim5$ quasars and measure their redshifts at high accuracy (1-2 %). It is also a cost-effective way to understand the contribution of quasars to the cosmic re-ionization history.

To date, most of the luminous quasars known at $z\sim6$ have been found to be in maximal accretion with the Eddington ratios, $\lambda_{\rm{Edd}}\sim1$, suggesting enhanced nuclear activities in the early universe. However, this may not be the whole picture of supermassive black hole (SMBH) growth since previous studies have not reached on faint quasars that are more likely to harbor SMBHs with low $\lambda_{\rm{Edd}}$. To gain a better understanding on the accretion activities in quasars in the early universe, we obtained a deep near-infrared (NIR) spectrum of a quasar, IMS J220417.92+011144.8 (hereafter IMS J2204+0112), one of the faintest quasars that have been identified at $z\sim6$. From the redshifted C IV $\lambda 1549$ emission line in the NIR spectrum, we find that IMS J2204+0112 harbors a SMBH with about a billion solar mass and $\lambda_{\rm{Edd}} \sim 0.1$, but with a large uncertainty in both quantities (0.41 dex). IMS J2204+0112 has one of the lowest Eddington ratios among quasars at $z\sim6$, but a common value among quasars at $z\sim2$. Its low $\lambda_{\rm{Edd}}$ can be explained with two scenarios; the SMBH growth from a stellar mass black hole through short-duration super-Eddington accretion events or from a massive black hole seed ($\sim10^{5}\,M_{\odot}$) with Eddington-limited accretion. NIR spectra of more faint quasars are needed to better understand the accretion activities of SMBHs at $z \sim 6$.

We present the variability and time lag measurements of PG 0934+013 based on a photometric and spectroscopic monitoring campaign over a two year period. We obtained 46 epochs of data from the spectroscopic campaign, which was carried out using the Southern African Large Telescope with $\sim$1 week cadence over two sets of 4 month-long observing period, while we obtained 80 epochs of \textit{B}-band imaging data using a few 1-m class telescopes. Due to the seven month gap between the two observing periods, we separately measured the time lags of broad emission lines including H$\beta$, by comparing the emission line light curve with the \textit{B}-band continuum light curve using the cross-correlation function techniques. We determined the H$\beta$ lag, $\tau_{\rm cent} = 8.46^{+2.08}_{-2.14}$ days in the observed-frame based on Year 2 data, while the time lag from Year 1 data was not reliably determined. Using the rms spectrum of Year 2 data, we measured the \Hb\ line dispersion \sigmaline = 668 $\pm$ 44 \kms\ after correcting for the spectral resolution. Adopting a virial factor f = 4.47 from Woo et al. 2015, we determined the black hole mass M$_{BH}$ = $3.13 ^{+0.91} _{-0.93} \times 10^{6}$ \msun based on the \Hb\ time lag and velocity.

We present our first results of the survey for high redshift quasars at $5 \lesssim {\rm z} \lesssim 5.7$. The search for quasars in this redshift range has been known to be challenging due to limitations of filter sets used in previous studies. We conducted a quasar survey for two specific redshift ranges, 4.60 $\leq$ z $\leq$ 5.40 and 5.50 $\leq$ z $\leq$ 6.05, using multi-wavelength data that include observations using custom-designed filters, $is$ and $iz$. Using these filters and a new selection technique, we were able to reduce the fraction of interlopers. Through optical spectroscopy, we confirmed seven quasars at 4.7 $\leq$ z $\leq$ 5.4 with $-27.4 < M_{1450} < -26.4$ which were discovered independently by another group recently. We estimated black hole masses and Eddington ratios of four of these quasars from optical and near-infrared spectra, and found that these quasars are undergoing nearly Eddington-limited accretion which is consistent with the rapid growth of supermassive black holes in luminous quasars at z $\sim$ 5.

The types of instability in the interacting binary stars are reviewed. The project "Inter-Longitude Astronomy" is a series of smaller projects on concrete stars or groups of stars. It has no special funds, and is supported from resources and grants of participating organizations, when informal working groups are created. Totally we studied 1900+ variable stars of different types. The characteristic timescale is from seconds to decades and (extrapolating) even more. The monitoring of the first star of our sample AM Her was initiated by Prof. V.P. Tsesevich (1907-1983). Since more than 358 ADS papers were published. Some highlights of our photometric and photo-polarimetric monitoring and mathematical modelling of interacting binary stars of different types are presented: classical, asynchronous, intermediate polars and magnetic dwarf novae (DO Dra) with 25 timescales corresponding to different physical mechanisms and their combinations (part "Polar"); negative and positive superhumpers in nova-like and many dwarf novae stars ("Superhumper"); eclipsing "non-magnetic" cataclysmic variables; symbiotic systems ("Symbiosis"); super-soft sources (SSS, QR And); spotted (and not spotted) eclipsing variables with (and without) evidence for a current mass transfer ("Eclipser") with a special emphasis on systems with a direct impact of the stream into the gainer star's atmosphere, or V361 Lyr-type stars. Other parts of the ILA project are "Stellar Bell" (interesting pulsating variables of different types and periods - M, SR, RV Tau, RR Lyr, Delta Sct) and "Novice"(="New Variable") discoveries and classification with a subsequent monitoring for searching and studying possible multiple components of variability. Special mathematical methods have been developed to create a set of complementary software for statistically optimal modelling of variable stars of different types.

A spectroscopic and photometric analysis has been carried out for thirty-two candidate SX_Phe variable blue straggler stars in the Kepler-field (Balona & Nemec 2012). Radial velocities (RVs), space motions (U,V,W), projected rotation velocities (v sin i), spectral types, and atmospheric characteristics (T_eff, log g, [Fe/H], xi_t, zeta_RT, etc.) are presented for 30 of the 32 stars. Although several stars are metal-weak with extreme halo orbits, the mean [Fe/H] of the sample is near solar, thus the stars are more metal-rich than expected for a typical sample of Pop.II stars, and more like halo metal-rich A-type stars (Perry 1969). Two thirds of the stars are fast rotators with v sin i > 50 km/s, including four stars with v sin i > 200 km/s. Three of the stars have (negative) RVs > 250 km/s, five have retrograde space motions, and 21 have total speeds (relative to the LSR) > 400 km/s. All but one of the 30 stars have positions in a Toomre diagram consistent with the kinematics of bona fide halo stars (the exception being a thick-disk star). Observed Roemer time delays, pulsation frequency modulations and light curves suggest that at least one third of the stars are in binary (or triple) systems with orbital periods ranging from 2.3 days to more than four years.

We describe the characteristics and performance of a camera system, Spectral energy distribution Camera for Quasars in Early Universe (SQUEAN). It was developed to measure SEDs of high redshift quasar candidates (z $\gtrsim$ 5) and other targets, e.g., young stellar objects, supernovae, and gamma-ray bursts, and to trace the time variability of SEDs of objects such as active galactic nuclei (AGNs). SQUEAN consists of an on-axis focal plane camera module, an auto-guiding system, and mechanical supporting structures. The science camera module is composed of a focal reducer, a customizable filter wheel, and a CCD camera on the focal plane. The filter wheel uses filter cartridges that can house filters with different shapes and sizes, enabling the filter wheel to hold twenty filters of 50 mm $\times$ 50 mm size, ten filters of 86 mm $\times$ 86 mm size, or many other combinations. The initial filter mask was applied to calibrate the filter wheel with high accuracy and we verified that the filter position is repeatable at much less than one pixel accuracy. We installed and tested 50 nm medium bandwidth filters of 600 -- 1,050 nm and other filters at the commissioning observation in 2015 February. We found that SQUEAN can reach limiting magnitudes of 23.3 - 25.3 AB mag at 5-$\sigma$ in a 1-hour total integration time. - 25.3 AB mag at 5-$\sigma$ in a 1-hour total integration time.

The $\Lambda$CDM cosmological model successfully reproduces many aspects of the galaxy and structure formation of the universe. However, the growth of large-scale structures (LSSs) in the early universe is not well tested yet with observational data. Here, we have utilized wide and deep optical--near-infrared data in order to search for distant galaxy clusters and superclusters ($0.8<z<1.2$). From the spectroscopic observation with the Inamori Magellan Areal Camera and Spectrograph (IMACS) on the Magellan telescope, three massive clusters at $z\sim$0.91 are confirmed in the SSA22 field. Interestingly, all of them have similar redshifts within $\Delta z\sim$0.01 with velocity dispersions ranging from 470 to 1300 km s$^{-1}$. Moreover, as the maximum separation is $\sim$15 Mpc, they compose a supercluster at $z\sim$0.91, meaning that this is one of the most massive superclusters at this redshift to date. The galaxy density map implies that the confirmed clusters are embedded in a larger structure stretching over $\sim$100 Mpc. $\Lambda$CDM models predict about one supercluster like this in our surveyed volume, consistent with our finding so far. However, there are more supercluster candidates in this field, suggesting that additional studies are required to determine if the $\Lambda$CDM cosmological model can successfully reproduce the LSSs at high redshift.

Multiple color selection techniques have been successful in identifying quasars from wide-field broad-band imaging survey data. Among the quasars that have been discovered so far, however, there is a redshift gap at $5 \lesssim {\rm z} \lesssim 5.7$ due to the limitations of filter sets in previous studies. In this work, we present a new selection technique of high redshift quasars using a sequence of medium-band filters: nine filters with central wavelengths from 625 to 1025 nm and bandwidths of 50 nm. Photometry with these medium-bands traces the spectral energy distribution (SED) of a source, similar to spectroscopy with resolution R $\sim$ 15. By conducting medium-band observations of high redshift quasars at 4.7 $\leq$ z $\leq$ 6.0 and brown dwarfs (the main contaminants in high redshift quasar selection) using the SED camera for QUasars in EArly uNiverse (SQUEAN) on the 2.1-m telescope at the McDonald Observatory, we show that these medium-band filters are superior to multi-color broad-band color section in separating high redshift quasars from brown dwarfs. In addition, we show that redshifts of high redshift quasars can be determined to an accuracy of $\Delta{\rm z}/(1+{\rm z}) = 0.002$ -- $0.026$. The selection technique can be extended to z $\sim$ 7, suggesting that the medium-band observation can be powerful in identifying quasars even at the re-ionization epoch.

Recent studies suggest that faint active galactic nuclei may be responsible for the reionization of the universe. Confirmation of this scenario requires spectroscopic identification of faint quasars ($M_{1450}>-24$ mag) at $z \gtrsim6$, but only a very small number of such quasars have been spectroscopically identified so far. Here, we report the discovery of a faint quasar IMS J220417.92+011144.8 at z~6 in a 12.5 deg$^{2}$ region of the SA22 field of the Infrared Medium-deep Survey (IMS). The spectrum of the quasar shows a sharp break at $\sim8443~\rm{\AA}$, with emission lines redshifted to $z=5.944 \pm 0.002$ and rest-frame ultraviolet continuum magnitude $M_{1450}=-23.59\pm0.10$ AB mag. The discovery of IMS J220417.92+011144.8 is consistent with the expected number of quasars at z~6 estimated from quasar luminosity functions based on previous observations of spectroscopically identified low-luminosity quasars . This suggests that the number of $M_{1450}\sim-23$ mag quasars at z~6 may not be high enough to fully account for the reionization of the universe. In addition, our study demonstrates that faint quasars in the early universe can be identified effectively with a moderately wide and deep near-infrared survey such as the IMS.

We study the host galaxy properties of the tidal disruption object, Swift J164449.3+573451 using long-term optical to near-infrared (NIR) data. First, we decompose the galaxy surface brightness distribution and analyze the morphology of the host galaxy using high resolution \emph{HST} WFC3 images. We conclude that the host galaxy is a bulge-dominant galaxy that is well described by a single S\'{e}rsic model with S\'{e}rsic index $n=3.43\pm0.05$. Adding a disk component, the bulge to total host galaxy flux ratio (B/T) is $0.83\pm0.03$, which still indicates a bulge-dominant galaxy. Second, we estimate multi-band fluxes of the host galaxy through long-term light curves. Our long-term NIR light curves reveal the pure host galaxy fluxes $\sim500$ days after the burst. We fit spectral energy distribution (SED) models to the multi-band fluxes from the optical to NIR of the host galaxy and determine its properties. The stellar mass, the star formation rate, and the age of stellar population are $\log(M_{\star}/M_{\odot}) = 9.14^{+0.13}_{-0.10}$, $0.03^{+0.28}_{-0.03}\, M_{\odot}$/yr, and $0.63^{+0.95}_{-0.43}$ Gyr. Finally, we estimate the mass of the central super massive black hole which is responsible for the tidal disruption event. The black hole mass is estimated to be $10^{6.7\pm0.4}\, M_{\odot}$ from $M_{\mathrm{BH}}$ - $M_{\star,\mathrm{bul}}$ and $M_{\mathrm{BH}}$ - $L_{\mathrm{bul}}$ relations for the $K$ band, although a smaller value of $\sim10^5\, M_{\odot}$ cannot be excluded convincingly if the host galaxy harbors a pseudobulge.

We study the dependence of angular two-point correlation functions on stellar mass ($M_{*}$) and specific star formation rate (sSFR) of $M_{*}>10^{10}M_{\odot}$ galaxies at $z\sim1$. The data from UKIDSS DXS and CFHTLS covering 8.2 deg$^{2}$ sample scales larger than 100 $h^{-1}$Mpc at $z\sim1$, allowing us to investigate the correlation between clustering, $M_{*}$, and star formation through halo modeling. Based on halo occupation distributions (HODs) of $M_{*}$ threshold samples, we derive HODs for $M_{*}$ binned galaxies, and then calculate the $M_{*}/M_{\rm halo}$ ratio. The ratio for central galaxies shows a peak at $M_{\rm halo}\sim10^{12}h^{-1}M_{\odot}$, and satellites predominantly contribute to the total stellar mass in cluster environments with $M_{*}/M_{\rm halo}$ values of 0.01--0.02. Using star-forming galaxies split by sSFR, we find that main sequence galaxies ($\rm log\,sSFR/yr^{-1}\sim-9$) are mainly central galaxies in $\sim10^{12.5} h^{-1}M_{\odot}$ haloes with the lowest clustering amplitude, while lower sSFR galaxies consist of a mixture of both central and satellite galaxies where those with the lowest $M_{*}$ are predominantly satellites influenced by their environment. Considering the lowest $M_{\rm halo}$ samples in each $M_{*}$ bin, massive central galaxies reside in more massive haloes with lower sSFRs than low mass ones, indicating star-forming central galaxies evolve from a low $M_{*}$--high sSFR to a high $M_{*}$--low sSFR regime. We also find that the most rapidly star-forming galaxies ($\rm log\,sSFR/yr^{-1}>-8.5$) are in more massive haloes than main sequence ones, possibly implying galaxy mergers in dense environments are driving the active star formation. These results support the conclusion that the majority of star-forming galaxies follow secular evolution through the sustained but decreasing formation of stars.

We present the J and H-band source catalog covering the AKARI North Ecliptic Pole field. Filling the gap between the optical data from other follow-up observations and mid-infrared (MIR) data from AKARI, our near-infrared (NIR) data provides contiguous wavelength coverage from optical to MIR. For the J and H-band imaging, we used the FLoridA Multi-object Imaging Near-ir Grism Observational Spectrometer (FLAMINGOS) on the Kitt Peak National Observatory 2.1m telescope covering a 5.1 deg2 area down to a 5 sigma depth of ~21.6 mag and ~21.3 mag (AB) for J and H-band with an astrometric accuracy of 0.14" and 0.17" for 1 sigma in R.A. and Decl. directions, respectively. We detected 208,020 sources for J-band and 203,832 sources for H-band. This NIR data is being used for studies including analysis of the physical properties of infrared sources such as stellar mass and photometric redshifts, and will be a valuable dataset for various future missions.

We carried out photometric observations for HBC 722 in SDSS r, i and z bands from 2011 April to 2013 May with a Camera for Quasars in Early uNiverse attached to the 2.1m Otto Struve telescope at McDonald Observatory. The post-outburst phenomena were classified into five phases according to not only brightness but also color variations, which might be caused by physical changes in the emitting regions of optical and near-infrared bands. A series of spectral energy distribution (SED) is presented to support color variations and track the time evolution of SED in optical/near-infrared bands after the outburst. Given two years of data, possible periodicities of r, i and z bands were checked. We found out three families of signals around ~6, ~10 and ~1 days in three bands, which is broadly consistent with Green et al. (2013). We also examined short term variability (intra-day and day scales) to search for evidences of flickering by using the micro-variability method. We found clear signs of day scale variability and weak indications of intra-day scale fluctuations, which implies that the flickering event occurs in HBC 722 after outburst.

We present a wide dataset of gamma-ray, X-ray, UVOIR, and radio observations of the Swift GRB100814A. At the end of the slow decline phase of the X-ray and optical afterglow, this burst shows a sudden and prominent rebrightening in the optical band only, followed by a fast decay in both bands. The optical rebrightening also shows chromatic evolution. Such a puzzling behaviour cannot be explained by a single component model. We discuss other possible interpretations, and we find that a model that incorporates a long-lived reverse shock and forward shock fits the temporal and spectral properties of GRB100814A the best.

We present the analysis of four first overtone RR Lyrae stars observed with the Kepler space telescope, based on data obtained over nearly 2.5yr. All four stars are found to be multiperiodic. The strongest secondary mode with frequency f_2 has an amplitude of a few mmag, 20 - 45 times lower than the main radial mode with frequency f_1. The two oscillations have a period ratio of P_2/P_1 = 0.612 - 0.632 that cannot be reproduced by any two radial modes. Thus, the secondary mode is nonradial. Modes yielding similar period ratios have also recently been discovered in other variables of the RRc and RRd types. These objects form a homogenous group and constitute a new class of multimode RR Lyrae pulsators, analogous to a similar class of multimode classical Cepheids in the Magellanic Clouds. Because a secondary mode with P_2/P_1 ~ 0.61 is found in almost every RRc and RRd star observed from space, this form of multiperiodicity must be common. In all four Kepler RRc stars studied, we find subharmonics of f_2 at ~1/2 f_2 and at ~3/2 f_2. This is a signature of period doubling of the secondary oscillation, and is the first detection of period doubling in RRc stars. The amplitudes and phases of f_2 and its subharmonics are variable on a timescale of 10 - 200d. The dominant radial mode also shows variations on the same timescale, but with much smaller amplitude. In three Kepler RRc stars we detect additional periodicities, with amplitudes below 1mmag, that must correspond to nonradial g-modes. Such modes never before have been observed in RR Lyrae variables.

Activity at the centers of galaxies, during which the central supermassive black hole is accreting material, is nowadays accepted to be rather ubiquitous and most probably a phase of every galaxy's evolution. It has been suggested that galactic mergers and interactions may be the culprits behind the triggering of nuclear activity. We use near-infrared data from the new Infrared Medium-Deep Survey (IMS) and the Deep eXtragalactic Survey (DXS) of the VIMOS-SA22 field and radio data at 1.4 GHz from the FIRST survey and a deep VLA survey to study the environments of radio-AGN over an area of ~25 sq. degrees and down to a radio flux limit of 0.1 mJy and a J-band magnitude of 23 mag AB. Radio-AGN are predominantly found in environments similar to those of control galaxies at similar redshift, J-band magnitude, and U-R rest-frame absolute color. However, a sub-population of radio-AGN is found in environments up to 100 times denser than their control sources. We thus preclude merging as the dominant triggering mechanism of radio-AGN. Through the fitting of the broadband spectral energy distribution of radio-AGN in the least and most dense environments, we find that those in the least dense environments show higher radio-loudness, higher star formation efficiencies, and higher accretion rates, typical of the so-called high-excitation radio-AGN. These differences tend to disappear at z>1. We interpret our results in terms of a different triggering mechanism for these sources that is driven by mass-loss through winds of young stars created during the observed ongoing star formation.

Context. The Vista Variables in the V\'ia L\'actea (VVV) ESO Public Survey is a variability survey of the Milky Way bulge and an adjacent section of the disk carried out from 2010 on ESO Visible and Infrared Survey Telescope for Astronomy (VISTA). VVV will eventually deliver a deep near-IR atlas with photometry and positions in five passbands (ZYJHK_S) and a catalogue of 1-10 million variable point sources - mostly unknown - which require classifications. Aims. The main goal of the VVV Templates Project, that we introduce in this work, is to develop and test the machine-learning algorithms for the automated classification of the VVV light-curves. As VVV is the first massive, multi-epoch survey of stellar variability in the near-infrared, the template light-curves that are required for training the classification algorithms are not available. In the first paper of the series we describe the construction of this comprehensive database of infrared stellar variability. Methods. First we performed a systematic search in the literature and public data archives, second, we coordinated a worldwide observational campaign, and third we exploited the VVV variability database itself on (optically) well-known stars to gather high-quality infrared light-curves of several hundreds of variable stars. Results. We have now collected a significant (and still increasing) number of infrared template light-curves. This database will be used as a training-set for the machine-learning algorithms that will automatically classify the light-curves produced by VVV. The results of such an automated classification will be covered in forthcoming papers of the series.

Although more than eight decades have passed since P. Th. Oosterhoff drew attention to differences in the properties of RR Lyrae variables in globular clusters, the origin and significance of the Oosterhoff groups remain unclear. Nonetheless, the accumulation of extensive new observations of RR Lyrae stars in globular clusters of the Milky Way and Local Group galaxies allows a fresh look at the phenomenon. Insights come not only from surveys of variables within the original Oosterhoff groups I and II but also from recent observations of the Oosterhoff-intermediate systems found especially in smaller Local Group galaxies. We will compare properties of RR Lyrae in several systems to investigate what they reveal about system-to-system differences of transition temperature between fundamental-mode and first overtone pulsators and of horizontal branch luminosity. Both transition temperature and horizontal branch luminosity have at various times been credited as playing roles in the creation of the Oosterhoff dichotomy.

Several lines of argument support the existence of a link between activity at the nuclei of galaxies, in the form of an accreting supermassive black hole, and star-formation activity in these galaxies. The exact nature of this link is still under debate. Radio jets have long been argued to be an ideal mechanism that allows AGN to interact with their host galaxy and regulate star-formation. In this context, we are using a sample of radio sources in the North Ecliptic Pole (NEP) field to study the nature of the putative link between AGN activity and star-formation. This is done by means of spectral energy distribution (SED) fitting. We use the excellent spectral coverage of the AKARI infrared space telescope together with the rich ancillary data available in the NEP to build SEDs extending from UV to far-IR wavelengths. Through SED fitting we constrain both the AGN and host galaxy components. We find a significant AGN component in our sample of relatively faint radio-sources ($<$mJy), that increases in power with increasing radio-luminosity. At the highest radio-luminosities, the presence of powerful jets dominates the radio emission of these sources. A positive correlation is found between the luminosity of the AGN component and that of star-formation in the host galaxy, independent of the radio luminosity. By contrast, for a given redshift and AGN luminosity, we find that increasing radio-luminosity leads to a decrease in the specific star-formation rate. The most radio-loud AGN are found to lie on the main sequence of star-formation for their respective redshifts. For the first time, such a two-sided feedback process is seen in the same sample. We conclude that radio jets do suppress star-formation in their host galaxies but appear not to totally quench it. Our results therefore support the maintenance nature of "radio-mode" feedback from radio-AGN jets.

We present the discovery and extensive early-time observations of the Type Ic supernova (SN) PTF12gzk. Our finely sampled light curves show a rise of 0.8mag within 2.5hr. Power-law fits [f(t)\sim(t-t_0)^n] to these data constrain the explosion date to within one day. We cannot rule out the expected quadratic fireball model, but higher values of n are possible as well for larger areas in the fit parameter space. Our bolometric light curve and a dense spectral sequence are used to estimate the physical parameters of the exploding star and of the explosion. We show that the photometric evolution of PTF12gzk is slower than that of most SNe Ic, and its high ejecta velocities (~30,000km/s four days after explosion) are closer to the observed velocities of broad-lined SNe Ic associated with gamma-ray bursts (GRBs) than to the observed velocities in normal Type Ic SNe. The high velocities are sustained through the SN early evolution, and are similar to those of GRB-SNe when the SN reach peak magnitude. By comparison with the spectroscopically similar SN 2004aw, we suggest that the observed properties of PTF12gzk indicate an initial progenitor mass of 25-35 solar mass and a large (5-10E51 erg) kinetic energy, close to the regime of GRB-SN properties. The host-galaxy characteristics are consistent with GRB-SN hosts, and not with normal SN Ic hosts as well, yet this SN does not show the broad lines over extended periods of time that are typical of broad-line Type Ic SNe.

This is the third in a series of papers studying the variable stars in old globular clusters in the Large Magellanic Cloud. The primary goal of this series is to look at how the characteristics and behavior of RR Lyrae stars in Oosterhoff-intermediate systems compare to those of their counterparts in Oosterhoff-I/II systems. In this paper we present the results of our new time-series BVI photometric study of the globular cluster Reticulum. We found a total of 32 variables stars (22 RRab, 4 RRc, and 6 RRd stars) in our field of view. We present photometric parameters and light curves for these stars. We also present physical properties, derived from Fourier analysis of light curves, for some of the RR Lyrae stars. We discuss the Oosterhoff classification of Reticulum and use our results to re-derive the distance modulus and age of the cluster.

We present ten medium-resolution, high signal-to-noise ratio near-infrared (NIR) spectra of SN 2011fe from SpeX on the NASA Infrared Telescope Facility (IRTF) and Gemini Near-Infrared Spectrograph (GNIRS) on Gemini North, obtained as part of the Carnegie Supernova Project. This data set constitutes the earliest time-series NIR spectroscopy of a Type Ia supernova (SN Ia), with the first spectrum obtained at 2.58 days past the explosion and covering -14.6 to +17.3 days relative to B-band maximum. C I {\lambda}1.0693 {\mu}m is detected in SN 2011fe with increasing strength up to maximum light. The delay in the onset of the NIR C I line demonstrates its potential to be an effective tracer of unprocessed material. For the first time in a SN Ia, the early rapid decline of the Mg II {\lambda}1.0927 {\mu}m velocity was observed, and the subsequent velocity is remarkably constant. The Mg II velocity during this constant phase locates the inner edge of carbon burning and probes the conditions under which the transition from deflagration to detonation occurs. We show that the Mg II velocity does not correlate with the optical light-curve decline rate {\Delta}m15. The prominent break at ~1.5 {\mu}m is the main source of concern for NIR k-correction calculations. We demonstrate here that the feature has a uniform time evolution among SNe Ia, with the flux ratio across the break strongly correlated with {\Delta}m15. The predictability of the strength and the onset of this feature suggests that the associated k-correction uncertainties can be minimized with improved spectral templates.

We present our successful Chandra program designed to identify, with sub-arcsecond accuracy, the X-ray afterglow of the short GRB 111117A, which was discovered by Swift and Fermi. Thanks to our rapid target of opportunity request, Chandra clearly detected the X-ray afterglow, though no optical afterglow was found in deep optical observations. The host galaxy was clearly detected in the optical and near-infrared band, with the best photometric redshift of z=1.31_{-0.23}^{+0.46} (90% confidence), making it one of the highest known short GRB redshifts. Furthermore, we see an offset of 1.0 +- 0.2 arcseconds, which corresponds to 8.4 +- 1.7 kpc, between the host and the afterglow position. We discuss the importance of using Chandra for obtaining sub-arcsecond X-ray localizations of short GRB afterglows to study GRB environments.

We present the detection of day-timescale periodic variability in the r-band lightcurve of newly outbursting FU Orionis-type object HBC 722, taken from > 42 nights of observation with the CQUEAN instrument on the McDonald Observatory 2.1m telescope. The optical/near-IR lightcurve of HBC 722 shows a complex array of periodic variability, clustering around 5.8 day (0.044 mag amplitude) and 1.28 day (0.016 mag amplitude) periods, after removal of overall baseline variation. We attribute the unusual number of comparable strength signals to a phenomenon related to the temporary increase in accretion rate associated with FUors. We consider semi-random "flickering", magnetic braking/field compression and rotational asymmetries in the disk instability region as potential sources of variability. Assuming the 5.8 day period is due to stellar rotation and the 1.28 day period is indicative of Keplerian rotation at the inner radius of the accretion disk (at 2 R(star)), we derive a B-field strength of 2.2-2.7 kG, slightly larger than typical T Tauri stars. If instead the 5.8 day signal is from a disk asymmetry, the instability region has an outer radius of 5.4 R(star), consistent with models of FUor disks. Further exploration of the time domain in this complicated source and related objects will be key to understanding accretion processes.

We present a photometric catalog of infrared (IR) sources based on the North Ecliptic PoleWide field (NEP-Wide) survey of AKARI, which is an infrared space telescope launched by Japan. The NEP-Wide survey covered 5.4 deg2 area, a nearly circular shape centered on the North Ecliptic Pole, using nine photometric filter-bands from 2 - 25 {\mu}m of the Infrared Camera (IRC). Extensive efforts were made to reduce possible false objects due to cosmic ray hits, multiplexer bleeding phenomena around bright sources, and other artifacts. The number of detected sources varied depending on the filter band: with about 109,000 sources being cataloged in the near-IR bands at 2 - 5 {\mu}m, about 20,000 sources in the shorter parts of the mid-IR bands between 7 - 11 {\mu}m, and about 16,000 sources in the longer parts of the mid-IR bands, with \sim 4,000 sources at 24 {\mu}m. The estimated 5? detection limits are approximately 21 magnitude (mag) in the 2 - 5 {\mu}m bands, 19.5 - 19 mag in the 7 - 11 {\mu}m, and 18.8 - 18.5 mag in the 15 - 24 {\mu}m bands in the AB magnitude scale. The completenesses for those bands were evaluated as a function of magnitude: the 50% completeness limits are about 19.8 mag at 3 {\mu}m, 18.6 mag at 9 {\mu}m, and 18 mag at 18 {\mu}m band, respectively. To construct a reliable source catalog, all of the detected sources were examined by matching them with those in other wavelength data, including optical and ground-based near-IR bands. The final band-merged catalog contains about 114,800 sources detected in the IRC filter bands. The properties of the sources are presented in terms of the distributions in various color-color diagrams.

New CCD photometric observations of BX Dra were obtained for 26 nights from 2009 April to 2010 June. The long-term photometric behaviors of the system are presented from detailed studies of the period and light variations, based on the historical data and our new observations. All available light curves display total eclipses at secondary minima and inverse O'Connell effects with Max I fainter than Max II, which are satisfactorily modeled by adding the slightly time-varying hot spot on the primary star. A total of 87 times of minimum light spanning over about 74 yrs, including our 22 timing measurements, were used for ephemeris computations. Detailed analysis of the O-C diagram showed that the orbital period has changed in combinations with an upward parabola and a sinusoidal variation. The continuous period increase with a rate of +5.65 \times 10^-7 d yr^-1 is consistent with that calculated from the Wilson-Devinney synthesis code. It can be interpreted as a mass transfer from the secondary to the primary star at a rate of 2.74 \times 10^-7 M\odot yr^-1, which is one of the largest rates for contact systems. The most likely explanation of the sinusoidal variation with a period of 30.2 yrs and a semi-amplitude of 0.0062 d is a light-traveltime effect due to the existence of a circumbinary object. We suggest that BX Dra is probably a triple system, consisting of a primary star with a spectral type of F0, its secondary component of spectral type F1-2, and an unseen circumbinary object with a minimum mass of M3 = 0.23 M\odot.

We present the gamma-ray, X-ray, optical and radio data for GRB100814A. At the end of the slow decline phase of the X-ray and optical afterglow, a sudden and prominent rebrightening in the optical band occurs followed by a fast decay in both bands. This optical rebrightening is accompanied by possible chromatic variations. We discuss possible interpretations, such as double component scenarios and internal dissipation mechanism, with their virtues and drawbacks. We also compare GRB100814A with other Swift bursts that show optical rebrightenings with similar properties.

We describe the overall characteristics and the performance of an optical CCD camera system, Camera for QUasars in EArly uNiverse (CQUEAN), which is being used at the 2.1 m Otto Struve Telescope of the McDonald Observatory since 2010 August. CQUEAN was developed for follow-up imaging observations of red sources such as high redshift quasar candidates (z >= 5), Gamma Ray Bursts, brown dwarfs, and young stellar objects. For efficient observations of the red objects, CQUEAN has a science camera with a deep depletion CCD chip which boasts a higher quantum efficiency at 0.7 - 1.1 um than conventional CCD chips. The camera was developed in a short time scale (~ one year), and has been working reliably. By employing an auto-guiding system and a focal reducer to enhance the field of view on the classical Cassegrain focus, we achieve a stable guiding in 20 minute exposures, an imaging quality with FWHM >= 0.6" over the whole field (4.8' * 4.8'), and a limiting magnitude of z = 23.4 AB mag at 5-sigma with one hour total integration time.

Y -band is a broad passband that is centered at ~ 1 micron. It is becoming a new, popular window for extragalactic study especially for observation of red objects thanks to recent CCD technology developments. In order to better understand the general characteristics of objects in Y -band, and to investigate the promise of Y -band observations with small telescopes, we carried out imaging observation of several extragalactic fields, brown dwarfs and high redshift quasars with Y -band filter at the Mt. Lemmon Optical Astronomy Observatory and the Maidanak observatory. From our observations, we constrain the bright end of the galaxy and the stellar number counts in Y -band. Also, we test the usefulness of high redshift quasar (z > 6) selection via i-z-Y color-color diagram, to demonstrate that the i-z-Y color-color diagram is effective for the selection of high redshift quasars even with a conventional optical CCD camera installed at a 1-m class telescope.

We present the mid-infrared (MIR) properties of galaxies within a supercluster in the North Ecliptic Pole region at z?0.087 observed with the AKARI satellite. We use data from the AKARI NEP-Wide (5.4 deg2) IR survey and the CLusters of galaxies EVoLution studies (CLEVL) mission program. We show that near-IR (3 {\mu}m)-mid- IR (11 {\mu}m) color can be used as an indicator of the specific star formation rate and the presence of intermediate age stellar populations. From the MIR observations, we find that red-sequence galaxies consist not only of passively evolving red early-type galaxies, but also of 1) "weak-SFG" (disk-dominated star-forming galaxies which have star formation rates lower by \sim 4 \times than blue-cloud galaxies), and 2) "intermediate- MXG" (bulge-dominated galaxies showing stronger MIR dust emission than normal red early-type galaxies). Those two populations can be a set of transition galaxies from blue, star-forming, late-type galaxies evolving into red, quiescent, early-type ones. We find that the weak-SFG are predominant at intermediate masses (1010M\odot < M\star < 1010.5M\odot) and are typically found in local densities similar to the outskirts of galaxy clusters. As much as 40% of the supercluster member galaxies in this mass range can be classified as weak-SFGs, but their proportion decreases to < 10% at larger masses (M\star > 1010.5 M\odot) at any galaxy density. The fraction of the intermediate-MXG among red- sequence galaxies at 1010M\odot < M\star < 1011M\odot also decreases as the density and mass increase. In particular, \sim42% of the red-sequence galaxies with early-type morphologies are classified as intermediate-MXG at intermediate densities. These results suggest that the star formation activity is strongly dependent on the stellar mass, but that the morphological transformation is mainly controlled by the environment.

This paper summarizes the main results of our recent study of the non-Blazhko RR Lyrae stars observed with the Kepler space telescope. These stars offer the opportunity for studying the stability of the pulsations of RR Lyrae stars and for providing a reference against which the Blazhko RR Lyrae stars can be compared. Of particular interest is the stability of the low-dispersion (sigma < 1mmag) light curves constructed from ~18,000 long-cadence (30-min) and (for FN Lyr and AW Dra) the ~150,000 short-cadence (1-min) photometric data points. Fourier-based [Fe/H] values and other physical characteristics are also derived. When the observed periods are compared with periods computed with the Warsaw non-linear convective pulsation code better agreement is achieved assuming pulsational L and M values rather than the (higher) evolutionary L and M values.

Nineteen of the ~40 RR Lyr stars in the Kepler field have been identified as candidate non-Blazhko (or unmodulated) stars. In this paper we present the results of Fourier decomposition of the time-series photometry of these stars acquired during the first 417 days of operation (Q0-Q5) of the Kepler telescope. Fourier parameters based on ~18400 long-cadence observations per star (and ~150000 short-cadence observations for FN Lyr and for AW Dra) are derived. None of the stars shows the recently discovered `period-doubling' effect seen in Blazhko variables; however, KIC 7021124 has been found to pulsate simultaneously in the fundamental and second overtone modes with a period ratio P2/P0 ~ 0.59305 and is similar to the double-mode star V350 Lyr. Period change rates are derived from O-C diagrams spanning, in some cases, over 100 years; these are compared with high-precision periods derived from the Kepler data alone. Extant Fourier correlations by Kovacs, Jurcsik et al. (with minor transformations from the V to the Kp passband) have been used to derive underlying physical characteristics for all the stars. This procedure seems to be validated through comparisons of the Kepler variables with galactic and LMC RR Lyr stars. The most metal-poor star in the sample is NR Lyr, with [Fe/H]=-2.3 dex; and the four most metal-rich stars have [Fe/H] ranging from -0.6 to +0.1 dex. Pulsational luminosities and masses are found to be systematically smaller than L and mass values derived from stellar evolution models, and are favoured over the evolutionary values when periods are computed with the Warsaw linear hydrodynamics code. Finally, the Fourier parameters are compared with theoretical values derived using the Warsaw non-linear convective pulsation code.

Massive black holes are believed to reside at the centres of most galaxies. They can be- come detectable by accretion of matter, either continuously from a large gas reservoir or impulsively from the tidal disruption of a passing star, and conversion of the gravitational energy of the infalling matter to light. Continuous accretion drives Active Galactic Nuclei (AGN), which are known to be variable but have never been observed to turn on or off. Tidal disruption of stars by dormant massive black holes has been inferred indirectly but the on- set of a tidal disruption event has never been observed. Here we report the first discovery of the onset of a relativistic accretion-powered jet in the new extragalactic transient, Swift J164449.3+573451. The behaviour of this new source differs from both theoretical models of tidal disruption events and observations of the jet-dominated AGN known as blazars. These differences may stem from transient effects associated with the onset of a powerful jet. Such an event in the massive black hole at the centre of our Milky Way galaxy could strongly ionize the upper atmosphere of the Earth, if beamed towards us.

Massive black holes are believed to reside at the centres of most galaxies. They can be- come detectable by accretion of matter, either continuously from a large gas reservoir or impulsively from the tidal disruption of a passing star, and conversion of the gravitational energy of the infalling matter to light. Continuous accretion drives Active Galactic Nuclei (AGN), which are known to be variable but have never been observed to turn on or off. Tidal disruption of stars by dormant massive black holes has been inferred indirectly but the on- set of a tidal disruption event has never been observed. Here we report the first discovery of the onset of a relativistic accretion-powered jet in the new extragalactic transient, Swift J164449.3+573451. The behaviour of this new source differs from both theoretical models of tidal disruption events and observations of the jet-dominated AGN known as blazars. These differences may stem from transient effects associated with the onset of a powerful jet. Such an event in the massive black hole at the centre of our Milky Way galaxy could strongly ionize the upper atmosphere of the Earth, if beamed towards us.

We report the results of our high resolution optical spectroscopic monitoring campaign ($\lambda$ = 3800 -- 8800 A, R = 30000 -- 45000) of the new FU Orionis-type object HBC 722. We observed HBC 722 with the BOES 1.8-m telescope between 2010 November 26 and 2010 December 29 and FU Orionis itself on 2011 January 26. We detect a number of previously unreported high-resolution K I and Ca II lines beyond 7500 A. We resolve the H$\alpha$ and Ca II line profiles into three velocity components, which we attribute to both disk and outflow. The increased accretion during outburst can heat the disk to produce the relatively narrow absorption feature and launch outflows appearing as high velocity blue and redshifted broad features.

We present and interpret simultaneous new photometric and spectroscopic observations of the peculiar magnetic white dwarf WD1953-011. The flux in the V-band filter and intensity of the Balmer spectral lines demonstrate variability with the rotation period of about 1.45 days. According to previous studies, this variability can be explained by the presence of a dark spot having a magnetic nature, analogous to a sunspot. Motivated by this idea, we examine possible physical relationships between the suggested dark spot and the strong-field magnetic structure (magnetic "spot", or "tube") recently identified on the surface of this star. Comparing the rotationally-modulated flux with the variable spectral observables related to the magnetic "spot" we establish their correlation, and therefore their physical relationship. Modeling the variable photometric flux assuming that it is associated with temperature variations in the stellar photosphere, we argue that the strong-field area and dark, low-temperature spot are comparable in size and located at the same latitudes, essentially overlapping each other with a possible slight longitudinal shift. In this paper we also present a new, improved value of the star's rotational period and constrain the characteristics of the thermal inhomogeneity over the degenerate's surface.

We present ground-based and HST optical observations of the optical transients (OTs) of long-duration Gamma Ray Bursts (GRBs) 060729 and 090618, both at a redshift of z = 0.54. For GRB 060729, bumps are seen in the optical light curves (LCs), and the late-time broadband spectral energy distributions (SEDs) of the OT resemble those of local type Ic supernovae (SNe). For GRB 090618, the dense sampling of our optical observations has allowed us to detect well-defined bumps in the optical LCs, as well as a change in colour, that are indicative of light coming from a core-collapse SN. The accompanying SNe for both events are individually compared with SN1998bw, a known GRB-supernova, and SN1994I, a typical type Ic supernova without a known GRB counterpart, and in both cases the brightness and temporal evolution more closely resemble SN1998bw. We also exploit our extensive optical and radio data for GRB 090618, as well as the publicly-available SWIFT -XRT data, and discuss the properties of the afterglow at early times. In the context of a simple jet-like model, the afterglow of GRB 090618 is best explained by the presence of a jet-break at t-to > 0.5 days. We then compare the rest-frame, peak V -band absolute magnitudes of all of the GRB and X-Ray Flash (XRF)-associated SNe with a large sample of local type Ibc SNe, concluding that, when host extinction is considered, the peak magnitudes of the GRB/XRF-SNe cannot be distinguished from the peak magnitudes of non-GRB/XRF SNe.

We present ground-based and HST optical observations of the optical transients (OTs) of long-duration Gamma Ray Bursts (GRBs) 060729 and 090618, both at a redshift of z = 0.54. For GRB 060729, bumps are seen in the optical light curves (LCs), and the late-time broadband spectral energy distributions (SEDs) of the OT resemble those of local type Ic supernovae (SNe). For GRB 090618, the dense sampling of our optical observations has allowed us to detect well-defined bumps in the optical LCs, as well as a change in colour, that are indicative of light coming from a core-collapse SN. The accompanying SNe for both events are individually compared with SN1998bw, a known GRB-supernova, and SN1994I, a typical type Ic supernova without a known GRB counterpart, and in both cases the brightness and temporal evolution more closely resemble SN1998bw. We also exploit our extensive optical and radio data for GRB 090618, as well as the publicly-available SWIFT -XRT data, and discuss the properties of the afterglow at early times. In the context of a simple jet-like model, the afterglow of GRB 090618 is best explained by the presence of a jet-break at t-to > 0.5 days. We then compare the rest-frame, peak V -band absolute magnitudes of all of the GRB and X-Ray Flash (XRF)-associated SNe with a large sample of local type Ibc SNe, concluding that, when host extinction is considered, the peak magnitudes of the GRB/XRF-SNe cannot be distinguished from the peak magnitudes of non-GRB/XRF SNe.

We present our analysis of the long cadence Kepler data for the well-studied Blazhko star RR Lyr, gathered during the first two quarters of the satellite's observations and covering a total of 127d. Besides being of great importance for our understanding of RR Lyrae stars in general, these RR Lyr data can be regarded as a case study for observations of bright stars with Kepler. Kepler can perform high-precision photometry on targets like RR Lyr, as the saturated flux is conserved to a very high degree. The Kepler data on RR Lyr are revolutionary in several respects. Even with long-cadence sampling (one measurement per 29.4 min), the unprecedented precision (< mmag) of the Kepler photometry allows the study of the star's extreme light curve variations in detail. The multiplet structures at the main frequency and its harmonics, typical for Blazhko stars, are clearly detected up to the quintuplets. For the first time, photometric data of RR Lyr reveal the presence of half-integer frequencies, linked to a period doubling effect. This phenomenon may be connected to the still unexplained Blazhko modulation. Moreover, with three observed Blazhko cycles at our disposal, we observe that there is no exact repetition in the light curve changes from one modulation cycle to the next for RR Lyr. This may be due to additional periodicities in the star, or to transient or quasi-periodic changes.

We present the source catalog and the properties of the $B-, R-$, and $I-$band images obtained to support the {\it AKARI} North Ecliptic Pole Wide (NEP-Wide) survey. The NEP-Wide is an {\it AKARI} infrared imaging survey of the north ecliptic pole covering a 5.8 deg$^2$ area over 2.5 -- 6 $\micron$ wavelengths. The optical imaging data were obtained at the Maidanak Observatory in Uzbekistan using the Seoul National University 4k $\times$ 4k Camera on the 1.5m telescope. These images cover 4.9 deg$^2$ where no deep optical imaging data are available. Our $B-, R-$, and $I-$band data reach the depths of $\sim$23.4, $\sim$23.1, and $\sim$22.3 mag (AB) at 5$\sigma$, respectively. The source catalog contains 96,460 objects in the $R-$band, and the astrometric accuracy is about 0.15$\arcsec$ at 1$\sigma$ in each RA and Dec direction. These photometric data will be useful for many studies including identification of optical counterparts of the infrared sources detected by {\it AKARI}, analysis of their spectral energy distributions from optical through infrared, and the selection of interesting objects to understand the obscured galaxy evolution.

GRB 090426 is a short duration burst detected by Swift ($T_{90}\sim 1.28$ s in the observer frame, and $T_{90}\sim 0.33$ s in the burst frame at $z=2.609$). Its host galaxy properties and some $\gamma$-ray related correlations are analogous to those seen in long duration GRBs, which are believed to be of a massive-star origin (so-called Type II GRBs). We present the results of its early optical observations with the 0.8-m TNT telescope at Xinglong observatory, and the 1-m LOAO telescope at Mt. Lemmon Optical Astronomy Observatory in Arizona. Our well-sampled optical afterglow lightcurve covers from $\sim 90$ seconds to $\sim 10^4$ seconds post the GRB trigger. It shows two shallow decay episodes that are likely due to energy injection, which end at $\sim 230$ seconds and $\sim 7100$ seconds, respectively. The decay slopes post the injection phases are consistent with each other ($\alpha\simeq 1.22$). The X-ray afterglow lightcurve appears to trace the optical, although the second energy injection phase was missed due to visibility constraints introduced by the {\em Swift} orbit. The X-ray spectral index is $\beta_X\sim 1.0$ without temporal evolution. Its decay slope is consistent with the prediction of the forward shock model. Both X-ray and optical emission is consistent with being in the same spectral regime above the cooling frequency ($\nu_c$). The fact that $\nu_c$ is below the optical band from the very early epoch of the observation provides a constraint on the burst environment, which is similar to that seen in classical long duration GRBs. We therefore suggest that death of a massive star is the possible progenitor of this short burst.

The first detection of the period doubling phenomenon is reported in the Kepler RR Lyrae stars RR Lyr, V808 Cyg and V355 Lyr. Interestingly, all these pulsating stars show Blazhko modulation. The period doubling manifests itself as alternating maxima and minima of the pulsational cycles in the light curve, as well as through the appearance of half-integer frequencies located halfway between the main pulsation period and its harmonics in the frequency spectrum. The effect was found to be stronger during certain phases of the modulation cycle. We were able to reproduce the period doubling bifurcation in our nonlinear RR Lyrae models computed by the Florida-Budapest hydrocode. This enabled us to trace the origin of this instability in RR Lyrae stars to a resonance, namely a 9:2 resonance between the fundamental mode and a high-order (9th) radial overtone showing strange-mode characteristics. We discuss the connection of this new type of variation to the mysterious Blazhko effect and argue that it may give us fresh insights to solve this century-old enigma.

We have observed a new cataclysmic variable (CV) SDSS J080434.20+510349.2 and study the origin of a long-term variability found in its light curve. Multi-longitude time-resolved photometric observations were carried out to analyze the uncommon behavior also found recently in two newly discovered CVs. This study of SDSS J080434.20+510349.2 mainly concerns the understanding of the nature of the observed double-humped light curve and its relation to a cyclic brightening occurring during quiescence. The observations were obtained early in 2007, when the object was at about V~17.1, 0.4 mag brighter than the pre-outburst magnitude. The light curve shows a sinusoidal variability with an amplitude of about 0.07 mag and a periodicity of 42.48 min, which is half of the orbital period of the system. In addition, we have observed two "mini-outbursts" of the system up to 0.6 mag, with a duration of about 4 days each. The "mini-outburst" had a symmetric profile and repeated in about 32 days. Subsequent monitoring of the system shows a cyclical behaviour of such "mini-outbursts" with a similar recurrence period. The origin of the double-humped light curve and the periodic brightening is discussed in the light of the evolutionary state of SDSS J080434.20+510349.2.

Results of a multi-color study of the variability of the magnetic cataclysmic variable BY Cam are presented. The observations were obtained at the Korean 1.8m and Ukrainian 2.6m, 1.2m and 38-cm telescopes in 2003-2005, 56 observational runs cover 189 hours. The variations of the mean brightness in different colors are correlated with a slope dR/dV=1.29(4), where the number in brackets denotes the error estimates in the last digits. For individual runs, this slope is much smaller ranging from 0.98(3) to 1.24(3), with a mean value of 1.11(1). Near the maximum, the slope becomes smaller for some nights, indicating more blue spectral energy distribution, whereas the night-to-night variability has an infrared character. For the simultaneous UBVRI photometry, the slopes increase with wavelength from dU/dR=0.23(1) to dI/dR=1.18(1). Such wavelength dependence is opposite to that observed in non-magnetic cataclysmic variables, in an agreement to the model of cyclotron emission. The principal component analysis shows two (with a third at the limit of detection) components of variablitity with different spectral energy distribution, which possibly correspond to different regions of emission. The scalegram analysis shows a highest peak corresponding to the 200-min spin variability, its quarter and to the 30-min and 8-min QPOs. The amplitudes of all these components are dependent on wavelength and luminosity state. The light curves were fitted by a statistically optimal trigonometrical polynomial (up to 4-th order) to take into account a 4-hump structure. The dependences of these parameters on the phase of the beat period and on mean brightness are discussed. The amplitude of spin variations increases with an increasing wavelength and with decreasing brightness

V-band time-series CCD photometric observations of the intermediate-age open cluster M11 were performed to search for variable stars. Using these time-series data, we carefully examined light variations of all stars in the observing field. A total of 82 variable stars were discovered, of which 39 stars had been detected recently by Hargis et al. (2005). On the basis of observational properties such as variable period, light curve shape, and position on a color-magnitude diagram, we classified their variable types as 11 delta Scuti-type pulsating stars, 2 gamma Doradus-type pulsating stars, 40 W UMa-type contact eclipsing binaries, 13 Algol-type detached eclipsing binaries, and 16 eclipsing binaries with long period. Cluster membership for each variable star was deduced from the previous proper motion results (McNamara et al. 1977) and position on the color-magnitude diagram. Many pulsating stars and eclipsing binaries in the region of M11 are probable members of the cluster.

Time-series CCD photometric observations of the intermediate-age open cluster NGC 2099 were performed to search for variable stars. We also carried out BV photometry to study physical properties of variables in the cluster. Using V-band time-series data, we carefully examined light variations of about 12,000 stars in the range of 10 < V < 22 mag. A total of 24 variable stars have been identified; seven stars are previously known variables and 17 stars are newly identified. On the basis of observational properties such as light curve shape, period, and amplitude, we classified the new variable stars as nine delta Scuti-type pulsating stars, seven eclipsing binaries, and one peculiar variable star. Judging from the position of delta Scuti-type stars in the color-magnitude diagram, only two stars are likely to have the cluster membership. One new variable KV10 shows peculiar light variations with a delta Scuti-type short period of about 0.044 day as well as a long period of 0.417 day.

The Angstrom Project is using a distributed network of two-meter class telescopes to conduct a high cadence pixel-lensing survey of the bulge of the Andromeda Galaxy (M31). With the expansion of global telescope network, the detection efficiency of pixel-lensing surveys is rapidly improving. In this paper, we estimate the detection rate of binary lens events expected from high-cadence pixel-lensing surveys toward M31 such as the Angstrom Project based on detailed simulation of events and application of realistic observational conditions. Under the conservative detection criteria that only high signal-to-noise caustic-crossing events with long enough durations between caustic crossings can be firmly identified as binary lens events, we estimate that the rate would be $\Gamma_{\rm b}\sim (7-15)f_{\rm b}(N/50)$ per season, where $f_{\rm b}$ is the fraction of binaries with projected separations of $10^{-3} {\rm AU}<\tilde{d}<10^3 {\rm AU}$ out of all lenses and $N$ is the rate of stellar pixel-lensing events. We find that detected binaries would have mass ratios distributed over a wide range of $q\gtrsim 0.1$ but with separations populated within a narrow range of $1 {\rm AU}\lesssim \tilde{d}\lesssim 5 {\rm AU}$. Implementation of an alert system and subsequent follow-up observations would be important not only for the increase of the binary lens event rate but also for the characterization of lens matter.

The Angstrom Project is undertaking an optical survey of stellar microlensing events across the bulge region of the Andromeda Galaxy (M31) using a distributed network of two-meter class telescopes. The Angstrom Project Alert System (APAS) has been developed to identify in real time candidate microlensing and transient events using data from the Liverpool and Faulkes North robotic telescopes. This is the first time that real-time microlensing discovery has been attempted outside of the Milky Way and its satellite galaxies. The APAS is designed to enable follow-up studies of M31 microlensing systems, including searches for gas giant planets in M31. Here we describe the APAS and we present a few example light curves obtained during its commissioning phase which clearly demonstrate its real-time capability to identify microlensing candidates as well as other transient sources.

Microlensing is one of the most powerful methods that can detect extrasolar planets and a future space-based survey with a high monitoring frequency is proposed to detect a large sample of Earth-mass planets. In this paper, we examine the sensitivity of the future microlensing survey to Earth-mass planets located in the habitable zone. For this, we estimate the fraction of Earth-mass planets that will be located in the habitable zone of their parent stars by carrying out detailed simulation of microlensing events based on standard models of the physical and dynamic distributions and the mass function of Galactic matter. From this investigation, we find that among the total detectable Earth-mass planets from the survey, those located in the habitable zone would comprise less than 1% even under a less-conservative definition of the habitable zone. We find the main reason for the low sensitivity is that the projected star-planet separation at which the microlensing planet detection efficiency becomes maximum (lensing zone) is in most cases substantially larger than the median value of the habitable zone. We find that the ratio of the median radius of the habitable zone to the mean radius of the lensing zone is roughly expressed as $d_{\rm HZ}/r_{\rm E}\sim 0.2(m/0.5 M_\odot)^{1/2}$.

We present BV CCD photometry for the open clusters Czernik 24 and Czernik 27. These clusters have never been studied before, and we provide, for the first time, the cluster parameters; reddening, distance, metallicity and age. Czernik 24 is an old open cluster with age 1.8 +/- 0.2 Gyr, metallicity [Fe/H]=-0.41 +/- 0.15 dex, distance modulus (m-M)_0 = 13.1 +/- 0.3 mag (d=4.1 +/- 0.5 kpc), and reddening E(B-V) = 0.54 +/- 0.12 mag. The parameters for Czernik 27 are estimated to be age = 0.63 +/- 0.07 Gyr, [Fe/H]= -0.02 +/- 0.10 dex, (m-M)_0 = 13.8 +/- 0.2 mag (d=5.8 +/- 0.5 kpc), and E(B-V) = 0.15 +/- 0.05 mag. The metallicity and distance values for Czernik 24 are consistent with the relation between the metallicity and the Galactocentric distance of other old open clusters. We find the metallicity gradient of 51 old open clusters including Czernik 24 to be Delta [Fe/H]/Delta R_gc= -0.064 +/- 0.009 dex/kpc.

The Angstrom Project is using a global network of 2m-class telescopes to conduct a high cadence pixel microlensing survey of the bulge of the Andromeda Galaxy (M31), with the primary aim of constraining its underlying bulge mass distribution and stellar mass function. Here we investigate the feasibility of using such a survey to detect planets in M31. We estimate the efficiency of detecting signals for events induced by planetary systems as a function of planet/star mass ratio and separation, source type and background M31 surface brightness. We find that for planets of a Jupiter-mass or above that are within the lensing zone (~1 -3 AU) detection is possible above 3 $\sigma$, with detection efficiencies ~3% for events associated with giant stars, which are the typical source stars of pixel-lensing surveys. A dramatic improvement in the efficiency of ~40 -- 60% is expected if follow-up observations on an 8m telescope are made possible by a real-time alert system.

Results of a CCD study of the variability of the cataclysmic variable BG CMi obtained at the Korean 1.8m telescope in 2002-2005 are presented. The "multi-comparison star" method had been applied for better accuracy estimates. The linear ephemeris based on 19 mean maxima for 2002--2005 is HJD 2453105.31448(6)+0.01057257716(198)(E-764707). The period differs from that predicted by the quadratic ephemeris by Pych et al. (1996) leading to a possible cycle miscount. The statistically optimal ephemeris is a fourth-order polynomial, as a quadratic or even a cubic ephemeris leads to unaceptably large residuals: Min.HJD=$ 2445020.28095(28)+0.0105729609(57)E -1.58(32)\cdot10^{-13}E^2-5.81(64)\cdot10^{-19}E^3+4.92(41)\cdot10^{-25}E^4.$ Thus the rate of the spin-up of the white dwarf is decreasing. An alternative explanation is that the spin-up has been stopped during recent years. The deviations between the amplutudes of the spin variability in V and R, as well as between phases are not statistically significant. However, the orbital light curves exhibit distinct difference; the corresponding color index shows a nearly sinusoidal shape with a maximum at orbital phase $\sim0.2.$ The variations of the amplitude of spin waves shows a short maximum at the phase of the orbital dip. The corrected ephemeris for orbital minima is Min.HJD=$2448368.7225(12)+0\fd13474841(6)\cdot(E-24849)$ with a narrow dip occuring 0.07P later. The rate of the spin period variation seems to be changed, justifying the necessity of regular observations of intermediate polars.

To maximize the number of planet detections, current microlensing follow-up observations are focusing on high-magnification events which have a higher chance of being perturbed by central caustics. In this paper, we investigate the properties of central caustics and the perturbations induced by them. We derive analytic expressions of the location, size, and shape of the central caustic as a function of the star-planet separation, $s$, and the planet/star mass ratio, $q$, under the planetary perturbative approximation and compare the results with those based on numerical computations. While it has been known that the size of the planetary caustic is \propto \sqrt{q}, we find from this work that the dependence of the size of the central caustic on $q$ is linear, i.e., \propto q, implying that the central caustic shrinks much more rapidly with the decrease of $q$ compared to the planetary caustic. The central-caustic size depends also on the star-planet separation. If the size of the caustic is defined as the separation between the two cusps on the star-planet axis (horizontal width), we find that the dependence of the central-caustic size on the separation is \propto (s+1/s). While the size of the central caustic depends both on $s$ and q, its shape defined as the vertical/horizontal width ratio, R_c, is solely dependent on the planetary separation and we derive an analytic relation between R_c and s. Due to the smaller size of the central caustic combined with much more rapid decrease of its size with the decrease of q, the effect of finite source size on the perturbation induced by the central caustic is much more severe than the effect on the perturbation induced by the planetary caustic. Abridged.

The Andromeda Galaxy Stellar Robotic Microlensing Project (The Angstrom Project) aims to use stellar microlensing events to trace the structure and composition of the inner regions of the Andromeda Galaxy (M31). We present microlensing rate and timescale predictions and spatial distributions for stellar and sub-stellar lens populations in combined disk and barred bulge models of M31. We show that at least half of the stellar microlenses in and around the bulge are expected to have characteristic durations between 1 and 10 days, rising to as much as 80% for brown-dwarf dominated mass functions. These short-duration events are mostly missed by current microlensing surveys that are looking for Macho candidates in the M31 dark matter halo. Our models predict that an intensive monitoring survey programme such as Angstrom, which will be able to detect events of durations upwards of a day, could detect around 30 events per season within ~5 arcminutes of the M31 centre, due to ordinary low-mass stars and remnants. This yield increases to more than 60 events for brown-dwarf dominated mass functions. The overall number of events and their average duration are sensitive diagnostics of the bulge mass, in particular the contribution of low-mass stars and brown dwarfs. The combination of an inclined disk, an offset bar-like bulge, and differences in the bulge and disk luminosity functions results in a four-way asymmetry in the number of events expected in each quadrant defined by the M31 disk axes. The asymmetry is sensitive to the bar prolongation, orientation and mass.

We carried out wide-field BVI CCD photometric observations of the GCs M3 and M13 using the BOAO 1.8 m telescope equipped with a 2K CCD. We present CMDs of M3 and M13. We have found AGB bumps at V = 14.85 for M3 at V = 14.25 for M13. It is found that AGB stars in M3 are more concentrated near the bump, while those in M13 are scattered along the AGB sequence. We identified the RGB bump of M3 at V = 15.50 and that of M13 at V = 14.80. We have estimated the ratios R and R2 for M3 and M13 and found that of R for M3 is larger than that for M13 while R2's for M3 and M13 are similar when only normal HB stars are used in R and R2 for M13. However, we found that R's for M3 and M13 are similar while R2 for M3 is larger than that for M13 when all the HB stars are included in R and R2 for M13. We have compared the observed RGB LFs of M3 and M13 with the theoretical RGB LF of Bergbusch & VandenBerg at the same radial distances from the cluster centers as used in R and R2 for M3 and M13. We found "extra stars" belonging to M13 in the comparison of the observed RGB LF of M13 and the theoretical RGB LF of Bergbusch & VandenBerg. In the original definition of R of Buzzoni et al., N(HB) corresponds to the lifetime of HB stars in the RR Lyrae instability strip at log T_eff = 3.85. So, the smaller R value resulting for M13 compared with that for M3 in the case where only normal HB stars are included in R and R2 for M13 may be partially caused by "extra stars", and the similar R's for M3 and M13 in the case where the all HB stars are included in R and R2 for M13 may be caused by "extra stars" in the upper RGB of M13. If "extra stars" in the upper RGB of M13 are caused by an effective "deep mixing" these facts support the contention that an effective "deep mixing" could lead to different HB morphologies between M3 and M13 and subsequent sequences.

Through time-series CCD photometry of the globular cluster NGC 5466, we have detected nine SX Phoenicis stars including three new ones. All the SX Phoenicis stars are located in the blue straggler region in the color-magnitude diagram of NGC 5466. Five of them show clearly double-radial mode features, the periods of which are well consistent with the theoretical ratio of the first overtone mode to the fundamental mode (P$_{1H}$/P$_F$). Normally, it has not been easy to secure a P-L relation of the SX Phoenicis stars, because determination of the pulsational mode of the SX Phoenicis stars has been difficult. The existence of five SX Phoenicis stars in NGC5466 with double-radial modes allows us to derive reliably a P-L relation for the fundamental mode of the SX Phoenicis stars. Using seven SX Phoenicis stars including five stars with double-radial modes, we derive a P-L relation for the fundamental mode in NGC 5466, $<V> = - 3.25 (\pm 0.46) Log P + 14.70 (\pm 0.06), (\sigma = \pm0.04) $, corresponding to $<M_V> = - 3.25 (\pm0.46)Log P - 1.30 (\pm0.06) $ for an adopted distance modulus of $(m-M)_0=16.00$ and zero reddening.

We present deep wide field $VI$ CCD photometry of the Sextans dwarf spheroidal galaxy (dSph) in the Local Group, covering a field of 42' x 28' located at the center of the galaxy. Color-magnitude diagrams of the Sextans dSph show well-defined red giant branch (RGB), blue horizontal branch (BHB), prominent red horizontal branch (RHB), asymptotic giant branch (AGB), about 120 variable star candidates including RR Lyraes and anomalous Cepheids, about 230 blue stragglers (BSs), and main sequence (MS) stars. The distance to the galaxy is derived using the $I$-band magnitude of the tip of the RGB at I(TRGB)=15.95+/-0.04: (m-M)_0=19.90+/-0.06 for an adopted reddening of E(B-V)=0.01. The mean metallicity of the RGB is estimated from the (V-I) color: [Fe/H]=-2.1+/-0.1(statistical error) +/-0.2(standard calibration error) dex, with a dispersion of sigma[Fe/H]=0.2 dex. The age of the MSTO of the main old population is estimated to be similar to that of M92, and there are seen some stellar populations with younger age. There is found to be one RGB bump at V=19.95+/-0.05 mag (M_V=0.03 mag), and a weak brighter bump at V=19.35+/-0.05 mag(Mv=-0.58 mag) which is probably an AGB bump. The V-band luminosity function of the RGB and MS stars is in general similar to that of M92, with a slight excess of stars in the magnitude range brighter than the MSTO with respect to that of M92. The bright BSs are more centrally concentrated than the faint BSs. Significant radial gradients are seen for several kinds of populations: the RHB, the red RGB, the red subgiant branch (SGB), and the bright BSs are more centrally concentrated toward the center of the galaxy, compared with the BHB, the blue RGB, the blue SGB, and the faint BSs, respectively.

Through time-series CCD photometry of the metal-poor globular cluster M53, we have discovered eight new SX Phoenicis type stars (labeled from SXP1 to SXP8). All the new SX Phoenicis stars are located in the blue straggler star region of a color-magnitude diagram of M53. One of these stars (SXP2) is found to have very closely separated pulsation frequencies: $f_1/f_2 = 0.9595$ where $f_1$ and $f_2$ are primary and secondary frequencies. This may be due to excitation of non-radial modes. Six of these SX Phoenicis stars are considered to be pulsating in the fundamental mode. They show a tight linear correlation between the period and luminosity. We derive a period - luminosity relation for the fundamental mode for the period range of $-1.36 < Log P[d]< -1.15$ : $<V>=-3.010(\pm0.262)Log P + 15.310(\pm0.048)$ with an rms scatter of 0.038, corresponding to $<M_V>=-3.010 Log P - 1.070$ for an adopted distance modulus of $(m-M)_V=16.38$ (Harris 1996).

We report on an observing campaign in March 2001 to monitor the brightness of the later arriving Q0957+561 B image in order to compare with the previously published brightness observations of the (first arriving) A image. The 12 participating observatories provided 3543 image frames which we have analyzed for brightness fluctuations. From our classical methods for time delay determination, we find a 417.09 +/- 0.07 day time delay which should be free of effects due to incomplete sampling. During the campaign period, the quasar brightness was relatively constant and only small fluctuations were found; we compare the structure function for the new data with structure function estimates for the 1995--6 epoch, and show that the structure function is statistically non-stationary. We also examine the data for any evidence of correlated fluctuations at zero lag. We discuss the limits to our ability to measure the cosmological time delay if the quasar's emitting surface is time resolved, as seems likely.

We have initiated a long-term project, the BOAO photometric survey of open clusters, to enlarge our understanding of galactic structure using UBVI CCD photometry of open clusters which have been little studied before. This is the second paper of the project in which we present the photometry of 12 open clusters. We have determined the cluster parameters by fitting the Padova isochrones to the color-magnitude diagrams of the clusters. All the clusters except for Be 0 and NGC 1348 are found to be intermediate-age to old (0.2 - 4.0 Gyrs) open clusters with a mean metallicity of [Fe/H] = 0.0.

A new SX Phe star (labelled SXP1) found from $BV$ CCD photometry is the first to be discovered in the globular cluster M15. It is a blue straggler and is located $102\arcsec.8$ north and $285\arcsec.6$ west of the center of M15 \citep{har96}. Mean magnitudes of SXP1 are $<$$B$$>$ = 18$\fm$671 and $<$$V$$>$ = 18$\fm$445. The amplitude of variability of SXP1 is measured to be $\Delta V \approx 0.15$. From multiple-frequency analysis based on the Fourier decomposition method, we detect two very closely separated pulsating frequencies: the primary frequency at $f_1=24.630$ c/d for both $B$- and $V$-bands, and the secondary frequency at $f_2=24.338$ c/d for the $B$-band and 24.343 c/d for the $V$-band. This star is the second among known SX Phe stars found to pulsate with very closely separated frequencies ($f_2/f_1\ge0.95$). These frequencies may be explained by excitation of nonradial modes; however, we have an incomplete understanding of this phenomenon in the case of SX Phe stars with relatively high amplitudes. The metallicity-period and the variability amplitude-period relations for SXP1 in M15 are found to be consistent with those for SX Phe stars in other globular clusters.

Open clusters are useful tools to investigate the structure and evolution of the Galactic disk. We have started a long-term project to obtain UBVI CCD Photometry of open clusters which were little studied before, using the Doyak 1.8 m telescope of Bohyunsan Optical Astronomy Observatory in Korea. The primary goals of this project are (1) to make a catalog of UBVI photometry of open clusters, (2) to make an atlas of open clusters, and (3) to survey and monitor variable stars in open clusters. Here we describe this project and report the first results based on preliminary analysis of the data on four open clusters in the survey sample: Be 14, Cr 74, Biu 9, and NGC 2355. Isochrone fitting of the color-magnitude diagrams of the clusters shows that all of them are intermediate age to old (0.3 - 1.6 Gyrs) open clusters with moderate metallicity.