Science Verification (SV) data of solar observations with Atacama Large Millimeter-submillimeter Array (ALMA) telescope were released to the scientific community. Understanding the centre to limb brightness function is necessary to compare features in full disk images. Our goals are to find the empirical centre to limb brightness functions in two available spectral bands and create flattened images with centre to limb brightness variations removed. We used second-order polynomial fit of the cosine of incidence angle to data points as a function of radial distance to the centre of the solar disk. The method also includes iterative removal of outliers based on the interquartile range. Fitting functions for all available images proved to adequately describe the data with comparatively small errors in the fitting coefficients. In both bands we found brightening towards the limb which is a consequence of increase in electron temperatures with radial distance in this region of the solar atmosphere. This study found that the $T_{b}$ of an active region has about 180 K difference between with and without the limb brightening at radial distance $\approx$0.75$R_{\odot}$ in Band 6. We also made flattened images with limb brightening removed. The limb brightening effect in ALMA images is significant enough (of the order of 10% for Band 3 and about 15% in Band 6) that it can not be neglected in further analyses. Since the effect of the side lobes was not included in this study, these values probably represent the lower limit of the limb brightening. The shape of the limb brightening function can also be used to constrain electron densities and temperatures in various layers of the solar atmosphere.

Reliable determination of the basic physical properties of hot emission-line binaries with Roche-lobe filling secondaries is important for developing the theory of mass exchange in binaries. It is not easy, however, due to the presence of circumstellar matter. Here, we report the first detailed investigation of a new representative of this class of binaries, HD~81357, based on the analysis of spectra and photometry from several observatories. HD~81357 was found to be a double-lined spectroscopic binary and an ellipsoidal variable seen under an intermediate orbital inclination of $\sim(63\pm5)^\circ$, having an orbital period of 33\fd77445(41) and a~circular orbit. From an automated comparison of the observed and synthetic spectra, we estimate the component's effective temperatures to be 12930(540)~K and 4260(24)~K. The combined light-curve and orbital solutions, also constrained by a very accurate Gaia Data Release 2 parallax, give the following values of the basic physical properties: masses $3.36\pm0.15$ and $0.34\pm0.04$~\Mnom, radii $3.9\pm0.2$ and 13.97\pm0.05$~\Rnom, and a~mass ratio $10.0\pm0.5$. Evolutionary modelling of the system including the phase of mass transfer between the components indicated that HD~81357 is a~system observed in the final slow phase of the mass exchange after the mass-ratio reversal. Contrary to what has been seen for similar binaries like AU~Mon, no cyclic light variations were found on a~time scale an~order of magnitude longer than the orbital period. 243,1 15%

In this work we use solar observations with the ALMA radio telescope at the wavelength of 1.21 mm. The aim of the analysis is to improve understanding of the solar chromosphere, a dynamic layer in the solar atmosphere between the photosphere and corona. The study has an observational and a modeling part. In the observational part full-disc solar images are analyzed. Based on a modified FAL atmospheric model, radiation models for various observed solar structures are developed. Finally, the observational and modeling results are compared and discussed.

Sunspot position data obtained from Kanzelh\"{o}he Observatory for Solar and Environmental Research (KSO) sunspot drawings and white light images in the period 1964 to 2016 were used to calculate the rotational and meridional velocities of the solar plasma. Velocities were calculated from daily shifts of sunspot groups and an iterative process of calculation of the differential rotation profiles was used to discard outliers. We found a differential rotation profile and meridional motions in agreement with previous studies using sunspots as tracers and conclude that the quality of the KSO data is appropriate for analysis of solar velocity patterns. By analysing the correlation and covariance of meridional velocities and rotation rate residuals we found that the angular momentum is transported towards the solar equator. The magnitude and latitudinal dependence of the horizontal component of the Reynolds stress tensor calculated is sufficient to maintain the observed solar differential rotation profile. Therefore, our results confirm that the Reynolds stress is the dominant mechanism responsible for transport of angular momentum towards the solar equator.

Various solar features can be seen on maps of the Sun in the mm and sub-mm wavelength range. The recently installed Atacama Large Millimeter/submillimeter Array (ALMA) is capable of observing the Sun in that wavelength range with an unprecedented spatial, temporal and spectral resolution. To interpret solar observations with ALMA the first important step is to compare ALMA maps with simultaneous images of the Sun recorded in other spectral ranges. First we identify different structures in the solar atmosphere seen in the optical, IR and EUV parts of the spectrum (quiet Sun (QS), active regions (AR), prominences on the disc, magnetic inversion lines (IL), coronal holes (CH) and coronal bright points (CBPs)) in a full disc solar ALMA image. The second aim is to measure the intensities (brightness temperatures) of those structures and compare them with the corresponding QS level. A full disc solar image at 1.21 mm obtained on December 18, 2015 during a CSV-EOC campaign with ALMA is calibrated and compared with full disc solar images from the same day in H\alpha, in He I 1083 nm core, and with SDO images (AIA at 170 nm, 30.4 nm, 21.1 nm, 19.3 nm, and 17.1 nm and HMI magnetogram). The brightness temperatures of various structures are determined by averaging over corresponding regions of interest in the ALMA image. Positions of the QS, ARs, prominences on the disc, ILs, CHs and CBPs are identified in the ALMA image. At 1.21 mm ARs appear as bright areas (but sunspots are dark), while prominences on the disc and CHs are not discernible from the QS background, although having slightly less intensity than surrounding QS regions. ILs appear as large, elongated dark structures and CBPs correspond to ALMA bright points. These results are in general agreement with sparse earlier measurements at similar wavelengths. The identification of CBPs represents the most important new result.

The sunspot position from Greenwich Photoheliographic Results (GPR), US Air Force Solar Optical Observing Network and National Oceanic and Atmospheric Administration (USAF/NOAA), and Debrecen Photoheliographic Data (DPD) data bases in the period 1874 to 2016 were used to calculate yearly values of the solar differential-rotation parameters $A$ and $B$. The calculated differential-rotation parameters were compared with the solar-activity level. We found that the Sun rotates more differentially at the minimum than at the maximum of activity during the 1977 - 2016 epoch. An inverse correlation between equatorial rotation and solar activity was found using the recently revised sunspot number. The secular decrease of equatorial rotation rate accompanying the increase of activity stopped in the last part of the 20th century. It was noted that when a significant peak of equatorial rotation velocity is observed during minimum of activity, the strength of the next maximum is smaller than the previous one.

The main aim of this work is to determine the solar differential rotation by tracing sunspot groups during the period 1964-2016, using the Kanzelh\"ohe Observatory for Solar and Environmental Research (KSO) sunspot drawings and white light images. Two procedures for the determination of the heliographic positions were applied: an interactive procedure on the KSO sunspot drawings (1964 - 2008, solar cycles nos. 20 - 23) and an automatic procedure on the KSO white light images (2009 - 2016, solar cycle no. 24). For the determination of the synodic angular rotation velocities two different methods have been used: a daily shift (DS) method and a robust linear least-squares fit (rLSQ) method. Afterwards, the rotation velocities had to be converted from synodic to sidereal, which were then used in the least-squares fitting for the solar differential rotation law. For the test data from 2014, we found the rLSQ method for calculating rotational velocities to be more reliable than the DS method. The best fit solar differential rotation profile for the whole time period is $\omega(b)$ = (14.47 $\pm$ 0.01) - (2.66 $\pm$ 0.10) $\sin^2b$ (deg/day) for the DS method and $\omega(b)$ = (14.50 $\pm$ 0.01) - (2.87 $\pm$ 0.12) $\sin^2b$ (deg/day) for the rLSQ method. A barely noticeable north - south asymmetry is observed for the whole time period 1964 - 2016 in the present paper. Rotation profiles, using different data sets (e.g. Debrecen Photoheliographic Data, Greenwich Photoheliographic Results), presented by other authors for the same time periods and the same tracer types, are in good agreement with our results. Therefore, the KSO data set is suitable for the investigation of the long-term variabilities in the solar rotation profile.

The Debrecen Photoheliographic Data catalogue is a continuation of the Greenwich Photoheliographic Results providing daily positions of sunspots and sunspot groups. We analyse the data for sunspot groups focusing on meridional motions and transfer of angular momentum towards the solar equator. Velocities are calculated with a daily shift method including an automatic iterative process of removing the outliers. Apart from the standard differential rotation profile, we find meridional motion directed towards the zone of solar activity. The difference in measured meridional flow in comparison to Doppler measurements and some other tracer measurements is interpreted as a consequence of different flow patterns inside and outside of active regions. We also find a statistically significant dependence of meridional motion on rotation velocity residuals confirming the transfer of angular momentum towards the equator. Analysis of horizontal Reynolds stress reveals that the transfer of angular momentum is stronger with increasing latitude up to about 40$^{\circ}$ where there is a possible maximum in absolute value.

Context. It is possible to detect and track coronal bright points (CBPs) in SDO/AIA images. Combination of high resolution and high cadence provides a wealth of data that can be used to determine velocity flows on the solar surface with very high accuracy. Aims. We derived a very accurate solar rotation profile and investigated meridional flows, torsional oscillations and horizontal Reynolds stress based on $\approx$6 months of SDO/AIA data. Methods. We used a segmentation algorithm to detect CBPs in SDO/AIA images. We also used invariance of the solar rotation profile with central meridian distance (CMD) to determine the height of CBPs in 19.3 nm channel. Results. Best fit solar rotation profile is given by $\omega(b)=(14.4060\pm0.0051 + (-1.662\pm0.050)\sin^{2}b + (-2.742\pm0.081)\sin^{4}b)${\degr} day$^{-1}$. Height of CBPs in SDO/AIA 19.3 nm channel was found to be $\approx$6500 km. Meridional motion is predominantly poleward for all latitudes, while solar velocity residuals show signs of torsional oscillations. Horizontal Reynolds stress was found to be small compared to similar works, but still showing transfer of angular momentum towards the solar equator. Conclusion. Most of the results are consistent with Doppler measurements rather than tracer measurements. Fairly small calculated value of horizontal Reynolds stress might be due to the particular phase of the solar cycle. Accuracy of the calculated rotation profile indicates that it is possible to measure changes in the profile as the solar cycle evolves. Analysis of further SDO/AIA CBP data will also help in better understanding of the temporal behaviour of the rotation velocity residuals, meridional motions and Reynolds stress.

Predicting transit times of Coronal Mass Ejections (CMEs) from their initial parameters is a very important subject, not only from the scientific perspective, but also because CMEs represent a hazard for human technology. We used a neural network to analyse transit times for 153 events with only two input parameters: initial velocity of the CME, $v$, and Central Meridian Distance, CMD, of its associated flare. We found that transit time dependence on $v$ is showing a typical drag-like pattern in the solar wind. The results show that the speed at which acceleration by drag changes to deceleration is $v\approx$500 km s$^{-1}$. Transit times are also found to be shorter for CMEs associated with flares on the western hemisphere than those originating on the eastern side of the Sun. We attribute this difference to the eastward deflection of CMEs on their path to 1 AU. The average error of the NN prediction in comparison to observations is $\approx$12 hours which is comparable to other studies on the same subject.

Coronal bright points (CBP) are ubiquitous small brightenings in the solar corona associated with small magnetic bipoles. We derive the solar differential rotation profile by tracing the motions of CBPs detected by the Atmospheric Imaging Assembly (AIA) instrument aboard the Solar Dynamics Observatory (SDO). We also investigate problems related to detection of coronal bright points resulting from instrument and detection algorithm limitations. To determine the positions and identification of coronal bright points we used a segmentation algorithm. A linear fit of their central meridian distance and latitude versus time was utilised to derive velocities. We obtained 906 velocity measurements in a time interval of only 2 days. The differential rotation profile can be expressed as $\omega_{rot} = (14.47\pm 0.10 + (0.6\pm 1.0)\sin^{2}(b) + (-4.7\pm 1.7)\sin^{4}(b))$\degr day$^{-1}$. Our result is in agreement with other work and it comes with reasonable errors in spite of the very short time interval used. This was made possible by the higher sensitivity and resolution of the AIA instrument compared to similar equipment as well as high cadence. The segmentation algorithm also played a crucial role by detecting so many CBPs, which reduced the errors to a reasonable level. Data and methods presented in this paper show a great potential to obtain very accurate velocity profiles, both for rotation and meridional motion and, consequently, Reynolds stresses. The amount of coronal bright point data that could be obtained from this instrument should also provide a great opportunity to study changes of velocity patterns with a temporal resolution of only a few months. Other possibilities are studies of evolution of CBPs and proper motions of magnetic elements on the Sun.

The main objective of the study is to determine the probability distributions of the geomagnetic Dst index as a function of the coronal mass ejection (CME) and solar flare parameters for the purpose of establishing a probabilistic forecast tool for the geomagnetic storm intensity. Several CME and flare parameters as well as the effect of successive-CME occurrence in changing the probability for a certain range of Dst index values, were examined. The results confirm some of already known relationships between remotely-observed properties of solar eruptive events and geomagnetic storms, namely the importance of initial CME speed, apparent width, source position, and the associated solar flare class. In this paper we quantify these relationships in a form to be used for space weather forecasting in future. The results of the statistical study are employed to construct an empirical statistical model for predicting the probability of the geomagnetic storm intensity based on remote solar observations of CMEs and flares.

Thanks to its long and rich observational history and rapid apsidal motion, the massive eclipsing binary Y Cyg represents one of the cornestones to critical tests of stellar evolution theory for massive stars. Yet, the determination of the basic physical properties is less accurate than it could be given the existing number of spectral and photometric observations. Our goal is to analyze all these data simultaneously with the new dedicated series of our own spectral and photometric observations from observatories widely separated in longitude. We obtained new series of UBV observations at three observatories separated in local time to obtain complete light curves of Y Cyg for its orbital period close to 3 days. This new photometry was reduced and carefully transformed to the standard UBV system using the HEC22 program. We also obtained new series of red spectra secured at two observatories and re-analyzed earlier obtained blue electronic spectra. Our analyses provide the most accurate so far published value of the apsidal period of 47.805 +/- 0.030 yrs and the following physical elements: M1=17.72+/-0.35$ Msun, M2=17.73+/-0.30 Msun, R1=5.785+/-0.091 Rsun, and R2=5.816+/-0.063 Rsun. The disentangling thus resulted in the masses, which are somewhat higher than all previous determinations and virtually the same for both stars, while the light curve implies a slighly higher radius and luminosity for star 2. The above empirical values imply the logarithm of the internal structure constant log k2 = -1.937. A comparison with Claret's stellar interior models implies an age close to 2 millions yrs for both stars. The claimed accuracy of modern element determination of 1-2 per cent seems still a bit too optimistic and obtaining new high-dispersion and high-resolution spectra is desirable.

The goal of this paper is to investigate Reynolds stresses and to check if it is plausible that they are responsible for angular momentum transfer toward the solar equator. We also analysed meridional velocity, rotation velocity residuals and correlation between the velocities. We used sunspot groups position measurements from GPR (Greenwich Photographic Result) and SOON/USAF/NOAA (Solar Observing Optical Network/United States Air Force/National Oceanic and Atmospheric Administration) databases covering the period from 1878 until 2011. In order to calculate velocities we used daily motion of sunspot groups. The sample was also limited to $\pm$58\degr in Central Meridian Distance (CMD) in order to avoid solar limb effects. We mainly investigated velocity patterns depending on solar cycle phase and latitude. We found that meridional motion of sunspot groups is toward the centre of activity from all available latitudes and in all phases of the solar cycle. The range of meridional velocities is $\pm10$ m s$^{-1}$. Horizontal Reynolds stress is negative at all available latitudes and indicates that there is a minimum value ($q\approx$ - 3000 m$^2$ s$^{-2}$) located at $b\approx\pm$30$^{0}$. In our convention this means that angular momentum is transported toward the solar equator in agreement with the observed rotational profile of the Sun.

We describe the construction of a highly reliable sample of approximately 7,000 optically faint periodic variable stars with light curves obtained by the asteroid survey LINEAR across 10,000 sq.deg of northern sky. Majority of these variables have not been cataloged yet. The sample flux limit is several magnitudes fainter than for most other wide-angle surveys; the photometric errors range from ~0.03 mag at $r=15$ to ~0.20 mag at r=18. Light curves include on average 250 data points, collected over about a decade. Using SDSS-based photometric recalibration of the LINEAR data for about 25 million objects, we selected ~200,000 most probable candidate variables and visually confirmed and classified approximately 7,000 periodic variables using phased light curves. The reliability and uniformity of visual classification across eight human classifiers was calibrated and tested using a SDSS Stripe 82 region variable star catalog, and verified using an unsupervised machine learning approach. The resulting sample of periodic LINEAR variables is dominated by 3,900 RR Lyrae stars and 2,700 eclipsing binary stars of all subtypes, and includes small fractions of relatively rare populations such as asymptotic giant branch stars and SX Phoenicis stars. We discuss the distribution of these mostly uncataloged variables in various diagrams constructed with optical-to-infrared SDSS, 2MASS and WISE photometry, and with LINEAR light curve features. An interesting side result is a robust and precise quantitative description of a strong correlation between the light-curve period and color/spectral type for close and contact eclipsing binary stars. These large samples of robustly classified variable stars will enable detailed statistical studies of the Galactic structure and physics of binary and other stars, and we make them publicly available.

Using our new photometric and spectroscopic observations as well as all available published data, we present a new interpretation of the properties of the peculiar emission-line binary UX Mon. We conclude that this binary is in a rare phase of fast mass transfer between the binary components prior to the mass ratio reversal. We firmly establish that the orbital period is secularly decreasing at a rate of $\dot{P}=-0.260$ seconds per year. From several lines of reasoning, we show that the mass ratio of the component losing mass to the mass-gaining component $q$ must be larger than 1 and find our most probable value to be $q=1.15\pm0.1$. The BINSYN suite of programs and the steepest descent method were used to perform the final modeling. We modeled the star as a W Ser star with a thick disk around its primary. Although the remaining uncertainties in some of the basic physical elements describing the system in our model are not negligible, the model is in fair agreement with available observations. Only the nature of the light variations outside the primary eclipse remains unexplained.

We present a new observational campaign, DWARF, aimed at detection of circumbinary extrasolar planets using the timing of the minima of low-mass eclipsing binaries. The observations will be performed within an extensive network of relatively small to medium-size telescopes with apertures of ~20-200 cm. The starting sample of the objects to be monitored contains (i) low-mass eclipsing binaries with M and K components, (ii) short-period binaries with sdB or sdO component, and (iii) post-common-envelope systems containing a WD, which enable to determine minima with high precision. Since the amplitude of the timing signal increases with the orbital period of an invisible third component, the timescale of project is long, at least 5-10 years. The paper gives simple formulas to estimate suitability of individual eclipsing binaries for the circumbinary planet detection. Intrinsic variability of the binaries (photospheric spots, flares, pulsation etc.) limiting the accuracy of the minima timing is also discussed. The manuscript also describes the best observing strategy and methods to detect cyclic timing variability in the minima times indicating presence of circumbinary planets. First test observation of the selected targets are presented.

The A-type star HR 6412 = V2368 Oph was used by several investigators as a photometric comparison star for the known eclipsing binary U Oph but was found to be variable by three independent groups, including us. By analysing series of new spectral and photometric observations and a critical compilation of available radial velocities, we were able to find the correct period of light and radial-velocity variations and demonstrate that the object is an eclipsing and double-lined spectroscopic binary moving in a highly eccentric orbit. We derived a linear ephemeris T min.I = HJD (2454294.67 +/- 0.01) + (38.32712 +/- 0.00004)d x E and estimated preliminary basic physical properties of the binary. The dereddened UBV magnitudes and effective temperatures of the primary and secondary, based on our light- and velocity-curve solutions, led to distance estimates that agree with the Hipparcos distance within the errors. We find that the mass ratio must be close to one, but the limited number and wavelength range of our current spectra does not allow a truly precise determination of the binary masses. Nevertheless, our results show convincingly that both binary components are evolved away from the main sequence, which makes this system astrophysically very important. There are only a few similarly evolved A-type stars among known eclipsing binaries. Future systematic observations and careful analyses can provide very stringent tests for the stellar evolutionary theory.

A series of 353 red electronic spectra obtained between 1994 and 2010, and of 171 UBV photometric observations of the 2010 eclipse, were analyzed in an effort to better understand the eclipsing binary eps Aur. The main results follow. (1) We attempted to recover a spectrum of the companion by disentangling the observed spectra of the eps Aur binary failed, but we were able to disentangle the spectrum of telluric lines and obtain a mean spectrum of the F-type primary star. The latter was then compared to a grid of synthetic spectra for a number of plausible values of T(eff) and log(g), but a reasonably good match was not found. However, we conclude that the observed spectrum is that of a low-gravity star. (2) We examined changes in the complex H-alpha line profiles over the past 16 years, with particular emphasis on the 2009-2011 eclipse period, by subtracting a mean out-of-eclipse H-alpha profile (appropriately shifted in radial velocity) from the observed spectra. We find that the dark disk around the unseen companion has an extended 'atmosphere' that manifests itself via blueshifted and redshifted H-alpha 'shell' absorptions seen projected against the F star. Significantly, the H-alpha shell line first appeared three years before first contact of the optical eclipse when the system was not far past maximum separation. (3) Analyses of radial velocities and central intensities of several strong, unblended spectral lines, as well as UBV photometry, demonstrated that these observables showed apparent multiperiodic variability during eclipse. The dominant period of 66.21 was common to all the observables, but with different phase shifts between these variables. This result strongly supports our earlier suggestion that the photometric variability seen during eclipse is intrinsic to the F star, and therefore, the idea of a central brightening due to a hole in the disk should be abandoned.

We collected rich series of RV measurements covering last 110 years and photometric observations from the past 6 primary eclipses, complemented them by our new observations and derived a new precise ephemeris and an orbital solution of epsilon Aur.