Source author record

Jae Woo Lee

Jae Woo Lee appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

49works
13topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

49 published item(s)

preprint2022arXiv

Stability and selective extinction in complex mutualistic networks

We study species abundance in the empirical plant-pollinator mutualistic networks exhibiting broad degree distributions, with uniform intra-group competition assumed, by the Lotka-Volterra equation. The stability of a fixed point is found to be identified by the signs of its non-zero components and those of its neighboring fixed points. Taking the annealed approximation, we derive the non-zero components to be formulated in terms of degrees and the rescaled interaction strengths, which lead us to find different stable fixed points depending on parameters, and we obtain the phase diagram. The selective extinction phase finds small-degree species extinct and effective interaction reduced, maintaining stability and hindering the onset of instability. The non-zero minimum species abundances from different empirical networks show data collapse when rescaled as predicted theoretically.

preprint2022arXiv

The Pre-He White Dwarfs in Eclipsing Binaries. III. WASP 1625-04

1SWASP J162545.15-043027.9 (WASP 1625-04) has been announced as one of EL CVn candidates showing total primary eclipses and ellipsoidal variations. This paper presents the absolute properties of the binary star, based on our high-resolution spectroscopy conducted from 2015 through 2020. From the spectral analysis, the radial velocities (RVs) for both components were obtained with the effective temperature $T_{\rm eff,1} = 8990 \pm 200$ K and the rotational rate $v_1\sin$$i=53\pm5$ km s$^{-1}$ for the more massive primary. The RV measurements were analyzed with archival WASP photometry. From the modelling we obtained: $M_1 = 1.745 \pm 0.013$ M$_\odot$, $M_2 = 0.187 \pm 0.002$ M$_\odot$, $R_1 = 1.626 \pm 0.008$ M$_\odot$, $R_2 = 0.290 \pm 0.003$ M$_\odot$, $L_1 = 15.5 \pm 1.4$ L$_\odot$, and $L_2 = 1.84 \pm 0.16$ L$_\odot$. In the Hertzsprung-Russell diagram, WASP 1625-04 A lies on the zero-age main sequence and its companion accords well with the helium-core white dwarf models of 0.19 M$_\odot$ in the constant luminosity phase. Our improved results demonstrate that WASP 1625-04 is a typical EL CVn-type binary with a low mass ratio and $M_2$ combination in the thin-disk population and is the product of the stable, non-conservative mass transfer of the precursor binary.

preprint2022arXiv

The Pre-He White Dwarfs in Eclipsing Binaries. IV. WASP 1814+48 with Multiperiodic Pulsations

For the EL CVn candidate 1SWASPJ181417.43+481117.0 (WASP 1814+48), we secured the first spectroscopic observations between 2015 April and 2021 March. Using the echelle spectra, the radial velocities (RVs) of the primary star were measured with its atmospheric parameters of $T_{\rm eff,1}=7770\pm130$ K and $v_1$$\sin$$i=47\pm6$ km s$^{-1}$. We fitted our single-lined RVs and the TESS light curve simultaneously. From the binary modeling, we determined the following fundamental parameters for each component: $M_1=1.659\pm0.048$ $M_\odot$, $R_1=1.945\pm0.027$ $R_\odot$, and $L_1=12.35\pm0.90$ $L_\odot$ for WASP 1814+48 A, and $M_2=0.172\pm0.005$ $M_\odot$, $R_2=0.194\pm0.005$ $R_\odot$, and $L_2=0.69\pm0.07$ $L_\odot$ for WASP 1814+48 B. The surface gravity of $\log g_2=5.098\pm0.026$ obtained from $M_2$ and $R_2$ is concurrent with 5.097$\pm$0.025 computed directly from the observable quantities. WASP 1814+48 B is well-matched with the 0.176 M$_\odot$ white dwarf (WD) evolutionary model for $Z=0.01$. The metallicity and our Galactic kinematics indicate that the program target is a thin-disk star. The whole light residuals after removal of the binary trend were analyzed and found to oscillate at a total of 52 frequencies. Among these, most of the low frequencies below 24 day$^{-1}$ are aliases and orbital harmonics. The five significant frequencies between 32 and 36 day$^{-1}$ are the pulsation modes of WASP 1814+48 A located in the $δ$ Sct domain on ZAMS, and the high frequencies of 128$-$288 day$^{-1}$ arise from WASP 1814+48 B in the pre-He WD instability strip. Our results reveal that WASP 1814+48 is the fifth EL CVn star that is composed of a $δ$ Sct-type primary and a pre-ELMV (extremely low-mass pre-He WD variable).

preprint2022arXiv

V608 Cassiopeiae: A W UMa-type eclipsing binary with two possible circumbinary companions

We present the photometric properties of V608 Cas from detailed studies of light curves and eclipse timings. The light curve synthesis indicates that the eclipsing pair is an overcontact binary with parameters of $ΔT$ = 155 K, $q$ = 0.328, and $f$ = 26 %. We detected the third light $\ell_{3}$, which corresponds to about 8 % and 5 % of the total systemic light in $V$ and $R$ bands, respectively. Including our 6 timing measurements, a total of 38 times of minimum light were used for a period study. It was found that the orbital period of V608 Cas has varied in some combination of an upward parabola and two periodic variations. The continuous period increase with a rate of $+$3.99$\times$10$^{-7}$ d yr$^{-1}$ can be interpreted as a mass transfer from the secondary component to the primary star at a rate of 1.51$\times$10$^{-7}$ M$_\odot$ yr$^{-1}$. The periods and semi-amplitudes of the two periodic variations are about $P_3$ = 16.0 yr and $P_4$ = 26.3 yr, and $K_3$ = 0.0341 d and $K_4$ = 0.0305 d, respectively. The most likely explanation of both cycles is a pair of light-traveling time effects operated by the possible presence of third and fourth components with estimated masses of $M_3$ = 2.20 M$_\odot$ and $M_4$ = 1.27 M$_\odot$ in eccentric orbits of $e_3$ = 0.66 and $e_4$ = 0.52. Because the contribution of $\ell_{3}$ is very low compared to the estimated masses of two circumbinary objects, they can be inferred as very faint compact objects.

preprint2021arXiv

TESS Eclipsing Binary Stars. I. Short cadence observations of 4584 eclipsing binaries in Sectors 1-26

In this paper we present a catalog of 4584 eclipsing binaries observed during the first two years (26 sectors) of the TESS survey. We discuss selection criteria for eclipsing binary candidates, detection of hither-to unknown eclipsing systems, determination of the ephemerides, the validation and triage process, and the derivation of heuristic estimates for the ephemerides. Instead of keeping to the widely used discrete classes, we propose a binary star morphology classification based on a dimensionality reduction algorithm. Finally, we present statistical properties of the sample, we qualitatively estimate completeness, and discuss the results. The work presented here is organized and performed within the TESS Eclipsing Binary Working Group, an open group of professional and citizen scientists; we conclude by describing ongoing work and future goals for the group. The catalog is available from http://tessEBs.villanova.edu and from MAST.

preprint2020arXiv

Estimation of Infection Rate and Prediction of Initial Infected Individuals of COVID-19

We consider the pandemic spreading of COVID-19 for some selected countries after the outbreak of the coronavirus in Wuhan City, China. We estimated the infection rate and the initial infected individuals of COVID-19 by using the officially reported data at the early stage of the epidemic for the susceptible (S), infectable (I), quarantined (Q), and the cofirmed recovered (Rk) population model, so called SIQRk model. In the reported data we know the quarantined cases and the recovered cases. We can not know the recovered cases from the asymptomatic cases. In the SIQRk model we can estimated the model parameters and the initial infecting cases (confirmed ans asymtomatic cases) from the data fits. We obtained the infection rate in the range between 0.233 and 0.462, the basic reproduction number Ro in the range between 1.8 and 3.5, and the initial number of infected individuals in the range betwee 10 and 8409 for some selected countries. By using fitting parameters we estimated the maximum time of the infection for Germany when the government are performing the quarantine policy. The disease is undergoing to the calm state about six months after first patients were identified.

preprint2020arXiv

TESS Photometry of the Eclipsing $δ$ Scuti Star AI Hydrae

AI Hya has been known as an eclipsing binary with a monoperiodic $δ$ Sct pulsator. We present the results from its {\it TESS} photometry observed during Sector 7. Including our five minimum epochs, the eclipse timing diagram displays the apsidal motion with a rate of $\dotω$ = 0.075$\pm$0.031 deg year$^{-1}$, which corresponds to an apsidal period of U = 4800$\pm$2000 years. The binary star model represents that the smaller, less massive primary component is 427 K hotter than the pulsating secondary, and our distance of 612$\pm$36 pc is in good agreement with the $Gaia$ distance of 644$\pm$26 pc. We subtracted the binary effects from the observed {\it TESS} data and applied a multifrequency analysis to these residuals. The result reveals that AI Hya is multiperiodic in its pulsation. Of 14 signals detected, four ($f_1$, $f_2$, $f_3$, $f_6$) may be considered independent pulsation frequencies. The period ratios of $P_{\rm pul}/P_{\rm orb}$ = 0.012$-$0.021 and the pulsation constants of $Q$ = 0.30$-$0.52 days correspond to $δ$ Sct pulsations in binaries. We found that the secondary component of AI Hya pulsates in both radial fundamental $F$ modes ($f_2$ and $f_3$) and non-radial $g_1$ modes with a low degree of $\ell$ = 2 ($f_1$ and $f_6$).

preprint2020arXiv

The Pre-He White Dwarfs in Eclipsing Binaries. I. WASP 0131+28

We report the first $BV$ light curves and high-resolution spectra of the post-mass transfer binary star WASP 0131+28 to study the absolute properties of extremely low-mass white dwarfs. From the observed spectra, the double-lined radial velocities were derived, and the effective temperature and rotational velocity of the brighter, more massive primary were found to be $T_{\rm eff,1} = 10,000 \pm 200$ K and $v_1\sin$$i$ = 55 $\pm$ 10 km s$^{-1}$, respectively. The combined analysis of the {\it TESS} archive data and ours yielded the accurate fundamental parameters of the program target. The masses were derived to about 1.0 \% accuracy and the radii to 0.6 \%, or better. The secondary component's parameters of $M_2 = 0.200 \pm 0.002$ M$_\odot$, $R_2 = 0.528 \pm 0.003$ R$_\odot$, $T_{\rm eff,2}$ = 11,186 $\pm$ 235 K, and $L_2 = 3.9 \pm 0.3$ L$_\odot$ are in excellent agreement with the evolutionary sequence for a helium-core white dwarf of mass 0.203 M$_\odot$, and indicates that this star is halfway through the constant luminosity phase. The results presented in this article demonstrate that WASP 0131+28 is an EL CVn eclipsing binary in a thin disk, which is formed from the stable Roche-lobe overflow channel and composed of a main-sequence dwarf with a spectral type A0 and a pre-He white dwarf.

preprint2020arXiv

Tidally Excited Modes and $δ$ Scuti Pulsations in the Eclipsing Triple Star IM Persei

IM Per is a triple star system whose eclipsing pair masses and radii are accurate to within 1 \%. The TESS light curve of the program target exhibits partial eclipses and multiple oscillations with mmag-level amplitudes. It is found that the oscillations affect eclipse timing measurements. Binary modeling of the high-quality TESS data indicates that the eclipsing components of the triple system are twin dwarfs with parameters of $M_2$/$M_1$ = 0.995, $R_2$/$R_1$ = 0.901, and $Δ$($T_{\rm eff,1}$--$T_{\rm eff,2}$) = 12 K in an eccentric ($e$ = 0.049), detached configuration. The third light of $l_3$ = 0.054 may mostly come from a G-type tertiary companion. Our predicted parallax of 1.52 $\pm$ 0.09 mas is concurrent with the Gaia measurement of 1.52 $\pm$ 0.05 mas. Multifrequency analysis of the outside-eclipse residuals reveals 22 significant pulsation signals: four in the gravity-mode region (0.03$-$2.22 day$^{-1}$) and 18 in the pressure-mode region (9.19$-$25.12 day$^{-1}$). Of the low frequencies, $f_{11}$ and $f_{14}$ are orbital harmonics that can be identified as tidally excited modes. The pulsation periods and constants for the high frequencies, and the position in the Cepheid instability strip demonstrate that the pulsating component of IM Per is a $δ$ Sct variable.

preprint2019arXiv

Absolute Properties of the Eclipsing $γ$ Dor Star V404 Lyrae

We present the first high-resolution spectra for the eclipsing binary V404 Lyr showing $γ$ Dor pulsations, which we use to study its absolute properties. By fitting models to the disentangling spectrum of the primary star, we found that it has an effective temperature of $T_{\rm eff,1}=7,330 \pm 150$ K and a rotational velocity of $v_1\sin$$i=148\pm18$ km s$^{-1}$. The simultaneous analysis of our double-lined radial velocities and the pulsation-subtracted ${\it Kepler}$ data gives us accurate stellar and system parameters of V404 Lyr. The masses, radii, and luminosities are $M_1$ = 2.17$\pm$0.06 M$_\odot$, $R_1$ = 1.91$\pm$0.02 R$_\odot$, and $L_1$ = 9.4$\pm$0.8 L$_\odot$ for the primary, and $M_2$ = 1.42$\pm$0.04 M$_\odot$, $R_2$ = 1.79$\pm$0.02 R$_\odot$, and $L_2$ = 2.9$\pm$0.2 L$_\odot$ for the secondary. The tertiary component orbiting the eclipsing pair has a mass of $M_{\rm 3b}$ = 0.71$\pm$0.15 $M_\odot$ in an orbit of $P_{\rm 3b}$ = 642$\pm$3 days, $e_{\rm 3b}$ = 0.21$\pm$0.04, and $a_{\rm 3b}$ = 509$\pm$2 R$_\odot$. The third light of $l_3=4.1\pm0.2\%$ could be partly attributable to the K-type circumbinary object. By applying a multiple frequency analysis to the eclipse-subtracted light residuals, we detected 45 frequencies with signal to noise amplitude ratios larger than 4.0. Identified as independent pulsation modes, seven frequencies ($f_1-f_6$, $f_9$), their new pulsation constants, and the location in the Hertzsprung-Russell diagram indicate that the pulsating primary is a $γ$ Dor-type variable star.

preprint2016arXiv

$γ$ Doradus Pulsations in the Eclipsing Binary Star KIC 6048106

We present the ${\it Kepler}$ photometry of KIC 6048106 exhibiting O'Connell effect and multiperiodic pulsations. Including a starspot on either of the components, light-curve synthesis indicates that this system is a semi-detached Algol with a mass ratio of 0.211, an orbital inclination of 73.9 deg, and a large temperature difference of 2,534 K. To examine in detail both spot variations and pulsations, we separately analyzed the {\it Kepler} time-series data at the interval of an orbital period by an iterative way. The results reveal that the variable asymmetries of the light maxima can be interpreted as the changes of a magnetic cool spot on the secondary component with time. Multiple frequency analyses were performed in the outside-eclipse light residuals after removal of the binarity effects from the observed {\it Kepler} data. We detected 30 frequencies with signal to noise amplitude ratios larger than 4.0, of which six ($f_2$--$f_6$ and $f_{10}$) can be identified as high-order (17 $\le n \le$ 25) low-degree ($\ell$ = 2) gravity-mode pulsations that were stable during the observing run of 200 d. In contrast, the other frequencies may be harmonic and combination terms. For the six frequencies, the pulsation periods and pulsation constants are in the ranges of 0.352$-$0.506 d and 0.232$-$0.333 d, respectively. These values and the position on the HR diagram demonstrate that the primary star is a $γ$ Dor variable. The evolutionary status and the pulsation nature of KIC 6048106 are discussed.

preprint2016arXiv

KIC 6220497: A New Algol-type Eclipsing Binary with Multiperiodic Pulsations

We present both binarity and pulsation of KIC 6220497 from the {\it Kepler} observations. The light curve synthesis shows that the eclipsing system is a semi-detached Algol with parameters of $q$ = 0.243$\pm$0.001, $i$ = 77.3$\pm$0.3 deg, and $ΔT$ = 3,372$\pm$58 K, in which the detached primary component fills its Roche lobe by $\sim$87\%. A multiple frequency analysis of the eclipse-subtracted light residuals reveals 33 frequencies in the range of 0.75$-$20.22 d$^{-1}$ with amplitudes between 0.27 and 4.56 mmag. Among these, four are pulsation frequencies in fundamental ($f_1$, $f_5$) and $p$ ($f_2$, $f_7$) modes, and six are orbital frequency ($f_8$, $f_{31}$) and its harmonics ($f_6$, $f_{11}$, $f_{20}$, $f_{24}$), which can be attributed to tidally excited modes. For the pulsation frequencies, the pulsation constants of 0.16$-$0.33 d and the period ratios of $P_{\rm pul}/P_{\rm orb}$ = 0.042$-$0.089 indicate that the primary component is a $δ$ Sct pulsating star and, thus, KIC 6220497 is an oscillating eclipsing Algol (oEA) star. The dominant pulsation period of 0.1174051$\pm$0.0000004 d is significantly longer than that expected from empirical relations that link the pulsation period with the orbital period. The surface gravity of $\log g_1$ = 3.78$\pm$0.03 is clearly smaller than those of the other oEA stars with similar orbital periods. The pulsation period and the surface gravity of the pulsating primary demonstrate that KIC 6220497 would be the more evolved EB, compared with normal oEA stars.

preprint2016arXiv

The Varying Light Curve and Timings of the Ultra-short Period Contact Binary KIC 9532219

KIC 9532219 is a W UMa-type eclipsing binary with an orbital period of 0.1981549 d that is below the short-period limit ($\sim$0.22 d) of the period distribution for contact binaries. The {\it Kepler} light curve of the system exhibits striking light changes in both eclipse depths and light maxima. Applying third-body and spot effects, the light-curve synthesis indicates that the eclipsing pair is currently in a marginal contact stage with a mass ratio of $q$=1.20, an orbital inclination of $i$=66.0 deg, a temperature difference of $Δ$ ($T_{1}$--$T_{2}$)=172 K, and a third light of $l_3$=75.9 \%. To understand the light variations with time, we divided up the light curve into 312 segments and separately analyzed them. The results reveal that variation of eclipse depth is primarily caused by changing amounts of contamination due to the nearby star KIC9532228 between the {\it Kepler} Quarters and that the variable O'Connell effect originates from the starspot activity on the less massive primary component. Based on our light-curve timings, a period study of KIC 9532219 indicates that the orbital period has varied as a combination of a downward parabola and a light-travel-time (LTT) effect due to a third body, with a period of 1196 d and a minimum mass of 0.0892 $M_\odot$ in an eccentric orbit of 0.150. The parabolic variation could be a small part of a second LTT orbit due to a fourth component in a wider orbit, instead of either mass transfer or angular momentum loss.

preprint2016arXiv

Time-series spectroscopy of the pulsating eclipsing binary XX Cephei

Oscillating Algol-type eclipsing binaries (oEA) are very interesting objects that have three observational features of eclipse, pulsation, and mass transfer. Direct measurement of their masses and radii from the double-lined radial velocity data and photometric light curves would be the most essential for understanding their evolutionary process and for performing the asteroseismological study. We present the physical properties of the oEA star XX Cep from high-resolution time-series spectroscopic data. The effective temperature of the primary star was determined to be 7,946 $\pm$ 240 K by comparing the observed spectra and the Kurucz models. We detected the absorption lines of the secondary star, which had never been detected in previous studies, and obtained the radial velocities for both components. With the published $BVRI$ light curves, we determined the absolute parameters for the binary via Wilson-Devinney modeling. The masses and radii are $M_{1} = 2.49 \pm 0.06$ $M_\odot$, $M_{2} = 0.38 \pm 0.01$ $M_\odot$, $R_{1} = 2.27 \pm 0.02$ $R_\odot$, and $R_{2} = 2.43 \pm 0.02$ $R_\odot$, respectively. The primary star is about $45 \%$ more massive and $60 \%$ larger than the zero-age main sequence (ZAMS) stars with the same effective temperature. It is probably because XX Cep has experienced a very different evolutionary process due to mass transfer, contrasting with the normal main sequence stars. The primary star is located inside the theoretical instability strip of $δ$ Sct-type stars on HR diagram. We demonstrated that XX Cep is an oEA star, consisting of a $δ$ Sct-type pulsating primary component and an evolved secondary companion.

preprint2015arXiv

KIC 4739791: A New R CMa-type Eclipsing Binary with a Pulsating Component

The {\it Kepler} light curve of KIC 4739791 exhibits partial eclipses, inverse O'Connell effect, and multiperiodic pulsations. Including a starspot on either of the binary components, the light-curve synthesis indicates that KIC 4739791 is in detached or semi-detached configurations with both a short orbital period and a low mass ratio. Multiple frequency analyses were performed in the light residuals after subtracting the binarity effects from the original {\it Kepler} data. We detected 14 frequencies: six in the low-frequency region (0.1$-$2.3 d$^{-1}$) and eight in the high-frequency region (18.2$-$22.0 d$^{-1}$). Among these, six high frequencies with amplitudes of 0.62$-$1.97 mmag were almost constant over time for 200 d. Their pulsation periods and pulsation constants are in the ranges of 0.048$-$0.054 d and 0.025$-$0.031 d, respectively. In contrast, the other frequencies may arise from the alias effects caused by the orbital frequency or combination frequencies. We propose that KIC 4739791 is a short-period R CMa binary with the lowest mass ratio in the known classical Algols and that its primary component is a $δ$ Sct pulsating star. Only four R CMa stars have been identified, three of which exhibit $δ$ Sct-type oscillations. These findings make KIC 4739791 an attractive target for studies of stellar interior structure and evolution.

preprint2015arXiv

State and group dynamics of world stock market by principal component analysis

We study the dynamic interactions and structural changes in global financial indices in the years 1998-2012. We apply a principal component analysis (PCA) to cross-correlation coefficients of the stock indices. We calculate the correlations between principal components (PCs) and each asset, known as PC coefficients. A change in market state is identified as a change in the first PC coefficients. Some indices do not show significant change of PCs in market state during crises. The indices exposed to the invested capitals in the stock markets are at the minimum level of risk. Using the first two PC coefficients, we identify indices that are similar and more strongly correlated than the others. We observe that the European indices form a robust group over the observation period. The dynamics of the individual indices within the group increase in similarity with time, and the dynamics of indices are more similar during the crises. Furthermore, the group formation of indices changes position in two-dimensional spaces due to crises. Finally, after a financial crisis, the difference of PCs between the European and American indices narrows.

preprint2015arXiv

The Physical Nature and Orbital Behavior of the Eclipsing System DK Cygni

New CCD photometry is presented for the hot overcontact binary DK Cyg, together with reasonable explanations for the light and period variations. Historical light and velocity curves from 1962 to 2012 were simultaneously analyzed with the Wilson-Devinney (W-D) synthesis code. The brightness disturbances were satisfactorily modeled by applying a magnetic cool spot on the primary star. Based on 261 times of minimum light including 116 new timings and spanning more than 87 yrs, a period study reveals that the orbital period has varied due to a periodic oscillation superposed on an upward parabola. The period and semi-amplitude of the modulation are about 78.1 yrs and 0.0037 d, respectively. This detail is interpreted as a light-travel-time effect due to a circumbinary companion with a minimum mass of $M_3$=0.065 $M_\odot$, within the theoretical limit of $\sim$0.07 M$_\odot$ for a brown dwarf star. The observed period increase at a fractional rate of $+$2.74 $\times $10$^{-10}$ is in excellent agreement with that calculated from our W-D synthesis. Mass transfer from the secondary to the primary component is mainly responsible for the secular period change. We examined the evolutionary status of the DK Cyg system from the absolute dimensions.

preprint2014arXiv

1SWASP J093010.78+533859.5: A Possible Hierarchical Quintuple System

We present the observational results of this kind of rare object 1SWASP J093010.78+533859.5, for which the doubly eclips- ing feature had been detected previously from the SuperWASP photometric archive. Individual PSF photometry for two objects with a separation of about 1.9 arcsec was performed for the first time in this study. Our time-series photometric data confirms the finding of Lohr et al. (2013) that the bright object A is an Algol-type detached eclipsing binary and the fainter B is a W UMa-type contact eclipsing. Using the high- resolution optical spectra, we obtained well-defined radial velocity variations of system A. Furthermore, stationary spectral lines were detected that must have originated from a further, previously unrecognized stellar component. It was confirmed by the third object contribution from the light curve analysis. No spectral feature of the system B was detected, probably due to motion blur by long exposure time. We obtained the binary parameters and the absolute dimensions of the systems A and B from light curve synthesis with and without radial velocities, respectively. The primary and secondary components of system A have a spectral type of K1 and K5 main sequences, respec- tively. Two components of system B have nearly the same type of K3 main sequence. Light variations for both binaries are satisfactorily modeled by using two-spot models with one starspot on each component. We estimated the distances to systems A and B individually. Two systems may have similar distances of about 70 pc and seem to be gravitationally bound with a separation of about 130 AU. In conclusion, we suggest that 1SWASP J093010.78+533859.5 is a quintuple stellar system with a hierarchical structure of a triple system A(ab)c and a binary system B(ab).

preprint2014arXiv

Correlation and Network Topologies in Global and Local Stock Indices

This study examined how the correlation and network structure of 30 global indices and 145 local Korean indices belonging to the KOSPI 200 have changed during the 13-year period, 2000-2012. The correlations among the indices were calculated. The results showed that although the average correlations of the global indices increased with time, the local indices showed a decreasing trend except for drastic changes during crises. The average correlation of the local indices exceeded the global indices during the crises from 2000-2002, implying a strong correlation structure among the local indices during this period due to the detrimental effect of the dot-com bubble. The threshold networks (TN) were constructed in the observation time window by assigning a threshold value and determining the network topologies. A significant change in the network topologies was observed due to the financial crises in both markets. The Jaccard similarities were also determined using the common links of TNs. The TNs of the financial network were not consistent with the evolution of the time, and the successive TNs of the global indices were more similar than those of the successive local indices. Finally, the Jaccard similarities identified the change in the market state due to a crisis in both markets.

preprint2014arXiv

On the RZ Draconis Sub-stellar Circumbinary Companions Stability Study of the Proposed Sub-stellar Circumbinary System

In this work we revisit the proposed multi-circumbinary system RZ Dra. We find the proposed system to be highly unstable. We attempt to find a best-fit light-travel time model rendering the orbits to follow stable orbits. We found a best-fit solution, but the resulting orbits also exhibit short-term orbital instabilities. We therefore conclude that the observed timing variation must be due to other astrophysical effects causing a change in the binary orbital period. Follow-up observations of this system is encouraged.

preprint2014arXiv

Physical Nature and Timing Variations of the Eclipsing System V407 Pegasi

New multiband CCD photometry is presented for V407 Peg; the $R_{\rm C}$ light curves are the first ever compiled. Our light curves, displaying a flat bottom at secondary minimum and an O'Connell effect, were simultaneously analyzed with the radial-velocity (RV) curves given by Rucinski et al. (2008). The light changes of the system are best modeled by using both a hot spot on the secondary star and a third light. The model represents historical light curves also. All available minimum epochs, including our six timing measurements, have been examined and indicate that the eclipse timing variation is mainly caused by light asymmetries due to the spot activity detected in the light-curve synthesis. The hot spot may be produced as a result of the impact of the gas stream from the primary star. Our light and velocity solutions indicate that V407 Peg is a totally-eclipsing A-type overcontact binary with values of $q$=0.251, $i$=87$^\circ.6$, $ΔT$=496 K, $f$=61 $\%$, and $l_3$=11$\sim$16 $\%$. Individual masses and radii of both components are determined to be $M_1$=1.72 M$_\odot$, $M_2$=0.43 M$_\odot$, $R_1$=2.15 R$_\odot$, and $R_2$=1.21 R$_\odot$. These results are very different from previous ones, which is probably caused by the light curves with distorted and inclined eclipses used in those other analyses. The facts that there are no objects optically related with the system and that the seasonal RVs show a large discrepancy in systemic velocity indicate that the third light source most likely arises from a tertiary component orbiting the eclipsing pair.

preprint2014arXiv

Self-organized criticality of a simplified integrate-and-fire neural model on random and small-world network

We consider the criticality for firing structures of a simplified integrate-and-fire neural model on the regular network, small-world network, and random networks. We simplify an integrate-and-fire model suggested by Levina, Herrmann and Geisel (LHG). In our model we set up the synaptic strength as a constant value. We observed the power law behaviors of the probability distribution of the avalanche size and the life time of the avalanche. The critical exponents in the small-world network and the random network were the same as those in the fully connected network. However, in the regular one-dimensional ring, the model does not show the critical behaviors. In the simplified LHG model, the short-cuts are crucial role in the self-organized criticality. The simplified LHG model in three types of networks such as the fully connected network, the small-world network, and random network belong to the same universality class.

preprint2014arXiv

The Eclipsing System V404 Lyr: Light-Travel Times and Gamma Doradus Pulsations

We present the physical properties of V404 Lyr exhibiting eclipse timing variations and multiperiodic pulsations from all historical data including the Kepler observations. Detailed analyses of 2,922 minimum epochs showed that the orbital period has varied through a combination with an upward-opening parabola and two sinusoidal variations, with periods of $P_3$=649 d and $P_4$=2,154 d and semi-amplitudes of $K_3$=193 s and $K_4$=49 s, respectively. The secular period increase could be interpreted as a combination of the secondary to primary mass transfer and angular momentum loss. The most reasonable explanation for both sinusoids is a pair of light-travel-time effects due to two circumbinary objects with projected masses of $M_3$=0.47 M$_\odot$ and $M_4$=0.047 M$_\odot$. The third-body parameters are consistent with those calculated using the W-D binary code. For the orbital inclinations higher than 43$^\circ$, the fourth component has a mass within the hydrogen-burning limit of $\sim$0.07 M$_\odot$, which implies that it is a brown dwarf. A satisfactory model for the Kepler light curves was obtained through applying a cool spot to the secondary component. The results demonstrate that the close eclipsing pair is in a semi-detached, but near-contact, configuration; the primary fills approximately 93% of its limiting lobe and is larger than the lobe-filling secondary. Multiple frequency analyses were applied to the light residuals after subtracting the synthetic eclipsing curve from the Kepler data. This revealed that the primary component of V404 Lyr is a $γ$ Dor type pulsating star, exhibiting seven pulsation frequencies in the range of 1.85$-$2.11 d$^{-1}$ with amplitudes of 1.38$-$5.72 mmag and pulsation constants of 0.24$-$0.27 d. The seven frequencies were clearly identified as high-order low-degree gravity-mode oscillations which might be excited through tidal interaction.

preprint2014arXiv

The Kepler Eclipsing System KIC 5621294 and its Substellar Companion

We present the physical properties of KIC 5621294 showing light and timing variations from the ${\it Kepler}$ photometry. Its light curve displays partial eclipses and O'Connell effect with Max II fainter than Max I, which was fitted quite well by applying third-body and spot effects to the system. The results indicate that the eclipsing pair is a classical Algol-type system with parameters of $q$=0.22, $i$=76$^\circ$.8, and $Δ$($T_{1}$--$T_{2}$)=4,235 K, in which the detached primary component fills about 77\% of its limiting lobe. Striking discrepancies exist between the primary and secondary eclipse times obtained with the method of Kwee \& van Woerden. These are mainly caused by surface inhomogeneities due to spot activity detected in our light-curve synthesis. The 1,253 light-curve timings from the Wilson-Devinney code were used for a period study. It was found that the orbital period of KIC 5621294 has varied due to a periodic variation overlaid on a downward parabola. The sinusoidal variation with a period of 961 d and a semi-amplitude of 22.5 s most likely arise from a light-time effect due to a third component with a mass of $M_3 \sin i_3$=46.9 M$\rm_{Jup}$, which is in good agreement with that calculated from the light curve itself. If its orbital inclination is larger than about 40$^\circ$, the mass of the circumbinary object would possibly match a brown dwarf. The parabolic variation could not be fully explained by either a mass transfer between the binary components or an angular momentum via magnetic braking. It is possible that the parabola may only be the observed part of a period modulation caused by the presence of another companion in a wider orbit.

preprint2014arXiv

The Pulsating sdB+M Eclipsing System NY Virginis and its Circumbinary Planets

We searched for circumbinary planets orbiting NY Vir in historical eclipse times including our long-term CCD data. Sixty-eight times of minimum light with accuracies better than 10 s were used for the ephemeris computations. The best fit to those timings indicated that the orbital period of NY Vir has varied due to a combination of two sinusoids with periods of $P_3$=8.2 yr and $P_4$=27.0 yr and semi-amplitudes of $K_3$=6.9 s and $K_4$=27.3 s, respectively. The periodic variations most likely arise from a pair of light-time effects due to the presence of third and fourth bodies that are gravitationally bound to the eclipsing pair. We have derived the orbital parameters and the minimum masses, $M_3 \sin i_3$ = 2.8 M$\rm_{Jup}$ and $M_4 \sin i_4$ = 4.5 M$\rm_{Jup}$, of both objects. A dynamical analysis suggests that the outer companion is less likely to orbit the binary on a circular orbit. Instead we show that future timing data might push its eccentricity to moderate values for which the system exhibits long-term stability. The results demonstrate that NY Vir is probably a star-planet system, which consists of a very close binary star and two giant planets. The period ratio $P_{3}/P_{4}$ suggests that a long-term gravitational interaction between them would result in capture into a nearly 3:10 mean motion resonance. When the presence of the circumbinary planets is verified and understood more comprehensively, the formation and evolution of this planetary system should be advanced greatly.

preprint2013arXiv

Network Topologies of Financial Market During the Global Financial Crisis

We consider the effects of the global financial crisis through a local Korean financial market around the 2008 crisis. We analyze 185 individual stock prices belonging to the KOSPI (Korea Composite Stock Price Index), cosidering three time periods: the time before, during, and after the crisis. The complex networks generate from the fully connected correlation network by using the cross-correlation coefficients among the stock price time series of the companies. We generate the threshold networks (TN), the minimal spanning tees (MST), and the hierarchical network (HN) from the fully connected cross-correlation networks. By assigning a threshold value of the cross-correlation coefficient, we obtain the threshold networks. We observe the power law of the degree distribution in the limited range of the threshold. The degree distribution of the largest cluster in the threshold networks during the crisis is fatter than other periods. The clustering coefficient of the threshold networks follows the power law in the scaling range. We also generate the minimal spanning trees from the fully connected correlation networks. The MST during the crisis period shrinks in comparison to the periods before and after the crisis. The cophenetic correlation coefficient increases during the crisis, indicating that the hierarchical structure increases during this period. When the crisis hit the market, the companies behave synchronously and their correlations become stronger than the normal period.

preprint2013arXiv

Random Matrix Theory and Cross-correlations in Global Financial Indices and Local Stock Market Indices

We analyzed cross-correlations between price fluctuations of global financial indices (20 daily stock indices over the world) and local indices (daily indices of 200 companies in the Korean stock market) by using random matrix theory (RMT). We compared eigenvalues and components of the largest and the second largest eigenvectors of the cross-correlation matrix before, during, and after the global financial the crisis in the year 2008. We find that the majority of its eigenvalues fall within the RMT bounds [λ_, λ+], where λ_- and λ_+ are the lower and the upper bounds of the eigenvalues of random correlation matrices. The components of the eigenvectors for the largest positive eigenvalues indicate the identical financial market mode dominating the global and local indices. On the other hand, the components of the eigenvector corresponding to the second largest eigenvalue are positive and negative values alternatively. The components before the crisis change sign during the crisis, and those during the crisis change sign after the crisis. The largest inverse participation ratio (IPR) corresponding to the smallest eigenvector is higher after the crisis than during any other periods in the global and local indices. During the global financial the crisis, the correlations among the global indices and among the local stock indices are perturbed significantly. However, the correlations between indices quickly recover the trends before the crisis.

preprint2013arXiv

Revisiting the proposed circumbinary multi-planet system NSVS14256825

In this work we carry out an analysis of the observed times of primary and secondary eclipses of the post-common envelope binary NSVS14256825. Recently, \cite{Almeida2013} proposed that two circumbinary companions orbit this short-period eclipsing binary, in order to explain observed variations in the timing of mutual eclipses between the two binary components. Using a standard weighted least-squares minimisation technique, we have extensively explored the topology of $χ^2$ parameter space of a single planet model. We find the data set to be insufficient to reliably constrain a one-companion model. Various models, each with similar statistical significance, result in substantially different orbital architectures for the additional companion. No evidence is seen for a second companion of planetary nature. We suspect insufficient coverage (baseline) of timing data causing the best-fit parameters to be unconstrained.

preprint2013arXiv

The Eclipsing System EP Andromedae and its Circumbinary Companions

We present new long-term CCD photometry for EP And acquired during the period 2007 to 2012. The light curves display total eclipses at primary minima and season-to-season light variability. Our synthesis for all available light curves indicates that the eclipsing pair is a W-type overcontact binary with parameters of q=2.578, i=83.3 deg, $ΔT$=27 K, f=28 %, and l3=2~3 %. The asymmetric light curves in 2007 were satisfactorily modeled by a cool spot on either of the eclipsing components from a magnetic dynamo. Including our 95 timing measurements, a total of 414 times of minimum light spanning about 82 yr were used for a period study. A detailed analysis of the eclipse timing diagram revealed that the orbital period of EP And has varied as a combination of an upward-opening parabola and two periodic variations, with cycle lengths of P3=44.6 yr and P4=1.834 yr and semi-amplitudes of K3=0.0100 d and K4=0.0039 d, respectively. The observed period increase at a fractional rate of +1.39x10^{-10} is in excellent agreement with that calculated from the W-D code and can be plausibly explained by some combination of mass transfer from the primary to the secondary star and angular momentum loss due to magnetic braking. The most reasonable explanation for both cycles is a pair of light-travel-time effects driven by the possible existence of a third and fourth component with projected masses of M3=0.25 Msun and M4=0.90 Msun. The more massive companion could be revealed using high-resolution spectroscopic data extending over the course of a few years and could also be a binary itself. It is possible that the circumbinary objects may have played an important role in the formation and evolution of the eclipsing pair, which would cause it to have a short initial orbital period and thus evolve into an overcontact configuration by angular momentum loss.

preprint2012arXiv

Interspecific competition underlying mutualistic networks

The architecture of bipartite networks linking two classes of constituents is affected by the interactions within each class. For the bipartite networks representing the mutualistic relationship between pollinating animals and plants, it has been known that their degree distributions are broad but often deviate from power-law form, more significantly for plants than animals. Here we consider a model for the evolution of the mutualistic networks and find that their topology is strongly dependent on the asymmetry and non-linearity of the preferential selection of mutualistic partners. Real-world mutualistic networks analyzed in the framework of the model show that a new animal species determines its partners not only by their attractiveness but also as a result of the competition with pre-existing animals, which leads to the stretched-exponential degree distributions of plant species.

preprint2012arXiv

Physical Properties of the Low-Mass Eclipsing Binary NSVS 02502726

NSVS 02502726 has been known as a double-lined, detached eclipsing binary that consists of two low-mass stars. We obtained $BVRI$ photometric follow-up observations in 2009 and 2011 to measure improved physical properties of the binary star. Each set of light curves, including the 2008 data given by \cCakirli et al., was simultaneously analyzed with the previously published radial-velocity curves using the Wilson-Devinney binary code. The conspicuous seasonal light variations of the system are satisfactorily modelled by a two-spot model with one starspot on each component and by changes of the spot parameters with time. Based on 23 eclipse timings calculated from the synthetic model and one ephemeris epoch, an orbital period study of NSVS 02502726 reveals that the period has experienced a continuous decrease of $-5.9\times10^{-7}$ d yr$^{-1}$ or a sinusoidal variation with a period and semi-amplitude of 2.51 yrs and 0.0011 d, respectively. The timing variations could be interpreted as either the light-travel-time effect due to the presence of an unseen third body, or as the combination of this effect and angular momentum loss via magnetic stellar wind braking. Individual masses and radii of both components are determined to be $M_1$=0.689$\pm$0.016 M$_\odot$, $M_2$=0.341$\pm$0.009 M$_\odot$, $R_1$=0.707$\pm$0.007 R$_\odot$, and $R_2$=0.657$\pm$0.008 R$_\odot$. The results are very different from those of \cCakirli et al. with the primary's radius (0.674$\pm$0.006 R$_\odot$) smaller the secondary's (0.763$\pm$0.007 R$_\odot$). We compared the physical parameters presented in this paper with current low-mass stellar models and found that the measured values of the primary star are best fitted to a 79-Myr isochrone. The primary is in good agreement with the empirical mass-radius relation from low-mass binaries, but the secondary is oversized by about 85%.

preprint2012arXiv

Scaling of nestedness in complex networks

Nestedness characterizes the linkage pattern of networked systems, indicating the likelihood that a node is linked to the nodes linked to the nodes with larger degrees than it. Networks of mutualistic relationship between distinct groups of species in ecological communities exhibit such nestedness, which is known to support the network robustness. Despite such importance, quantitative characteristics of nestedness is little understood. Here we take graph-theoretic approach to derive the scaling properties of nestedness in various model networks. Our results show how the heterogeneous connectivity patterns enhance nestedness. Also we find that the nestedness of bipartite networks depend sensitively on the fraction of different types of nodes, causing nestedness to scale differently for nodes of different types.

preprint2012arXiv

The Light and Period Variations of the Eclipsing Binary BX Draconis

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.

preprint2012arXiv

The sub-Saturn Mass Transiting Planet HAT-P-12b

We present new photometric data of the transiting planet HAT-P-12b observed in 2011. Our three transit curves are modelled using the JKTEBOP code and adopting the quadratic limb-darkening law. Including our measurements, 18 transit times spanning about 4.2 yr were used to determine the improved ephemeris with a transit epoch of 2,454,187.85560$\pm$0.00011 BJD and an orbital period of 3.21305961$\pm$0.00000035 d. The physical properties of the star-planet system are computed using empirical calibrations from eclipsing binary stars and stellar evolutionary models, combined with both our transit parameters and previously-known spectroscopic results. We found that the absolute dimensions of the host star are $M_{\rm A}$=0.73$\pm$0.02 M$_\odot$, $R_{\rm A}$=0.70$\pm$0.01 R$_\odot$, log $g_{\rm A}$=4.61$\pm$0.02, $ρ_{\rm A}$=2.10$\pm$0.09 $ρ_{\odot}$, and $L_{\rm A}$=0.21$\pm$0.01 L$_\odot$. The planetary companion has $M_{\rm b}$=0.21$\pm$0.01 M$_{\rm Jup}$, $R_{\rm b}$=0.94$\pm$0.01 R$_{\rm Jup}$, log $g_{\rm b}$=2.77$\pm$0.02, $ρ_{\rm b}$=0.24$\pm$0.01 $ρ_{\rm Jup}$, and $T_{\rm eq}$=960$\pm$14 K. Our results agree well with standard models of irradiated gas giants with a core mass of 11.3 M$_{\oplus}$.

preprint2012arXiv

The Triply Eclipsing Hierarchical Triple Star KIC002856960

In a recent study, Armstrong et al. presented an eclipsing binary star of about 6.2 h period with transit-like tertiary signals occurring every 204.2 d in the {\it Kepler} public data of KIC002856960 and proposed three possible hierarchical structures: (AB)b, (AB)C, and A(BC). We analyzed the {\it Kepler} light curve by including a third light source and one starspot on each binary component. The results represent that the close eclipsing pair is in a low-mass eccentric-orbit, detached configuration. Based on 123 eclipse timings calculated from the Wilson-Devinney binary model, a period study of the close binary reveals that the orbital period has experienced a sinusoidal variation with a period and a semi-amplitude of 205$\pm$2 d and 0.0021$\pm$0.0002 d, respectively. The period variation would be produced by the light-travel-time effect due to a gravitationally-bound third body with a minimum mass of $M_3 \sin i_3$=0.76 M$_\odot$ in an eccentric orbit of $e_3$=0.61. This is consistent with the presence of third light found in our light-curve solution and the tertiary signal of 204.2 d period most likely arises from the K-type star crossed by the close eclipsing binary. Then, KIC002856960 is a triply eclipsing hierarchical system, A(BC), consisting of a close binary with two M-type dwarfs and a more massive K-type component. The presence of the third star may have played an important role in the formation and evolution of the close pair, which may ultimately evolve into a contact system by angular momentum loss.

preprint2012arXiv

Trends of Papers Published from 2006 to 2010 in Journals Nature and Science

We present an analysis of the papers published in the journals Nature and Science in the years from 2006 to 2010. During this period, a total of 7788 papers were published in the two journals. This includes 544 astronomy papers that comprise 7.0% of the papers in `all' research fields and 18.9% of those in the fields of `physical sciences'. The sub-fields of research of the astronomy papers are distributed, in descending order of number of papers, in Solar System, stellar astronomy, galaxies and the universe, the Milky Way Galaxy, and exoplanets. The observational facilities used for the studies are mainly ground-based telescopes (31.1%), spacecrafts (27.0%), and space telescopes (22.8%), while 16.0% of papers did not use any noticeable facilities and 1.7% used other facilities. Korean scientists have published 86 papers (33 in Nature and 53 in Science), which is 1.10% of all the papers (N=7788) in the two journals. The share of papers by Korean astronomers among the scientific papers by Koreans is 8.14%, slightly higher than the contribution of astronomy papers (7.0%) in both journals.

preprint2011arXiv

New Light-Travel Time Models and Orbital Stability Study of the Proposed Planetary System HU Aquarii

In this work we propose a new orbital architecture for the two proposed circumbinary planets around the polar eclipsing binary HU Aquarii. We base the new two-planet, light-travel time model on the result of a Monte Carlo simulation driving a least-squares Levenberg-Marquardt minimisation algorithm on the observed eclipse egress times. Our best-fitting model with $χ_{r}^2=1.43$ resulted in high final eccentricities for the two companions leading to an unstable orbital configuration. From a large ensemble of initial guesses we examined the distribution of final eccentricities and semi-major axes for different $χ_{r}^2$ parameter intervals and encountered qualitatively a second population of best-fitting parameters. The main characteristic of this population is described by low-eccentric orbits favouring long-term orbital stability of the system. We present our best-fitting model candidate for the proposed two-planet system and demonstrate orbital stability over one million years using numerical integrations.

preprint2011arXiv

The Algol System SZ Herculis: Physical Nature and Orbital Behavior

Multiband CCD photometric observations of SZ Her were obtained between 2008 February and May. The light curve was completely covered and indicated a significant temperature difference between both components. The light-curve synthesis presented in this paper indicates that the eclipsing binary is a classical Algol-type system with parameters of $q$=0.472, $i$=87$^\circ$.57, and $Δ$($T_{1}$--$T_{2}$)=2,381 K; the primary component fills approximately 77% of its limiting lobe and is slightly larger than the lobe-filling secondary. More than 1,100 times of minimum light spanning more than one century were used to study an orbital behavior of the binary system. It was found that the orbital period of SZ Her has varied due to a combination of two periodic variations with cycle lengths of $P_3$=85.8 yr and $P_4$=42.5 yr and semi-amplitudes of $K_3$=0.013 d and $K_4$=0.007 d, respectively. The most reasonable explanation for them is a pair of light-time effects (LITEs) driven by the possible existence of two M-type companions with minimum masses of $M_3$=0.22 M$_\odot$ and $M_4$=0.19 M$_\odot$, that are located close to the 2:1 mean motion resonance. If two additional bodies exist, then the overall dynamics of the multiple system may provide a significant clue to the formation and evolution of the eclipsing pair.

preprint2011arXiv

WZ Cygni: a Marginal Contact Binary in a Triple System?

We present new multiband CCD photometry for WZ Cyg made on 22 nights in two observing seasons of 2007 and 2008. Our light-curve synthesis indicates that the system is in poor thermal contact with a fill-out factor of 4.8% and a temperature difference of 1447 K. Including our 40 timing measurements, a total of 371 times of minimum light spanning more than 112 yr were used for a period study. Detailed analysis of the $O$--$C$ diagram showed that the orbital period has varied by a combination with an upward parabola and a sinusoid. The upward parabola means the continuous period increase and indicates that some stellar masses are thermally transferred from the less to the more massive primary star at a rate of about 5.80$\times10^{-8}$ M$_\odot$ yr$^{-1}$. The sinusoidal variation with a period of 47.9 yr and a semi-amplitude of 0.008 d can be interpreted most likely as the light-travel-time effect due to the existence of a low-mass M-type tertiary companion with a projected mass of $M_3 \sin i_3$=0.26 M$_\odot$. We examined the evolutionary status of WZ Cyg from the absolute dimensions of the eclipsing pair. It belongs to the marginal contact binary systems before broken-contact phase, consisting of a massive primary star with spectral type of F4 and a secondary with the type K1.

preprint2010arXiv

Physical Properties of the Transiting Planetary System TrES-3

We present four new transits of the planetary system TrES-3 observed between 2009 May and 2010 June. Among these, the third transit by itself indicates possible evidence for brightness disturbance, which might be the result of the planet blocking a cool starspot on the stellar surface. A total of 109 transit times, including our measurements, were used to determine the improved ephemeris with a transit epoch of 2454185.910944$\pm$0.000072 HJED and an orbital period of 1.30618700$\pm$0.00000015 d. We analyzed the transit light curves using the JKTEBOP code and adopting the quadratic limb-darkening law. In order to derive the physical properties of the TrES-3 system, the transit parameters are combined with the empirical relations from eclipsing binary stars and stellar evolutionary models. The stellar mass and radius obtained from a calibration using $T_A$, log $ρ_{\rm A}$ and [Fe/H] are consistent with those from the isochrone analysis. We found that the exoplanet TrES-3b has a mass of 1.93$\pm$0.07 M$_{\rm Jup}$, a radius of 1.30$\pm$0.04 R$_{\rm Jup}$, a surface gravity of log $g_{\rm b}$=3.45$\pm$0.02, a density of 0.82$\pm$0.06 $ρ_{\rm Jup}$, and an equilibrium temperature of 1641$\pm$23 K. The results are in good agreement with theoretical models for gas giant planets.

preprint2010arXiv

The First Comprehensive Photometric Study of the Algol-type System CL Aurigae

We present the first extensive photometric results of CL Aur from our BVRI CCD photometry made on 22 nights from 2003 November through 2005 February. Fifteen new timings of minimum light were obtained. During the past 104 years, the orbital period has varied due to a periodic oscillation superposed on a continuous period increase. The period and semi-amplitude of the oscillation are about 21.6 yrs and 0.0133 d, respectively. This detail is interpreted as a light-travel-time effect due to a low-luminosity K-type star gravitationally bound to the CL Aur close system. Our photometric study indicates that CL Aur is a relatively short-period Algol-type binary with values of q=0.602 and i=88$^\circ$.2. Mass transfer from the secondary to the primary eclipsing component is at least partly responsible for the observed secular period change with a rate of dP/dt = +1.4$\times10^{-7}$ d yr$^{-1}$. A cool spot model has been calculated but we think that an alternative hot-spot model resulting from a gas stream impact on the hot star is more reasonable despite two difficulties with the explanation. Absolute dimensions of the eclipsing system are deduced and its present state is compared with tracks for single star and conservative close binary evolution. Finally, we examine the possible reconciliation of two different calculations of the luminosity of the hot spot and a re-interpretation of the secular term of the period variability.

preprint2010arXiv

The Light and Period Variations of the Eclipsing Binary AA Ursae Majoris

We present new multiband CCD photometry for AA UMa made on 8 nights between January and March 2009; the $R$ light curves are the first ever compiled. Historical light curves, as well as ours, display partial eclipses and inverse O'Connell effects with Max I fainter than Max II. Among possible spot models, a cool spot on either of the component stars and its variability with time permit good light-curve representations for the system. A total of 194 eclipse timings over 81 yrs, including our five timings, were used for ephemeris computations. We found that the orbital period of the system has varied due to a periodic oscillation overlaid on an upward parabolic variation. The continuous period increase at a fractional rate of $+$1.3$\times$10$^{-10}$ is consistent with that calculated from the W-D code and can be interpreted as a thermal mass transfer from the less to the more massive secondary star at a rate of 6.6$\times$10$^{-8}$ M$_\odot$ yr$^{-1}$. The periodic component is in satisfactory accord with a light-time effect due to an unseen companion with a period of 28.2 yrs, a semi-amplitude of 0.007 d, and a minimum mass of $M_3 \sin i_3$=0.25 $M_\odot$ but this period variation could also arise from magnetic activity.

preprint2009arXiv

Long-term Photometric Behavior of the Eclipsing Binary GW Cephei

New CCD photometry during 4 successive years from 2005 is presented for the eclipsing binary GW Cep, together with reasonable explanations for the light and period variations. All historical light curves, obtained over a 30-year interval, display striking light changes, and are best modeled by the simultaneous existence of a cool spot and a hot spot on the more massive cool component star. The facts that the system is magnetically active and that the hot spot has consistently existed on the inner hemisphere of the star indicate that the two spots are formed by (1) magnetic dynamo-related activity on the cool star and (2) mass transfer from the primary to the secondary component. Based on 38 light-curve timings from the Wilson-Devinney code and all other minimum epochs, a period study of GW Cep reveals that the orbital period has experienced a sinusoidal variation with a period and semi-amplitude of 32.6 yrs and 0.009 d, respectively. In principle, these may be produced either by a light-travel-time effect due to a third body or by an active magnetic cycle of at least one component star. Because we failed to find any connection between luminosity variability and the period change, that change most likely arises from the existence of an unseen third companion star with a minimum mass of 0.22 $M_\odot$ gravitationally bound to the eclipsing pair.

preprint2009arXiv

The Period Variation of and a Spot Model for the Eclipsing Binary AR Bootis

New CCD photometric observations of the eclipsing system AR Boo were obtained from February 2006 to April 2008. The star's photometric properties are derived from detailed studies of the period variability and of all available light curves. We find that over about 56 years the orbital period of the system has varied due to a combination of an upward parabola and a sinusoid rather than in a monotonic fashion. Mass transfer from the less massive primary to the more massive secondary component is likely responsible for at least a significant part of the secular period change. The cyclical variation with a period of 7.57 yrs and a semi-amplitude of 0.0015 d can be produced either by a light-travel-time effect due to an unseen companion with a scaled mass of $M_3 \sin i_3$=0.081 $M_\odot$ or by a magnetic period modulation in the secondary star. Historical light curves of AR Boo, as well as our own, display season-to-season light variability, which are best modeled by including both a cool spot and a hot one on the secondary star. We think that the spots express magnetic dynamo-related activity and offer limited support for preferring the magnetic interpretation of the 7.57-year cycle over the third-body understanding. Our solutions confirm that AR Boo belongs to the W-subtype contact binary class, consisting of a hotter, less massive primary star with a spectral type of G9 and a companion of spectral type K1.

preprint2009arXiv

The Solar-Type Contact Binary BX Pegasi Revisited

We present the results of new CCD photometry for the contact binary BX Peg, made during three successive months beginning on September 2008. As do historical light curves, our observations display an O'Connell effect and the November data by themselves indicate clear evidence for very short-time brightness disturbance. For these variations, model spots are applied separately to the two data set of Group I (Sep.--Oct.) and Group II (Nov.). The former is described by a single cool spot on the secondary photosphere and the latter by a two-spot model with a cool spot on the cool star and a hot one on either star. These are generalized manifestations of the magnetic activity of the binary system. Twenty light-curve timings calculated from Wilson-Devinney code were used for a period study, together with all other minimum epochs. The complex period changes of BX Peg can be sorted into a secular period decrease caused dominantly by angular momentum loss due to magnetic stellar wind braking, a light-travel-time (LTT) effect due to the orbit of a low-mass third companion, and a previously unknown short-term oscillation. This last period modulation could be produced either by a second LTT orbit with a period of about 16 yr due to the existence of a fourth body or by the effect of magnetic activity with a cycle length of about 12 yr.

preprint2004arXiv

Reversible random sequential adsorption on a one-dimensional lattice

We consider the reversible random sequential adsorption of line segments on a one-dimensional lattice. Line segments of length $l \geq 2$ adsorb on the lattice with a adsorption rate $K_a$, and leave with a desorption rate $K_d$. We calculate the coverage fraction, and steady-state jamming limits by a Monte Carlo method. We observe that coverage fraction and jamming limits do not follow mean-field results at the large $K=K_a/K_d >>1$. Jamming limits decrease when the length of the line segment $l$ increases. However, jamming limits increase monotonically when the parameter $K$ increases. The distribution of two consecutive empty sites is not equivalent to the square of the distribution of isolated empty sites.

preprint1997arXiv

Effects of Impurities in Random Sequential Adsorption on a One-Dimensional Substrate

We have solved the kinetics of random sequential adsorption of linear $k$-mers on a one-dimensional disordered substrate for the random sequential adsorption initial condition and for the random initial condition. The jamming limits $θ(\infty, k', k)$ at fixed length of linear $k$-mers have a minimum point at a particular density of the linear $k'$-mers impurity for both cases. The coverage of the surface and the jamming limits are compared to the results for Monte Carlo simulation. The Monte Carlo results for the jamming limits are in good agreement with the analytical results. The continuum limits are derived from the analytical results on lattice substrates.

preprint1997arXiv

Irreversible Deposition of Line Segment Mixtures on a Square Lattice: Monte Carlo Study

We have studied kinetics of random sequential adsorption of mixtures on a square lattice using Monte Carlo method. Mixtures of linear short segments and long segments were deposited with the probability $p$ and $1-p$, respectively. For fixed lengths of each segment in the mixture, the jamming limits decrease when $p$ increases. The jamming limits of mixtures always are greater than those of the pure short- or long-segment deposition. For fixed $p$ and fixed length of the short segments, the jamming limits have a maximum when the length of the long segment increases. We conjectured a kinetic equation for the jamming coverage based on the data fitting.

preprint1997arXiv

Three-Species Diffusion-Limited Reaction with Continuous Density-Decay Exponents

We introduce a model of three-species two-particle diffusion-limited reactions A+B -> A or B, B+C -> B or C, and C+A -> C or A, with three persistence parameters (survival probabilities in reaction) of the hopping particle. We consider isotropic and anisotropic diffusion (hopping with a drift) in 1d. We find that the particle density decays as a power-law for certain choices of the persistence parameter values. In the anisotropic case, on one symmetric line in the parameter space, the decay exponent is monotonically varying between the values close to 1/3 and 1/2. On another, less symmetric line, the exponent is constant. For most parameter values, the density does not follow a power-law. We also calculated various characteristic exponents for the distance of nearest particles and domain structure. Our results support the recently proposed possibility that 1d diffusion-limited reactions with a drift do not fall within a limited number of distinct universality classes.