Over the course of its history, the Milky Way has ingested multiple smaller satellite galaxies. While these accreted stellar populations can be forensically identified as kinematically distinct structures within the Galaxy, it is difficult in general to precisely date the age at which any one merger occurred. Recent results have revealed a population of stars that were accreted via the collision of a dwarf galaxy, called \textit{Gaia}-Enceladus, leading to a substantial pollution of the chemical and dynamical properties of the Milky Way. Here, we identify the very bright, naked-eye star $\nu$\,Indi as a probe of the age of the early in situ population of the Galaxy. We combine asteroseismic, spectroscopic, astrometric, and kinematic observations to show that this metal-poor, alpha-element-rich star was an indigenous member of the halo, and we measure its age to be $11.0 \pm 0.7$ (stat) $\pm 0.8$ (sys)$\,\rm Gyr$. The star bears hallmarks consistent with it having been kinematically heated by the \textit{Gaia}-Enceladus collision. Its age implies that the earliest the merger could have begun was 11.6 and 13.2 Gyr ago at 68 and 95% confidence, respectively. Input from computations based on hierarchical cosmological models tightens (i.e. reduces) slightly the above limits.

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 apply the Rayleigh-Ritz Variational Scheme for studying the transverse oscillations of a magnetic flux tube. The flux tube is considered in low-$\beta$ solar coronal condition. The perturbations decomposed into irrotational and solenoidal components and MHD equations are reduced as a matrix eigenvalue problem. In the case of longitudinally stratified thin flux tube, the fundamental and higher order kink frequencies are computed. The results are in good agreements with previous studies.

The distributions of solar flare energies and waiting times have not been described simultaneously by a single physical model, yet. In this research, we investigate whether recent avalanche models can describe the distributions for both the released energies and waiting times of flares in an active region. Flaring events are simulated using the modified Lu and Hamilton model (Lu and Hamilton (1991), ApJ, 380, 89) and also the optimized model (Farhang et al. (2018), ApJ, 859, 41). Applying a probability balance equation approach, we study the statistics of the simulated flaring events and investigate the origin of the observed power law in the flare frequency-size distribution. The results indicate that the power law originates in the distribution of transition rates (the distribution of the probabilities of transitions between different energies) rather than the distribution of the energy of the active region. It is also observed that the waiting-time distribution of simulated flaring events follows a q-exponential function which approximates a simple Poisson distribution.

The collective transverse coronal loop oscillations seem to be detected in the observational studies. In this regard, Luna et al. (2009, ApJ, 692, 1582) modeled the collective kinklike normal modes of several cylindrical loops system using the T-matrix theory. This paper investigates the effects of longitudinal density stratification along the loop axis, on the collective kinklike modes of system of coronal loops. The coronal loops system is modeled as cylinders of parallel flux tubes, with two ends of each loop at the dense photosphere. The flux tubes are considered as uniform magnetic fields, with stratified density along the loop axis which changes discontinuously at the lateral surface of each cylinder. The MHD equations are reduced to solve a set of two coupled dispersion relations for frequencies and wave numbers, in the presence of stratification parameter. The fundamental and first overtone frequencies and longitudinal wave numbers are computed. The previous results are verified for unstratified coronal loops system. Finally we would conclude that increased longitudinal density stratification parameter will result in increase of the frequencies. The frequencies ratios, first overtones to fundamentals, are very sensitive functions of density scale height parameter. Therefore, stratification should be included in dynamics of coronal loops systems. For the unstratified coronal loops system, these ratios are the same as monoloop ones.

Modeling fast solar wind based on the kinetic theory is an important task for scientists. In this paper, we present a two-fluid model for fast solar wind with anisotropic Kappa-Maxwellian electrons and Bi-Maxwellian protons. In the simulation, the energy exchange between the plasma particles and low-frequency Alfv\'en waves is considered. A set of eleven coupled equations is derived by applying the zeroth- to fourth-order moments of the Vlasov equation and the modified electromagnetic Maxwell equations. A characteristic of the Kappa distribution (indicated by $\kappa$ index) is explicit in the equation for the parallel component of the electron heat flux (parallel to the ambient magnetic field line) and differs from the equation derived for the proton heat flux due to the different nature of the distributions. Within the large $\kappa$ index, the equations for the two-fluid model tend to the equations obtained by the Maxwellian distribution. Using an iterated Crank-Nicolson method, the coupled equations are numerically solved for the fast solar wind conditions. We show that at (0.3 - 1) AU from the Sun, the electron density, components of temperature, and components of heat flux follow the power-law behavior. We also showed that near the Earth, the flow speed (electron or proton) increases with decreasing $\kappa$. We concluded that applying the small $\kappa$ index (the non-Maxwellian distribution), the extraordinary nature of the solar atmosphere, with its temperature of several million kelvin temperature for electrons, has been captured.

We present the first results from the Transiting Exoplanet Survey Satellite (TESS) on the rotational and pulsational variability of magnetic chemically peculiar A-type stars. We analyse TESS 2-min cadence data from sectors 1 and 2 on a sample of 83 stars. Five new rapidly oscillating Ap (roAp) stars are announced. One of these pulsates with periods around 4.7 min, making it the shortest period roAp star known to date. Four out of the five new roAp stars are multiperiodic. Three of these, and the singly-periodic one show the presence of rotational mode splitting. Individual frequencies are provided in all cases. In addition, seven previously known roAp stars are analysed. Additional modes of oscillation are found in some stars, while in others we are able to distinguish the true pulsations from possible aliases present in the ground-based data. We find that the pulsation amplitude in the TESS filter is typically a factor 6 smaller than that in the $B$ filter which is usually used for ground-based observations. For four roAp stars we set constraints on the inclination angle and magnetic obliquity, through the application of the oblique pulsator model. We also confirm the absence of roAp-type pulsations down to amplitude limits of 6 and 13 micromag, respectively, in two of the best characterised non-oscillating Ap (noAp) stars. We announce 27 new rotational variables along with their rotation periods, and provide different rotation periods for seven other stars. Finally, we discuss how these results challenge state-of-the-art pulsation models for roAp stars.

Solar coronal loops represent the variety of fast, intermediate, and slow normal mode oscillations. In this study, the transverse oscillations of the loops with a few-minutes period and also with damping caused by the resonant absorption were analyzed using extreme ultraviolet (EUV) images of the Sun. We employed the 171 $\AA$ data recorded by Solar Dynamic Observatory (SDO)/Atmospheric Imaging Assembly (AIA) to analyze the parameters of coronal loop oscillations such as period, damping time, loop length, and loop width. For the loop observed on 11 October 2013, the period and the damping of this loop are obtained to be 19 and 70 minutes, respectively. The damping quality, the ratio of the damping time to the period, is computed about 3.6. The period and damping time for the extracted loop recorded on 22 January 2013 are about 81 and 6.79 minutes, respectively. The damping quality is also computed as 12. It can be concluded that the damping of the transverse oscillations of the loops is in the strong damping regime, so resonant absorption would be the main reason for the damping.

Aims: The statistics of the photospheric granulation pattern are investigated using continuum images observed by Solar Dynamic Observatory (SDO)/Helioseismic and Magnetic Imager (HMI) taken at 6713~\AA. Methods: The supergranular boundaries can be extracted by tracking photospheric velocity plasma flows. The local ball-tracking method is employed to apply on the HMI data gathered over the years 2011-2015 to estimate the boundaries of the cells. The edge sharpening techniques are exerted on the output of ball-tracking to precisely identify the cells borders. To study the fractal dimensionality (FD) of supergranulation, the box counting method is used. Results: We found that both the size and eccentricity follow the log-normal distributions with peak values about 330 Mm$^2$ and 0.85, respectively. The five-year mean value of the cells number appeared in half-hour sequences is obtained to be about 60 $\pm$ 6 within an area of $350^{\prime\prime}\times350^{\prime\prime}$. The cells orientation distribution presents the power-law behavior. Conclusions: The orientation of supergranular cells ($O$) and their size ($S$) follows a power-law function as $|O| \propto S^{9.5}$. We found that the non-roundish cells with smaller and larger sizes than 600 Mm$^2$ are aligned and perpendicular with the solar rotational velocity on the photosphere, respectively. The FD analysis shows that the supergranular cells form the self-similar patterns.

Solar flares are an abrupt release of magnetic energy in the Sun's atmosphere due to reconnection of the coronal magnetic field. This occurs in response to turbulent flows at the photosphere which twist the coronal field. Similar to earthquakes, solar flares represent the behavior of a complex system, and expectedly their energy distribution follows a power law. We present a statistical model based on the principle of minimum energy in a coronal loop undergoing magnetic reconnection, which is described as an avalanche process. We show the distribution of peaks for the flaring events in this self-organized critical system is scale-free. The obtained power law index 1.84 $ \pm $ 0.02 for the peaks is in good agreement with satellite observations of soft X-ray flares. The principle of minimum energy can be applied for general avalanche models to describe many other phenomena.

We investigate the characteristics of the solar flares complex network. The limited predictability, non-linearity, and self-organized criticality of the flares allow us to study systems of flares in the field of the complex systems. Both the occurrence time and the location of flares detected from January 1, 2006 to July 21, 2016 are used to design the growing flares network. The solar surface is divided into cells with equal areas. The cells, which include flare(s), are considered as nodes of the network. The related links are equivalent to sympathetic flaring. The extracted features present that the network of flares follows quantitative measures of complexity. The power-law nature of the connectivity distribution with a degree exponent greater than three reveals that flares form a scale-free and small-world network. The great value of the clustering coefficient, small characteristic path length, and slowly change of the diameter are all characteristics of the flares network. We show that the degree correlation of the flares network has the characteristics of a disassortative network. About 11% of the large energetic flares (M and X types in GOES classification) that occurred in the network hubs cover 3% of the solar surface.

Magnetic elements of the solar surface are studied in magnetograms recorded with the high-resolution Solar Dynamics Observatory / Helioseismic and Magnetic Imager . To extract some statistical and physical properties of these elements (e.g., filling factors, magnetic flux, size, lifetimes), the Yet Another Feature Tracking Algorithm (YAFTA), a region-based method, is employed. An area with 400$^{\prime\prime}\times$400$^{\prime\prime}$ was selected to investigate the magnetic characteristics during the year 2011. The correlation coefficient between filling factors of negative and positive polarities is 0.51. A broken power law fit was applied to the frequency distribution of size and flux. Exponents of the power-law distributions for sizes smaller and greater than 16 arcsec$^2$ were found to be -2.24 and -4.04, respectively. The exponents of power$-$law distributions for fluxes smaller and greater than 2.63$\times$10$^{19}$ Mx were found to be -2.11 and -2.51, respectively. The relationship between the size ($S$) and flux ($F$) of elements can be expressed by a power-law behavior in the form of $S\propto F~^{0.69}$. The lifetime and its relationship with the flux and size of quiet-Sun (QS) elements are studied during three days. The code detected patches with lifetimes of about 15 hours, which we call long-duration events. It is found that more than 95\% of the magnetic elements have lifetimes of less than 100 minutes. About 0.05\% of the elements were found with lifetimes of more than 6 hours. The relationships between the size (S), lifetime (T), and the flux (F) for patches in the QS, indicate the power$-$law relationships $S\propto T~^{0.25}$ and $F\propto T~^{0.38}$, respectively. Executing a detrended fluctuation analysis of the time series of new emerged magnetic elements, we find a Hurst exponent of 0.82, which implies long-range temporal correlation in the system.

Here, we developed a complex network of solar active regions (ARs) to study various local and global properties of the network. The values of the Hurst exponent ($0.8-0.9$) were evaluated by both the detrended fluctuation analysis and the rescaled range analysis applied on the time series of the AR numbers. The findings suggest that ARs can be considered as a system of self-organized criticality. We constructed a growing network based on locations, occurrence times, and the lifetimes of 4,227 ARs recorded from 1 January 1999 to 14 April 2017. The behaviour of the clustering coefficient shows that the ARs network is not a random network. The logarithmic behaviour of the length scale has the characteristics of a so-called \textquotedblleft small-world network\textquotedblright. It is found that the probability distribution of the node degrees for undirected networks follows the power-law with exponents of about 3.7 to 4.2. This indicates the scale-free nature of the ARs network. The scale-free and small-world properties of the ARs network confirm that the system of ARs forms a system of self-organized criticality. Our results show that the occurrence probability of flares (classified by GOES class C $> 5$, M, and X flares) in the position of the ARs network hubs take values greater than that obtained for other nodes.

KIC 10553698A is a hot pulsating subdwarf B (sdB) star observed by the Kepler satellite. It exhibits dipole (l = 1) and quadrupole (l = 2) gravity modes with a clear period spacing structure. The seismic properties of the KIC 10553698A provide a test of stellar evolution models, and offer a unique opportunity to determine mixing processes. We consider mixing due to convective overshooting beyond the boundary of the helium burning core. Very small overshooting ( f = 10^{-6} ) results in a progressive increase in the size of convective core. However, moderate ( f = 10^{-2} ) and small ( f = 10^{-5} ) overshooting both lead to the occurrence of inert outer convective shells in the near-core region. We illustrate that the chemical stratifications induced by convective shells are able to change the g-mode period spacing pattern of a sdB star appreciably. The mean period spacing and trapping of the gravity modes in the model with moderate and small core overshooting are fully consistent with the period spacing trends observed in KIC 10553698A. Atomic diffusion driven by gravitational settling as well as thermal and chemical gradients is applied to reach a better match with the observed period spacings. Models that include small or very small overshooting with atomic diffusion have a decreased lifetime of the extreme horizontal branch phase and produce chemical stratification induced by convective shells during helium burning phase. In addition of being consistent with asteroseismology, their calculated values of the R2 parameter are more compatible with the observed R2 values.

Prediction of solar flares is an important task in solar physics. The occurrence of solar flares is highly dependent on the structure and the topology of solar magnetic fields. A new method for predicting large (M and X class) flares is presented, which uses machine learning methods applied to the Zernike moments of magnetograms observed by the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO) for a period of six years from 2 June 2010 to 1 August 2016. Magnetic field images consisting of the radial component of the magnetic field are converted to finite sets of Zernike moments and fed to the Support Vector Machine (SVM) classifier. Zernike moments have the capability to elicit unique features from any 2-D image, which may allow more accurate classification. The results indicate whether an arbitrary active region has the potential to produce at least one large flare. We show that the majority of large flares can be predicted within 48 hours before their occurrence, with only 10 false negatives out of 385 flaring active region magnetograms, and 21 false positives out of 179 non-flaring active region magnetograms. Our method may provide a useful tool for prediction of solar flares which can be employed alongside other forecasting methods.

In this study, the properties of wavefunctions of the MHD oscillations for a single and a system of straight flux tubes are investigated. Magnetic flux tubes with a straight magnetic field and longitudinal density stratification were considered in zero-\beta approximation. A single three-dimensional wave equation (eigen value problem) is solved for longitudinal component of the perturbed magnetic field using the finite element method (FEM). Wavefunctions (eigenfunction of wave equation) of the MHD oscillations are categorized into sausage, kink, helical kink and fluting modes. Concerning the amplitude location of the waves that arise from the flux tube, those waves are identified as body, surface, and leaky waves, appeared in the case of both a single and a system of flux tubes. Exact recognition of the wavefunctions can be used in coronal seismology and also helps to the future high resolution instruments that would be designed for studying the properties of the solar loop oscillations in details.

The solar corona is the origin of very dynamic events that are mostly produced in active regions (AR) and coronal holes (CH). The exact location of these large-scale features can be determined by applying image-processing approaches to extreme-ultraviolet (EUV) data. We here investigate the problem of segmentation of solar EUV images into ARs, CHs, and quiet-Sun (QS) images in a firm Bayesian way. On the basis of Bayes' rule, we need to obtain both prior and likelihood models. To find the prior model of an image, we used a Potts model in non-local mode. To construct the likelihood model, we combined a mixture of a Markov-Gauss model and non-local means. After estimating labels and hyperparameters with the Gibbs estimator, cellular learning automata were employed to determine the label of each pixel. We applied the proposed method to a Solar Dynamics Observatory/ Atmospheric Imaging Assembly (SDO/AIA) dataset recorded during 2011 and found that the mean value of the filling factor of ARs is 0.032 and 0.057 for CHs. The power-law exponents of the size distribution of ARs and CHs were obtained to be -1.597 and -1.508, respectively, with the maximum likelihood estimator method. When we compare the filling factors of our method with a manual selection approach and the SPoCA algorithm, they are highly compatible.

In this study, we propose methods for the automatic detection of photospheric features (bright points and granules) from ultra-violet (UV) radiation, using a feature-based classifier. The methods use quiet-Sun observations at 214 nm and 525 nm images taken by Sunrise on 9 June 2009. The function of region growing and mean shift procedure are applied to segment the bright points (BPs) and granules, respectively. Zernike moments of each region are computed. The Zernike moments of BPs, granules, and other features are distinctive enough to be separated using a support vector machine (SVM) classifier. The size distribution of BPs can be fitted with a power-law slope -1.5. The peak value of granule sizes is found to be about 0.5 arcsec^2. The mean value of the filling factor of BPs is 0.01, and for granules it is 0.51. There is a critical scale for granules so that small granules with sizes smaller than 2.5 arcsec^2 cover a wide range of brightness, while the brightness of large granules approaches unity. The mean value of BP brightness fluctuations is estimated to be 1.2, while for granules it is 0.22. Mean values of the horizontal velocities of an individual BP and an individual BP within the network were found to be 1.6 km/s and 0.9 km/s, respectively. We conclude that the effect of individual BPs in releasing energy to the photosphere and maybe the upper layers is stronger than what the individual BPs release into the network.

Solar mini-dimmings can be detect in the Extreme Ultra-Violet coro- nal eruptions. Here, sequences of 171 A images taken by Solar Dynamic Observa- tory/Atmospheric Imaging Assembaly on 13 June 2010 are used. In this special day, both of coronal mass ejection and prominence were observed. The average velocities and accelerations of 500 mini-dimmings which were detected using on feature based classifier (Alipour et al 2012) are studied. The large number of mini-dimmings shows positive accelerations in the beginning times as similar as Coronal Mass Ejections. On the the start time of prominence eruptions, the number of mini-dimmings is increased to a maximum values. There is a positive correlation between the kinematics of mini- dimmings and both CME and prominences. This study can be extended to understand the exact relationship of CMEs and mini-dimmings.

The effect of temperature inhomogeneity on the periods, their ratios (fundamental vs. first overtone), and the damping times of the standing slow modes in gravitationally stratified solar coronal loops are studied. The effects of optically thin radiation, compressive viscosity, and thermal conduction are considered. The linearized one-dimensional magnetohydrodynamic (MHD) equations (under low-$\beta$ condition) were reduced to a fourth--order ordinary differential equation for the perturbed velocity. The numerical results indicate that the periods of non-isothermal loops (i.e. temperature increases from the loop base to apex) are smaller compared to those of isothermal loops. In the presence of radiation, viscosity, and thermal conduction, an increase in the temperature gradient is followed by a monotonic decrease in the periods (compared with the isothermal case), while the period ratio turns out to be a sensitive function of the gradient of the temperature and the loop lengths. We verify that radiative dissipation is not a main cooling mechanism of both isothermal and non-isothermal hot coronal loops and has a small effect on the periods. Thermal conduction and compressive viscosity are primary mechanisms in the damping of slow modes of the hot coronal loops. The periods and damping times in the presence of compressive viscosity and/or thermal conduction dissipation are consistent with the observed data in specific cases. By tuning the dissipation parameters, the periods and the damping times could be made consistent with the observations in more general cases.

Nanoflares, the basic unit of impulsive energy release may produce much of the solar background emission. Extrapolation of the energy frequency distribution of observed microflares, which follows a power law to lower energies can give an estimation of the importance of nanoflares for heating the solar corona. If the power law index is greater than 2, then the nanoflare contribution is dominant. We model time series of extreme ultraviolet emission radiance, as random flares with a power law exponent of the flare event distribution. The model is based on three key parameters, the flare rate, the flare duration and the power law exponent of the flare intensity frequency distribution. We use this model to simulate emission line radiance detected in 171 \AA, observed by STEREO/EUVI and SDO/AIA. The Observed light curves are matched with simulated light curves using an Artificial Neural Network and parameter values are determined across regions of active region, quiet sun, and coronal hole. The damping rate of nanoflares is compared with radiative losses cooling time. The effect of background emission, data cadence, and network sensitivity on the key parameters of model is studied. Most of the observed light curves have a power law exponent, $\alpha$, greater than the critical value 2. At these sites nanoflare heating could be significant.

Small-scale extreme ultraviolet (EUV) dimming often surrounds sites of energy release in the quiet Sun. This paper describes a method for the automatic detection of these small-scale EUV dimmings using a feature based classifier. The method is demonstrated using sequences of 171 A images taken by STEREO/EUVI on 13 June 2007 and by SDO/AIA on 27 August 2010. The feature identification relies on recognizing structure in sequences of space-time 171\AA\ images using the Zernike moments of the images. The Zernike moments space-time slices with events and non-events are distinctive enough to be separated using a Support Vector Machine (SVM) classifier. The SVM is trained using 150 event and 700 non-event space-time slices. We find a total of 1217 events in the EUVI images and 2064 events in the AIA images on the days studied. Most of the events are found between latitudes -35 degree and +35 degree. The sizes and expansion speeds of central dimming regions are extracted using a region grow algorithm. The histograms of the sizes in both EUVI and AIA follow a steep power law with slope about -5. The AIA slope extends to smaller sizes before turning over. The mean velocity of 1325 dimming regions seen by AIA is found to be about 14 km/s.

We study the resonant absorption of MHD waves in magnetized flux tubes with a radial density inhomogeneity. Within the approximation that resistive and viscous processes are operative in thin layers surrounding the singularities of the MHD equations, we give the full spectrum of the eigenfrequencies and damping rates of the MHD quasi modes of the tube. Both surface and body modes are analyzed. \keywords{Sun -- corona: magnetohydrodynamics (MHD) -- Sun: magnetic fields -- Sun: oscillations

Context. Nanoflares are small impulsive bursts of energy that blend with and possibly make up much of the solar background emission. Determining their frequency and energy input is central to understanding the heating of the solar corona. One method is to extrapolate the energy frequency distribution of larger individually observed flares to lower energies. Only if the power law exponent is greater than 2, is it considered possible that nanoflares contribute significantly to the energy input. Aims. Time sequences of ultraviolet line radiances observed in the corona of an active region are modelled with the aim of determining the power law exponent of the nanoflare energy distribution. Methods. A simple nanoflare model based on three key parameters (the flare rate, the flare duration time, and the power law exponent of the flare energy frequency distribution) is used to simulate emission line radiances from the ions Fe XIX, Ca XIII, and Si iii, observed by SUMER in the corona of an active region as it rotates around the east limb of the Sun. Light curve pattern recognition by an Artificial Neural Network (ANN) scheme is used to determine the values. Results. The power law exponents, alpha 2.8, 2.8, and 2.6 for Fe XIX, Ca XIII, and Si iii respectively. Conclusions. The light curve simulations imply a power law exponent greater than the critical value of 2 for all ion species. This implies that if the energy of flare-like events is extrapolated to low energies, nanoflares could provide a significant contribution to the heating of active region coronae.

The influence of longitudinal structuring on the fast kink modes of coronal loops is investigated. Analytical dispersion relations and mode profiles are derived for the second-order ordinary differential equation governing the z- component of the perturbation in the magnetic field, $\delta B_z$. All other components are given in terms of $\delta B_z$. Deviations from the frequencies and mode profiles of homogenous loops are given as functions of the density scale height. The effects of the fixed and variable column masses, negative scale heights, and density contrasts inside and outside of the loops are studied. The frequency ratios, mode profiles, and the antinode shifts from those of the sine profiles of the homogenous loops, are tools to estimate solar photospheric parameters. To this end, we have expanded the relevant factors up to the second order in the stratification parameter. In particular, we verify that the first overtone antinode shifts are in the Mm range and are within the reach of the resolutions of the present day observations.

We investigate the standing kink modes of a cylindrical model of coronal loops. The density is stratified along the loop axis and changes discontinuously at the surface of the cylinder. The periods and mode profiles are studies with their deviation from those of the unstratified loops. The aim is to extract information on the density scale heights prevailing in the solar corona. The problem is reduced to solving a single second-order partial differential equation for $\delta B_z(r,z)$, the longitudinal component of the Eulerian perturbations of the magnetic field. This equation, in turn, is separated into two second-order ordinary, differential equations in $r$ and $z$ that are, however, connected through a dispersion relation between the frequencies and the longitudinal wave numbers. In the thin tube approximation, the eigensolutions are obtained by a perturbation technique, where the perturbation parameter is the density stratification parameter. Otherwise the problem is solved numerically.