As an illustration of the value of historical measures we present some examples of the dynamical evolution of multiple systems resembling the Orion Trapezium. We constructed models by combining carefully selected historical measures of the separations among components of young massive stellar systems with modern observations. By computing large numbers of fictitious systems resembling real trapezia we were able to simulate the dynamical evolution of such systems. Our results on the dynamical fate of the Orion Trapezium and of ten additional young clusters resembling the Orion Trapezium show extremely short dynamical lifetimes for these systems.

Preliminary results on the ongoing spectroscopic study of the six brightest Orion Trapezium stars is presented here. The main purpose of this work is to better understand the multiplicity and stability of each of these subsystems and the dynamical future of the group. So far the most interesting results reached are: 1) The orbit of the secondary star of the eclipsing Component A (V1016 Ori) is highly inclined with respect to the equatorial plane of its primary star. 2) The also eclipsing binary BM Ori (Trapezium Component B) does have a tertiary member with period about 3.5 years, as proposed by Vitrichenko & Klochkova (2004), and is the same as the companion recently found by the GRAVITY collaboration et al. (2018}. 3) Component D is indeed a spectroscopic and interferometric double star with a relatively high-mass companion ($q=M_2/M_1=0.5$) and period $52.90\pm0.05\,d$. 4) Component F, is a CP star (B7.5 p Si); its radial velocity, $23.2\pm4.2\,km\,s^{-1}$, is smaller than that of all other Trapezium members and, possibly, the evolutionary stage of the star is more advanced than that of members with similar mass. Consequently, Component F is probably not physically related to the Trapezium. Several evidences point to the extreme youth of this stellar group; its further study, most likely, will shed light on the formation processes of massive stars.

We analyse new optical spectroscopic, direct-image and X-ray observations of the recently discovered a high proper motion cataclysmic variable V1838 Aql. The data were obtained during its 2013 superoutburst and its subsequent quiescent state. An extended emission around the source was observed up to 30 days after the peak of the superoutburst, interpreted it as a bow--shock formed by a quasi-continuous outflow from the source in quiescence. The head of the bow--shock is coincident with the high--proper motion vector of the source ($v_{\perp}=123\pm5$ km s$^{-1}$) at a distance of $d=202\pm7$ pc. The object was detected as a weak X-ray source ($0.015\pm0.002$ counts s$^{-1}$) in the plateau of the superoutburst, and its flux lowered by two times in quiescence (0.007$\pm$0.002 counts s$^{-1}$). Spectroscopic observations in quiescence we confirmed the orbital period value $P_{\rm{orb}}=0.0545\pm 0.0026$ days, consistent with early-superhump estimates, and the following orbital parameters: $\gamma= -21\pm3$ km s$^{-1}$ and $K_1 = 53\pm3$ km s$^{-1}$. The white dwarf is revealed as the system approaches quiescence, which enables us to infer the effective temperature of the primary $T_{eff}=11,600\pm400$K. The donor temperature is estimated $\lesssim 2200$K and suggestive of a system approaching the period minimum. Doppler maps in quiescence show the presence of the hot spot in HeI line at the expected accretion disc-stream shock position and an unusual structure of the accretion disc in H$\alpha$.

We present an in-depth photometric study of the 2013 superoutburst of the recently discovered cataclysmic variable V1838 Aql and subsequent photometry near its quiescent state. A careful examination of the development of the superhumps is presented. Our best determination of the orbital period is Porb = 0.05698(9) days, based on the periodicity of early superhumps. Comparing the superhump periods at stages A and B with the early superhump value we derive a period excess of {\epsilon} = 0.024(2) and a mass ratio of q = 0.10(1). We suggest that V1838 Aql is approaching the orbital period minimum and thus has a low-mass star as a donor instead of a sub-stellar object.

We have selected archival observational data for several O and B trapezia in the Milky Way. For each of the main components of the trapezia we obtained transverse velocities from the historical separation data. With this information, and with the stellar masses of the main components, we studied the dynamical evolution of ensembles of multiple systems mimicking each one of the trapezia. For this purpose we conducted numerical $N-$body integrations using the best available values for the masses, the observed positions and {\bf transverse} velocities, randomly generated radial velocities, and random line-of-sight ($z$) positions for all components. Random perturbations were assigned to the observed quantities, compatible with the observational errors. A large fraction of the simulated systems ({\bf between 70 and 90 percent}) turned out to be unbound. The properties of the evolving systems are studied at different values of the evolution time. We find that the dynamical lifetimes of both the bound and unbound seems to be quite short, of less than 10~000 yr for the unbound systems, and of 10~000 to 20~000 yr for the bound systems. The end result of the simulations is usually a binary {\bf with semiaxes of a few hundred AU} sometimes a triple of hierarchical on non-hierarchical type. Non-hierarchical triples formed during the integrations (which are dynamically unstable) were found to have much longer lifetimes, of 250~000 to 500~000 yr. The frequency distributions of the major semi-axes and eccentricities of the resulting binaries are discussed and compared with observational properties of binary systems from the literature.

We present the second paper of a series of publications aiming at obtaining a better understanding regarding the nature of type Ia supernovae (SNIa) progenitors by studying a large sample of detached F, G and K main sequence stars in close orbits with white dwarf companions (i.e. WD+FGK binaries). We employ the LAMOST (Large Sky Area Multi-Object Fibre Spectroscopic Telescope) data release 4 spectroscopic data base together with GALEX (Galaxy Evolution Explorer) ultraviolet fluxes to identify 1,549 WD+FGK binary candidates (1,057 of which are new), thus doubling the number of known sources. We measure the radial velocities of 1,453 of these binaries from the available LAMOST spectra and/or from spectra obtained by us at a wide variety of different telescopes around the globe. The analysis of the radial velocity data allows us to identify 24 systems displaying more than 3sigma radial velocity variation that we classify as close binaries. We also discuss the fraction of close binaries among WD+FGK systems, which we find to be ~10 per cent, and demonstrate that high-resolution spectroscopy is required to efficiently identify double-degenerate SNIa progenitor candidates.

We discuss recent observational data on the transverse and radial velocities, as well as on the masses of the main components of the Orion Trapezium. Based on the most reliable values of these quantities we study the dynamical evolution of ensembles of multiple systems mimicking the Orion Trapezium. To this end we conduct numerical $N-$body integrations using the observed masses, planar positions and velocities, radial velocities, and random line-of-sight ($z$) positions for all components. We include perturbations in these quantities compatible with the observational errors. We find the dynamical lifetimes of such systems to be quite short, of the order of 10 to 50 thousand years. The end result of the simulations is usually a tight binary, or sometimes a hierarchical triple. The properties of the evolved systems are studied at different values of the crossing times. The frequency distributions of the major semiaxes and eccentricities of the resulting binaries are discussed and compared with observations.

We have obtained simultaneous photometric and spectroscopic observations of the cataclysmic variable 1RXS J064434.5+334451. We have calibrated the spectra for slit losses using the simultaneous photometry allowing to construct reliable Doppler images from H$\alpha$ and HeII 4686 emission lines. We have improved the ephemeris of the object based on new photometric eclipse timings, obtaining $HJD = 2453403.759533 + 0.26937446E$. Some eclipses present a clear internal structure which we attribute to a central HeII emission region surrounding the white dwarf, a finding supported by the Doppler tomography. This indicates that the system has a large inclination angle $i=78 \pm 2^{\circ}$. We have also analysed the radial velocity curve from the emission lines to measure its semi--amplitude, $K_1$, from H$\alpha$ and HeII 4686 and derive the masses of the components: $M_1=0.82\pm0.06$ M$_{\odot}$, $M_2=0.78\pm0.04$ M$_{\odot}$ and their separation $a=2.01\pm0.06$ $R_{\odot}$. The Doppler tomography and other observed features in this nova-like system strongly suggests that this is a SW Sex-type system.

Using the novel Diffracto-Astrometry technique, we analyze 44 Hubble Space Telescope Wide Field Planetary Camera 2 images of the Orion Trapezium (OT) taken over a span of 12 yr (1995-2007). We measure the relative positions of the six brighter OT components (A-F) and supplement these results with measurements of the relative separations and position angles taken from the literature, thus extending our analysis time base to ~200 yr. For every pair of components we find the relative rate of separation as well as the temporal rate of change of their position angles, which enable us to determine the relative kinematics of the system. Component E shows a velocity larger than the OT's escape velocity, thus confirming that it is escaping from the gravitational pull of this system.

We used extensive ground-based multisite and archival spectroscopy to derive observational constraints for a seismic modelling of the magnetic beta Cep star V2052 Ophiuchi. The line-profile variability is dominated by a radial mode (f_1=7.14846 d^{-1}) and by rotational modulation (P_rot=3.638833 d). Two non-radial low-amplitude modes (f_2=7.75603 d^{-1} and f_3=6.82308 d^{-1}) are also detected. The four periodicities that we found are the same as the ones discovered from a companion multisite photometric campaign (Handler et al. 2012) and known in the literature. Using the photometric constraints on the degrees l of the pulsation modes, we show that both f_2 and f_3 are prograde modes with (l,m)=(4,2) or (4,3). These results allowed us to deduce ranges for the mass (M \in [8.2,9.6] M_o) and central hydrogen abundance (X_c \in [0.25,0.32]) of V2052 Oph, to identify the radial orders n_1=1, n_2=-3 and n_3=-2, and to derive an equatorial rotation velocity v_eq \in [71,75] km s^{-1}. The model parameters are in full agreement with the effective temperature and surface gravity deduced from spectroscopy. Only models with no or mild core overshooting (alpha_ov \in [0,0.15] local pressure scale heights) can account for the observed properties. Such a low overshooting is opposite to our previous modelling results for the non-magnetic beta Cep star theta Oph having very similar parameters, except for a slower surface rotation rate. We discuss whether this result can be explained by the presence of a magnetic field in V2052 Oph that inhibits mixing in its interior.

We present a theoretical analysis of the measuring algorithm we use when applying the Diffracto-Astrometry technique to Hubble Space Telescope Wide Field Planetary Camera 2 (WFPC2) saturated stellar images. Theoretical Point Spread Functions (PSFs) were generated using the Tiny Tim software for the four CCDs in the WFPC2 and for some of the available filters. These images were then measured with our Diffracto-Astrometry measuring algorithm using only their diffraction pattern, and positions for the simulated PSFs on each generated CCD-frame were obtained. The measuring algorithm recovers the original positions reasonably well ($\pm 0.1\ \rm{to} \ \pm 0.4$ pixels). However, slight deviations from the original values are observed. These also vary with position over the entire surface of the CCD. We adjust the difference between the real and the measured position with a quadratic function of the coordinates. The transformation coefficients also present a slight correlation with the filter effective wavelength. Application of these transformation coefficients allows us to determine the position of a stellar image with a precision of a few hundredths of a pixel.

We have made a search for common proper motion (CPM) companions to the wide binaries in the solar vicinity. We found that the binary GJ 282AB has a very distant CPM companion (NLTT 18149) at a separation $s=1.09 \arcdeg$. Improved spectral types and radial velocities are obtained, and ages determined for the three components. The Hipparcos trigonometric parallaxes and the new radial velocities and ages turn out to be very similar for the three stars, and provide strong evidence that they form a physical system. At a projected separation of 55733AU from GJ 282AB, NLTT 18149 ranks among the widest physical companions known.

We report Hubble Space Telescope imaging, obtained 155 and 449 days after the 2006 outburst of the recurrent nova RS Ophiuchi, together with ground-based spectroscopic observations, obtained from the Observatorio Astron\'omico Nacional en San Pedro M\'artir, Baja California, M\'exico and at the Observatorio Astrof\'isico Guillermo Haro, at Cananea, Sonora, M\'exico. The observations at the first epoch were used as inputs to model the geometry and kinematic structure of the evolving RS Oph nebular remnant. We find that the modeled remnant comprises two distinct co-aligned bipolar components; a low-velocity, high-density innermost (hour glass) region and a more extended, high-velocity (dumbbell) structure. This overall structure is in agreement with that deduced from radio observations and optical interferometry at earlier epochs. We find that the asymmetry observed in the west lobe is an instrumental effect caused by the profile of the HST filter and hence demonstrate that this lobe is approaching the observer. We then conclude that the system has an inclination to the line of sight of 39$^{+1}_{-10}$ degrees. This is in agreement with the inclination of the binary orbit and lends support to the proposal that this morphology is due to the interaction of the outburst ejecta with either an accretion disk around the central white dwarf and/or a pre-existing red giant wind that is significantly denser in the equatorial regions of the binary than at the poles. The second epoch HST observation was also modeled. However, as no spectra were taken at this epoch, it is more difficult to constrain any model. Nevertheless, we demonstrate that between the two HST epochs the outer dumbbell structure seems to have expanded linearly.

High-dispersion time-resolved spectroscopy of the unique magnetic cataclysmic variable AE Aqr is presented. A radial velocity analysis of the absorption lines yields K_2 = 168.7+/- 1 km/s. Substantial deviations of the radial velocity curve from a sinusoid are interpreted in terms of intensity variations over the secondary star's surface. A complex rotational velocity curve as a function of orbital phase is detected which has a modulation frequency of twice the orbital frequency, leading to an estimate of the binary inclination angle that is close to 70^o. The minimum and maximum rotational velocities are used to indirectly derive a mass ratio of q= 0.6 and a radial velocity semi-amplitude of the white dwarf of K_1 = 101+/-3 km/s. We present an atmospheric temperature indicator, based on the absorption line ratio of Fe I and Cr I lines, whose variation indicates that the secondary star varies from K0 to K4 as a function of orbital phase. The ephemeris of the system has been revised, using more than one thousand radial velocity measurements, published over nearly five decades. From the derived radial velocity semi-amplitudes and the estimated inclination angle, we calculate that the masses of the stars are M_1 = 0.63+/-0.05M_sun; M_2 = 0.37+/-0.04 M_sun, and their separation is a = 2.33+/-0.02R_sun. Our analysis indicates the presence of a late-type star whose radius is larger, by a factor of nearly two, than the radius of a normal main sequence star of its mass. Finally we discuss the possibility that the measured variations in the rotational velocity, temperature, and spectral type of the secondary star as functions of orbital phase may, like the radial velocity variations, be attributable to regions of enhanced absorption on the star's surface.

High-resolution spectroscopy of U Gem was obtained during quiescence. We did not find a hot spot or gas stream around the outer boundaries of the accretion disk. Instead, we detected a strong narrow emission near the location of the secondary star. We measured the radial velocity curve from the wings of the double-peaked H$\alpha$ emission line, and obtained a semi-amplitude value that is in excellent agreement with the obtained from observations in the ultraviolet spectral region by Sion et al. (1998). We present also a new method to obtain K_2, which enhances the detection of absorption or emission features arising in the late-type companion. Our results are compared with published values derived from the near-infrared NaI line doublet. From a comparison of the TiO band with those of late type M stars, we find that a best fit is obtained for a M6V star, contributing 5 percent of the total light at that spectral region. Assuming that the radial velocity semi-amplitudes reflect accurately the motion of the binary components, then from our results: K_em = 107+/-2 km/s; K_abs = 310+/-5 km/s, and using the inclination angle given by Zhang & Robinson(1987); i = 69.7+/-0.7, the system parameters become: M_WD = 1.20+/-0.05 M_sun,; M_RD = 0.42+/-0.04 M_sun; and a = 1.55+/- 0.02 R_sun. Based on the separation of the double emission peaks, we calculate an outer disk radius of R_out/a ~0.61, close to the distance of the inner Lagrangian point L_1/a~0.63. Therefore we suggest that, at the time of observations, the accretion disk was filling the Roche-Lobe of the primary, and that the matter leaving the L_1 point was colliding with the disc directly, producing the hot spot at this location.

High-dispersion spectroscopy of EY Cyg obtained from data spanning twelve years show, for the first time, the radial velocity curves from both emission and absorption line systems, yielding semi-amplitudes K_{em}=24+/- 4 km s^-1 and K_{abs}=54+/- 2 km s^-1. The orbital period of this system is found to be 0.4593249(1)d. The masses of the stars, their mass ratio and their separation are found to be M_1 sin^3 i = 0.015+/-0.002 M_sun, M_2 sin^3 i = 0.007+/-0.002 M_sun, q = K_1/K_2 = M_2/M_1 = 0.44+/-0.02 and a sin i = 0.71+/-0.04 R_sun. We also found that the spectral type of the secondary star is around K0,consistent with an early determination by Kraft(1962). From the spectral type of the secondary star and simple comparisons with single main sequence stars, we conclude that the radius of the secondary star is about 30 per cent larger than a main sequence star of the same mass. We also present VRI CCD photometric observations, some of them simultaneous with the spectroscopic runs. The photometric data shows several light modulations, including a sinusoidal behaviour with twice the frequency of the orbital period, characteristic of the modulation coming from an elongated, irradiated secondary star. Low and high states during quiescence are also detected and discussed. From several constrains, we obtain tight limits for the inclination angle of the binary system between 13 and 15 degrees, with a best value of 14 degrees obtained from the sinusoidal light curve analysis. From the above results we derive masses M_1 = 1.10+/-0.09 M_sun, M_2 = 0.49+/-0.09 M_sun, and a binary separation a = 2.9+/- 0.1 R_sun.