EDP Sciences
Free Access
Volume 511, February 2010
Article Number L7
Number of page(s) 4
Section Letters
DOI https://doi.org/10.1051/0004-6361/200913832
Published online 10 March 2010
A&A 511, L7 (2010)


HR 7355 - another rapidly braking He-strong CP star?[*]

Z. Mikulásek1,5 - J. Krticka1 - G. W. Henry2 - S. N. de Villiers3 - E. Paunzen4 - M. Zejda1

1 - Department of Theoretical Physics and Astrophysics, Masaryk University, Brno, Czech Republic
2 - Center of Excellence in Information Systems, Tennessee State University, Nashville, Tennessee, USA
3 - Private Observatory, 61 Dick Burton Road, Plumstead, Cape Town, South Africa
4 - Institute for Astronomy of the University of Vienna, Vienna, Austria
5 - Observatory and Planetarium of J. Palisa, VSB - Technical University, Ostrava, Czech Republic

Received 9 December 2009 / Accepted 4 February 2010

Context. Strong meridional mixing induced by rapid rotation is one reason why all hot main-sequence stars are not chemically peculiar. However, the finding that the He-strong CP star HR 7355 is a rapid rotator complicates this concept.
Aims. Our goal is to explain the observed behaviour of HR 7355 based on period analysis of all available photometry.
Methods. Over two years, we acquired 114 new BV observations of HR 7355 at observatories in Arizona, U.S.A and Cape Town, South Africa. We performed period analyses of the new observations along with new analyses of 732 archival measurements from the Hipparcos and ASAS projects.
Results. We find that the light curves of HR 7355 in various filters are quite similar, with amplitudes 0.035(4), 0.036(4), and 0.038(3) mag in B, Hp, and V, respectively. The light curves are double-peaked, with unevenly deep minima. We substantially refine the rotational period to be $P=0\hbox{$.\!\!^{\rm d}$ }5214410(4)$, indicating that HR 7355 is the most rapidly rotating CP star known. Our period analyses reveal a possible lengthening of the rotational period with $\dot{P}/P=2.4(8)\times 10^{-6}~\rm {yr}^{-1}$.
Conclusions. We conclude that the shape and amplitude of HR 7355 light curves are typical of magnetic He-strong CP stars, for which light variations are the result of photospheric spots on the surface of a rotating star. We hypothesise that the light variations are caused mainly by an uneven distribution of overabundant helium on the star's surface. We briefly describe and discuss the cause of the rapid rotational braking of the star.

Key words: stars: chemically peculiar - stars: variables: general - stars: individual: HR 7355 - stars: rotation

1 Introduction

Chemically peculiar (CP) stars are an important class of stars that occupy the upper main sequence, where radiative diffusion and gravitational settling result in atmospheric chemical abundances that differ remarkably from the Sun's. According to their patterns of chemical anomalies, CP stars are classified into several subclasses that also follow a temperature sequence. Among them, SrCrEu, Si, He-weak, and He-strong stars have strong global magnetic fields. These ``magnetic'' CP (mCP) stars also exhibit synchronous variability in their spectra and brightness with periods longer than one half of a day. Their photometric amplitudes are a few hundredths of magnitude.

mCP stars have inhomogeneous surface distributions of chemical elements as determined from their rotationally modulated spectral-line variability (e.g., Lüftinger et al. 2010b,a). The uneven surface distribution of various elements, together with rotation, has been presumed to be the main cause of these stars' light variability. Line blanketing caused by overabundant elements, mainly iron-group metals and b-f transitions, may induce the flux redistribution providing the mechanism for the light variability (Lanz et al. 1996). A strong magnetic field may also influence the light curves (e.g., Kochukhov et al. 2005).

Krticka et al. (2007,2009) used surface abundance maps to simulate successfully the light curves of the He-strong star HD 37776 and the Si star HR 7224. They demonstrated that the inhomogeneous surface distribution of silicon, iron, and helium, along with the b-f and b-b transitions of these elements, accounted for most of the light variability in these CP stars.

The light curves of magnetic CP stars are stable on a timescale of decades or more and so indicate persistency of their spectroscopic and photometric spots. We can observe period changes in their light curves in only a few cases (Mikulásek et al. 2008b, and references therein).

It is generally expected that the meridional currents induced by the rapid rotation of hot MS stars are able to erase the effects of the slow diffusion of chemical elements and thus prevent the formation of CP abundance anomalies. We study therefore the light variability of HR 7355, which is one of the most rapidly rotating mCP stars.

2 The star

HR 7355 (HD 182180, HIP 95480) is a bright (V=6.02 mag) but poorly studied southern B2V star. It is known to be a rapid rotator, Oksala et al. (2010) evaluating its projected equatorial velocity to be $v\sin i=(300\pm15)~{\rm km~s}^{-1}$. The Hipparcos catalog (ESA 1997) classified the star as an ``unsolved variable'', but Koen & Eyer (2002) reanalysed the Hipparcos photometry and found a frequency of variation of $1/0\hbox{$.\!\!^{\rm d}$ }26072$.

The first detailed study of HR 7355 was performed by Rivinius et al. (2008), who describe the star as a helium-strong CP star with Balmer emission. They argue that the true period, corresponding to the star's rotational period, is twice that derived by Koen & Eyer (2002): $P=0\hbox{$.\!\!^{\rm d}$ }521428(6)$. This suggests that HR 7355 is the most rapidly rotating He-strong CP star known.

In the present study, we analyse our new and archival photometry of HR 7355 to determine the star's periodicity and light curve shape at various epochs between 1990 and 2009.

2.1 Parameters of the star

We estimate the mass and age of HR 7355 using the revised version of the Hipparcos catalogue (van Leeuwen 2007), the Geneva and Strömgren photometry (from GCPD[*]), and the isochrones of the Padova group (Marigo et al. 2008).

The Geneva colours (Cramer 1999) are consistent with a B2 V star. The Z index shows no peculiarity beyond what we would expect for a classical CP4 star (Paunzen et al. 2005). From the various Geneva indices and the available Strömgren photometry, we deduce a reddening of E(b-y) = 0.05 mag and an effective temperature of 18 000 K. The photometric estimates for the absolute magnitude, MV, vary between -1.75 and -1.95 mag. The parallax gives a distance of 273(26) pc and an absolute magnitude of  MV=-1.38(0.20) mag (assuming V = 6.02 mag). This result is fainter than those estimated from the photometric indices, which may be a result of the star's rapid rotation.

We estimate the age of HR 7355 to be between 15 and 25 Myr, as derived from the parallax-based and the photometry-based values of MV, respectively. We note that Westin (1985) lists $\log
t=7.379$ (23.9 Myr) for HR 7355. Consequently, HR 7355 appears to be near the middle of its main-sequence lifetime. We derive a mass of 6.3(0.3) $M_{\odot}$ from the parallax-based MV. The photometry-based absolute magnitude would result in a higher mass.

From its distance and the Galactic coordinates, we derive the Galactic [X,Y] distances/coordinates of HR 7355 to be [+254,+48] pc. Comparing these values and their uncertainties with the borders of the Scorpius-Centaurus association (Preibisch & Mamajek 2008), we find that HR 7355 is probably beyond the outer edge of the association. The mean proper motion of the association is [-25,-10] mas yr-1 (de Zeeuw et al. 1999), compared to [+12,-15] mas yr-1 for HR 7355. Given its motion direction and apparent age, we suggest that HR 7355 does not originate in the Scorpius-Centaurus association.

3 Photometric data of HR 7355

3.1 New photometry

Most of our new photometry of HR 7355 was acquired with the T3 0.4 m automatic photoelectric telescope (APT) at Fairborn Observatory in southern Arizona during 2008 May and June. The T3 APT uses a temperature-stabilised EMI 9924B photomultiplier tube to measure photon count rates sequentially through B and V filters. For additional details about the collection and reduction of the APT data, see cf. Mikulásek et al. (2008b) and references therein. The comparison and check star used by the APT were HD 179520 and HD 181240, respectively. The southern declination of HR 7355 ( $\delta=-28^{\circ}$) causes airmass values of $\sim$2 for all of our APT measurements. The standard deviation of the difference between variable and comparison stars is $\sim$0.01 mag, twice the typical scatter with this telescope. Thus we acquired 43 B and 38 V measurements.

We also obtained seven continuous hours of BV photometry with a photoelectric photometer attached to a 0.28 m reflector on the night of 2009 July 26 UT at the private observatory of one of us (SNdV) in Cape Town, South Africa. We acquired 17 and 16 measurements in Vand B. The comparison and check stars were the same as used by the APT; the scatter in the observations was 0.01 mag. When combined with the APT measurements, these continuous single-night observations helped us to minimise the aliasing inherent in single-sight data sets.

The newly acquired photometric data are available at the CDS[*] or the On-line database of photometric observations of mCP stars[*].

3.2 Archival photometry

The Hipparcos photometry (ESA 1997) of HR 7355 consists of 57 measurements in $B_{{\rm T}}$ (scatter $\sigma $ of 23 mmag), 45 in Hp ($\sigma=9$ mmag), and 56 in $V_{{\rm T}}$($\sigma=28$ mmag). The data cover the period from 1990 March to 1993 March.

The continuously updated archive of the All Sky Automated Survey (ASAS) (Pojmanski 2002) is a useful source of photometric observations of stars near the equator. Unfortunately, the precision for relatively bright stars is poor - the data have significant non-Gaussian scatter with plenty of outliers. Furthermore, the data contain a long-term trend with a range of 0.05 mag.

Nonetheless, the ASAS data are the only photometric data of HR 7355 available in the period between the end of the Hipparcos mission (1993) and the start of our own photometry (2008). We extracted only the individual V measurements with a quality of ``A'' or ``B'' and removed the long-term trend. This resulted in an ASAS data set of 574 V observations covering the time interval from 2001 February to 2009 October. These data have a scatter of 50 mmag.

Therefore, we consider in our period analysis a total of 846 measurements spread over 20 years.

4 The ephemeris

4.1 The period

The periodogram of all available HR 7355 photometry displays only two prominent peaks - one at a period of $0\hbox{$.\!\!^{\rm d}$ }52144$ and the second at $0\hbox{$.\!\!^{\rm d}$ }26072$ - that agree well with the previous determinations of Koen & Eyer (2002) and Rivinius et al. (2008). The shorter period is the consequence of the double-peaked nature of the light curve (see Fig. 1). This $0\hbox{$.\!\!^{\rm d}$ }26072$ period is shorter than the critical period of centrifugal disruption of main-sequence stars of the same spectral type and so cannot represent HR 7355's rotation period, which we take to be $0\hbox{$.\!\!^{\rm d}$ }52144$. The longer period is also compatible with the period of the well documented spectral and spectropolarimetric variations (Rivinius et al. 2010; Oksala et al. 2010).

\par\includegraphics[width=6.7cm,clip]{13832fg1.eps} \end{figure} Figure 1:

BHpV light curves plotted versus the quadratic phase. $\Box $ - Hipparcos, $\circ $ - Cape Town, $\triangle $ - APT, $\blacklozenge $ - ASAS (each point represents the mean of 30 adjacent ASAS measurements).

Open with DEXTER

4.2 The light variations model

The double-waved light curves in B, Hp, and V are similar (see Fig. 1), differing only slightly in their effective amplitudes Ac (for the definition see Mikulásek et al. 2007), where the subscript c denotes the filter colour. Light-curve magnitudes can then be expressed as

\end{displaymath} (1)

where mcj(t) is the magnitude in colour c observed by j-th observer, $\overline{m}_{cj}$ is the mean magnitude, which can be variable over the long term. This was the case for ASAS measurements, where we had to assume a cubic trend. The function $F(\vartheta)$ is the simplest normalised periodic function that represents the observed photometric variations of HR 7355 in detail. The phase of maximum brightness is defined to be 0.0, and the effective amplitude is defined to be 1.0. The function, being the sum of three terms, is described by two parameters $\beta_1$ and $\beta_2$, where $\beta_1$ quantifies the difference in depths of the primary and the secondary minima and $\beta_2$expresses any asymmetry in the light curve
                  $\displaystyle F(\vartheta,~ \beta_1,~\beta_2)$ = $\displaystyle \beta_1~\cos(2~\pi~\vartheta)+
\sqrt{1\!-\!\beta_1^2\!-\!\beta_2^2}\ \cos(4~\pi~\vartheta)$  
    $\displaystyle +~\beta_2\left[\textstyle{\frac{2}{\sqrt{5}}}~\sin(2~\pi~\vartheta)
-\textstyle{\frac{1}{\sqrt{5}}}\sin(4~\pi~ \vartheta)\right],$ (2)

where $\vartheta$ is the phase function of time described in Mikulásek et al. (2008b). The O-C diagram (see Fig. 2) indicates that the function may be parabolic; thus, we assumed it to have the form

\begin{displaymath}\vartheta \cong
...1-\alpha_2\right);\quad \vartheta_1=\left(t-M_{01}\right)/P_1,
\end{displaymath} (3)

where M01 is the maximum nearest the weighted center of the observations, P1 is the period of the linear fit, $\dot{P}$ is the time derivative of the period (assumed to be constant here), and $\alpha_1$ and $\alpha_2$ are constants determined by the data distribution in time so that the first and second terms in Eq. (3) are orthogonal, such that $\alpha_{1,2}=\Theta_1 \pm
\sqrt{\Theta_1^2+\overline{\vartheta_1^2}}$, where $\Theta_1=

All 19 model parameters were computed simultaneously by a weighted non-linear LSM regression applied to the complete observational material without any artificial divisions.

\end{figure} Figure 2:

Time variation in mutual phase shifts of observed light curves in days computed with a linear ephemeris. $\Box $ - Hipparcos, $\diamondsuit $ - ASAS, $\triangle $ - Arizona, and $\circ $ - Cape Town. The Hipparcos and ASAS data sets are divided into two pieces with equal numbers of observations and plotted separately. This division was made for only illustrative reasons and did not affect the period analysis in any way. The quadratic time function is denoted by the solid line, while the 1 $\sigma $-uncertainties are plotted with dashed lines.

Open with DEXTER

4.3 Results

We determined the model parameters $\alpha_1=4420,\ \alpha_2=-7416$, $M_{01}=2~452~367.1221(24),\ P_1=0\hbox{$.\!\!^{\rm d}$ }5~214~410(4),\
\dot{P}=3.4(1.1)\times 10^{-9}=0.107(35)$  s y-1, and $\dot{P}/P_1=2.4(8)\times 10^{-6}~{\rm yr}^{-1}$. Time of the primary maximum being given by

                     $\displaystyle {\rm HJD}_{{\rm maxI}}(E)$ $\textstyle \cong$ $\displaystyle 2~452~367.1221(24)+0\hbox{$.\!\!^{\rm d}$ }
5~214~410(4)\times E$  
    $\displaystyle +8\hbox{$.\!\!^{\rm d}$ }8(2\hbox{$.\!\!^{\rm d}$ }9)\times10^{-10}\ (E+7416)~(E-4420),$ (4)

where E is an integer. The effective amplitudes Ac are almost identical and equal to 0.035(4), 0.036(4), and 0.038(3) mag in B, Hp, and V, respectively, $\beta_1=0.08(6)$, and $\beta_2=-0.26(6)$.

The reliability of the LSM procedure results depends critically on the appropriateness of the chosen model for the description of the reality - in our case light variations over the past twenty years. Inspecting Fig. 1, we can conclude that the model of light curves described by Eqs. (1) and (2) is fully adequate.

Adequacy of the quadratic model used for the phase function Eq. (3) can be tested by changes in mutual phase shifts of observed light curves. The shift (O-C)k can be evaluated for any subset k of observational data by the relation: (O-C) $_k=-\overline{\Delta
m_{k}~\dot{m}_{k}}/\overline{\dot{m}^2_{k}}$, where $\Delta m_{k}$are deviations of observed magnitudes from the fit, and $\dot{m}_{k}$ are time derivatives of the fitted function (Mikulásek et al. 2006). From Fig. 2 containing shifts of 6 observational data subsets with respect to the linear ephemeris, we conclude that our quadratic model of the phase function in Eq. (3) is acceptable.

To be sure that non-zero $\dot{P}$ is not a mere conjunction of the non-Gaussian noise in ASAS data, we treated the observational set excluding all the ASAS data. Even in this case, we obtained $\dot
P=2.8(1.6)\times 10^{-9}$, supporting the reality of the period lengthening.

5 Discussion

5.1 HR 7355 - the most rapidly rotating CP star

The $P_1=0\hbox{$.\!\!^{\rm d}$ }52144$ rotational period of HR 7355 is among the shortest known to date for CP stars: HD 164429 ( $P=0\hbox{$.\!\!^{\rm d}$ }51899$, Adelman 1999); HD 124224 = CU Vir ( $P=0\hbox{$.\!\!^{\rm d}$ }52070$, Sokolov 2000); HD 92385 ( $P=0\hbox{$.\!\!^{\rm d}$ }54909$, ESA 1997). Among these fast rotating CP stars, HR 7355 is the hottest, most massive and, consequently, the largest thus, we conclude that HR 7355 has the largest equatorial rotational velocity ( $V_{{\rm eq}}\approx 370 \pm 80 \
{\rm km~s}^{-1}$ using the parallax-based MV for a B2 V star and the period) and also the highest ratio of equatorial velocity to critical equatorial velocity $(0.75\pm0.25)$ among all known CP stars. Therefore, we consider HR 7355 to be the most rapid rotator of all known chemically peculiar stars. Because of its rapid rotation, it can serve as an important benchmark for theories describing the influence of rotational mixing on chemical peculiarity.

5.2 Nature of the light variation

Our analysis of the B, Hp, and V light curves demonstrates that all three passbands have the same effective amplitude. This suggests that the variability mechanism at optical wavelengths could be unique (Mikulásek et al. 2008a). We suggest that these light variations are caused by the uneven distribution of optically active, overabundant elements on the surface of the star. Unfortunately, we do not have maps of the abundance distribution or spectrograms suitable for their creation.

We consider the possibility that the broad-band optical light variability in HR 7355 arises from variations in the spectral lines of $He {\sc i}$, $Si {\sc iii}$, and $C {\sc ii}$, as depicted for HR 7355 in Fig. 1 of Rivinius et al. (2008). They compare the line profiles of two spectra obtained at quadrature (HJD  $=2~451~385.507,\ \vartheta=-1882.445(20)$ and HJD  $=2~453~191.879,\ \vartheta=1581.731(15)$, as defined for Eq. (3)); their first spectrum was thus taken shortly after secondary light maximum, while the second was acquired at the secondary light minimum. The only overabundant element whose line intensity reaches maximum in the optical is helium (see the relevant analysis of He-strong HD 37776 in Krticka et al. 2007). Lines of the other ions mentioned above are relatively weak and nearly constant. The photometric effect of the weak emission in H$\alpha$ is also negligible.

We modelled the light curves of HR 7355 with the code described by Krticka et al. (2007) assuming the inclination angle of $\mathit{i}=54^\circ \pm18^\circ$ ( $v \sin i=300\pm 15$ km s-1, Oksala et al. 2010) and found that two circular, helium-rich spots with [He/H]=1.4 and radii $60^\circ$ on opposite hemispheres of HR 7355 reproduce the observed light variations in all three passbands. We therefore conclude that light variations in HR 7355 may be the result of the uneven distribution of helium.

This hypothesis should be tested with Doppler tomography, which we plan to do in the near future.

5.3 Rotational braking and its nature

An increase in the rotational period with $\dot{P}/P_1=2.4(8)\times 10^{-6}~{\rm yr}^{-1}$ is observed at the $3.1\sigma$level of certainty. During the past 20 years, the period has increased by 2.1 s! Unfortunately, no observations were taken in the interval 1994-2000. This adds some uncertainty to our quadratic fit of the O-C residuals. If the period change were real, then HR 7355 would be the fourth known CP star displaying an increase in its rotational period.

The star is similar in several aspects to the most rapidly braking He-strong CP star HD 37776 with a well-determined $\overline{\dot{P}/P}=4.01(17)\times10^{-6}$ year-1(Mikulásek et al. 2008b). The period increase in HR 7355 could also be interpreted as a deceleration of the rotation of its surface layers due to momentum loss by a magnetically confined stellar wind (see also ud-Doula et al. 2009). However, a possible change in the rate of period change for HD 37776 ( $\ddot{P}=-29(13)\times10^{-13}$ d-1) suggests that the process of the rotational braking need not be fully monotonic!

Both stars have strong magnetic fields; that of HR 7355 is approximately dipolar ranging in strength from -2 kG to +2.5 kG (Oksala et al. 2010), while the field of HD 37776 is dominated by a quadruple component (Thompson & Landstreet 1985).However, the main difference between the stars discussed here is their age. While HD 37776 is a very young CP star with an age of around 1 Myr, HR 7355 is about twenty times older (see Sect. 2.1). Because the characteristic braking time of HR 7355 ( $P/\dot{P}\approx 4 \times 10^{5}$ yr) is about fifty times shorter than its stellar age, it seems probable that variations in rotational period are confined to the outer layers of the star. The period variations could be cyclic with a long interval of period increase followed by a rapid period decrease. The period decrease may be caused by greater friction between the surface and inner layers of the star caused by, e.g., reconnection events. The length of this cycle can be roughly estimated as the time when the surface layers lag behind the inner part of the star by one revolution (i.e., $T_{{\rm cycle}}\sim
P~\sqrt{~2/\dot{P}}$). This produces estimates of 35 yr and 45 yr for HR 7355 and HD 37776, respectively. The hypothesis should be tested a few decades hence with additional observations.

6 Conclusions

We have improved the accuracy of the rotational period of the most rapidly rotating CP star HR 7355 to $P=0\hbox{$.\!\!^{\rm d}$ }5214410(4)$ and determined its rate of period change to be $\dot{P}/P=2.4(8)\times 10^{-6}~\rm {yr}^{-1}$. We propose that the period variations could be cyclic on a timescale of a few decades. Observed light variations may be caused by the uneven surface distribution of overabundant helium. HR 7355 remains a very appealing target for continued photometric and spectroscopic observations, as well as for the modelling of its unusual behaviour.


Supports of grants GAAV IAA301630901 and MEB 060807/WTZ CZ-11/2008 are acknowledged. G.W.H. acknowledges long-term support from NASA, NSF, Tennessee State University, and the state of Tennessee through its Centers of Excellence program. E.P. acknowledges support by the City of Vienna (Hochschuljubiläumsstiftung project: H-1930/2008). We thank to dr. D. A. Bohlender for his inspiring refereeing of the article.



... star?[*]
Photometric data are only available in electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr ( or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/511/L7
... GCPD[*]
... CDS[*]
... stars[*]
Copyright ESO 2010

All Figures

\par\includegraphics[width=6.7cm,clip]{13832fg1.eps} \end{figure} Figure 1:

BHpV light curves plotted versus the quadratic phase. $\Box $ - Hipparcos, $\circ $ - Cape Town, $\triangle $ - APT, $\blacklozenge $ - ASAS (each point represents the mean of 30 adjacent ASAS measurements).

Open with DEXTER
In the text

\end{figure} Figure 2:

Time variation in mutual phase shifts of observed light curves in days computed with a linear ephemeris. $\Box $ - Hipparcos, $\diamondsuit $ - ASAS, $\triangle $ - Arizona, and $\circ $ - Cape Town. The Hipparcos and ASAS data sets are divided into two pieces with equal numbers of observations and plotted separately. This division was made for only illustrative reasons and did not affect the period analysis in any way. The quadratic time function is denoted by the solid line, while the 1 $\sigma $-uncertainties are plotted with dashed lines.

Open with DEXTER
In the text

Copyright ESO 2010

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.