The Committee on the Status of Women in Astronomy (CSWA) is calling on federal science funding agencies, in their role as the largest sources of funding for astronomy in the United States, to take actions that will advance career development and improve workplace conditions for women and minorities in astronomy. Funding agencies can and should lead the charge to produce excellent diversity and inclusion outcomes in our field by the 2030 Astrophysics Decadal Survey. Anecdotal and quantitative evidence, gathered both by the CSWA and other groups, shows that many structural barriers to success remain in our community. We acknowledge the success of programs like NSF-INCLUDES and NSF-ADVANCE and endorse their continued work. We also recommend further action to remove barriers to success for women and minority astronomers. Key recommendations are: - Federal agencies should push academic institutions to reconsider their support systems for scientists by requiring a high standard of pay and benefits. - Federal agencies should fund research and programs that explore and implement strategies for improving the experiences of scientists. - Federal agencies should require Principal Investigators (PIs) to provide plans for creating inclusive work environments and mentoring young scientists.

The Committee on the Status of Women in Astronomy (CSWA) is calling on federal science funding agencies, in their role as the largest sources of funding for astronomy in the United States, to take actions that will end harassment, particularly sexual harassment, in astronomical workplaces. Funding agencies can and should lead the charge to end harassment in astronomy by the 2030 Astrophysics Decadal Survey. Anecdotal and quantitative evidence, gathered both by the CSWA and other groups, shows that harassment is prevalent and damaging for women and minority astronomers and those in related fields. Actions recommended herein will increase the rate of reporting of harassment to agencies and improve their ability to investigate and take action against harassers. We also recommend that agencies participate in harassment prevention by creating and implementing the best anti-harassment education possible. Key recommendations are: - Federal agencies should improve their ethics policies by making harassment a form of scientific misconduct. - Federal agencies should mandate that institutions report to them when a funded Principal Investigator (PI) or co-Principal Investigator (co-PI) is found to be a perpetrator of harassment. - Federal funding agencies should provide online guides to help scientists identify harassment and connect them to the right resources for making confidential or official reports. - Federal agencies should create and ensure the implementation of anti-harassment trainings by making them a requirement of receiving grant funding.

Organizations that support science (astronomy) such as federal agencies, research centers, observatories, academic institutions, societies, etc. employ advisory committees and boards as a mechanism for reviewing their activities and giving advice on practices, policies and future directions. As with any scientific endeavor, there is concern over complementing these committees with enough members who have as broad a range of expertise and understanding as possible, so that bias is mitigated. However, for a number of reasons (logistical, practical, financial, etc.), committees can also not be infinitely large and thus trade-offs must be made. It is often recognized that conflicts of interest must be acknowledged within these committees, but what is not often recognized it the potential for unmitigated biases and "group think" that can be introduced as part of these committees. In this white paper, we recommend that advisory committees that collect community input, (e.g., the Decadal Survey review committee), also collect, compile and review input demographic data before finalizing reports, (e.g., the final 2020 Decadal Survey Report). A summary of these data should be released alongside the final survey report. This information would enable the committee to understand potential "blind spots" and biases of the data collection phase and inform future data collections of any barriers that affect the omission of perspectives from various demographics.

The late-type Be star $\beta$ CMi is remarkably stable compared to other Be stars that have been studied. This has led to a realistic model of the outflowing Be disk by Klement et al. These results showed that the disk is likely truncated at a finite radius from the star, which Klement et al.~suggest is evidence for an unseen binary companion in orbit. Here we report on an analysis of the Ritter Observatory spectroscopic archive of $\beta$ CMi to search for evidence of the elusive companion. We detect periodic Doppler shifts in the wings of the H$\alpha$ line with a period of 170 d and an amplitude of 2.25 km s$^{-1}$, consistent with a low-mass binary companion ($M\approx 0.42 M_\odot$). We then compared the small changes in the violet-to-red peak height changes ($V/R$) with the orbital motion. We find weak evidence that it does follow the orbital motion, as suggested by recent Be binary models by Panoglou et al. Our results, which are similar to those for several other Be stars, suggest that $\beta$ CMi may be a product of binary evolution where Roche lobe overflow has spun up the current Be star, likely leaving a hot subdwarf or white dwarf in orbit around the star. Unfortunately, no direct sign of this companion star is found in the very limited archive of {\it International Ultraviolet Explorer} spectra.

We analyze a 162 ks HETG Chandra observation of the O7.5 III(n)((f)) star xi Per, together with contemporaneous H alpha observations. The X-ray spectrum of this star is similar to other single O stars, and not pathological in any way. Its UV wind lines are known to display cyclical time variability, with a period of 2.086 days, which is thought to be associated with co-rotating interaction regions (CIRs). We examine the Chandra and H alpha data for variability on this time scale. We find that the X-rays vary by about 15% over the course of the observations and that this variability is out of phase with variable absorption on the blue wing of the H alpha profiles (assumed to be a surrogate for the UV absorption associated with CIRs). While not conclusive, both sets of data are consistent with models where the CIRs are either a source of X-rays or modulate them.

The work is aimed at a study of the circumstellar disk of the bright classical binary Be star {\pi} Aqr. We analysed variations of a double-peaked profile of the H{\alpha} emission line in the spectrum of {\pi} Aqr that was observed in many phases during ~40 orbital cycles in 2004--2013. We applied the Discrete Fourier Transform (DFT) method to search for periodicity in the peak intensity (V/R) ratio. Doppler tomography was used to study the structure of the disk around the primary. The dominant frequency in the power spectrum of the H{\alpha} V/R ratio is 0.011873 day^-1 that correspond to a period of 84.2(2) days and is in agreement with the earlier determined orbital period of the system, Porb=84.1 days. The V/R ratio shows a sinusoidal variation phase-locked with the orbital period. Doppler maps of all our spectra show a non-uniform structure of the disk around the primary: a ring with the inner and outer radii at Vin~ 450 km/s and Vout~ 200km/s, respectively, along with an extended stable region (spot) at V_x ~ 225 km/s and V_y~100 km/s. The disk radius of ~ 65 Rsun = 0.33 AU was estimated assuming Keplerian motion of a particle on a circular orbit at the disk outer edge.

We present the first high angular resolution observations in the nearinfrared H-band (1.6 microns) of the Luminous Blue Variable star P Cygni. We obtained six-telescope interferometric observations with the CHARA Array and the MIRC beam combiner. These show that the spatial flux distribution is larger than expected for the stellar photosphere. A two component model for the star (uniform disk) plus a halo (two-dimensional Gaussian) yields an excellent fit of the observations, and we suggest that the halo corresponds to flux emitted from the base of the stellar wind. This wind component contributes about 45% of the H-band flux and has an angular FWHM = 0.96 mas, compared to the predicted stellar diameter of 0.41 mas. We show several images reconstructed from the interferometric visibilities and closure phases, and they indicate a generally spherical geometry for the wind. We also obtained near-infrared spectrophotometry of P Cygni from which we derive the flux excess compared to a purely photospheric spectral energy distribution. The H-band flux excess matches that from the wind flux fraction derived from the two component fits to the interferometry. We find evidence of significant near-infrared flux variability over the period from 2006 to 2010 that appears similar to the variations in the H-alpha emission flux from the wind. Future interferometric observations may be capable of recording the spatial variations associated with temporal changes in the wind structure.

We present the results of a time-series analysis of 130 echelle spectra of epsilon Ori (B0 Ia), acquired over seven observing seasons between 1998 and 2006 at Ritter Observatory. The equivalent widths of Halpha (net) and He I 5876 were measured and radial velocities were obtained from the central absorption of He I 5876. Temporal variance spectra (TVS) revealed significant wind variability in both Halpha and He I 5876. The He I TVS have a double-peaked profile consistent with radial velocity oscillations. A periodicity search was carried out on the equivalent width and radial velocity data, as well as on wavelength-binned spectra. This analysis has revealed several periods in the variability with time scales of 2-7 d. Many of these periods exhibit sinusoidal modulation in the associated phase diagrams. Several of these periods were present in both Halpha and He I, indicating a possible connection between the wind and the photosphere. Due to the harmonic nature of these periods, stellar pulsations may be the origin of some of the observed variability. Periods on the order of the rotational period were also detected in the He I line in the 98-99 season and in both lines during the 04-05 season. These periods may indicate rotational modulation due to structure in the wind.

We report on the preliminary results of a three-year monitoring campaign of Galactic and Magellanic LBVs with both moderate and high resolution spectroscopy. We have collected more than 500 moderate-resolution spectra of 20 stars in the red portion of the optical spectrum, which includes the wind-sensitive transitions of H$\alpha$ and He I 5876 and 6678 Angstroms. We summarize our detailed study of 15 years of high resolution H$\alpha$ spectroscopy of the prototypical luminous blue variable, P Cygni. We report on the discovery of discrete absorption components in P Cygni's H$\alpha$ profile, and we discuss their relationship to structure in the wind. These results are compared to our recent high resolution interferometric imaging of the $H$-band emitting region surrounding the star. We discuss recent observations of $\eta$ Car, highlighting its unusual recovery from the previous periastron passage. Our results on HDE 326823 indicate that the star has a short period binary orbit, with Roche Lobe overflow onto an unseen massive companion.

P Cygni is a prototype of the Luminous Blue Variables (or S Doradus variables), and the star displays photometric and emission line variability on a timescale of years (known as the "short S Doradus phase" variations). Here we present new high resolution H-alpha spectroscopy of P Cyg that we combine with earlier spectra and concurrent V-band photometry to document the emission and continuum flux variations over a 24 y time span. We show that the emission and continuum fluxes vary in concert on timescales of 1.6 y and longer, but differ on shorter timescales. The H-alpha profile shape also varies on the photometric timescales, and we describe the observed co-variations of the emission peak and absorption trough properties. We argue that the episodes of photometric and emission brightening are caused by increases in the size of the emission region that are related to variations in wind mass loss rate and outflow speed. We find evidence of blueward accelerating, Discrete Absorption Components (DACs) in the absorption trough of the H-alpha profile, and these features have slower accelerations and longer durations than those observed in other lines. The DAC strengths also appear to vary on the photometric timescales, and we suggest that the propagation of the DAC-related wind structures is closely related to changes in the overall wind mass loss rate and velocity.

Deneb is often considered the prototypical A-type supergiant, and is one of the visually most luminous stars in the Galaxy. A-type supergiants are potential extragalactic distance indicators, but the variability of these stars needs to be better characterized before this technique can be considered reliable. We analyzed 339 high resolution echelle spectra of Deneb obtained over the five-year span of 1997 through 2001 as well as 370 Stromgren photometric measurements obtained during the same time frame. Our spectroscopic analysis included dynamical spectra of the H-alpha profile, H-alpha equivalent widths, and radial velocities measured from Si II 6347, 6371. Time-series analysis reveals no obvious cyclic behavior that proceeds through multiple observing seasons, although we found a suspected 40 day period in two, non-consecutive observing seasons. Some correlations are found between photometric and radial velocity data sets, and suggest radial pulsations at two epochs. No correlation is found between the variability of the H-alpha profiles and that of the radial velocities or the photometry. Lucy (1976) found evidence that Deneb was a long period single-lined spectroscopic binary star, but our data set shows no evidence for radial velocity variations caused by a binary companion.

We provide a quantitative analysis of time-variable phenomena in the photospheric, near-star, and outflow regions of the late-B supergiant (SG) HD 199478. The analysis is based primarily on optical spectroscopic datasets secured between 1999 and 2000 from the Bulgarian NAO, Tartu, and Ritter Observatories. The temporal behaviour of HD 199478 is characterised by three key empirical properties: (i) systematic central velocity shifts in the photospheric absorption lines, including C II and He I, over a characteristic time-scale of abou 20 days; (ii) extremely strong, variable H alpha emission with no clear modulation signal, and (iii) the occurrence in 2000 of a (rare) high-velocity absorption (HVA) event in H alpha, which evolved over about 60 days, showing the clear signature of mass infall and outflows. In these properties HD 199478 resembles few other late-B SGs with peculiar emission and HVAs in H alpha (HD 91619, HD 34085, HD 96919). Non-LTE line synthesis modelling is conducted using FASTWIND for these late-B SGs to constrain and compare their fundamental parameters within the context of extreme behaviour in the H alpha lines. Our analysis indicate that at the cooler temperature edge of B SGs, there are objects whose wind properties, as traced by H alpha, are inconsistent with the predictions of the smooth, spherically symmetric wind approximation. This discordance is still not fully understood and may highlight the role of a non-spherical, disk-like, geometry, which may result from magnetically-driven equatorial compression of the gas. Ordered dipole magnetic fields may also lead to confined plasma held above the stellar surface, which ultimately gives rise to transient HVA events.

We present the results of an optical spectroscopic study of the massive Wolf-Rayet binary WR137. These data cover the dust-formation maximum in 1997. Combining all available measurements of radial velocities, we derive, for the first time, a spectroscopic orbit with period 13.05 +/- 0.18 years. The resulting masses, adopting i=67degrees, are M(O)= 20 +/- 2 Mo and M(WR)= 4.4+/- 1.5 Mo. These appear, respectively, around normal and on the low side for the given spectral types. Analysis of the intense multi-site spectroscopic monitoring in 1999 shows that the CIII5696 and CIV5802/12 lines have the highest intrinsic variability levels. The periodogram analysis yields a small-amplitude modulation in the absorption troughs of the CIV5802/12 and HeI5876 lines with a period of 0.83 days, which could be related either to pulsations or large-scale rotating structures as seen in the WN4 star EZ Canis Majoris (WR6).Wavelet analysis of the strong emission lines of CIII5696 and CIV5802/12 enabled us to isolate and follow for several hours small structures associated with density enhancements within the wind of the Wolf-Rayet star. Cross-correlating the variability patterns seen in different lines, we find a weak but significant correlation between the varability in emission lines with different ionization potential, i.e. in lines formed at different distances from the WR stellar core. Adopting a beta wind-velocity law, from the motion of individual subpeaks we find beta around 5, which is significantly larger than the canonical value of 1 usually found in O-star winds.

We present an overview of the spectral variability of the peculiar F-type hypergiant Rho Cas, obtained from our long-term monitoring campaigns over the past 8.5 years with four spectrographs in the northern hemisphere. Between 2000 June and September an exceptional variability phase occurred when the V-brightness dimmed by about a full magnitude. The star recovered from this deep minimum by 2001 April. It is the third outburst of Rho Cas on record in the last century. We observe TiO absorption bands in high-resolution near-IR spectra obtained with the Utrecht Echelle Spectrograph during the summer of 2000. TiO formation in the outer atmosphere occurred before the deep brightness minimum. Atmospheric models reveal that the effective temperature decreases by at least 3000 K, and the TiO shell is driven supersonically with Mdot ~= 5.4 10^{-2} Msun/yr. Strong episodic mass loss and TiO have also been observed during the outbursts of 1945-47 and 1985-86. A detailed analysis of the exceptional outburst spectra is provided, by comparing with high-resolution optical spectra of the early M-type supergiants Mu Cep (Ia) and Betelgeuse (Iab). The outburst spectra indicate the formation of a low-temperature, optically thick circumstellar gas shell of 3 10^{-2} Msun during ~200 d, caused by dynamic instability of the upper atmosphere of this pulsating massive supergiant near the Eddington luminosity limit. We observe that the mass-loss rate during the outburst is of the same order of magnitude as has been proposed for the outbursts of Eta Carinae. We present calculations that correctly predict the outburst time-scale, whereby the shell ejection is driven by the release of hydrogen ionization-recombination energy.

The spectroscopic and spectropolarimetric variability of the peculiar variable V838 Monocerotis during the brighter phases of its multiple outbursts in 2002 is presented. Significant line profile variability of H$\alpha$ and Si II 6347.10\AA & 6371.36\AA occurred in spectra obtained between 2002 February 5 and 2002 March 14, and a unique secondary absorption component was observed near the end of this time period. Our observations also suggest that multiple shifts in ionization states occurred during the outbursts. Spectropolarimetric observations reveal that V838 Mon exhibited both intrinsic and interstellar polarization components during the initial stages of the second outburst, indicating the presence of an asymmetric geometry; however, the intrinsic component had significantly declined by February 14. We determine the interstellar polarization to be $P_{max} = 2.746 \pm 0.011 %$, $\lambda_{max} = 5790 \pm 37\AA$, $PA = 153.43 \pm 0.12 ^{\circ} $, and we find the integrated intrinsic V band polarization on February 5 to be $P = 0.983 \pm 0.012 %$ at a position angle of $127.0 \pm 0.5^{\circ}$. The implications of these observations for the nature of V838 Monocerotis, its distance, and its ejecta are discussed.

Between 2000 June and September an exceptional variability phase occurred in the peculiar F-type hypergiant Rho Cas, when the V-brightness dimmed by at least a full magnitude. The star recovered from this deep minimum by 2001 April. It is the third outburst of Rho Cas on record in the last century. We observe TiO absorption bands in high-resolution optical and near-IR spectra obtained with the Utrecht Echelle Spectrograph during the summer of 2000. TiO formation in the outer atmosphere occurred before the deep brightness minimum.Atmospheric models reveal that the effective temperature decreases by at least 3000 K, and the TiO shell is driven supersonically with Mdot =~5.4 10^-2 Msun/yr. Strong episodic mass loss and TiO have also been observed during the outbursts of 1945-47 and 1985-86. An analysis of the exceptional outburst spectra of 2000-01 is provided, by comparing with high-resolution optical spectra of the early M-type supergiants Mu Cep (Ia) and Betelgeuse (Iab). During the outburst, central emission appears above the local continuum level in the split Na D lines. The outburst spectra indicate the formation of a low- temperature, optically thick circumstellar gas shell of 3 10^-2 Msun during \~200 d, caused by dynamic instability of the upper atmosphere of this pulsating massive supergiant near the Eddington luminosity limit. We present an equation that correctly predicts the outburst time-scale, whereby the shell ejection is driven by the release of hydrogen recombination energy. We observe that the mass-loss rate during the outburst is of the same order of magnitude as has been proposed for the outbursts of Eta Carinae. The research results in this paper are described in further detail in Lobel et al. (2003). A spectral movie sequence of the outburst is available at http://cfa-www.harvard.edu/~alobel/

The bright star delta Sco has been considered a typical B0-type object for many years. Spectra of the star published prior to 1990 showed no evidence of emission, but only of short-term line profile variations attributed to nonradial pulsations. Speckle interferometric observations show that delta Sco is a binary system with a highly-eccentric orbit and a period of 10.6 years. Weak emission in the H-alpha line was detected in its spectrum for the first time during a periastron passage in 1990. Shortly before the next periastron passage in the summer of 2000, the binary entered a strong H-alpha emission and enhanced mass loss phase. We monitored the spectroscopic development of the Be outburst from July 2000 through March 2001. In this paper we present results from our spectroscopy, refine elements of the binary orbit, and discuss possible mechanisms for the mass loss.