We observed the afterglow of GRB 130427A with the RATIR instrument on the 1.5-m Harold L. Johnson telescope of the Observatorio Astron\'omico Nacional on Sierra San Pedro M\'artir. Our homogenous $griZYJH$ photometry extends from the night of burst to three years later. We fit a model for the afterglow. There is a significant positive residual which matches the behavior of SN 1998bw in the $griZ$ filters; we suggest that this is a photometric signature of the supernova SN 2013cq associated with the GRB. The peak absolute magnitude of the supernova is $M_r=-18.43\pm0.11$.

We report on the first open-use based Atacama Large Millimeter/submm Array (ALMA) 345-GHz observation for the late afterglow phase of GRB131030A. The ALMA observation constrained a deep limit at 17.1 d for the afterglow and host galaxy. We also identified a faint submillimeter source (ALMAJ2300-0522) near the GRB131030A position. The deep limit at 345 GHz and multifrequency observations obtained using {\it Swift} and RATIR yielded forward shock modeling with a two-dimensional relativistic hydrodynamic jet simulation and described X-ray excess in the afterglow. The excess was inconsistent with the synchrotron self-inverse Compton radiation from the forward shock. The host galaxy of GRB131030A and optical counterpart of ALMAJ2300-0522 were also identified in the SUBARU image. Based on the deep ALMA limit for the host galaxy, the 3-$\sigma$ upper limits of IR luminosity and the star formation rate (SFR) is estimated as $L_{IR}<1.11\times10^{11} L_{\odot}$ and SFR$<18.7$ ($M_{\odot}$~yr$^{-1}$), respectively. Although the separation angle from the burst location (3.5 arcsec) was rather large, ALMAJ2300-0522 may be one component of the GRB131030A host galaxy, according to previous host galaxy cases.

We present optical and near-infrared light curves and optical spectra of SN 2013dx, associated with the nearby (redshift 0.145) gamma-ray burst GRB 130702A. The prompt isotropic gamma-ray energy released from GRB 130702A is measured to be $E_{\gamma,iso}=6.4_{-1.0}^{+1.3}\times10^{50}$erg (1keV-10MeV in the rest frame), placing it intermediate between low-luminosity GRBs like GRB 980425/SN 1998bw and the broader cosmological population. We compare the observed $g'r'i'z'$ light curves of SN 2013dx to a SN 1998bw template, finding that SN 2013dx evolves ~20% faster (steeper rise time), with a comparable peak luminosity. Spectroscopically, SN 2013dx resembles other broad-lined Type Ic supernovae, both associated with (SN 2006aj and SN 1998bw) and lacking (SN 1997ef, SN 2007I, and SN 2010ah) gamma-ray emission, with photospheric velocities around peak of ~21,000km s$^{-1}$. We construct a quasi-bolometric ($g'r'i'z'yJ$) light curve for SN 2013dx, only the fifth GRB-associated SN with extensive NIR coverage and the third with a bolometric light curve extending beyond $\Delta t>40$d. Together with the measured photospheric velocity, we derive basic explosion parameters using simple analytic models. We infer a $^{56}$Ni mass of $M_{\mathrm{Ni}}=0.37\pm0.01$M$_{\odot}$, an ejecta mass of $M_{\mathrm{ej}}=3.1\pm0.1$M$_{\odot}$, and a kinetic energy of $E_{\mathrm{K}}=(8.2\pm0.43)\times10^{51}$ erg (statistical uncertainties only), consistent with previous GRB-associated SNe. When considering the ensemble population of GRB-associated SNe, we find no correlation between the mass of synthesized $^{56}$Ni and high-energy properties, despite clear predictions from numerical simulations that $M_{\mathrm{Ni}}$ should correlate with the degree of asymmetry. On the other hand, $M_{\mathrm{Ni}}$ clearly correlates with the kinetic energy of the supernova ejecta across a wide range of core-collapse events.

We present our extensive observational campaign on the Swift-discovered GRB141121A, al- most ten years after its launch. Our observations covers radio through X-rays, and extends for more than 30 days after discovery. The prompt phase of GRB 141121A lasted 1410 s and, at the derived redshift of z = 1.469, the isotropic energy is E{\gamma},iso = 8.0x10^52 erg. Due to the long prompt duration, GRB141121A falls into the recently discovered class of UL-GRBs. Peculiar features of this burst are a flat early-time optical light curve and a radio-to-X-ray rebrightening around 3 days after the burst. The latter is followed by a steep optical-to-X-ray decay and a much shallower radio fading. We analyze GRB 141121A in the context of the standard forward-reverse shock (FS,RS) scenario and we disentangle the FS and RS contributions. Finally, we comment on the puzzling early-time (t ~3 d) behavior of GRB 141121A, and suggest that its interpretation may require a two-component jet model. Overall, our analysis confirms that the class of UL-GRBs represents our best opportunity to firmly establish the prominent emission mechanisms in action during powerful GRB explosions, and future missions (like SVOM, XTiDE, or ISS-Lobster) will provide many more of such objects.

We report here the discovery by the Intermediate Palomar Transient Factory (iPTF) of iPTF14yb, a luminous ($M_{r}\approx-27.8$ mag), cosmological (redshift 1.9733), rapidly fading optical transient. We demonstrate, based on probabilistic arguments and a comparison with the broader population, that iPTF14yb is the optical afterglow of the long-duration gamma-ray burst GRB 140226A. This marks the first unambiguous discovery of a GRB afterglow prior to (and thus entirely independent of) an associated high-energy trigger. We estimate the rate of iPTF14yb-like sources (i.e., cosmologically distant relativistic explosions) based on iPTF observations, inferring an all-sky value of $\Re_{\mathrm{rel}}=610$ yr$^{-1}$ (68% confidence interval of 110-2000 yr$^{-1}$). Our derived rate is consistent (within the large uncertainty) with the all-sky rate of on-axis GRBs derived by the Swift satellite. Finally, we briefly discuss the implications of the nondetection to date of bona fide "orphan" afterglows (i.e., those lacking detectable high-energy emission) on GRB beaming and the degree of baryon loading in these relativistic jets.

We present optical and near-infrared (NIR) photometry of 28 gamma-ray bursts (GRBs) detected by the \textit{Swift} satellite and rapidly observed by the Reionization and Transients Infrared/Optical (RATIR) camera. We compare the optical flux at fiducial times of 5.5 and 11 hours after the high-energy trigger to that in the X-ray regime to quantify optical darkness. 46$\pm$9 per cent (13/28) of all bursts in our sample and 55$\pm$10 per cent (13/26) of long GRBs are optically dark, which is statistically consistently with previous studies. Fitting RATIR optical and NIR spectral energy distributions (SEDs) of 19 GRBs, most (6/7) optically dark GRBs either occur at high-redshift ($z>4.5$) or have a high dust content in their host galaxies ($A_{\rm V} > 0.3$). Performing K-S tests, we compare the RATIR sample to those previously presented in the literature, finding our distributions of redshift, optical darkness, host dust extinction and X-ray derived column density to be consistent. The one reported discrepancy is with host galaxy dust content in the BAT6 sample, which appears inconsistent with our sample and other previous literature. Comparing X-ray derived host galaxy hydrogen column densities to host galaxy dust extinction, we find that GRBs tend to occur in host galaxies with a higher metal-to-dust ratio than our own Galaxy, more akin to the Large and Small Magellanic Clouds. Finally, to mitigate time evolution of optical darkness, we measure $\beta_{\rm OX,rest}$ at a fixed rest frame time, $t_{\rm rest}=1.5$ hours and fixed rest frame energies in the X-ray and optical regimes. Choosing to evaluate optical flux at $\lambda_{\rm rest}=0.25~\mu$m, we remove high-redshift as a source of optical darkness, demonstrating that optical darkness must result from either high-redshift, dust content in the host galaxy along the GRB sight line, or a combination of the two.

We present results from the first attempts to derive various physical characteristics of the dusty Wolf-Rayet star WR 48a based on a multi-wavelength view of its observational properties. This is done on the basis of new optical and near-infrared spectral observations and on data from various archives in the optical, radio and X-rays. The optical spectrum of WR 48a is acceptably well represented by a sum of two spectra: of a WR star of the WC8 type and of a WR star of the WN8h type. The strength of the interstellar absorption features in the optical spectra of WR 48a and the near-by stars D2-3 and D2-7 (both members of the open cluster Danks 2) indicates that WR 48a is located at a distance of ~4 kpc from us. WR 48a is very likely a thermal radio source and for such a case and smooth (no clumps) wind its radio emission suggests a relatively high mass-loss rate of this dusty WR star (dM/dt = a few x 10^(-4) solar masses per year). Long timescale (years) variability of WR 48a is established in the optical, radio and X-rays. Colliding stellar winds likely play a very important role in the physics of this object. However, some LBV-like (luminous blue variable) activity could not be excluded as well.

We present optical and near-infrared (NIR) photometry for three gamma-ray burst supernovae (GRB-SNe): GRB 120729A, GRB 130215A / SN 2013ez and GRB 130831A / SN 2013fu. In the case of GRB 130215A / SN 2013ez, we also present optical spectroscopy at t-t0=16.1 d, which covers rest-frame 3000-6250 Angstroms. Based on Fe II (5169) and Si (II) (6355), our spectrum indicates an unusually low expansion velocity of 4000-6350 km/s, the lowest ever measured for a GRB-SN. Additionally, we determined the brightness and shape of each accompanying SN relative to a template supernova (SN 1998bw), which were used to estimate the amount of nickel produced via nucleosynthesis during each explosion. We find that our derived nickel masses are typical of other GRB-SNe, and greater than those of SNe Ibc that are not associated with GRBs. For GRB 130831A / SN 2013fu, we use our well-sampled R-band light curve (LC) to estimate the amount of ejecta mass and the kinetic energy of the SN, finding that these too are similar to other GRB-SNe. For GRB 130215A, we take advantage of contemporaneous optical/NIR observations to construct an optical/NIR bolometric LC of the afterglow. We fit the bolometric LC with the millisecond magnetar model of Zhang & Meszaros (2001), which considers dipole radiation as a source of energy injection to the forward shock powering the optical/NIR afterglow. Using this model we derive an initial spin period of P=12 ms and a magnetic field of B=1.1 x 10^15 G, which are commensurate with those found for proposed magnetar central engines of other long-duration GRBs.

We present a template fitting algorithm for determining photometric redshifts, $z_{\rm phot}$, of candidate high-redshift gamma-ray bursts (GRBs). Using afterglow photometry, obtained by the Reionization And Transients InfraRed (RATIR) camera, this algorithm accounts for the intrinsic GRB afterglow spectral energy distribution (SED), host dust extinction and the effect of neutral hydrogen (local and cosmological) along the line of sight. We present the results obtained by this algorithm and RATIR photometry of GRB 130606A, finding a range of best fit solutions $5.6 < z_{\rm phot} < 6.0$ for models of several host dust extinction laws (none, MW, LMC and SMC), consistent with spectroscopic measurements of the redshift of this GRB. Using simulated RATIR photometry, we find our algorithm provides precise measures of $z_{\rm phot}$ in the ranges $4 < z_{\rm phot} \lesssim 8$ and $9 < z_{\rm phot} < 10$ and can robustly determine when $z_{\rm phot}>4$. Further testing highlights the required caution in cases of highly dust extincted host galaxies. These tests also show that our algorithm does not erroneously find $z_{\rm phot} < 4$ when $z_{\rm sim}>4$, thereby minimizing false negatives and allowing us to rapidly identify all potential high-redshift events.

We present spatially- and velocity-resolved echelle spectroscopy for NGC 7009 obtained with the UVES spectrograph at the European Southern Observatory's Very Large Telescope. Our objective is to analyze the kinematics of emission lines excited by recombination and collisions with electrons to determine whether similarities or differences could be useful in elucidating the well-known abundance discrepancy derived from them. We construct position-velocity maps for recombination, fluorescence, charge transfer, and collisionally-excited lines. We find a plasma component emitting in the C II, N II, O II, and Ne II recombination lines whose kinematics are discrepant: They are incompatible with the ionization structure derived from all other evidence and the kinematics derived from all of these lines are unexpectedly very similar. We find direct evidence for a recombination contribution to [N II] 5755. Once taken into account, the electron temperatures from [N II], [O III], and [Ne III] agree at a given position and velocity. The electron densities derived from [O II] and [Ar IV] are consistent with direct imaging and the distribution of hydrogen emission. The kinematics of the C II, N II, O II, and Ne II lines does not coincide with the kinematics of the [O III] and [Ne III] forbidden emission, indicating that there is an additional plasma component to the recombination emission that arises from a different volume from that giving rise to the forbidden emission from the parent ions within NGC 7009. Thus, the chemical abundances derived from either type of line are correct only for the plasma component from which they arise. Apart from [N II] 5755, we find no anomaly with the forbidden lines usually used to determine chemical abundances in ionized nebulae, so the abundances derived from them should be reliable for the medium from which they arise.

Aims. We report new results for the massive evolved and main sequence members of the young galactic cluster DBS2003 179. We determine the physical parameters and investigate the high-mass stellar content of the cluster, as well as of its close vicinity. Methods. Our analysis is based on ISAAC/VLT moderate-resolution (R\approx4000) infrared spectroscopy of the brightest cluster members. We derive stellar parameters for sixteen of the stellar members, using full non-LTE modeling of the obtained spectra. Results. The cluster contains three late WN or WR/LBV stars (Obj 4, Obj 15, and Obj 20:MDM32) and at least 5 OIf and 5 OV stars. According to the Hertzsprung-Russell diagram for DBS2003 179, the WR stars show masses above 85Msun, the OIf stars are between 40 and 80Msun, and the main sequence O stars are >20Msun. There are indications of binarity for Obj 4 and Obj 11, and Obj 3 shows a variable spectrum. The cluster is surrounded by a continuous protostar formation region most probably triggered by DBS2003 179. Conclusions. We confirm that DBS2003 179 is young massive cluster (2.5 10^4Msun) very close to the Galactic center at the distance of 7.9+-0.8 kpc.

Two observational campaigns were carried out during the eclipses of EE Cep in 2003 and 2008/9 to verify whether the eclipsing body in the system is indeed a dark disk and to understand the observed changes in the depth and durations of the eclipses. Multicolour photometric data and spectroscopic observations at both low and high resolution were collected. We numerically modelled the variations in brightness and colour during the eclipses. We tested models with different disk structure. We considered the possibility of disk precession. The complete set of observational data collected during the last three eclipses are made available to the astronomical community. Two blue maxima in the colour indices were detected during these two eclipses, one before and one after the photometric minimum. The first (stronger) blue maximum is simultaneous with a "bump" that is very clear in all the UBVRI light curves. Variations in the spectral line profiles seem to be recurrent during each cycle. NaI lines always show at least three absorption components during the eclipse minimum and strong absorption is superimposed on the H_alpha emission. These observations confirm that the eclipsing object in EE Cep system is indeed a dark, dusty disk around a low luminosity object. The primary appears to be a rapidly rotating Be star that is strongly darkened at the equator and brightened at the poles. Some of the conclusions of this work require verification in future studies: (i) a complex, possibly multi-ring structure of the disk in EE Cep; (ii) our explanation of the "bump" observed during the last two eclipses in terms of the different times of obscuration of the hot polar regions of the Be star by the disk; and (iii) our suggested period of the disk precession (~11-12 P_orb) and predicted depth of about 2 mag the forthcoming eclipse in 2014.

We compare calculated intensities of lines of CII, NI, NII, OI and OII with a published deep spectroscopic survey of planetary nebula IC 418. Our calculations use a self-consistent nebular model and a synthetic spectrum of the central star atmosphere to take into account line excitation by continuum fluorescence and electron recombination. We found that the NII spectrum of s, p and most d states is excited by fluorescence due to the low excitation conditions of the nebula. Many CII and OII lines have significant amounts of excitation by fluorescence. Recombination excites all the lines from f and g states and most OII lines. In the neutral-ionized boundary the NI quartet and OI triplet dipole allowed lines are excited by fluorescence, while the quintet OI lines are excited by recombination. Electron excitation produces the forbidden optical lines of OI, and continuum fluorescence enhances the NI forbidden line intensities. Lines excited by fluorescence of light below the Lyman limit thus suggest a new diagnostic to explore the inner boundary of the photodissociation region of the nebula.

We predict intensities of lines of CII, NI, NII, OI and OII and compare them with a deep spectroscopic survey of IC 418 to test the effect of excitation of nebular emission lines by continuum fluorescence of starlight. Our calculations use a nebular model and a synthetic spectrum of its central star to take into account excitation of the lines by continuum fluorescence and recombination. The NII spectrum is mostly produced by fluorescence due to the low excitation conditions of the nebula, but many CII and OII lines have more excitation by fluorescence than recombination. In the neutral envelope, the NI permitted lines are excited by fluorescence, and almost all the OI lines are excited by recombination. Electron excitation produces the forbidden optical lines of OI, but continuum fluorescence excites most of the NI forbidden line intensities. Lines excited by fluorescence of light below the Lyman limit thus suggest a new diagnostic to explore the photodissociation region of a nebula.

Over the past 40 years, the massive LBV/WR system HD 5980 in the Small Magellanic Cloud has undergone a long-term S Doradus type variability cycle and two brief and violent eruptions in 1993 and 1994. In this paper we analyze a collection of UV and optical spectra obtained between 1979 and 2009 and perform CMFGEN model fits to spectra of 1994, 2000, 2002 and 2009. The results are as follows: a) The long term S Dor-type variability is associated with changes of the hydrostatic radius; b) The 1994 eruption involved changes in its bolometric luminosity and wind structure; c) the emission-line strength, the wind velocity and the continuum luminosity underwent correlated variations in the sense that a decreasing V$_\infty$ is associated with increasing emission line and continuum levels; and d) The spectrum of the third star in the system ({\it Star C}) is well-fit by a T$_{eff}$=32 K model atmosphere with SMC chemical abundances. For all epochs, the wind of the erupting star is optically thick at the sonic point and is thus driven mainly by the continuum opacity. We speculate that the wind switches between two stable regimes driven by the "hot" (during the eruption) and the "cool" (post-eruption) iron opacity bumps as defined by Lamers & Nugis (2002) and Gr\"afener and Hamann (2008), and thus the wind may undergo a bi-stability jump of a different nature from that which occurs in OB-stars.

We present an analysis of high resolution spectral observations of the symbiotic star StHa 190. A 30 days period has been derived from radial velocities of the G-type absorption lines and the HeII 4686A emission line. The main aim of this work was to look for explanation of the very wide absorption lines of the yellow giant. The very low mass function obtained from the absorption lines radial velocities suggests that the observed changes probably do not corespond to the orbital motion of this star.

We present a coherent stellar and nebular model reproducing the observations of the Planetary Nebula IC418. We want to test whether a stellar model obtained by fitting the stellar observations is able to satisfactory ionize the nebula and reproduce the nebular observations, which is by no mean evident. This allows us to determine all the physical parameters of both the star and the nebula, including the abundances and the distance. We used all the observational material available (FUSE, IUE, STIS and optical spectra) to constrain the stellar atmosphere model performed using the CMFGEN code. The photoionization model is done with Cloudy_3D, and is based on CTIO, Lick, SPM, IUE and ISO spectra as well as HST images. More than 140 nebular emission lines are compared to the observed intensities. We reproduce all the observations for the star and the nebula. The 3D morphology of the gas distribution is determined. The effective temperature of the star is 36.7kK. Its luminosity is 7700 solar luminosity. We describe an original method to determine the distance of the nebula using evolutionary tracks. No clumping factor is need to reproduce the age-luminosity relation. The distance of 1.25 kpc is found in very good agreement with recent determination using parallax method. The chemical composition of both the star and the nebula are determined. Both are Carbon-rich. The nebula presents evidence of depletion of elements Mg, Si, S, Cl (0.5 dex lower than solar) and Fe (2.9 dex lower than solar). This is the first self-consistent stellar and nebular model for a Planetary Nebula that reproduces all the available observations ranging from IR to UV, showing that the combined approach for the modeling process leads to more restrictive constraints and, in principle, more trustworthy results.

We are proposing to conduct a multicolor, synoptic infrared (IR) imaging survey of the Northern sky with a new, dedicated 6.5-meter telescope at San Pedro M\'artir (SPM) Observatory. This initiative is being developed in partnership with astronomy institutions in Mexico and the University of California. The 4-year, dedicated survey, planned to begin in 2017, will reach more than 100 times deeper than 2MASS. The Synoptic All-Sky Infrared (SASIR) Survey will reveal the missing sample of faint red dwarf stars in the local solar neighborhood, and the unprecedented sensitivity over such a wide field will result in the discovery of thousands of z ~ 7 quasars (and reaching to z > 10), allowing detailed study (in concert with JWST and Giant Segmented Mirror Telescopes) of the timing and the origin(s) of reionization. As a time-domain survey, SASIR will reveal the dynamic infrared universe, opening new phase space for discovery. Synoptic observations of over 10^6 supernovae and variable stars will provide better distance measures than optical studies alone. SASIR also provides significant synergy with other major Astro2010 facilities, improving the overall scientific return of community investments. Compared to optical-only measurements, IR colors vastly improve photometric redshifts to z ~ 4, enhancing dark energy and dark matter surveys based on weak lensing and baryon oscillations. The wide field and ToO capabilities will enable a connection of the gravitational wave and neutrino universe - with events otherwise poorly localized on the sky - to transient electromagnetic phenomena.

The next decade of survey astronomy has the potential to transform our knowledge of variable stars. Stellar variability underpins our knowledge of the cosmological distance ladder, and provides direct tests of stellar formation and evolution theory. Variable stars can also be used to probe the fundamental physics of gravity and degenerate material in ways that are otherwise impossible in the laboratory. The computational and engineering advances of the past decade have made large-scale, time-domain surveys an immediate reality. Some surveys proposed for the next decade promise to gather more data than in the prior cumulative history of astronomy. The actual implementation of these surveys will have broad implications for the types of science that will be enabled. We examine the design considerations for an optimal time-domain photometric survey dedicated to variable star science, including: observing cadence, wavelength coverage, photometric and astrometric accuracy, single-epoch and cumulative depth, overall sky coverage, and data access by the broader astronomical community. The best surveys must combine aspects from each of these considerations to fully realize the potential for the next decade of time-domain science.

Recent near- and mid-infrared surveys have brought evidence that the Milky Way continues to form massive clusters. We carry out a program to determine the basic physical properties of the new massive cluster candidate [DBS2003]179. Medium-resolution K-band spectra and deep near-infrared images of [DBS2003]179 were used to derive the spectral types of eight member stars, and to estimate the distance and reddening to the cluster. Seven of ten stars with spectra show emission lines. Comparison with template spectra indicated that they are early O-type stars. The mean radial velocity of the cluster is Vrad=-77+-6 km/s. Knowing the spectral types of the members and the color excesses, we determined extinction Av~16.6 and distance modulus (m-M)0~14.5 mag (D~7.9 kpc). The presence of early O-stars and a lack of red supergiants suggests a cluster age of 2-5Myr. The total cluster mass is approximated to 0.7x10^4 Msun and it is not yet dynamically relaxed. The candidate [DBS2003]179 further increases the family of the massive young clusters in the Galaxy, although it appears less massive than the prototypical starburst clusters.

We analyze the chemical composition of the central star of the planetary nebula NGC 6543 based upon a detailed NLTE model of its stellar wind. The logarithmic abundances by number are H=12.00, He=11.00, C=9.03, N=8.36, O=9.02, Si=8.19, P=5.53, S=7.57 and Fe=7.24. Compared with the solar abundances, most of the elements have solar composition with respect to hydrogen except C which is overabundant by 0.28 dex and Fe which is depleted by $\sim 0.2$ dex. Contrary to most previous work, we find that the star is not H-poor and has a normal He composition. These abundances are compared with those found in the diffuse X-ray plasma and the nebular gas. Compared to the plasma emitting in diffuse X-rays, the stellar wind is much less depleted in iron. Since the iron depletions in the nebular gas and X-ray plasma are similar, we conclude that the plasma emitting diffuse X-rays is derived from the nebular gas rather than the stellar wind. Excellent agreement is obtained between the abundances in the stellar wind and the nebular recombination line abundances for He, C, and O relative to H. On the other hand, the derived stellar N abundance is smaller than the nebular N abundance derived from recombination lines and agrees with the abundance found from collisionally-excited lines. The mean temperature variation determined by five different methods indicates that the difference in the nebular abundances between the recombination lines and collisionally excited lines can be explained as due to the temperature variations in a chemically homogeneous medium.

We calculate the hydrogen and helium ionization in B[e] envelopes and explore their dependence on mass-loss and effective temperature. We also present simulated observations of the Halpha emission line and the C IV 1550 doublet, and study their behavior. This paper reports our first results in an ongoing study of B[e] supergiants, and provides a glimpse on the ionization of the most important elements in self-consistent numerical simulations. Our newly developed 2D stellar atmosphere code, ASTAROTH, was used for the numerical simulations. The code self-consistently solves for the continuum radiation, non-LTE level populations, and electron temperature in axi-symmetric stellar envelopes. Observed profiles were calculated by an auxiliary program developed separately from ASTAROTH. In all but one of our models, H remained fully ionized. Due to ionizations from excited states it is much more difficult to get a H neutral disk than indicated by previous analytical calculations. Near the poles, the ionization is high in all models, while helium recombined in the equatorial regions for all but our lowest mass-loss rate. Although the model parameters were not adjusted to provide fits to any particular star, the theoretical profiles show some features seen in the profiles of R126. These include the partially resolved double peaked profile of Halpha, and the weak emission associated with the UV C IV resonance line.

We present the results of optical wavelength observations of the unusual SMC eclipsing binary system HD 5980 obtained in 1999 and 2004--2005. Radial velocity curves for the erupting LBV/WR object (star A) and its close WR-like companion (star B) are obtained by deblending the variable emission-line profiles of N IV and N V lines under the simplistic assumption that these lines originate primarily in the winds of star A and star B. The derived masses M_A=58--79 Mo and M_B=51--67 Mo, are more consistent with the stars' location near the top of the HRD than previous estimates. The presence of a wind-wind interaction region is inferred from the orbital phase-dependent behavior of He I P Cygni absorption components. The emission-line intensities continued with the declining trend previously seen in UV spectra. The behavior of the photospheric absorption lines is consistent with the results of Schweickhardt (2002) who concludes that the third object in the combined spectrum, star C, is also a binary system with P(starC)~96.5 days, e=0.83. The data used in this paper will be made publicly available for further analysis.

Young massive clusters are perfect astrophysical laboratories for study of massive stars. Clusters with Wolf-Rayet (WR) stars are of special importance, since this enables us to study a coeval WR population at a uniform metallicity and known age. GLIMPSE30 (G30) is one of them. The cluster is situated near the Galactic plane (l=298.756deg, b=-0.408deg) and we aimed to determine its physical parameters and to investigate its high-mass stellar content and especially WR stars. Our analysis is based on SOFI/NTT JsHKs imaging and low resolution (R~2000) spectroscopy of the brightest cluster members in the K atmospheric window. For the age determination we applied isochrone fits for MS and Pre-MS stars. We derived stellar parameters of the WR stars candidates using a full nonLTE modeling of the observed spectra. Using a variety of techniques we found that G30 is very young cluster, with age t~4Myr. The cluster is located in Carina spiral arm, it is deeply embedded in dust and suffers reddening of Av~10.5+-1.1mag. The distance to the object is d=7.2+-0.9kpc. The mass of the cluster members down to 2.35Msol is ~1600Msol. Cluster's MF for the mass range of 5.6 to 31.6Msol shows a slope of Gamma=-1.01+-0.03. The total mass of the cluster obtained by this MF down to 1Msol is about 3x10^3Msol. The spectral analysis and the models allow us to conclude that in G30 are at least one Ofpe/WN and two WR stars. The WR stars are of WN6-7 hydrogen rich type with progenitor masses more than 60Msol. G30 is a new member of the exquisite family of young Galactic clusters, hosting WR stars. It is a factor of two to three less massive than some of the youngest super-massive star clusters like Arches, Quintuplet and Central cluster and is their smaller analog.

We present extensive optical spectroscopy of the early-type magnetic star HD 191612 (O6.5f?pe-O8fp). The Balmer and HeI lines show strongly variable emission which is highly reproducible on a well-determined 538-d period. Metal lines and HeII absorptions (including many selective emission lines but excluding He II 4686A emission) are essentially constant in line strength, but are variable in velocity, establishing a double-lined binary orbit with P(orb) = 1542d, e=0.45. We conduct a model-atmosphere analysis of the primary, and find that the system is consistent with a O8: giant with a B1: main-sequence secondary. Since the periodic 538-d changes are unrelated to orbital motion, rotational modulation of a magnetically constrained plasma is strongly favoured as the most likely underlying `clock'. An upper limit on the equatorial rotation is consistent with this hypothesis, but is too weak to provide a strong constraint.

The X-ray source 2S0114+650=LSI+65 010 is a binary system containing a B-type primary and a low mass companion believed to be a neutron star. The system has three reported periodicities: the orbital period, P{orb}~11.6 d, X-ray flaring with P{flare}~2.7 hr, and a "superorbital" X-ray periodicity P{super}~30.7 d. The objective of this paper is to show that the puzzling periodicities in the system may be explained in the context of scenarios in which tidal interactions drive oscillations in the B-supergiant star. We calculate the solution of the equations of motion for one layer of small surface elements distributed along the equator of the star, as they respond to the forces due to gas pressure, centrifugal, coriolis, viscous forces, and the gravitational forces of both stars. This calculation provides variability timescales that can be compared with the observations. In addition, we use observational data obtained at the Observatorio Astron\'omico Nacional en San Pedro M\'artir (OAN/SPM) between 1993-2004 to determine which of the periodicities may be present in the optical region. We suggest that the tidal oscillations lead to a structured stellar wind which, when fed to the neutron star, produces the X-ray modulations. The connection between the stellar oscillations and the modulation of the mass ejection may lie in the shear energy dissipation generated by the tangential motions that are produced by the tidal interaction, particularly in the tidal bulge region. The tidal oscillation scenario weakens the case for 2S0114+650 containing a magnetar descendent.

We have searched for the emission from \fex and \fexiv that is expected from the gas emitting in diffuse X-rays in \bd, NGC 6543, NGC 7009, and NGC 7027. Neither line was detected in any object. Models that fit the X-ray spectra of these objects indicate that the \fex emission should be below our detection thresholds, but the predicted \fexiv emission exceeds our observed upper limits (one sigma) by factors of at least 3.5 to 12. The best explanation for the absence of \fexiv is that the X-ray plasma is depleted in iron. In principle, this result provides a clear chemical signature that may be used to determine the origin of the X-ray gas in either the nebular gas or the stellar wind. At present, though various lines of evidence appear to favour a nebular origin, the lack of atmospheric and nebular iron abundances in the objects studied here precludes a definitive conclusion.

We present a new radiative transfer code for axi-symmetric stellar atmospheres and compare test results against 1D and 2D models with and without velocity fields. The code uses the short characteristic method with modifications to handle axi-symmetric and non-monotonic 3D wind velocities, and allows for distributed calculations. The formal solution along a characteristic is evaluated with a resolution that is proportional to the velocity gradient along the characteristic. This allows us to accurately map the variation of the opacities and emissivities as a function of frequency and spatial coordinates, but avoids unnecessary work in low velocity regions. We represent a characteristic with an impact-parameter vector p (a vector that is normal to the plane containing the characteristic and the origin) rather than the traditional unit vector in the direction of the ray. The code calculates the incoming intensities for the characteristics by a single latitudinal interpolation without any further interpolation in the radiation angles. Using this representation also provides a venue for distributed calculations since the radiative transfer can be done independently for each p.

Terminal wind velocity and mass loss rate are the most fundamental parameters of stellar winds. Unfortunately, their determination calls for high resolution spectroscopy in a range of wavelengths spanning from the ultraviolet to the infrared. For weak and/or distant objects, this becomes unfeasible. Nevertheless, it is possible to obtain an accurate estimate of these parameters through a simplified study of the formation processes of resonant lines which show P Cyg profiles. In this case, the line profile is a complex function of 6 parameters, with the radiative transport treatment becoming relatively inexpensive. However, preforming a 6 dimensional parameter fit raises a number of problems if one seeks an objective and automatic procedure to yield the optimal values, from which wind velocities and mass loss rates can be estimated. Useing a Likelihood function to construct a well defined statistical estimator of the goodness of fit which corresponds to a given model, we turn to a genetic algorithm through which we find the global maximum of the 6 dimensional Likelihood hyper-surface. We here present the implementation of the method, its successful testing with synthetic line profiles, where the answer is known in advance, together with first results of its application to real data

We discuss work toward developing a 2.5D non-LTE radiative transfer code. Our code uses the short characteristic method with modifications to handle realistic 3D wind velocities. We also designed this code for parallel computing facilities by representing the characteristics with an impact parameter vector p. This makes the interpolation in the radiation angles unnecessary and allows for an independent calculation for characteristics represented by different p vectors. The effects of the velocity field are allowed for by increasing, as needed, the number of grid points along a short characteristic. This allows us to accurately map the variation of the opacities and emissivities as a function of frequency and spatial coordinates. In the future we plan to use this transfer code with a 2D non-LTE stellar atmosphere program to self-consistently solve for level populations, the radiation field and temperature structure for stars with winds and without spherical symmetry.

We present a catalogue of OB associations in IC 1613. Using an automatic and objective method (Battinelli's 1991 technique) 60 objects were found. The size distribution reveals a significant peak at about 60 parsecs if a distance modulus of 24.27 mag is assumed. Spatial distributions of the detected associations and H II regions are strongly correlated.

We present high resolution spectroscopic observations of the massive X-ray binary system LS I+65 010=2S 0114+650 in the optical wavelength region. A correlation between equivalent width and radial velocity of photospheric absorption lines is found. The systemic velocity, inferred from the weaker lines is v_helio = -31 +/- 5 km s^-1, which, if attributed solely to the Galactic rotation curve, implies that LS I+65 010 lies within 3 kpc from the Sun. The ISM Na I D lines display 2 resolved high velocity components at v_helio = -70, -48 km s^-1, possibly associated with gas surrounding the binary system, in addition to the -24, -8 km s^-1 ISM features due to the Orion and the Perseus arm regions. Strong photospheric line profile variability is present on a night to night timescale, with He I 5875 A displaying an additional blue-shifted absorption in some of the spectra, similar to what is observed in the optical counterpart of Vela X-1. A connection between the extended blue wing and X-ray maximum is suggested. Short timescale variations in line profiles are detected on only two nights, but the evidence that these variations occur on the 2.78 hour X-ray flaring period is marginal.

We present Fabry-Perot and multi-object spectroscopy of the galaxies in Hickson compact group 31 (HCG31). Based upon our H$\alpha$ data cubes, galaxies A and C are a single entity, showing no discontinuity in their kinematics. Kinematically, galaxy E is probably a component of the A+C complex; otherwise it is a recently detached fragment. Galaxy F appears, both kinematically and chemically, to have formed from material tidally removed from the A+C complex. Galaxies B and G are kinematically distinct from this complex. Galaxy Q also has a radial velocity compatible with group membership. Galaxies A, B, C, and F have nearly identical oxygen abundances, despite spanning a luminosity range of 5$ $mag. Galaxy B's oxygen abundance is normal for its luminosity, while galaxy F's abundance is that expected given its origin as a tidal fragment of the A+C complex. The oxygen abundances in galaxies A and C are also understandable if the A+C complex is a late-type spiral suffering strong gas inflow and star formation as a result of a tidal interaction. Given the kinematics of both the galaxies and the \ion{H}{i} gas, the oxygen abundances, and the position of galaxy G, we propose that an interaction of galaxy G with the A+C complex, rather than a merger of galaxies A and C, is a more complete explanation for the tidal features and other properties of HCG31. In this case, the A+C complex need not be a merger in progress, though this is not ruled out.

A detailed kinematic analysis of ionized gas in the nearby irregular galaxy NGC 4449 is presented. Observations are conducted in the spectral lines of Halpha and [SII]. Our scanning Fabry--Perot interferometric observations are presented from both a global as well as a local perspective. We have analysed the global velocity field, the spatially extended diffuse gaseous component (DIG), the HII region populations, and, furthermore, have determined the rotation curve based on the heliocentric radial velocities of the global Halpha spatial distribution. Our results for NGC 4449 show that the optical velocity field has a decreasing value in radial velocity along the optical bar from NE to SW. The DIG component that permeates the entire galaxy was analysed (up to a limiting surface brightness of ~ 3.165x10^-5 ergs cm$^-2 s^-1 steradian^-1). We find that the diffuse gas component presents peculiar kinematical features such as abrupt velocity gradients and highly supersonic velocity dispersions (sigma~4 times the values of the nearest HII regions) but that its kinematical and dynamical influence is important on both global and local scales. The optical rotation curve of this nearby irregular shows that the NE sector rotates like a solid body (V_rot~40 km s^-1 at R=2 kpc). For the SW side, our results are not conclusive; the behavior of the gas at those locations is chaotic. We conclude that the origin of such complex kinematics and dynamics is undoubtedly related to the aftermath of an interaction experienced by this galaxy in the past (abridged).

King, Modjaz, & Li (1999) discovered Nova 1999 in IC1613 at Lick Observatory. Both Fugazza et al. (2000) and Borissova et al. (2000) questioned this classification, because they were able to detect the star on images obtained in previous years. In infrared frames taken on Oct. 15, 1998, the Nova 1999 has (J-K) = 1.14 and K = 14.69 mag. Our light curve study, based primarily on 92 unfiltered Lick images, suggested that the object could be a Mira-type variable with a period of 640.7 days. This period is very close to that obtained by Fugazza et al. (2000) - 631 days. The star is overluminous with respect to the period-luminosity (PL) relation derived by Feast et al. (1989) for Mira variables in the LMC. At longer periods (P > 400 ~ days), many LMC Miras show such behavior and the PL relation appears to break down. It is possible that the situation in IC1613 is similar. An optical spectrum obtained with the Keck-II telescope shows features typical of M3Ie or M3IIIe stars. We conclude that the star is a normal long-period M-type Mira variable, the first such star confirmed in IC1613.

We present a detailed kinematical analysis of two selected nebulae in the Local Group irregular galaxy IC 1613. The nebulae are: S8, the only known supernova remnant in this galaxy, and S3, a Wolf-Rayet nebula associated with the only WO star in this galaxy. For S8, we have obtained and analyzed its radial velocity field, where we found complex profiles which can be fitted by several velocity components. These profiles also show the presence of high velocity, low density gas. From this, we have obtained the expansion velocity, estimated the preshock density and calculated the basic kinematical parameters of this SNR. We suggest that in S8 we are seing a SNR partially hidden by dust. This suggestion comes from the fact that the SNR is located between two superbubbles where a ridge of obscured material unveils the existence of dust. Moreover, we show that this hypothesis prevails when energetic arguments are taken into account. In the case of S3, this nebula shows bipolar structure. By means of its kinematics, we have analyzed its two lobes, the ``waist'', as well as its relation with the nearest superbubbles. For the first time we are able to see closed the NW lobe, showing a clover leaf shape. This fact allows a better quantitative knowledge of the nebula as a whole. Furthermore, we found evidence of an expansion motion in the NW lobe. In the light of our results, we can express that these nebulae are the product of very massive stellar evolution. It is surprising the influence these stars still have in shaping their surrounding gas, and on the energy liberation towards the interstellar medium of this galaxy.

We present optical spectroscopy of H {\sc ii} regions in the Local Group galaxy IC 10 and $UBVR$ photometry of foreground stars in three fields towards this galaxy. From these data, we find that the foreground reddening due to the Milky Way is $E(B-V)=0.77\pm 0.07 \mathrm{mag}$. We find that IC 10 contains considerable internal dust, which qualitatively explains the variety of reddening values found by studies of its different stellar populations. Based upon our foreground reddening, IC 10 has intrinsic photometric properties like those of a blue compact dwarf galaxy, and not those of a dwarf irregular. This result is consistent with much evidence that IC 10 is in the throes of a starburst that began at least $10 \mathrm{Myr}$ ago. We also report the discovery of a new WR star in the H {\sc ii} region HL111c.

We present the light curve of Luminous Blue Variable candidate star GR290 (Romano's star) in M33. The photographic photometry was made in photographic plates taken in B band of the M33 galaxy and cover an eight year period, 1982 - 1990. Twenty five plates, separated in seven groups, have been used. CCD B magnitude of the star is also presented. The analysis of our data together with the Romano's magnitudes (1978) shows "normal" eruptions with amplitude of more than 1 mag and timescale of about 20 years and smaller oscillations with amplitude 0.5 mag and a period of about 320 days. This is a typical photometrical behavior for LBVs.

We present near infrared JHK photometry for the central area of the irregular dwarf galaxy IC10. Comparing the red supergiants of IC10 and IC1613 we have determined the reddening of E(B-V)=1.05+-0.10 and the dereddened distance modulus (m-M)_0=23.86+-0.12 mag of Population I stars in IC10. Detection of six Brgamma emission structures clearly outline two star forming regions in the central area of IC10. There is no H_2 emissions down to 4.65 $10^{-16}$ $\rm erg s^{-1} cm^{-2} arcsec^{-2}$ indicating that the dominant excitation mechanism of the molecular gas is from UV radiation from hot young stars. From the comparison of Brgamma and Halpha fluxes, we derive average extinction toward th e star forming regions E(B-V)_ Halpha= 1.8+-0.2. The summarized SFR derived from six Halpha and Brgamma emission structures in our field of view is 0.6 $M_\odot {\rm year^{-1}}$.

We present Halpha and [SII] observations for the Local Group irregular galaxy IC 1613 using the PUMA scanning Fabry-Perot interferometer. Our goal is to analyze the kinematics of the ionized gas in the complex sample of superbubbles located in the whole extension of our field (10 arcmin), which includes most of the optical emission of this galaxy, and to study the inter-relationship between young stellar associations and nebulae based on a previous study that we have made on the stellar associations of the central region of this galaxy. The ionized gas in this galaxy is distributed in classical HII regions and in a series of superbubbles (also called giant shells) covering a large fraction of the optical extent of the galaxy. We present a catalog of kinematical properties of both the HII regions of this galaxy and the superbubbles. We have also compared the kinematics of the ionized gas in HII regions to search for possible dynamic differences between neutral and ionized gas.

We present a moderately-deep JK photometry for three selected areas of the dwarf irregular galaxy IC1613. The color-magnitude diagrams contain a mixture of red supergiants, asymptotic giant branch stars and the brightest red giant stars. The red supergiants are massive (20 - 25 M_\odot) and young -with ages between 8 and 25 Myr. The most important result is the evidence of the decreasing density of the intermediate age AGB population in the vicinity of the HII regions in the galaxy. We also find age differences between AGB stars in the main body of the galaxy and those near the HII regions in the North-East. The former span a range in ages between 1 and 10 Gyr, while the latter are younger than 1 Gyr. Using the period-luminosity relation derived by Madore & Freedman (1991) and JK magnitudes of the Cepheid variable V20, we calculated (m-M)_K = 24.37+-0.2. The recently discovered Nova (King et al. 1999) was identified in Field III. Its presence of our images and its brightness questioned its classification as a nova.

In recent years, much studies have focused on determining the origin of the large-scale line-profile and/or photometric patterns of variability displayed by some apparently single Wolf-Rayet stars, with the existence of an unseen (collapsed?) companion or of spatially extended wind structures as potential candidates. We present observations of WR 1 which highlight the unusual character of the variations in this object. Our narrowband photometric observations reveal a gradual increase of the stellar continuum flux amounting to Delta v = 0.09 mag followed by a decline on about the same timescale (3-4 days). Only marginal evidence for variability is found during the 11 following nights. Strong, daily line-profile variations are also observed but they cannot be easily linked to the photometric variations. Similarly to the continuum flux variations, coherent time-dependent changes are observed in 1996 in the centroid, equivalent width, and skewness of He II 4686. Despite the generally coherent nature of the variations, we do not find evidence in our data for the periods claimed in previous studies. While the issue of a cyclical pattern of variability in WR 1 is still controversial, it is clear that this object might constitute in the future a cornerstone for our understanding of the mechanisms leading to the formation of largely anisotropic outflows in Wolf-Rayet stars.