We have identified 230 Tycho-2 Spectral Catalog stars that exhibit 8 micron mid-infrared extraphotospheric excesses in the MidCourse Space Experiment (MSX) and Spitzer Space Telescope Galactic Legacy MidPlane Survey Extraordinaire (GLIMPSE) surveys. Of these, 183 are either OB stars earlier than B8 in which the excess plausibly arises from a thermal bremsstrahlung component or evolved stars in which the excess may be explained by an atmospheric dust component. The remaining 47 stars have spectral classifications B8 or later and appear to be main sequence or late pre-main-sequence objects harboring circumstellar disks. Six of the 47 stars exhibit multiple signatures characteristic of pre-main-sequence circumstellar disks, including emission lines, near-infrared K-band excesses, and X-ray emission. Approximately one-third of the remaining 41 sources have emission lines suggesting relative youth. Of the 25 GLIMPSE stars with SST data at >24 microns, 20 also show an excess at 24 microns. Three additional objects have 24 micron upper limits consistent with possible excesses, and two objects have photospheric measurements at 24 microns. Six MSX sources had a measurement at wavelengths >8 microns. We modeled the excesses in 26 stars having two or more measurements in excess of the expected photospheres as single-component blackbodies. We determine probable disk temperatures and fractional infrared luminosities in the range 191 < T < 787 and 3.9x10^-4 < L_IR/L_* < 2.7x10^-1. We estimate a lower limit on the fraction of Tycho-2 Spectral Catalog main-sequence stars having mid-IR, but not near-IR, excesses to be 1.0+-0.3%.

We investigate the diffuse absolute calibration of the InfraRed Array Camera on the Spitzer Space Telescope at 8.0microns using a sample of 43 HII regions with a wide range of morphologies near GLON=312deg. For each region we carefully measure sky-subtracted,point-source- subtracted, areally-integrated IRAC 8.0-micron fluxes and compare these with Midcourse Space eXperiment (MSX) 8.3-micron images at two different spatial resolutions, and with radio continuum maps. We determine an accurate median ratio of IRAC 8.0-micron/MSX\8.3-micron fluxes, of 1.55+/-0.15. From robust spectral energy distributions of these regions we conclude that the present 8.0-micron diffuse calibration of the SST is 36% too high compared with the MSX validated calibration, perhaps due to scattered light inside the camera. This is an independent confirmation of the result derived for the diffuse calibration of IRAC by the Spitzer Science Center (SSC). From regression analyses we find that 843-MHz radio fluxes of HII regions and mid-infrared (MIR) fluxes are linearly related for MSX at 8.3-microns and Spitzer at 8.0 microns, confirming the earlier MSX result by Cohen & Green. The median ratio of MIR/843-MHz diffuse continuum fluxes is 600 times smaller in nonthermal than thermal regions, making it a sharp discriminant. The ratios are largely independent of morphology up to a size of ~24 arcsec. We provide homogeneous radio and MIR morphologies for all sources. MIR morphology is not uniquely related to radio structure. Compact regions may have MIR filaments and/or diffuse haloes, perhaps infrared counter- parts to weakly ionized radio haloes found around compact HII regions. We offer two IRAC colour-colour plots as quantitative diagnostics of diffuse HII regions.

The GLIMPSE (Galactic Legacy Mid-Plane Survey Extraordinaire) Point Source Catalog of ~ 30 million mid-infrared sources towards the inner Galaxy, 10 < |l| < 65 degrees and |b| < 1 degree, was used to determine the distribution of stars in Galactic longitude, latitude, and apparent magnitude. The counts versus longitude can be approximated by the modified Bessel function N=N_0*(l/l_0)*K_1(l/l_0), where l_0 is insensitive to limiting magnitude, band choice, and side of Galactic center: l_0= 17-30 degrees with a best fit value in the the 4.5 micron band of l_0=24 +/- 4 degrees. Modeling the source distribution as an exponential disk yields a radial scale length of H= 3.9 +/- 0.6 kpc. There is a pronounced north-south asymmetry in source counts for |l| < 30 degrees, with ~ 25% more stars in the north. For l=10-30 degrees, there is a strong enhancement of stars of m= 11.5-13.5 mag. A linear bar passing through the Galactic center with half-length R_bar=4.4 +/- 0.5 kpc, tilted by phi=44 +/- 10 degrees to the Sun-Galactic Center line, provides the simplest interpretation of this data. We examine the possibility that enhanced source counts at l=26-28 degrees, 31.5-34 degrees, and 306-309 degrees are related to Galactic spiral structure. Total source counts are depressed in regions where the counts of red objects (m_K-m_[8.0] >3) peak. In these areas, the counts are reduced by extinction due to molecular gas and/or high diffuse backgrounds associated with star formation.

We apply our three dimensional photoionization code to model Wisconsin H-alpha Mapper observations of the H II region surrounding the O9.5V star Zeta Oph. Our models investigate the porosity of the interstellar medium around zeta Oph and the effects of 3D densities on the H-alpha surface brightness and variation in the [N II]/H-alpha line ratio. The Zeta Oph H II region has a well characterized ionizing source, so it is an excellent starting point for 3D models of diffuse ionized gas. We investigate various hierarchically clumped density structures, varying the overall smoothness within the clumping algorithm. By simulating the observations, we can estimate the porosity of the medium in the vicinity of Zeta Oph and find that within the context of our hierarchically clumped models, around 50% to 80% of the volume is occupied by clumps surrounded by a low density smooth medium. We also conclude that in order for O stars to ionize the diffuse Warm Ionized Medium, the O star environment must be more porous than that surrounding Zeta Oph, with clumps occupying less than one half of the interstellar volume. Our clumpy models have irregular boundaries, similar to observed H II regions. However, in observed H II regions it is difficult to identify the precise location of the boundary because of the foreground and/or background emission from the widespread Warm Ionized Medium. This complicates the interpretation of the predicted rapid rise of some emission line ratios near the edge of uniform density H II regions and combined with the three dimensional clumpy nature of the interstellar medium may explain the apparent lack of distinctive emission line ratios near H I -- H II interfaces.

We present an extinction analysis of 9 paths through the LMC and SMC based on FUSE observations. To date, just two LMC sight lines have probed dust grain composition and size distributions in the Clouds using spectra including wavelengths as short as 950 A. We supplement these with results from 4 regions distinguished by their IR through UV extinction curves and grouped as LMCAvg, LMC2, SMC bar and SMC wing. Despite the distinct characters of Milky Way and Magellanic Cloud extinction, our results are generally analogous to those found for Galactic curves in that the FUSE portions of each extinction curve are described reasonably well by FM curves fitted only to longer wavelength data and lack any dramatic new extinction features, and any deviations from the CCM formalism continue into FUV wavelengths. An MEM analysis of these curves suggests that LMCAvg and SMC wing sight lines require more silicon and/or carbon in dust than current abundance measurements would allow, while the requirements for LMC2 and SMC bar sight lines do not fully tax the available reservoirs. An intermediate product of this analysis is the measurement of new H_2 abundances in the Magellanic Clouds. Collectively considering Cloud sight lines that possess significant H_2 column densities, E(B-V)/N(HI) ratios are reduced by significant factors relative to the Galactic mean, whereas the corresponding E(B-V)/N(H_2) values more closely resemble their Galactic counterpart. These trends reflect the fact that among these sight lines f(H_2)-values are lower than those common in the Milky Way for paths with similar degrees of reddening.

Spitzer IRAC 3.6-8 micron photometry obtained as part of the GLIMPSE survey has revealed mid-infrared excesses for 33 field stars with known spectral types in a 1.2 sq. degree field centered on the southern Galactic HII region RCW49. These stars comprise a subset of 184 stars with known spectral classification, most of which were pre-selected to have unusually red IR colors. We propose that the mid-IR excesses are caused by circumstellar dust disks that are either very late remnants of stellar formation or debris disks generated by planet formation. Of these 33 stars, 29 appear to be main-sequence stars based on optical spectral classifications. Five of the 29 main-sequence stars are O or B stars with excesses that can be plausibly explained by thermal bremsstrahlung emission, and four are post main-sequence stars. The lone O star is an O4V((f)) at a spectrophotometric distance of 3233+ 540- 535 pc and may be the earliest member of the Westerlund 2 cluster. Of the remaining 24 main-sequence stars, 18 have SEDs that are consistent with hot dusty debris disks, a possible signature of planet formation. Modeling the excesses as blackbodies demonstrates that the blackbody components have fractional bolometric disk-to-star luminosity ratios, LIR/L*, ranging from 10^-3 to 10^-2 with temperatures ranging from 220 to 820 K. The inferred temperatures are more consistent with asteroid belts rather than the cooler temperatures expected for Kuiper belts. Mid-IR excesses are found in all spectral types from late B to early K.

We have compared Monte Carlo photoionization models of H II regions with a uniform density distribution with models with the same central stars and chemical compositions but with 3-D hierarchical clumps. We compare the abundances of He, N, O, Ne, and S obtained from emission line strengths and [O III] and [N II] temperatures to those in our models. We consider stellar temperatures in the range 37.5 -- 45kK and ionizing luminosities from 10^{48} to 10^{51} photons/s. Clumped models have different ionic abundances than uniform. For hot stars, He^0/He^+ is 2 -- 3%, much larger than with uniform models. This amount of He I is independent of metallicity and so impacts the determination of the primordial abundance of He. The total abundances of O, Ne, and S obtained by the usual methods of analysis, using T([OIII) for high stages of ionization and T([NII]) for low, are about as accurate for clumped models as for uniform and within about 20% of the true values. If T([OIII]) is used for analyzing all ions, the derived (O/H) is 40 to 60% too large for cool stars but is good for hot stars. Uniform models have similar errors, so the clumping does not change the accuracy of abundance analysis. The physical causes of the ionic abundance errors are present in real nebulae. In clumped models, helium ionizing radiation from zones of high ionization (low He^0 and low UV opacity) can penetrate nearby regions near the edge of the ionized zone. This effect allows He^0 to absorb more stellar photons than in uniform or radially symmetrical geometries. In turn, these absorptions compete with O+, etc., for those energetic stellar photons.

We present a new planetary nebula, first identified in images from the Australia Telescope Compact Array, although not recognized at that time. Recent observations with the Spitzer Space Telescope during the GLIMPSE Legacy program have rediscovered the object. The high-resolution radio and infrared images enable the identification of the central star or its wind, the recognition of the radio emission as thermal, and the probable presence of polycylic aromatic hydrocarbons in and around the source. These lead to the conclusion that G313.3+00.3 is a planetary nebula. This object is of particular interest because it was discovered solely through radio and mid-infrared imaging, without any optical (or near-infrared) confirmation, and acts as a proof of concept for the discovery of many more highly extinguished planetary nebulae. G313.3+00.3 is well-resolved by both the instruments with which it was identified, and suffers extreme reddening due to its location in the Scutum-Crux spiral arm.

We present extinction curves that include data down to far ultraviolet wavelengths (FUV; 1050 - 1200 A) for nine Galactic sight lines. The FUV extinction was measured using data from the Far Ultraviolet Spectroscopic Explorer. The sight lines were chosen for their unusual extinction properties in the infrared through the ultraviolet; that they probe a wide range of dust environments is evidenced by the large spread in their measured ratios of total-to-selective extinction, R_V = 2.43 - 3.81. We find that extrapolation of the Fitzpatrick & Massa relationship from the ultraviolet appears to be a good predictor of the FUV extinction behavior. We find that predictions of the FUV extinction based upon the Cardelli, Clayton & Mathis (CCM) dependence on R_V give mixed results. For the seven extinction curves well represented by CCM in the infrared through ultraviolet, the FUV extinction is well predicted in three sight lines, over-predicted in two sight lines, and under-predicted in 2 sight lines. A Maximum Entropy Method analysis using a simple three component grain model shows that seven of the nine sight lines in the study require a larger fraction of grain materials to be in dust when FUV extinction is included in the models. Most of the added grain material is in the form of small (radii < 200 A) grains.

We determine and tabulate A(lambda)/A(K), the wavelength dependence of interstellar extinction, in the Galactic plane for 1.25um<lambda<8.0um along two lines of sight: l~42 and l~284. The first is a relatively quiescent and unremarkable region; the second contains the giant HII region RCW 49 as well as a ``field'' region unrelated to the cluster and nebulosity. Areas near these Galactic longitudes were imaged at J, H, and K bands by 2MASS and at 3--8um by Spitzer for the GLIMPSE Legacy program. We measure the mean values of the color excess ratios (A(lambda)-A(K))/(A(J)-A(K)) directly from the color distributions of observed stars. The extinction ratio between two of the filters, e.g. A(J)/A(K), is required to calculate A(lambda)/A(K) from those measured ratios. We use the apparent JHK magnitudes of giant stars along our two sightlines, and fit the reddening as a function of magnitude (distance) to determine A(J)/kpc, A(K)/kpc, and A(J)/A(K). Our values of A(lambda)/A(K) show a flattening across the 3--8um wavelength range, roughly consistent with the Lutz 1996 extinction measurements derived for the sightline toward the Galactic center.

Spitzer Space Telescope imaging as part of the Galactic Legacy Mid-Plane Survey Extraordinaire (GLIMPSE) reveals a previously unidentified low-latitude rich star cluster near l=31.3 degrees, b=-0.1 degrees. Near-infrared JHK' photometry from the Wyoming Infrared Observatory indicates an extinction of A_V ~ 15+/-3 mag for cluster members. Analysis of 13CO features along the same sightline suggests a probable kinematic distance of 3.1 -- 5.2 kpc. The new cluster has an angular diameter of ~1-2 pc, a total magnitude m_{K_0}=2.1 corrected for extinction, and a luminosity of M_K ~ -10.3 at 3.1 kpc. In contrast to young massive Galactic clusters with ages <100 Myr, the new cluster has no significant radio emission. Comparison to theoretical K-band luminosity functions indicates an age of at least several Gyr and a mass of at least 10^5 solar masses. Unlike known old open clusters, this new cluster lies in the inner Galaxy at R_{GC} ~ 6.1 kpc. We designate this object ``GLIMPSE-C01'' and classify it as a Milky Way globular cluster passing through the Galactic disk. We also identify a region of star formation and fan-shaped outflows from young stellar objects in the same field as the cluster. The cluster's passage through the Galactic molecular layer may have triggered this star formation activity.

GLIMPSE imaging using the Infrared Array Camera (IRAC) on the Spitzer Space Telescope indicates that star formation is ongoing in the RCW 49 giant H II region. A photometric comparison of the sources in RCW 49 to a similar area to its north finds that at least 300 stars brighter than 13th magnitude in band [3.6] have infrared excesses inconsistent with reddening due to foreground extinction. These are likely young stellar objects (YSOs) more massive than 2.5 M_sun suggesting that thousands more low-mass stars are forming in this cloud. Some of the YSOs are massive (B stars) and therefore very young, suggesting that a new generation of star formation is occurring, possibly triggered by stellar winds and shocks generated by the older (2-3 Myr) central massive cluster. The Spitzer IRAC camera has proven to be ideally suited for distinguishing young stars from field stars, and the GLIMPSE survey of the Galactic Plane will likely find thousands of new star formation regions.

We have calculated the radiation field, dust grain temperatures, and far infrared emissivity of numerical models of turbulent molecular clouds. When compared to a uniform cloud of the same mean optical depth, most of the volume inside the turbulent cloud is brighter, but most of the mass is darker. There is little mean attenuation from center to edge, and clumping causes the radiation field to be somewhat bluer. There is also a large dispersion, typically by a few orders of magnitude, of all quantities relative to their means. However, despite the scatter, the 850 micron emission maps are well correlated with surface density. The fraction of mass as a function of intensity can be reproduced by a simple hierarchical model of density structure.

We demonstrate that the observed increase of some nebular line ratios with height above the midplane in the diffuse ionized gas (DIG) in the Milky Way and other galaxies is a natural consequence of the progressive hardening of the radiation field far from the midplane ionizing sources. To obtain increasing temperatures and line ratios away from the midplane, our photoionization simulations of a multi-component interstellar medium do not require as much additional heating (over and above that from photoionization) as previous studies that employed one dimensional, spherically averaged models. Radiation leaking into the DIG from density bounded H II regions is generally harder in the H-ionizing continuum and has its He-ionizing photons suppressed compared to the ionizing source of the H II region. In line with other recent investigations, we find that such leaky H II region models can provide elevated temperatures and line ratios, and a lower He+ fraction in the DIG. For a composite model representing the relative spectral types of O stars in the solar neighbourhood, we find that additional heating less than 1.0E-26 n_e ergs/s/cm^3 can reproduce the observed elevated line ratios in the DIG. This additional heating is considerably less than previous estimates due to the natural hardening of the radiation field reaching large heights in our simulations.

We have developed a three dimensional Monte Carlo photoionization code tailored for the study of Galactic H II regions and the percolation of ionizing photons in diffuse ionized gas. We describe the code, our calculation of photoionization, heating & cooling, and the approximations we have employed for the low density H II regions we wish to study. Our code gives results in agreement with the Lexington H II region benchmarks. We show an example of a 2D shadowed region and point out the very significant effect that diffuse radiation produced by recombinations of helium has on the temperature within the shadow.

GLIMPSE (Galactic Legacy Infrared Mid-Plane Survey Extraordinaire), a SIRTF Legacy Science Program, will be a fully sampled, confusion-limited infrared survey of the inner two-thirds of the Galactic disk with a pixel resolution of \~1.2" using the Infrared Array Camera (IRAC) at 3.6, 4.5, 5.8, and 8.0 microns. The survey will cover Galactic latitudes |b| <1 degree and longitudes |l|=10 to 65 degrees (both sides of the Galactic center). The survey area contains the outer ends of the Galactic bar, the Galactic molecular ring, and the inner spiral arms. The GLIMPSE team will process these data to produce a point source catalog, a point source data archive, and a set of mosaicked images. We summarize our observing strategy, give details of our data products, and summarize some of the principal science questions that will be addressed using GLIMPSE data. Up-to-date documentation, survey progress, and information on complementary datasets are available on the GLIMPSE web site: www.astro.wisc.edu/glimpse.

By radiation transfer models with a realistic power spectra of the projected density distributions, we show that the optical properties of grains are poorly constrained by observations of reflection nebulae. The ISM is known to be hierarchically clumped from a variety of observations (molecules, H I, far-infrared). Our models assume the albedo and phase parameter of the dust, the radial optical depth of the sphere averaged over all directions, and random distributions of the dust within the sphere. The outputs are the stellar extinction, optical depth, and flux of scattered light as seen from various viewing angles. Observations provide the extinction and scattered flux from a particular direction. Hierarchical geometry has a large effect on the flux of scattered light emerging from a nebula for a particular extinction of the exciting star. There is a very large spread in both scattered fluxes and extinctions for any distribution of dust. Consequently, an observed stellar extinction and scattered flux can be fitted by a wide range of albedos. With hierarchical geometry it is not completely safe to determine even relative optical constants from multiwavelength observations of the same reflection nebula. The geometry effectively changes with wavelength as the opacity of the clumps varies. Limits on the implications of observing the same object in various wavelengths are discussed briefly. Henry (2002) uses a recipe to determine the scattered flux from a star with a given extinction. It is claimed to be independent of the geometry. It provides considerably more scattering than our models, probably leading to an underestimate of the grain albedos from the UV Diffuse Galactic Light.

Based on the ISO spectral catalogue of compact HII regions by Peeters et al. (2001), we present a first analysis of the hydrogen recombination and atomic fine-structure lines originated in the ionized gas. The sample consists of 34 HII regions located at galactocentric distances between Rgal = 0 and 15 kpc. The SWS HI recombination lines between 2 and 8 mum are used to estimate the extinction law at these wavelengths for 14 HII regions. An extinction in the K band between 0 and $\sim$ 3 mag. has been derived. The fine-structure lines of N, O, Ne, S and Ar are detected in most of the sources. Most of these elements are observed in two different ionization stages probing a range in ionization potential up to 41 eV. The ISO data, by itself or combined with radio data taken from the literature, is used to derive the elemental abundances relative to hydrogen. The present data thus allow us to describe for each source its elemental abundance, its state of ionization and to constrain the properties of the ionizing star(s).

Observations of the "Warm Ionized Medium" (or, equivalently, the "Diffuse Ionized Gas") of the local ISM, the Perseus arm in the Milky Way, and also in several other galaxies show strong [NII]6563 (~H-alpha in some cases) and [SII]6717/[NII]6583 = 0.6 - 0.7 in all locations and objects. Other line ratios (e.g., [O III]5007/H-beta) vary considerably. Simple photoionization models reproduce the observed spectra, providing extra heating beyond that supplied by photoionization is assumed (Reynolds, Haffner, & Tufte 1999). With observed gas-phase abundances (not solar), the line ratios in the local arm at b = 0 deg are fitted with no extra heating and (S/H) = 13 ppm (solar is 20 ppm). Local gas observed at b = -35 deg requires extra heating of about gamma = 0.75, where gamma is the extra heating in units of 10^{-25} erg H^{-1} s^{-1}. In the Perseus arm, there are similar results, with a domposition consistent with the Galactic abundance gradient. The requirements for NGC 891 are similar to the Perseus arm: little or no extra heating at |z| = 1 kpc and gamma 3 at 2 kpc. In NGC 891 there is also an increase of 5007/H-alpha with |z| that can only come about if most of the ionizing radiation is supplied by stars with T~50000 K. Either their radiation must propagate from the plane to high |z| through very little intervening matter, or else the stars are located at high |z|. The total power requirement of the extra heating is <15% of the photoionization power. [O~II]3727/H-beta can serve as a useful diagnostic of extra heating, but [S~III] 9065,9531/H-alpha is not useful in this regard.

Sightlines through high- and intermediate-velocity clouds allow measurements of ionic gas phase abundances, A, at very low values of HI column density, N(HI). Present observations cover over 4 orders of magnitude in N(HI). Remarkably, for several ions we find that the A vs N(HI) relation is the same at high and low column density and that the abundances have a relatively low dispersion (factors of 2-3) at any particular N(HI). Halo gas tends to have slightly higher values of A than disk gas at the same N(HI), suggesting that part of the dispersion may be attributed to the environment. We note that the dispersion is largest for NaI; using NaI as a predictor of N(HI) can lead to large errors. Important implications of the low dispersions regarding the physical nature of the ISM are: (a) because of clumping, over sufficiently long pathlengths N(HI) is a reasonable measure of the_local_ density of_most_ of the H atoms along the sight line; (b) the destruction of grains does not mainly take place in catastrophic events such as strong shocks, but is a continuous function of the mean density; (c) the cycling of the ions becoming attached to grains and being detached must be rapid, and the two rates must be roughly equal under a wide variety of conditions; (d) in gas that has a low average density the attachment should occur within denser concentrations.

Through a program of narrow band imaging, we have observed the changing structure of the H alpha emission line around Nova Cyg 1992 (V1974 Cyg) at regular intervals from 1996 to 1999. Between 1994 and 1996, the nebular boundary advanced to the southwest at nearly the speed of light, implying that the nebula was created by an expanding wave of radiation originating in the explosion interacting with surrounding material. The expansion speed dropped to 0.35c during 1996-1999. We have taken spectra of the nebula in 1998 and 1999. Only Balmer lines are detected, no He I, [OIII], [OII], [NII], or [SII]. There is also no trace of the high excitation nova lines (He II, NeV, etc). These spectra show that the nebula is NOT a reflection nebula, a conventional HII region, or a shock involving motions of the gas. The integrated H alpha luminosity of the nebula between 1996 and 1999 is in the range =~1.3-2.2 x 10(35) erg s(-1). The density is poorly determined, but is probably very large (10(4) cm(-3)) in order to explain the brightest region of the nebula. The dynamical timescale is about a year and the recombination timescale of the same order. Bright patches are observed to fade in these times. The energy required to ionize the nebula is the bolometric luminosity of the nova for 30 days, smaller than the time during which the temperature of the nova photosphere was in the right range to produce the ionizing photons. We have also undertaken sensitive surveys of H alpha nebulae around recent novae but find no evidence of other such nebulae, so this type of object must be rare.

Interstellar grains may be composite collections of particles of distinct materials, including voids, agglomerated together. We determine the various optical cross sections of such composite grains, given the optical properties of each constituent, using an approximate model of the composite grain. We assume it consists of many concentric spherical layers of the various materials, each with a specified volume fraction. In such a case the usual Mie theory can be generalized and the extinction, scattering, and other cross sections determined exactly. We find that the ordering of the materials in the layering makes some difference to the derived cross sections, but averaging over the various permutations of the order of the materials provides rapid convergence as the number of shells (each of which is filled by all of the materials proportionately to their volume fractions) is increased. Three shells, each with one layer of a particular constituent material, give a very satisfactory estimate of the average cross section produced by larger numbers of shells. We give the formulae for the Rayleigh limit (small size parameter) for multi-layered spheres and use it to propose an ``Effective Medium Theory'' (EMT), in which an average optical constant is taken to represent the ensemble of materials. Multi-layered models are used to compare the accuracies of several EMTs already in the literature.