Showing 1–50 of 76 results for author: Massa, D

Wind line variability and intrinsic errors in observational mass loss rates

Authors: Derck Massa, Raman K. Prinja, Lidia Oskinova

Abstract: UV wind line variability in OB stars appears to be universal. To quantify this variation and to estimate its effect on a mass loss rate determined from a single observation, we use the IUE archive to identify non-peculiar OB stars with well developed but unsaturated Si IV 1400 doublets and at least 10 independent observations. This resulted in 1699 spectra of 25 stars. We use SEI modelling to tran… ▽ More UV wind line variability in OB stars appears to be universal. To quantify this variation and to estimate its effect on a mass loss rate determined from a single observation, we use the IUE archive to identify non-peculiar OB stars with well developed but unsaturated Si IV 1400 doublets and at least 10 independent observations. This resulted in 1699 spectra of 25 stars. We use SEI modelling to translate the profile variations into optical depth variations and, hence, variations in measured mass loss rates. These variations quantify the intrinsic error inherent in any mass loss rate derived from a single observation. The derived rates have an overall RMS variation of about 22%, but this differs with effective temperature, being as small at 8% for the hottest stars and up to 45% for the cooler ones. Furthermore, any single determination can differ from the mean by a factor of 2 or more. Our results also imply that mass loss rates determined from non-simultaneous observations (such as UV and ground based data) need not agree. We also use our results to examine the nature of the wind structures responsible for the variability. Our findings suggest that the optical depth variations result from optically very thick structures occulting more or less of the line of sight to the stellar disk. Further, the smaller optical depth variations in the hottest stars suggests that the responsible structures are disrupted in their more powerful winds. △ Less

Submitted 16 June, 2024; originally announced July 2024.

Comments: 23 pages, 8 figures, Accepted by ApJ

In-Context Learning of Physical Properties: Few-Shot Adaptation to Out-of-Distribution Molecular Graphs

Authors: Grzegorz Kaszuba, Amirhossein D. Naghdi, Dario Massa, Stefanos Papanikolaou, Andrzej Jaszkiewicz, Piotr Sankowski

Abstract: Large language models manifest the ability of few-shot adaptation to a sequence of provided examples. This behavior, known as in-context learning, allows for performing nontrivial machine learning tasks during inference only. In this work, we address the question: can we leverage in-context learning to predict out-of-distribution materials properties? However, this would not be possible for struct… ▽ More Large language models manifest the ability of few-shot adaptation to a sequence of provided examples. This behavior, known as in-context learning, allows for performing nontrivial machine learning tasks during inference only. In this work, we address the question: can we leverage in-context learning to predict out-of-distribution materials properties? However, this would not be possible for structure property prediction tasks unless an effective method is found to pass atomic-level geometric features to the transformer model. To address this problem, we employ a compound model in which GPT-2 acts on the output of geometry-aware graph neural networks to adapt in-context information. To demonstrate our model's capabilities, we partition the QM9 dataset into sequences of molecules that share a common substructure and use them for in-context learning. This approach significantly improves the performance of the model on out-of-distribution examples, surpassing the one of general graph neural network models. △ Less

Submitted 3 June, 2024; originally announced June 2024.

Comments: 12 pages, 4 figures

Expanded Sample of Small Magellanic Cloud Ultraviolet Dust Extinction Curves: Correlations between the 2175 A bump, q_pah, UV extinction shape, and N(HI)/A(V)

Authors: Karl D. Gordon, E. L. Fitzpatrick, Derck Massa, Ralph Bohlin, Jeremy Chastenet, Claire E. Murray, Geoffrey C. Clayton, Daniel J. Lennon, Karl A. Misselt, Karin Sandstrom

Abstract: The Small Magellanic Cloud (SMC) shows a large variation in ultraviolet (UV) dust extinction curves, ranging from Milky Way-like (MW) to significantly steeper curves with no detectable 2175 A bump. This result is based on a sample of only nine sightlines. From HST/STIS and IUE spectra of OB stars, we have measured UV extinction curves along 32 SMC sightlines where eight of these curves were publis… ▽ More The Small Magellanic Cloud (SMC) shows a large variation in ultraviolet (UV) dust extinction curves, ranging from Milky Way-like (MW) to significantly steeper curves with no detectable 2175 A bump. This result is based on a sample of only nine sightlines. From HST/STIS and IUE spectra of OB stars, we have measured UV extinction curves along 32 SMC sightlines where eight of these curves were published previously. We find 16 sightlines with steep extinction with no detectable 2175 A bump, four sightlines with MW-like extinction with a detectable 2175 A bump, two sightlines with fairly flat UV extinction and weak/absent 2175 A bumps, and 10 sightlines with unreliable curves due to low SMC dust columns. Our expanded sample shows that the sightlines with and without the 2175 A bump are located throughout the SMC and not limited to specific regions. The average extinction curve of the 16 bumpless sightlines is very similar to the previous average based on four sightlines. We find no correlation between dust column and the strength of the 2175 A bump. We test the hypothesis that the 2175 A bump is due to the same dust grains that are responsible for the mid-infrared carbonaceous (PAH) emission features and find they are correlated, confirming recent work in the MW. Overall, the slope of the UV extinction increases as the amplitudes of the 2175 A bump and far-UV curvature decrease. Finally, the UV slope is correlated with $N(HI)/A(V)$ and the 2175 A bump and nonlinear far-UV rise amplitudes are anti-correlated with $N(HI)/A(V)$. △ Less

Submitted 21 May, 2024; originally announced May 2024.

Comments: 24 pages, 12 figures, ApJ, accepted

Material Informatics through Neural Networks on Ab-Initio Electron Charge Densities: the Role of Transfer Learning

Authors: Dario Massa, Stefanos Papanikolaou, Piotr Sankowski

Abstract: In this work, the dynamic realms of Materials Science and Computer Science advancements meet the critical challenge of identifying efficient descriptors capable of capturing the essential features of physical systems. Such task has remained formidable, with solutions often involving ad-hoc scalar and vectorial sets of materials properties, making optimization and transferability challenging. We ex… ▽ More In this work, the dynamic realms of Materials Science and Computer Science advancements meet the critical challenge of identifying efficient descriptors capable of capturing the essential features of physical systems. Such task has remained formidable, with solutions often involving ad-hoc scalar and vectorial sets of materials properties, making optimization and transferability challenging. We extract representations directly from ab-initio differential electron charge density profiles using Neural Networks, highlighting the pivotal role of transfer learning in such task. Firstly, we demonstrate significant improvements in regression of a specific defected-materials property with respect to training a deep network from scratch, both in terms of predictions and their reproducibilities, by considering various pre-trained models and selecting the optimal one after fine-tuning. The remarkable performances obtained confirmed the transferability of the existent pre-trained Convolutional Neural Networks (CNNs) on physics domain data, very different from the original training data. Secondly, we demonstrate a saturation in the regression capabilities of computer vision models towards properties of an extensive variety of undefected systems, and how it can be overcome with the help of large language model (LLM) transformers, with as little text information as composition names. Finally, we prove the insufficiency of open-models, like GPT-4, in achieving the analogous tasks and performances as the proposed domain-specific ones. The work offers a promising avenue for enhancing the effectiveness of descriptor identification in complex physical systems, shedding light over the power of transfer learning to easily adapt and combine available models, with different modalities, to the physics domain, at the same time opening space to a benchmark for LLMs capabilities in such domain. △ Less

Submitted 17 January, 2024; originally announced January 2024.

Unbiased Atomistic Predictions of Crystal Dislocation Dynamics using Bayesian Force Fields

Authors: Cameron J. Owen, Amirhossein D. Naghdi, Anders Johansson, Dario Massa, Stefanos Papanikolaou, Boris Kozinsky

Abstract: Crystal dislocation dynamics, especially at high temperatures, represents a subject where experimental phenomenological input is commonly required, and parameter-free predictions, starting from quantum methods, have been beyond reach. This is especially true for phenomena like stacking faults and dislocation cross-slip, which are computationally intractable with methods like density functional the… ▽ More Crystal dislocation dynamics, especially at high temperatures, represents a subject where experimental phenomenological input is commonly required, and parameter-free predictions, starting from quantum methods, have been beyond reach. This is especially true for phenomena like stacking faults and dislocation cross-slip, which are computationally intractable with methods like density functional theory, as $\sim 10^5-10^6$ atoms are required to reliably simulate such systems. Hence, this work extends quantum-mechanical accuracy to mesoscopic molecular dynamics simulations and opens unprecedented possibilities in material design for extreme mechanical conditions with direct atomistic insight at the deformation mesoscale. To accomplish this, we construct a Bayesian machine-learned force field (MLFF) from ab initio quantum training data, enabling direct observations of high-temperature and high-stress dislocation dynamics in single-crystalline Cu with atomistic resolution. In doing so, a generalizable training protocol is developed for construction of MLFFs for dislocation kinetics, with wide-ranging applicability to other single element systems and alloys. The resulting FLARE MLFF provides excellent predictions of static bulk elastic properties, stacking fault widths and energies, dynamic evolutions and mobilities of edge and screw dislocations, as well as cross-slip energy barriers for screw dislocations, all of which are compared to available experimental measurements. This work marks the first reliable quantitative determination of dislocation mobilities and cross-slip barriers, demonstrating a substantial advancement over previous empirical and machine learned force field models. △ Less

Submitted 9 January, 2024; originally announced January 2024.

Comments: 15 pages, 7 figures, 3 tables

Optically-thick Structure in Early B Type Supergiant Stellar Winds at Low Metallicities

Authors: Timothy N. Parsons, Raman K. Prinja, Matheus Bernini-Peron, Alex W. Fullerton, Derck L. Massa, Lidia M. Oskinova, Daniel Pauli, Andreas A. C. Sander, Matthew J. Rickard

Abstract: Accurate determination of mass-loss rates from massive stars is important to understanding stellar and galactic evolution and enrichment of the interstellar medium. Large-scale structure and variability in stellar winds have significant effects on mass-loss rates. Time-series observations provide direct quantification of such variability. Observations of this nature are available for some Galactic… ▽ More Accurate determination of mass-loss rates from massive stars is important to understanding stellar and galactic evolution and enrichment of the interstellar medium. Large-scale structure and variability in stellar winds have significant effects on mass-loss rates. Time-series observations provide direct quantification of such variability. Observations of this nature are available for some Galactic early supergiant stars but not yet for stars in lower metallicity environments such as the Magellanic Clouds. We utilise ultraviolet spectra from the Hubble Space Telescope ULLYSES program to demonstrate that the presence of structure in stellar winds of supergiant stars at low metallicities may be discerned from single-epoch spectra. We find evidence that, for given stellar luminosities and mean stellar wind optical depths, structure is more prevalent at higher metallicities. We confirm, at Large Magellanic Cloud (0.5 Z_solar), Small Magellanic Cloud (0.2 Z_solar) and lower (0.14 -- 0.1 Z_solar) metallicities, earlier Galactic results that there does not appear to be correlation between the degree of structure in stellar winds of massive stars and stellar effective temperature. Similar lack of correlation is found with regard to terminal velocity of stellar winds. Additional and revised values for radial velocities of stars and terminal velocities of stellar winds are presented. Direct evidence of temporal variability, on timescales of several days, in stellar wind at low metallicity is found. We illustrate that narrow absorption components in wind-formed profiles of Galactic OB stellar spectra remain common in early B supergiant spectra at low metallicities, providing means for better constraining hot, massive star mass-loss rates. △ Less

Submitted 2 January, 2024; originally announced January 2024.

Comments: 18 pages plus Appendices, 9 figures plus 37 figures in Appendices. Accepted by MNRAS

Neural Network Interatomic Potentials For Open Surface Nano-mechanics Applications

Authors: Amirhossein D. Naghdi, Franco Pellegrini, Emine Küçükbenli, Dario Massa, F. Javier Dominguez Gutierrez, Efthimios Kaxiras, Stefanos Papanikolaou

Abstract: Material characterization in nano-mechanical tests requires precise interatomic potentials for the computation of atomic energies and forces with near-quantum accuracy. For such purposes, we develop a robust neural-network interatomic potential (NNIP), and we provide a test for the example of molecular dynamics (MD) nanoindentation, and the case of body-centered cubic crystalline molybdenum (Mo).… ▽ More Material characterization in nano-mechanical tests requires precise interatomic potentials for the computation of atomic energies and forces with near-quantum accuracy. For such purposes, we develop a robust neural-network interatomic potential (NNIP), and we provide a test for the example of molecular dynamics (MD) nanoindentation, and the case of body-centered cubic crystalline molybdenum (Mo). We employ a similarity measurement protocol, using standard local environment descriptors, to select ab initio configurations for the training dataset that capture the behavior of the indented sample. We find that it is critical to include generalized stacking fault (GSF) configurations, featuring a dumbbell interstitial on the surface, to capture dislocation cores, and also high-temperature configurations with frozen atom layers for the indenter tip contact. We develop a NNIP with distinct dislocation nucleation mechanisms, realistic generalized stacking fault energy (GSFE) curves, and an informative energy landscape for the atoms on the sample surface during nanoindentation. We compare our NNIP results with nanoindentation simulations, performed with three existing potentials -- an embedded atom method (EAM) potential, a gaussian approximation potential (GAP), and a tabulated GAP (tabGAP) potential -- that predict different dislocation nucleation mechanisms, and display the absence of essential information on the shear stress at the sample surface in the elastic region. We believe that these features render specialized NNIPs essential for simulations of nanoindentation and nano-mechanics with near-quantum accuracy. △ Less

Submitted 17 December, 2023; originally announced December 2023.

Alloy Informatics through Ab Initio Charge Density Profiles: Case Study of Hydrogen Effects in Face-Centered Cubic Crystals

Authors: Dario Massa, Efthimios Kaxiras, Stefanos Papanikolaou

Abstract: Materials design has traditionally evolved through trial-error approaches, mainly due to the non-local relationship between microstructures and properties such as strength and toughness. We propose 'alloy informatics' as a machine learning based prototype predictive approach for alloys and compounds, using electron charge density profiles derived from first-principle calculations. We demonstrate t… ▽ More Materials design has traditionally evolved through trial-error approaches, mainly due to the non-local relationship between microstructures and properties such as strength and toughness. We propose 'alloy informatics' as a machine learning based prototype predictive approach for alloys and compounds, using electron charge density profiles derived from first-principle calculations. We demonstrate this framework in the case of hydrogen interstitials in face-centered cubic crystals, showing that their differential electron charge density profiles capture crystal properties and defect-crystal interaction properties. Radial Distribution Functions (RDFs) of defect-induced differential charge density perturbations highlight the resulting screening effect, and, together with hydrogen Bader charges, strongly correlate to a large set of atomic properties of the metal species forming the bulk crystal. We observe the spontaneous emergence of classes of charge responses while coarse-graining over crystal compositions. Nudge-Elastic-Band calculations show that RDFs and charge features also connect to hydrogen migration energy barriers between interstitial sites. Unsupervised machine-learning on RDFs supports classification, unveiling compositional and configurational non-localities in the similarities of the perturbed densities. Electron charge density perturbations may be considered as bias-free descriptors for a large variety of defects. △ Less

Submitted 24 September, 2023; originally announced September 2023.

Substitutional Alloying Using Crystal Graph Neural Networks

Authors: Dario Massa, Daniel Cieśliński, Amirhossein Naghdi, Stefanos Papanikolaou

Abstract: Materials discovery, especially for applications that require extreme operating conditions, requires extensive testing that naturally limits the ability to inquire the wealth of possible compositions. Machine Learning (ML) has nowadays a well established role in facilitating this effort in systematic ways. The increasing amount of available accurate DFT data represents a solid basis upon which new… ▽ More Materials discovery, especially for applications that require extreme operating conditions, requires extensive testing that naturally limits the ability to inquire the wealth of possible compositions. Machine Learning (ML) has nowadays a well established role in facilitating this effort in systematic ways. The increasing amount of available accurate DFT data represents a solid basis upon which new ML models can be trained and tested. While conventional models rely on static descriptors, generally suitable for a limited class of systems, the flexibility of Graph Neural Networks (GNNs) allows for direct learning representations on graphs, such as the ones formed by crystals. We utilize crystal graph neural networks (CGNN) to predict crystal properties with DFT level accuracy, through graphs with encoding of the atomic (node/vertex), bond (edge), and global state attributes. In this work, we aim at testing the ability of the CGNN MegNet framework in predicting a number of properties of systems previously unseen from the model, obtained by adding a substitutional defect in bulk crystals that are included in the training set. We perform DFT validation to assess the accuracy in the prediction of formation energies and structural features (such as elastic moduli). Using CGNNs, one may identify promising paths in alloy discovery. △ Less

Submitted 19 June, 2023; originally announced June 2023.

Comments: 34 pages, 14 figures, Supplementary Material

One Relation for All Wavelengths: The Far-Ultraviolet to Mid-Infrared Milky Way Spectroscopic R(V) Dependent Dust Extinction Relationship

Authors: Karl D. Gordon, Geoffrey C. Clayton, Marjorie Decleir, E. L. Fitzpatrick, Derck Massa, Karl A. Misselt, Erik J. Tollerud

Abstract: Dust extinction is one of the fundamental measurements of dust grain sizes, compositions, and shapes. Most of the wavelength dependent variations seen in Milky Way extinction are strongly correlated with the single parameter R(V)=A(V)/E(B-V). Existing R(V) dependent extinction relationships use a mixture of spectroscopic and photometry observations, hence do not fully capture all the important dus… ▽ More Dust extinction is one of the fundamental measurements of dust grain sizes, compositions, and shapes. Most of the wavelength dependent variations seen in Milky Way extinction are strongly correlated with the single parameter R(V)=A(V)/E(B-V). Existing R(V) dependent extinction relationships use a mixture of spectroscopic and photometry observations, hence do not fully capture all the important dust features nor continuum variations. Using four existing samples of spectroscopically measured dust extinction curves, we consistently measure the R(V) dependent extinction relationship spectroscopically from the far-ultraviolet to mid-infrared for the first time. Linear fits of A(lambda)/A(V) dependent on R(V) are done using a method that fully accounts for their significant and correlated uncertainties. These linear parameters are fit with analytic wavelength dependent functions to determine the smooth R(V) (2.3-5.6) and wavelength (912 A-32 micron) dependent extinction relationship. This relationship shows that the far-UV rise, 2175 A bump, and the three broad optical features are dependent on R(V), but the 10 and 20 micron features are not. Existing literature relationships show significant deviations compared to this relationship especially in the far-ultraviolet and infrared. Extinction curves that clearly deviate from this relationship illustrate that this relationship only describes the average behavior versus R(V). We find tentative evidence that the relationship may not be linear with 1/R(V) especially in the ultraviolet. For the first time, this relationship provides measurements of dust extinction that spectroscopically resolve the continuum and features in the ultraviolet, optical, and infrared as a function of R(V) enabling detailed studies of dust grains properties and full spectroscopic accounting for the effects of dust extinction on astrophysical objects. △ Less

Submitted 6 April, 2023; v1 submitted 4 April, 2023; originally announced April 2023.

Comments: 16 pages, 9 figures, ApJ, in press

Caught in the Act: A Metal-Rich High-Velocity Cloud in the Inner Galaxy

Authors: Frances H. Cashman, Andrew J. Fox, Bart P. Wakker, Trisha Ashley, Derck Massa, Edward B. Jenkins, Dhanesh Krishnarao, Felix J. Lockman, Robert A. Benjamin, Rongmon Bordoloi, Tae-Sun Kim

Abstract: We characterize the chemical and physical conditions in an outflowing high-velocity cloud in the inner Galaxy. We report a super-solar metallicity of [O/H] = $+0.36\pm0.12$ for the high-velocity cloud at $v_\mathrm{LSR}$ = 125.6 km s$^{-1}$ toward the star HD 156359 ($l$ = 328.$^{\circ}$7, $b$ = $-$14.$^{\circ}$5, $d$ = 9 kpc, $z$ = $-$2.3 kpc). Using archival observations from FUSE, HST STIS, and… ▽ More We characterize the chemical and physical conditions in an outflowing high-velocity cloud in the inner Galaxy. We report a super-solar metallicity of [O/H] = $+0.36\pm0.12$ for the high-velocity cloud at $v_\mathrm{LSR}$ = 125.6 km s$^{-1}$ toward the star HD 156359 ($l$ = 328.$^{\circ}$7, $b$ = $-$14.$^{\circ}$5, $d$ = 9 kpc, $z$ = $-$2.3 kpc). Using archival observations from FUSE, HST STIS, and ESO FEROS we measure high-velocity absorption in H I, O I, C II, N II, Si II, Ca II, Si III, Fe III, C IV, Si IV, N V, and O VI. We measure a low H I column density of log $N$(H I) = $15.54\pm0.05$ in the HVC from multiple unsaturated H I Lyman series lines in the FUSE data. We determine a low dust depletion level in the HVC from the relative strength of silicon, iron, and calcium absorption relative to oxygen, with [Si/O]=$-0.33\pm0.14$, [Fe/O]=$-0.30\pm0.20$, and [Ca/O] =$-0.56\pm0.16$. Analysis of the high-ion absorption using collisional ionization models indicates that the hot plasma is multi-phase, with the C IV and Si IV tracing 10$^{4.9}$ K gas and N V and O VI tracing 10$^{5.4}$ K gas. The cloud's metallicity, dust content, kinematics, and close proximity to the disk are all consistent with a Galactic wind origin. As the HD 156359 line of sight probes the inner Galaxy, the HVC appears to be a young cloud caught in the act of being entrained in a multi-phase Galactic outflow and driven out into the halo. △ Less

Submitted 20 October, 2022; originally announced October 2022.

Comments: Submitted to ApJ, 18 pages, 13 figures, comments welcome

Phase-resolved spectroscopic analysis of the eclipsing black hole X-ray binary M33 X-7: System properties, accretion, and evolution

Authors: V. Ramachandran, L. M. Oskinova, W. -R. Hamann, A. A. C. Sander, H. Todt, D. Pauli, T. Shenar, J. M. Torrejón, K. A. Postnov, J. M. Blondin, E. Bozzo, R. Hainich, D. Massa

Abstract: M33 X-7 is the only known eclipsing black hole high mass X-ray binary. The system is reported to contain a very massive O supergiant donor and a massive black hole in a short orbit. The high X-ray luminosity and its location in the metal-poor galaxy M33 make it a unique laboratory for studying the winds of metal-poor donor stars with black hole companions and it helps us to understand the potentia… ▽ More M33 X-7 is the only known eclipsing black hole high mass X-ray binary. The system is reported to contain a very massive O supergiant donor and a massive black hole in a short orbit. The high X-ray luminosity and its location in the metal-poor galaxy M33 make it a unique laboratory for studying the winds of metal-poor donor stars with black hole companions and it helps us to understand the potential progenitors of black hole mergers. Using phase-resolved simultaneous HST- and XMM-Newton-observations, we traced the interaction of the stellar wind with the black hole. Our comprehensive spectroscopic investigation of the donor star (X-ray+UV+optical) yields new stellar and wind parameters for the system that differs significantly from previous estimates. In particular, the masses of the components are considerably reduced to 38 for the O-star donor and 11.4 for the black hole. The O giant is overfilling its Roche lobe and shows surface He enrichment. The donor shows a densely clumped wind with a mass-loss rate that matches theoretical predictions. We investigated the wind-driving contributions from different ions and the changes in the ionization structure due to X-ray illumination. Toward the black hole, the wind is strongly quenched due to strong X-ray illumination. For this system, the standard wind-fed accretion scenario alone cannot explain the observed X-ray luminosity, pointing toward an additional mass overflow, which is in line with our acceleration calculations. The X-ray photoionization creates an He II emission region emitting $10^{47}$ ph/s. We computed binary evolutionary tracks for the system using MESA. Currently, the system is transitioning toward an unstable mass transfer phase, resulting in a common envelope of the black hole and donor. Since the mass ratio is q~3.3 and the period is short, the system is unlikely to survive the common envelope, but will rather merge. △ Less

Submitted 16 August, 2022; originally announced August 2022.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 667, A77 (2022)

Wind line variability and intrinsic errors in observational mass loss rates

Authors: Derck Massa, Raman K. Prinja

Abstract: UV wind line variability in OB stars appears to be universal. We review the evidence that the variability is due to large, dense, optically thick structures rooted in or near the photosphere. Using repeated bservations and a simple model we translate observed profile variations into optical depth variations and, consequently, variations in measured mass loss rates. Although global rates may be sta… ▽ More UV wind line variability in OB stars appears to be universal. We review the evidence that the variability is due to large, dense, optically thick structures rooted in or near the photosphere. Using repeated bservations and a simple model we translate observed profile variations into optical depth variations and, consequently, variations in measured mass loss rates. Although global rates may be stable, measured rates vary. Consequently, profile variations infer how mass loss rates determined from UV wind lines vary. These variations quantify the intrinsic error inherent in any mass loss rate derived from a single observation. These derived rates can differ by factors of 3 or more. Our results also imply that rates from non-simultaneous observations (such as UV and ground based data) need not agree. Finally, we use our results to examine the nature of the structures responsible for the variability. △ Less

Submitted 28 July, 2022; originally announced July 2022.

Comments: 6 pages, 5 figures, paper presented at IAU Symp 361 "Massive Stars Near and Far" (eds. Nicole St-Louis, Jorick Vink, Jonathan Mackey)

Relations between mid-IR dust emission and UV extinction

Authors: Derck Massa, Karl D. Gordon, E. L. Fitzpatrick

Abstract: We analyze low resolution Spitzer infrared (IR) 5-14 micron spectra of the diffuse emission toward a carefully selected sample of stars. The sample is composed of sight lines toward stars that have well determined ultraviolet (UV) extinction curves and which are shown to lie beyond effectively all of the extinguishing and emitting dust along their lines of sight. Our sample includes sight lines wh… ▽ More We analyze low resolution Spitzer infrared (IR) 5-14 micron spectra of the diffuse emission toward a carefully selected sample of stars. The sample is composed of sight lines toward stars that have well determined ultraviolet (UV) extinction curves and which are shown to lie beyond effectively all of the extinguishing and emitting dust along their lines of sight. Our sample includes sight lines whose UV curve extinction curves exhibit a wide range of curve morphology and which sample a variety of interstellar environments. As a result, this unique sample enabled us to study the connection between the extinction and emission properties of the same grains, and to examine their response to different physical environments. We quantify the emission features in terms of the PAH model given by Draine & Li (2007) and a set of additional features, not known to be related to PAH emission. We compare the intensities of the different features in the Spitzer mid-IR spectra with the Fitzpatrick & Massa (2007) parameters which describe the shapes of UV to near-IR extinction curves. Our primary result is that there is a strong correlation between the area of the 2175 A UV bump in the extinction curves of the program stars and the strengths of the major PAH emission features in the mid-IR spectra for the same lines of sight. △ Less

Submitted 5 December, 2021; originally announced December 2021.

Comments: 29 pages, 17 figures, accepted to the ApJ

DNS Covert Channel Detection via Behavioral Analysis: a Machine Learning Approach

Authors: Salvatore Saeli, Federica Bisio, Pierangelo Lombardo, Danilo Massa

Abstract: Detecting covert channels among legitimate traffic represents a severe challenge due to the high heterogeneity of networks. Therefore, we propose an effective covert channel detection method, based on the analysis of DNS network data passively extracted from a network monitoring system. The framework is based on a machine learning module and on the extraction of specific anomaly indicators able to… ▽ More Detecting covert channels among legitimate traffic represents a severe challenge due to the high heterogeneity of networks. Therefore, we propose an effective covert channel detection method, based on the analysis of DNS network data passively extracted from a network monitoring system. The framework is based on a machine learning module and on the extraction of specific anomaly indicators able to describe the problem at hand. The contribution of this paper is two-fold: (i) the machine learning models encompass network profiles tailored to the network users, and not to the single query events, hence allowing for the creation of behavioral profiles and spotting possible deviations from the normal baseline; (ii) models are created in an unsupervised mode, thus allowing for the identification of zero-days attacks and avoiding the requirement of signatures or heuristics for new variants. The proposed solution has been evaluated over a 15-day-long experimental session with the injection of traffic that covers the most relevant exfiltration and tunneling attacks: all the malicious variants were detected, while producing a low false-positive rate during the same period. △ Less

Submitted 4 October, 2020; originally announced October 2020.

Comments: This is a pre-print of an article published in the proceedings of 14th International Conference on Malicious and Unwanted Software (MALWARE), 2019

MSC Class: 63P30; 68T99 ACM Class: C.2.0; I.5.0

An Analysis of the Shapes of Interstellar Extinction Curves. VIII. The Optical Extinction Structure

Authors: D. Massa, E. L. Fitzpatrick, K. D. Gordon

Abstract: New HST/STIS optical spectra were obtained for a sample of early type stars with existing IUE UV spectra. These data were used to construct optical extinction curves whose general properties are discussed elsewhere. In this paper, we identify extinction features in the curves that are wider than diffuse interstellar bands (DIBs) but narrower than the well known broad band variability. This interme… ▽ More New HST/STIS optical spectra were obtained for a sample of early type stars with existing IUE UV spectra. These data were used to construct optical extinction curves whose general properties are discussed elsewhere. In this paper, we identify extinction features in the curves that are wider than diffuse interstellar bands (DIBs) but narrower than the well known broad band variability. This intermediate scale structure, or ISS, contains distinct features whose peaks can contribute a few percent to 20% of the total extinction. Most of the ISS variation can be captured by three principal components. We model the ISS with three Drude profiles and show that their strengths and widths vary from one sight line to another, but their central positions are stable, near 4370, 4870 and 6300A. The Very Broad Structure, VBS, in optical curves appears to be a minimum between the 4870 and 6300A absorption peaks. We find relations among the fit parameters and provide a physical interpretation of them in terms of a simplistic grain model. Finally, we note that the strengths of the 4370 and 4870A features are correlated to the strength of the 2175A UV bump, but that the 6300A feature is not, and that none of the ISS features are related to R(V). However, we verify that the broad band curvature of the continuous optical extinction is strongly related to R(V). △ Less

Submitted 29 January, 2020; originally announced January 2020.

Comments: 13 pages, 10 Figures, accepted by ApJ

An Analysis of the Shapes of Interstellar Extinction Curves. VII. Milky Way Spectrophotometric Optical-through-Ultraviolet Extinction and Its R-Dependence

Authors: Edward L. Fitzpatrick, Derck Massa, Karl D. Gordon, Ralph Bohlin, Geoffrey C. Clayton

Abstract: We produce a set of 72 NIR through UV extinction curves by combining new HST/STIS optical spectrophotometry with existing IUE spectrophotometry (yielding gapless coverage from 1150 to 10000 Angstroms) and NIR photometry. These curves are used to determine a new, internally consistent, NIR through UV Milky Way Mean Curve and to characterize how the shapes of the extinction curves depend on R(V). We… ▽ More We produce a set of 72 NIR through UV extinction curves by combining new HST/STIS optical spectrophotometry with existing IUE spectrophotometry (yielding gapless coverage from 1150 to 10000 Angstroms) and NIR photometry. These curves are used to determine a new, internally consistent, NIR through UV Milky Way Mean Curve and to characterize how the shapes of the extinction curves depend on R(V). We emphasize that while this dependence captures much of the curve variability, there remains considerable variation which is independent of R(V). We use the optical spectrophotometry to verify the presence of structure at intermediate wavelength scales in the curves. The fact that the optical through UV portions of the curves are sampled at relatively high resolution makes them very useful for determining how extinction affects different broad band systems, and we provide several examples. Finally, we compare our results to previous investigations. △ Less

Submitted 19 October, 2019; originally announced October 2019.

Comments: To appear in The Astrophysical Journal

Coordinated UV and X-ray spectroscopic observations of the O-type giant xi Per: the connection between X-rays and large-scale wind structure

Authors: Derck Massa, Lida Oskinova, Raman Prinja, Richard Ignace

Abstract: We present new, contemporaneous HST STIS and XMM observations of the O7 III(n)((f)) star xi Per. We supplement the new data with archival IUE spectra, to analyze the variability of the wind lines and X-ray flux of xi Pper. The variable wind of this star is known to have a 2.086 day periodicity. We use a simple, heuristic spot model which fits the low velocity (near surface) IUE wind line variabili… ▽ More We present new, contemporaneous HST STIS and XMM observations of the O7 III(n)((f)) star xi Per. We supplement the new data with archival IUE spectra, to analyze the variability of the wind lines and X-ray flux of xi Pper. The variable wind of this star is known to have a 2.086 day periodicity. We use a simple, heuristic spot model which fits the low velocity (near surface) IUE wind line variability very well, to demonstrate that the low velocity absorption in the new STIS spectra of N IV 1718 and Si IV 1402 vary with the same 2.086 day period. It is remarkable that the period and amplitude of the STIS data agree with those of the IUE spectra obtained 22 years earlier. We also show that the time variability of the new XMM fluxes are also consistent with the 2.086 day period. Thus, our new, multi-wavelength coordinated observations demonstrate that the mechanism which causes the UV wind line variability is also responsible for a significant fraction of the X-rays in single O stars. The sequence of events for the multi-wavelength light curve minima is: Si IV 1402, N IV 1718, and X-ray flux, each separated by a phase of about 0.06 relative to the 2.086 day period. Analysis of the X-ray fluxes shows that they become softer as they weaken. This is contrary to expectations if the variability is caused by periodic excess absorption. Further, the high resolution X-ray spectra suggest that the individual emission lines at maximum are more strongly blue shifted. If we interpret the low velocity wind line light curves in terms of our model, it implies that there are two bright regions, i.e., regions with less absorption, separated by 180 deg, on the surface of the star. We note that the presence and persistent of two spots separated by 180 deg suggests that a weak dipole magnetic field is responsible for the variability of the UV wind line absorption and X-ray flux in xi Per. △ Less

Submitted 28 January, 2019; originally announced January 2019.

Comments: 23 pages, 14 figures

A Study of H$α$ Line Profile Variations in $β$ Lyr

Authors: Richard Ignace, Sharon Gray, Macon Magno, Gary Henson, Derck Massa

Abstract: We examine over 160 archival H$α$ spectra from the Ritter Observatory for the interacting binary $β$~Lyr obtained between 1996 and 2000. The emission is characteristically double-peaked, but asymmetric, and with an absorption feature that is persistently blueshifted. Using a set of simplifying assumptions, phase varying emission line profiles are calculated for H$α$ formed entirely in a Keplerian… ▽ More We examine over 160 archival H$α$ spectra from the Ritter Observatory for the interacting binary $β$~Lyr obtained between 1996 and 2000. The emission is characteristically double-peaked, but asymmetric, and with an absorption feature that is persistently blueshifted. Using a set of simplifying assumptions, phase varying emission line profiles are calculated for H$α$ formed entirely in a Keplerian disk, and separately for the line formed entirely from an off-center bipolar flow. However, a dynamic spectrum of the data indicate the blueshifted feature is not always present, and the data are even suggestive of a drift of the feature in velocity shift. We explore whether a circumbinary envelope, hot spot on the accretion disk, or accretion stream could explain the observations. While none are satisfactory, an accretion stream explanation is somewhat promising. △ Less

Submitted 11 July, 2018; originally announced July 2018.

Comments: accepted to Astronomical Journal

Fast Flux Service Network Detection via Data Mining on Passive DNS Traffic

Authors: Pierangelo Lombardo, Salvatore Saeli, Federica Bisio, Davide Bernardi, Danilo Massa

Abstract: In the last decade, the use of fast flux technique has become established as a common practice to organise botnets in Fast Flux Service Networks (FFSNs), which are platforms able to sustain illegal online services with very high availability. In this paper, we report on an effective fast flux detection algorithm based on the passive analysis of the Domain Name System (DNS) traffic of a corporate n… ▽ More In the last decade, the use of fast flux technique has become established as a common practice to organise botnets in Fast Flux Service Networks (FFSNs), which are platforms able to sustain illegal online services with very high availability. In this paper, we report on an effective fast flux detection algorithm based on the passive analysis of the Domain Name System (DNS) traffic of a corporate network. The proposed method is based on the near-real-time identification of different metrics that measure a wide range of fast flux key features; the metrics are combined via a simple but effective mathematical and data mining approach. The proposed solution has been evaluated in a one-month experiment over an enterprise network, with the injection of pcaps associated with different malware campaigns, that leverage FFSNs and cover a wide variety of attack scenarios. An in-depth analysis of a list of fast flux domains confirmed the reliability of the metrics used in the proposed algorithm and allowed for the identification of many IPs that turned out to be part of two notorious FFSNs, namely Dark Cloud and SandiFlux, to the description of which we therefore contribute. All the fast flux domains were detected with a very low false positive rate; a comparison of performance indicators with previous works show a remarkable improvement. △ Less

Submitted 18 September, 2018; v1 submitted 17 April, 2018; originally announced April 2018.

Comments: This is a pre-print of an article published in the proceedings of 21st International Conference, ISC 2018, Guildford, UK, September 9-12, 2018. The final authenticated version is available online at: https://doi.org/10.1007/978-3-319-99136-8_25

MSC Class: 62-07; 62P30 ACM Class: C.2.0

Probing the Outflowing Multiphase Gas ~1 kpc Below the Galactic Center

Authors: Blair D. Savage, Tae-Sun Kim, Andrew J. Fox, Derck Massa, Rongmon Bordoloi, Edward. B. Jenkins, Nicolas Lehner, Joss Bland-Hawthorn, Felix J. Lockman, Svea Hernandez, Bart P. Wakker

Abstract: Comparison of ISM absorption in the UV spectrum of LS 4825, a B1 Ib-II star d=21+/-5 kpc from the Sun toward l = 1.67 deg and b = -6.63 deg, with ISM absorption toward an aligned foreground star at d < 7.0+/-1.7 kpc, allows us to isolate and study gas associated with the Milky Way nuclear wind. Spectra from the Space Telescope Imaging Spectrograph (STIS) show low ion absorption out to d < 7 kpc (… ▽ More Comparison of ISM absorption in the UV spectrum of LS 4825, a B1 Ib-II star d=21+/-5 kpc from the Sun toward l = 1.67 deg and b = -6.63 deg, with ISM absorption toward an aligned foreground star at d < 7.0+/-1.7 kpc, allows us to isolate and study gas associated with the Milky Way nuclear wind. Spectra from the Space Telescope Imaging Spectrograph (STIS) show low ion absorption out to d < 7 kpc ( e.g., O I, C II, Mg II, Si II, Fe II, S II) only between 0 and 40 km/s, while absorption at d > 7 kpc, ~1 kpc below the galactic plane, is complex and spans -290 to + 94 km/ s. The intermediate and high ions Si III, C IV, Si IV and N V show extremely strong absorption with multiple components from -283 to 107 km/ s implying that the ISM ~1 kpc below the galactic center has a substantial reservoir of plasma, and more gas containing C IV and N V than in the Carina OB1 association at z = 0 kpc. Abundances and physical conditions are presented for many absorption components. The high ion absorption traces cooling transition temperature plasma probably driven by the outflowing hot gas, while the extraordinary large thermal pressure, p/k ~ 10^5 cm^-3 K^-1 in an absorption component at -114 km/s probably arises from the ram pressure of the outflowing hot gas. The observations are consistent with a flow whose ionization structure in the high ions can be understood through a combination of non-equilibrium radiative cooling and turbulent mixing. △ Less

Submitted 1 November, 2017; v1 submitted 21 July, 2017; originally announced July 2017.

Comments: pdf file contains 56 pages, 14 tables and 16 figures. Manuscript accepted for 2017, ApJS, 232, 25

Journal ref: Savage, B. D., Kim, T-S, Fox, A. J., Massa, D., Bordoloi, R., Jenkins, E. B., Lehner, N., Bland-Hawthorn, J., Lockman, F. J, Hernandez, S., & Wakker, B. 2017, ApJS, 232, 25

Mass loss rates from mid-IR excesses in LMC and SMC O stars

Authors: D. Massa, A. W. Fullerton, R. K. Prinja

Abstract: We use a combination of BVJHK and Spitzer [3.6], [5.8] and [8.0] photometry to determine IR excesses for a sample of 58 LMC and 46 SMC O stars. This sample is ideal for determining IR excesses because the very small line of sight reddening minimizes uncertainties due to extinction corrections. We use the core-halo model developed by Lamers & Waters (1984a) to translate the excesses into mass loss… ▽ More We use a combination of BVJHK and Spitzer [3.6], [5.8] and [8.0] photometry to determine IR excesses for a sample of 58 LMC and 46 SMC O stars. This sample is ideal for determining IR excesses because the very small line of sight reddening minimizes uncertainties due to extinction corrections. We use the core-halo model developed by Lamers & Waters (1984a) to translate the excesses into mass loss rates and demonstrate that the results of this simple model agree with the more sophisticated CMFGEN models to within a factor of 2. Taken at face value, the derived mass loss rates are larger than those predicted by Vink et al. (2001), and the magnitude of the disagreement increases with decreasing luminosity. However, the IR excesses need not imply large mass loss rates. Instead, we argue that they probably indicate that the outer atmospheres of O stars contain complex structures and that their winds are launched with much smaller velocity gradients than normally assumed. If this is the case, it could affect the theoretical and observational interpretations of the "weak wind" problem, where classical mass loss indicators suggest that the mass loss rates of lower luminosity O stars are far less than expected. △ Less

Submitted 8 June, 2017; originally announced June 2017.

Comments: 15 pages, 10 figures. Accepted for publication in MNRAS

The highly ionized, high velocity gas in NGC 6231

Authors: Derck Massa

Abstract: It is well known that clusters of massive stars are influenced by the presence of strong winds, that they are sources of diffuse X-rays from shocked gas, and that this gas can be vented into the surrounding region or the halo through the champagne effect. However, the details of how these different environments interact and evolve are far from complete. This paper attributes the broad C IV 1500 ab… ▽ More It is well known that clusters of massive stars are influenced by the presence of strong winds, that they are sources of diffuse X-rays from shocked gas, and that this gas can be vented into the surrounding region or the halo through the champagne effect. However, the details of how these different environments interact and evolve are far from complete. This paper attributes the broad C IV 1500 absorption features (extending to -1900 \kms) that are seen in the spectra of main sequence B stars in NGC 6231 to gas in the cluster environment and not the B stars themselves. It is shown that the presence of a WC star, WR 79, in the cluster makes this gas detectable because its wind enriches the cluster gas with carbon. Given the available data, it is not clear whether the absorbing gas is simply the far wind of WR 79 or a collective cluster wind enriched by carbon from the wind of WR 79. If it is simply due to the wind, then this wind must flow, unimpeded for more than 2 pc, suggesting that the inner region of the cluster is nearly devoid of obstructing material. If it is actually a collective wind from the cluster, then we could be witnessing an important stage of galactic feedback. In either case, the observations provide a unique and significant piece to the puzzle of how massive, open clusters evolve. △ Less

Submitted 25 October, 2016; originally announced October 2016.

Comments: 9 pages, 6 figures

New Ultraviolet Extinction Curves for Interstellar Dust in M31

Authors: Geoffrey C. Clayton, Karl D. Gordon, Luciana C. Bianchi, Derck L. Massa, Edward L. Fitzpatrick, R. C. Bohlin, Michael J. Wolff

Abstract: New low-resolution UV spectra of a sample of reddened OB stars in M31 were obtained with HST/STIS to study the wavelength dependence of interstellar extinction and the nature of the underlying dust grain populations. Extinction curves were constructed for four reddened sightlines in M31 paired with closely matching stellar atmosphere models. The new curves have a much higher S/N than previous stud… ▽ More New low-resolution UV spectra of a sample of reddened OB stars in M31 were obtained with HST/STIS to study the wavelength dependence of interstellar extinction and the nature of the underlying dust grain populations. Extinction curves were constructed for four reddened sightlines in M31 paired with closely matching stellar atmosphere models. The new curves have a much higher S/N than previous studies. Direct measurements of N(H I) were made using the Ly$α$ absorption lines enabling gas-to-dust ratios to be calculated. The sightlines have a range in galactocentric distance of 5 to 14 kpc and represent dust from regions of different metallicities and gas-to-dust ratios. The metallicities sampled range from Solar to 1.5 Solar. The measured curves show similarity to those seen in the Milky Way and the Large Magellanic Cloud. The Maximum Entropy Method was used to investigate the dust composition and size distribution for the sightlines observed in this program finding that the extinction curves can be produced with the available carbon and silicon abundances if the metallicity is super-Solar. △ Less

Submitted 23 October, 2015; originally announced October 2015.

Comments: ApJ, in press, 9 pages, 5 figures

On the Origin of Wind Line Variability in O Stars

Authors: D. Massa, R. K. Prinja

Abstract: We analyze 10 UV time series for 5 stars which fulfill specific sampling and spectral criteria to constrain the origin of large-scale wind structure in O stars. We argue that excited state lines must arise close to the stellar surface and are an excellent diagnostic complement to resonance lines which, due to radiative transfer effects, rarely show variability at low velocity. Consequently, we spl… ▽ More We analyze 10 UV time series for 5 stars which fulfill specific sampling and spectral criteria to constrain the origin of large-scale wind structure in O stars. We argue that excited state lines must arise close to the stellar surface and are an excellent diagnostic complement to resonance lines which, due to radiative transfer effects, rarely show variability at low velocity. Consequently, we splice dynamic spectra of the excited state line, N IV1718, at low velocity to those of 1393 component of the Si IV 1400 doublet at high velocity in order to examine the temporal evolution of wind line features. These spliced time series reveal that nearly all of the features observed in the time series originate at or very near the stellar surface. Further, we positively identify the observational signature of equatorial co-rotating interaction regions in two of the five stars and possibly two others. In addition, we see no evidence for features originating further out in the wind. We use our results to consolidate the fact that the features seen in dynamic spectra must be huge, in order to remain in the line of sight for days, persisting to very large velocity and that the photospheric footprint of the features must also be quite large, ~ 15 - 20% of the stellar diameter. △ Less

Submitted 22 June, 2015; originally announced June 2015.

Comments: Accepted by ApJ

CIR Modulation of the X-ray Flux from the O7.5 III(n)((f)) Star xi Persei?

Authors: D. Massa, L. Oskinova, A. W. Fullerton, R. K. Prinja, D. A. Bohlender, N. D. Morrison, M. Blake, W. Pych

Abstract: 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… ▽ More 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. △ Less

Submitted 21 March, 2014; originally announced March 2014.

Comments: Accepted by MNRAS. 9 pages, 9 figures

UV diagnostic of porosity-free mass-loss estimates in B stars

Authors: Raman Prinja, Derck Massa

Abstract: We seek to establish evidence in UV P Cygni line profiles that the signs of wind clumping and porosity vary with velocity. We aim to demonstrate empirically that while at most wind velocities optically thick clumps cover only a fraction of the stellar surface, close to the terminal velocity where narrow absorption components (NACs) appear in UV lines the covering factor is approximately unity. SEI… ▽ More We seek to establish evidence in UV P Cygni line profiles that the signs of wind clumping and porosity vary with velocity. We aim to demonstrate empirically that while at most wind velocities optically thick clumps cover only a fraction of the stellar surface, close to the terminal velocity where narrow absorption components (NACs) appear in UV lines the covering factor is approximately unity. SEI line-synthesis models are used to determine the radial optical depths of blue and red components of the SiIV 1400 resonance line doublet in a sample of 12 B0 to B4 supergiants. We focus on stars with well developed NACs and relatively low terminal velocity so that the SiIV doublet components can be treated as radiatively decoupled and formed independently. For all 12 stars the mean optical depth ratio of the blue to red components is closer to ~ 2 (i.e. the ratio of oscillator strengths) in the NACs than at intermediate and lower velocities. The product of mass-loss rate and Si^3+ ion fraction calculated from the NAC optical depths is a factor of ~ 2 to 9 higher compared to mass-loss values sampled at ~ 0.4 to 0.6 of the terminal velocity. Since the wind effectively becomes `smooth' at the high NAC velocities and the column density is uniformly distributed over the stellar disk, the optical depths of the NACs are not seriously affected by porosity and this feature thus provides the most reliable measurement of mass-loss rate in the UV lines. Applications of this result to the weak-wind problem of late O-dwarf stars and the "PV mass loss discordance" in early O supergiants are discussed. △ Less

Submitted 9 September, 2013; originally announced September 2013.

Comments: Accepted for publication in Astronomy and Astrophysics; 7 pages; 3 figures

Looking deep into the Cat's Eye: Structure and rotation in the fast wind of the PN central star of NGC6543

Authors: R. K. Prinja, D. L. Massa, M. Cantiello

Abstract: We present HST/STIS time-series spectroscopy of the central star of the "Cat's Eye" planetary nebula NGC 6543. Intensive monitoring of the UV lines over a 5.8 hour period reveals well defined details of large-scale structure in the fast wind, which are exploited to provide new constraints on the rotation rate of the central star. We derive characteristics of the line profile variability that suppo… ▽ More We present HST/STIS time-series spectroscopy of the central star of the "Cat's Eye" planetary nebula NGC 6543. Intensive monitoring of the UV lines over a 5.8 hour period reveals well defined details of large-scale structure in the fast wind, which are exploited to provide new constraints on the rotation rate of the central star. We derive characteristics of the line profile variability that support a physical origin due to co-rotating interaction regions (CIRs) that are rooted at the stellar surface. The recurrence time of the observed spectral signatures of the CIRs is used to estimate the rotation period of the central star and, adopting a radius between 0.3 and 0.6 \Rsun constrains the rotational velocity to the range 54 \leq v_{rot} \leq 108 \kms. The implications of these results for single star evolution are discussed based on models calculated here for low-mass stars. Our models predict a sub-surface convective layer in NGC 6543 which we argue to be causally connected to the occurrence of structure in the fast wind. △ Less

Submitted 1 October, 2012; originally announced October 2012.

Comments: Accepted for publication in ApJ Letters

PN fast winds: Temporal structure and stellar rotation

Authors: R. K. Prinja, D. L. Massa, M. A. Urbaneja, R. -P. Kudritzki

Abstract: To diagnose the time-variable structure in the fast winds of central stars of planetary nebulae (CSPN), we present an analysis of P Cygni line profiles in FUSE satellite far-UV spectroscopic data. Archival spectra are retrieved to form time-series datasets for the H-rich CSPN NGC 6826, IC 418, IC 2149, IC 4593 and NGC 6543. Despite limitations due to the fragmented sampling of the time-series, we… ▽ More To diagnose the time-variable structure in the fast winds of central stars of planetary nebulae (CSPN), we present an analysis of P Cygni line profiles in FUSE satellite far-UV spectroscopic data. Archival spectra are retrieved to form time-series datasets for the H-rich CSPN NGC 6826, IC 418, IC 2149, IC 4593 and NGC 6543. Despite limitations due to the fragmented sampling of the time-series, we demonstrate that in all 5 CSPN the UV resonance lines are variable primarily due to the occurrence of blueward migrating discrete absorption components (DACs). Empirical (SEI) line-synthesis modelling is used to determine the range of fluctuations in radial optical depth, which are assigned to the temporal changes in large-scale wind structures. We argue that DACs are common in CSPN winds, and their empirical properties are akin to those of similar structures seen in the absorption troughs of massive OB stars. Constraints on PN central star rotation velocities are derived from Fast-Fourier Transform analysis of photospheric lines for our target stars. Favouring the causal role of co-rotating interaction regions, we explore connections between normalised DAC accelerations and rotation rates of PN central stars and O stars. The comparative properties suggest that the same physical mechanism is acting to generate large-scale structure in the line-driven winds in the two different settings. △ Less

Submitted 29 February, 2012; originally announced February 2012.

Comments: Accepted for publication in MNRAS; 10 pages, 5 figures

Characterizing the magnetic fields of the first tau Sco analogues

Authors: V. Petit, O. Kochukhov, D. L. Massa, W. L. F. Marcolino, G. A. Wade, R. Ignace

Abstract: The B0.2 V magnetic star tau Sco stands out from the larger population of massive OB stars due to its high X-ray activity, peculiar wind diagnostics and complex magnetic field. Recently, Petit et al. 2011 presented the discovery of the first two tau Sco analogues -- HD 66665 and HD 63425, identified by the striking similarity of their UV spectra to that of tau Sco. ESPaDOnS and Narval spectropolar… ▽ More The B0.2 V magnetic star tau Sco stands out from the larger population of massive OB stars due to its high X-ray activity, peculiar wind diagnostics and complex magnetic field. Recently, Petit et al. 2011 presented the discovery of the first two tau Sco analogues -- HD 66665 and HD 63425, identified by the striking similarity of their UV spectra to that of tau Sco. ESPaDOnS and Narval spectropolarimetric observations were obtained by the Magnetism in Massive Stars CFHT and TBL Large Programs, in order to characterize the stellar and magnetic properties of these stars. A magnetic field of similar surface strength was found on both stars, reinforcing the connection between the presence of a magnetic field and wind peculiarities. We present additional phase-resolved observations secured by the MiMeS collaboration for HD 66665 in order to measure its magnetic geometry, and correlate that geometry with diagnostics of mass-loss. △ Less

Submitted 7 November, 2011; v1 submitted 3 November, 2011; originally announced November 2011.

Comments: to appear in the proceedings of Stellar polarimetry: From birth to death, Madison, USA, June 27 - June 30, 2011 (version 2 was updated to correct a typo in the arXiv metadata)

The Cosmic Origins Spectrograph: On-Orbit Instrument Performance

Authors: S. Osterman, J. Green, C. Froning, S. Béland, E. Burgh, K. France, S. Penton, T. Delker, D. Ebbets, D. Sahnow, J. Bacinski, R. Kimble, J. Andrews, E. Wilkinson, J. McPhate, O. Siegmund, T. Ake, A. Aloisi, C. Biagetti, R. Diaz, W. Dixon, S. Friedman, P. Ghavamian, P. Goudfrooij, G. Hartig , et al. (9 additional authors not shown)

Abstract: The Cosmic Origins Spectrograph (COS) was installed in the Hubble Space Telescope in May, 2009 as part of Servicing Mission 4 to provide high sensitivity, medium and low resolution spectroscopy at far- and near-ultraviolet wavelengths (FUV, NUV). COS is the most sensitive FUV/NUV spectrograph flown to date, spanning the wavelength range from 900Å to 3200Å with peak effective area approaching 3000… ▽ More The Cosmic Origins Spectrograph (COS) was installed in the Hubble Space Telescope in May, 2009 as part of Servicing Mission 4 to provide high sensitivity, medium and low resolution spectroscopy at far- and near-ultraviolet wavelengths (FUV, NUV). COS is the most sensitive FUV/NUV spectrograph flown to date, spanning the wavelength range from 900Å to 3200Å with peak effective area approaching 3000 cm^2. This paper describes instrument design, the results of the Servicing Mission Orbital Verification (SMOV), and the ongoing performance monitoring program. △ Less

Submitted 28 December, 2010; originally announced December 2010.

Comments: 9 pages, 7 figures, submitted to Astrophysics and Space Science

Observing with HST below 1150Å: Extending the Cosmic Origins Spectrograph Coverage to 900Å

Authors: Steve Osterman, Steven V. Penton, Kevin France, Stéphane Béland, Stephan McCandliss, Jason McPhate, Derck Massa

Abstract: The far-ultraviolet (FUV) channel of the Cosmic Origins Spectrograph (COS) is designed to operate between 1130Å and 1850Å, limited at shorter wavelengths by the reflectivity of the MgF2 protected aluminum reflective surfaces on the Optical Telescope Assembly and on the COS FUV diffraction gratings. However, because the detector for the FUV channel is windowless, it was recognized early in the desi… ▽ More The far-ultraviolet (FUV) channel of the Cosmic Origins Spectrograph (COS) is designed to operate between 1130Å and 1850Å, limited at shorter wavelengths by the reflectivity of the MgF2 protected aluminum reflective surfaces on the Optical Telescope Assembly and on the COS FUV diffraction gratings. However, because the detector for the FUV channel is windowless, it was recognized early in the design phase that there was the possibility that COS would retain some sensitivity at shorter wavelengths due to the first surface reflection from the MgF2 coated optics. Preflight testing of the flight spare G140L grating revealed ~5% efficiency at 1066Å, and early on-orbit observations verified that the COS G140L/1230 mode was sensitive down to at least the Lyman limit with 10-20 cm^2 effective area between 912Å and 1070Å, and rising rapidly to over 1000 cm2 beyond 1150Å. Following this initial work we explored the possibility of using the G130M grating out of band to provide coverage down to 900Å. We present calibration results and ray trace simulations for these observing modes and explore additional configurations that have the potential to increase spectroscopic resolution, signal to noise, and observational efficiency below 1130Å. △ Less

Submitted 28 December, 2010; originally announced December 2010.

Comments: 6 pages, 5 figures. To appear in 2010 HST Calibration Workshop Proceedings

Discovery of the first tau Sco analogues: HD 66665 and HD 63425

Authors: V. Petit, D. L. Massa, W. L. F. Marcolino, G. A. Wade, R. Ignace, The MiMeS Collaboration

Abstract: The B0.2 V magnetic star tau Sco stands out from the larger population of massive OB stars due to its high X-ray activity, peculiar wind diagnostics and highly complex magnetic field. This paper presents the discovery of the first two tau Sco analogues - HD 66665 and HD 63425, identified by the striking similarity of their UV spectra to that of tau Sco. ESPaDOnS spectropolarimetric observations we… ▽ More The B0.2 V magnetic star tau Sco stands out from the larger population of massive OB stars due to its high X-ray activity, peculiar wind diagnostics and highly complex magnetic field. This paper presents the discovery of the first two tau Sco analogues - HD 66665 and HD 63425, identified by the striking similarity of their UV spectra to that of tau Sco. ESPaDOnS spectropolarimetric observations were secured by the Magnetism in Massive Stars CFHT Large Program, in order to characterize the stellar and magnetic properties of these stars. CMFGEN modelling of optical ESPaDOnS spectra and archived IUE UV spectra showed that these stars have stellar parameters similar to those of tau Sco. A magnetic field of similar surface strength is found on both stars, reinforcing the connection between the presence of a magnetic field and wind peculiarities. However, additional phase-resolved observations will be required in order to assess the potential complexity of the magnetic fields, and verify if the wind anomalies are linked to this property. △ Less

Submitted 20 December, 2010; originally announced December 2010.

Comments: 6 pages, 2 tables, 3 figures. Accepted for publication in MNRAS. The definitive version will be available at www.blackwel-synergy.com

Tau Sco: discovery of the clones

Authors: V. Petit, D. L. Massa, W. L. F. Marcolino, G. A. Wade, R. Ignace, the MiMeS Collaboration

Abstract: The B0.2 V magnetic star tau Sco stands out from the larger population of massive magnetic OB stars due to its high X-ray activity and remarkable wind, apparently related to its peculiar magnetic field - a field which is far more complex than the mostly-dipolar fields usually observed in magnetic OB stars. tau Sco is therefore a puzzling outlier in the larger picture of stellar magnetism - a star… ▽ More The B0.2 V magnetic star tau Sco stands out from the larger population of massive magnetic OB stars due to its high X-ray activity and remarkable wind, apparently related to its peculiar magnetic field - a field which is far more complex than the mostly-dipolar fields usually observed in magnetic OB stars. tau Sco is therefore a puzzling outlier in the larger picture of stellar magnetism - a star that still defies interpretation in terms of a physically coherent model. Recently, two early B-type stars were discovered as tau Sco analogues, identified by the striking similarity of their UV spectra to that of tau Sco, which was - until now - unique among OB stars. We present the recent detection of their magnetic fields by the MiMeS collaboration, reinforcing the connection between the presence of a magnetic field and wind anomalies (Petit et al. 2010). We will also present ongoing observational efforts undertaken to establish the precise magnetic topology, in order to provide additional constrains for existing models attempting to reproduce the unique wind structure of tau Sco-like stars. △ Less

Submitted 10 December, 2010; originally announced December 2010.

Comments: To appear in the proceedings of Magnetic Stars 2010

Mass loss rates inferred from mid-IR color excesses of LMC and SMC O stars

Authors: D. Massa, A. Fullerton, D. Lennon, R. K. Prinja

Abstract: We use a combination of VJHK and Spitzer} [3.6], [5.8] and [8.0] photometry, to determine IR excesses in a sample of LMC and SMC O stars. This sample is ideal for determining excesses because: 1) the distances to the stars, and hence their luminosities, are well-determined, and; 2) the very small line of sight reddenings minimize the uncertainties introduced by extinction corrections. We find IR e… ▽ More We use a combination of VJHK and Spitzer} [3.6], [5.8] and [8.0] photometry, to determine IR excesses in a sample of LMC and SMC O stars. This sample is ideal for determining excesses because: 1) the distances to the stars, and hence their luminosities, are well-determined, and; 2) the very small line of sight reddenings minimize the uncertainties introduced by extinction corrections. We find IR excesses much larger than expected from Vink et al. (2001) mass loss rates. This is in contrast to previous wind line analyses for many of the LMC stars which suggest mass loss rates much less than the Vink et al. predictions. ogether, these results indicate that the winds of the LMC and SMC O stars are strongly structured (clumped). △ Less

Submitted 16 October, 2010; originally announced October 2010.

Comments: 2 pages, 2 figures; to appear in the proceedings of IAU Symposium 272

tau Sco: The discovery of the clones

Authors: V. Petit, D. L. Massa, W. L. F. Marcolino, G. A. Wade, R. Ignace, the MiMeS Collaboration

Abstract: The B0.2 V magnetic star tau Sco stands out from the larger population of massive magnetic OB stars due to its remarkable, superionized wind, apparently related to its peculiar magnetic field - a field which is far more complex than the mostly-dipolar fields usually observed in magnetic OB stars. tau Sco is therefore a puzzling outlier in the larger picture of stellar magnetism - a star that still… ▽ More The B0.2 V magnetic star tau Sco stands out from the larger population of massive magnetic OB stars due to its remarkable, superionized wind, apparently related to its peculiar magnetic field - a field which is far more complex than the mostly-dipolar fields usually observed in magnetic OB stars. tau Sco is therefore a puzzling outlier in the larger picture of stellar magnetism - a star that still defies interpretation in terms of a physically coherent model. Recently, two early B-type stars were discovered as tau Sco analogues, identified by the striking similarity of their UV spectra to that of tau Sco, which was - until now - unique among OB stars. We present the recent detection of their magnetic fields by the MiMeS collaboration, reinforcing the connection between the presence of a magnetic field and a superionized wind. We will also present ongoing observational efforts undertaken to establish the precise magnetic topology, in order to provide additional constrains for existing models attempting to reproduce the unique wind structure of tau Sco-like stars. △ Less

Submitted 11 October, 2010; originally announced October 2010.

Comments: Proceedings of IAUS272: Active OB stars

Infrared properties of Active OB stars in the Magellanic Clouds from the Spitzer SAGE Survey

Authors: A. Z. Bonanos, D. J. Lennon, D. L. Massa, M. Sewilo, F. Köhlinger, N. Panagia, J. Th. van Loon, C. J. Evans, L. J. Smith, M. Meixner, K. Gordon, the SAGE teams

Abstract: We present a study of the infrared properties of 4922 spectroscopically confirmed massive stars in the Large and Small Magellanic Clouds, focusing on the active OB star population. Besides OB stars, our sample includes yellow and red supergiants, Wolf-Rayet stars, Luminous Blue Variables (LBVs) and supergiant B[e] stars. We detect a distinct Be star sequence, displaced to the red, and find a highe… ▽ More We present a study of the infrared properties of 4922 spectroscopically confirmed massive stars in the Large and Small Magellanic Clouds, focusing on the active OB star population. Besides OB stars, our sample includes yellow and red supergiants, Wolf-Rayet stars, Luminous Blue Variables (LBVs) and supergiant B[e] stars. We detect a distinct Be star sequence, displaced to the red, and find a higher fraction of Oe and Be stars among O and early-B stars in the SMC, respectively, when compared to the LMC, and that the SMC Be stars occur at higher luminosities. We also find photometric variability among the active OB population and evidence for transitions of Be stars to B stars and vice versa. We furthermore confirm the presence of dust around all the supergiant B[e] stars in our sample, finding the shape of their spectral energy distributions (SEDs) to be very similar, in contrast to the variety of SED shapes among the spectrally variable LBVs. △ Less

Submitted 1 October, 2010; originally announced October 2010.

Comments: 5 pages, 1 figure, to appear in the proceedings of the IAUS 272 on "Active OB stars: structure, evolution, mass loss and critical limits" (Paris, July 19-23, 2010), Cambridge University Press. Editors C. Neiner, G. Wade, G. Meynet and G. Peters

The Infrared Properties of Massive Stars in the Magellanic Clouds

Authors: A. Z. Bonanos, D. J. Lennon, D. L. Massa, M. Sewilo, F. Köhlinger, N. Panagia, J. Th. van Loon, C. J. Evans, L. J. Smith, M. Meixner, K. Gordon, the SAGE teams

Abstract: We present results of our study of the infrared properties of massive stars in the Large and Small Magellanic Clouds, which are based on the Spitzer SAGE surveys of these galaxies. We have compiled catalogs of spectroscopically confirmed massive stars in each galaxy, as well as photometric catalogs for a subset of these stars that have infrared counterparts in the SAGE database, with uniform photo… ▽ More We present results of our study of the infrared properties of massive stars in the Large and Small Magellanic Clouds, which are based on the Spitzer SAGE surveys of these galaxies. We have compiled catalogs of spectroscopically confirmed massive stars in each galaxy, as well as photometric catalogs for a subset of these stars that have infrared counterparts in the SAGE database, with uniform photometry from 0.3 to 24 microns in the UBVIJHKs+IRAC+MIPS24 bands. These catalogs enable a comparative study of infrared excesses of OB stars, classical Be stars, yellow and red supergiants, Wolf-Rayet stars, Luminous Blue Variables and supergiant B[e] stars, as a function of metallicity, and provide the first roadmaps for interpreting luminous, massive, resolved stellar populations in nearby galaxies at infrared wavelengths. △ Less

Submitted 1 October, 2010; originally announced October 2010.

Comments: 5 pages, 1 figure, contribution to 39th Liege Astrophysical Colloquium, 'The multi-wavelength view of hot, massive stars'

Signature of wide-spread clumping in B supergiant winds

Authors: Raman Prinja, Derck Massa

Abstract: We seek to establish additional observational signatures of the effects of clumping in OB star winds. The action of clumping on strategic wind-formed spectral lines is tested to steer the development of models for clumped winds and thus improve the reliability of mass-loss determinations for massive stars.The SiIV 1400 resonance line doublets of B0 to B5 supergiants are analysed using empirical li… ▽ More We seek to establish additional observational signatures of the effects of clumping in OB star winds. The action of clumping on strategic wind-formed spectral lines is tested to steer the development of models for clumped winds and thus improve the reliability of mass-loss determinations for massive stars.The SiIV 1400 resonance line doublets of B0 to B5 supergiants are analysed using empirical line-synthesis models. The focus is on decoding information on wind clumping from measurements of ratios of the radial optical depths (tau_(rad)(w)) of the red and blue components of the SiIV doublet. We exploit in particular the fact that the two doublet components are decoupled and formed independently for targets with relatively low wind terminal velocities. Line-synthesis analyses reveal that the mean ratio of tau_(rad)(w) of the blue to red SiIV components are rarely close to the canonical value of ~ 2 (expected from atomic constants), and spread instead over a range of values between ~1 and 2. These results are interpreted in terms of a photosphere that is partially obscured by optically thick structures in the outflowing gas.The spectroscopic signatures established in this study demonstrate the wide-spread existence of wind clumping in B supergiants. The additional information in unsaturated doublet profiles provides a means to quantify the porosity of the winds. △ Less

Submitted 16 July, 2010; originally announced July 2010.

Comments: Accepted for publication in A&A Letters

A massive runaway star from 30 Doradus

Authors: C. J. Evans, N. R. Walborn, P. A. Crowther, V. Henault-Brunet, D. Massa, W. D. Taylor, I. D. Howarth, H. Sana, D. J. Lennon, J. Th. van Loon

Abstract: We present the first ultraviolet (UV) and multi-epoch optical spectroscopy of 30 Dor 016, a massive O2-type star on the periphery of 30 Doradus in the Large Magellanic Cloud. The UV data were obtained with the Cosmic Origins Spectrograph on the Hubble Space Telescope as part of the Servicing Mission Observatory Verification program after Servicing Mission 4, and reveal #016 to have one of the fas… ▽ More We present the first ultraviolet (UV) and multi-epoch optical spectroscopy of 30 Dor 016, a massive O2-type star on the periphery of 30 Doradus in the Large Magellanic Cloud. The UV data were obtained with the Cosmic Origins Spectrograph on the Hubble Space Telescope as part of the Servicing Mission Observatory Verification program after Servicing Mission 4, and reveal #016 to have one of the fastest stellar winds known. From analysis of the CIV 1548-51 doublet we find a terminal velocity, v_infty=3450 +/- 50km/s. Optical spectroscopy is from the VLT-FLAMES Tarantula Survey, from which we rule out a massive companion (with 2d<P<1yr) to a confidence of 98%. The radial velocity of #016 is offset from the systemic value by -85km/s, suggesting that the star has traveled the 120pc from the core of 30 Doradus as a runaway, ejected via dynamical interactions. △ Less

Submitted 29 April, 2010; originally announced April 2010.

Comments: 6 pages, ApJL, in press

Spitzer SAGE-SMC Infrared Photometry of Massive Stars in the Small Magellanic Cloud

Authors: A. Z. Bonanos, D. J. Lennon, F. Köhlinger, J. Th. van Loon, D. L. Massa, M. Sewilo, C. J. Evans, N. Panagia, B. L. Babler, M. Block, S. Bracker, C. W. Engelbracht, K. D. Gordon, J. L. Hora, R. Indebetouw, M. R. Meade, M. Meixner, K. A. Misselt, T. P. Robitaille, B. Shiao, B. A. Whitney

Abstract: We present a catalog of 5324 massive stars in the Small Magellanic Cloud (SMC), with accurate spectral types compiled from the literature, and a photometric catalog for a subset of 3654 of these stars, with the goal of exploring their infrared properties. The photometric catalog consists of stars with infrared counterparts in the Spitzer, SAGE-SMC survey database, for which we present uniform phot… ▽ More We present a catalog of 5324 massive stars in the Small Magellanic Cloud (SMC), with accurate spectral types compiled from the literature, and a photometric catalog for a subset of 3654 of these stars, with the goal of exploring their infrared properties. The photometric catalog consists of stars with infrared counterparts in the Spitzer, SAGE-SMC survey database, for which we present uniform photometry from 0.3-24 um in the UBVIJHKs+IRAC+MIPS24 bands. We compare the color magnitude diagrams and color-color diagrams to those of the Large Magellanic Cloud (LMC), finding that the brightest infrared sources in the SMC are also the red supergiants, supergiant B[e] (sgB[e]) stars, luminous blue variables, and Wolf-Rayet stars, with the latter exhibiting less infrared excess, the red supergiants being less dusty and the sgB[e] stars being on average less luminous. Among the objects detected at 24 um are a few very luminous hypergiants, 4 B-type stars with peculiar, flat spectral energy distributions, and all 3 known luminous blue variables. We detect a distinct Be star sequence, displaced to the red, and suggest a novel method of confirming Be star candidates photometrically. We find a higher fraction of Oe and Be stars among O and early-B stars in the SMC, respectively, when compared to the LMC, and that the SMC Be stars occur at higher luminosities. We estimate mass-loss rates for the red supergiants, confirming the correlation with luminosity even at the metallicity of the SMC. Finally, we confirm the new class of stars displaying composite A & F type spectra, the sgB[e] nature of 2dFS1804 and find the F0 supergiant 2dFS3528 to be a candidate luminous blue variable with cold dust. △ Less

Submitted 4 July, 2010; v1 submitted 6 April, 2010; originally announced April 2010.

Comments: 23 pages, 17 figures, 5 tables, accepted for publication in the Astronomical Journal

Journal ref: Astron.J.140:416-429,2010

A FUSE View of the Stellar Winds of Planetary Nebula Central Stars

Authors: Martin A. Guerrero, Gerardo Ramos-Larios, Derck Massa

Abstract: Since the IUE satellite produced a vast collection of high-resolution UV spectra of central stars of planetary nebulae (CSPNe), there has not been any further systematic study of the stellar winds of these stars. The high spectral resolution, sensitivity and large number of archival observations in the FUSE archive allow the study of the stellar winds of CSPNe in the far UV domain where lines of… ▽ More Since the IUE satellite produced a vast collection of high-resolution UV spectra of central stars of planetary nebulae (CSPNe), there has not been any further systematic study of the stellar winds of these stars. The high spectral resolution, sensitivity and large number of archival observations in the FUSE archive allow the study of the stellar winds of CSPNe in the far UV domain where lines of species spanning a wide excitation range can be observed. We present here a preliminary analysis of the P Cygni profiles of a sample of 60 CSPNe observed by FUSE. P Cygni profiles evidencing fast stellar winds with velocities between 200 and 4,300 km/s have been found in 40 CSPNe. In many cases, this is the first time that fast stellar winds have been reported for these PNe. A detailed study of these far-UV spectra is on-going. △ Less

Submitted 1 December, 2009; originally announced December 2009.

Comments: 4 pages, 3 figures; accepted, to appear in Publications of the Astronomical Society of Australia

The dense and asymmetric central star wind of the young PN He 2-138

Authors: R. K. Prinja, S. E. Hodges, M. A. Urbaneja, D. L. Massa

Abstract: We present optical ESO time-series and UV archival (FUSE, HST, IUE) spectroscopy of the H-rich central star of He 2-138. Our study targets the central star wind in a very young planetary nebula, and explores physical conditions that may provide clues to the nature of the preceding post-AGB super-wind phases of the star. We provide evidence for a dense, slowly accelerating outflow that is variabl… ▽ More We present optical ESO time-series and UV archival (FUSE, HST, IUE) spectroscopy of the H-rich central star of He 2-138. Our study targets the central star wind in a very young planetary nebula, and explores physical conditions that may provide clues to the nature of the preceding post-AGB super-wind phases of the star. We provide evidence for a dense, slowly accelerating outflow that is variable on time-scales of hours. Line-synthesis modelling (SEI and CMFGEN) of low and high ionization UV and optical lines is interpreted in terms of an asymmetric, two-component outflow, where high-speed high-ionization gas forms mostly in the polar region. Slower, low ionization material is then confined primarily to a cooler equatorial component of the outflow. A dichotomy is also evident at photospheric levels. We also document temporal changes in the weak photospheric lines of He 2-138, with tentative evidence for a 0.36-day modulation in blue-to-red migrating features in the absorption lines. These structures may betray 'wave-leakage' of prograde non-radial pulsations of the central star. These multi-waveband results on the aspherical outflow of He 2-138 are discussed in the context of current interest in understanding the origin of axi- and point-symmetric planetary nebulae. △ Less

Submitted 26 October, 2009; originally announced October 2009.

Comments: 10 pages, 9 figures, MNRAS in press

Spitzer SAGE Infrared Photometry of Massive Stars in the Large Magellanic Cloud

Authors: A. Z. Bonanos, D. L. Massa, M. Sewilo, D. J. Lennon, N. Panagia, L. J. Smith, M. Meixner, B. L. Babler, S. Bracker, M. R. Meade, K. D. Gordon, J. L. Hora, R. Indebetouw, B. A. Whitney

Abstract: We present a catalog of 1750 massive stars in the Large Magellanic Cloud, with accurate spectral types compiled from the literature, and a photometric catalog for a subset of 1268 of these stars, with the goal of exploring their infrared properties. The photometric catalog consists of stars with infrared counterparts in the Spitzer SAGE survey database, for which we present uniform photometry fr… ▽ More We present a catalog of 1750 massive stars in the Large Magellanic Cloud, with accurate spectral types compiled from the literature, and a photometric catalog for a subset of 1268 of these stars, with the goal of exploring their infrared properties. The photometric catalog consists of stars with infrared counterparts in the Spitzer SAGE survey database, for which we present uniform photometry from 0.3-24 microns in the UBVIJHKs+IRAC+MIPS24 bands. The resulting infrared color-magnitude diagrams illustrate that the supergiant B[e], red supergiant and luminous blue variable (LBV) stars are among the brightest infrared point sources in the Large Magellanic Cloud, due to their intrinsic brightness, and at longer wavelengths, due to dust. We detect infrared excesses due to free-free emission among ~900 OB stars, which correlate with luminosity class. We confirm the presence of dust around 10 supergiant B[e] stars, finding the shape of their spectral energy distributions (SEDs) to be very similar, in contrast to the variety of SED shapes among the spectrally variable LBVs. The similar luminosities of B[e] supergiants (log L/Lo>=4) and the rare, dusty progenitors of the new class of optical transients (e.g. SN 2008S and NGC 300 OT), plus the fact that dust is present in both types of objects, suggests a common origin for them. We find the infrared colors for Wolf-Rayet stars to be independent of spectral type and their SEDs to be flatter than what models predict. The results of this study provide the first comprehensive roadmap for interpreting luminous, massive, resolved stellar populations in nearby galaxies at infrared wavelengths. △ Less

Submitted 11 August, 2009; v1 submitted 8 May, 2009; originally announced May 2009.

Comments: 57 pages, 19 figures, 5 tables, accepted for publication in the Astronomical Journal (v3 corrects typos)

Journal ref: Astron.J.138:1003-1021,2009

An Analysis of the Shapes of Interstellar Extinction Curves. VI. The Near-IR Extinction Law

Authors: Edward L. Fitzpatrick, Derck Massa

Abstract: We combine new HST/ACS observations and existing data to investigate the wavelength dependence of NIR extinction. Previous studies suggest a power-law form, with a "universal" value of the exponent, although some recent observations indicate that significant sight line-to-sight line variability may exist. We show that a power-law model provides an excellent fit to most NIR extinction curves, but… ▽ More We combine new HST/ACS observations and existing data to investigate the wavelength dependence of NIR extinction. Previous studies suggest a power-law form, with a "universal" value of the exponent, although some recent observations indicate that significant sight line-to-sight line variability may exist. We show that a power-law model provides an excellent fit to most NIR extinction curves, but that the value of the power, beta, varies significantly from sight line-to-sight line. Therefore, it seems that a "universal NIR extinction law" is not possible. Instead, we find that as beta decreases, R(V) [=A(V)/E(B-V)] tends to increase, suggesting that NIR extinction curves which have been considered "peculiar" may, in fact, be typical for different R(V) values. We show that the power law parameters can depend on the wavelength interval used to derive them, with the beta increasing as longer wavelengths are included. This result implies that extrapolating power law fits to determine R(V) is unreliable. To avoid this problem, we adopt a different functional form for NIR extinction. This new form mimics a power law whose exponent increases with wavelength, has only 2 free parameters, can fit all of our curves over a longer wavelength baseline and to higher precision, and produces R(V) values which are consistent with independent estimates and commonly used methods for estimating R(V). Furthermore, unlike the power law model, it gives R(V)'s that are independent of the wavelength interval used to derive them. It also suggests that the relation R(V) = -1.36 E(K-V)/E(B-V) - 0.79 can estimate R(V) to +/-0.12. Finally, we use model extinction curves to show that our extinction curves are in accord with theoretical expectations. △ Less

Submitted 1 May, 2009; originally announced May 2009.

Comments: To appear in the Astrophysical Journal

Journal ref: Astrophys.J.699:1209-1222,2009

Massive Star Multiplicity: The Cepheid W Sgr

Authors: Nancy Remage Evans, Derck Massa, Charles Proffitt

Abstract: We have obtained spectra of the W Sgr system with the STIS spectrograph on the Hubble Space Telescope. The spectra resolve the system into a distant companion B which is the hottest star in the system and the spectroscopic binary (A = Aa + Ab). A and B are separated by 0.16". We have extracted the spectra of both of these. We see no flux in the Aa + Ab spectrum which cannot be accounted for by t… ▽ More We have obtained spectra of the W Sgr system with the STIS spectrograph on the Hubble Space Telescope. The spectra resolve the system into a distant companion B which is the hottest star in the system and the spectroscopic binary (A = Aa + Ab). A and B are separated by 0.16". We have extracted the spectra of both of these. We see no flux in the Aa + Ab spectrum which cannot be accounted for by the Cepheid, and put an upper limit on the spectral type and mass of the companion Ab of F5 V and $\leq$1.4Msun. Using the orbit from HST FGS measurements from Benedict, et al., this results in an upper limit to the mass of the Cepheid of $\leq$5.4Msun. We also discuss two possible distant companions. Based on photometry from the 2MASS Point Source Catalog, they are not physical companions of the W Sgr system. △ Less

Submitted 18 February, 2009; originally announced February 2009.

Comments: accepted for AJ

Quantitative studies of the optical and UV spectra of Galactic early B supergiants I. Fundamental parameters

Authors: Samantha C. Searle, Raman K. Prinja, Derck Massa, Robert Ryans

Abstract: We undertake an optical and ultraviolet spectroscopic analysis of a sample of 20 Galactic B0 - B5 supergiants of luminosity classes Ia, Ib, Iab and II. Fundamental stellar parameters are obtained from optical diagnostics and a critical comparison of the model predictions to observed UV spectral features is made. These parameters are derived for individual stars using CMFGEN, the nLTE, line-blank… ▽ More We undertake an optical and ultraviolet spectroscopic analysis of a sample of 20 Galactic B0 - B5 supergiants of luminosity classes Ia, Ib, Iab and II. Fundamental stellar parameters are obtained from optical diagnostics and a critical comparison of the model predictions to observed UV spectral features is made. These parameters are derived for individual stars using CMFGEN, the nLTE, line-blanketed model atmosphere code of Hillier et al., 1998. The B supergiant temperature scale derived here shows a reduction of 1000 - 3000 K compared to previous results obtained using unblanketed codes. Mass loss rate estimates are in good agreement with predicted theoretical values and all of the 20 B0 - B5 supergiants analysed show evidence for CNO processing. The observed WLR values calculated for B0 - B0.7 supergiants are larger than predicted values, whereas the reverse is true for B1 - B5 supergiants. This means that the discrepancy between observed and theoretical values cannot be resolved by adopting clumped (i.e., lower) mass loss rates, as for O stars. The most surprising result is that, although CMFGEN succeeds in reproducing the optical spectrum accurately, it fails to reproduce key UV diagnostics, such as NV and CIV P Cygni profiles, precisely. This problem arises because the models are not ionised enough and fail to reproduce the full extent of the observed absorption trough of the P Cygni profiles. These findings add further support to the need to revise the standard model of massive star winds. △ Less

Submitted 28 January, 2008; originally announced January 2008.

Comments: A&A in press; 35 pages, 33 figures

New insights into the nature of the SMC WR/LBV binary HD 5980

Authors: C. Foellmi, G. Koenigsberger, L. Georgiev, O. Toledano, S. V. Marchenko, P. Massey, T. H. Dall, A. F. J. Moffat, N. Morrell, M. Corcoran, A. Kaufer, Y. Naze, J. Pittard, N. St. -Louis, A. Fullerton, D. Massa, A. M. T. Pollock

Abstract: 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 primaril… ▽ More 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. △ Less

Submitted 29 November, 2007; originally announced November 2007.

Comments: 48 pages, 26 figures

Journal ref: Rev.Mex.Astron.Astrofis.44:3-27,2008

Structure and clumping in the fast wind of NGC6543

Authors: R. K. Prinja, S. E. Hodges, D. L. Massa, A. W. Fullerton, A. W. Burnley

Abstract: Far-UV spectroscopy from the FUSE satellite is analysed to uniquely probe spatial structure and clumping in the fast wind of the central star of the H-rich planetary nebula NGC6543 (HD164963). Time-series data of the unsaturated PV 1118, 1128 resonance line P Cygni profiles provide a very sensitive diagnostic of variable wind conditions in the outflow. We report on the discovery of episodic and… ▽ More Far-UV spectroscopy from the FUSE satellite is analysed to uniquely probe spatial structure and clumping in the fast wind of the central star of the H-rich planetary nebula NGC6543 (HD164963). Time-series data of the unsaturated PV 1118, 1128 resonance line P Cygni profiles provide a very sensitive diagnostic of variable wind conditions in the outflow. We report on the discovery of episodic and recurrent optical depth enhancements in the PV absorption troughs, with some evidence for a 0.17-day modulation time-scale. SEI line-synthesis modelling is used to derive physical properties, including the optical depth evolution of individual `events'. The characteristics of these features are essentially identical to the `discrete absorption components' (DACs) commonly seen in the UV lines of massive OB stars. We have also employed the unified model atmosphere code CMFGEN to explore spectroscopic signatures of clumping, and report in particular on the clear sensitivity of the PV lines to the clump volume filling factor. The results presented here have implications for the downward revision of mass-loss rates in PN central stars. We conclude that the temporal structures seen in the PV lines of NGC6543 likely have a physical origin that is similar to that operating in massive, luminous stars, and may be related to near-surface perturbations caused by stellar pulsation and/or magnetic fields. △ Less

Submitted 17 August, 2007; originally announced August 2007.

Comments: 11 pages, 11 figures. Accepted for publication in MNRAS

The effects of clumping on wind line variability

Authors: D. Massa, R. K. Prinja, A. W. Fullerton

Abstract: We review the effects of clumping on the profiles of resonance doublets. By allowing the ratio of the doublet oscillator strenghts to be a free parameter, we demonstrate that doublet profiles contain more information than is normally utilized. In clumped (or porous) winds, this ratio can lies between unity and the ratio of the f-values, and can change as a function of velocity and time, dependin… ▽ More We review the effects of clumping on the profiles of resonance doublets. By allowing the ratio of the doublet oscillator strenghts to be a free parameter, we demonstrate that doublet profiles contain more information than is normally utilized. In clumped (or porous) winds, this ratio can lies between unity and the ratio of the f-values, and can change as a function of velocity and time, depending on the fraction of the stellar disk that is covered by material moving at a particular velocity at a given moment. Using these insights, we present the results of SEI modeling of a sample of B supergiants, zeta Pup and a time series for a star whose terminal velocity is low enough to make the components of its Si IV 1400 doublet independent. These results are interpreted within the framework of the Oskinova et al. (2007) model, and demonstrate how the doublet profiles can be used to extract infromation about wind structure. △ Less

Submitted 6 August, 2007; originally announced August 2007.

Comments: 3 pages, to appear in Clumping in Hot Star Winds, W.-R. Hamann, A. Feldmeier & L. Oskinova, eds., Potsdam: Univ.-Verl., 2007, URN: http://nbn-resolving.de/urn:nbn:de:kobv:517-opus-13981