PHOEBE 2 is a Python package for modeling the observables of eclipsing star systems, but until now has focused entirely on the forward-model -- that is, generating a synthetic model given fixed values of a large number of parameters describing the system and the observations. The inverse problem, obtaining orbital and stellar parameters given observational data, is more complicated and computationally expensive as it requires generating a large set of forward-models to determine which set of parameters and uncertainties best represent the available observational data. The process of determining the best solution and also of obtaining reliable and robust uncertainties on those parameters often requires the use of multiple algorithms, including both optimizers and samplers. Furthermore, the forward-model of PHOEBE has been designed to be as physically robust as possible, but is computationally expensive compared to other codes. It is useful, therefore, to use whichever code is most efficient given the reasonable assumptions for a specific system, but learning the intricacies of multiple codes presents a barrier to doing this in practice. Here we present the 2.3 release of PHOEBE (publicly available from http://phoebe-project.org) which introduces a general framework for defining and handling distributions on parameters, and utilizing multiple different estimation, optimization, and sampling algorithms. The presented framework supports multiple forward-models, including the robust model built into PHOEBE itself.

The Kepler mission has provided a treasure trove of eclipsing binaries (EBs), observed at extremely high photometric precision, nearly continuously for several years. We are carrying out a survey of ~100 of these EBs to derive dynamical masses and radii with precisions of 3% or better. We use multiplexed near-infrared H band spectroscopy from the SDSS-III and -IV APOGEE instrument and optical spectroscopy from the HET High-Resolution Spectrograph to derive double-lined spectroscopic orbits and dynamical mass-ratios for the EB sample. This information is combined with Kepler photometry to derive orbital inclination, dynamical masses of the components, radii and temperatures. These measurements are directly applicable for benchmarking stellar models that are integrating the next generation of improvements, such as the magnetic suppression of convection efficiency, updated opacity tables, and fine-tuned equations of state. We selected our EB sample to include systems with low-mass ($M \lt 0.8 \; M_{\odot}$) primary or secondary components, many expected to populate the relatively sparse parameter space below $0.5 \; M_{\odot}$. In this paper, we describe our EB sample and the analysis techniques we are utilizing, and present masses and radii for two systems that inhabit particularly underpopulated regions of mass-radius-period space: KIC 2445134 and KIC 3003991. Our analysis of KIC 2445134 ($q=0.411 \pm 0.001$) yields masses and radii of $M_A = 1.29 \pm 0.03 \; M_{\odot}$, $M_B = 0.53 \pm 0.01 \; M_{\odot}$, $R_A = 1.42 \pm 0.01 \; R_{\odot}$, $R_B = 0.510 \pm 0.004 \; R_{\odot}$, and a temperature ratio of $T_B/T_A = 0.635 \pm 0.001$; our analysis of KIC 3003991 ($q=0.298 \pm 0.006$) yields $M_A = 0.74 \pm \; 0.04 M_{\odot}$, $M_B = 0.222 \pm \; 0.007 M_{\odot}$, $R_A = 0.84 \pm 0.01 \; R_{\odot}$, $R_B = 0.250 \pm 0.004 \; R_{\odot}$, and a temperature ratio of $T_B/T_A= 0.662 \pm 0.001$.

Tidal forces in eccentric binary stars known as heartbeat stars excite detectable oscillations that shed light on the processes of tidal synchronization and circularization. We examine the pulsation phases of tidally excited oscillations (TEOs) in heartbeat binary systems. The target list includes four published heartbeat binaries and four additional systems observed by {\it Kepler}. To the first order, the pulsation phases of TEOs can be explained by the geometric effect of the dominant $l=2$, $m=0$, or $\pm 2$ modes assuming pulsations are adiabatic. We found that this simple theoretical interpretation can account for more than half of the systems on the list, assuming their spin and orbit axes are aligned. We do find significant deviations from the adiabatic predictions for some other systems, especially for the misaligned binary KIC 8164262. The deviations can potentially help to probe the non-adiabaticity of pulsation modes as well as resonances in the tidal forcing.

We present the characterization of KIC 4142768, an eclipsing binary with two evolved A-type stars in an eccentric orbit with a period of 14 days. We measure the fundamental parameters of the two components ($M_1=2.05M_{\odot}, R_1=2.96R_{\odot}$ and $M_2=2.05M_{\odot}, R_2=2.51R_{\odot}$) by combining {\it Kepler} photometry and spectra from {\it Keck} HIRES. The measured surface rotation rates are only one-fifth of the pseudo-synchronous rate of the eccentric orbit. Fourier spectrum of the light curve reveals hybrid pulsations of $\delta$ Scuti and $\gamma$ Doradus type, with pulsation frequencies at about $15-18$ day$^{-1}$ for p modes and about $0.2-1.2$ day$^{-1}$ for low-frequency g modes. Some of the g modes are exact orbital harmonics and are likely tidally excited. Their pulsation amplitudes and phases both agree with predictions from the linear tidal theory for $l=2, m=2$ prograde modes. We examine the period spacing patterns in the free oscillating g modes and identify them mostly as prograde sectoral dipole modes. The unstable frequency range and frequency spacing of p modes and the inferred asymptotic g-mode period spacings both agree with the stellar model for the primary star evolved to a late stage of the main sequence. The inferred rotation rate of the convective core boundary is very slow, similar to the small surface rotation rate inferred from the spectroscopy. The measured surface and near-core rotation rates provide constraints for testing the mechanism of angular momentum transfer and tidal synchronization in evolved eccentric binary star systems.

Theoretical scaling relations for solar-like oscillators and red giants are widely used to estimate fundamental stellar parameters. The accuracy and precision of these relations have often been questioned in the literature, with studies often utilising binarity for model-independent validation. However, it has not been tested if the photometric effects of binarity introduce a systematic effect on the extraction of the seismic properties of the pulsating component(s). In this paper, we present an estimation of the impact of a contaminating photometric signal with a distinct background profile on the global asteroseismic parameter $\nu_{\text{max}}$ through the analysis of synthetic red-giant binary light curves. We generate the pulsational and granulation parameters for single red giants with different masses, radii and effective temperatures from theoretical scaling relations and use them to simulate single red-giant light curves with the characteristics of Kepler long-cadence photometric data. These are subsequently blended together according to their light ratio to generate binary red-giant light curves of various configurations. We then perform a differential analysis to characterise the systematic effects of binarity on the extraction of $\nu_{\text{max}}$. We find that the $\nu_{\text{max}}$ extraction for red-giant power spectra featuring overlapping power excesses is unreliable if unconstrained priors are used. Outside of this scenario, we obtain results that are nearly identical to single-star case. We conclude that i) the photometric effects of binarity on the extraction of $\nu_{\text{max}}$ are largely negligible as long as the power excesses of the individual components do not overlap, and that ii) there is minimal advantage to using more than two super-Lorentzian components to model the granulation signal of a binary red-giant.

Binary systems where the axis of rotation (spin) of one or both components is tilted w.r.t. the axis of revolution are called misaligned binary systems. The angle of misalignment, obliquity, has been measured for a handful of stars and extrasolar planets to date. Here we present a mathematical framework for a complete and rigorous treatment of misalignment and introduce an extension to the public PHOEBE code that implements this framework. We discuss misalignment for the Roche geometry and introduce methods for computing stellar shapes, equilibrium (generalized Lagrange) points of the potential and minimal requirements for lobe existence. Efficient parametrization of misalignment is proposed in the plane-of-sky coordinates and implementation details in PHOEBE are given alongside the proof-of-concept toy model, comparison with a known misaligned binary DI Her, and comparison with a misaligned planetary system Kepler-13. We provide important mathematical details of the model in the Appendix. This paper accompanies the release of PHOEBE 2.1, which will be available soon from its website http://phoebe-project.org.

The main aim of this experiment is to provide a complete census of old (t > 10 Gyr, RR Lyrae, type II Cepheids, red horizontal branch), intermediate age (red clump, Miras) and young (classical Cepheids) stellar tracers across the Galactic Bulge. To fully exploit the unique photometric quality of LSST images, we plan to perform a Shallow minisurvey (ugrizy, -20 < l < 20 deg, -15 < b < 10 deg) and a Deep minisurvey (izy, -20 < l < 20 deg, -3 < b < 3 deg). The former one is aimed at constraining the 3D structure of the galactic Bulge across the four quadrants, and in particular, the transition between inner and outer Bulge. The u,g,r,i,z,y LSST bands provide fundamental diagnostics to constrain the evolutionary properties of low and intermediate-mass stars when moving from a metal-poor to a metal-rich regime. The deep minisurvey is aimed at tracing RR Lyrae, Red Clump stars, Miras and classical Cepheids in highly reddened regions of the Galactic center. These images will allow us to investigate the role that baryonic mass and dark matter played in the early formation and evolution of the MW.

Regular sampling in the time domain results in aliasing in the frequency domain that complicates the accurate determination of the periods of astrophysical variables. We propose to actively break the regularity of this sampling by providing an additional consideration for the scheduler that weights fields according to when observations will contribute the least to aliasing. The current aliases for each field can be computed during daytime from the history of observations. We can then calculate the times when additional observations would worsen or alleviate these aliases for different fields. The scheduler should give preference to observation epochs that lessen the effect of aliasing, while still meeting all other cadence requirements.

Tidal dissipation in binary star and planetary systems is poorly understood. Fortunately, eccentric binaries known as heartbeat stars often exhibit tidally excited oscillations, providing observable diagnostics of tidal circularization mechanisms and timescales. We apply tidal theories to observations of the heartbeat star KIC 8164262, which contains an F-type primary in a very eccentric orbit that exhibits a prominent tidally excited oscillation. We demonstrate that the prominent oscillation is unlikely to result from a chance resonance between tidal forcing and a stellar oscillation mode. However, the oscillation has a frequency and amplitude consistent with the prediction of resonance locking, a mechanism in which coupled stellar and orbital evolution maintain a stable resonance between tidal forcing and a stellar oscillation mode. The resonantly excited mode produces efficient tidal dissipation (corresponding to an effective tidal quality factor $Q \sim 5 \times 10^4$), such that tidal orbital decay/circularization proceeds on a stellar evolution time scale.

We present the analysis of KIC 8164262, a heartbeat star with a high-amplitude (~1 mmag), tidally resonant pulsation (a mode in resonance with the orbit) at 229 times the orbital frequency and a plethora of tidally induced g-mode pulsations (modes excited by the orbit). The analysis combines Kepler light curves with follow-up spectroscopic data from the Keck telescope, KPNO (Kitt Peak National Observatory) 4-m Mayal telescope and the 2.7-m telescope at the McDonald observatory. We apply the binary modelling software, PHOEBE, to the Kepler light curve and radial velocity data to determine a detailed binary star model that includes the prominent pulsation and Doppler boosting, alongside the usual attributes of a binary star model (including tidal distortion and reflection). The results show that the system contains a slightly evolved F star with an M secondary companion in a highly eccentric orbit (e = 0.886). We use the results of the binary star model in a companion paper (Fuller et al., 2017) where we show that the prominent pulsation can be explained by a tidally excited oscillation mode held near resonance by a resonance locking mechanism.

Eccentric binaries known as heartbeat stars experience strong dynamical tides as the stars pass through periastron, providing a laboratory to study tidal interactions. We measure the rotation periods of 24 heartbeat systems, using the Kepler light curves to identify rotation peaks in the Fourier transform. Where possible, we compare the rotation period to the pseudosynchronization period derived by Hut 1981. Few of our heartbeat stars are pseudosynchronized with the orbital period. For four systems, we were able to identify two sets of rotation peaks, which we interpret as the rotation from both stars in the binary. The majority of the systems have a rotation period that is approximately 3/2 times the pseudosynchronization period predicted by Hut 1981, suggesting that other physical mechanisms influence the stars' rotation, or that stars typically reach tidal spin equilibrium at a rotation period slightly longer than predicted.

The precision of photometric and spectroscopic observations has been systematically improved in the last decade, mostly thanks to space-borne photometric missions and ground-based spectrographs dedicated to finding exoplanets. The field of eclipsing binary stars strongly benefited from this development. Eclipsing binaries serve as critical tools for determining fundamental stellar properties (masses, radii, temperatures and luminosities), yet the models are not capable of reproducing observed data well either because of the missing physics or because of insufficient precision. This led to a predicament where radiative and dynamical effects, insofar buried in noise, started showing up routinely in the data, but were not accounted for in the models. PHOEBE (PHysics Of Eclipsing BinariEs; http://phoebe-project.org) is an open source modeling code for computing theoretical light and radial velocity curves that addresses both problems by incorporating missing physics and by increasing the computational fidelity. In particular, we discuss triangulation as a superior surface discretization algorithm, meshing of rotating single stars, light time travel effect, advanced phase computation, volume conservation in eccentric orbits, and improved computation of local intensity across the stellar surfaces that includes photon-weighted mode, enhanced limb darkening treatment, better reflection treatment and Doppler boosting. Here we present the concepts on which PHOEBE is built on and proofs of concept that demonstrate the increased model fidelity.

Heartbeat stars (HB stars) are a class of eccentric binary stars with close periastron passages. The characteristic photometric HB signal evident in their light curves is produced by a combination of tidal distortion, heating, and Doppler boosting near orbital periastron. Many HB stars continue to oscillate after periastron and along the entire orbit, indicative of the tidal excitation of oscillation modes within one or both stars. These systems are among the most eccentric binaries known, and they constitute astrophysical laboratories for the study of tidal effects. We have undertaken a radial velocity (RV) monitoring campaign of Kepler HB stars in order to measure their orbits. We present our first results here, including a sample of 21 Kepler HB systems, where for 19 of them we obtained the Keplerian orbit and for 3 other systems we did not detect a statistically significant RV variability. Results presented here are based on 218 spectra obtained with the Keck/HIRES spectrograph during the 2015 Kepler observing season, and they have allowed us to obtain the largest sample of HB stars with orbits measured using a single instrument, which roughly doubles the number of HB stars with an RV measured orbit. The 19 systems measured here have orbital periods from 7 to 90 d and eccentricities from 0.2 to 0.9. We show that HB stars draw the upper envelope of the eccentricity - period distribution. Therefore, HB stars likely represent a population of stars currently undergoing high eccentricity migration via tidal orbital circularization, and they will allow for new tests of high eccentricity migration theories.

Heartbeat stars are eccentric (e > 0.2) ellipsoidal variables whose light curves resemble a cardiogram. We present the observations and corresponding model of KIC 3749404, a highly eccentric (e = 0.66), short period (P = 20.3 d) heartbeat star with tidally induced pulsations. A binary star model was created using PHOEBE, which we modified to include tidally induced pulsations and Doppler boosting. The morphology of the photometric periastron variation (heartbeat) depends strongly on the eccentricity, inclination and argument of periastron. We show that the inclusion of tidally induced pulsations in the model significantly changes the parameter values, specifically the inclination and those parameters dependent on it. Furthermore, we determine the rate of apsidal advance by modelling the periastron variation at the beginning and end of the 4-yr Kepler data set and dividing by the elapsed time. We compare the model with the theoretical expectations for classical and general relativistic apsidal motion and find the observed rate to be two orders of magnitude greater than the theoretical rate. We find that the observed rate cannot be explained by tidally induced pulsations alone and consequently hypothesise the presence of a third body in the system.

We explore the phenomenology of a class of models where the dark matter particle can inelastically up-scatter to a heavier excited state via off-diagonal dipolar interactions with the interstellar plasma (gas or free electrons). The heavier particle then rapidly decays back to the dark matter particle plus a quasi-monochromatic photon. For the process to occur at appreciable rates, the mass splitting between the heavier state and the dark matter must be comparable to, or smaller than, the kinetic energy of particles in the plasma. As a result, the predicted photon line falls in the soft X-ray range, or, potentially, at arbitrarily lower energies. We explore experimental constraints from cosmology and particle physics, and present accurate calculations of the dark matter thermal relic density and of the flux of monochromatic X-rays from thermal plasma excitation. We find that the model provides a natural explanation for the observed 3.5 keV line from clusters of galaxies and from the Galactic center, and is consistent with null detections of the line from dwarf galaxies. The unique line shape, which will be resolved by future observations with the Hitomi (formerly Astro-H) satellite, and the predicted unique morphology and target-temperature dependence will enable easy discrimination of this class of models versus other scenarios for the generation of the 3.5 keV line or of any other unidentified line across the electromagnetic spectrum.

The primary Kepler Mission provided nearly continuous monitoring of ~200,000 objects with unprecedented photometric precision. We present the final catalog of eclipsing binary systems within the 105 square degree Kepler field of view. This release incorporates the full extent of the data from the primary mission (Q0-Q17 Data Release). As a result, new systems have been added, additional false positives have been removed, ephemerides and principal parameters have been recomputed, classifications have been revised to rely on analytical models, and eclipse timing variations have been computed for each system. We identify several classes of systems including those that exhibit tertiary eclipse events, systems that show clear evidence of additional bodies, heartbeat systems, systems with changing eclipse depths, and systems exhibiting only one eclipse event over the duration of the mission. We have updated the period and galactic latitude distribution diagrams and included a catalog completeness evaluation. The total number of identified eclipsing and ellipsoidal binary systems in the Kepler field of view has increased to 2878, 1.3% of all observed Kepler targets. An online version of this catalog with downloadable content and visualization tools is maintained at http://keplerEBs.villanova.edu.

The unparalleled photometric data obtained by NASA's Kepler Space Telescope has led to improved understanding of red-giant stars and binary stars. We discuss the characterization of known eccentric system, containing a solar-like oscillating red-giant primary component. We also report several new binary systems that are candidates for hosting an oscillating companion. A powerful approach to study binary stars is to combine asteroseimic techniques with light curve fitting. Seismology allows us to deduce the properties of red giants. In addition, by modeling the ellipsoidal modulations we can constrain the parameters of the binary system. An valuable independent source are ground-bases, high-resolution spectrographs.

KIC 8569819 is an eclipsing binary star with an early F primary and G secondary in a 20.85-d eccentric orbit. The primary is a delta Sct - gamma Dor star pulsating in both p modes and g modes. Using 4 years of Kepler Mission photometric data we independently model the light curve using the traditional technique with the modelling code PHOEBE, and we study the orbital characteristics using the new frequency modulation (FM) technique. We show that both methods provide the equivalent orbital period, eccentricity and argument of periastron, thus illustrating and validating the FM technique. In the amplitude spectrum of the p-mode pulsations we also discovered an FM signal compatible with a third body in the system, a low mass M dwarf in an 861-d orbit around the primary pair. However, the eclipses show no timing variations, indicating that the FM signal is a consequence of the intrinsic change in pulsation frequency, thus providing a cautionary tale. Our analysis shows the potential of the FM technique using {\it Kepler} data, and we discuss the prospects to detect planets and brown dwarfs in Kepler data for A and F stars even in the absence of transits and with no spectroscopic radial velocity curves. This opens the possibility of finding planets orbiting hotter stars that cannot be found by traditional techniques.

KIC 8560861 (HD 183648) is a marginally eccentric (e=0.05) eclipsing binary with an orbital period of P_orb=31.973d, exhibiting mmag amplitude pulsations on time scales of a few days. We present the results of the complex analysis of high and medium-resolution spectroscopic data and Kepler Q0 -- Q16 long cadence photometry. The iterative combination of spectral disentangling, atmospheric analysis, radial velocity and eclipse timing variation studies, separation of pulsational features of the light curve, and binary light curve analysis led to the accurate determination of the fundamental stellar parameters. We found that the binary is composed of two main sequence stars with an age of 0.9\+-0.2 Gyr, having masses, radii and temperatures of M_1=1.93+-0.12 M_sun, R_1=3.30+-0.07 R_sun, T_eff1=7650+-100 K for the primary, and M_2=1.06+-0.08 M_sun, R_2=1.11+-0.03 R_sun, T_eff2=6450+-100 K for the secondary. After subtracting the binary model, we found three independent frequencies, two of which are separated by twice the orbital frequency. We also found an enigmatic half orbital period sinusoidal variation that we attribute to an anomalous ellipsoidal effect. Both of these observations indicate that tidal effects are strongly influencing the luminosity variations of HD 183648. The analysis of the eclipse timing variations revealed both a parabolic trend, and apsidal motion with a period of (P_apse)_obs=10,400+-3,000 y, which is three times faster than what is theoretically expected. These findings might indicate the presence of a distant, unseen companion.

The unparalleled photometric data obtained by NASA's Kepler space telescope led to an improved understanding of red giant stars and binary stars. Seismology allows us to constrain the properties of red giants. In addition to eclipsing binaries, eccentric non-eclipsing binaries, exhibiting ellipsoidal modulations, have been detected with Kepler. We aim to study the properties of eccentric binary systems containing a red giant star and derive the parameters of the primary giant component. We apply asteroseismic techniques to determine masses and radii of the primary component of each system. For a selected target, light and radial velocity curve modelling techniques are applied to extract the parameters of the system. The effects of stellar on the binary system are studied. The paper presents the asteroseismic analysis of 18 pulsating red giants in eccentric binary systems, for which masses and radii were constrained. The orbital periods of these systems range from 20 to 440days. From radial velocity measurements we find eccentricities between e=0.2 to 0.76. As a case study we present a detailed analysis of KIC5006817. From seismology we constrain the rotational period of the envelope to be at least 165 d, roughly twice the orbital period. The stellar core rotates 13 times faster than the surface. From the spectrum and radial velocities we expect that the Doppler beaming signal should have a maximum amplitude of 300ppm in the light curve. Through binary modelling, we determine the mass of the secondary component to be 0.29$\pm$0.03\,$M_\odot$. For KIC5006817 we exclude pseudo-synchronous rotation of the red giant with the orbit. The comparison of the results from seismology and modelling of the light curve shows a possible alignment of the rotational and orbital axis at the 2$\sigma$ level. Red giant eccentric systems could be progenitors of cataclysmic variables and hot subdwarf B stars.

We present Kepler photometry and ground based spectroscopy of KIC 4544587, a short-period eccentric eclipsing binary system with self-excited pressure and gravity modes, tidally excited modes, tidally influenced p modes, and rapid apsidal motion of 182 y per cycle. The primary and secondary components of KIC 4544587 reside within the delta Scuti and gamma Dor instability region of the Hurtzsprung-Russell diagram, respectively. By applying the binary modelling software PHOEBE to prewhitened Kepler photometric data and radial velocity data obtained using the William Herschel Telescope and 4-m Mayall telescope at KPNO, the fundamental parameters of this important system have been determined, including the stellar masses, 1.98+/-0.07 MSun and 1.60+/-0.06 MSun, and radii, 1.76+/-0.03 RSun and 1.42+/-0.02 RSun, for the primary and secondary components, respectively. Frequency analysis of the residual data revealed 31 modes, 14 in the gravity mode region and 17 in the pressure mode region. Of the 14 gravity modes 8 are orbital harmonics: a signature of tidal resonance. While the measured amplitude of these modes may be partially attributed to residual signal from binary model subtraction, we demonstrate through consideration of the folded light curve that these frequencies do in fact correspond to tidally excited pulsations. Furthermore, we present an echelle diagram of the pressure mode frequency region (modulo the orbital frequency) and demonstrate that the tides are also influencing the p modes. A first look at asteroseismology hints that the secondary component is responsible for the p modes, which is contrary to our expectation that the hotter star should pulsate in higher radial overtone, higher frequency p modes.

We have discovered a class of eccentric binary systems within the Kepler data archive that have dynamic tidal distortions and tidally-induced pulsations. Each has a uniquely shaped light curve that is characterized by periodic brightening or variability at time scales of 4-20 days, frequently accompanied by shorter period oscillations. We can explain the dominant features of the entire class with orbitally-varying tidal forces that occur in close, eccentric binary systems. The large variety of light curve shapes arises from viewing systems at different angles. This hypothesis is supported by spectroscopic radial velocity measurements for five systems, each showing evidence of being in an eccentric binary system. Prior to the discovery of these 17 new systems, only four stars, where KOI-54 is the best example, were known to have evidence of these dynamic tides and tidally-induced oscillations. We perform preliminary fits to the light curves and radial velocity data, present the overall properties of this class and discuss the work required to accurately model these systems.

Encoded within the morphological structure of galaxies are clues related to their formation and evolutionary history. Recent advances pertaining to the statistics of galaxy morphology include sophisticated measures of concentration (C), asymmetry (A), and clumpiness (S). In this study, these three parameters (CAS) have been applied to a suite of simulated galaxies and compared with observational results inferred from a sample of nearby galaxies. The simulations span a range of late-type systems, with masses between ~1e10 Msun and ~1e12 Msun, and employ star formation density thresholds between 0.1 cm^-3 and 100 cm^-3. We have found that the simulated galaxies possess comparable concentrations to their real counterparts. However, the results of the CAS analysis revealed that the simulated galaxies are generally more asymmetric, and that the range of clumpiness values extends beyond the range of those observed. Strong correlations were obtained between the three CAS parameters and colour (B-V), consistent with observed galaxies. Furthermore, the simulated galaxies possess strong links between their CAS parameters and Hubble type, mostly in-line with their real counterparts.