{'fraction_non_alphanumeric': 0.03054448871181939, 'fraction_numerical': 0.03380920761398849, 'mean_word_length': 4.308415332647966, 'pattern_counts': {'":': 0, '<': 2, '<?xml version=': 0, '>': 0, 'https://': 1, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Spots on Am stars

We investigate the light variations of 15 Am stars using four years of high-precision photometry from the Kepler spacecraft and an additional 14 Am stars from the K2 Campaign 0 field. We find that most of the Am stars in the Kepler field have light curves characteristic of rotational modulation due to star spots. Of the 29 Am stars observed, 12 are {\delta} Scuti variables and one is a {\gamma} Doradus star. One star is an eclipsing binary and another was found to be a binary from time-delay measurements. Two Am stars show evidence for flares which are unlikely to be due to a cool companion. The fact that 10 out of 29 Am stars are rotational variables and that some may even flare strongly suggests that Am stars possess significant magnetic fields. This is contrary to the current understanding that the enhanced metallicity in these stars is due to diffusion in the absence of a magnetic field. The fact that so many stars are {\delta} Scuti variables is also at odds with the prediction of diffusion theory. We suggest that a viable alternative is that the metal enhancement could arise from accretion.


INTRODUCTION
The Am stars (metallic-lined A stars) are a group of A-and early F-type stars in which the CaII K line is relatively weak, and/or metallic lines relatively strong, compared with the spectral type indicated by the Balmer lines. More detailed analysis shows that Am stars are characterized by an underabundance of Ca (and/or Sc) and/or an over-abundance of Fe and the Fe-group elements. The anomalous abundances in Am and Ap stars are thought to be a result of the interplay between gravitational settling and radiative acceleration of different atomic species. As a result, some elements diffuse upwards and others settle downwards provided there is little mixing. Mixing due to convection or meridional circulation due to rotation destroys the natural segregation of elements due to diffusion that would otherwise occur. Convection in A stars is confined to a very thin sub-surface layer which means that the diffusion process should proceed unhindered in slowly-rotating A stars.
The difference between Ap and Am stars is attributed to the fact that Ap stars have strong magnetic fields, whereas Am stars do not. Aurière et al. (2010) have confirmed that there is at least an order of magnitude difference in magnetic field strength between Ap/Bp stars and Am and HgMn stars. The presence of a large-scale strong dipole magnetic field in Ap stars modifies the diffusion rate in accordance with the geometry of the field and is responsible for the abundance patches in Ap stars. On the other hand, the weak (or perhaps tangled) magnetic fields in Am stars leads to a relatively homogeneous abundance distribution across the stellar surface.
Early calculations predict a much larger overabundance of heavy elements in Am stars than actually observed. Richer et al. (2000) found that it is important to consider the diffusion of heavy elements as well as helium. As the heavy elements settle, enrichment of the deeper iron bump opacity region of the iron-group causes a convective zone to develop in that region which mixes with the convective envelope by overshooting. The result is a dilution of the abundance anomalies, leading to abundances more in line with what is observed. More recent work shows that it is also important to include the destabilizing effect (thermohaline instability) caused by the composition gradient produced by diffusion (µ gradient) .
It is thought that stars in binaries with orbital periods of 2-10 d become Am stars because their rotational velocities, v, have been reduced by tidal interactions below v ≈ 120 km s −1 , a requirement for diffusion to act. Both Am and normal stars occur in binaries with orbital periods between 10-100 d, but why this is the case is not understood. A recent search for eclipsing binaries among Am stars suggests that around 60-70 percent of Am stars are spectroscopic binaries (Smalley et al. 2014), which is consistent with radial velocity studies (Abt & Levy 1985;Carquillat & Prieur 2007). Clearly, a substantial number of Am stars are probably single or in wide binaries and it is no longer clear whether binarity is a requirement for the Am phenomenon.
The effect of diffusion on pulsation of Am stars is to reduce the width of the instability strip, with the blue edge shifting towards the red edge, eventually leading to the disappearance of instability when helium is sufficiently depleted from the HeII ionization zone . The amount of driving is related to the amount of He still present in the driving zone. The cooler the star, the deeper the mixing and the greater the He abundance in the driving zone. Models therefore predict that Am stars should not pulsate unless they are sufficiently cool. More recent calculations suggest that rotational mixing is important in maintaining the He abundance in the driving region (Talon et al. 2006).
These ideas can be tested if we can obtain sufficiently accurate physical parameters to place pulsating and nonpulsating Am stars in the H-R diagram. Smalley et al. (2011) studied over 1600 Am stars and found that around 200 Am stars are pulsating δ Sct and γ Dor stars. It seems that Am stars, in general, seem to pulsate at somewhat lower amplitudes than normal δ Sct stars. However, Smalley et al. (2011) could find no real difference between the location of pulsating Am stars and δ Sct stars in the instability strip, except for a slight tendency of pulsating Am stars to be cooler than normal δ Sct stars. Catanzaro & Balona (2012) also concluded that there is a tendency for pulsating Am stars to be cooler than normal δ Sct stars and that the percentage of pulsating Am stars is about the same as the percentage of δ Sct stars in the instability strip. Similar conclusions were obtained by Balona et al. (2011c). It is clear that the prediction of diffusion theory regarding pulsation needs to be re-examined. The diffusion process, in fact, seems to play only a minor role in determining the pulsational stability of A and early F stars. Even Ap stars pulsate as δ Sct stars at low amplitudes (Balona et al. 2011a).
Recent high-precision photometric observations of A stars using Kepler have further undermined the general consensus that these stars have stable envelopes where diffusion has a strong competitive advantage. It seems that about 2 percent of A stars observed by Kepler show flares (Balona 2012(Balona , 2014c. It can be easily demonstrated that the flares are associated with the A star and not with a presumed cool companion. Furthermore, the Kepler data suggest that at least 40 percent of A stars are rotational variables (Balona 2013), presumably due to starspots. Since it is generally believed that Am stars do not have magnetic fields, the incidence of rotational modulation should be absent or small in Am stars compared to normal A stars.
The previous study of Kepler Am stars (Balona et al. 2011c) used only data for approximately the first 50 d and was mainly aimed at detecting δ Sct pulsations in these stars. The low-amplitude light variations due to rotational modulation were suspected but could not be studied. There are now four years of almost continuous high-precision photometry for these stars, making it easier to detect rotational modulation. In this paper we investigate rotational modulation in Am stars using Kepler data. Our aim is to establish whether the incidence of rotational modulation or other activity is lower in Am stars than in normal A stars. The answer might explain why there are many slowly-rotating A stars which are not Am stars and provide further clues as to the nature of the Am phenomenon.

THE DATA
There are 15 known Am stars in the Kepler field (Balona et al. 2011c). We also include 14 Am stars observed by the Kepler K2 Campaign 0 mission. These data are far less favourable to detecting rotational modulation because only about 45 d of data are available and the photometric precision is severely compromised by the telescope motion. However, The K2 data are useful in detecting δ Sct pulsations in these stars.
For the vast majority of stars, Kepler photometry is available only in long-cadence (LC) mode which consist of practically uninterrupted 30-min exposures covering a time span of over three years. For a few thousand stars, including all the Am stars discussed here, short-cadence (SC) 1min exposures are also available, but these usually cover only one or two months. Characteristics of LC data are described in Jenkins et al. (2010b) while Gilliland et al. (2010) describe the characteristics of SC data. The science processing pipeline is described by Jenkins et al. (2010a). The Kepler data are divided into quarters of approximately 90d duration which are numbered Q0, Q1, Q2, etc., Q0 being the initial 10-d commissioning run. In this study we use all available LC data (Q0 -Q17) from JD 2454953 to JD 2456424 (1471 d). These data are publicly available on the Barbara A. Mikulski Archive for Space Telescopes (MAST, archive.stsci.edu).
The light curve files contain simple aperture photometry (SAP) values and a more processed version of SAP with artifact mitigation included called Pre-search Data Conditioning (PDC) flux. The SAP photometry contains discontinuities, outliers, systematic trends and other instrumental signatures. The PDC tries to remove these errors while preserving astrophysically interesting signals by using a subset of highly correlated and quiet stars to generate a cotrending basis vector set. The PDC flux uses the best fit that simultaneously removes systematic effects while reducing the signal distortion and noise injection. Details of how this is accomplished are described in Stumpe et al. (2012) and Smith et al. (2012). Because we are interested in lowfrequency variations, we have used the PDC fluxes in this analysis.
The Kepler SAP data contains drifts and jumps. There is a zero-point difference between data in different quarters and the instrumental drift varies from quarter to quarter. Although the PDC data are designed to remove these effects, it is possible that residual long-term instrumental effects are still present at a low level. To understand how correction of data may effect the low-frequency end of the periodogram, we performed a simple simulation where we applied our own correction to the SAP data of a few stars which appeared to be the least variable among 20000 Kepler stars. To correct the SAP data we fitted and removed polynomials of low degree (typically degree 3 or 5) to individual quarters. The polynomials do not join smoothly between any two quarters, but we ignored these small jumps. To these data we added sinusoidal variations comprising of 100 frequencies of equal amplitudes and random phases. These frequencies start at zero frequency and are equally spaced by 0.1 d −1 . The periodogram recovers all frequencies, including the lowest frequency of 0.1 d −1 . The same simulation using a signal of 0.02 d −1 recovers a lowest frequency of 0.08 d −1 . From this simulation, we may deduce that frequencies higher than about 0.08 d −1 are not affected by the crude correction procedure. The corrections used to obtain the PDC data are, of course, more sophisticated and we can safely assume that frequencies higher than 0.08 d −1 are easily recovered.
The light curves from the extended Kepler K2 Campaign 0 mission were obtained from the raw pixel data files using a suitable aperture. Our own software, keplc, was used for this purpose. For fields where there were few or faint stars in close proximity to the target, we used the largest possible aperture in order to minimize the instrumental light variation due to image motion. The telescope motion can be found by tracking the centroid of the star. There is a 4.08 d −1 periodic component in the image motion which is reflected in the light curve. We therefore removed this frequency and all its harmonics from the light curve so as to enhance the actual stellar signal. Long-term drifts and jumps were also removed. The periodogram of the processed light curve was then calculated. Due to image motion and long-term drifts, the periodogram noise level for the K2 data increases strongly towards zero frequency. As a result, it is generally not possible to detect low frequencies in the star unless these are of high amplitude.

KEPLER LIGHT CURVES
A list of known Am stars in the Kepler and K2 Campaign 0 fields is shown in Table 1. The identifications are from the Kepler Input Catalogue (KIC, Brown et al. (2011)) and from the K2 Ecliptic Plane Input Catalog (EPIC, https://archive.stsci.edu/k2/KSCI-19082-002.pdf). The spectral classification for Am stars in the Kepler field are mostly quite recent and detailed metal abundances are available for most of the stars. The Am nature of these stars is hardly in doubt. The spectral classification for Am stars in the K2 field are far less secure. All of them are listed as Am stars by Renson & Manfroid (2009), but the original reference is often missing. In Table 1 we have referenced what we believe to be the original source, but sometimes the Am nature is not even mentioned. Part of the problem may be in the way the recognition of the Am class has evolved with time.
Examination of the light curves shows that of the 29 Am stars, there are 12 δ Sct variables, 1 γ Dor variable and 2 detached eclipsing binaries. Most of the others have light curves similar to those generally attributed to rotational modulation caused by starspots. Periodograms of the δ Sct Am stars in the Kepler field are shown in Fig. 1. Periodograms of the other stars in the Kepler field are shown in Fig. 2.
It is evident from Fig. 2 that highly significant, coherent, low-frequency variations and harmonics are present in all six non-δ Sct stars. In KIC 3836439 this can be attributed to proximity effects in a binary. Luyten (1936) found the star to be a spectroscopic binary with P orb = 1.54039 d, semi-amplitude K = 88.7 ± 2.8 km s −1 and a circular orbit. Batten et al. (1978)  f (m) = 0.112 and a1 = 1.88 × 10 6 km = 2.68R⊙. Kepler observations confirm that this is an eclipsing binary with P = 1.540407 d in a detached Algol-type system with sin i = 0.96 . The ratio of effective temperatures T2/T1 = 0.758 and the sum of radii R1 + R2 = 0.35a where a is the semi-major axis. Slightly different results are given by Prša et al. (2011): sin i = 0.95, T2/T1 = 0.48 and R1 + R2 = 0.37a. It is possible that the Am classification is a result of a composite spectrum. Table 1. List of Am stars in the Kepler field (top) and the K2 Campaign 0 field (bottom). The stars are identified by their KIC number and their EPIC number respectively while the second column is the HD number. The spectral classification is given in the third column. The apparent magnitude and the derived effective temperature, T eff , are mostly from the KIC except where noted. The surface gravity, log g and projected rotation velocity, v sin i, are from Uytterhoeven et al. (2011) except where noted. The variable star classification is estimated from the light curve and periodogram. The rotation period of the star (Prot in days) and the rotational light amplitude (A in ppm) is estimated from the periodogram. References are as follows: 1 - Abt (1981), 2 - Abt (1984), 3 - Ammons et al. (2006), 4 -Antoci et al. (2011), 5 -Bertaud (1960), 6 -Catanzaro et al. (2011), 7 -Catanzaro & Ripepi (2014), 8 -Catanzaro et al. (2015), 9 -Floquet (1970, 10 -Floquet (1975), 11 - (Halbedel 1985), 12 -McCuskey (1967), 13 -McCuskey (1959, 14 - Renson & Manfroid (2009). In the other stars, the low-frequency variations are unlikely to be due to proximity effects in a binary because the amplitudes often vary and the periodogram lines attributed to rotation are broadened in some stars. This is to be expected in a rotational variable because spots change in size and intensity with time. Differential rotation, which appears to be a general feature in A stars (Balona 2013), will also cause such broadening. As an example of rotational modulation, Fig. 3 shows part of the long-cadence PDC light curve of KIC 8323104. Quasi-regular light variability is present with a peak-to-peak amplitude of about 40 ppm and a timescale of around 6 d. This cannot be a result of proximity effects in a binary because of the irregularity of the variations. The most plausible explanation is that it is a result of rotational modulation due to stellar activity.
The second-last column of Table 1 shows the rotation period estimated from the periodogram. In the case of δ Sct stars there are often many other low-frequency peaks which are always present in any δ Sct star. In this case we identified the rotation peak because of its relatively high ampli-tude and the presence of an harmonic (Balona 2013). The presumed rotation period decreases with v sin i, as might be expected.
KIC 8703413 was first reported as an Am star by Mendoza (1974) and subsequently classified as kA2mF0 by Floquet (1975). The star was confirmed as Am by Catanzaro & Ripepi (2014). As can be seen in Fig. 4, the Kepler light curve shows distinct beating and traveling wave characteristic of star spots Balona 2013). The periodogram shown in the middle panel has a main peak at f1 = 0.1532 d −1 and its harmonics. In addition, there is another peak f2 = 0.1404 d −1 which looks to be a bit broader. It is not clear how this pattern may be interpreted, but it seems reasonable to suppose that either f1 or f2 could be the rotation frequency. Perhaps both frequencies could be ascribed to starspots if the star is in differential rotation. No other frequencies are present. The projected rotational velocity, v sin i = 15 ± 2 km s −1 , is compatible with the rotation frequency, giving a radius R sin i ≈ 1.96R⊙. There are seven radial velocities in the literature (Fehrenbach et al. 1997), which is not enough for an analysis. The velocity range is quite large and it is possible that the star may be a binary. KIC 8703413 is a flare star (Balona 2012) as can be seen in the bottom panel of Fig. 4. The flare intensity is about 0.7 ppt. The A star has a luminosity of about 100 times larger than a K dwarf, which means that if the flare is attributed to an unseen K dwarf companion, its relative intensity would have to be about 70 ppt. The mean flare intensities in Kepler field K-dwarfs are about 3.3 ppt (Walkowicz et al. 2011), which means that this particular flare would have to be of unusual intensity. In general, flares in A stars cannot be attributed to a cool companion since the average flare intensity is about two orders of magnitude larger than that in a typical K or M dwarf.
For KIC 3429637, Catanzaro et al. (2011) obtained T eff = 7300 ± 200 K, log g = 3.16 ± 0.25 and v sin i = 50 ± 5 km s −1 , values later confirmed by Murphy et al. (2012). The remarkable amplitude increase in the δ Sct pulsations of this star has been discussed by Murphy et al. (2012). Since that time, additional Kepler data show a more complex amplitude variation for the frequencies of largest amplitude (Fig. 5,  f3 = 10.936 d −1 show steady amplitude changes, the steady amplitude increase for f2 = 12.472 d −1 reverses and decreases sharply around JD 2455950. The latest observations indicate that the amplitude of f2 is, once again, increasing. The timescale seems to be too short for an evolutionary explanation, as suggested by Murphy et al. (2012). KIC 5219533 has a very rich δ Sct frequency spectrum, including a large number of modes in the γ Dor range. The projected rotational velocity is v sin i = 115 km s −1 . The principal peak at f1 = 10.2853 d −1 is clearly variable in frequency. The phase, φ, varies smoothly in the range 2.0 < φ < 4.0 radians over 400 d which corresponds to approximately 0.0008 d −1 .
A recent abundance analysis of KIC 9117875 by Catanzaro & Ripepi (2014) indicates T eff = 7400 ± 150 K, log g = 3.6 ± 0.1 dex, and v sin i = 58 ± 6 km s −1 . Both Ca and Sc show underabundances of about 1 dex, while the heavy elements are all overabundant: 0.5 dex for iron-peak elements and about 2 dex for Ba. The Am nature of this star is thus confirmed. The Kepler light curve shows beating and low frequencies resembling a γ Dor star and it was classified as such by Uytterhoeven et al. (2011). However, a more detailed examination of the periodogram shows a broad feature at about f1 = 0.672 d −1 with a very sharp neighbouring peak at f2 = 0.7157 d −1 . Several harmonics of f2 are visible, which is quite unlike the structure found in typical γ Dor stars. In fact, this mysterious structure is found in a considerable fraction of A stars, but is unexplained (Balona 2013). Balona (2014c) suggests that the sharp feature might be explained by a synchronously orbiting planet, not necessarily transiting. Assuming that f1 is the rotation frequency and using the above projected rotational velocity leads to a radius R sin i ≈ 1.69R⊙, which is reasonable for an A star.
There is a further mystery regarding this star in that it appears to flare. In Fig. 7 we show examples of possible flares in this star in the Kepler LC data. It should be noted that LC data is not ideal to detect flares since the exposure time of 30 min is a considerable fraction of the duration of the flare. Short-cadence data are available for this star, but only for about 10 d during Q0.
KIC 9204718 was confirmed as an Am star by Catanzaro & Ripepi (2014) who find T eff = 7600 ± 150 K, log g = 3.8 ± 0.1 dex and v sin i = 27 ± 3 km s −1 . There are clearly several low-frequency peaks in the periodogram, as shown in Fig. 8. The largest peak is at 0.1145 d −1 and its harmonic. These peaks are both rather broad suggesting amplitude or frequency modulation; in fact, the light variation with this frequency is clearly visible in the light curve (Fig. 8). This variation is suggestive of a starspot, in which case 0.1145 d −1 would be the rotation frequency. This is reasonable since we measured v sin i = 27 ± 3 km s −1 .
Even more interesting is the presence of several harmonics of 0.5474 d −1 . If we remove the 0.1145 d −1 and the δ Sct frequencies from the data, what is left is the 0.5474 d −1 variation and its harmonics. By binning the light curve and averaging the data in each bin, we obtained the phased light curve shown in the bottom panel of Fig. 8. There is no simple explanation for this variation. It cannot be the rotational frequency, since this has already been identified. The light curve resembles that of a contact binary, but if so it implies that the rotational period is not synchronized with the orbital period. There is a great deal of uncertainty regarding the nature of KIC 9272082. Macrae (1952) noted its possible peculiar spectrum, but did not give any details. Bertaud (1960) gives a classification of A5m, but Floquet (1970) classified it as a normal A7 star. The star is included as an Am star in the Renson & Manfroid (2009) catalogue. Catanzaro & Ripepi (2014) found it to be an A4 main sequence star with T eff = 8500±200 K, log g = 3.9±0.1 dex and v sin i = 75±7 km s −1 . Most of the elements in this star are somewhat overabundant by about 0.5-1.8 dex. However, Si, Cr, Ni, and Ba have normal abundances. Although the abundance pattern is not typical, it does not conform to that expected for Am stars.
The light curve (top panel of Fig. 9) shows distinct variations which are confirmed by the periodogram (middle panel). There is a constant frequency f1 = 1.03058 d −1 and its harmonic, suggestive of rotational modulation. The corrected LC data phased with this frequency shows a roughly sinusoidal variation (Fig. 9), which we assume to be rotational modulation.

BINARY MOTION FROM TIME DELAY
It is generally thought that most Am stars are binaries with orbital periods in the range 1-10 d (Abt 1967). It is of interest to determine whether any of the stars in our sample are spectroscopic binaries. A pulsating star in a binary system acts as a clock. If the star is a member of a binary system, its distance will vary as it orbits the barycentre. The changing distance leads to a variable phase of pulsation. The method has been discussed by Shibahashi & Kurtz (2012) using data transformed to the frequency domain. More recently, Murphy et al. (2014) has used the more direct analysis of phase variation. Application to δ Sct stars is difficult because of the many close frequencies which distort the phase variation. By direct minimization of the photometric error, it is possible to construct a diagram similar to the periodogram where the orbital semi-major axis is shown as a function of the orbital period (the binarogram). This technique allows relatively quick detection of binary candidates (Balona 2014a). KIC 3429637 has been discussed above in the context of its δ Sct pulsations. f1, shows a clear periodic variation with P = 704.7 d. There is a rather sudden change in phase of f2 (corresponding to a change in frequency) at a time soon after the sharp decrease in amplitude shown in the top panel. Close inspection reveals that the 704.7-d period is also present in the phase variation of f2. Short-cadence observations are available for this star, so there is no ambiguity regarding the frequencies.
We applied the time-delay algorithm described above to both f1 and f2. Because there is a strong sudden decrease of amplitude and frequency around JD 2455950, we decided to use only the data up to this date. There is still a longterm trend which was removed by fitting a Fourier series combined with a polynomial of third degree. we obtain the orbital solution shown in Table 4 and shown in the bottom panel of Fig. 5. For KIC 11445913 short cadence observations show that the dominant peaks are at f1 = 31.5578, f2 = 25.3771, f3 = 22.1330 and f4 = 37.8204 d −1 . As usual there are quite a number of low-frequency peaks, but none that could be identified as due to starspots. The value v sin i = 51 ± 1 km s −1 (Balona et al. 2011c). Time delay analysis of f1, f2 and f3 shows that they have the same phase variation, but the period, if it exists, is considerably longer than the almost 4 yr duration of the Kepler data. Fig. 10 shows the resulting radial velocities derived from the first derivative of the time delay variation. This star might be a binary with a period of around 10 yr or longer. Table 2. Binary elements of KIC 3429637 from the time delay of the Kepler light curve. The orbital period P orb , half-range of primary radial velocity variation K 1 , eccentricity e, longitude of periastron ω, JD of periastron Tper, projected semi-major axis of the primary a 1 sin i, and mass function f (M ).

K2 LIGHT CURVES
Owing to the failure of two of the four reaction wheels, the Kepler spacecraft has lost its pointing ability. The telescope drift is corrected at regular intervals, with the result that there is substantial image drift and a large increase in photometric error from about 20 ppm for the original Kepler field to about 300 ppm for the K2 field. Here we present results for Am stars observed in the first K2 scientific mission, Campaign 0. The observations occur over a timespan of about 77 d, but with a large gap of 26.5 d. The total ontarget time is about 47 d. This rather short time and the much increased scatter severely compromises detection of significant frequencies less than about 3 d −1 . Many of these stars were observed by the STEREO spacecraft. No significant variability could be found for EPIC 202059291, 202059336, 202061329, 202061333, 202061336, 202061338 and 202061363 from STEREO observations Paunzen et al. (2013). The remarkably low effective temperature for EPIC 202061363 reported in Ammons et al. (2006) is confirmed by our own analysis of available multicolour photometry. STEREO observations of EPIC 202061349 also do not show any significant variability (Wraight et al. 2012). Of these stars, we detected variability in EPIC 202059291, 202061333 and EPIC 202061338 (see below). Some of the stars are δ Sct variables. The periodograms of the δ Sct stars are shown in Fig. 11. The effective tem- perature of EPIC 202062436 was obtained through our own analysis of the multicolour photometry. EPIC 202061353 is a known δ Sct star from SuperWasp photometry . This is confirmed by the K2 data (Fig. 11). The main feature is a triplet of lowfrequency peaks at 1.452, 0.980 and 2.000 d −1 which are not quite equally spaced.
The periodogram of EPIC 202059291 shows only two significant peaks. Considering the very low amplitude, the variability is best understood as rotational modulation with period P = 2.984 d, but it could also be interpreted as a contact binary. One of us (GC) obtained a single spectrum of this star which appears to be a double-lined spectroscopic binary. It can probably be ruled out as an Am star. The periodogram and phased light curve is shown in Fig. 12.
The light curve of EPIC 202062447 (Fig. 13) shows that it is a γ Dor of the ASYM type (Balona et al. 2011b). The periodogram shows dominant periods of 0.912 d and 0.773 d.
The field of EPIC 202062450 contains two stars of almost equal brightness separated by 5 arcsec. The brighter star, FI Ori, is a detached eclipsing binary with P = 4.4478 d and spectral type K0III/IV (Halbedel 1985). The eclipses are very easily seen in the K2 combined light of the two stars. The fainter star, HD 256320 B has a spectral type kA5hA8mF5 (Halbedel 1985). We constructed an aperture which only included this star, but could not find any evidence of periodic variation with an amplitude exceeding 5 ppt.

DISCUSSION AND CONCLUSION
It is evident from Table 1 that most Am stars in the Kepler field exhibit rotational modulation. An idea of the spot size may be gained from the light amplitude of the modulation shown in Table 1. There are only 10 Am stars in the Kepler field known to be rotational variables and the amplitude range is large. The mean amplitude is 93±25 ppm. The mean rotational amplitude from 875 A stars in Balona (2013) is 541 ± 86 ppm. This seems to suggest that, in general, spots on Am stars are smaller than on normal A stars, but it is not possible to draw a definite conclusion owing to the small sample.
One or more flares appear to be seen in two confirmed Am stars, KIC 8703413 and KIC 9117875. The relatively large flare intensities imply that if the origin of the flare is a cool companion, the flare must be about two orders of magnitude more intense than on an isolated cool star owing to the much larger luminosity of the Am star. While this is not impossible, it seems very improbable. It is clear from a study of a much larger number of A stars that the flares cannot all be attributed to a cool companion (Balona 2012(Balona , 2013.
The presence of starspots and flares on Am stars suggests that Am stars may have significant magnetic fields similar in strength to normal A stars. The magnetic fields in both Am and normal A stars are sufficiently intense to form starspots and flares. Hence the idea that the superficial metal enrichment in Am stars is a result of the effects of diffusion and gravitational settling in the absence of a magnetic field may need to be revised. Furthermore, the relative number of pulsating Am stars is similar to the relative number of δ Sct stars among the A star population. Am stars which are δ Sct variables occupy the same region of the instability strip as normal δ Sct stars . All these observations are contrary to the view that Am stars are formed by diffusion of elements in the absence of a magnetic field and that δ Sct pulsation in Am stars only occurs near the red edge of the instability strip. Given these facts, one has to conclude either that the diffusion model for Am stars needs to be revised or that it is incorrect.
Thermohaline convection takes place when a layer of enhanced metal abundance is formed above layers of lighter mean molecular weight. This happens in the case of accretion of predominantly metal-rich material (Vauclair 2004), and also when such a layer is formed due to atomic diffusion. This important effect has only recently been studied. Théado et al. (2009) finds that accumulation of heavy elements due to atomic diffusion is attenuated when thermohaline mixing is taken into account, though not completely suppressed. The inclusion of the thermohaline instability in A-star models needs to be more fully investigated to see whether this can reconcile the diffusion model of Am stars with observations.
The presence of pulsations in Am stars is only one symptom of a more general problem which has been encountered since the advent of Kepler photometry. For example, it is not understood why less than half of stars in the δ Sct instability strip pulsate (Balona & Dziembowski 2011), or why lowfrequency pulsations are a general feature of all δ Sct stars (Balona 2014b). Further exploration of pulsational driving needs to be undertaken. For example, coherent pulsations may be driven by turbulent pressure in the hydrogen ionization zone (Antoci 2014).
It is well established that stars hosting planets are metal rich, but the reason for such metal enhancement is still a subject for debate. The enhancement may be primordial, the result of accretion or both. The argument against an accretion origin is related to the mass of the outer convective zones, which varies by more than one order of magnitude among the considered stars, while the observed over-abundances of metals are about the same. Vauclair (2004) found that thermohaline convection mixes the accreted matter and leaves only a very small µ-gradient at the end of the mixing process. The remaining µ gradient is too small to account for the observed overmetallicity (see also ). It should be noted, however, that all these calculations are confined to solar-type stars. The effect of thermohaline convection on intermediate mass A stars has not yet been determined.
Considering the problems encountered by standard diffusion theory in explaining the properties of Am stars, it may be appropriate to again consider the possibility that the Am phenomenon may be a combination of diffusion with accretion. For example, Böhm-Vitense (2006) has argued that accretion of interstellar material by A stars with tangled magnetic fields, which are weaker than those in peculiar A (Ap) stars, has the best chance of explaining the main characteristics of the peculiar heavy-element abundances in Am star photospheres. Another look at accretion of close-in planets should also be made. We know, for example, that accretion of comet-sized bodies is very likely occurring in the A6V star β Pic which is surrounded by a debris disk (Beust et al. 1996). There is no reason to suppose that this star is unique and we can expect infalling material, even of planets, to occur. As Pinsonneault et al. (2001) has pointed out, the mass of the outer convective zone decreases rapidly for hot stars, being very thin in A stars. Therefore one may expect infalling material, mixed only in the thin outer convective zone, to have a substantial effect on the observed metallicity. For example, the mass above the base of the He II convective zone in a typical mid-A star is less than the mass of the Earth. Since the original fraction of metals in the star is only about 2 percent of the mass of the zone, even a body much smaller than the Earth could increase the apparent metal abundance by a factor of two or more. It would be surprising if no such impacts ever occurred. Therefore one should not be surprised by the peculiar metal abundances in some A stars.
The chemically peculiar λ Boo stars have an unusually low superficial abundance of iron peak elements. It is thought that this may be a result of accretion of interstellar material as the star travels through a diffuse interstellar cloud (Kamp et al. 2008). However, the abundance pattern in λ Boo stars is very different from Am stars. This may reflect differences in abundance between the metal-poor interstellar medium and the metal-rich composition of infalling small planetary bodies.