Dark energy cosmological models with cosmic string

In this paper we have studied the anisotropic Kantowski-Sachs, locally rotationally symmetric (LRS) Bianchi type-I and LRS Bianchi type-III geometries filled with dark energy and one dimensional cosmic string in the Saez-Ballester theory of gravitation. To get physically valid solution we take hybrid expansion law of the average scale factor which is a product of power and exponential type of functions that results in time dependent deceleration parameter (\(q\)). The equation of state parameter of dark energy (\(\omega _{\mathit{de}}\)) has been discussed and we have observed that for the three models it crosses the phantom divide line (\(\omega _{\mathit{de}} = -1\)) and shows quintom like behavior. The density of dark energy (\(\rho _{\mathit{de}}\)) is an increasing function of redshift and remains positive throughout the evolution of the universe for the three models. Moreover in Kantowski-Sachs and LRS Bianchi type-I geometries the dark energy density dominates the string tension density (\(\lambda \)) and proper density (\(\rho \)) throughout the evolution of the universe. The physical and geometrical aspects of the statefinder parameters (\(r,s\)), squared speed of sound (\(v_{s}^{2} \)) and \(\omega _{\mathit{de}}\)–\(\omega ^{\prime }_{\mathit{de}}\) plane are also discussed.     

Forecasting Solar Flares Using Magnetogram-based Predictors and Machine Learning

We propose a forecasting approach for solar flares based on data from Solar Cycle 24, taken by the Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamics Observatory (SDO) mission. In particular, we use the Space-weather HMI Active Region Patches (SHARP) product that facilitates cut-out magnetograms of solar active regions (AR) in the Sun in near-realtime (NRT), taken over a five-year interval (2012?–?2016). Our approach utilizes a set of thirteen predictors, which are not included in the SHARP metadata, extracted from line-of-sight and vector photospheric magnetograms. We exploit several machine learning (ML) and conventional statistics techniques to predict flares of peak magnitude \({>}\,\mbox{M1}\) and \({>}\,\mbox{C1}\) within a 24 h forecast window. The ML methods used are multi-layer perceptrons (MLP), support vector machines (SVM), and random forests (RF). We conclude that random forests could be the prediction technique of choice for our sample, with the second-best method being multi-layer perceptrons, subject to an entropy objective function. A Monte Carlo simulation showed that the best-performing method gives accuracy \(\mathrm{ACC}=0.93(0.00)\), true skill statistic \(\mathrm{TSS}=0.74(0.02)\), and Heidke skill score \(\mathrm{HSS}=0.49(0.01)\) for \({>}\,\mbox{M1}\) flare prediction with probability threshold 15% and \(\mathrm{ACC}=0.84(0.00)\), \(\mathrm{TSS}=0.60(0.01)\), and \(\mathrm{HSS}=0.59(0.01)\) for \({>}\,\mbox{C1}\) flare prediction with probability threshold 35%.     

On the photo-gravitational restricted four-body problem with variable mass

This paper deals with the photo-gravitational restricted four-body problem (PR4BP) with variable mass. Following the procedure given by Gascheau (C. R. 16:393–394, 1843) and Routh (Proc. Lond. Math. Soc. 6:86–97, 1875), the conditions of linear stability of Lagrange triangle solution in the PR4BP are determined. The three radiating primaries having masses \(m_{1}\), \(m_{2}\) and \(m_{3}\) in an equilateral triangle with \(m_{2}=m_{3}\) will be stable as long as they satisfy the linear stability condition of the Lagrangian triangle solution. We have derived the equations of motion of the mentioned problem and observed that there exist eight libration points for a fixed value of parameters \(\gamma (\frac{m \ \text{at time} \ t}{m \ \text{at initial time}}, 0<\gamma\leq1 )\), \(\alpha\) (the proportionality constant in Jeans’ law (Astronomy and Cosmogony, Cambridge University Press, Cambridge, 1928), \(0\leq\alpha\leq2.2\)), the mass parameter \(\mu=0.005\) and radiation parameters \(q_{i}, (0< q_{i}\leq1, i=1, 2, 3)\). All the libration points are non-collinear if \(q_{2}\neq q_{3}\). It has been observed that the collinear and out-of-plane libration points also exist for \(q_{2}=q_{3}\). In all the cases, each libration point is found to be unstable. Further, zero velocity curves (ZVCs) and Newton–Raphson basins of attraction are also discussed.     

First integrals for the Kepler problem with linear drag

In this work we consider the Kepler problem with linear drag, and prove the existence of a continuous vector-valued first integral, obtained taking the limit as \(t\rightarrow +\infty \) of the Runge–Lenz vector. The norm of this first integral can be interpreted as an asymptotic eccentricity \(e_{\infty }\) with \(0\le e_{\infty } \le 1\). The orbits satisfying \(e_{\infty } <1\) approach the singularity by an elliptic spiral and the corresponding solutions \(x(t)=r(t)e^{i\theta (t)}\) have a norm r(t) that goes to zero like a negative exponential and an argument \(\theta (t)\) that goes to infinity like a positive exponential. In particular, the difference between consecutive times of passage through the pericenter, say \(T_{n+1} -T_n\), goes to zero as \(\frac{1}{n}\).     

Resonance tongues in the linear Sitnikov equation

In this paper, we deal with a Hill’s equation, depending on two parameters \(e\in [0,1)\) and \(\varLambda >0\), that has applications to some problems in Celestial Mechanics of the Sitnikov type. Due to the nonlinearity of the eccentricity parameter e and the coexistence problem, the stability diagram in the \((e,\varLambda )\)-plane presents unusual resonance tongues emerging from points \((0,(n/2)^2),\ n=1,2,\ldots \) The tongues bounded by curves of eigenvalues corresponding to \(2\pi \)-periodic solutions collapse into a single curve of coexistence (for which there exist two independent \(2\pi \)-periodic eigenfunctions), whereas the remaining tongues have no pockets and are very thin. Unlike most of the literature related to resonance tongues and Sitnikov-type problems, the study of the tongues is made from a global point of view in the whole range of \(e\in [0,1)\). Indeed, an interesting behavior of the tongues is found: almost all of them concentrate in a small \(\varLambda \)-interval [1, 9 / 8] as \(e\rightarrow 1^-\). We apply the stability diagram of our equation to determine the regions for which the equilibrium of a Sitnikov \((N+1)\)-body problem is stable in the sense of Lyapunov and the regions having symmetric periodic solutions with a given number of zeros. We also study the Lyapunov stability of the equilibrium in the center of mass of a curved Sitnikov problem.     

De Sitter universe described by a binary mixture with a variable cosmological constant $\lambda $

We have constructed a self-consistent system of Bianchi Type VI₀ cosmology, and mingling of perfect fluid and dark energy in five dimensions. The usual equation of state \(p = \gamma \rho \) with \(\gamma \in [0, 1]\) is chosen by the perfect fluid. The dark energy assumed to be chosen is taken into consideration to be either the quintessence or Chaplygin gas. The same solutions pertaining to the corresponding the field equations of Einstein are obtained as a quadrature. State parameter’s equations for dark energy \(\omega \) is found to be consistent enough with the recent observations of SNe Ia data (SNe Ia data with CMBR anisotropy) and galaxy clustering statistics. Here our models predict that the rate of expansion of Universe would increase with passage of time. The cosmological constant \(\varLambda \) is traced as a declining function of time and it gets nearer to a small positive value later on which is supported by the results from the current supernovae Ia observations. Also a detail discussion is made on the physical and geometrical aspects of the models.     

Geomagnetic Effects of Corotating Interaction Regions

We present an analysis of the geoeffectiveness of corotating interaction regions (CIRs), employing the data recorded from 25 January to 5 May 2005 and throughout 2008. These two intervals in the declining phase of Solar Cycle 23 are characterised by a particularly low number of interplanetary coronal mass ejections (ICMEs). We study in detail how four geomagnetic-activity parameters (the Dst, Ap, and AE indices, as well as the Dst time derivative, \(\mathrm{dDst}/\mathrm{d}t\)) are related to three CIR-related solar wind parameters (flow speed, \(V\), magnetic field, \(B\), and the convective electric field based on the southward Geocentric solar magnetospheric (GSM) magnetic field component, \(\mathit{VB}_{s}\)) on a three-hour time resolution. In addition, we quantify statistical relationships between the mentioned geomagnetic indices. It is found that Dst is correlated best to \(V\), with a correlation coefficient of \(\mathrm{cc}\approx0.6\), whereas there is no correlation between \(\mathrm{dDst}/\mathrm{d}t\) and \(V\). The Ap and AE indices attain peaks about half a day before the maximum of \(V\), with correlation coefficients ranging from \(\mathrm{cc}\approx0.6\) to \(\mathrm{cc}\approx0.7\), depending on the sample used. The best correlations of Ap and AE are found with \(\mathit{VB}_{s}\) with a delay of 3 h, being characterised by \(\mathrm{cc}\gtrsim 0.6\). The Dst derivative \(\mathrm{dDst}/\mathrm{d}t\) is also correlated with \(\mathit{VB}_{s}\), but the correlation is significantly weaker \(\mathrm{cc}\approx 0.4\)?–?0.5, with a delay of 0?–?3 h, depending on the employed sample. Such low values of correlation coefficients indicate that there are other significant effects that influence the relationship between the considered parameters. The correlation of all studied geomagnetic parameters with \(B\) are characterised by considerably lower correlation coefficients, ranging from \(\mathrm{cc}=0.3\) in the case of \(\mathrm{dDst}/\mathrm{d}t\) up to \(\mathrm{cc}=0.56\) in the case of Ap. It is also shown that peak values of geomagnetic indices depend on the duration of the CIR-related structures. The Dst is closely correlated with Ap and AE (\(\mathrm{cc}=0.7\)), Dst being delayed for about 3 h. On the other hand, \(\mathrm{dDst}/\mathrm{d}t\) peaks simultaneously with Ap and AE, with correlation coefficients of 0.48 and 0.56, respectively. The highest correlation (\(\mathrm{cc}=0.81\)) is found for the relationship between Ap and AE.     

Secular resonances between bodies on close orbits: a case study of the Himalia prograde group of jovian irregular satellites

The gravitational interaction between two objects on similar orbits can effect noticeable changes in the orbital evolution even if the ratio of their masses to that of the central body is vanishingly small. Christou (Icarus 174:215–229, 2005) observed an occasional resonant lock in the differential node \(\varDelta \varOmega \) between two members in the Himalia irregular satellite group of Jupiter in the N-body simulations (corresponding mass ratio \(\sim 10^{-9}\)). Using a semianalytical approach, we have reproduced this phenomenon. We also demonstrate the existence of two additional types of resonance, involving angle differences \(\varDelta \omega \) and \(\varDelta (\varOmega +\varpi )\) between two group members. These resonances cause secular oscillations in eccentricity and/or inclination on timescales \(\sim \)1 Myr. We locate these resonances in (a, e, i) space and analyse their topological structure. In subsequent N-body simulations, we confirm these three resonances and find a fourth one involving \(\varDelta \varpi \). In addition, we study the occurrence rates and the stability of the four resonances from a statistical perspective by integrating 1000 test particles for 100 Myr. We find \(\sim \)10 to 30 librators for each of the resonances. Particularly, the nodal resonance found by Christou is the most stable: 2 particles are observed to stay in libration for the entire integration.     

Relative equilibria of the restricted three-body problem in curved spaces

We use a formulation of the N-body problem in spaces of constant Gaussian curvature, \({\kappa }\in \mathbb {R}\), as widely used by A. Borisov, F. Diacu and their coworkers. We consider the restricted three-body problem in \(\mathbb {S}^2\) with arbitrary \({\kappa }>0\) (resp. \(\mathbb {H}^2\) with arbitrary \({\kappa }<0\)) in a formulation also valid for the case \({\kappa }=0\). For concreteness when \({\kappa }>0\) we restrict the study to the case of the three bodies at the upper hemisphere, to be denoted as \(\mathbb {S}^2_+\). The main goal is to obtain the totality of relative equilibria as depending on the parameters \({\kappa }\) and the mass ratio \(\mu \). Several general results concerning relative equilibria and its stability properties are proved analytically. The study is completed numerically using continuation from the \({\kappa }=0\) case and from other limit cases. In particular both bifurcations and spectral stability are also studied. The \(\mathbb {H}^2\) case is similar, in some sense, to the planar one, but in the \(\mathbb {S}^2_+\) case many differences have been found. Some surprising phenomena, like the coexistence of many triangular-like solutions for some values \(({\kappa },\mu )\) and many stability changes will be discussed.     

Density variation effect on multi-ions with kinetic Alfven wave around cusp region—a kinetic approach

The kinetic Alfven waves in the presence of homogeneous magnetic field plasma with multi-ions effect are investigated. The dispersion relation and normalised damping rate are derived for low-\(\beta\) plasma using kinetic theory. The effect of density variation of \(\text{H}^{+}\), \(\text{He}^{+}\) and \(\text{O}^{+}\) ions is observed on frequency and damping rate of the wave. The variation of frequency (\(\omega\)) and normalised damping rate (\(\gamma / \varOmega_{H^{ +}} \)) of the wave are studied with respect to \(k_{ \bot} \rho_{j}\), where \(k_{ \bot} \) is the perpendicular wave number, \(\rho_{j}\) is the ion gyroradius and \(j \) denotes \(\text{H}^{+}\), \(\text{He}^{+}\) and \(\text{O}^{+}\) ions. The variation with \(k_{ \bot} \rho_{j}\) is considered over wide range. The parameters appropriate to cusp region are used for the explanation of results. It is found that with hydrogen and helium ions gyration, the frequency of wave is influenced by the density variation of \(\text{H}^{+}\) and \(\text{He}^{+}\) ions but remains insensitive to the change in density of \(\text{O}^{+}\) ions. For oxygen ion gyration, the frequency of wave varies over a short range only for \(\text{O}^{+}\) ion density variation. The wave shows damping at lower altitude due to variation in density of lighter \(\text{H}^{+}\) and \(\text{He}^{+}\) ions whereas at higher altitude only heavy \(\text{O}^{+}\) ions contribute in wave damping. The damping of wave may be due to landau damping or energy transfer from wave to particles. The present study signifies that the both lighter and heavier ions dominate differently to change the characteristics of kinetic Alfven wave and density variation is also an important parameter to understand wave phenomena in cusp region.     

New holographic dark energy model with constant bulk viscosity in modified $$ gravity theory

The aim of this paper is to study new holographic dark energy (HDE) model in modified \(f(R,T)\) gravity theory within the framework of a flat Friedmann-Robertson-Walker model with bulk viscous matter content. It is thought that the negative pressure caused by the bulk viscosity can play the role of dark energy component, and drive the accelerating expansion of the universe. This is the motive of this paper to observe such phenomena with bulk viscosity. In the specific model \(f(R,T)=R+\lambda T\), where \(R\) is the Ricci scalar, \(T\) the trace of the energy-momentum tensor and \(\lambda \) is a constant, we find the solution for non-viscous and viscous new HDE models. We analyze new HDE model with constant bulk viscosity, \(\zeta =\zeta _{0}= \text{const.}\) to explain the present accelerated expansion of the universe. We classify all possible scenarios (deceleration, acceleration and their transition) with possible positive and negative ranges of \(\lambda \) over the constraint on \(\zeta _{0}\) to analyze the evolution of the universe. We obtain the solutions of scale factor and deceleration parameter, and discuss the evolution of the universe. We observe the future finite-time singularities of type I and III at a finite time under certain constraints on \(\lambda \). We also investigate the statefinder and \(\mathit{Om}\) diagnostics of the viscous new HDE model to discriminate with other existing dark energy models. In late time the viscous new HDE model approaches to \(\varLambda \mathit{CDM}\) model. We also discuss the thermodynamics and entropy of the model and find that it satisfies the second law of thermodynamics.     

Relating Solar Energetic Particle Event Fluences to Peak Intensities

Recently we (Kahler and Ling, Solar Phys.292, 59, 2017: KL) have shown that time–intensity profiles [\(I(t)\)] of 14 large solar energetic particle (SEP) events can be fitted with a simple two-parameter fit, the modified Weibull function, which is characterized by shape and scaling parameters [\(\alpha\) and \(\beta\)]. We now look for a simple correlation between an event peak energy intensity [\(I_{\mathrm{p}}\)] and the time integral of \(I(t)\) over the event duration: the fluence [\(F\)]. We first ask how the ratio of \(F/I_{\mathrm{p}}\) varies for the fits of the 14 KL events and then examine that ratio for three separate published statistical studies of SEP events in which both \(F\) and \(I_{\mathrm{p}}\) were measured for comparisons of those parameters with various solar-flare and coronal mass ejection (CME) parameters. The three studies included SEP energies from a 4?–?13 MeV band to \(E > 100~\mbox{MeV}\). Within each group of SEP events, we find a very robust correlation (\(\mathrm{CC} > 0.90\)) in log–log plots of \(F\)versus\(I_{\mathrm{p}}\) over four decades of \(I_{\mathrm{p}}\). The ratio increases from western to eastern longitudes. From the value of \(I_{\mathrm{p}}\) for a given event, \(F\) can be estimated to within a standard deviation of a factor of \({\leq}\,2\). Log–log plots of two studies are consistent with slopes of unity, but the third study shows plot slopes of \({<}\,1\) and decreasing with increasing energy for their four energy ranges from \(E > 10~\mbox{MeV}\) to \({>}\,100~\mbox{MeV}\). This difference is not explained.     

Landau damping of dust acoustic waves in the presence of hybrid nonthermal nonextensive electrons

Based on the kinetic theory, Landau damping of dust acoustic waves (DAWs) propagating in a dusty plasma composed of hybrid nonthermal nonextensive distributed electrons, Maxwellian distributed ions and negatively charged dust grains is investigated using Vlasov-Poisson’s equations. The characteristics of the DAWs Landau damping are discussed. It is found that the wave frequency increases by decreasing (increasing) the value of nonextensive (nonthermal) parameter, \(q\) (\(\alpha \)). It is recognized that \(\alpha \) plays a significant role in observing damping or growing DAW oscillations. For small values of \(\alpha \), damping modes have been observed until reaching a certain value of \(\alpha \) at which \(\omega _{i}\) vanishes, then a growing mode appears in the case of superextensive electrons. However, only damping DAW modes are observed in case of subextensive electrons. The present study is useful in the space situations where such distribution exists.     

Geodesics around oscillatons made of exponential scalar field potential

Some of the spherically symmetric solutions to the Einstein–Klein–Gordon (EKG) equations can describe the astronomical soliton objects made of a real time-dependent scalar fields. The solutions are known as oscillatons which are non-singular satisfying the flatness conditions asymptotically with periodic (separated) time-dependency. In this paper, we investigate the geodesic motion around an oscillaton. The Spherically Symmetric Geometry allows the bound orbits in the plan \(\theta=\frac{\pi}{2}\) under a given initial conditions. The potential for the scalar field \(\varPhi=\varPhi(r,t)\), is an exponential function of the form \(V(\varPhi)=V_{0}\exp(\lambda\sqrt{k_{0}}\varPhi)\).     

Anisotropic new holographic dark energy model in Saez–Ballester theory of gravitation

This work is devoted to the investigation of new holographic dark energy (infrared cutoff is the Hubble radius) in locally rotationally symmetric Bianchi type-\(I\) universe within the framework of Saez–Ballester (Phys. Lett. A 113:467, 1986) scalar–tensor theory of gravitation. We construct interacting and non-interacting dark energy models by solving the field equations using a relationship between the metric potentials. This leads to a variable deceleration parameter model which exhibits a transition of the universe from deceleration to acceleration. We evaluate various cosmological parameters of our models. We have observed that the energy density parameters, equation of state and important cosmological planes (\(\omega _{\mathit{de}} - \omega _{\mathit{de}}'\) and \(r - s\)) yield the results compatible with the modern observational data. We have, also, discussed the stability analysis of our models.     

An analysis of the convergence of Newton iterations for solving elliptic Kepler’s equation

In this note a study of the convergence properties of some starters \( E_0 = E_0(e,M)\) in the eccentricity–mean anomaly variables for solving the elliptic Kepler’s equation (KE) by Newton’s method is presented. By using a Wang Xinghua’s theorem (Xinghua in Math Comput 68(225):169–186, 1999) on best possible error bounds in the solution of nonlinear equations by Newton’s method, we obtain for each starter \( E_0(e,M)\) a set of values \( (e,M) \in [0, 1) \times [0, \pi ]\) that lead to the q-convergence in the sense that Newton’s sequence \( (E_n)_{n \ge 0}\) generated from \( E_0 = E_0(e,M)\) is well defined, converges to the exact solution \(E^* = E^*(e,M)\) of KE and further \( \vert E_n - E^* \vert \le q^{2^n -1}\; \vert E_0 - E^* \vert \) holds for all \( n \ge 0\). This study completes in some sense the results derived by Avendaño et al. (Celest Mech Dyn Astron 119:27–44, 2014) by using Smale’s \(\alpha \)-test with \(q=1/2\). Also since in KE the convergence rate of Newton’s method tends to zero as \( e \rightarrow 0\), we show that the error estimates given in the Wang Xinghua’s theorem for KE can also be used to determine sets of q-convergence with \( q = e^k \; \widetilde{q} \) for all \( e \in [0,1)\) and a fixed \( \widetilde{q} \le 1\). Some remarks on the use of this theorem to derive a priori estimates of the error \( \vert E_n - E^* \vert \) after n Kepler’s iterations are given. Finally, a posteriori bounds of this error that can be used to a dynamical estimation of the error are also obtained.     

Timing irregularities of PSR J1705–1906

Timing analysis of PSR J1705–1906 using data from Nanshan 25-m and Parkes 64-m radio telescopes, which span over fourteen years, shows that the pulsar exhibits significant proper motion, and rotation instability. We updated the astrometry parameters and the spin parameters of the pulsar. In order to minimize the effect of timing irregularities on measuring its position, we employ the Cholesky method to analyse the timing noise. We obtain the proper motion of \(-77(3)\) mas?yr^?1 in right ascension and \(-38(29)\) mas?yr^?1 in declination. The power spectrum of timing noise is analyzed for the first time, which gives the spectral exponent \(\alpha =-5.2\) for the power-law model indicating that the fluctuations in spin frequency and spin-down rate dominate the red noise. We detect two small glitches from this pulsar with fractional jump in spin frequency of \(\Delta \nu /\nu \sim 2.9 \times 10^{-10}\) around MJD 55199 and \(\Delta \nu /\nu \sim 2.7\times 10^{-10}\) around MJD 55953. Investigations of pulse profile at different time segments suggest no significant changes in the pulse profiles around the two glitches.     

A new holographic dark energy model in Brans-Dicke theory with logarithmic scalar field

We study a holographic dark energy model in the framework of Brans-Dicke (BD) theory with taking into account the interaction between dark matter and holographic dark energy. We use the recent observational data sets, namely SN Ia compressed Joint Light-Analysis (cJLA) compilation, Baryon Acoustic Oscillations (BAO) from BOSS DR12 and the Cosmic Microwave Background (CMB) of Planck 2015. After calculating the evolution of the equation of state as well as the deceleration parameters, we find that with a logarithmic form for the BD scalar field the phantom crossing can be achieved in the late time of cosmic evolution. Unlike the conventional theory of holographic dark energy in standard cosmology (\(\omega_{D}=0\)), our model results in a late time accelerated expansion. It is also shown that the cosmic coincidence problem may be resolved in the proposed model. We execute the statefinder and Om diagnostic tools and demonstrate that interaction term does not play a significant role. Based on the observational data sets used in this paper it seems that the best value with \(1\sigma \) and \(2\sigma \) confidence interval are \(\varOmega_{m}=0.268^{+0.008~+0.010}_{-0.007~-0.009}\), \(\alpha =3.361^{+0.332~+0.483} _{-0.401~-0.522}\), \(\beta =5.560^{+0.541~+0.780}_{-0.510~-0.729}\), \(c=0.777^{+0.023~+0.029}_{-0.017~-0.023}\) and \(b^{2} =0.045\), according to which we find that the proposed model in the presence of interaction is compatible with the recent observational data.     

<Emphasis Type="Italic">GALEX</Emphasis> far-ultraviolet observations of metal-poor subdwarfs

Far-ultraviolet photometry derived from the GALEX satellite observatory has been compiled for a sample of metal-poor subdwarfs with \(\mathrm{[Fe/H]} < -1.0\). The FUV properties of these subdwarfs are compared with those of a set of Population I dwarfs that are known to have low levels of chromospheric activity. Comparisons are made via a number of photometric plots, including an absolute FUV magnitude versus \((V-K_{s})\) diagram, two-colour diagrams involving both \((m_{ \mathrm{FUV}}-B)\) and \((m_{\mathrm{FUV}}-V)\) versus \(B-V\), and a two-colour diagram composed of \((m_{\mathrm{FUV}}-V)\) versus \((V-K_{s})\). The warmest subdwarfs with \((V-K_{s}) \sim1.2\mbox{--}1.4\) show FUV excesses ranging from \(\sim2\mbox{--}3~\mbox{mag}\) relative to the Population I dwarfs, with the amount of FUV enhancement decreasing among subdwarfs of decreasing effective temperature. The coolest dwarfs that are compared have \((V-K_{s}) \sim1.8\), and among these stars the subdwarfs with \(-2.0 \leq{\mathrm{[Fe/H]}} \leq-1.0\) approach the locus of low activity Population I dwarfs in the \((m_{\mathrm{FUV}}-V, V-K_{s})\) diagram. In the \((m_{\mathrm{FUV}}-B, B-V)\) diagram the subdwarfs in this metallicity range overlap the Population I dwarf sequence for \((B-V) > 0.6\). The behaviour of the subdwarfs is consistent with their FUV fluxes being determined by a combination of a photospheric FUV spectrum, the strength of which diminishes towards cooler effective temperatures, and a spectrum of emission lines arising from a chromosphere and/or transition region which are of comparable strength between the coolest dwarfs and subdwarfs.