Automated analysis of eclipsing binary light curves – II. Statistical analysis of OGLE LMC eclipsing binaries

In the first paper of this series, we presented EBAS – Eclipsing Binary Automated Solver, a new fully automated algorithm to analyse the light curves of eclipsing binaries, based on the ebop code. Here, we apply the new algorithm to the whole sample of 2580 binaries found in the Optical Gravitational Lensing Experiment (OGLE) Large Magellanic Cloud (LMC) photometric survey and derive the orbital elements for 1931 systems. To obtain the statistical properties of the short-period binaries of the LMC, we construct a well-defined subsample of 938 eclipsing binaries with main-sequence B-type primaries. Correcting for observational selection effects, we derive the distributions of the fractional radii of the two components and their sum, the brightness ratios and the periods of the short-period binaries. Somewhat surprisingly, the results are consistent with a flat distribution in log P between 2 and 10 d. We also estimate the total number of binaries in the LMC with the same characteristics, and not only the eclipsing binaries, to be about 5000. This figure leads us to suggest that  (0.7 ± 0.4)  per cent of the main-sequence B-type stars in the LMC are found in binaries with periods shorter than 10 d. This frequency is substantially smaller than the fraction of binaries found by small Galactic radial-velocity surveys of B stars. On the other hand, the binary frequency found by Hubble Space Telescope ( HST ) photometric searches within the late main-sequence stars of 47 Tuc is only slightly higher and still consistent with the frequency we deduced for the B stars in the LMC.     

LU Vel (GJ 375): A M3.5Ve Flare and Double-Lined Spectroscopic Binary

High resolution echelle spectroscopic observations taken with the FEROS spectrograph at the 2.2 m telescope ESO confirm the binary nature of the flare M3.5V star LU Vel (GJ 375, RE J0958-462) previously reported by Christian and Mathioudakis (2002). Emission of similar intensity from both components is detected in the Balmer, Na i D₁&D₂, He i D₃, Ca ii H&K, and Ca ii IRT lines. We have determined precise radial velocities by cross correlation with radial velocity standard stars, which have allowed us to obtain for the first time the orbital solution of the system. The binary consists of two near-equal M3.5V components with an orbital period shorter than 2 days. We have analyzed the behaviour of the chromospheric activity indicators (variability and possible flares). In addition, we have determined its rotational velocity and kinematics.     

Single star scidar: atmospheric parameters profiling using the simulated annealing algorithm

In this paper, we present a new method to estimate, for each turbulent layer labelled i , the horizontal wind speed   v ( h _i )  , the standard deviation of the horizontal wind speed fluctuations  σ _v ( h_i )  and the integrated value of   C ² _n   over the thickness  Δ h_i   of the turbulent layer   C ² _n ( h_i )Δ h_i   , where   h_i   is the altitude of the turbulent layer. These parameters are extracted from single star scintillation spatiotemporal cross-correlation functions of atmospheric speckles obtained within the generalized mode. This method is based on the simulated annealing algorithm to find the optimal solution required to solve the problem. Astrophysics parameters for adaptive optics are also calculated using   C ² _n ( h_i )  and   v ( h_i )  values. The results of other techniques support this new method.     

Smeared crack approach: back to the original track

This paper briefly reviews the formulations used over the last 40 years for the solution of problems involving tensile cracking, with both the discrete and the smeared crack approaches. The paper focuses on the smeared approach, identifying as its main drawbacks the observed mesh‐size and mesh‐bias spurious dependence when the method is applied ‘straightly’. A simple isotropic local damage constitutive model is considered, and the (exponential) softening modulus is regularized according to the material fracture energy and the element size. The continuum and discrete mechanical problems corresponding to both the weak discontinuity (smeared cracks) and the strong discontinuity (discrete cracks) approaches are analysed and the question of propagation of the strain localization band (crack) is identified as the main difficulty to be overcome in the numerical procedure. A tracking technique is used to ensure stability of the solution, attaining the necessary convergence properties of the corresponding discrete finite element formulation. Numerical examples show that the formulation derived is stable and remarkably robust. As a consequence, the results obtained do not suffer from spurious mesh‐size or mesh‐bias dependence, comparing very favourably with those obtained with other fracture and continuum mechanics approaches. Copyright © 2006 John Wiley & Sons, Ltd.