cataclysmic variables ( cvs ) are a diverse group of short - period binaries in which a late - type , roche - lobe filling secondary star transfers matter through an accretion disk onto a white dwarf primary .
the standard paradigm for the formation of cvs postulates that they start out as wide binaries ( on the order of 1 au ) , of moderate orbital period ( @xmath5 1 yr ) . as the more massive component ,
the white dwarf progenitor , evolves off of the main sequence , the secondary star suddenly finds itself orbiting within the red giant photosphere . during this `` common envelope '' ( ce ) phase ,
much of the angular momentum of the binary is believed to be shed due to interactions of the secondary star with the atmosphere of the red giant . during this process
, the binary period gets shortened , and the red giant envelope gets ejected . depending on the input parameters used in standard ce models ,
between 50 and 90% of all of these close binary stars merge ( c.f . , politano & weiler 2007 ) .
after emerging from the ce phase , an epoch of angular momentum loss via magnetic braking ( verbunt & zwaan 1981 ) , lasting @xmath6 10@xmath7 yr ( warner 1995 ) , is then required that leads to the formation of a semi - detached binary , and the mass transfer phase that signals the birth of a cv .    in this scenario , the vast majority of cv secondary stars should not show any signs of evolution , as there just has not been sufficient time for the low mass secondary stars ( m@xmath8 @xmath9 1 m@xmath10 ) in these systems to have begun to evolve off of the main sequence .
thus , the composition of the secondary stars should be relatively normal unless they have accreted matter with an unusual composition during the ce phase , or perhaps during classical novae eruptions . in an infrared spectroscopic survey of cvs , harrison et al .
( 2004 , 2005 ) found that the co features of the secondary stars in a number of cvs were much weaker than expected for their spectral types .
for u gem , the co features were extremely weak , yet the water vapor feature was normal for a spectral type of m4v .
this suggested a deficit of carbon .
hamilton et al .
( 2011 ) have tabulated the list of cvs with weaker than expected co features , and show that , when such data exists , analysis of uv spectra _ also _ find evidence for carbon deficiencies .
one issue with the infrared spectroscopy surveys was that they were qualitative , and were derived from lower resolution ( r @xmath0 2000 ) spectroscopy .
the line features in those spectra _
might _ be compromised due to veiling by a continuum source that arises from the accretion disk .
we have obtained moderate , and lower resolution , @xmath1-band spectroscopy of three of the brightest cvs : ss cyg , ru peg , and gk per .
previous near - infrared spectroscopy with spex ( harrison et al . 2004 , 2007 ) revealed that the secondary stars in these three systems had weaker than expected co features . in the following ,
we use a modified version of the moog spectral synthesis program to model moderate resolution @xmath1-band spectra to examine the atomic and molecular abundances in their secondary stars . as shown in hamilton et al .
( 2011 ) , higher resolution spectra dramatically reduce the effects of accretion disk veiling , allowing us to see secondary stars even in cv systems where the accretion luminosity totally dominates that of the secondary star . in the systems discussed here , however , the secondary star dominates the accretion disk flux in the @xmath1-band , and any veiling has to be quite small .
we then compare the best fit model spectra derived from the higher resolution data to lower resolution data .    in the next section
we discuss the data to be modeled , in section 3 we discuss the spectral synthesis procedure , in section 4 we present the modeling of the data , discuss the results in section 5 and draw our conclusions in section 6 .
ss cyg , ru peg , and gk per were observed using nirspec in high resolution mode on keck ii on 2004 august 28 .
the 0.432 x 12  slit was used providing a resolving power of r @xmath0 25000 , or a velocity resolution of @xmath11 @xmath0 12 km s@xmath12 .
we employed the two - nod script , and we used five minute exposure times for the three program cvs . to completely cover the @xmath1-band using nirspec requires two grating settings . due to time limitations , however , we only obtained data at a single grating position chosen to insure coverage of the first overtone bandhead of @xmath13co@xmath14 at 2.294 @xmath15 m .
this grating setting also covers the @xmath13co@xmath16 bandhead at 2.383 @xmath15 m .
in addition to the three program cvs , we also observed two late - type stars , @xmath17 eri ( k2v ) and @xmath18 psc ( g9iii ) , to provide templates for our synthetic spectral program given that the properties of those two objects are well known .
to correct for telluric absorption , we observed bright a0v stars located close to the program objects so as to minimize their relative differences in airmass .
these data were reduced using the idl routine `` @xmath19 '' , specially developed for nirspec package can be found here : http://www2.keck.hawaii.edu/inst/nirspec/redspec/index.html ] . in the @xmath1-band ,
the spectra of a0v stars are nearly featureless , except for the prominent h i brackett @xmath18 feature at 2.16 @xmath15 m .
the @xmath19 package does not attempt to correct for this feature , but can interpolate across such lines to reduce their impact upon division into the program star spectrum .
note that there is a weak telluric feature located very close to the brackett @xmath18 line , and thus the h i line profiles in spectra produced by the division of a `` patched '' a - star spectrum are slightly compromised .    just prior to our observing run on keck
, we observed the three program cvs using spex at the irtf on 2004 august 15 .
spex was used in cross - dispersed mode generating spectra spanning the range from 0.8 to 2.45 @xmath15 m . in the @xmath1-band ,
the dispersion is 5.4  / pix .
the raw data were reduced using spextool ( cushing et al . 2004 ) .
a log of all of our observations is presented in table 1 , and includes the orbital phase at the mid - point of the observational sequence for each of the three cvs .
the aavso light curve data base shows that ss cyg was in quiescence at the time of both the spex and nirspec observations , having reached visual maximum on 17 july .
ru peg was in outburst on 23 july , but was in quiescence for both observing runs . at the time of the infrared observations , gk per had @xmath20 = 13.1 , its normal minimum light value .
to calculate a synthetic spectrum , three basic building blocks are required : a stellar atmosphere model , a list of atomic and molecular lines , and a tool to calculate the corresponding flux at the stellar surface as a function of wavelength .
a full description of the modified version of moog and the grid generation and spectral matching python code ( called `` _ kmoog _ '' ) used below is described in hamilton ( 2013 ) , here we highlight the most important aspects .
there are several groups that provide precomputed atmosphere models , each with their own slightly different assumptions , simplifications , and basic input physics .
most of the models in general circulation , however , have the same underlying assumptions : the atmosphere is either modeled as a series of plane - parallel slabs ( 1d ) , or as a series of concentric shells ( 2d ) . the 1d and 2d atmosphere models
are calculated assuming local thermodynamic equilibrium ( lte ) conditions , in which one temperature can describe the particle velocity distribution as well as the population of atomic levels and ionization states around a given point in the stellar atmosphere .
it is also assumed that the atmosphere is time invariant , with the effects of dynamical processes such as mixing and convection parameterized .
bonifacio et al .
( 2012 ) provide an overview of the more detailed assumptions and differences between the available atmosphere computation methods and programs .
full 3d ( radiation - hydrodynamic ) stellar atmosphere models have been calculated across parameter grids ( e.g. , ludwig et al .
2009 , hauschildt and baron 2010 , magic et al .
2013a , 2013b ) . for this program ,
however , we have only used the 1d and 2d models .
the most commonly used atmosphere model codes are phoenix ( hauschildt et al .
( 1999 ) , allard et al .
2011 ) , marcs ( gustafsson et al . 1975 , 2008 ) , and atlas ( sbordone et al . 2004 ) .
phoenix is perhaps the most computationally sophisticated , sampling 42 million atomic and 550 million molecular lines when calculating opacity .
marcs model atmospheres are consistently updated through a web interface providing both 1d and 2d models . ultimately for this project
, we chose to use marcs atmospheres .
an important consideration for this choice was the existence of the interpolation code to create subgrid models , making it much simpler to construct a large grid of models quickly .
very recently , mszros et al .
( 2012 ) presented a new grid of atmospheres based off of atlas and marcs codes with an updated h@xmath21o line list as well as finer grid spacings in temperature , gravity , metallicity , and c and @xmath22 abundances .
given their recent availability , however , all analyses and spectra computed here use solely the original marcs grid .
the next crucial step to constructing a synthetic spectrum is to assemble a list of atomic and molecular transitions for which you wish to calculate an opacity at the stellar surface .
for our project we used the atlas line list with the addition of the co line list by goorvitch ( 1994 ) .
each line list was converted into a common format before combining , and , in the case of the co line list , only lines with @xmath23 and 2.28 @xmath24 2.405 @xmath15 m were kept .
the final list contains a total of 11,681 lines . for the purposes of this work , only computations involving diatomic molecules
were considered .
this eliminates lines of the important molecular species for the cooler m dwarfs commonly found in cvs , such as h@xmath21o . as described below , synthetic spectra for the sun and arcturus
were generated and compared to observations to provide the initial validation of this line list .
there are a number of codes that can readily make use of marcs and kurucz atmospheres .
moogchris / moog.html ] ( sneden 1973 ) has one of the longest track records , and given its ubiquity , we chose it to compute the synthetic spectra for this work . when creating a spectrum , moog
is internally limited to considering only 2500 lines in the line list . given the size of our line list , we modified moog ( see hamilton 2013 ) to allow for lists of 250,000 lines .
two additional changes were made to the moog source code .
first , a version was created to produce absolutely no output to the screen .
this prevents unnecessary lag and slowdowns , as well as greatly enhances usability when running a large number of models simultaneously .
second , and most importantly , a new binary output file type of moog was created .
significant savings in storage size were realized by creating a true fits output type , reducing the file size of a successfully created model by a factor of five or more .
the fits files created also have the benefit of being fully compliant and compatible with the standard cfitsio libraries , unlike the fits - like output already in the code .
the fits output created also has the benefit of storing both a full resolution spectrum as well as one convolved to lower resolutions , useful for quickly comparing against observations from different telescopes and instruments .      to validate our line list , we constructed synthetic spectra for the sun ( using the solar marcs model ) and arcturus ( using the parameters listed in table 2 , but with log@xmath25 = 3.0 ) to compare them to high resolution @xmath1-band spectra published by wallace et al .
( 1996 ) , and wallace & hinkle ( 1996 ) .
we concentrated on the reproduction of only the strongest spectral features observed in these two objects .
note that the arcturus comparison was only used to insure that we were not missing lines that become important at lower temperatures .
we then altered the oscillator strengths of the most highly discrepant lines to insure that they would both fit the observations better , and not bias the automated fitting process .
for example , there were a large number of sc i transitions that had oscillator strengths that were much too large .
unfortunately , oscillator strengths for many of the transitions in the @xmath1-band do not appear in the nist atomic spectra database , and thus we used trial and error adjustment of the tabulated oscillator strength in the line list to better reproduce the observed spectral feature ( for a primer on the identification of the strongest spectral features found in the @xmath1-band for k - type stars , see fig . 5 in harrison et al .
2004 ) .    to test our ability to reproduce observations ,
we have constructed synthetic spectra to compare to the nirspec observations of @xmath17 eri and @xmath18 psc , and to a sample lower resolution spectra of k dwarfs contained in the irtf spectral library ( cushing et al .
we summarize the spectroscopically derived parameters for these objects in table 2 , averaging the data compiled in the pastel catalog ( soubiran et al .
2010 ) for each object , constrained to measures made after 1970 .
the final column of table 2 denotes how many unique measurements went into the determination of each of the average properties for each star ( note that hd45977 only has a single measurement of its characteristics by sousa et al .
as we will see , this sample more than spans the range in metallicity observed for our program cvs .
note that the solar abundance pattern derived by grevesse et al .
( 2007 ) was used for all of the synthetic spectra , and only the abundance of carbon was altered for the investigation of the cvs .    in preparation for comparison to the models ,
the spectra were continuum normalized using a third order polynomial fit in iraf .
for the high resolution spectra , the process was quite trivial . for the full @xmath1-band spectra , however , we used an identical set of wavelength bins for identifying the continuum to be fit .
the continua were then fitted over the range 2.08 to 2.29 @xmath15 m , with bins that avoided all of the strongest spectral features .
the nirspec spectrum of @xmath17 eri is presented in fig . [ eps ] , where data for five of the seven spectral orders are shown ( the first order , spanning @xmath262.43 to 2.47 @xmath15 m , was too far red to obtain appreciable signal in our cv spectra , while the seventh order spectrum is dominated by the strong telluric feature at 2.05 @xmath15 m ) . as shown in table 2 , @xmath17 eri has a temperature near 5100 k , and a near - solar metallicity of [ fe / h ] = @xmath20.1 . in fig .
[ eps ] , we overplot a synthetic spectrum with t = 5100 k , and [ fe / h ] = @xmath20.125 ( interpolated between the model atmospheres with [ fe / h ] = @xmath20.25 , and [ fe / h ] = 0.0 ) .
the model fit to the strongest lines is excellent .
note that we have inserted an h i br@xmath18 line into our line list .
we do this to primarily show where this line is located .
the profiles of this feature as found in the infrared spectra for each of the objects are unreliable due to the reduction technique .
in addition , the profile of such a strong line can not be easily reproduced using moog .
@xmath18 psc allows us to test our ability to reproduce a much lower gravity object with a very low metallicity .
we used a marcs spherical , alpha - enhanced model atmosphere with log@xmath25 = 2.5 , and [ fe / h ] = @xmath20.5 to generate a synthetic spectrum for @xmath18 psc .
as noted in alves - brito et al .
( 2010 ) , @xmath18 psc has a modest alpha - enhancement of @xmath0 @xmath270.2 , and thus many of the metal lines are stronger in @xmath18 psc than they would be in a model that has a solar abundance pattern . since @xmath28 does not yet have the capabilities to adjust the abundances of any other species besides carbon , we ran the moog binary used by @xmath28 in a stand - alone mode , increasing the abundances of o , al and mg as detailed in alves - brito et al .
we overplot this model on the observed data in fig .
[ gamma ] .
the resulting synthetic spectrum is a reasonable match to the data .
given that the majority of cvs are too faint for moderate resolution infrared spectroscopy , most of the existing @xmath1-band data have been obtained at a lower resolution . to examine
how well we can model such data sets we have downloaded spex data for dwarf stars spanning the spectral types k0v to k7v .
the spectral type range k0v to k3v is shown in fig .
[ k0tok3 ] , while the range k4v to k7v is shown in fig .
[ k4tok7 ] .
we generated synthetic spectra for each of these stars using the parameters listed in table 2 .
the published metallicities produced excellent fits over most of the spectral range for all of the stars except hd145675 ( k0v ) .
we found that the best fitting model for hd145675 over the entire @xmath1-band required a lower metallicity , [ fe / h ] = @xmath270.25 , than tabulated for this object ( though within the observed error bars ) .
close inspection of fig .
[ k0tok3 ] shows that the synthetic spectra do not do a very good job of matching the strongest spectral features in the 2.10 to 2.15 @xmath15 m bandpass .
the origin of this discrepancy is not obvious .
the fit of the synthetic spectra over the rest of the @xmath1-band suggests the tabulated ( or derived ) metallicities are correct .
this might be a problem with continuum subtraction . to investigate this , we reduced and plot the spex data for the g9v star hd161198 observed during the same run as our program cvs [ note that hd161198 is a spectroscopic binary with a cool secondary star ( duquennoy et al .
1996 ) , and it is not a useful spectral template ] .
the model for hd145675 does a good job at fitting the spectrum of hd161198 in this range , indicating that continuum subtraction is not the issue .
this suggests that the problem could be due to either the data reduction process , or some issue with the telluric standard and/or object spectra .
the maximum amount of observed signal occurs in this bandpass , and the k0v , k1v , and k3v have @xmath1 @xmath9 4.7 ( hd3765 has @xmath1 = 5.16 , while hd145675 has @xmath1 = 5.55 ) , so perhaps it is a linearity issue .
the discrepancy between the models and observed spectra in the 2.10 to 2.15 @xmath15 m bandpass is not as severe for the k4v to k7v templates ( even though the k7v template is 61 cyg b , with @xmath1 = 2.32 ) .
clearly , however , our synthetic spectra do an excellent job at matching objects with known temperatures and metallicities . especially important for the current program
is the match of the model spectra to the observed co features over a wide range of temperature and metallicity .      in the preceding we have modeled the spectra of objects with known parameters . for modeling the cvs , where the temperatures and metallicities are poorly known , requires the generation of a large number of models , and a robust method to match the synthetic spectra to the observations . with the python wrapper
we have constructed to the modified version of the moog program described above , we are able to generate a large grid of models covering a user - defined range of the following parameters : t@xmath29 , log @xmath25 , [ fe / h ] , and [ c / fe ] .
microturbulence was fixed at 2 km s@xmath30 , and macroturbulence was not included since it is a small effect at the resolution of our observational data . in _
, if the synthetic spectra are designated to be stored as fits files , each model actually contains two synthetic spectra : one unsmoothed , full resolution model , and one model convolved with a rotationally broadened spectrum appropriate for the target at the observed spectral resolution ( which must be specified before the grid creation ) . a full grid covering the entire parameter space ( 3800 @xmath9 t@xmath29 @xmath9 5200 k , 3.0 @xmath9 log @xmath25 @xmath9 5.0
, @xmath31 @xmath9 [ fe / h ] @xmath9 @xmath271 , and @xmath31 @xmath9 [ c / fe ] @xmath9 @xmath271 ) requires about one day to create on a 32 core machine and occupies @xmath0 86 gbyte in fits output format .
obviously , smaller grids can be generated much more rapidly if the input parameter space can be narrowed .      with such large grids ,
it is easy to be overwhelmed by the quantity of data . given a suitable subset of parameters , a direct grid search is possible ; gather a list of valid models , and step through them one at a time comparing the model to observations .
given the sometimes low s / n of the observational data , however , a @xmath32 statistic is not calculated across the entire observational range , but restricted to the areas of greatest impact .
for example , including na i , ca i , mg i , al i , fe i lines , and the co bands .
this restricted region has the advantage of avoiding areas prone to increased levels of noise , and generally improves the quality of the fits to the observations .
another method of stepping through the parameter space involves the use of an evolutionary algorithm to find the best fit to the data .
evolutionary algorithms are based off the ideas of biological evolution , the simple premise that in some given population only the `` fittest individuals '' ( combination of all parameters ) will `` breed '' and influence the next generation of that population . if we imagine representing our parameters as 1d binary strings , then the behavior and the manipulation of the binary strings over the course of the fit is not unlike that of dna during reproduction .
sets of parameters can be mutated , introducing diversity into the population which helps to avoid getting stuck in a local maximum or minimum in the `` @xmath33 landscape '' .
mutations take the form of bit flips or swaps , analogous to transcription errors in dna replication .
we used a genetic algorithm framework called @xmath34 .
the algorithm proceeds as follows :    1 .   make an initial first guess at the parameters based on spectral shape 2 .
choose an appropriate maximum number of generations to evolve through as well as the initial population size 3 .
`` breed '' and `` mutate '' the `` fittest '' members in the population to create a new generation 4 .
repeat # 3 until the population converges to the desired level    in this manner , all parameters are varied continually as the population converges towards some minimum @xmath32 , at least a factor of four faster than the plain grid searching technique . using both grid search and genetic algorithm techniques ,
the grid created with moog can be readily compared with the appropriate subset of cv spectra to attempt to determine as many stellar parameters as possible .
both methods are available in _ kmoog _ , though we find the brute force is sufficiently rapid and as reliable as the genetic algorithm technique .
heat maps showing the best fitting models for any two parameters are then generated to allow the determination of the best fitting model(s ) .    as we will find below , determining the exact value t@xmath29 or
[ fe / h ] for an unconstrained object from a @xmath1-band spectrum is difficult .
many of the relevant oscillator strengths for transitions in this bandpass are unknown , and we have had to estimate a number of them by comparison of models to observations .
in addition , our line list is incomplete , and thus weaker spectra features can contribute in a negative fashion to the assessment of the fit of a model to the data .
the issues of telluric correction and continuum subtraction , as noted above , can introduce additional uncertainties into the goal of finding the best fitting solution .
thus , the @xmath32 analysis just described must be done carefully , focussing on the strongest spectral features , insuring that there is not a spectroscopic anomaly in that feature which can be attributed to random noise variation , or to the telluric correction process .
it is extremely difficult to produce a fully automated best fit regimen that can blindly solve for the parameters of a previously unconstrained object in the @xmath1-band , even when focused on specific spectral lines .
currently , we manually examine the @xmath32 analysis for each line / region to insure that it is reliable , and that no unexpected events have occurred that affect the line / region under analysis .
one such anomaly will be encountered shortly when we discuss the al i line ( at 2.116 @xmath15 m ) in ss cyg .
it has an unusual strength that defies obvious explanation . @xmath35
analysis of this line results in a super - solar abundance . we continue to work on the matching aspects of _ kmoog _ , but in its current state , a more interactive ( `` hands - on '' ) approach is required to insure the most robust results .
given that we have some information on the spectral types of the secondary stars in our cvs , we decided to use smaller grid than can be generated with the default settings of _
kmoog_. for example , tests in hamilton ( 2013 ) show that even the moderate resolution @xmath1-band spectra are inadequate for assessing gravity .
given the large rotational broadening of our targets , it will be impossible to constrain the surface gravity for the program cvs .
thus we fixed log @xmath25 = 4.5 for modeling ss cyg and ru peg , and log @xmath25 = 4.0 for gk per ( see below ) , to reduce parameter space , and increase the efficiency of the grid searching process .
we postpone the discussion of error bar estimation to section 4.4 .
ss cyg is the prototype long period ( p@xmath36 = 6.603 hr ) dwarf nova with eruptions that recur on a monthly basis .
harrison et al . ( 2004 ) estimate the spectral type of the secondary as k4/k5 , consistent with other published values .
they noted , however , that the co features were very weak for this spectral classification .
bitner et al . (
2007 ) have used moderate resolution visual band spectroscopy to examine the masses of both components , and found that the secondary star is 10 to 50% larger than a star of its mass and temperature should be .
they suggest that the secondary star in ss cyg is evolved , and may have a core that is depleted in hydrogen .
they also measured a rotational velocity of @xmath37sin@xmath38 = 89 km s@xmath12 for the secondary star .
the nirspec spectrum of ss cyg is presented in fig . [ sscyg ] .
the large rotational broadening of the spectral features is apparent .
our observations spanned @xmath39 = 0.05 . at the observed phase ( using the k@xmath8 = 162.5 km s@xmath12 value derived by bitner et al . ) , the total change in radial velocity over our set of observations would be @xmath11 = 11 km s@xmath12 , resulting in minimal `` orbital smearing . ''
concentrating on the sixth order bandpass ( 2.095 to 2.128 @xmath15 m ) , we convolved solar metallicity models having t = 4700 k with rotational velocities ranging from 70 to 110 km s@xmath12 .
the lowest @xmath33 values were for @xmath37sin@xmath38 = 90 @xmath40 10 km s@xmath12 .
for the modeling below , we adopt the bitner et al .
value of @xmath37sin@xmath38 = 89 km s@xmath12 for ss cyg .    to attempt to reproduce the nirspec data
, we generated a grid of models spanning 4200 k @xmath9 t@xmath29 @xmath9 4900 k , @xmath41 @xmath9 [ fe / h ] @xmath9 @xmath270.50 , and @xmath42 @xmath9 [ c / fe ] @xmath9 @xmath270.25 .
the resulting heat map for the sixth order , based on @xmath33 fits to _ only _ the first two of the three strong lines in the wavelength interval 2.104 and 2.118 @xmath15 m is shown in fig .
[ sshm02 ] .
the best fit occurs for t = 4700 k , and [ fe / h ] = @xmath20.25 .
the fourth order contains a number of metal lines from ti i and fe i , though this spectral region is noisier than the data at shorter wavelengths .
the resulting heat map , fig . [ sshm04 ] , is consistent with the result for the sixth order , though the value of [ fe / h ] is less constrained .
the results for the third order containing the @xmath13co@xmath14 bandhead is shown in fig .
[ sshm05 ] , and suggests [ c / fe ] @xmath9 @xmath20.5 .
the best fit synthetic spectrum is overplotted on the data in fig .
[ sscyg ] has these values , and fits the data set quite well .
there is only one apparent anomaly in the nirspec spectrum of ss cyg , and that is the al i line at 2.116 @xmath15 m .
we overplot the derived best fit model spectrum on the irtf spex observation in fig .
[ ssirtf ] .
the fit is quite good , given the s / n of those data .
the main differences between the model and the low resolution spectrum is again the al i line 2.116 @xmath15 m , and an mg i feature at 2.282 @xmath15 m .
the al i line was perfectly reproduced for the solar spectrum and for @xmath17 eri , but does appear to be slightly under - predicted for the irtf spectral sequence ( though note our earlier discussion of the anomalous continuum in this bandpass ) .
the line at 2.109 @xmath15 m is also due to al i , thus the abundance of aluminum is not unusual . it is odd that this line is perfectly reproduced for our other cvs .
there is an he i line at this position that is often seen weakly in emission in the @xmath1-band spectra of other cvs with more prominent accretion disks , such as tt ari , see fig .
[ ttari ] .
this emission line might be the cause of the peaks in the continuum on either side of this line in fig .
[ ssirtf ] , but the effect of any such emission line would be to lessen the depth of the feature , not increase its strength .
the mg i line at 2.282 @xmath15 m is well reproduced by _ kmoog _ for all of the template star spectra .
this line is very weak in the spex spectrum of ss cyg .
since the mg i line at 2.106 @xmath15 m is at the appropriate strength in both the moderate and lower resolution data sets , we conclude that this is simply a s / n or telluric correction issue .
the nirspec and spex data were obtained at quite different orbital phases , 0.20 and 0.40 , respectively .
the @xmath1-band light curve of ss cyg presented in harrison et al .
( 2007a ) shows no evidence for irradiation at @xmath43 = 0.5 , and the @xmath1 magnitudes at these two phases in this light curve are nearly identical .
thus , the fact that the same synthetic spectrum works for both epochs is not surprising .
this does prove , however , that irradiation of the secondary star is not the source that weakens the co absorption features , as the amount of the irradiated hemisphere seen at the dramatically different orbital phases of our spectra has no influence on the strength of the co features .
ru peg is a u gem type dwarf nova of very long period ( p@xmath36 = 8.99 hr ) , whose binary parameters were first examined by kraft ( 1962 ) , observing doubled ca ii lines at some points in the orbital phase and a spectrum broadly consistent with a g8 iv to k0 v. friend et al .
( 1990 ) observed the @xmath268190  na i doublet to measure radial velocity motion of ru peg , finding that the secondary star is cooler and less massive than expected , with the spectrum being consistent with a k3 dwarf .
they also note that the secondary star must be somewhat evolved , as the temperature , and masses derived are inconsistent with a main sequence dwarf .
harrison et al . ( 2004 ) presented irtf observations of ru peg , showing that the nir spectrum was consistent with a spectral type of k2 , but that ru peg must contain a subgiant since it would be over - luminous for a dwarf at the distance given by its parallax .
dunford et al . (
2012 ) recently presented spectra of ru peg analyzed with roche tomograms , deriving a much cooler spectral type of k5 .
along with the estimation of a cooler spectral type than found previously , dunford et al . used moderate resolution visual band spectra to estimate @xmath37sin@xmath38 = 89 km s@xmath12 .
as for ss cyg , we convolved a solar abundance model with t@xmath29 = 4700 k to make our own estimate of the rotation velocity and found a slightly higher value of @xmath37sin@xmath38 = 95 @xmath40 10 km s@xmath12 .
we have used this value in generating the synthetic spectra that follow . using the same size grid as that for ss cyg , the @xmath44 analysis for the sixth order , shown in fig .
[ ruhm04 ] , suggests solar metallicities , and 4600 @xmath9 t@xmath29 @xmath9 4700 k. the results for the fourth order ( fig .
[ ruhm02 ] ) indicate slightly subsolar metallicities , and 4500 @xmath9 t@xmath29 @xmath9 4900 k. we conclude that the secondary star has t@xmath29 = 4700 k , and [ fe / h ] = 0.0 .
analysis of the co features in the third spectral order ( fig .
[ ruhm05 ] ) finds [ c / fe ] = @xmath20.75 .
a synthetic spectrum with these values is overplotted on the nirspec data in fig .
[ rupeg ] , and provides an excellent fit to the data in all of the spectral orders .
the h i emission from the disk is lower than seen in ss cyg , and the underlying h
i br@xmath18 absorption from the secondary star leaves its imprint on the emission line profile .
note that the long period and orbital phase at which the observations occurred make for minimal orbital smearing for the duration of our data taking : @xmath11 = 20 km s@xmath12 .
the final model derived for the nirspec spectrum is overplotted on the spex data in fig .
[ rupegirtf ] .
this model fits the lower resolution data very well , including all of the other strong @xmath13co bandheads in the red part of the @xmath1-band .
sion et al .
( 2004 ) found that the carbon abundance in the photosphere of the white dwarf was 10% of solar .
thus , the carbon deficit observed in the uv can be directly traced to the secondary star .
sion et al . also found that silicon suffered a similar deficit as that of carbon .
the most prominent si i lines in the wavelength range plotted in fig .
[ rupegirtf ] are at 2.136 and 2.189 @xmath15 m , and while these lines are weak in our data , and the s / n not very high , they seem to be well represented by our solar abundance model .
gk per is an old classical nova , having experienced an extremely luminous nova eruption in 1901 ( see harrison et al .
it has since exhibited periodic dwarf nova outbursts every 2 - 3 years since 1966 , brightening by @xmath0 3 mag during these events .
kraft ( 1964 ) detected the mass donor star in the optical , adopting a spectral type of k2 iv and noting that the spectrum was highly variable with an equivalent width of sr ii 4077   consistent with a giant at one point and a dwarf the next .
more recently , morales - reuda et al .
( 2012 ) analyzed spectra of gk per , finding that the donor star appeared as k1 iv .
their analysis consisted of an iterative template subtraction method , however , not direct detection of the secondary star .
harrison et al . (
2007b ) presented a @xmath1-band spectrum of the mass donor star in gk per , and found that it has extremely weak co features .    the radial velocity study conducted by morales - reuda et al .
found @xmath37sin@xmath38 = 61.5 @xmath40 11.8 km s@xmath12 , k@xmath8 = 120.5 @xmath40 0.7 km s@xmath12 , and a mass ratio of @xmath45 = 0.55 @xmath40 0.21 .
given the luminosity of its classical nova eruption suggests a massive white dwarf .
if this is true , the radial velocity solution implies a surface gravity for the secondary star of log @xmath25 @xmath9 4 .
as noted above , it is not possible to constrain the surface gravity of the program objects given the limitations of our data set .
thus , in contrast to ss cyg and ru peg , we set the surface gravity to log @xmath25 = 4.0 for our modeling of gk per ( though this produced identical results to models with log @xmath25 = 4.5 ) . as for the preceding two objects
, we analyzed the rotational velocity apparent in our nirspec data and find @xmath37sin@xmath38 = 55 @xmath40 10 km s@xmath12 , slightly lower than that derived by morales - rueda et al . , but agreeing within the error bars of the two measurements .
we use this value in all of our modeling .
the nirspec data is plotted in fig .
[ gkper ] .
@xmath35 analysis of the model grid found these best fit values : t@xmath29 = 5100 k , [ fe / h ] = @xmath20.125 , and [ c / fe ] = @xmath21.0 .
gk per is fainter than our other two targets so the co region is very noisy , and thus the measurement of [ c / fe ] for gk per is slightly less certain than found for the other two targets .
it is interesting to note that the h i absorption of the underlying secondary star is sufficient to create the appearance of no h i br@xmath18 emission at the center of this feature .
the spex spectrum of gk per , fig .
[ gkperirtf ] , is even noisier than the nirspec data .
clearly , however , the same model does an excellent job at fitting the observations .      as demonstrated above , we are easily able to reproduce the spectra of objects with known parameters .
our examination of k stars in the irtf spectral library shows that if we know the value of [ fe / h ] , we can derive values of t@xmath29 , or vice versa , that are correct to within the published error bars for that object . to examine
how well we can determine t@xmath29 or [ fe / h ] if we do not know either value , we present a @xmath35 analysis for the irtf template hd36003 .
the heat map for the spectral region 2.18 @xmath9 @xmath26 @xmath9 2.37 @xmath15 m of this star is shown in fig .
[ k4vhm ] . while the models with parameters similar to the published values for both t@xmath29 ( 4615 @xmath40 29 k ) and [ fe / h ] ( @xmath20.14 @xmath40 0.08 ) do have low values of @xmath35 , there is a range of solutions with quantitatively similar fits to the data . as demonstrated in the plots shown in figs .
[ k0tok3 ] and [ k4tok7 ] , for k dwarfs , the strong na i doublet ( at 2.209 @xmath15 m ) , the ca i triplet ( at 2.26 @xmath15 m ) , and the co features all slowly get stronger as we decrease temperature .
thus , models can mimic the observed spectrum by either an enhanced abundance of [ fe / h ] and a hotter temperature , or a reduced value of [ fe / h ] and a lower temperature .
our @xmath35 analysis over such a large span in wavelength is certainly also affected by missing transitions for weak lines , the s / n of the data , and for this lower resolution spectrum , the continuum subtraction process .
the red end of the @xmath1-band is filled with strong telluric features , and even for high s / n data , small discrepancies can appear due to a slight mis - match of the telluric standard with the program object .
this is especially clear for the co features , as some bandheads in a particular spectrum are well fitted by the models , while others are not .    given these results , if we are able to constrain the limits on the temperature of the object , such as `` the secondary star in ss cyg has a spectral type near k5 '' , and not expect super - solar metallicities , we find that temperatures derived from these data are good to @xmath40 250 k , and [ fe / h ] is good to @xmath40 0.25 .
moderate resolution data that covered either the na i doublet or ca i triplet would almost certainly lead to more precise values of both quantities .    putting error bars on our values for [ c / fe ]
is more difficult in that none of our templates are known to have unusual values for this parameter .
it is clear that our models reproduce the co features in the template spectra over a large range in metallicity .
thus , altering just the abundance of carbon should not introduce any significant issues , as it simply acts to reduce the strength of the co features to mimic those seen in lower metallicity stars .
how a reduced carbon abundance affects the input stellar atmosphere is difficult to assess given the unknown evolutionary state of the cv secondary stars .
the presence of more than one co bandhead in all of the spectra also allows us to have more confidence in a particular value of [ c / fe ] . for the higher
s / n spectra , we expect that the precision of the values for [ c / fe ] are similar in scale to [ fe / h ] : @xmath40 0.25 .
for gk per the data are quite poor compared to the other two cvs , and the derived value for [ c / fe ] is less certain .
it is obvious , however , that the value of [ c / fe ] is lower in gk per than in the other two sources .
this result will be amplified if the gravity of gk per is significantly lower than our input value of log @xmath25 = 4.0 .
we have obtained moderate resolution @xmath1-band spectroscopy of three cvs that were previously noted to have weaker than expected co features : ss cyg , ru peg , and gk per . because those earlier data were obtained at lower resolution , r @xmath0 2,000 , questions were raised about the possibility of `` veiling '' of the co features to make them appear weaker . while we show below that this argument has no validity if realistic veiling scenarios are examined , it is the increased dispersion of the spectra obtained with nirspec that easily argues against such a scenario .
the nirspec data have a dispersion of 0.32  / pix , nearly 1/17 that of the spex data .
thus , the flux per pixel from any veiling component would be diminished by this factor .
model spectra fitted to the two different sets of data would give different values for the strength of the co features if there was some type of veiling confined to the co region .    to provide quantitative estimates of the relative carbon abundances for cvs
we have altered moog to allow us to model the entire @xmath1-band , and developed a software package that constructs large grids of synthetic spectra , and then searches these grids to find the best fitting solution to the observed data .
we specifically allowed for the abundance of carbon to be a separate input under the assumption that the @xmath13c abundance has been changed in the initial branches of the cno cycle , as argued in harrison et al .
obviously , weak co features could be explained by a deficit of oxygen , but a process that could cause such a deficit is not obvious . with this software we derived values of t@xmath29 , [ fe / h ] , and [ c / fe ] from the nirspec spectra for the three program objects .
we list these results in table 3 . except for the deficits of carbon , no other abundance anomalies are apparent , and the secondary stars of all three objects have solar , or slightly sub - solar abundances .
the total @xmath13c abundance for ss cyg and ru peg are similar , @xmath0 25% of solar , while for gk per it is @xmath0 10% of solar . if these peculiar values are the result of the cno cycle ( c.f .
, marks & sarna 1998 ) , than the abundance of @xmath46c should be enhanced , allowing for the detection of the bandheads of @xmath46co . the two strongest such bandheads with @xmath47 2.4 @xmath15 m are @xmath46co@xmath14 at 2.345 @xmath15 m , and @xmath46co@xmath48 at 2.374 @xmath15 m .
the first of these is not covered by our nirspec data , but dhillon et al .
( 2002 ) present a @xmath1-band spectrum of ss cyg which appears to show a strong @xmath46co@xmath14 absorption feature .
they attribute this feature to metal line absorption in star spots as indicated by anomalously strong na i absorption features ( at 2.3355 and 2.3386 @xmath15 m ) in their spectra .
these two na i lines are not anomalous in our spex data , nor is there a @xmath46co@xmath14 feature of the strength observed by dhillon et al .
the second of the @xmath46co features is covered by our nirspec observations , but there is a deep telluric feature at this wavelength that makes good correction quite difficult , and results in a noisy portion of the second order spectrum . even with this issue , there is certainly no sign of a strong absorption feature at this location in either ss cyg or ru peg .
the spex data covers both features , and while there are hints of both of these bandheads in the data for ss cyg and ru peg , the spectra are simply too noisy to draw any firm conclusions .
new , higher s / n moderate resolution spectroscopy will be required to properly investigate the @xmath13c/@xmath46c ratio .
harrison et al . ( 2004 ) found that the @xmath1-band continua of many long period cvs were flatter than expected if the secondary star was the sole source of emission in that bandpass . to more
closely examine the slope of the continuum of ss cyg , we compare the full @xmath1-band spectrum of ss cyg with that of the irtf template star that has the most similar temperature : hd45977 . in fig .
[ sscygcomp ] , we subtract the spectrum of hd45977 from that of ss cyg , both spectra having been normalized at 2.2 @xmath15 m . the result is a residual spectrum that is fairly flat , with f@xmath49 @xmath50 @xmath51 .
given the uncertainties in the temperatures of the two objects , and the differing sources of the data , such a flat spectrum is consistent with free - free emission .
the continuum subtraction process step required to allow us to directly compare the spex data to the synthetic spectra easily removes a weak continuum source such as this , and this is why the results for the nirspec and spex are simultaneously consistent with the same model .
we can also attempt to model the contribution of line emission sources that might preferentially veil the co features .
there are only two possible sources that we can envision that might do this : the h i pfund continuum , and co emission . as shown in fig .
[ sscygcomp ] , the h i br@xmath18 emission line in ss cyg is quite strong , so the pfund series of hydrogen will also be present .
the pfund series limit is at @xmath262.279 @xmath15 m , and the transitions pf23 and higher , will fall in our modeling bandpass .
the emissivities of these lines are very small ; hummer & storey ( 1987 ) list them as being @xmath9 1.7% that of br@xmath18 ( for `` case b '' conditions ) . fitting the br@xmath18 line with a gaussian profile
gives fwhm = 73  .
we will assume all of the h i lines have this profile . in fig .
[ sscygpfcomp ] , we present spectra where we added a rotationally broadened spectrum of the pfund series ( up to pf100 ) to the k4v spectral template , and compare it to the spectrum of ss cyg .
given their low emissivities , the effect of the h i lines is not detectable .
we also artificially multiplied the emissivities of the pfund lines by factors 10 and 100 and added them to the k4v template . only at very large multiplicative factors
can their presence be detected , and even then , they add a smooth continuum to the co region , leaving the absorption features relatively undisturbed , but producing a highly distorted continuum .
even with this simplistic model , it is clear that h i pfund emission can not provide the veiling necessary to explain weak co features .
the detection of co emission in wz sge ( howell et al . 2004 ) has led to the suggestion that co emission is the source of the veiling needed to produce weaker than expected co features seen in the program cvs . since the co emission must come from the accretion disk , a similar result as that for the pfund continuum emission attains . as discussed in scoville et al .
( 1980 ) , such emission is expected to be confined to regions where t @xmath52 3000 k , and where the densities are high , @xmath53 @xmath54 10@xmath55 @xmath56 . since ss cyg was in a minimum during our observations ,
presumably these conditions exist somewhere within its accretion disk .
as shown in bik et al .
( 2006 ) , optically thin co emission features are essentially the inverse of the co absorption features seen in late type stars .
this similarity allows us to construct crude models to investigate what co emission might look like if its origin was in a cv accretion disk .
if we assign the co emission to the disk , then the broadening of any spectral feature located there will be much greater than that of the slower rotational broadening of the secondary star ( 89 km s@xmath12 in the case of ss cyg ) .
for example , if the accretion disk extends to 75% of r@xmath57 ( the roche lobe radius of the white dwarf ) , using the stellar parameters for ss cyg in bitner et al .
( 2007 ) , the keplerian velocity at this radius will be @xmath0 400 km s@xmath12 . in fig .
[ sscygcocomp ] , we add a rotationally broadened co spectrum to the best fitting synthetic spectrum for ss cyg using the value for the h i br@xmath18 line ( fwhm = 73  ) , and a co emission spectrum broadened by half that value . as could be expected ,
the result is a dramatic change in the profiles of the co features , with an emission peak preceding each bandhead , a subtle redward shift for the location of the strongest absorption , and a dramatic narrowing of the co bandhead profiles .
obviously , the data for none of our objects resembles this result .
more complicated scenarios for the emission profile could be constructed , but to truly diminish the co features without distorting the resulting spectrum _ requires _ the addition of co emission with a broadening profile that is very similar to that of the secondary star .
there is no other location in a cv system besides the secondary star photosphere where such low velocities occur .
we have modified the spectral synthesis program moog and have demonstrated that we are able to model moderate and lower resolution @xmath1-band spectra for stars hotter than t@xmath29 @xmath54 4000 k.
we then used this program to investigate such spectra for three cvs : ss cyg , ru peg , and gk per .
it is clear that all three cvs in this sample have subsolar carbon abundances . that synthetic spectra with identical parameters simultaneously model data with dramatically different dispersions for each of the objects , and realistic co veiling scenarios
fail to reproduce the observed co features , proves that the carbon deficits are real .
it is important to note that these results represent the first direct measurements of abundances in the photosphere of the mass donor stars in cvs . while our determinations of t@xmath58 , [ fe / h ] , and [ c / fe ] for these three objects are not overly precise , they have no parallel .
there have not been many predictions for the behavior of the abundances in the atmospheres of cv donor stars as a function of age / evolutionary state .
marks & sarna ( 1998 ) explored how the photospheric abundances of cno elements in cv secondaries might evolve with time , including the possibilities of pollution by the common envelope phase , or through the accretion of classical novae ejecta .
they found that the latter two scenarios were inadequate in producing observable changes , and concluded that only pre - contact evolution could significantly alter the cno abundances and isotopic ratios in cv secondary star photospheres .
schenker et al .
( 2002 ) examined how cvs might evolve from supersoft binaries and , if so , could show dramatic changes in cno species . in light of our results
, it would be useful to revisit this research with newer generations of binary star evolution codes , such as starevol ( stancliffe & eldridge 2009 ) .
it is also obvious that moderate resolution spectroscopy of the brightest cvs covering the _ entire _ @xmath1-band are required to achieve more robust measures of [ fe / h ] , as well as allow for the proper examination of the co bands .
if the deficits of carbon arise from the cno cycle , than the isotopic abundances of carbon will be altered .
moderate resolution @xmath1-band spectroscopy covering the main @xmath46co features will allow for the confirmation of this inference .
both teh and rth were partially supported by a grant from the nsf ( ast-1209451 ) .
we acknowledge with thanks the variable star observations from the aavso international database contributed by observers worldwide and used in this research .
allard , f. , homeier , d. , & freytag , b. 2011 , aspc , 448 , 91 + alves - brito , a. , melendez , j. , asplund , m. , ramirez , i. , et al .
2010 , a&a , 513 , 35 + baron , e. , chen , b. , & hauschildt , p. h. 2009 , in aip conf .
, vol . 1171 , american institute of physics conference series , ed . i. hubeny , j. m. stone , k. macgregor , & k. werner , 148160 + bitner , m. a. , robinson , e. l. , & behr , b. b. 2007 , apj , 662 , 564 + bonifacio , p. , caffau , e. , ludwig , h .-
, & steffen , m. 2012 , iaus 282 , 213 + cushing , m. c. , rayner , j. t. , & vacca , w. d. 2005 , apj , 623 , 1115 + cushing , m. , vacca , w. d. , & rayner , j. t. 2004 , pasp , 116 , 362 + dhillon , v. s. , littlefair , s. p. , marsh , t. r. , sarna , m. j. , & boakes , e. h. 2002 , a&a , 393 , 611 + dunford , a. , watson , c. a. , & smith , r. c. 2012 , mnras , 422 , 3444 + duquennoy , a. , tokovinin , a. a. , leinert , ch .
, glindemann , et al .
1996 , a&a , 314 , 846 + friend , m. t. , martin , j. s. , connon - smith , r. , & jones , d. h. p. 1990 , mnras , 246 , 654 + goorvitch , d. 1994 , apjs , 95 , 535 + gray , d. f. , the observation and analysis of stellar photospheres ( cambridge : cambridge univ . press ) , 374 + gray , r. o. , & corbally , c. j. 1994 , aj , 107 , 742 + grevesse , n. , asplund , m. , & sauval , a. j. 2007 , space sci .
, 130 , 105 + gustafsson , b. , bell , r. a. , eriksson , k. , & nordlund , a. 1975 , a&a , 42 , 407 + gustafsson , b. , edvardsson , b. , eriksson , k. , j@xmath59rgensen , u. g. , nordlund , a. , & plez , b. 2008 , a&a , 486 , 951 + hamilton , r. t. 2013 , phd thesis , new mexico state university + hamilton , r. t. , harrison , t. e. tappert , c. , & howell , s. b. 2011 , apj , 728 , 16 + harrison , t. e. , bornak , j. , mcarthur , b. e. , & benedict , g. f. 2013 , apj , 767 , 7 + harrison , t. e. , howell , s. b. , szkody , p. , & cordova , f. a. 2007a , aj , 133 , 162 + harrison , t. e. , campbell , r. k. , howell , s. b. , cordova , f. a. , & schwope , a. d. 2007b , apj , 656 , 444 + harrison , t. e. , osborne , h. l. , & howell , s. b. 2005 , aj , 129 , 2400 + harrison , t. e. , osborne , h. l. , & howell , s. b. 2004 , aj , 127 , 3493 + hauschildt , p. h. , & baron , e. 2010 , a&a , 509 , 36 + hauschildt , p. h. , allard , f. , ferguson , j. , baron , e. , & alexander , d. r. 1999 , apj , 525 , 871 + howell , s. b. , harrison , t. e. , & szkody , p. 2004 ,
apj , 602 , l49 + hummer , d. g. , & storey , p. j. 1987 , mnras , 224 , 801 + kraft , r. r. 1964 , apj , 139 , 457 + kraft , r. p. 1962 ,
apj , 135 , 408 + ludwig , h. -g . , caffau , e. , steffen , m. , freytag , b. , et al .
2009 , mmsai , 80 , 711 + magic , z. , collett , r. , hayek , w. , & asplund , m. 2013 , a&a , 560 , 8 + magic , z. , collett , r. , asplund , m. , trampedach , r. , et al .
2013 , a&a , 557 , 26 + marks , p. b. , & sarna , m. j. 1998 , mnras , 301 , 699 + martnez - arniz , r. , maldonado , j. , montes , d. , eiroa , c. , & montesinos , b. 2010 , a&a , 520 , 79 + mszros , sz . ,
allende prieto , c. , edvardsson , b. , castelli , f. , et al .
2012 , aj , 144 , 120 + morales - reuda , l. , still , m. d. , roche , p. , wood , j. h. , & lockley , j. j. 2002 , mnras , 329 , 597 + politano , m. & weiler , k. p. 2007
, apj , 665 , 663 + reddy , r. r. , & viswanath , r. 1990 , journal of astrophysics and astronomy , 11 , 67 + sbordone , l. , bonifacio , p. , castelli , f. , & kurucz , r. l. 2004 , memorie della societa astronomica italiana supplementi , 5 , 93 + scoville , n. z. , krotkov , r. , & wang , d. 1980 , apj , 240 , 929 + schenker , k. , king , a. r. , kolb , u. , wynn , g. a. , & zhang , z. 2002 , mnras , 337 , 1105 + sion , e. m. , cheng , f. , godon , p. , urban , j. a. , & szkody , p. 2004a
, aj , 128 , 1834 + sneden , c. a. 1973 , phd thesis , university of texas + soubiran , c. , le campion , j. -f .
, cayrel de strobel , g. , & caillo , a. 2010 , a&a , 515 , 111 + sousa , s. g. , santos , n. c. , israelian , g. , mayor , m. , et al . 2011 , a&a , 533 , 141 + stancliffe , r. j. , & eldridge , j. j. 2009 , mnras , 396 , 1699 + verbunt , f. & zwaan , c. 1981 , a&a , 100 , 7 + wallace , k. , & hinkle , k. 1996 , apjs , 107 , 312 + wallace , k. , livingston , w. , hinkle , k. , & bernath , p. 1996 ,
apjs , 106 , 165 + warner , b. 1995 , cataclysmic variable stars ( cambridge : cambridge univ . press , 443 +    lccccc ss cyg & nirspec & 2004 aug .
28 & 12:31:51 & 4 @xmath60 300 & 0.20 + ru peg & nirspec & `` & 13:17:12 & 4 @xmath60 300 & 0.87 + gk per & nirspec & '' & 14:20:24 & 6 @xmath60 300 & 0.46 + @xmath18 psc & nirspec & `` & 13:45:17 & 8 @xmath60 1 & + @xmath17 eri & nirspec & '' & 15:06:18 & 24 @xmath60 1 & + ss cyg & spex & 2004 aug .
15 & 08:55:17 & 8 @xmath60 180 & 0.40 + ru peg & spex & `` & 10:25:06 & 8 @xmath60 240 & 0.85 + gk per & spex & '' & 13:53:54 & 6 @xmath60 240 & 0.94 +    lccccc sun & g2v & 5777 & 4.44 & 0.0 & + arcturus & k0iii & 4313 @xmath40 84 & 1.65 @xmath40 0.30 & @xmath20.54 @xmath40 0.10 & 34/31/31 + @xmath17 eri & k2v & 5088 @xmath40 72 & 4.54 @xmath40 0.16 & @xmath20.09 @xmath40 0.09 & 26/24/24 + @xmath18 psc & g9iii & 4853 @xmath40 45 & 2.59 @xmath40 0.26 & @xmath20.46 @xmath40 0.12 & 11/8/8 + hd145675 & k0v & 5320 @xmath40 114 & 4.45 @xmath40 0.07 & @xmath270.41 @xmath40 0.12 & 18/15/15 + hd10476 & k1v & 5189 @xmath40 53 & 4.46 @xmath40 0.12 & @xmath20.06 @xmath40 0.06 & 16/10/10 + hd3765 & k2v & 5023 @xmath40 66 & 4.53 @xmath40 0.16 & @xmath270.05 @xmath40 0.10 & 7/4/4 + hd219134 & k3v & 4837 @xmath40 138 & 4.51 @xmath40 0.13 & @xmath270.05 @xmath40 0.10 & 15/10/10 + hd45977 & k4v & 4689 @xmath40 174 & 4.30
@xmath40 0.39 & @xmath270.03 @xmath40 0.18 & 1/1/1 + hd36003 & k5v & 4615 @xmath40 29 & 4.35 @xmath40 0.06 & @xmath20.14 @xmath40 0.08 & 3/2/2 + hd201092 & k7v & 3964 @xmath40 175 & 4.50 @xmath40 0.18 & @xmath20.26 @xmath40 0.27 & 10/5/5 +    lcccc ss cyg & 4700 @xmath40 250 & 4.50 & @xmath20.25 @xmath61 0.25 & @xmath20.50 @xmath40 0.25 + ru peg & 4700 @xmath40 250 & 4.50 & @xmath270.00 @xmath40 0.25 & @xmath20.75 @xmath61 0.25 + gk per & 5100 @xmath40 250 & 4.00 & @xmath20.125 @xmath40 0.25 & @xmath21.00 @xmath40 0.50 +     eri ( black solid line ) .
the synthetic model spectrum ( green , dashed line ) generated for this object is overplotted . for the purposes of this paper
, the bottom panel will be defined as the second order , while the top panel will be labeled the sixth order
. the synthetic spectrum has been rotationally broaded by 4 km s@xmath12 as measured by martinez - arnaiz et al .
the s / n ratio ranged from 214 in the sixth order , to 95 in the second order.,title="fig:",width=566 ]      m are not especially well matched by the models .
the source of this discrepancy is hard to identify . in the top panel we also plot in red the spectrum of a g9v star ( hd161198 ) that we observed with spex during the 2004 august observing run , and the model spectrum generated for the k0v template fits the spectrum of hd161198 quite well.,title="fig:",width=566 ]         heat map for the order 6 spectral region for the nirspec spectrum of ss cyg . in this , and all subsequent heat maps , blue denotes @xmath62 @xmath63 0.7 ,
green is 0.7 @xmath9 @xmath62 @xmath63 0.9 , yellow is 0.9 @xmath9 @xmath62 @xmath63 1.5 , magenta is 1.5 @xmath9 @xmath62 @xmath63 2.5 , and red is @xmath62 @xmath54 2.5 . , title="fig:",width=566 ]       heat map for the order 3 spectral region for the nirspec spectrum of ss cyg . for this analysis ,
a grid of models with the metallicity fixed to [ fe / h ] = @xmath20.25 was constructed , while the carbon abundance was altered over the range @xmath21.0 @xmath9 [ c / fe ] @xmath9 @xmath270.25 .
, title="fig:",width=566 ]               heat map for the model grid fits to the third order spectrum of ru peg . like that for ss cyg ,
a grid of single metallicity models ( here with [ fe / h ] = 0.0 ) were generated with only the carbon abundance altered for this analysis.,title="fig:",width=566 ]             analysis of the irtf spectral template hd36003 for a large grid of model spectra .
the mean of the published values for this star , see table 2 , are t@xmath29 = 4550 k , and [ fe / h ] = @xmath20.17 .
the grid has temperature intervals of 100 k , and [ fe / h ] stepped by @xmath400.125.,title="fig:",width=491 ]       line of ss cyg . the spectrum below that has had the emissivities of the pfund lines increased by a factor of 10 .
the bottom spectrum has the pfund line emissivity increased by a factor of 100.,title="fig:",width=491 ]     line , added to the best fitting ( solar abundance ) synthetic spectrum .
the bottom spectrum is the result of the identical process , except the velocity broadening was one half that of br @xmath18 profile .
the co emission component was normalized so that when it was added to the model spectrum , the depth of the minima in the co features of ss cyg were reproduced .
dashed vertical lines indicate the wavelengths of the @xmath13co@xmath65 ( left ) and @xmath13co@xmath66 ( right ) bandheads.,title="fig:",width=642 ]