Since the end of 2018, the Transiting Exoplanet Survey Satellite (TESS) provides high-quality space data on stellar photometry to the astronomical community. We present the results of an analysis of TESS photometric data for known slowly rotating magnetic chemically peculiar (mCP) stars. In general, mCP stars show an inhomogeneous distribution of elements in their stellar atmospheres that leads to spectroscopic (line profile) and photometric (light curve) variations over the rotation period. In the frame of the oblique magnetic rotator (OMR) model, patches of enhanced chemical abundance on the stellar surface reveal the frequency of stellar rotation. Using this approach, we have compiled a list of slowly rotating mCP stars with rotation periods longer than two days from the analysis of the photometric data provided by TESS for the first eight sectors of observations. Slowly rotating mCP stars usually possess a hydrodynamically stable stellar atmosphere where a magnetic field can amplify the process of atomic diffusion and leads to the horizontal and vertical stratification of chemical abundances.

The new photometric data on pulsating Ap star HD~27463 obtained recently with the Transiting Exoplanet Survey Satellite (\textit{TESS}) are analysed to search for variability. Our analysis shows that HD~27463 exhibits two types of photometric variability. The low frequency variability with the period $P$ =~2.834274 $\pm$ 0.000008 d can be explained in terms of axial stellar rotation assuming the oblique magnetic rotator model and presence of surface abundance/brightness spots, while the detected high-frequency variations are characteristics of $\delta$~Scuti pulsations. From the analysis of Balmer line profiles visible in two FEROS spectra of HD~27463 we have derived its effective temperature and surface gravity, finding values that are close to those published for this star in the \textit{TESS} Input Catalogue (TIC). Knowing the rotation period and the v$\sin{i}$ value estimated from the fitting of Balmer line profiles we found that the rotational axis is inclined to the line of sight with an angle of $i=33\pm8\deg$. Our best-fitting model of the observed pulsation modes results in an overshoot parameter value $f_{ov} = 0.014$ and values of global stellar parameters that are in good agreement with the data reported in the TIC and with the data derived from fitting Balmer line profiles. This model indicates an age of 5.0 $\pm$~0.4 $\times 10^8$~yrs, which corresponds to a core hydrogen fraction of 0.33.

We present the first asteroseismic results for $\delta$ Scuti and $\gamma$ Doradus stars observed in Sectors 1 and 2 of the TESS mission. We utilise the 2-min cadence TESS data for a sample of 117 stars to classify their behaviour regarding variability and place them in the Hertzsprung-Russell diagram using Gaia DR2 data. Included within our sample are the eponymous members of two pulsator classes, $\gamma$ Doradus and SX Phoenicis. Our sample of pulsating intermediate-mass stars observed by TESS also allows us to confront theoretical models of pulsation driving in the classical instability strip for the first time and show that mixing processes in the outer envelope play an important role. We derive an empirical estimate of 74% for the relative amplitude suppression factor as a result of the redder TESS passband compared to the Kepler mission using a pulsating eclipsing binary system. Furthermore, our sample contains many high-frequency pulsators, allowing us to probe the frequency variability of hot young $\delta$ Scuti stars, which were lacking in the Kepler mission data set, and identify promising targets for future asteroseismic modelling. The TESS data also allow us to refine the stellar parameters of SX Phoenicis, which is believed to be a blue straggler.

In this paper we introduce the MOBSTER collaboration and lay out its scientific goals. We present first results based on the analysis of nineteen previously known magnetic O, B and A stars observed in 2-minute cadence in sectors 1 and 2 of the Transiting Exoplanet Survey Satellite (TESS) mission. We derive precise rotational periods from the newly obtained light curves and compare them to previously published values. We also discuss the overall photometric phenomenology of the known magnetic massive and intermediate-mass stars and propose an observational strategy to augment this population by taking advantage of the high-quality observations produced by TESS.

Uncertainties in stellar structure and evolution theory are largest for stars undergoing core convection on the main sequence. A powerful way to calibrate the free parameters used in the theory of stellar interiors is asteroseismology, which provides direct measurements of angular momentum and element transport. We report the detection and classification of new variable O and B stars using high-precision short-cadence (2-min) photometric observations assembled by the Transiting Exoplanet Survey Satellite (TESS). In our sample of 154 O and B stars, we detect a high percentage (90%) of variability. Among these we find 23 multiperiodic pulsators, 6 eclipsing binaries, 21 rotational variables, and 25 stars with stochastic low-frequency variability. Several additional variables overlap between these categories. Our study of O and B stars not only demonstrates the high data quality achieved by TESS for optimal studies of the variability of the most massive stars in the Universe, but also represents the first step towards the selection and composition of a large sample of O and B pulsators with high potential for joint asteroseismic and spectroscopic modeling of their interior structure with unprecedented precision.

We investigate the pulsation properties of stellar models representative of $\delta$ Scuti and $\gamma$ Doradus variables. We have calculated a grid of stellar models from 1.2 to 2.2 M$_{\odot}$, including the effects of both rotation and convective overshoot using MESA, and we investigate the pulsation properties of these models using GYRE. We discuss observable patterns in the frequency spacing for $p$ modes and the period spacings for g modes. Using the observable patterns in g mode period spacings, it may be possible to observationally constrain the convective overshoot and rotation of a model. We also calculate the pulsation constant (Q) for all models in our grid, and investigate the variation with convective overshoot and rotation. The variation in Q values of radial modes can be used to place constraints on the convective overshoot and rotation of stars in this region. As a test case, we apply this method to a sample of 22 high amplitude $\delta$ Scuti stars (HADS), and provide estimates for the convective overshoot of the sample.

We present multidimensional modeling of convection and oscillations in main-sequence stars somewhat more massive than the Sun, using three separate approaches: 1) Using the 3-D planar StellarBox radiation hydrodynamics code to model the envelope convection zone and part of the radiative zone. Our goals are to examine the interaction of stellar pulsations with turbulent convection in the envelope, excitation of acoustic modes, and the role of convective overshooting; 2) Applying the spherical 3-D MHD ASH (Anelastic Spherical Harmonics) code to simulate the core convection and radiative zone. Our goal is to determine whether core convection can excite low-frequency gravity modes, and thereby explain the presence of low frequencies for some hybrid gamma Doradus/delta Scuti variables for which the envelope convection zone is too shallow for the convective blocking mechanism to drive gravity modes; 3) Applying the ROTORC 2-D stellar evolution and dynamics code to calculate evolution with a variety of initial rotation rates and extents of core convective overshooting. The nonradial adiabatic pulsation frequencies of these nonspherical models are calculated using the 2-D pulsation code NRO. We present new insights into pulsations for stars of one to two solar masses gained by multidimensional modeling.

Among the most spectacular variable stars are the Luminous Blue Variables (LBVs), which can show three types of variability. The LBV phase of evolution is poorly understood, and the driving mechanisms for the variability are not known. The most common type of variability, the S Dor instability, occurs on timescales of tens of years. During an S Dor outburst, the visual magnitude of the star increases, while the bolometric magnitude stays approximately constant. In this work, we investigate pulsation as a possible trigger for the S Dor type outbursts. We calculate the pulsations of envelope models using a nonlinear hydrodynamics code including a time-dependent convection treatment. We initialize the pulsation in the hydrodynamic model based on linear non-adiabatic calculations. Pulsation properties for a full grid of models from 20 to 85 M$_{\odot}$ were calculated, and in this paper we focus on the few models that show either long-period pulsations or outburst-like behaviour, with photospheric radial velocities reaching 70-80 km/s. At the present time, our models cannot follow mass loss, so once the outburst event begins, our simulations are terminated. Our results show that pulsations alone are not able to drive enough surface expansion to eject the outer layers. However, the outbursts and long-period pulsations discussed here produce large variations in effective temperature and luminosity, which are expected to produce large variations in the radiatively driven mass-loss rates.

We have calculated the pulsations of massive stars using a nonlinear hydrodynamic code including time-dependent convection. The basic structure models are based on a standard grid published by Meynet et al. (1994). Using the basic structure, we calculated envelope models, which include the outer few percent of the star. These models go down to depths of at least 2 million K. These models, which range from 40 to 85 solar masses, show a range of pulsation behaviours. We find models with very long period pulsations ( $>$ 100 d), resulting in high amplitude changes in the surface properties. We also find a few models that show outburst-like behaviour. The details of this behaviour are discussed, including calculations of the resulting wind mass-loss rates.

The Luminous Blue Variable stars exhibit behavior ranging from light curve `microvariations' on timescales of tens of days, to `outbursts' accompanied by mass loss of up to 10e-03 solar masses per year, occurring decades apart, to `giant eruptions' such as seen in Eta Carinae ejecting one or more solar masses and recurring on timescales of centuries. Here we review the work of the Los Alamos group since 1993 to investigate pulsations and instabilities in massive stars using linear pulsation models and non-linear hydrodynamic models. The models predict pulsational variability that may be associated with the microvariations. Using a nonlinear pulsation hydrodynamics code with a time-dependent convection treatment, we show that, in some circumstances, the Eddington limit is exceeded periodically in the pulsation driving region of the stellar envelope, accelerating the outer layers, and perhaps initiating mass loss or LBV outbursts. We discuss how pulsations and mass loss may be responsible for the location of the Humphreys-Davidson Limit in the H-R diagram. The `giant eruptions', however, must involve much deeper regions in the stellar core to cause such large amounts of mass to be ejected. We review and suggest some possible explanations, including mixing from gravity modes, secular instabilities, the epsilon mechanism, or the SASI instability as proposed for Type II supernovae. We outline future work and required stellar modeling capabilities to investigate these possibilities.

The first stars are the key to the formation of primitive galaxies, early cosmological reionization and chemical enrichment, and the origin of supermassive black holes. Unfortunately, in spite of their extreme luminosities, individual Population III stars will likely remain beyond the reach of direct observation for decades to come. However, their properties could be revealed by their supernova explosions, which may soon be detected by a new generation of NIR observatories such as JWST and WFIRST. We present light curves and spectra for Pop III pair-instability supernovae calculated with the Los Alamos radiation hydrodynamics code RAGE. Our numerical simulations account for the interaction of the blast with realistic circumstellar envelopes, the opacity of the envelope, and Lyman absorption by the neutral IGM at high redshift, all of which are crucial to computing the NIR signatures of the first cosmic explosions. We find that JWST will detect pair-instability supernovae out to z > 30, WFIRST will detect them in all-sky surveys out to z ~ 15 - 20 and LSST and Pan-STARRS will find them at z ~ 7 - 8. The discovery of these ancient explosions will probe the first stellar populations and reveal the existence of primitive galaxies that might not otherwise have been detected.

The detection of Pop III supernovae could directly probe the primordial IMF for the first time, unveiling the properties of the first galaxies, early chemical enrichment and reionization, and the seeds of supermassive black holes. Growing evidence that some Pop III stars were less massive than 100 solar masses may complicate prospects for their detection, because even though they would have been more plentiful they would have died as core-collapse supernovae, with far less luminosity than pair-instability explosions. This picture greatly improves if the SN shock collides with a dense circumstellar shell ejected during a prior violent LBV type eruption. Such collisions can turn even dim SNe into extremely bright ones whose luminosities can rival those of pair-instability SNe. We present simulations of Pop III Type IIn SN light curves and spectra performed with the Los Alamos RAGE and SPECTRUM codes. Taking into account Lyman-alpha absorption in the early universe and cosmological redshifting, we find that 40 solar mass Pop III Type IIn SNe will be visible out to z ~ 20 with JWST and out to z ~ 7 with WFIRST. Thus, even low mass Pop III SNe can be used to probe the primeval universe.

A large number of massive stars are known to rotate, resulting in significant distortion and variation in surface temperature from the pole to the equator. Radiatively driven mass loss is temperature dependent, so rapid rotation produces variation in mass loss and angular momentum loss rates across the surface of the star, which is expected to affect the evolution of rapidly rotating massive stars. In this work, we investigate the two dimensional effects of rotation on radiatively driven mass loss and the associated loss of angular momentum in ZAMS models with solar metallicity. Using 2D stellar models, which give the variation in surface parameters as a function of co-latitude, we implement two different mass loss prescriptions describing radiatively driven mass loss. We find a significant variation in mass loss rates and angular momentum loss as a function of co-latitude. We find that the mass loss rate decreases as the rotation rate increases for models at constant initial mass, and derive scaling relations based on these models. When comparing 2D to 1D mass loss rates, we find that although the total angle integrated mass loss does not differ significantly, the 2D models predict less mass loss from the equator and more mass loss from the pole than the 1D predictions using von Zeipel's law. As a result, rotating models lose less angular momentum in 2D than in 1D, which will change the subsequent evolution of the star. The evolution of these models will be investigated in future work.

(abridged) Recent work on several beta Cephei stars has succeeded in constraining both their interior rotation profile and their convective core overshoot. In particular, a recent study focusing on theta$ Oph has shown that a convective core overshoot parameter of alpha = 0.44 is required to model the observed pulsation frequencies, significantly higher than for other stars of this type. We investigate the effects of rotation and overshoot in early type main sequence pulsators, and attempt to use the low order pulsation frequencies to constrain these parameters. This will be applied to a few test models and theta Oph. We use a 2D stellar evolution code and a 2D linear adiabatic pulsation code to calculate pulsation frequencies for 9.5 Msun models. We calculate low order p-modes for models with a range of rotation rates and convective core overshoot parameters. Using these models, we find that the convective core overshoot has a larger effect on the pulsation frequencies than the rotation, except in the most rapidly rotating models considered. When the differences in radii are accounted for by scaling the frequencies, the effects of rotation diminish, but are not entirely accounted for. We find that increasing the convective core overshoot decreases the large separation, while producing a slight increase in the small separations. We created a model frequency grid which spanned several rotation rates and convective core overshoot values. Using a modified chi^2 statistic, we are able to recover the rotation velocity and core overshoot for a few test models. Finally, we discuss the case of the beta Cephei star theta Oph. Using the observed frequencies and a fixed mass and metallicity, we find a lower overshoot than previously determined, with alpha = 0.28 +/- 0.05. Our determination of the rotation rate agrees well with both previous work and observations, around 30 km/s.

We have investigated the effects of uniform rotation and a specific model for differential rotation on the pulsation frequencies of 10 \Msun\ stellar models. Uniform rotation decreases the frequencies for all modes. Differential rotation does not appear to have a significant effect on the frequencies, except for the most extreme differentially rotating models. In all cases, the large and small separations show the effects of rotation at lower velocities than do the individual frequencies. Unfortunately, to a certain extent, differential rotation mimics the effects o f more rapid rotation, and only the presence of some specific observed frequencies with well identified modes will be able to uniquely constrain the internal rotation of pulsating stars.

Spectral energy distributions for models of arbitrarily rotating stars are computed using two dimensional rotating stellar models, NLTE plane parallel model atmospheres, and a code to integrate the appropriately weighted intensities over the visible surface of the stellar disk. The spectral energy distributions depend on the inclination angle between the observer and the rotation axis of the model. We use these curves to deduce what one would infer the model's luminosity and effective temperature to be assuming the object was nonrotating.

Radial and nonradial oscillations offer the opportunity to investigate the interior properties of stars. We use 2D stellar models and a 2D finite difference integration of the linearized pulsation equations to calculate non-radial oscillations. This approach allows us to directly calculate the pulsation modes for a distorted rotating star without treating the rotation as a perturbation. We are also able to express the finite difference solution in the horizontal direction as a sum of multiple spherical harmonics for any given mode. Using these methods, we have investigated the effects of increasing rotation and the number of spherical harmonics on the calculated eigenfrequencies and eigenfunctions and compared the results to perturbation theory. In slowly rotating stars, current methods work well, and we show that the eigenfunction can be accurately modelled using 2nd order perturbation theory and a single spherical harmonic. We use 10 Msun models with velocities ranging from 0 to 420 km/s (0.89 Omega_c) and examine low order p modes. We find that one spherical harmonic remains reasonable up to a rotation rate around 300km s^{-1} (0.69 Omega_c) for the radial fundamental mode, but can fail at rotation rates as low as 90 km/s (0.23 Omega_c) for the 2H mode or l = 2 p_2 mode, based on the eigenfrequencies alone. Depending on the mode in question, a single spherical harmonic may fail at lower rotation rates if the shape of the eigenfunction is taken into consideration. Perturbation theory, in contrast, remains valid up to relatively high rotation rates for most modes. We find the lowest failure surface equatorial velocity is 120 km/s (0.30 Omega_c) for the l = 2 p_2 mode, but failure velocities between 240 and 300 km/s (0.58-0.69 Omega_c)are more typical.

The spectral energy distribution (SED) of a non-spherical star could differ significantly from the SED of a spherical star with the same average temperature and luminosity. Calculation of the SED of a deformed star is often approximated as a composite of several spectra, each produced by a plane parallel model of given effective temperature and gravity. The weighting of these spectra over the stellar surface, and hence the inferred effective temperature and luminosity, will be dependent on the inclination of the rotation axis of the star with respect to the observer, as well as the temperature and gravity distribution on the stellar surface. Here we calculate the surface conditions of rapidly rotating stars with a 2D stellar structure and evolution code and compare the effective temperature distribution to that predicted by von Zeipel's law. We calculate the composite spectrum for a deformed star by interpolating within a grid of intensity spectra of plane parallel model atmospheres and integrating over the surface of the star. Using this method, we find that the deduced variation of effective temperature with inclination can be as much as 3000 K for an early B star, depending on the details of the underlying model.