The influence of macroclumps for free–free spectral energy distributions (SEDs) of ionized winds is considered. The goal is to emphasize distinctions between microclumping and macroclumping effects. Microclumping can alter SED slopes and flux levels if the volume filling factor of the clumps varies with radius; however, the modifications are independent of the clump geometry. To what extent does macroclumping alter SED slopes and flux levels? In addressing the question, two specific types of macroclump geometries are explored: shell fragments (pancake-shaped) and spherical clumps. Analytic and semi-analytic results are derived in the limiting case that clumps never obscure one another. Numerical calculations based on a porosity formalism is used when clumps do overlap. Under the assumptions of a constant expansion, isothermal, and fixed ionization wind, the fragment model leads to results that are essentially identical to the microclumping result. Mass-loss rate determinations are not affected by porosity effects for shell fragments. By contrast, spherical clumps can lead to a reduction in long-wavelength fluxes, but the reductions are only significant for extreme volume filling factors.

In this paper we investigate the multiwavelengths properties of the magnetic early B-type star HR7355. We present its radio light curves at several frequencies, taken with the Jansky Very Large Array, and X-ray spectra, taken with the XMM X-ray telescope. Modeling of the radio light curves for the Stokes I and V provides a quantitative analysis of the HR7355 magnetosphere. A comparison between HR7355 and a similar analysis for the Ap star CUVir, allows us to study how the different physical parameters of the two stars affect the structure of the respective magnetospheres where the non-thermal electrons originate. Our model includes a cold thermal plasma component that accumulates at high magnetic latitudes that influences the radio regime, but does not give rise to X-ray emission. Instead, the thermal X-ray emission arises from shocks generated by wind stream collisions close to the magnetic equatorial plane. The analysis of the X-ray spectrum of HR7355 also suggests the presence of a non-thermal radiation. Comparison between the spectral index of the power-law X-ray energy distribution with the non-thermal electron energy distribution indicates that the non-thermal X-ray component could be the auroral signature of the non-thermal electrons that impact the stellar surface, the same non-thermal electrons that are responsible for the observed radio emission. On the basis of our analysis, we suggest a novel model that simultaneously explains the X-ray and the radio features of HR7355 and is likely relevant for magnetospheres of other magnetic early type stars.

The β Cephei stars represent an important class of massive star pulsators that probe the evolution of B-type stars and the transition from main sequence to hydrogen-shell burning evolution. By understanding β Cep stars, we gain insights into the detailed physics of massive star evolution, including rotational mixing, convective core overshooting, magnetic fields, and stellar winds, all of which play important roles. Similarly, modeling their pulsation provides additional information into their interior structures. Furthermore, measurements of the rate of change of pulsation period offer a direct measure of β Cephei stellar evolution. In this work, we compute state-of-the-art stellar evolution models assuming different amounts of initial rotation and convective core overshoot and measure the theoretical rates of period change, that we compare to rates previously measured for a sample of β Cephei stars. The results of this comparison are mixed. For three stars, the rates are too low to infer any information from stellar evolution models, whereas for three other stars the rates are too high. We infer stellar parameters, such as mass and age, for two β Cephei stars: ξ1 CMa and δ Cet, which agree well with independent measurements. We explore ideas for why models may not predict the higher rates of period change. In particular, period drifts in β Cep stars can artificially lead to overestimated rates of secular period change.

Two equations in Ignace et al. (2019) have minor typos. The first correction is in Eq. (15). That expression has three lines. In the first line, there should be no multiplicative factor of R0 appearing. The correct equation is: (Formula Presented) The other two lines in Eq. (15) are correct. The second correction is in Eq. (19). As in Eq. (15), the appearance of R0 is in error. The correct expression is (Formula Presented) The mistakes appearing in Eqs. (15) and (19) are typos that have no bearing on the results appearing in the figures which were computed correctly, nor on the discussion or conclusions of the paper. (Equation Presented).

Within this past decade, advances in spectropolarimetric analyses allowed mapping surface characteristics of nearby non-resolved stars with unique characteristics. In this study, we attempt mapping of the magnetic field structure and surface distribution of elements for such a star, the magnetic Ap star phi Dra. Using high-precision photometric data from TESS satellite, we improve its rotation period. Spectrum synthesis calculations allow compilation of a list of lines present in its spectrum. The resulting synthetic spectrum and observed NARVAL spectra are used to re-estimate element abundances. Least-squares deconvolved (LSD) intensity and circular polarisation profiles are then computed from a sample of 1260 metal lines. To determine which element(s) are most suitable for in-depth analysis, variability of LSD profiles is studied qualitatively for Fe, Cr and Si. The longitudinal magnetic field of phi Dra is calculated from LSD circular polarisation profiles. Stellar magnetic field maps and distributions of Fe concentration are derived through Zeeman Doppler Imaging (ZDI). The resulting maps of this study show five areas of high Fe concentrations, in the Northern stellar hemisphere. The magnetic field topology of phi Dra resulting from our analysis is that of an offset dipole with small quadrupole contributions. Our abundance and magnetic maps suggest correlation between high concentrations of Fe and high magnetic field strength. The field is primarily radial in 4 out of 5 such regions, contrary to theoretical expectations.

We present new ATCA multiwavelength radio measurements (range 2.1–21.2 GHz) of the early-type magnetic star ρ Oph A, performed in 2019 March during three different observing sessions. These new ATCA observations evidence a clear rotational modulation of the stellar radio emission and the detection of coherent auroral radio emission from ρ Oph A at 2.1 GHz. We collected high-resolution optical spectra of ρ Oph A acquired by several instruments over a time span of about 10 yr. We also report new magnetic field measurements of ρ Oph A that, together with the radio light curves and the temporal variation of the equivalent width of the He I line (λ = 5015 Å), were used to constrain the rotation period and the stellar magnetic field geometry. The above results have been used to model the stellar radio emission, modelling that allowed us to constrain the physical condition of ρ Oph A’s magnetosphere. Past XMM–Newton measurements showed periodic X-ray pulses from ρ Oph A. We correlate the X-ray light curve with the magnetic field geometry of ρ Oph A. The already published XMM–Newton data have been re-analysed showing that the X-ray spectra of ρ Oph A are compatible with the presence of a non-thermal X-ray component. We discuss a scenario where the emission phenomena occurring at the extremes of the electromagnetic spectrum, radio and X-ray, are directly induced by the same plasma process. We interpret the observed X-ray and radio features of ρ Oph A as having an auroral origin.

The non-thermal radio emission of main-sequence early-type stars is a signature of stellar magnetism. We present multiwavelength (1.6-16.7 GHz) ATCA measurements of the early-type magnetic star ρ OphC, which is a flat-spectrum non-thermal radio source. The ρ OphC radio emission is partially circularly polarized with a steep spectral dependence: the fraction of polarized emission is about 60 at the lowest frequency sub-band (1.6 GHz) while is undetected at 16.7 GHz. This is clear evidence of coherent Auroral Radio Emission (ARE) from the ρ OphC magnetosphere. Interestingly, the detection of the ρ OphC's ARE is not related to a peculiar rotational phase. This is a consequence of the stellar geometry, which makes the strongly anisotropic radiation beam of the amplified radiation always pointed towards Earth. The circular polarization sign evidences mainly amplification of the ordinary mode of the electromagnetic wave, consistent with a maser amplification occurring within dense regions. This is indirect evidence of the plasma evaporation from the polar caps, a phenomenon responsible for the thermal X-ray aurorae. ρ OphC is not the first early-type magnetic star showing the O-mode dominated ARE but is the first star with the ARE always on view.

This thesis primarily uses submillimetre molecular line data from HARP, a heterodyne array on the James Clerk Maxwell Telescope (JCMT), to further investigate star formation in the Ophiuchus L1688 cloud. HARP was used to observe CO J = 3-2 isotopologues: 12CO, 13CO and C18O; and the dense gas tracer HCO+ J = 4-3. A method for calculating molecular line contamination in the SCUBA-2 450 and 850 μm dust continuum data was developed, which can be used to convert 12CO J =6-5and J =3-2 maps of integrated intensity (K km s−1) to molecular line flux (mJy beam−1) contaminating the continuum emission. Using HARP maps of 12CO J = 3-2, I quantified the amount of molecular line contamination found in the SCUBA-2 850 μm maps of three different regions, including NGC 1333 of Perseus and NGC 2071 and NGC 2024 of Orion B. Regions with ‘significant’ (i.e. > 20%) molecular line contamination correspond to molecular outflows. This method is now being used to remove molecular line contamination from regions with both SCUBA-2 dust continuum and HARP 12CO map coverage in the Gould Belt Legacy Survey (GBS). The Ophiuchus L1688 cloud was observed in all three CO J = 3-2 isotopologues. I carried out a molecular outflow analysis in the region on a list of 30 sources from the Spitzer ‘c2d’ survey [Evans et al., 2009]. Out of the 30 sources, 8 had confirmed bipolar outflows, 20 sources had ‘confused’ outflow detections and 2 sources did not have outflow detections. The Ophiuchus cloud was found to be gravitationally bound with the turbulent kinetic energy a factor of 7 lower than the gravitational binding energy. The high-velocity outflowing gas was found to be only 21% of the turbulence in the cloud, suggesting outflows are significant but not the dominant source of turbulence in the region. Other factors were found to influence the global high-velocity outflowing gas in addition to molecular outflows, including hot dust from nearby B-type stars, outflow remnants from less embedded sources and stellar winds from the Upper Scorpius OB association. To trace high density gas in the Ophiuchus L1688 cloud, HCO+ J = 4-3 was observed to further investigate the relationship between high column density and high density in the molecular cloud. Non-LTE codes RADEX and TORUS were used to develop density models corresponding to the HCO+ emission. The models involved both constant density and peaked density profiles. RADEX [van der Tak et al., 2007] models used a constant density model along the line-of-sight and indicated the HCO+ traced densities that were predominantly subthermally excited with den- sities ranging from 10^3–10^5 cm^−3. Line-of-sight estimates ranged from several parsecs to 90 pc, which was unrealistic for the Ophiuchus cloud. This lead to the implementation of peaked density profiles using the TORUS non-LTE radiative transfer code. Initial models used a ‘triangle’ density profile and a more complicated log-normal density probability density function (PDF) profile was subsequently implemented. Peaked density models were relatively successful at fitting the HCO+ data. Triangle models had density fits ranging from 0.2–2.0×10^6 cm^−3 and 0.1–0.3×10^6 cm^−3 for the 0.2 and 0.3 pc cloud length models re- spectively. Log-normal density models with constant-σ had peak density ranges from 0.2–1.0 ×10^5 cm^−3 and 0.6–2.0×10^5 cm^−3 for 0.2 and 0.3 pc models respectively. Similarly, log-normal models with varying-σ had lower and upper density limits corresponding to the range of FWHM velocities. Densities (lower and upper limits) ranged from 0.1–1.0 ×10^6 and 0.5–3.0 ×10^5 cm^-3 for the 0.2 and 0.3 pc models respectively. The result of the HCO+ density modelling indicated the distributions of starless, prestellar and protostellar cores do not have a preference for higher densities with respect to the rest of the cloud. This is contrary to past research suggesting the probability of finding a submillimetre core steeply rises as a function of column density (i.e. density; Belloche et al. 2011; Hatchell et al. 2005). Since the majority of sources are less embedded (i.e Class II/III), it is possible the evolutionary state of Ophiuchus is the main reason the small sample of Class 0/I protostars do not appear to have a preference for higher densities in the cloud.

The notion that low- to intermediate-mass young stellar objects (YSOs) gain mass at a constant rate during the early stages of their evolution appears to be challenged by observations of YSOs suffering sudden increases of the rate at which they gain mass from their circumstellar discs. Also, this idea that stars spend most of their lifetime with a low accretion rate and gain most of their final mass during short-lived episodes of high accretion bursts, helps to solve some long-standing problems in stellar evolution. The original classification of eruptive variables divides them in two separate subclasses known as FU Orionis stars (FUors) and EX Lupi stars (EXors). In this classical view FUors are at an early evolutionary stage and are still gaining mass from their parent envelopes, whilst EXors are thought to be older objects only surrounded by an accretion disc. The problem with this classical view is that it excludes younger protostars which have higher accretion rates but are too deeply embedded in circumstellar matter to be observed at optical wavelengths. Optically invisible protostars have been observed to display large variability in the near-infrared. These and some recent discoveries of new eruptive variables, show characteristics that can be attributed to both of the optically-defined subclasses of eruptive variables. The new objects have been proposed to be part of a new class of eruptive variables. However, a more accepted scenario is that in fact the original classes only represent two extremes of the same phenomena. In this sense eruptive variability could be explained as arising from one physical mechanism, i.e. unsteady accretion, where a variation in the parameters of such mechanism can cause the different characteristics observed in the members of this class. With the aim of studying the incidence of episodic accretion among young stellar objects, and to characterize the nature of these eruptive variables we searched for high amplitude variability in two multi-epoch infrared surveys: the UKIDSS Galactic Plane Survey (GPS) and the Vista Variables in the Via Lactea (VVV). In order to further investigate the nature of the selected variable stars, we use photometric information arising from public surveys at near- to farinfrared wavelengths. In addition we have performed spectroscopic and photometric follow-up for a large subset of the samples arising from GPS and VVV. We analyse the widely separated two-epoch K-band photometry in the 5th, 7th and 8th data releases of the UKIDSS Galactic Plane Survey. We find 71 stars with _K > 1 mag, including 2 previously known OH/IR stars and a Nova. Even though the mid-plane is mostly excluded from the dataset, we find the majority (66%) of our sample to be within known star forming regions (SFRs), with two large concentrations in the Serpens OB2 association (11 stars) and the Cygnus-X complex (27 stars). The analysis of the multi-epoch K-band photometry of 2010-2012 data from VVV covering the Galactic disc at |b| < 1◦ yields 816 high amplitude variables, which include known variables of different classes such as high mass X-ray binaries, Novae and eclipsing binaries among others. Remarkably, 65% of the sample are found concentrated towards areas of star formation, similar to the results from GPS. In both surveys, sources in SFRs show spectral energy distributions (SEDs) that support classification as YSOs. This indicates that YSOs dominate the Galactic population of high amplitude infrared variable stars at low luminosities and therefore likely dominate the total high amplitude population. Spectroscopic follow-up allows us to confirm the pre-main sequence nature of several GPS and VVV Objects. Most objects in both samples show spectroscopic signatures that can be attributed to YSOs undergoing high states of accretion, such as veiling of photospheric features and CO emission, or show FUor-like spectra. We also find a large fraction of objects with 2.12 μm H2 emission that can be explained as arising from shock-excited emission caused by molecular outflows. Whether these molecular outflows are related to outbursts events cannot be confirmed from our data. Adding the GPS and VVV spectroscopic results, we find that between 6 and 14 objects are new additions to the FUor class from their close resemblance to the near-infrared spectra of FUors, and at least 23 more objects are new additions to the eruptive variable class. For most of these we are unable to classify them into any of the original definitions for this variable class. In any case, we are adding up to 37 new stars to the eruptive variable class which would double the current number of known objects. We note that most objects are found to be deeply embedded optically invisible stars, thus increasing the number of objects belonging to this subclass by a much larger factor. In general, objects in our samples which are found to be likely eruptive variable stars show a mixture of characteristics that can be attributed to both of the optically-defined classes. This agrees well with the recent discoveries in the literature. Finally, we are able to derive a first rough estimate on the incidence of episodic accretion among class I YSOs in the star-forming complex G305. We find that _ 9% of such objects are in a state of high accretion. This number is in agreement with previous theoretical and observational estimates among class I YSOs.

The blue hypergiant Cyg OB2 12 (B3Ia+) is a representative member of the class of very massive stars in a poorly understood evolutionary stage. We obtained its high-resolution X-ray spectrum using the Chandra observatory. PoWR model atmospheres were calculated to provide realistic wind opacities and to establish the wind density structure. We find that collisional de-excitation is the dominant mechanism depopulating the metastable upper levels of the forbidden lines of the He-like ions Si xivand Mg xii. Comparison between the model and observations reveals that X-ray emission is produced in a dense plasma, which could reside only at the photosphere or in a colliding wind zone between binary components. The observed X-ray spectra are well-fitted by thermal plasma models, with average temperatures in excess of 10 MK. The wind speed in Cyg OB2 12 is not high enough to power such high temperatures, but the collision of two winds in a binary system can be sufficient. We used archival data to investigate the X-ray properties of other blue hypergiants. In general, stars of this class are not detected as X-ray sources. We suggest that our new Chandra observations of Cyg OB2 12 can be best explained if Cyg OB2 12 is a colliding wind binary possessing a late O-type companion. This makes Cyg OB2 12 only the second binary system among the 16 known Galactic hypergiants. This low binary fraction indicates that the blue hypergiants are likely products of massive binary evolution during which they either accreted a significant amount of mass or already merged with their companions.