Source author record

C. Ian Short

C. Ian Short appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

14works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

14 published item(s)

preprint2022arXiv

Chroma+GAS: An Expedited Solution for the Chemical Equilibrium for Cool Stellar Atmospheres

We describe a unique approach to economizing the solution to the general chemical equilibrium and equation-of-state problem for late-type stars, including diatomic and polyatomic molecules, that is fast, accurate, and suitable for responsive approximate data modelling applications, and to more intensive modelling approaches in which the calculation of the gas equilibrium must be expedited to allow other aspects to be treated more realistically. The method, based on a novel economization of the Newton's method of solution of the linearized Saha and conservation equations, has been implemented in Python and made available as a stand-alone package, GASPy, and has been integrated into the interactive Python atmosphere and spectrum modelling code ChromaStarPy. As a result, ChromaStarPy now computes the state of the gas, the number density of absorbers, and the surface flux spectrum, with consistent inclusion of 105 chemical species, including 34 diatomic, and 16 polyatomic, neutral molecules, as well as H$^-$ and H$_2^+$, as well as many neutral and ionized atomic species. The economized method converges very rapidly and greatly improves the code's relevance to late-type stellar and brown dwarf spectrum modelling. We provide a brief overview of the GAS methodology, and present some illustrative results for the chemical equilibrium and spectrum for an M-type bright giant and dwarf, and a comparison to results of the PHOENIX/PPRESS package. All codes are available from the OpenStars www site: www.ap.smu.ca/OpenStars.

preprint2020arXiv

In situ exo-planet transit lightcurve modelling with the Chroma+ suite

We have added to the Chroma+ suite of stellar atmosphere and spectrum modelling codes the ability to synthesize the exo-planet transit lightcurve for planets of arbitrary size up to 10% of the host stellar radius, and arbitrary planetary and stellar mass and orbital radius (thus determining orbital velocity) and arbitrary orbital inclination. The lightcurves are computed in situ, integrated with the radiative transfer solution for the radiation field emerging from the stellar surface, and there is no limb-darkening parameterization. The lightcurves are computed for the Johnson-Bessel photometric system UBVRIHJK. We describe our method of computing the transit path, and the reduction in flux caused by occultation, and compare our lightcurve to an analytic solution with a four-parameter limb-darkening parameterization for the case of an edge-on transit of the Sun by Earth. This capability has been added to all ports and variations, including the Python port, ChromaStarPy, and the version that interpolates among the fully line-blanketed ATLAS9 surface intensity distributions, ChromaStarAtlas. All codes may be accessed at www.ap.smu.ca/OpenStars and at GitHub (github.com/sevenian3).

preprint2016arXiv

GrayStarServer: Server-side spectrum synthesis with a browser-based client-side user interface

I present GrayStarServer (GSS), a stellar atmospheric modeling and spectrum synthesis code of pedagogical accuracy that is accessible in any web browser on commonplace computational devices and that runs on a time-scale of a few seconds. The addition of spectrum synthesis annotated with line identifications extends the functionality and pedagogical applicability of GSS beyond that of its predecessor, GrayStar3 (GS3). The spectrum synthesis is based on a line list acquired from the NIST atomic spectra database, and the GSS post-processing and user interface (UI) client allows the user to inspect the plain text ASCII version of the line list, as well as to apply macroscopic broadening. Unlike GS3, GSS carries out the physical modeling on the server side in Java, and communicates with the JavaScript and HTML client via an asynchronous HTTP request. I also describe other improvements beyond GS3 such as more realistic modeling physics and use of the HTML <canvas> element for higher quality plotting and rendering of results, and include a comparison to Phoenix modeling. I also present LineListServer, a Java code for converting custom ASCII line lists in NIST format to the byte data type file format required by GSS so that users can prepare their own custom line lists. I propose a standard for marking up and packaging model atmosphere and spectrum synthesis output for data transmission and storage that will facilitate a web-based approach to stellar atmospheric modeling and spectrum synthesis. I describe some pedagogical demonstrations and exercises enabled by easily accessible, on-demand, responsive spectrum synthesis. GSS may serve as a research support tool by providing quick spectroscopic reconnaissance. GSS may be found at www.ap.smu.ca/~ishort/OpenStars/.

preprint2015arXiv

grayStar3 - gray no more: More physical realism and a more intuitive interface - all still in a WWW browser

The goal of the openStar project is to turn any WWW browser, running on any platform, into a virtual star equipped with parameter knobs and instrumented with output displays that any user can experiment with using any device for which a browser is available. grayStar3 (gS3) is a major improvement upon GrayStar 2.0 (GS2), both in the physical realism of the modeling and the intuitiveness of the user interface. The code integrates scientific modeling in JavaScript with output visualization HTML. The user interface is adaptable so as to be appropriate for a large range of audiences from the high-school to the introductory graduate level. The modeling is physically based and all outputs are determined entirely and directly by the results of in situ modeling, giving the code significant generality and credibility for pedagogical applications. gS3 also models and displays the circumstellar habitable zone (CHZ) and allows the user to adjust the greenhouse effect and albedo of the planet. In its default mode the code is guaranteed to return a result within a few second of wall-clock time on any device. The more advanced user has the option of turning on more realistic physics modules that address more advanced topics in stellar astrophysics. gS3 is a public domain, open source project and the code is available from www.ap.smu.ca/~ishort/grayStar3/ and is on GitHub. gS3 effectively serves as a public library of generic JavaScript+HTML plotting routines that may be recycled by the community.

preprint2015arXiv

NLTE and LTE Lick indices for red giants from [M/H] 0.0 to -6.0 at SDSS and IDS spectral resolution

We investigate the dependence of the complete system of 22 Lick indices on overall metallicity scaled from solar abundances, [M/H], from the solar value, 0.0, down to the extremely-metal-poor (XMP) value of -6.0, for late-type giant stars (MK luminosity class III, log(g)=2.0) of MK spectral class late-K to late-F (3750 < Teff < 6500 K) of the type that are detected as "fossils" of early galaxy formation in the Galactic halo and in extra-galactic structures. Our investigation is based on synthetic index values, I, derived from atmospheric models and synthetic spectra computed with PHOENIX in LTE and Non-LTE (NLTE), where the synthetic spectra have been convolved to the spectral resolution, R, of both IDS and SDSS (and LAMOST) spectroscopy. We identify nine indices, that we designate "Lick-XMP", that remain both detectable and significantly [M/H]-dependent down to [M/H] values of at least ~-5.0, and down to [M/H] ~ -6.0 in five cases, while also remaining well-behaved . For these nine, we study the dependence of I on NLTE effects, and on spectral resolution. For our LTE I values for spectra of SDSS resolution, we present the fitted polynomial coefficients, C_n, from multi-variate linear regression for I with terms up to third order in the independent variable pairs (Teff, [M/H]), and (V-K, [M/H]), and compare them to the fitted C_n values of Worthey et al. (1994) at IDS spectral resolution.

preprint2014arXiv

GrayStar: A Web application for pedagogical stellar atmosphere and spectral line modelling and visualisation

GrayStar is a stellar atmospheric and spectral line modelling, post-processing, and visualisation code, suitable for classroom demonstrations and laboratory-style assignments, that has been developed in Java and deployed in JavaScript and HTML. The only software needed to compute models and post-processed observables, and to visualise the resulting atmospheric structure and observables, is a common Web browser. Therefore, the code will run on any common PC or related X86 (-64) computer of the type that typically serves classroom data projectors, is found in undergraduate computer laboratories, or that students themselves own, including those with highly portable form-factors such as net-books and tablets. The user requires no experience with compiling source code, reading data files, or using plotting packages. More advanced students can view the JavaScript source code using the developer tools provided by common Web browsers. The code is based on the approximate gray atmospheric solution and runs quickly enough on current common PCs to provide near-instantaneous results, allowing for real time exploration of parameter space. I describe the user interface and its inputs and outputs and suggest specific pedagogical applications and projects. Therefore, this paper may serve as a GrayStar user manual for both instructors and students. In an accompanying paper, I describe the computational strategy and methodology as necessitated by Java and JavaScript. I have made the application itself, and the HTML, CSS, JavaScript, and Java source files available to the community. The Web application and source files may be found at www.ap.smu.ca/~ishort/GrayStar.

preprint2014arXiv

GrayStar: A Web application for pedagogical stellar atmosphere and spectral line modelling and visualisation II: Methods

GrayStar is a stellar atmospheric and spectral line modelling, post-processing, and visualisation code, suitable for classroom demonstrations and laboratory-style assignments, that has been developed in Java and deployed in JavaScript and HTML. The only software needed to compute models and post-processed observables, and to visualise the resulting atmospheric structure and observables, is a common Web browser. Therefore, the code will run on any common PC or related X86 (-64) computer of the type that typically serves classroom data projectors, is found in undergraduate computer laboratories, or that students themselves own, including those with highly portable form-factors such as net-books and tablets. The user requires no experience with compiling source code, reading data files, or using plotting packages. More advanced students can view the JavaScript source code using the developer tools provided by common Web browsers. The code is based on the approximate gray atmospheric solution and runs quickly enough on current common PCs to provide near-instantaneous results, allowing for real time exploration of parameter space. I describe the computational strategy and methodology as necessitated by Java and JavaScript. In an accompanying paper, I describe the user interface and its inputs and outputs and suggest specific pedagogical applications and projects. I have made the application itself, and the HTML, CSS, JavaScript, and Java source files available to the community. The Web application and source files may be found at www.ap.smu.ca/~ishort/GrayStar.

preprint2014arXiv

NLTE 1.5D Modeling of Red Giant Stars

Spectra for 2D stars in the 1.5D approximation are created from synthetic spectra of 1D non-local thermodynamic equilibrium (NLTE) spherical model atmospheres produced by the PHOENIX code. The 1.5D stars have the spatially averaged Rayleigh-Jeans flux of a K3-4 III star, while varying the temperature difference between the two 1D component models ($ΔT_{\mathrm{1.5D}}$), and the relative surface area covered. Synthetic observable quantities from the 1.5D stars are fitted with quantities from NLTE and local thermodynamic equilibrium (LTE) 1D models to assess the errors in inferred $T_{\mathrm{eff}}$ values from assuming horizontal homogeneity and LTE. Five different quantities are fit to determine the $T_{\mathrm{eff}}$ of the 1.5D stars: UBVRI photometric colors, absolute surface flux SEDs, relative SEDs, continuum normalized spectra, and TiO band profiles. In all cases except the TiO band profiles, the inferred $T_{\mathrm{eff}}$ value increases with increasing $ΔT_{\mathrm{1.5D}}$. In all cases, the inferred $T_{\mathrm{eff}}$ value from fitting 1D LTE quantities is higher than from fitting 1D NLTE quantities and is approximately constant as a function of $ΔT_{\mathrm{1.5D}}$ within each case. The difference between LTE and NLTE for the TiO bands is caused indirectly by the NLTE temperature structure of the upper atmosphere, as the bands are computed in LTE. We conclude that the difference between $T_{\mathrm{eff}}$ values derived from NLTE and LTE modelling is relatively insensitive to the degree of the horizontal inhomogeneity of the star being modeled, and largely depends on the observable quantity being fit.

preprint2013arXiv

Modeling the near-UV band of GK stars, Paper III: Dependence on abundance pattern

We extend the grid of NLTE models presented in Paper II to explore variations in abundance pattern in two ways: 1) The adoption of the Asplund et al. (2009) (GASS10) abundances, 2) For stars of metallicity, [M/H], of -0.5, the adoption of a non-solar enhancement of alpha-elements by +0.3 dex. Moreover, our grid of synthetic spectral energy distributions (SEDs) is interpolated to a finer numerical resolution in both T_eff (Delta T_eff = 25 K) and log g (Delta log g = 0.25). We compare the values of T_eff and log g inferred from fitting LTE and Non-LTE SEDs to observed SEDs throughout the entire visible band, and in an ad hoc "blue" band. We compare our spectrophotometrically derived T_eff values to a variety of T_eff calibrations, including more empirical ones, drawn from the literature. For stars of solar metallicity, we find that the adoption of the GASS10 abundances lowers the inferred T_eff value by 25 - 50 K for late-type giants, and NLTE models computed with the GASS10 abundances give T_eff results that are marginally in better agreement with other T_eff calibrations. For stars of [M/H]=-0.5 there is marginal evidence that adoption of alpha-enhancement further lowers the derived T_eff value by 50 K. Stellar parameters inferred from fitting NLTE models to SEDs are more dependent than LTE models on the wavelength region being fitted, and we find that the effect depends on how heavily line blanketed the fitting region is, whether the fitting region is to the blue of the Wien peak of the star's SED, or both.

preprint2012arXiv

Matching the Spectral Energy Distribution and p Mode Oscillation Frequencies of the Rapidly Rotating Delta Scuti Star ? Ophiuchi with a 2D Rotating Stellar Model

Spectral energy distributions are computed using 2D rotating stellar models and NLTE plane parallel model atmospheres. A rotating, 2D stellar model has been found which matches the observed ultraviolet and visible spectrum of ? Oph. The SED match occurs for the interferometrically deduced surface shape and inclination, and is different from the SED produced by spherical models. The p mode oscillation frequencies in which the latitudinal variation is modelled by a linear combination of eight Legendre polynomials were computed for this model. The five highest and seven of the nine highest amplitude modes show agreement between computed axisymmetric, equatorially symmetric mode frequencies and the mode frequencies observed by MOST to within the observational error. Including nonaxisymmetric modes up through |m| = 2 and allowing the possibility that the eight lowest amplitude modes could be produced by modes which are not equatorially symmetric produces matches for 24 out of the 35 MOST modes to within the observational error and another eight modes to within twice the observational error. The remaining three observed modes can be fit within 4.2 times the observational error, but even these may be fit to within the observational error if the criteria for computed modes are expanded.

preprint2012arXiv

Modeling the near-UV band of GK stars, Paper II: NLTE models

We present a grid of atmospheric models and synthetic spectral energy distributions (SEDs) for late-type dwarfs and giants of solar and 1/3 solar metallicity with many opacity sources computed in self-consistent Non-Local Thermodynamic Equilibrium (NLTE), and compare them to the LTE grid of Short & Hauschildt (2010) (Paper I). We describe, for the first time, how the NLTE treatment affects the thermal equilibrium of the atmospheric structure (T(tau) relation) and the SED as a finely sampled function of Teff, log g, and [A/H] among solar metallicity and mildly metal poor red giants. We compare the computed SEDs to the library of observed spectrophotometry described in Paper I across the entire visible band, and in the blue and red regions of the spectrum separately. We find that for the giants of both metallicities, the NLTE models yield best fit Teff values that are ~30 to 90 K lower than those provided by LTE models, while providing greater consistency between \log g values, and, for Arcturus, Teff values, fitted separately to the blue and red spectral regions. There is marginal evidence that NLTE models give more consistent best fit Teff values between the red and blue bands for earlier spectral classes among the solar metallicity GK giants than they do for the later classes, but no model fits the blue band spectrum well for any class. For the two dwarf spectral classes that we are able to study, the effect of NLTE on derived parameters is less significant.

preprint2011arXiv

Spectral Type and Radial Velocity Variations in Three SRC Variables

SRC variables are M supergiants, precursors to Type II supernovae, that vary in brightness with moderately regular periods of order 100-1000 days. Although identified as pulsating stars that obey their own period-luminosity relation, few have been examined in enough detail to follow the temperature and spectral changes that they undergo during their long cycles. The present study examines such changes for several SRC variables revealed by CCD spectra obtained at the Dominion Astrophysical Observatory (DAO) during 2005-2009, as well as by archival spectra from the DAO (and elsewhere) for some stars from the 1960s to 1980s, and Cambridge radial velocity spectrometer measures for Betelgeuse. Described here is our classification procedure and information on the spectral type and radial velocity changes in three of the stars. The results provide insights into the pulsation mechanism in M supergiants.

preprint2010arXiv

Modeling the near-UV band of GK stars, Paper I: LTE models

We present a grid of LTE atmospheric models and synthetic spectra that cover the spectral class range from mid-G to mid-K, and luminosity classes from V to III, that is dense in Teff sampling (Delta Teff=62.5 K), for stars of solar metallicity and moderately metal poor scaled solar abundance ([A/H]=0.0 and -0.5). All models have been computed with two choices of atomic line list: a) the "big" line lists of Kurucz (1992) that best reproduce the broad-band solar blue and near UV flux level, and b) the "small" lists of Kurucz & Peytremann (1975) that provide the best fit to the high resolution solar blue and near-UV spectrum. We compare our model SEDs to a sample of stars carefully selected from the large catalog of uniformly re-calibrated spectrophotometry of Burnashev (1985) with the goal of determining how the quality of fit varies with stellar parameters, especially in the historically troublesome blue and near-UV bands. We confirm that our models computed with the "big" line list recover the derived Teff values of the PHOENIX NextGen grid, but find that the models computed with the "small" line list provide greater internal self-consistency among different spectral bands, and closer agreement with the empirical Teff scale of Ramirez & Melendez (2005), but not to the interferometrically derived Teff values of Baines et al. (2010). We find no evidence that the near UV band discrepancy between models and observations for Arcturus (alpha Boo) reported by Short & Hauschildt (2003 and 2009) is pervasive, and that Arcturus may be peculiar in this regard.

preprint2008arXiv

Non-LTE modeling of the near UV band of late-type stars

We investigate the ability of both LTE and Non-LTE models to fit the near UV band absolute flux distribution and individual spectral line profiles of three standard stars for which high quality spectrophotometry and high resolution spectroscopy are available: The Sun (G2 V), Arcturus (K2 III), and Procyon (F5 IV-V). We investigate 1) the effect of the choice of atomic line list on the ability of NLTE models to fit the near UV band flux level, 2) the amount of a hypothesized continuous thermal absorption extinction source required to allow NLTE models to fit the observations, and 3) the semi-empirical temperature structure required to fit the observations with NLTE models and standard continuous near UV extinction. We find that all models that are computed with high quality atomic line lists predict too much flux in the near UV band for Arcturus, but fit the warmer stars well. The variance among independent measurements of the solar irradiance in the near UV is sufficiently large that we cannot definitely conclude that models predict too much near UV flux, in contrast to other recent results. We surmise that the inadequacy of current atmospheric models of K giants in the near UV band is best addressed by hypothesizing that there is still missing continuous thermal extinction, and that the missing near UV extinction becomes more important with decreasing effective temperature for spectral classes later than early G, suggesting a molecular origin.