We present modifications to the Athena++ framework to enable use of general equations of state (EOS). Part of our motivation for doing so is to model transient astrophysics phenomena, as these types of events are often not well approximated by an ideal gas. This necessitated changes to the Riemann solvers implemented in Athena++. We discuss the adjustments made to the HLLC, and HLLD solvers and EOS calls required for arbitrary EOS. We demonstrate the reliability of our code in a number of tests which utilize a relatively simple, but non-trivial EOS based on hydrogen ionization, appropriate for the transition from atomic to ionized hydrogen. Additionally, we perform tests using an electron-positron Helmholtz EOS, appropriate for regimes where nuclear statistical equilibrium is a good approximation. These new complex EOS tests overall show that our modifications to Athena++ accurately solve the Riemann problem with linear convergence and linear-wave tests with quadratic convergence. We provide our test solutions as a means to check the accuracy of other hydrodynamic codes. Our tests and additions to Athena++ will enable further research into (magneto)hydrodynamic problems where realistic treatments of the EOS are required.

We present new methods to solve the Riemann problem both exactly and approximately for general equations of state (EoS) to facilitate realistic modeling and understanding of astrophysical flows. The existence and uniqueness of the new exact general EoS Riemann solution can be guaranteed if the EoS is monotone regardless of the physical validity of the EoS. We confirm that: (1) the solution of the new exact general EoS Riemann solver and the solution of the original exact Riemann solver match when calculating perfect gas Euler equations; (2) the solution of the new Harten-Lax-van Leer-Contact (HLLC) general EoS Riemann solver and the solution of the original HLLC Riemann solver match when working with perfect gas EoS; and (3) the solution of the new HLLC general EoS Riemann solver approaches the new exact solution. We solve the EoS with two methods, one is to interpolate 2D EoS tables by the bi-linear interpolation method, and the other is to analytically calculate thermodynamic variables at run-time. The interpolation method is more general as it can work with other monotone and realistic EoS while the analytic EoS solver introduced here works with a relatively idealized EoS. Numerical results confirm that the accuracy of the two EoS solvers is similar. We study the efficiency of these two methods with the HLLC general EoS Riemann solver and find that analytic EoS solver is faster in the test problems. However, we point out that a combination of the two EoS solvers may become favorable in some specific problems. Throughout this research, we assume local thermal equilibrium.

We show that recent observations of the compact binary, AM CVn type system, ES Ceti are fully consistent with theoretical predictions of stable mass transfer moderated by angular momentum loss due to gravitational-wave radiation. One of the main predictions of this model (for degenerate donors) is a widening of the binary. The mass transfer rate inferred from the observed rate of change in the orbital frequency is consistent with that inferred from the observed flux using the recent Gaia DR2 parallax

We review the abundance and diversity of terrestrial rock hosted life, the environments it inhabits, the evolution of its metabolisms, and its fossil biomarkers to provide guidance in the search for the biomarkers of rock hosted life on Mars. Key finds are metabolic pathways for chemolithotrophic microorganisms evolved much earlier in Earth history than those of surface dwelling phototrophic microorganisms,the emergence of the former occurred at a time when Mars was habitable, whereas that of the latter occurred at a time when the martian surface would have been uninhabitable, subsurface biomass do not correlate with organic carbon and tends to be highest at interfaces where chemical redox gradients are most pronounced, deep subsurface metabolic activity does not rely upon the respiration of organic photosynthate but upon the flux of inorganic energy and the abiotic and biotic recycling of metabolic waste products, and the rock record reveals examples of subsurface life back to 3.45 Ga with several examples of good preservation potential in rock types that are quite different from those preserving the photospheric supported biosphere.These findings suggest that rock hosted life would have likely to emerge and be preserved in a martian context. We thus propose a Mars exploration strategy that scales spatially, focusing initially on identifying rocks with evidence for groundwater flow and low temperature mineralization, then identifying redox and permeability interfaces preserved within rock outcrops, and finally focusing on finding minerals associated with redox reactions and traces of carbon and diagnostic biosignatures. The lessons from Earth show that ancient rock hosted life is preserved in the fossil record and confirmable via a suite of morphologic, organic, mineralogical and isotopic fingerprints and microscopic textures.

We present the results of local, vertically stratified, radiation magnetohydrodynamic shearing box simulations of magnetorotational instability (MRI) turbulence for a (hydrogen poor) composition applicable to accretion disks in AM CVn type systems. Many of these accreting white dwarf systems are helium analogues of dwarf novae (DNe). We utilize frequency-integrated opacity and equation of state tables appropriate for this regime to accurately portray the relevant thermodynamics. We find bistability of thermal equilibria in the effective temperature, surface mass density plane typically associated with disk instabilities. Along this equilibrium curve (i.e. the S-curve) we find that the stress to thermal pressure ratio $\alpha$ varied with peak values of $\sim 0.15$ near the tip of the upper branch. Similar to DNe, we found enhancement of $\alpha$ near the tip of the upper branch caused by convection; this increase in $\alpha$ occurred despite our choice of zero net vertical magnetic flux. Two notable differences we find between DN and AM CVn accretion disk simulations are that AM CVn disks are capable of exhibiting persistent convection in outburst, and ideal MHD is valid throughout quiescence for AM CVns. In contrast, DNe simulations only show intermittent convection, and non-ideal MHD effects are likely important in quiescence. By combining our previous work with these new results, we also find that convective enhancement of the MRI is anticorrelated with mean molecular weight.

Vertically stratified shearing box simulations of magnetorotational turbulence commonly exhibit a so-called butterfly diagram of quasi-periodic azimuthal field reversals. However, in the presence of hydrodynamic convection, field reversals no longer occur. Instead, the azimuthal field strength fluctuates quasi-periodically while maintaining the same polarity, which can either be symmetric or antisymmetric about the disc midplane. Using data from the simulations of Hirose et al. (2014), we demonstrate that the lack of field reversals in the presence of convection is due to hydrodynamic mixing of magnetic field from the more strongly magnetized upper layers into the midplane, which then annihilate field reversals that are starting there. Our convective simulations differ in several respects from those reported in previous work by others, in which stronger magnetization likely plays a more important role than convection.

The phenomenological Disc Instability Model has been successful in reproducing the observed light curves of dwarf nova outbursts by invoking an enhanced Shakura-Sunyaev $\alpha$ parameter $\sim0.1-0.2$ in outburst compared to a low value $\sim0.01$ in quiescence. Recent thermodynamically consistent simulations of magnetorotational (MRI) turbulence with appropriate opacities and equation of state for dwarf nova accretion discs have found that thermal convection enhances $\alpha$ in discs in outburst, but only near the hydrogen ionization transition. At higher temperatures, convection no longer exists and $\alpha$ returns to the low value comparable to that in quiescence. In order to check whether this enhancement near the hydrogen ionization transition is sufficient to reproduce observed light curves, we incorporate this MRI-based variation in $\alpha$ into the Disc Instability Model, as well as simulation-based models of turbulent dissipation and convective transport. These MRI-based models can successfully reproduce observed outburst and quiescence durations, as well as outburst amplitudes, albeit with different parameters from the standard Disc Instability Models. The MRI-based model lightcurves exhibit reflares in the decay from outburst, which are not generally observed in dwarf novae. However, we highlight the problematic aspects of the quiescence physics in the Disc Instability Model and MRI simulations that are responsible for this behavior.

We present the results of local, vertically stratified, radiation MHD shearing box simulations of MRI turbulence appropriate for the hydrogen ionizing regime of dwarf nova and soft X-ray transient outbursts. We incorporate the frequency-integrated opacities and equation of state for this regime, but neglect non-ideal MHD effects and surface irradiation, and do not impose net vertical magnetic flux. We find two stable thermal equilibrium tracks in the effective temperature versus surface mass density plane, in qualitative agreement with the S-curve picture of the standard disk instability model. We find that the large opacity at temperatures near $10^4$K, a corollary of the hydrogen ionization transition, triggers strong, intermittent thermal convection on the upper stable branch. This convection strengthens the magnetic turbulent dynamo and greatly enhances the time-averaged value of the stress to thermal pressure ratio $\alpha$, possibly by generating vertical magnetic field that may seed the axisymmetric magnetorotational instability, and by increasing cooling so that the pressure does not rise in proportion to the turbulent dissipation. These enhanced stress to pressure ratios may alleviate the order of magnitude discrepancy between the $\alpha$-values observationally inferred in the outburst state and those that have been measured from previous local numerical simulations of magnetorotational turbulence that lack net vertical magnetic flux.

Based on a deep imaging survey, we present the first homogeneous star formation history (SFH) of the Fornax dwarf spheroidal (dSph) galaxy. We have obtained two-filter photometry to a depth of B ~ 23 over the entire surface of Fornax, the brightest dSph associated with the Milky Way, and derived its SFH using a CMD-fitting technique. We show that Fornax has produced the most complex star formation and chemical enrichment histories of all the Milky Way dSphs. This system has supported multiple epochs of star formation. A significant number of stars were formed in the early Universe, however the most dominant population are the intermediate age stars. This includes a strong burst of star formation approximately 3 to 4 Gyr ago. Significant population gradients are also evident. Similar to other dSphs, we have found that recent star formation was concentrated towards the centre of the system. Furthermore, we show that the central region harboured a faster rate of chemical enrichment than the outer parts of Fornax. At the centre, the ancient stars (age > 10 Gyr) display a mean metallicity of [Fe/H] ~ -1.4, with evidence for three peaks in the metallicity distribution. Overall, enrichment in Fornax has been highly efficient: the most recent star formation burst has produced stars with close to solar metallicity. Our results support a scenario in which Fornax experienced an early phase of rapid chemical enrichment producing a wide range of abundances. Star formation gradually decreased until ~4 Gyr ago, when Fornax experienced a sudden burst of strong star formation activity accompanied by substantial chemical enrichment. Weaker star forming events followed, and we have found tentative evidence for stars with ages less than 100 Myr.

We have used the two-degree field (2dF) multi-fibre spectrograph of the Anglo-Australian Telescope to search for candidate members of the unusual globular cluster omega Centauri at and beyond the cluster tidal radius. Velocities with an accuracy of ~10 kms-1 were obtained for 4105 stars selected to lie in the vicinity of the lower giant branch in the cluster colour-magnitude diagram and which cover an area on the sky of ~2.4x3.9 deg2 centered on the cluster. Within the velocity interval 190-270 kms-1, the cluster member candidates have a steeply declining surface density distribution consistent with the adopted tidal radius of 57`. For stars in the sample beyond the tidal radius, an analysis of line-strengths from the spectra, as well as radial velocities, identifies only six stars as possible candidates for extra-tidal association with the cluster. If all six of these stars are indeed related to the cluster, then a maximum of 0.7 +/- 0.2 per cent of the total cluster mass is contained in the region between one and two tidal radii. Given this limit, we conclude that there is no compelling evidence for any significant extra-tidal population in omega Cen. The effects of tidal shocks on the outer parts of the cluster are consistent with this limit. Theories for the origin of omega Cen frequently suggest that the cluster is the former nucleus of a tidally stripped dwarf galaxy. Our results require that the stripping process must have been largely complete at early epochs, consistent with dynamical models of the process. The stripped former dwarf galaxy stars are therefore now widely distributed around the Galaxy.

We present deep, wide-field g and r photometry of the transition type dwarf galaxy Leo T, obtained with the blue arm of the Large Binocular Telescope. The data confirm the presence of both very young (<1 Gyr) as well as much older (>5 Gyr) stars. We study the structural properties of the old and young stellar populations by preferentially selecting either population based on their color and magnitude. The young population is significantly more concentrated than the old population, with half-light radii of 104+-8 and 148+-16 pc respectively, and their centers are slightly offset. Approximately 10% of the total stellar mass is estimated to be represented by the young stellar population. Comparison of the color-magnitude diagram (CMD) with theoretical isochrones as well as numerical CMD-fitting suggest that star formation began over 10 Gyr ago and continued in recent times until at least a few hundred Myr ago. The CMD-fitting results are indicative of two distinct star formation bursts, with a quiescent period around 3 Gyr ago, albeit at low significance. The results are consistent with no metallicity evolution and [Fe/H] ~ -1.5 over the entire age of the system. Finally, the data show little if any sign of tidal distortion of Leo T.

We present the first deep color-magnitude diagram of the Canes Venatici I (CVnI) dwarf galaxy from observations with the wide field Large Binocular Camera on the Large Binocular Telescope. Reaching down to the main-sequence turnoff of the oldest stars, it reveals a dichotomy in the stellar populations of CVnI: it harbors an old (> 10 Gyr), metal-poor ([Fe/H] ~ -2.0) and spatially extended population along with a much younger (~ 1.4-2.0 Gyr), 0.5 dex more metal-rich, and spatially more concentrated population. These young stars are also offset by 64_{-20}^{+40} pc to the East of the galaxy center. The data suggest that this young population, which represent ~ 3-5 % of the stellar mass of the galaxy within its half-light radius, should be identified with the kinematically cold stellar component found by Ibata et al. (2006). CVnI therefore follows the behavior of the other remote MW dwarf spheroidals which all contain intermediate age and/or young populations: a complex star formation history is possible in extremely low-mass galaxies.

We explore how well crowded field point-source photometry can be accomplished with SDSS data: We present a photometric pipeline based on DoPhot, and tuned for analyzing crowded-field images from the SDSS. Using Monte Carlo simulations we show that the completeness of source extraction is above 80% to i < 21 (AB) and a stellar surface density of about 200 sq.amin. Hence, a specialized data pipeline can efficiently be used for e.g. nearby resolved galaxies in SDSS images, where the standard SDSS photometric package Photo, when applied in normal survey mode, gives poor results. We apply our pipeline to an area of about 3.55sq.deg. around the dwarf spheroidal galaxy (dSph) Leo I, and construct a high S/N star-count map of Leo I via an optimized filter in color-magnitude space (g,r,i). Although the radial surface-density profile of the dwarf deviates from the best fit empirical King model towards outer radii, we find no evidence for tidal debris out to a stellar surface-density of 4*10^(-3) of the central value. We determine the total luminosity of Leo I, and model its mass using the spherical and isotropic Jeans equation. Assuming that 'mass follows light' we constrain a lower limit of the total mass of the dSph to be (1.7+/-0.2)*10^7 Msol. Contrary, if the mass in Leo I is dominated by a constant density dark-matter (DM) halo, then the mass within the central 12' is (2+/-0.6)*10^8 Msol. This leads to a mass-to-light ratio of >>6 (Ic_sol), and possibly >75 if the DM halo dominates the mass and extends further out than 12'. In summary, our results show that Leo I is a symmetric, relaxed and bound system; this supports the idea that Leo I is a dark-matter dominated system.

This contribution describes photometry for two Galactic dSphs obtained with the Large Binocular Telescope to a magnitude of ~25.5. Using the Large Binocular Camera, a purpose-built wide-field imager for the LBT, we have examined the structure and star formation histories of two newly-discovered Local Group members, the Hercules dSph and the Leo T dSph/dIrr system. We have constructed a structural map for the Hercules system using three-filter photometry to V ~ 25.5. This is the first deep photometry for this system, and it indicates that Hercules is unusually elongated, possibly indicating distortion due to the Galactic tidal field. We have also derived the first star formation history for the Leo T system, and find that its oldest population of stars (age ~ 13 Gyr) were relatively metal-rich, with [Fe/H] ~ -1.5.

Using SDSS I data, we have analysed the stellar distribution of the Leo II dwarf spheroidal galaxy (distance of 233 kpc) to search for evidence of tidal deformation. The existing SDSS photometric catalogue contains gaps in regions of high stellar crowding, hence we filled the area at the centre of Leo II using the DAOPHOT algorithm applied to the SDSS images. The combined DAOPHOT-SDSS dataset contains three-filter photometry over a 4x4 square degree region centred on Leo II. By defining a mask in three-filter colour-magnitude space, we removed the majority of foreground field stars. We have measured the following Leo II structural parameters: a core radius of r_c = 2.64 +/- 0.19 arcmin (178 +/- 13 pc), a tidal radius of r_t = 9.33 +/- 0.47 arcmin (632 +/- 32 pc) and a total V-band luminosity of L_V = (7.4 +/- 2.0) times 10^5 L_sun (M_V = -9.9 +/- 0.3). Our comprehensive analysis of the Leo II structure did not reveal any significant signs of tidal distortion. The internal structure of this object contains only mild isophotal twisting. A small overdensity was discovered appoximately 4.5 tidal radii from the Leo II centre, however we conclude it is unlikely to be material tidally stripped from Leo II based on its stellar population, and is most likely a foreground overdensity of stars. Our results indicate that the influence of the Galactic graviational field on the structure of Leo II has been relatively mild. We rederived the mass-to-light ratio of this system using existing kinematic data combined with our improved structural measurements, and favour the scenario in which Leo II is strongly dominated by dark matter with (M/L)_V ~ 100 in solar units.

We present a deep, wide-field photometric survey of the newly-discovered Hercules dwarf spheroidal galaxy, based on data from the Large Binocular Telescope. Images in B, V and r were obtained with the Large Binocular Camera covering a 23' times 23' field of view to a magnitude of ~25.5 (5 sigma). This permitted the construction of colour-magnitude diagrams that reach approximately 1.5 magnitudes below the Hercules main sequence turnoff. Three-filter photometry allowed us to preferentially select probable Hercules member stars, and examine the structure of this system at a previously unattained level. We find that the Hercules dwarf is highly elongated (3:1), considerably more so than any other dSph satellite of the Milky Way except the disrupting Sagittarius dwarf. While we cannot rule out that the unusual structure is intrinsic to Hercules as an equilibrium system, our results suggest tidal disruption as a likely cause of this highly elliptical structure. Given the relatively large Galactocentric distance of this system (132 +/- 12 kpc), signs of tidal disruption would require the Hercules dwarf to be on a highly eccentric orbit around the Milky Way.

We have imaged the recently discovered stellar overdensity located approximately one core radius from the center of the Fornax dwarf spheroidal galaxy using the Magellan Clay 6.5m telescope with the Magellan Instant Camera (MagIC). Superb seeing conditions allowed us to probe the stellar populations of this overdensity and of a control field within Fornax to a limiting magnitude of R=26. The color-magnitude diagram of the overdensity field is virtually identical to that of the control field with the exception of the presence of a population arising from a very short (less than 300 Myr in duration) burst of star formation 1.4 Gyr ago. Coleman et al. have argued that this overdensity might be related to a shell structure in Fornax that was created when Fornax captured a smaller galaxy. Our results are consistent with this model, but we argue that the metallicity of this young component favors a scenario in which the gas was part of Fornax itself.

The results of a wide-field survey of the Sculptor dwarf spheroidal galaxy are presented. Our aims were to obtain an accurate map of the outer structure of Sculptor, and to determine the level of interaction between this system and the Galaxy. Photometry was obtained in two colours down to the magnitude limits of V=20 and I=19, covering a 3.1 times 3.1 square deg area centred on Sculptor. The resulting colour-magnitude data were used as a mask to select candidate horizontal branch and red giant branch stars for this system. Previous work has shown that the red horizontal branch (HB) stars are more concentrated than the blue HB stars. We have determined the radial distributions of these two populations and show that the overall Sculptor density profile is well described by a two component model based on a combination of these radial distributions. Additionally, spectra of the Ca ii triplet region were obtained for over 700 candidate red giant stars over the 10 square deg region using the 2dF instrument on the Anglo-Australian Telescope. These spectra were used to remove foreground Galactic stars based on radial velocity and Ca ii triplet strength. The final list of Sculptor members contained 148 stars, seven of which are located beyond the nominal tidal radius. Both the photometric and spectroscopic datasets indicate no significant extra-tidal structure. These results support at most a mild level of interaction between this system and the Galaxy, and we have measured an upper mass limit for extra-tidal material to be 2.3 +/- 0.6% of the Sculptor luminous mass. This lack of tidal interaction indicates that previous velocity dispersion measurements (and hence the amount of dark matter detected) in this system are not strongly influenced by the Galactic tidal field.

The results of a photometric survey of the Fornax dwarf spheroidal galaxy are presented herein. Accurate photometry in two colours (V and I) was collected over a 10 square degree area centred on the Fornax system with the aim of searching for extra-tidal structure. The data were complete to a magnitude of V=20, or approximately one magnitude brighter than the Fornax red clump stars. Stars were selected with a colour and magnitude lying near the Fornax red giant branch, thereby reducing contamination from the field star population. We were thus able to probe the outer structure of this dwarf galaxy. Immediately visible was a shell-like structure located 1.3 degrees NW from the centre of Fornax, approximately 30 arcmin beyond the nominal tidal radius at this position angle. We have measured the absolute visual magnitude of this feature to be M_V ~ -7. The feature is aligned with a previously reported shell (age approximately 2 Gyr) located near the core radius of Fornax. A statistical analysis of the extra-tidal stellar distribution further revealed two lobes situated on the Fornax minor axis which are aligned with the two shell-like features. The two-lobed structure combined with the two shells provide strong evidence that Fornax has experienced a merger event in the relatively recent past.

In the search for tidal structure in Galactic satellite systems, we have conducted a photometric survey over a 10 square degree area centred on the Fornax dSph galaxy. The survey was made in two colours, and the resulting colour-magnitude data were used as a mask to select candidate Fornax RGB stars, thereby increasing the contrast of Fornax stars to background sources in the outer regions. Previously, we reported the presence of a shell (age 2 Gyr) located towards the centre of Fornax. In this contribution we reveal a second shell, significantly larger than the first, located 1.3 degrees NW from the centre of Fornax, outside the nominal tidal radius. Moreover, the distribution of Fornax RGB stars reveals two lobes extending to the spatial limit of our survey, and aligned with the minor axis and with the two shells. These results support the hypothesis of a merger between Fornax and a gas-rich companion approximately 2 Gyr ago.

Recent years have seen a series of large-scale photometric surveys with the aim of detecting substructure in nearby dwarf galaxies. Some of these objects display a varying distribution of each stellar population, reflecting their star formation histories. Also, dwarf galaxies are dominated by dark matter, therefore luminous substructure may represent a perturbation in the underlying dark material. Substructure can also be the effect of tidal interaction, such as the disruption of the Sagittarius dSph by the Milky Way. Therefore, substructure in dwarf galaxies manifests the stellar, structural and kinematic evolution of these objects.

We present initial results from a wide field photometric survey of the Fornax dwarf spheroidal (dSph) galaxy. The focus here is on a small overdensity of stars located approximately 17 arcmin (670 pc in projection) from the centre of Fornax. Based on imaging data in both V and I bands down to approximately V = 21, the dimensions of the feature are approximately 1.7 arcmin by 3.2 arcmin (68 pc by 128 pc) with an average surface brightness of 25.8 mag/arcsec^2 in V. A follow-up analysis using the deep B and R band data obtained by Stetson et al. (1998) indicates that the overdensity is apparently dominated by a relatively young stellar population with an age of approximately 2 Gyr, though the current data do not rule out the presence of older stars. Our preferred interpretation is that the overdensity represents shell structure in the Fornax dSph galaxy, a phenomenon previously unseen in dwarf galaxies. The shell may be the remnant of a merger with Fornax of a smaller, gas-rich system that occured approximately 2 Gyr ago. The comparitively recent star formation seen in some dSph and dE galaxies may then also be the result of infall of small relatively gas-rich systems.

We obtained B and R-band CCD images for the dwarf elliptical (dE) galaxies DDO44, UGC4998, KK98_77, DDO71, DDO113, and UGC7356 at the NOT. Using Fourier analysis technique we measure stellar R-band surface brightness fluctuations (SBFs) and magnitudes in 29 different fields of the galaxies. Independent tip of the red giant branch distances for DDO44, KK98_77, DDO71 are used to convert their set of apparent into absolute SBF magnitudes. The results are combined with the local (B-R) colours and compared with the (B-R)-\bar{M}_R relation for mainly old, metal-poor stellar populations as predicted by Worthey's population synthesis models using Padova isochrones. While the colour dependency of the theoretical relation is confirmed by the empirical data, we find a systematic zero point offset between observations and theory in the sense that models are too faint by 0.13+-0.02 mag. Based on these findings we establish a new semiempirical calibration of the SBF method as distance indicator for dE galaxies with an estimated uncertainty of ~10 percent. Taking first advantage of the improved calibration, we determine SBF distances for the other three early-type dwarfs. Although found in the M81 group region, previous velocity measurements suggested UGC4998 is in the background. This picture is confirmed by our SBF distance of 10.5+-0.9 Mpc. We can further identify DDO113 as a faint stellar system at the near side of the Canes Venatici I (CVn I) cloud at a distance of 3.1+-0.3 Mpc. The second CVn I member in our sample, UGC7356, lies at 6.7+-0.6 Mpc and spatially close to M106. We derive BR surface brightness profiles and colour gradients for all dwarfs and determine photometric and Sersic parameters. Finally, we discuss two objects in DDO71 and UGC7356 which may resemble GCs.