Massive exoplanets on extremely tight orbits, such as WASP-12 b, induce equilibrium tides in their host stars. Following the orbital motion of the planet, the tidal fluid flow in the star can be detected with the radial velocity method. Its signature manifests as the second harmonics of the orbital frequency that mimics a non-zero orbital eccentricity. Using the new radial velocity measurements acquired with the HARPS-N spectrograph at the Telescopio Nazionale Galileo and combining them with the literature data, we show that the apparent eccentricity of WASP-12 b's orbit is non-zero at a 5.8 sigma level, and the longitude of periastron of this apparently eccentric orbit is close to 270 degrees. This orbital configuration is compatible with a model composed of a circular orbit and a signature of tides raised in the host star. The radial velocity amplitude of those tides was found to be consistent with the equilibrium tide approximation. The tidal deformation is predicted to produce a flux modulation with an amplitude of 80 ppm which could be detected using space-born facilities.

We present the latest results of our spectroscopic observations and refined modelling of a sample of detached eclipsing binaries (DEBs), selected from the $Kepler$} Eclipsing Binary Catalog, that are also double-lined spectroscopic binaries (SB2). New high resolution spectra obtained with the HIDES spectrograph, attached to the 1.88-m telescope of the Okayama Astrophysical Observatory supplemented the previous observations, allowing to improve physical parameters of some systems, and characterize three previously omitted. All the obtained radial velocities (RVs) were combined with $Kepler$ photometry, in order to derive a full set of orbital and physical parameters. Ten out of eleven SB2s have their component spectra tomographically disentangled, and spectral analysis was performed with iSpec, in order to derive effective temperatures of components and metallicities of the systems. By comparing our results with theoretical models, we assess the age and evolutionary status of the studied objects. We find a good match to all but one systems. We have derived distances from thus determined parameters, and compare them with those from the $Gaia$ Data Release 2. For systems previously studied by other authors, our new results lead to better consistency between observations and models.

We present results of the combined photometric and spectroscopic analysis of three detached eclipsing binaries, which secondary components are not visible or very hard to identify in the optical spectra - ASAS J052743-0359.7, ASAS J065134-2211.5, and ASAS J073507-0905.7. The first one is a known visual binary ADS 4022, and we found that it is a quadruple system, composed of two spectroscopic binaries, one of which shows eclipses. None of the systems was previously recognized as a spectroscopic binary. We collected a number of high-resolution optical and IR spectra to calculate the radial velocities (RVs) and later combined them with MITSuME and ASAS photometry. The IR spectra were crucial for secure identification of the cooler components' lines. RV measurements were done with the TODCOR technique, and RV curves modelled with our own procedure V2FIT. Light curve modelling was performed with JKTEBOP and PHOEBE codes. Temperatures and metallicities of two systems were estimated from spectra. For the ADS 4022 system we also used the archival WDS data and new SOAR observations in order to derive the orbit of the visual pair for the first time. Ages were estimated by comparing our results with PARSEC isochrones. The eclipsing pair A052743 A ($P=5.27$ d) is composed of a 1.03(6) M$_\odot$, 1.03(2) R$_\odot$ primary and a 0.60(2) M$_\odot$, 0.59(2) R$_\odot$ secondary. The components of the $P=21.57$ d non-eclipsing pair B likely have masses in between the two eclipsing components, and both pairs are on a $\sim$188 yr orbit around their common centre of mass. The system A065134 ($P=8.22$ d) consists of a 0.956(12) M$_\odot$, 0.997(4) R$_\odot$ primary and a 0.674(5) M$_\odot$, 0.690(7) R$_\odot$ secondary. Finally, A073507 ($P=1.45$ d), which consists of a 1.452(34) M$_\odot$, 1.635(12) R$_\odot$ primary and a 0.808(13) M$_\odot$, 0.819(11) R$_\odot$ secondary, is likely a PMS system.

We have obtained high-resolution spectra of four eclipsing binary systems (FM Leo, NN Del, V963 Cen and AI Phe) with the view to gaining insight into the relative orientation of their stellar spin axes and orbital axes. The so called Rossiter-McLaughlin (RM) effect, i.e. the fact that the broadening and the amount of blue- or redshift in the spectra during an eclipse depend on the tilt of the spin axis of the background star, has the potential of reconciling observations and theoretical models if such a tilt is found. We analyse the RM effect by disentangling the spectra, removing the front component and measuring the remaining, distorted lines with a broadening function (BF) obtained from single value decomposition (SVD), weighting by the intensity centre of the BF in the eclipse. All but one of our objects show no significant misalignment, suggesting that aligned systems are dominant. We provide stellar as well as orbital parameters for our systems. With five measured spin-orbit angles we significantly increase (from 9 to 14) the number of stars for which it has been measured. The spin-orbit angle $\beta$ calculated for AI Phe's secondary component shows a misalignment of 87$\pm$17 degrees. NN Del, with a large separation of components and a long dynamical timescale for circularisation and synchronisation, is an example of a close to primordial spin-orbit angle measurement.

The spectroscopic binary system HD 66051 (V414 Pup) consists of a highly peculiar CP3 (HgMn) star and an A-type component. It also shows out-of-eclipse variability that is due to chemical spots. This combination allows the derivation of tight constraints for the testing of time-dependent diffusion models. We analysed radial velocity and photometric data using two different methods to determine astrophysical parameters and the orbit of the system. Appropriate isochrones were used to derive the age of the system. The orbital solution and the estimates from the isochrones are in excellent agreement with the estimates from a prior spectroscopic study. The system is very close to the zero-age main sequence and younger than 120 Myr. HD 66051 is a most important spectroscopic binary system that can be used to test the predictions of the diffusion theory explaining the peculiar surface abundances of CP3 stars.

We present the design and commissioning of Project Solaris, a global network of autonomous observatories. Solaris is a Polish scientific undertaking aimed at the detection and characterization of circumbinary exoplanets and eclipsing binary stars. To accomplish this, a network of four fully autonomous observatories has been deployed in the Southern Hemisphere: Solaris-1 and Solaris-2 in the South African Astronomical Observatory in South Africa; Solaris-3 in Siding Spring Observatory in Australia; and Solaris-4 in Complejo Astronomico El Leoncito in Argentina. The four stations are nearly identical and are equipped with 0.5-m Ritchey-Cr\'etien (f/15) or Cassegrain (f/9, Solaris-3) optics and high-grade 2 K x 2 K CCD cameras with Johnson and Sloan filter sets. We present the design and implementation of low-level security; data logging and notification systems; weather monitoring components; all-sky vision system, surveillance system; and distributed temperature and humidity sensors. We describe dedicated grounding and lighting protection system design and robust fiber data transfer interfaces in electrically demanding conditions. We discuss the outcomes of our design, as well as the resulting software engineering requirements. We describe our system's engineering approach to achieve the required level of autonomy, the architecture of the custom high-level industry-grade software that has been designed and implemented specifically for the use of the network. We present the actual status of the project and first photometric results; these include data and models of already studied systems for benchmarking purposes (Wasp-4b, Wasp-64b, and Wasp-98b transits, PG 1663-018, an eclipsing binary with a pulsator) as well J024946-3825.6, an interesting low-mass binary system for which a complete model is provided for the first time.

We present the results of our analysis of KIC 4150611 (HD 181469) - an interesting, bright quintuple system that includes a hybrid $\delta$ Sct/$\gamma$ Dor pulsator. Four periods of eclipses - 94.2, 8.65, 1.52 and 1.43 d - have been observed by the Kepler satellite, and three point sources (A, B, and C) are seen in high angular resolution images. From spectroscopic observations made with the HIDES spectrograph attached to the 1.88-m telescope of the Okayama Astrophysical Observatory (OAO), for the first time we calculated radial velocities (RVs) of the component B - a pair of G-type stars - and combined them with Kepler photometry in order to obtain absolute physical parameters of this pair. We also managed to directly measure RVs of the pulsator, also for the first time. Additionally, we modelled the light curves of the 1.52 and 1.43-day pairs, and measured their eclipse timing variations (ETVs). We also performed relative astrometry and photometry of three sources seen on the images taken with the NIRC2 camera of the Keck II telescope. Finally, we compared our results with theoretical isochrones. The brightest component Aa is the hybrid pulsator, transited every 94.2 days by a pair of K/M-type stars (Ab1+Ab2), which themselves form a 1.52-day eclipsing binary. The components Ba and Bb are late G-type stars, forming another eclipsing pair with a 8.65 day period. Their masses and radii are $M_{Ba}=0.894\pm0.010$ M$_\odot$, $R_{Ba}=0.802\pm0.044$ R$_\odot$ for the primary, and $M_{Bb}=0.888\pm0.010$ M$_\odot$, $R_{Bb}=0.856\pm0.038$ R$_\odot$ for the secondary. The remaining period of 1.43 days is possibly related to a faint third star C, which itself is most likely a background object. The system's properties are well-represented by a 35 Myr isochrone. There are also hints of additional bodies in the system.

We present the results of our spectroscopic observations of eight detached eclipsing binaries (DEBs), selected from the $Kepler$ Eclipsing Binary Catalog. Radial velocities (RVs) were calculated from high resolution spectra obtained with the HIDES pectrograph, attached to the 1.88-m telescope of the Okayama Astrophysical Observatory, and were used to characterize the targets in combination with the $Kepler$ light curves. For each binary we obtained a full set of orbital and physical parameters, reaching precision below 3% in masses and radii for 5 pairs. By comparing our results with theoretical models, we assess the distance, age and evolutionary status of the researched objects. We also study eclipse timing variations of selected objects, and identify a new system with a $\gamma$Dor pulsator. Two systems are triples, and show lines coming from three components. In one case the motion of the outer star and the perturbation in the RVs of the inner binary are clearly visible and periodical, which allows us to directly calculate the mass of the third star, and inclination of the outer orbit. In the second case we only see a clear motion of the tertiary, and investigate two scenarios: that it is a linear trend coming from the orbital motion around the inner binary, and that it is caused by a planetary mass companion. When possible, we also compare our results with the literature, and conclude that only by combining photometry with RVs it is possible to obtain correct physical parameters of both components of a DEB.

We present results of our spectroscopic observations of nine detached eclipsing binaries (DEBs), selected from the $Kepler$ Eclipsing Binary Catalog, that only show one set of spectral lines. Radial velocities (RVs) were calculated from the high resolution spectra obtained with the HIDES instrument, attached to the 1.88-m telescope at the Okayama Astrophysical Observatory, and from the public APOGEE archive. In our sample we found five single-lined binaries, with one component dominating the spectrum. The orbital and light curve solutions were found for four of them, and compared with isochrones, in order to estimate absolute physical parameters and evolutionary status of the components. For the fifth case we only update the orbital parameters, and estimate the properties of the unseen star. Two other systems show orbital motion with a period known from the eclipse timing variations (ETVs). For these we obtained parameters of outer orbits, by translating the ETVs to RVs of the centre of mass of the eclipsing binary, and combining with the RVs of the outer star. Of the two remaining ones, one is most likely a blend of a faint background DEB with a bright foreground star, which lines we see in the spectra, and the last case is possibly a quadruple bearing a sub-stellar mass object. Where possible, we compare our results with literature, especially with results from asteroseismology. We also report possible detections of solar-like oscillations in our RVs.

We present spectroscopic and photometric solutions for three spotted systems with red giant components. Absolute physical and orbital parameters for these double-lined detached eclipsing binary stars are presented for the first time. These were derived from the V-, and I-band ASAS and WASP photometry, and new radial velocities calculated from high quality optical spectra we obtained with a wide range of spectrographs and using the two-dimensional cross-correlation technique (TODCOR). All of the investigated systems (ASAS J184949-1518.7, BQ Aqr, and V1207 Cen) show the differential evolutionary phase of their components consisting of a main sequence star or a subgiant and a red giant, and thus constitute very informative objects in terms of testing stellar evolution models. Additionally, the systems show significant chromospheric activity of both components. They can be also classified as classical RS CVn-type stars. Besides the standard analysis of radial velocities and photometry, we applied spectral disentangling to obtain separate spectra for both components of each analysed system which allowed for a more detailed spectroscopic study. We also compared the properties of red giant stars in binaries that show spots, with those that do not, and found that the activity phenomenon is substantially suppressed for stars with Rossby number higher than $\sim$1 and radii larger than $\sim$20 R$_\odot$.

We present scientific results obtained with a recently commissioned \'{e}chelle spectrograph on the 0.5-m Solaris-1 telescope in the South African Astronomical Observatory. BACHES is a low-cost slit \'{e}chelle spectrograph that has a resolution of 21,000 at 5,500 \AA. The described setup is fully remotely operated and partly automated. Custom hardware components have been designed to allow both spectroscopic and photometric observations. The setup is controlled via dedicated software. The throughput of the system allows us to obtain spectra with an average SNR of 22 at 6375 {\AA} for a 30-min exposure of a $V=10$ mag target. The stability of the instrument is influenced mainly by the ambient temperature changes. We have obtained radial velocity RMS values for a bright (V = 5.9 mag) spectroscopic binary as good as 0.59 km s$^{-1}$ and 1.34 km s$^{-1}$ for a $V = 10.2$ mag eclipsing binary. Radial velocity measurements have been combined with available photometric light curves. We present models of six eclipsing binary systems, and for previously known targets, we compare our results with those available in the literature. Masses of binary components have been determined with 3% errors for some targets. We confront our results with benchmark values based on measurements from the HARPS and UCLES spectrographs on 4-m class telescopes and find very good agreement. The described setup is very efficient and well suited for a spectroscopic survey. We can now spectroscopically characterize about 300 eclipsing binary stars per year up to 10.2 mag assuming typical weather conditions at SAAO without a single observing trip.

We present our results of a combined analysis of radial velocity and light curves of a double-lined spectroscopic and eclipsing binary KIC 09246715, observed photometrically by the $Kepler$ satellite, and spectroscopically with the OAO-1.88m telescope with the HIgh-Dispertion Echelle Spectrograph (HIDES). The target was claimed to be composed of two red giants, one of which is showing solar-like oscillations. We have found that the mass and radius of the primary are $M_1=2.169\pm0.024$ M$_\odot$ and $R_1=8.47\pm0.13$ R$_\odot$, and of the secondary: $M_2=2.143\pm0.025$ M$_\odot$ and $R_2=8.18\pm0.09$ R$_\odot$, which confirms the double-giant status. Our secondary is the star to which the oscillations were attributed. Results of its previous asteroseismic analysis are in agreement with ours, only significantly less precise, but the subsequent light-curve-based study failed to derive correct mass and radius of our primary. KIC 09246715 is one of the rare cases where asteroseismic parameters of a solar-like oscillator were confirmed by an independent method, and only the third example of a Galactic double-giant eclipsing binary with masses and radii measured with precision below 2%.

For more than three years now we have been conducting a spectroscopic survey of detached eclipsing binaries (DEBs) from the All-Sky Automated Survey (ASAS) database. Thousands of high-resolution spectra of over 300 systems were secured, and used for radial velocity measurements and spectral analysis. In our sample we found a zoo of multiple systems, such as spectroscopic triples and quadruples, visual binaries with eclipsing components, and circumbinary low-mass companions, including sub-stellar-mass candidates

We present the orbital and physical parameters of the detached eclipsing binary V1200~Centauri (ASAS~J135218-3837.3) from the analysis of spectroscopic observations and light curves from the \textit{All Sky Automated Survey} (ASAS) and SuperWASP database. The radial velocities were computed from the high-resolution spectra obtained with the OUC 50-cm telescope and PUCHEROS spectrograph and with 1.2m Euler telescope and CORALIE spectrograph using the cross-correlation technique \textsc{todcor}. We found that the absolute parameters of the system are $M_1= 1.394\pm 0.030$ M$_\odot$, $M_2= 0.866\pm 0.015$ M$_\odot$, $R_1= 1.39\pm 0.15$ R$_\odot$, $R_2= 1.10\pm 0.25$ R$_\odot$. We investigated the evolutionary status and kinematics of the binary and our results indicate that V1200~Centauri is likely a member of the Hyades moving group, but the largely inflated secondary's radius may suggest that the system may be even younger, around 30 Myr. We also found that the eclipsing pair is orbited by another, stellar-mass object on a 351-day orbit, which is unusually short for hierarchical triples. This makes V1200 Cen a potentially interesting target for testing the formation models of multiple stars.

We present the first full orbital and physical analysis of HD 187669, recognized by the All-Sky Automated Survey (ASAS) as the eclipsing binary ASAS J195222-3233.7. We combined multi-band photometry from the ASAS and SuperWASP public archives and 0.41-m PROMPT robotic telescopes with our high-precision radial velocities from the HARPS spectrograph. Two different approaches were used for the analysis: 1) fitting to all data simultaneously with the WD code, and 2) analysing each light curve (with JKTEBOP) and RVs separately and combining the partial results at the end. This system also shows a total primary (deeper) eclipse, lasting for about 6 days. A spectrum obtained during this eclipse was used to perform atmospheric analysis with the MOOG and SME codes in order to constrain physical parameters of the secondary. We found that ASAS J195222-3233.7 is a double-lined spectroscopic binary composed of two evolved, late-type giants, with masses of $M_1 = 1.504\pm0.004$ and $M_2=1.505\pm0.004$ M$_\odot$, and radii of $R_1 = 11.33\pm0.28$ and $R_2=22.62\pm0.50$ R$_\odot$, slightly less metal abundant than the Sun, on a $P=88.39$ d orbit. Its properties are well reproduced by a 2.38 Gyr isochrone, and thanks to the metallicity estimation from the totality spectrum and high precision in masses, it was possible to constrain the age down to 0.1 Gyr. It is the first so evolved galactic eclipsing binary measured with such a good accuracy, and as such is a unique benchmark for studying the late stages of stellar evolution.

We evaluate a pre-production BACHES \'{e}chelle spectrograph in terms of its usefulness for radial velocity surveys of binary stars with small telescopes in a remote and autonomous. We use the Solaris-4 observatory located in Casleo, Argentina, that is part of a global network of autonomous observatories as the test-bed for the instrument. The setup is designed in such a way that spectroscopy and photometry can be carried out using the same telescope without the need to mechanically modify the imaging train. We observe single spectroscopic standard stars as well as binary stars up to 9.75 mag. We present results of mechanical tests of the instruments and spectroscopic observations carried out between Nov 26th and Dec 8th 2013. We conclude that BACHES is a very compact and capable spectrograph well suited for remote and autonomous operation. Coupled to a 0.5-m telescope it is capable of obtaining spectra of 10 mag targets with a SNR of 20 for 30-minute exposures. This is a very good result considering the price and size of the instrument. It brings new possibilities to the scientific community and opens a whole new range of research opportunities available to new and existing observatories.

We present the results of the combined photometric and spectroscopic analysis of a bright (V=9.14), nearby (d=31 pc), late-type detached eclipsing binary AK Fornacis. This P=3.981 d system has not been previously recognised as a double-lined spectroscopic binary, and this is the first full physical model of this unique target. With the FEROS, CORALIE and HARPS spectrographs we collected a number of high-resolution spectra in order to calculate radial velocities of both components of the binary. Measurements were done with our own disentangling procedure and the TODCOR technique, and were later combined with the photometry from the ASAS and SuperWASP archives. We also performed an atmospheric analysis of the component spectra with the Spectroscopy Made Easy (SME) package. Our analysis shows that AK For consists of two active, cool dwarfs having masses of $M_1=0.6958 \pm 0.0010$ and $M_2=0.6355 \pm 0.0007$ M$_\odot$ and radii of $R_1=0.687 \pm 0.020$ and $R_2=0.609 \pm 0.016$ R$_\odot$, slightly less metal abundant than the Sun. Parameters of both components are well reproduced by the models. AK For is the brightest system among the known eclipsing binaries with K or M type stars. Its orbital period is one of the longest and rotational velocities one of the lowest, which allows us to obtain very precise radial velocity measurements. The precision in physical parameters we obtained places AK For among the binaries with the best mass measurements in the literature. It also fills the gap in our knowledge of stars in the range of 0.5-0.8 M$_\odot$, and between short and long-period systems. All this makes AK For a unique benchmark for understanding the properties of low-mass stars.

We present absolute physical and orbital parameters for three double-lined detached eclipsing binary systems from the All Sky Automated Survey (ASAS) catalogue with subgiant and giant components. These parameters were derived from archival V-band ASAS photometry and new radial velocities. The radial velocities were calculated from high quality optical spectra we obtained with the 8.2 m Subaru/HDS, ESO 3.6 m/HARPS, 1.9 m Radcliffe/GIRAFFE, CTIO 1.5 m/CHIRON, and 1.2 m Euler/CORALIE using the two-dimensional cross-correlation technique (TODCOR) and synthetic template spectra chosen for every system separately as references. The physical and orbital parameters of the systems were derived with the PHOEBE and JKTEBOP codes. We checked the evolutionary status of the systems with several sets of isochrones and determined distances for each system. The derived uncertainties for individual masses of ASAS J010538-8003.7, ASAS J182510-2435.5, and V1980 Sgr components vary from 0.7% to 3.6% while the radii are in the range 1% to 24%. For all of the investigated systems such a detailed orbital and physical analysis is presented for the first time.

Current precision in radial velocity (RV) measurements of binary stars reaches $\sim$2 ms$^{-1}$. This level of precision means that RV models have to take into account additional non-Keplerian effects such as tidal and rotational distortion of the components of a binary star, relativistic effects and orbital precession. We generate synthetic binaries using Yonsei-Yale stellar models. For typical representatives we investigate the impact of various orbital orientations and different non-Keplerian effects on the RV curves. To this end we simulate RV observations with an added white noise of different scale. Subsequently we try to reconstruct the input orbital parameters and their errors by fitting a model using a standard least-squares method. In particular we investigate the connection between the tidal distortion of the shape of the stars and the best-fit orbital eccentricity, the possibility of deriving orbital inclination of a non-eclipsing binary star by exploiting relativistic effects and the circumstances in which the orbital precession can be detected. We confirm that the method proposed by to obtain orbital inclination with use of the relativistic effect does work in favourable cases and that it can be used even for orbital configurations far from an edge-on orientation. We show that the RV variations imposed by tidally distorted stars can mimic non-zero eccentricity in some binaries. The scale of such an effect depends on the RV accuracy. Finally, we demonstrate that the apsidal precession can be easily detected with precision RVs. In particular we can detect orbital precession of $10^{-4}$ rad yr$^{-1}$, $10^{-3}$ rad yr$^{-1}$ for precision of RVs of 1 ms$^{-1}$ and 10 ms$^{-1}$ respectively.

We present the orbital and physical parameters of a newly discovered low-mass detached eclipsing binary from the All-Sky Automated Survey (ASAS) database: ASAS J011328-3821.1 A - a member of a visual binary system with the secondary component separated by about 1.4 seconds of arc. The radial velocities were calculated from the high-resolution spectra obtained with the 1.9-m Radcliffe/GIRAFFE, 3.9-m AAT/UCLES and 3.0-m Shane/HamSpec telescopes/spectrographs on the basis of the TODCOR technique and positions of H_alpha emission lines. For the analysis we used V and I band photometry obtained with the 1.0-m Elizabeth and robotic 0.41-m PROMPT telescopes, supplemented with the publicly available ASAS light curve of the system. We found that ASAS J011328-3821.1 A is composed of two late-type dwarfs having masses of M_1 = 0.612 +/- 0.030 M_sun, M_2 = 0.445 +/- 0.019 M_sun and radii of R_1 = 0.596 +/- 0.020 R_sun, R_2 = 0.445 +/- 0.024 R_sun, both show a substantial level of activity, which manifests in strong H_alpha and H_beta emission and the presence of cool spots. The influence of the third light on the eclipsing pair properties was also evaluated and the photometric properties of the component B were derived. Comparison with several popular stellar evolution models shows that the system is on its main sequence evolution stage and probably is more metal rich than the Sun. We also found several clues which suggest that the component B itself is a binary composed of two nearly identical ~0.5 M_sun stars.

The Light-Time Effect (LTE) is observed whenever the distance between the observer and any kind of periodic event changes in time. The usual cause of this distance change is the reflex motion about the system's barycenter due to the gravitational influence of one or more additional bodies. We analyze 5032 eclipsing contact (EC) and detached (ED) binaries from the All Sky Automated Survey (ASAS) catalogue to detect variations in the times of eclipses which possible can be due to the LTE effect. To this end we use an approach known from the radio pulsar timing where a template radio pulse of a pulsar is used as a reference to measure the times of arrivals of the collected pulses. In our analysis as a template for a photometric time series from ASAS, we use a best-fitting trigonometric series representing the light curve of a given EC or ED. Subsequently, an O-C diagram is built by comparing the template light curve with light curves obtained from subsets of a given time series. Most of the variations we detected in O-Cs correspond to a linear period change. Three show evidence of more than one complete LTE-orbit. For these objects we obtained preliminary orbital solutions. Our results demonstrate that the timing analysis employed in radio pulsar timing can be effectively used to study large data sets from photometric surveys.

We present the most precise to date orbital and physical parameters of the well known short period (P=5.975 d), eccentric (e=0.3) double-lined spectroscopic binary BY Draconis, a prototype of a class of late-type, active, spotted flare stars. We calculate the full spectroscopic/astrometric orbital solution by combining our precise radial velocities (RVs) and the archival astrometric measurements from the Palomar Testbed Interferometer (PTI). The RVs were derived based on the high resolution echelle spectra taken between 2004 and 2008 with the Keck I/HIRES, Shane/CAT/HamSpec and TNG/SARG telescopes/spectrographs using our novel iodine-cell technique for double-lined binary stars. The RVs and available PTI astrometric data spanning over 8 years allow us to reach 0.2-0.5% level of precision in M sin3(i) and the parallax but the geometry of the orbit (i=154 deg) hampers the absolute mass precision to 3.3%, which is still an order of magnitude better than for previous studies. We compare our results with a set of Yonsei-Yale theoretical stellar isochrones and conclude that BY Dra is probably a main sequence system more metal-rich than the Sun. Using the orbital inclination and the available rotational velocities of the components, we also conclude that the rotational axes of the components are likely misaligned with the orbital angular momentum. Given BY Dra's main sequence status, late spectral type and the relatively short orbital period, its high orbital eccentricity and probable spin-orbit misalignment are not in agreement with the tidal theory. This disagreement may possibly be explained by smaller rotational velocities of the components and the presence of a substellar mass companion to BY Dra AB.

We present the results of our spectroscopic and photometric analysis of two newly discovered low-mass detached eclipsing binaries found in the All-Sky Automated Survey (ASAS) catalogue: ASAS J093814-0104.4 and ASAS J212954-5620.1. Using the GIRAFFE instrument on the 1.9-m Radcliffe telescope at SAAO and the UCLES spectrograph on the 3.9-m Anglo-Australian Telescope, we obtained high-resolution spectra of both objects and derived their radial velocities (RVs) at various orbital phases. The RVs of both objects were measured with the TODCOR technique using synthetic template spectra as references. We also obtained V and I band photometry using the 1.0-m Elizabeth telescope at SAAO and the 0.4-m PROMPT instruments located at the CTIO. The orbital and physical parameters of the systems were derived with PHOEBE and JKTEBOP codes. We compared our results with several sets of widely-used isochrones. Our multi-epoch photometric observations demonstrate that both objects show significant out-of-eclipse modulations, which vary in time. We believe that this effect is caused by stellar spots, which evolve on time scales of tens of days. For this reason, we constructed our models on the basis of photometric observations spanning short time scales (less than a month). Our modeling indicates that (1) ASAS-09 is a main sequence active system with nearly-twin components with masses of M1 = 0.771(33) Msun, M2 = 0.768(21) Msun and radii of R1 = 0.772(12) Rsun and R2 = 0.769(13) Rsun. (2) ASAS-21 is a main sequence active binary with component masses of M1 = 0.833(17) Msun, M2 = 0.703(13) Msun and radii of R1 = 0.845(12) Rsun and R2 = 0.718(17) Rsun. Both systems confirm the characteristic of active low-mass stars, for which the observed radii are larger and the temperatures lower than predicted by evolutionary models. Other parameters agree within errors with the models of main sequence stars.

The Palomar High-precision Astrometric Search for Exoplanet Systems monitored 51 subarcsecond binary systems to evaluate whether tertiary companions as small as Jovian planets orbited either the primary or secondary stars, perturbing their otherwise smooth Keplerian motions. Six binaries are presented that show evidence of substellar companions orbiting either the primary or secondary star. Of these six systems, the likelihoods of two of the detected perturbations to represent real objects are considered to be "high confidence", while the remaining four systems are less certain and will require continued observations for confirmation.

(Abridged) Differential astrometry measurements from the Palomar High-precision Astrometric Search for Exoplanet Systems (PHASES) are used to constrain the astrometric orbit of the previously known lesssim 2 day subsystem in the triple system 63 Gem A and have detected a previously unknown 2 year Keplerian wobble superimposed on the visual orbit of the much longer period (213 years) binary system HR 2896. The very small astrometric perturbation caused by the inner pair in 63 Gem A stretches the limits of current astrometric capabilities, but PHASES observations are able to constrain the orientation of the orbit. The two bright stars comprising the HR 2896 long period (213 year) system have a combined spectral type of K0III and the newly detected object's mass estimate places it in the regime of being a M dwarf. The motion of the stars are slow enough that their spectral features are always blended, preventing Doppler studies. The PHASES measurements and radial velocities (when available) have been combined with lower precision single-aperture measurements covering a much longer timeframe to improve the characterization of the long period orbits in both binaries. The visual orbits of the short and long period systems are presented for both systems, and used to calculate two possible values of the mutual inclinations between inner and outer orbits of 152 pm 12 degrees or a less likely value of 31 pm 11 degrees for 63 Gem A and 10.2 pm 2.4 degrees or 171.2 pm 2.8 degrees for HR 2896. The first is not coplanar, whereas the second is either nearly coplanar or anti-coplanar.

The Palomar High-precision Astrometric Search for Exoplanet Systems (PHASES) monitored 51 subarcsecond binary systems to evaluate whether tertiary companions as small as Jovian planets orbited either the primary or secondary stars, perturbing their otherwise smooth Keplerian motions. Twenty-one of those systems were observed 10 or more times and show no evidence of additional companions. A new algorithm is presented for identifying astrometric companions and establishing the (companion mass)-(orbital period) combinations that can be excluded from existence with high confidence based on the PHASES observations, and the regions of mass-period phase space being excluded are presented for 21 PHASES binaries.

Differential astrometry measurements from the Palomar High-precision Astrometric Search for Exoplanet Systems have been combined with lower precision single-aperture measurements covering a much longer timespan (from eyepiece measurements, speckle interferometry, and adaptive optics) to determine improved visual orbits for 20 binary stars. In some cases, radial velocity observations exist to constrain the full three-dimensional orbit and determine component masses. The visual orbit of one of these binaries---alpha Com (HD 114378)---shows that the system is likely to have eclipses, despite its very long period of 26 years. The next eclipse is predicted to be within a week of 2015 January 24.

The Palomar High-precision Astrometric Search for Exoplanet Systems (PHASES) monitored 51 sub-arcsecond binary systems to determine precision binary orbits, study the geometries of triple and quadruple star systems, and discover previously unknown faint astrometric companions as small as giant planets. PHASES measurements made with the Palomar Testbed Interferometer (PTI) from 2002 until PTI ceased normal operations in late 2008 are presented. Infrared differential photometry of several PHASES targets were measured with Keck Adaptive Optics and are presented.

We present photometric and spectroscopic analysis of the bright detached eclipsing binary BG Ind. The masses of the components are found to be 1.428 +- 0.008 and 1.293 +- 0.008 Msun and the radii to be 2.290+-0.017 and 1.680+-0.038 Rsun for primary and secondary stars, respectively. Spectra- and isochrone-fitting coupled with color indices calibration yield [Fe/H]=-0.2+-0.1. At an age of 2.65+-0.20 Gyr BG Ind is well advanced in the main-sequence evolutionary phase - in fact, its primary is at TAMS or just beyond it. Together with three similar systems (BK Peg, BW Aqr and GX Gem) it offers an interesting opportunity to test the theoretical description of overshooting in the critical mass range 1.2-1.5 Msun.

Aims. We present the results of our detailed spectroscopic and photometric analysis of two previously unknown < 1 M{_\sun} detached eclipsing binaries: ASAS J045304-0700.4 and ASAS J082552-1622.8. Methods. With the HIgh Resolution Echelle Spectrometer (HIRES) on the Keck-I telescope, we obtained spectra of both objects covering large fractions of orbits of the systems. We also obtained V and I band photometry with the 1.0-m Elizabeth telescope of the South African Astronomical Observatory (SAAO). The orbital and physical parameters of the systems were derived with the PHOEBE and JKTEBOP codes. We investigated the evolutionary status of both binaries with several sets of widely-used isochrones. Results. Our modelling indicates that (1) ASAS J045304-0700.4 is an old, metal-poor, active system with component masses of M1 = 0.8452 \pm 0.0056 M$_\odot$, M2 = 0.8390 \pm 0.0056 M$_\odot$ and radii of R1 = 0.848 \pm 0.005 R$_\odot$ and R2 = 0.833 \pm 0.005 R$_\odot$, which places it at the end of the Main Sequence evolution - a stage rarely observed for this type of stars. (2) ASAS J082552-1622.8 is a metal-rich, active binary with component masses of M1 = 0.703 \pm 0.003 M$_\odot$, M2 = 0.687 \pm 0.003 M$_\odot$ and radii of R1 = 0.694(+0.007-0.011) R$_\odot$ and R2 = 0.699(+0.011-0.014) R$_\odot$. Both systems show significant out-of-eclipse variations, probably owing to large, cold spots. we also investigated the influence of a third light in the second system.

The presence of a body in an orbit around a close eclipsing binary star manifests itself through the light time effect influencing the observed times of eclipses as the close binary and the circumbinary companion both move around the common centre of mass. This fact combined with the periodicity with which the eclipses occur can be used to detect the companion. Given a sufficient precision of the times of eclipses, the eclipse timing can be employed to detect substellar or even planetary mass companions. The main goal of the paper is to investigate the potential of the photometry based eclipse timing of binary stars as a method of detecting circumbinary planets. In the models we assume that the companion orbits a binary star in a circular Keplerian orbit. We analyze both the space and ground based photometry cases. In particular, we study the usefulness of the on-going COROT and Kepler missions in detecting circumbinary planets. We also explore the relations binding the planet discovery space with the physical parameters of the binaries and the geometrical parameters of their light curves. We carry out detailed numerical simulations of the eclipse timing by employing a relatively realistic model of the light curves of eclipsing binary stars. We study the influence of the white and red photometric noises on the timing precision. We determine the sensitivity of the eclipse timing technique to circumbinary planets for the ground and space based photometric observations. We provide suggestions for the best targets, observing strategies and instruments for the eclipse timing method. Finally, we compare the eclipse timing as a planet detection method with the radial velocities and astrometry.

First spectroscopic and new photometric observations of the eclipsing binary FM Leo are presented. The main aims were to determine orbital and stellar parameters of two components and their evolutionary stage. First spectroscopic observations of the system were obtained with DDO and PST spectrographs. The results of the orbital solution from radial velocity curves are combined with those derived from the light-curve analysis (ASAS-3 photometry and supplementary observations of eclipses with 1 m and 0.35 m telescopes) to derive orbital and stellar parameters. JKTEBOP, Wilson-Devinney binary modelling codes and a two-dimensional cross-correlation (TODCOR) method were applied for the analysis. We find the masses to be M_1 = 1.318 $\pm$ 0.007 and M_2 = 1.287 $\pm$ 0.007 M_sun, the radii to be R_1 = 1.648 $\pm$ 0.043 and R_2 = 1.511 $\pm$ 0.049 R_sun for primary and secondary stars, respectively. The evolutionary stage of the system is briefly discussed by comparing physical parameters with current stellar evolution models. We find the components are located at the main sequence, with an age of about 3 Gyr.

We present high precision radial velocities (RVs) of double-lined spectroscopic binary stars HD78418, HD123999, HD160922, HD200077 and HD210027. They were obtained based on the high resolution echelle spectra collected with the Keck I/Hires, Shane/CAT/Hamspec and TNG/Sarge telescopes/spectrographs over the years 2003-2008 as a part of TATOOINE search for circumbinary planets. The RVs were computed using our novel iodine cell technique for double-line binary stars. The precision of the RVs is of the order of 1-10 m/s. Our RVs combined with the archival visibility measurements from the Palomar Testbed Interferometer allow us to derive very precise spectroscopic/astrometric orbital and physical parameters of the binaries. In particular, we derive the masses, the absolute K and H band magnitudes and the parallaxes. The masses together with the absolute magnitudes in the K and H bands enable us to estimate the ages of the binaries. These RVs allow us to obtain some of the most accurate mass determinations of binary stars. The fractional accuracy in m*sin(i) only and hence based on the RVs alone ranges from 0.02% to 0.42%. When combined with the PTI astrometry, the fractional accuracy in the masses ranges in the three best cases from 0.06% to 0.5%. Among them, the masses of HD210027 components rival in precision the mass determination of the components of the relativistic double pulsar system PSRJ0737-3039. In the near future, for double-lined eclipsing binary stars we expect to derive masses with a fractional accuracy of the order of up to ~0.001% with our technique. This level of precision is an order of magnitude higher than of the most accurate mass determination for a body outside the Solar System - the double neutron star system PSRB1913+16.

We present preliminary results of the first and on-going radial velocity survey for circumbinary planets. With a novel radial velocity technique employing an iodine absorption cell we achieve an unprecedented RV precision of up to 2 m/s for double-lined binary stars. The high resolution spectra collected with the Keck I/Hires, TNG/Sarg and Shane/CAT/Hamspec telescopes/spectrographs over the years 2003-2008 allow us to derive RVs and compute planet detection limits for ten double-lined binary stars. For this initial sample of targets, we can rule out planets on dynamically stable orbits with masses as small as ~0.3 to 3 MJup for the orbital periods of up to ~5.3 years. Even though the presented sample of stars is too small to make any strong conclusions, it is clear that the search for circumbinary planets is now technique-wise possible and eventually will provide new constraints for the planet formation theories.

Using the adaptive optics facilities at the 200-in Hale and 10-m Keck II, we observed in the near infrared a sample of 12 binary and multiple stars and one open cluster. We used the near diffraction limited images of these systems to measure the relative separations and position angles between their components. In this paper, we investigate and correct for the influence of the differential chromatic refraction and chip distortions on our relative astrometric measurements. Over one night, we achieve an astrometric precision typically well below 1 miliarcsecond and occasionally as small as 40 microarcseconds. Such a precision is in principle sufficient to astrometrically detect planetary mass objects around the components of nearby binary and multiple stars. Since we have not had sufficiently large data sets for the observed sample of stars to detect planets, we provide the limits to planetary mass objects based on the obtained astrometric precision.

We derive the absolute physical and orbital parameters for a sample of 18 detached eclipsing binaries from the \emph{All Sky Automated Survey} (ASAS) database based on the available photometry and our own radial velocity measurements. The radial velocities (RVs) are computed using spectra we collected with the 3.9-m Anglo-Australian Telescope and its \emph{University College London Echelle Spectrograph} and the 1.9-m SAAO Radcliffe telescope and its \emph{Grating Instrument for Radiation Analysis with a Fibre Fed Echelle}. In order to obtain as precise RVs as possible, most of the systems were observed with an iodine cell available at the AAT/UCLES and/or analyzed using the two-dimensional cross-correlation technique (TODCOR). The RVs were measured with TODCOR using synthetic template spectra as references. However, for two objects we used our own approach to the tomographic disentangling of the binary spectra to provide observed template spectra for the RV measurements and to improve the RV precision even more. For one of these binaries, AI Phe, we were able to the obtain an orbital solution with an RV $rms$ of 62 and 24 m s$^{-1}$ for the primary and secondary respectively. For this system, the precision in $M \sin^3{i}$ is 0.08%. For the analysis, we used the photometry available in the ASAS database. We combined the RV and light curves using PHOEBE and JKTEBOP codes to obtain the absolute physical parameters of the systems. Having precise RVs we were able to reach $\sim$0.2 % precision (or better) in masses in several cases but in radii, due to the limited precision of the ASAS photometry, we were able to reach a precision of only 1% in one case and 3-5 % in a few more cases. For the majority of our objects, the orbital and physical analysis is presented for the first time.

V838 Mon erupted at the beginning of 2002 becoming an extremely luminous star. Among various scenarios proposed to explain the nature of the outburst the most promising is a stellar merger event. In this paper we investigate the observational properties of the star and its surroundings in the post outburst phase. We have obtained a high resolution optical spectrum of V838 Mon in October 2005 using the Keck I telescope.We have identified numerous atomic features and molecular bands present in the spectrum and provided an atlas of those features. In order to improve the spectrum interpretation we have performed simple modeling of the molecular bands. Our analysis indicates that the spectrum is dominated by molecular absorption features arising in photospheric regions with temperatures of ~2400 K and in colder outer layers, where the temperature decreases down to ~500 K. A number of resonance lines of neutral alkali metals are observed to show P-Cyg profiles. Particularly interesting are numerous prominent emission lines of [FeII]. All of them show practically the same profile, which can be well described by a Lorentzian profile. In the blue part of the spectrum photospheric signatures of the B-type companion are easily seen. We have fitted the observed spectrum with a synthetic one and the obtained parameters are consistent with the B3V type. We have also estimated radial and rotational velocities of the companion.

We report the discovery of one or more planet-mass companions to the K0-giant HD 102272 with the Hobby-Eberly Telescope. In the absence of any correlation of the observed periodicities with the standard indicators of stellar activity, the observed radial velocity variations are most plausibly explained in terms of a Keplerian motion of at least one planet-mass body around the star. With the estimated stellar mass of 1.9M$_\odot$, the minimum mass of the confirmed planet is 5.9M$_J$. The planet's orbit is characterized by a small but nonzero eccentricity of $e$=0.05 and the semi-major axis of 0.61 AU, which makes it the most compact one discovered so far around GK-giants. This detection adds to the existing evidence that, as predicted by theory, the minimum size of planetary orbits around intermediate-mass giants is affected by both planet formation processes and stellar evolution. The currently available evidence for another planet around HD 102272 is insufficient to obtain an unambiguous two-orbit solution.

We present the results of our study of astrometric stability of 200-in Hale (Mt. Palomar) and 10-m Keck II (Mauna Kea) telescopes, both with Adaptive Optics (AO) facilities. A group of nearby visual binaries and multiples was observed in near infrared, relative separations and position angles measured. We have also checked the influence of some systematic effects (e.g. atmospherical refraction, varying plate scale factor) on result and precision of astrometric measurements. We conclude that in visual binaries astrometrical observations it is possible to achieve much better precision than 1 miliarcsecond, which in many cases allows detection of the astrometrical signal produced by planetary-mass object.

Precision radial velocity (RV) measurements of the Sun-like dwarf 14 Herculis published by Naef et. al (2004), Butler et. al (2006) and Wittenmyer et al (2007) reveal a Jovian planet in a 1760 day orbit and a trend indicating the second distant object. On the grounds of dynamical considerations, we test a hypothesis that the trend can be explained by the presence of an additional giant planet. We derive dynamical limits to th orbital parameters of the putative outer Jovian companion in an orbit within ~13 AU. In this case, the mutual interactions between the Jovian planets are important for the long-term stability of the system. The best self-consistent and stable Newtonian fit to an edge-on configuration of Jovian planets has the outer planet in 9 AU orbit with a moderate eccentricity ~0.2 and confined to a zone spanned by the low-order mean motion resonances 5:1 and 6:1. This solution lies in a shallow minimum of \Chi and persists over a wide range of the system inclination. Other stable configurations within 1\sigma confidence interval of the best fit are possible for the semi-major axis of the outer planet in the range of (6,13) AU and the eccentricity in the range of (0,0.3). The orbital inclination cannot yet be determined but when it decreases, both planetary masses approach ~10 Jupiter masses and for ~30 deg the hierarchy of the masses is reversed.

We have used high precision differential astrometry from the Palomar High-precision Astrometric Search for Exoplanet Systems (PHASES) project and radial velocity measurements covering a time-span of 20 years to determine the orbital parameters of the 88 Tau A system. 88 Tau is a complex hierarchical multiple system comprising a total of six stars; we have studied the brightest 4, consisting of two short-period pairs orbiting each other with an 18-year period. We present the first orbital solution for one of the short-period pairs, and determine the masses of the components and distance to the system to the level of a few percent. In addition, our astrometric measurements allow us to make the first determination of the mutual inclinations of the orbits. We find that the sub-systems are not coplanar.

mu Orionis was identified by spectroscopic studies as a quadruple star system. Seventeen high precision differential astrometry measurements of mu Ori have been collected by the Palomar High-precision Astrometric Search for Exoplanet Systems (PHASES). These show both the motion of the long period binary orbit and short period perturbations superimposed on that caused by each of the components in the long period system being themselves binaries. The new measurements enable the orientations of the long period binary and short period subsystems to be determined. Recent theoretical work predicts the distribution of relative inclinations between inner and outer orbits of hierarchical systems to peak near 40 and 140 degrees. The degree of coplanarity of this complex system is determined, and the angle between the planes of the A-B and Aa-Ab orbits is found to be 136.7 +/- 8.3 degrees, near the predicted distribution peak at 140 degrees; this result is discussed in the context of the handful of systems with established mutual inclinations. The system distance and masses for each component are obtained from a combined fit of the PHASES astrometry and archival radial velocity observations. The component masses have relative precisions of 5% (component Aa), 15% (Ab), and 1.4% (each of Ba and Bb). The median size of the minor axes of the uncertainty ellipses for the new measurements is 20 micro-arcseconds. Updated orbits for delta Equulei, kappa Pegasi, and V819 Herculis are also presented.

Searches for planets in close binary systems explore the degree to which stellar multiplicity inhibits or promotes planet formation. There is a degeneracy between planet formation models when only systems with single stars are studied--several mechanisms appear to be able to produce such a final result. This degeneracy is lifted by searching for planets in binary systems; the resulting detections (or evidence of non-existence) of planets in binaries isolates which models may contribute to how planets form in nature. In this chapter, we consider observational efforts to detect planetary companions to binary stars in two types of hierarchical planet-binary configurations: first ``S-type'' planets which orbit just one of the stars, with the binary period being much longer than the planet's; second, ``P-type'' or circumbinary planets, where the planet simultaneously orbits both stars, and the planetary orbital period is much longer than that of the binary. The S-type planet finding techniques are different for binaries that can or cannot be spatially resolved. For wider systems, techniques reviewed include dualstar interferometric differential astrometry and precision radial velocities. Alternatively, unresolved binaries can be studied using modified dualstar "PHASES-style" differential astrometry or a modification of the radial velocity technique for composite spectra. Should a fortunately aligned--but still long period--binary be found, eclipse timing can also reveal the presence of S-type planets. Methods for detecting P-type planets include the composite-spectra variant of the radial velocity technique and eclipse timing.

We report the discovery of a substellar-mass companion to the K0-giant HD 17092 with the Hobby-Eberly Telescope. In the absence of any correlation of the observed 360-day periodicity with the standard indicators of stellar activity, the observed radial velocity variations are most plausibly explained in terms of a Keplerian motion of a planetary-mass body around the star. With the estimated stellar mass of 2.3Msun, the minimum mass of the planet is 4.6MJ. The planet's orbit is characterized by a mild eccentricity of e=0.17 and a semi-major axis of 1.3 AU. This is the tenth published detection of a planetary companion around a red giant star. Such discoveries add to our understanding of planet formation around intermediate-mass stars and they provide dynamical information on the evolution of planetary systems around post-main sequence stars.

A new observing mode for the Palomar Testbed Interferometer was developed in2002-2003 which enables differential astrometry at the level of 20 micro-arcseconds for binary systems with separations of several hundred milli-arcseconds (mas). This phase-referenced mode is the basis of the Palomar High-precision Astrometric Search for Exoplanet Systems (PHASES), a search for giant planets orbiting either the primary or secondary star in fifty binary systems. We present the first science results from the PHASES search. The properties of the stars comprising binary systems are determined to high precision. The mutual inclinations of several hierarchical triple star systems have been determined. We will present upper limits constraining the the existence of giant planets in a few of the target systems.

Two nearby stars, HD 128311 and HD 82943, are believed to host pairs of Jupiter-like planets involved in a strong first order 2:1 mean motion resonance (MMR). In this work we reanalyze available radial velocity (RV) measurements and demonstrate that it is also possible to explain the observed RV variations of the parents stars as being induced by a pair of Trojan planets (i.e., in a 1:1 MMR). We show that these Trojan configurations reside in extended zones of stability in which such systems may easily survive in spite of large masses of the planets, large eccentricities and nonzero mutual inclinations of their orbits. We also show that HD 82943 could harbor a previously unknown third planet about of 0.5 Jupiter masses in 2.1 AU orbit.

We perform a dynamical analysis of the recently published radial velocity (RV) measurements of a few solar type stars which host multiple Jupiter-like planets. In particular, we re-analyze the data for HD 202206, 14 Her, HD 37124 and HD 108874. We derive dynamically stable configurations which reproduce the observed RV signals using our method called GAMP (an acronym of the Genetic Algorithm with MEGNO Penalty). The GAMP relies on the N-body dynamics and makes use of genetic algorithms merged with a stability criterion. For this purpose, we use the maximal Lyapunov exponent computed with the dynamical fast indicator MEGNO. Through a dynamical analysis of the phase-space in a neighborhood of the obtained best-fit solutions, we derive meaningful limits on the parameters of the planets. We demonstrate that GAMP is especially well suited for the analysis of the RV data which only partially cover the longest orbital period and/or correspond to multi-planet configurations involved in low-order mean motion resonances (MMRs). In particular, our analysis reveals a presence of a second Jupiter-like planet in the 14 Her system (14 Her c) involved in a 3:1 or 6:1 MMR with the known companion b. We also show that the dynamics of the HD 202206 system may be qualitatively different when coplanar and mutually-inclined orbits of the companions are considered. We demonstrate that the two outer planets in the HD 37124 system may reside in a close neighborhood of the 5:2 MMR. Finally, we found a clear indication that the HD 108874 system may be very close to, or locked in an exact 4:1 MMR.

V819 Herculis is a well-studied triple star system consisting of a ``wide'' pair with 5.5 year period, one component of which is a 2.2-day period eclipsing single-line spectroscopic binary. Differential astrometry measurements from the Palomar High-precision Astrometric Search for Exoplanet Systems (PHASES) are presented and used to determine a relative inclination between the short- and long-period orbits of 23.6 +- 4.9 degrees. This represents only the sixth unambiguous determination of the mutual inclination of orbits in a hierarchical triple system. This result is combined with those for the other five systems for analysis of the observed distribution of mutual inclinations in nearby triple systems. It is found that this distribution is different than that which one would expect from random orientations with statistical significance at the 94% level; implications for studying the spatial distribution of angular momentum in star forming regions is discussed.

kappa Pegasi is a well-known, nearby triple star system. It consists of a ``wide'' pair with semi-major axis 235 milli-arcseconds, one component of which is a single-line spectroscopic binary (semi-major axis 2.5 milli-arcseconds). Using high-precision differential astrometry and radial velocity observations, the masses for all three components are determined and the relative inclinations between the wide and narrow pairs' orbits is found to be 43.8 +/- 3.0 degrees, just over the threshold for the three body Kozai resonance. The system distance is determined to 34.60 +/- 0.21 parsec, and is consistent with trigonometric parallax measurements.

delta Equulei is among the most well-studied nearby binary star systems. Results of its observation have been applied to a wide range of fundamental studies of binary systems and stellar astrophysics. It is widely used to calibrate and constrain theoretical models of the physics of stars. We report 27 high precision differential astrometry measurements of delta Equulei from the Palomar High-precision Astrometric Search for Exoplanet Systems (PHASES). The median size of the minor axes of the uncertainty ellipses for these measurements is 26 micro-arcseconds. These data are combined with previously published radial velocity data and other previously published differential astrometry measurements using other techniques to produce a combined model for the system orbit. The distance to the system is determined to within a twentieth of a parsec and the component masses are determined at the level of a percent. The constraints on masses and distance are limited by the precisions of the radial velocity data; we outline plans improve this deficiency and discuss the outlook for further study of this binary.

A spectroscopic technique employing an iodine absorption cell (I_2) to superimpose a reference spectrum onto a stellar spectrum is currently the most widely adopted approach to obtain precision radial velocities of solar-type stars. It has been used to detect ~80 extrasolar planets out of ~130 know. Yet in its original version, it only allows us to measure precise radial velocities of single stars. In this paper, we present a novel method employing an I_2 absorption cell that enables us to accurately determine radial velocities of both components of double-lined binaries. Our preliminary results based on the data from the Keck I telescope and HIRES spectrograph demonstrate that 20-30 m/s radial velocity precision can be routinely obtained for "early" type binaries (F3-F8). For later type binaries, the precision reaches ~10 m/s. We discuss applications of the technique to stellar astronomy and searches for extrasolar planets in binary systems. In particular, we combine the interferometric data collected with the Palomar Testbed Interferometer with our preliminary precision velocities of the spectroscopic double-lined binary HD 4676 to demonstrate that with such a combination one can routinely obtain masses of the binary components accurate at least at the level of 1.0%.

We report a transiting extrasolar giant planet around the star OGLE-TR-10 (orbital period = 3.1 days), which was uncovered as a candidate by the OGLE team in their photometric survey towards the Galactic center (Udalski et al, 2002a). We observed OGLE-TR-10 spectroscopically over a period of two years (2002-2004) using the HIRES instrument with an iodine cell on the Keck I telescope, and measured small radial velocity variations that are consistent with the presence of a planetary companion. This confirms the earlier identification of OGLE-TR-10b by our team and also recently by Bouchy et al (2004a) as a possible planet. Additionally, in this paper we are able to rule out a blend scenario as an alternative explanation. From an analysis combining all available radial velocity measurements with the OGLE light curve we find that OGLE-TR-10b has a mass of 0.57+/-0.12 MJup and a radius of 1.24+/-0.09 RJup. These parameters bear close resemblance to those of the first known transiting extrasolar planet, HD 209458b.

We perform a dynamical analysis of the recently updated set of the radial velocity (RV) measurements of the HD160691 (mu Arae). The kinematic, 2-Keplerian model of the measurements leads to the best-fit solution in which the eccentricity of the outer planet is about 0.7 and its semi-major axis is about 4AU. The parameters of the inner planet are well determined. The eccentricity is about 0.3 and the semi-major axis is about 1.65AU. The best 2-Keplerian model leads to a catastrophically unstable configuration, self disrupting in less than 20,000yr. To derive dynamically stable configurations which are simultaneously consistent with the RV data, we use the so called GAMP (Genetic Algorithm with MEGNO Penalty). Using this method, we derive meaningful limits on the parameters of the outer planet which provide a stable behavior of the system. The best-fit solutions are located in a shallow valley of Chi^2, in the (a_c,e_c)-plane, extending over 2AU (for the formal 1sigma confidence interval of the best fit). We find two equally good best-fit solutions leading to the qualitatively different orbital configurations. One of them corresponds to the center of the 5:1 mean motion resonance (MMR) and the second one describes a configuration between the 6:1 and 17:2MMRs. The mu Arae system can be found in the zone confined to other low-order MMRs of the type p:1 with p>5. Our results support the classification of the mu Arae as a hierarchical planetary system, dynamically similar to other known multi-planet systems around HD12661 and HD169830. The results of the GAMP on the extended data set are fully consistent with our previous conclusions concerning a much shorter observational window.

We study the long-term dynamics of the PSR 1257+12 planetary system. Using the recently determined accurate initial condition by Konacki & Wolszczan (2003) who derived the orbital inclinations and the absolute masses of the planets B and C, we investigate the system stability by long-term, 1Gyr direct integrations. No secular changes of the semi-major axes eccentricities and inclinations appear during such an interval. This stable behavior is confirmed with the fast indicator MEGNO. The analysis of the orbital stability in the neighborhood of the nominal initial condition reveals that the PSR 1257+12 system is localized in a wide stable region of the phase space but close to a few weak 2 and 3-body mean motion resonances. The long term stability is additionally confirmed by a negligible exchange of the Angular Momentum Deficit between the innermost planet A and the pair of the outer planets B and C. An important feature of the system that helps sustain the stability is the secular apsidal resonance (SAR) between the planets B and C with the center of libration about 180deg. We find useful limits on the elements of the innermost planet A which are otherwise unconstrained by the observations. We also perform a preliminary study of the short-term dynamics of massless particles in the system.

We re-examine the photometric and spectroscopic evidence available for the star OGLE-TR-56, recently discovered to harbor a giant planet presenting transits and orbiting with a period of 1.21 days. We investigate the possibility that the observational signatures reported might be the result of "blending" with the light of an eclipsing binary along the same line of sight. Using techniques developed earlier we perform fits to the light curve under a variety of blend configurations, subject to all observational constraints, and we infer further properties of these possible blends. We then carry out realistic end-to-end simulations based on those properties in order to quantify the spectral line asymmetries and radial velocity variations expected from such scenarios. We confront these calculations with the observations. The results from these models are clearly inconsistent with the measured radial velocity and bisector span variations, ruling out blends and confirming the planetary nature of the companion. The example of OGLE-TR-56 serves to illustrate the sort of tests that can and should be performed on transiting planet candidates to eliminate the possibility of "false positives", and in particular, line-of-sight contamination from eclipsing binaries.

The Palomar Testbed Interferometer has observed several binary star systems whose separations fall between the interferometric coherence length (a few hundredths of an arcsecond) and the typical atmospheric seeing limit of one arcsecond. Using phase-referencing techniques we measure the relative separations of the systems to precisions of a few tens of micro-arcseconds. We present the first scientific results of these observations, including the astrometric detection of the faint third stellar component of the kappa Pegasi system.

We report high-resolution spectroscopic follow-up observations of the faint transiting planet candidate OGLE-TR-33 (V = 14.7), located in the direction of the Galactic center. Small changes in the radial velocity of the star were detected that suggested initially the presence of a large planet or brown dwarf in orbit. However, further analysis revealed spectral line asymmetries that change in phase with the 1.95-day period, casting doubt on those measurements. These asymmetries make it more likely that the transit-like events in the light curve are the result of contamination from the light of an eclipsing binary along the same line of sight (referred to as a "blend"). We performed detailed simulations in which we generated synthetic light curves resulting from such blend scenarios and fitted them to the measured light curve. Guided by these fits and the inferred properties of the stars, we uncovered a second set of lines in our spectra that correspond to the primary of the eclipsing binary and explain the asymmetries. Using all the constraints from spectroscopy we were then able to construct a model that satisfies all the observations, and to characterize the three stars based on model isochrones. OGLE-TR-33 is fully consistent with being a hierarchical triple system composed of a slightly evolved F6 star (the brighter object) near the end of its main-sequence phase, and an eclipsing binary with a K7-M0 star orbiting an F4 star. The application to OGLE-TR-33 of the formalism developed to fit light curves of transit candidates illustrates the power of such simulations for predicting additional properties of the blend and for guiding further observations that may serve to confirm that scenario, thereby ruling out a planet. Tests such as this can be very important for validating faint candidates.

We report the independent discovery of a new extrasolar transiting planet around OGLE-TR-113, a candidate star from the Optical Gravitational Lensing Experiment. Small radial-velocity variations have been detected based on observations conducted with the MIKE spectrograph on the Magellan I (Baade) telescope at the Las Campanas Observatory (Chile) during 2003. We have also carried out a light-curve analysis incorporating new photometry and realistic physical parameters for the star. OGLE-TR-113b has an orbital period of only 1.43 days, a mass of 1.08 +/- 0.28 M_Jup, and a radius of 1.09 +/- 0.10 R_Jup. Similar parameters have been obtained very recently in an independent study by Bouchy et al., from observations taken a year later. The orbital period of OGLE-TR-113b, and also that of the previously announced planet OGLE-TR-56b (P_orb = 1.21 days) -- the first two found photometrically -- are much shorter than the apparent cutoff of close-in giant planets at 3-4-day periods found from high-precision radial velocities surveys. Along with a third case reported by Bouchy et al. (OGLE-TR-132b, P_orb = 1.69 days), these objects appear to form a new class of ``very hot Jupiters'' that pose very interesting questions for theoretical study.

We present optical interferometric observations of two double-lined spectroscopic binaries, HD 6118 and HD 27483, taken with the Palomar Testbed Interferometer (PTI) in the K band. HD 6118 is one of the most eccentric spectroscopic binaries and HD 27483 a spectroscopic binary in the Hyades open cluster. The data collected with PTI in 2001-2002 allow us to determine astrometric orbits and when combined with the radial velocity measurements derive all physical parameters of the systems. The masses of the components are 2.65 +/- 0.27 M_Sun and 2.36 +/- 0.24 M_Sun for HD 6118 and 1.38 +/- 0.13 M_Sun and 1.39 +/- 0.13 M_Sun for HD 27483. The apparent semi-major axis of HD 27483 is only 1.2 mas making it the closest binary successfully observed with an optical interferometer.

The new 2-planetary system around HD169830 has been announced during the XIX-th IAP Colloquium "Extrasolar Planets: Today & Tomorrow" (Paris, June 30 - July 4, 2003) by the Geneva Extrasolar Planet Search team. We study the orbital dynamics of this system in the framework of the $N$-body problem. The analysis of its orbital stability is performed using the long-term integrations and the fast indicators, the Mean Exponential Growth factor of Nearby Orbits and the Frequency Map Analysis. The HD169830 appears to be located in a wide stable region of the phase space. The ratio of the mean motions of the planets HD169830b and c is between low-order mean motion resonances, 9:1 and 10:1. The long-term integration of the coplanar configurations, conducted over 1Gyr, reveals that the eccentricities of the companions vary with a large amplitude about 0.4-0.5 but there is no sign of instability. The orbital parameters of the planets resemble those of another 2-planetary system, around HD12661. Both of them can be classified as hierarchical planetary systems. We investigate whether these two exosystems are dynamically similar. Such similarities may be important for finding out if the formation and subsequent orbital evolution of exoplanetary systems obey common rules.

We report new spectroscopic observations of the recently discovered transiting planet OGLE-TR-56b with the Keck/HIRES instrument. Our radial velocity measurements with errors of ~100 m/s show clear variations that are in excellent agreement with the phasing (period and epoch) derived from the OGLE transit photometry, confirming the planetary nature of the companion. The new data combined with measurements from the previous season allow an improved determination of the mass of the planet, M_p = 1.45 +/- 0.23 M_Jup. All available OGLE photometry, including new measurements made this season, have also been analyzed to derive an improved value for the planetary radius of R_p = 1.23 +/- 0.16 R_Jup. We discuss the implications of these results for the theory of extrasolar planets.

We present new results of timing and single pulse measurements for 18 radio pulsars discovered in 1993 - 1997 by the Penn State/NRL declination-strip survey conducted with the 305-m Arecibo telescope at 430 MHz. Long-term timing measurements have led to significant improvements of the rotational and the astrometric parameters of these sources, including the millisecond pulsar, PSR J1709+2313, and the pulsar located within the supernova remnant S147, PSR J0538+2817. Single pulse studies of the brightest objects in the sample have revealed an unusual "bursting" pulsar, PSR J1752+2359, two new drifting subpulse pulsars, PSR J1649+2533 and PSR J2155+2813, and another example of a pulsar with profile mode changes, PSR J1746+2540. PSR J1752+2359 is characterized by bursts of emission, which appear once every 3-5 min. and decay exponentially on a ~45 sec timescale. PSR J1649+2533 spends ~30% of the time in a null state with no detectable radio emission.

We report the results of our campaign to follow-up spectroscopically several candidate extrasolar transiting planets from the OGLE-III survey in the direction of the Galactic center, announced in 2001 (Udalski et al 2002a,b). All of these objects present shallow and periodic dips in brightness that may be due to planetary companions. Our high-resolution Keck I/HIRES observations have revealed two interesting cases (OGLE-TR-10 with a period of 3.1 days, and OGLE-TR-58 with a period of 4.3 days) that show no radial velocity variations at the level of 100--200 m/s. If orbited by companions, their masses would be similar to Jupiter. With the information in hand (including the light curves) we are not able to rule out that these candidates are instead the result of contamination from an eclipsing binary in the same line of sight (a "blend"). We discuss also the case of OGLE-TR-56 that was recently reported by Konacki et al (2003) to have a Jupiter-size companion, based on an earlier analysis of our data, and we present supporting information. Two other candidates, OGLE-TR-3 and OGLE-TR-33, show clear evidence that they are blends.

We present measurements of the true masses and orbital inclinations of the two Earth-mass planets in the PSR B1257+12 system, based on the analysis of their mutual gravitational perturbations detectable as microsecond variations of the arrival times of radio pulses from the pulsar. The 6.2-millisecond pulsar, PSR B1257+12, has been regularly timed with the Arecibo telescope since late 1990. Assuming the standard pulsar mass of 1.4 M_Sun, the derived masses of planets B and C are 4.3 +/- 0.2 M_Earth and 3.9 +/- 0.2 M_Earth, respectively. The corresponding orbital inclinations of 53.4 and 47.3 deg (or 127 and 133 deg) imply that the two orbits are almost coplanar. This result, together with the known near 3:2 resonance between the orbits of the two planets, strongly supports the hypothesis of a disk origin of the PSR B1257+12 planetary system. The system's long-term stability is guaranteed by the low, Earth-like masses of planets B and C.

The set of radial velocity measurements of the HD 160691 has been recently published by Jones et.al 2002 (MNRAS). It reveals a linear trend that indicates a presence of the second planet in this system. The preliminary double-Keplerian orbital fit to the observations, announced by the discovery team, describes a highly unstable, self-disrupting configuration. Because the observational window of the HD 160691 system is narrow, the orbital parameters of the hypothetical second companion are unconstrained. In this paper we try to find out whether a second giant planet can exist up to the distance of Jupiter and search for the dynamical constraints on its orbital parameters. Our analysis employs a combination of fitting algorithms and simultaneous examination of the dynamical stability of the obtained orbital fits. It reveals that if the semi-major axis of the second planet is smaller than $\simeq 5.2$ AU, the observations are consistent with quasi-periodic, regular motions of the system confined to the islands of various low-order mean motion resonances, e.g., 3:1, 7:2, 4:1, 5:1, or to their vicinity. In such cases the second planet has smaller eccentricity $\simeq 0.2-0.5$ than estimated in the previous works. We show that the currently available Doppler data rather preclude the 2:1 mean motion resonance expected by some authors to be present in the HD 160691 planetary system. We also demonstrate that the MEGNO-penalty method, developed in this paper, which is a combination of the genetic minimization algorithm and the MEGNO analysis, can be efficiently used for predicting stable planetary configurations when only a limited number of observations is given or the data do not provide tight constraints on the orbital elements.

We have measured radial velocities of a star, OGLE-TR-56, which shows a 1.2-day transit-like light curve found photometrically by Udalski et al.(2002ab). Here we show that the velocity changes we detect are probably induced by an object of 0.9 Jupiter masses - a very close-in gas-giant planet only 0.023 AU from its star, with a planetary radius of 1.3 Jupiter radii and a mean density of 0.5 g/cm3. At its small orbital distance, OGLE-TR-56b is hotter than any known planet, approaching 1900K, but it is stable against long-term evaporation or tidal disruption.

We present a new analysis of the dynamics of the planetary system around the pulsar B1257+12. A semi-analytical theory of perturbation between terrestrial-mass planets B and C is developed and applied to improve multi-orbit timing formula for this object. We use numerical simulations of the pulse arrival times for PSR B1257+12 to demonstrate that our new timing model can serve as a toll to determine the masses of the two planets.

We present theoretical analysis of the astrometric searches for extrasolar planets with the Space Interferometry Mission (SIM). Particularly, we derive a model for the future measurements with SIM and discuss the problem of reliable estimation of orbital elements of planets. For this purpose we propose a new method of data analysis and present a numerical test of its application on simulated SIM astrometric measurements of $\upsilon$ Andromedae planetary system. We demonstrate that our approach allows successfull determination of its orbital elements.

Since it has become possible to discovery planets orbiting nearby solar-type stars through very precise Doppler-shift measurements, the role of methods used to analyze such observations has grown significantly. The widely employed model-dependent approach based on the least-squares fit of the Keplerian motion to the radial-velocity variations can be, as we show, unsatisfactory. Thus, in this paper, we propose a new method that may be easily and successfully applied to the Doppler-shift measurements. This method allows us to analyze the data without assuming any specific model and yet to extract all significant features of the observations. This very simple idea, based on the subsequent subtraction of all harmonic components from the data, can be easily implemented. We show that our method can be used to analyze real 16 Cygni B Doppler-shift observations with a surprising but correct result which is substantially different from that based on the least-squares fit of a Keplerian orbit. Namely, using frequency analysis we show that with the current accuracy of this star's observations it is not possible to determine the value of the orbital eccentricity which is claimed to be as high as 0.6.

In this paper we present a new method that can be used for analysis of time of arrival of a pulsar pulses (TOAs). It is designated especially to detect quasi-periodic variations of TOAs. We apply our method to timing observations of PSR B1257+12 and demonstrate that using it it is possible to detect not only first harmonics of a periodic variations, but also the presence of a resonance effect. The resonance effect detected, independently of its physical origin, can appear only when there is a non-linear interaction between two periodic modes. The explanation of TOAs variations as an effect of the existence of planets is, till now, the only known and well justified. In this context, the existence of the resonance frequency in TOAs is the most significant signature of the gravitational interaction of planets.