We report the Transiting Exoplanet Survey Satellite ($TESS$) detection of a multi-planet system orbiting the $V=10.9$ K0 dwarf TOI 125. We find evidence for up to five planets, with varying confidence. Three high signal-to-noise transit signals correspond to sub-Neptune-sized planets ($2.76$, $2.79$, and $2.94\ R_{\oplus}$), and we statistically validate the planetary nature of the two inner planets ($P_b = 4.65$ days, $P_c = 9.15$ days). With only two transits observed, we report the outer object ($P_{.03} = 19.98$ days) as a high signal-to-noise ratio planet candidate. We also detect a candidate transiting super-Earth ($1.4\ R_{\oplus}$) with an orbital period of only $12.7$ hours and a candidate Neptune-sized planet ($4.2\ R_{\oplus}$) with a period of $13.28$ days, both at low signal-to-noise. This system is amenable to mass determination via radial velocities and transit timing variations, and provides an opportunity to study planets of similar size while controlling for age and environment. The ratio of orbital periods between TOI 125 b and c ($P_c/P_b = 1.97$) is slightly smaller than an exact 2:1 commensurability and is atypical of multiple planet systems from $Kepler$, which show a preference for period ratios just $wide$ of first-order period ratios. A dynamical analysis refines the allowed parameter space through stability arguments and suggests that, despite the nearly commensurate periods, the system is unlikely to be in resonance.

We present results of the combined photometric and spectroscopic analysis of three detached eclipsing binaries, which secondary components are not visible or very hard to identify in the optical spectra - ASAS J052743-0359.7, ASAS J065134-2211.5, and ASAS J073507-0905.7. The first one is a known visual binary ADS 4022, and we found that it is a quadruple system, composed of two spectroscopic binaries, one of which shows eclipses. None of the systems was previously recognized as a spectroscopic binary. We collected a number of high-resolution optical and IR spectra to calculate the radial velocities (RVs) and later combined them with MITSuME and ASAS photometry. The IR spectra were crucial for secure identification of the cooler components' lines. RV measurements were done with the TODCOR technique, and RV curves modelled with our own procedure V2FIT. Light curve modelling was performed with JKTEBOP and PHOEBE codes. Temperatures and metallicities of two systems were estimated from spectra. For the ADS 4022 system we also used the archival WDS data and new SOAR observations in order to derive the orbit of the visual pair for the first time. Ages were estimated by comparing our results with PARSEC isochrones. The eclipsing pair A052743 A ($P=5.27$ d) is composed of a 1.03(6) M$_\odot$, 1.03(2) R$_\odot$ primary and a 0.60(2) M$_\odot$, 0.59(2) R$_\odot$ secondary. The components of the $P=21.57$ d non-eclipsing pair B likely have masses in between the two eclipsing components, and both pairs are on a $\sim$188 yr orbit around their common centre of mass. The system A065134 ($P=8.22$ d) consists of a 0.956(12) M$_\odot$, 0.997(4) R$_\odot$ primary and a 0.674(5) M$_\odot$, 0.690(7) R$_\odot$ secondary. Finally, A073507 ($P=1.45$ d), which consists of a 1.452(34) M$_\odot$, 1.635(12) R$_\odot$ primary and a 0.808(13) M$_\odot$, 0.819(11) R$_\odot$ secondary, is likely a PMS system.

We have obtained high-resolution spectra of four eclipsing binary systems (FM Leo, NN Del, V963 Cen and AI Phe) with the view to gaining insight into the relative orientation of their stellar spin axes and orbital axes. The so called Rossiter-McLaughlin (RM) effect, i.e. the fact that the broadening and the amount of blue- or redshift in the spectra during an eclipse depend on the tilt of the spin axis of the background star, has the potential of reconciling observations and theoretical models if such a tilt is found. We analyse the RM effect by disentangling the spectra, removing the front component and measuring the remaining, distorted lines with a broadening function (BF) obtained from single value decomposition (SVD), weighting by the intensity centre of the BF in the eclipse. All but one of our objects show no significant misalignment, suggesting that aligned systems are dominant. We provide stellar as well as orbital parameters for our systems. With five measured spin-orbit angles we significantly increase (from 9 to 14) the number of stars for which it has been measured. The spin-orbit angle $\beta$ calculated for AI Phe's secondary component shows a misalignment of 87$\pm$17 degrees. NN Del, with a large separation of components and a long dynamical timescale for circularisation and synchronisation, is an example of a close to primordial spin-orbit angle measurement.

We present the design and commissioning of Project Solaris, a global network of autonomous observatories. Solaris is a Polish scientific undertaking aimed at the detection and characterization of circumbinary exoplanets and eclipsing binary stars. To accomplish this, a network of four fully autonomous observatories has been deployed in the Southern Hemisphere: Solaris-1 and Solaris-2 in the South African Astronomical Observatory in South Africa; Solaris-3 in Siding Spring Observatory in Australia; and Solaris-4 in Complejo Astronomico El Leoncito in Argentina. The four stations are nearly identical and are equipped with 0.5-m Ritchey-Cr\'etien (f/15) or Cassegrain (f/9, Solaris-3) optics and high-grade 2 K x 2 K CCD cameras with Johnson and Sloan filter sets. We present the design and implementation of low-level security; data logging and notification systems; weather monitoring components; all-sky vision system, surveillance system; and distributed temperature and humidity sensors. We describe dedicated grounding and lighting protection system design and robust fiber data transfer interfaces in electrically demanding conditions. We discuss the outcomes of our design, as well as the resulting software engineering requirements. We describe our system's engineering approach to achieve the required level of autonomy, the architecture of the custom high-level industry-grade software that has been designed and implemented specifically for the use of the network. We present the actual status of the project and first photometric results; these include data and models of already studied systems for benchmarking purposes (Wasp-4b, Wasp-64b, and Wasp-98b transits, PG 1663-018, an eclipsing binary with a pulsator) as well J024946-3825.6, an interesting low-mass binary system for which a complete model is provided for the first time.

We present the results of our analysis of KIC 4150611 (HD 181469) - an interesting, bright quintuple system that includes a hybrid $\delta$ Sct/$\gamma$ Dor pulsator. Four periods of eclipses - 94.2, 8.65, 1.52 and 1.43 d - have been observed by the Kepler satellite, and three point sources (A, B, and C) are seen in high angular resolution images. From spectroscopic observations made with the HIDES spectrograph attached to the 1.88-m telescope of the Okayama Astrophysical Observatory (OAO), for the first time we calculated radial velocities (RVs) of the component B - a pair of G-type stars - and combined them with Kepler photometry in order to obtain absolute physical parameters of this pair. We also managed to directly measure RVs of the pulsator, also for the first time. Additionally, we modelled the light curves of the 1.52 and 1.43-day pairs, and measured their eclipse timing variations (ETVs). We also performed relative astrometry and photometry of three sources seen on the images taken with the NIRC2 camera of the Keck II telescope. Finally, we compared our results with theoretical isochrones. The brightest component Aa is the hybrid pulsator, transited every 94.2 days by a pair of K/M-type stars (Ab1+Ab2), which themselves form a 1.52-day eclipsing binary. The components Ba and Bb are late G-type stars, forming another eclipsing pair with a 8.65 day period. Their masses and radii are $M_{Ba}=0.894\pm0.010$ M$_\odot$, $R_{Ba}=0.802\pm0.044$ R$_\odot$ for the primary, and $M_{Bb}=0.888\pm0.010$ M$_\odot$, $R_{Bb}=0.856\pm0.038$ R$_\odot$ for the secondary. The remaining period of 1.43 days is possibly related to a faint third star C, which itself is most likely a background object. The system's properties are well-represented by a 35 Myr isochrone. There are also hints of additional bodies in the system.

We present the results of our spectroscopic observations of eight detached eclipsing binaries (DEBs), selected from the $Kepler$ Eclipsing Binary Catalog. Radial velocities (RVs) were calculated from high resolution spectra obtained with the HIDES pectrograph, attached to the 1.88-m telescope of the Okayama Astrophysical Observatory, and were used to characterize the targets in combination with the $Kepler$ light curves. For each binary we obtained a full set of orbital and physical parameters, reaching precision below 3% in masses and radii for 5 pairs. By comparing our results with theoretical models, we assess the distance, age and evolutionary status of the researched objects. We also study eclipse timing variations of selected objects, and identify a new system with a $\gamma$Dor pulsator. Two systems are triples, and show lines coming from three components. In one case the motion of the outer star and the perturbation in the RVs of the inner binary are clearly visible and periodical, which allows us to directly calculate the mass of the third star, and inclination of the outer orbit. In the second case we only see a clear motion of the tertiary, and investigate two scenarios: that it is a linear trend coming from the orbital motion around the inner binary, and that it is caused by a planetary mass companion. When possible, we also compare our results with the literature, and conclude that only by combining photometry with RVs it is possible to obtain correct physical parameters of both components of a DEB.

We present scientific results obtained with a recently commissioned \'{e}chelle spectrograph on the 0.5-m Solaris-1 telescope in the South African Astronomical Observatory. BACHES is a low-cost slit \'{e}chelle spectrograph that has a resolution of 21,000 at 5,500 \AA. The described setup is fully remotely operated and partly automated. Custom hardware components have been designed to allow both spectroscopic and photometric observations. The setup is controlled via dedicated software. The throughput of the system allows us to obtain spectra with an average SNR of 22 at 6375 {\AA} for a 30-min exposure of a $V=10$ mag target. The stability of the instrument is influenced mainly by the ambient temperature changes. We have obtained radial velocity RMS values for a bright (V = 5.9 mag) spectroscopic binary as good as 0.59 km s$^{-1}$ and 1.34 km s$^{-1}$ for a $V = 10.2$ mag eclipsing binary. Radial velocity measurements have been combined with available photometric light curves. We present models of six eclipsing binary systems, and for previously known targets, we compare our results with those available in the literature. Masses of binary components have been determined with 3% errors for some targets. We confront our results with benchmark values based on measurements from the HARPS and UCLES spectrographs on 4-m class telescopes and find very good agreement. The described setup is very efficient and well suited for a spectroscopic survey. We can now spectroscopically characterize about 300 eclipsing binary stars per year up to 10.2 mag assuming typical weather conditions at SAAO without a single observing trip.

We evaluate a pre-production BACHES \'{e}chelle spectrograph in terms of its usefulness for radial velocity surveys of binary stars with small telescopes in a remote and autonomous. We use the Solaris-4 observatory located in Casleo, Argentina, that is part of a global network of autonomous observatories as the test-bed for the instrument. The setup is designed in such a way that spectroscopy and photometry can be carried out using the same telescope without the need to mechanically modify the imaging train. We observe single spectroscopic standard stars as well as binary stars up to 9.75 mag. We present results of mechanical tests of the instruments and spectroscopic observations carried out between Nov 26th and Dec 8th 2013. We conclude that BACHES is a very compact and capable spectrograph well suited for remote and autonomous operation. Coupled to a 0.5-m telescope it is capable of obtaining spectra of 10 mag targets with a SNR of 20 for 30-minute exposures. This is a very good result considering the price and size of the instrument. It brings new possibilities to the scientific community and opens a whole new range of research opportunities available to new and existing observatories.

We present a new multi-fluid, grid MHD code PIERNIK, which is based on the Relaxing TVD scheme (Jin & Xin, 1995). The original scheme (see Trac & Pen (2003) and Pen et al. (2003)) has been extended by an addition of dynamically independent, but interacting fluids: dust and a diffusive cosmic ray gas, described within the fluid approximation, with an option to add other fluids in an easy way. The code has been equipped with shearing-box boundary conditions, and a selfgravity module, Ohmic resistivity module, as well as other facilities which are useful in astrophysical fluid-dynamical simulations. The code is parallelized by means of the MPI library. In this paper we present an extension of PIERNIK, which is designed for simulations of diffusive propagation of the Cosmic-Ray (CR) component in the magnetized ISM.

We present new developments on the Cosmic--Ray driven, galactic dynamo, modeled by means of direct, resistive CR--MHD simulations, performed with ZEUS and PIERNIK codes. The dynamo action, leading to the amplification of large--scale galactic magnetic fields on galactic rotation timescales, appears as a result of galactic differential rotation, buoyancy of the cosmic ray component and resistive dissipation of small--scale turbulent magnetic fields. Our new results include demonstration of the global--galactic dynamo action driven by Cosmic Rays supplied in supernova remnants. An essential outcome of the new series of global galactic dynamo models is the equipartition of the gas turbulent energy with magnetic field energy and cosmic ray energy, in saturated states of the dynamo on large galactic scales.

We conduct global galactic--scale magnetohydrodynamical (MHD) simulations of the cosmic--ray driven dynamo. We assume that exploding stars deposit small--scale, randomly oriented, dipolar magnetic fields into the differentially rotating ISM, together with a portion of cosmic rays, accelerated in supernova shocks. Our simulations are performed with the aid of a new parallel MHD code PIERNIK. We demonstrate that dipolar magnetic fields supplied on small SN--remnant scales, can be amplified exponentially by the CR--driven dynamo to the present equipartition values, and transformed simultaneously to large galactic--scales by an inverse cascade promoted by resistive processes.

We present a new multi-fluid, grid MHD code PIERNIK, which is based on the Relaxing TVD scheme (Jin & Xin, 1995). The original scheme (see Trac & Pen (2003) and Pen et al. (2003)) has been extended by an addition of dynamically independent, but interacting fluids: dust and a diffusive cosmic ray gas, described within the fluid approximation, with an option to add other fluids in an easy way. The code has been equipped with shearing-box boundary conditions, and a selfgravity module, Ohmic resistivity module, as well as other facilities which are useful in astrophysical fluid-dynamical simulations. The code is parallelized by means of the MPI library. In this paper we present Ohmic resistivity extension of the original Relaxing TVD MHD scheme, and show examples of magnetic reconnection in cases of uniform and current-dependent resistivity prescriptions.

We present a new multi-fluid, grid MHD code PIERNIK, which is based on the Relaxing TVD scheme (Jin and Xin, 1995). The original scheme (see Trac & Pen (2003) and Pen et al. (2003)) has been extended by an addition of dynamically independent, but interacting fluids: dust and a diffusive cosmic ray gas, described within the fluid approximation, with an option to add other fluids in an easy way. The code has been equipped with shearing-box boundary conditions, and a selfgravity module, Ohmic resistivity module, as well as other facilities which are useful in astrophysical fluid-dynamical simulations. The code is parallelized by means of the MPI library. In this paper we introduce the multifluid extension of Relaxing TVD scheme and present a test case of dust migration in a two-fluid disk composed of gas and dust. We demonstrate that due to the difference in azimuthal velocities of gas and dust and the drag force acting on both components dust drifts towards maxima of gas pressure distribution.

We present a new multi-fluid, grid MHD code PIERNIK, which is based on the Relaxing TVD scheme. The original scheme has been extended by an addition of dynamically independent, but interacting fluids: dust and a diffusive cosmic ray gas, described within the fluid approximation, with an option to add other fluids in an easy way. The code has been equipped with shearing-box boundary conditions, and a selfgravity module, Ohmic resistivity module, as well as other facilities which are useful in astrophysical fluid-dynamical simulations. The code is parallelized by means of the MPI library. In this paper we shortly introduce basic elements of the Relaxing TVD MHD algorithm, following Trac & Pen (2003) and Pen et al. (2003), and then focus on the conservative implementation of the shearing box model, constructed with the aid of the Masset's (2000) method. We present results of a test example of a formation of a gravitationally bounded object (planet) in a self-gravitating and differentially rotating fluid.