the yuan - tseh lee array for microwave background anisotropy ( amiba ) is a platform - mounted 7-element interferometer operating at 3-mm wavelength to study the structure of the cosmic microwave background ( cmb ) radiation .
it is being constructed as part of the cosmology and particle astrophysics ( cospa ) project , funded by the taiwan ministry of education initiative on academic excellence .
this excellence initiative was aimed at stimulating interdisciplinary research and large scale integration of independent research programs .
cospa is designed to jump start a program of research in cosmology , with both theory and experimental projects , while incorporating research in high energy physics , development of infrastructure for optical astronomy in taiwan , as well as accessing observing time on a 4-m class optical telescope .
amiba is a collaboration between principally the academia sinica institute of astronomy and astrophysics ( asiaa ) , the national taiwan university ( ntu ) physics and electrical engineering departments , and the australia telescope national facility ( atnf ) .
the project was started in 2000 .
a two - element prototype was deployed in 2002 to mauna loa ( elevation 3396@xmath0 m ) in hawaii for testing of design concepts .
site development was completed in 2004 .
the amiba mount was delivered and installed in 2004 , while the platform was delivered and integrated in 2005 . with the integration of the first seven elements of the array and successful first light ,
the amiba was dedicated in october 2006 , and named after then academia sinica president yuan tseh lee for his important contributions in promoting the growth of astronomy in taiwan .
figure 1 shows the amiba at the dedication ceremony .
the aim of amiba is to study spatial structure in the cmb radiation ( ho et al .
2008 ; wu et al .
2008a ) , which carries imprints of various physical processes in early epochs of the universe . since its initial detection by penzias & wilson ( 1965 ) ,
the cmb has been recognized as the definitive signature of the big bang which began the expansion of the universe .
subsequent studies have established the properties of this relic radiation after its decoupling from the matter in the early universe around @xmath1 : a mean temperature of @xmath2k ( present ) with minute fractional anisotropies at the level of @xmath3 ( cobe , mather et al .
1990 ; smoot et al .
1992 ) , and polarization at the level of a few to 10% of temperature fluctuations ( dasi , kovac et al .  2002 : wmap , kogut et al .
2003 ; page et al .  2007 : nolta et al .
2008 : cbi , readhead et al .
2006 : quad , pryke et al .
in particular , the cmb structures seen on various angular scales by cobe and then wmap ( bennett et al .  2003 ; spergel et al .  2007 ;
komatsu et al .
2008 ) demonstrated that the angular power spectrum of cmb anisotropies is a powerful probe of our cosmological model of the universe .
amiba is built to sample the angular range from @xmath4 to @xmath5 , corresponding to spherical harmonic multipoles @xmath6@xmath7 , at a wavelength of 3@xmath0 mm , with full polarization .
these capabilities complement existing , on - going , and planned experiments .
the angular scales sampled by amiba address the higher - order acoustic peaks of the cmb structures to further constrain cosmological models .
amiba also aims to search for , and study , distant high - redshift clusters of galaxies whose hot intracluster gas will distort the cmb spectrum via the thermal sunyaev zeldovich effect ( hereafter sze , sunyaev & zeldovich 1970 , 1972 ; rephaeli 1995 ; birkinshaw 1999 ; carlstrom , holder , & reese 2002 ) . the optical and x - ray surface brightnesses of clusters of galaxies decrease rapidly with increasing redshift due to cosmological redshift dimming , while the detectable sze is close to being independent of redshift because it is a spectral distortion of the cmb radiation which itself increases in intensity with increasing redshift as @xmath8 in a standard cosmological model .
thus sze measurements are potentially more sensitive than x - ray observations for finding clusters of galaxies beyond a redshift @xmath9 , and will be an important probe for the matter distribution in the high - redshift universe .    in this paper , we describe the design and construction of amiba , the first observational results , and the scientific potential of this instrument . throughout this paper ,
we adopt a concordance @xmath10cdm cosmology with @xmath11 , @xmath12 , and @xmath13 .
the basic characteristics of amiba are summarized in table 1 .        with the funding of the cospa / amiba projects in 2000 ,
a workshop was held to define the scientific objectives and the design of amiba ( lo et al .
2001 ; liang 2002 ) . _
the first design criterion was to operate at 3-mm wavelength .
_ this was to take advantage of the sweet spot at 3@xmath0 mm where the fractional sze decrement with respect to the primary cmb is close to its maximum ( see figure 1 of zhang et al .
2002 ) and the sze signal is minimally - contaminated by the galactic synchrotron emission , dust foregrounds , and the population of cluster and background radio sources .
operations at 3@xmath0 mm also complement the wavelength coverage of other existing and planned cmb instruments : interferometers such as cbi at 30@xmath0ghz ( padin et al .
2002 ) , ami at 15@xmath0ghz ( scaife et al .
2008 ) , sza ( mroczkowski et al .  2008 ) at 30 and 90@xmath0ghz , and vsa at 30@xmath0ghz ( watson et al .
2002 ) ; bolometer arrays such as act , apex - sz ( halverson et al .  2008 ) , and spt . _
the second design criterion was to choose interferometry .
_ this was a somewhat difficult choice as many new cmb projects were then planning to use bolometer arrays , which held the promise of a greater inherent sensitivity because of the broad wavelength coverage and a greater speed because of the multiple elements of the detector arrays .
the choice of interferometry was based on the desire to utilize cross correlations to suppress systematic effects , since the ability for bolometers to integrate down to theoretical noise was unknown at that point . furthermore , asiaa had the experience of working in millimeter wavelength interferometry from being a partner on the submillimeter array ( sma ) project ( ho et al .
amiba was seen as an extension and application of the technical capabilities within asiaa .
interferometry is also a natural way to sample simultaneously the spatial structures on various scales and to construct a map by fourier inversion .
imaging the entire primary beam of the individual elements of an interferometer is equivalent to a multi - element detector array .
_ the third design criterion concerns the angular sizes to be pursued .
_ amiba was specifically designed to sample structures in the cmb on small angular scales at multipoles @xmath6@xmath7 , complementing in angular scales earlier large - sky cmb experiments ( e.g. , cobe , maxima , boomerang ) and the dasi experiment on degree / sub - degree angular scales ( @xmath14@xmath15 ; leitch et al .
2002a , 2002b ) . in the left panel of figure 2
, we show the initial compact configuration of seven 0.6@xmath0 m antennas ( _ small solid circles _ ) on the 6@xmath0 m platform ( _ solid outer circle _ ) , along with the distribution of holes ( _ crosses _ ) for mounting receiver packages . at each of the frequency channels centered at about 90 and 98@xmath0ghz ,
this compact configuration provides 21 simultaneous baselines with three baseline lengths of @xmath16 , 1.05 , and 1.21@xmath0 m , corresponding to angular multipoles @xmath17 of @xmath18 at @xmath19 . also shown in the right panel of figure 2
is the sensitivity of the 7-element compact array as a function of @xmath20 at two frequency channels , represented by window functions ( see eq .
[ 14 ] of white et al .  1999 ) for the three different baseline - lengths . for the initial configuration of amiba at 3@xmath0 mm ( see figures 2 and 3 ) , we are sensitive to the multipole range of @xmath21 , which roughly overlaps with the cbi experiment at 1 cm ( padin et al .  2002 ; mason et al .  2003 ; pearson et al .
2003 ) but is complementary in observing frequency .
hence amiba aimed at angular scales from @xmath22 to @xmath23arcmin , in order to extend the coverage of angular sizescales by one more order of magnitude . the choice of 3-mm wavelength meant the required baselines were 0.6@xmath0 m to 6@xmath0 m .
a maximum baseline of 6@xmath0 m immediately suggested that this interferometer is small enough to be mounted on a platform , as pioneered by cbi and dasi .
_ the fourth design criterion was therefore to choose a carbon fiber platform for weight and stiffness considerations . _
the platform was made with multiple holes for mounting the receiver packages , to accommodate multiple baseline configurations .
the layout of the receiver ports has a hexagonal pattern in order to allow close packing , and also to utilize the triangular patterns used by small interferometers such as the sma , cbi , and dasi , in order to achieve the most uniform @xmath24-coverage .
an asiaa team led by robert n. martin and philippe raffin designed the platform , which was manufactured by composite mirrors applications .
a single rigid platform has the advantages of stable differential pointing , a stable baseline solution without differential delay tracking , no mutual shadowing by individual reflectors , and a single drive system .
however , a platform also means the interferometer does not have different projected baselines because of earth rotation .    _
the fifth design criterion was the choice of a hexapod mount , which provides six degrees of freedom ( which are tightly constrained ) in driving the telescope and rotating the platform .
_ by rotating the reflectors with respect to the celestial source , additional @xmath24-spacings are sampled .
the rotation of the platform , described by the so - called polarization angle of the platform , also allows us to set different orientations of the receivers relative to the ground , which is useful for discriminating between various environmental effects and checking for ground pick - up .
these capabilities were introduced both by cbi and dasi , and we have incorporated them . the hexapod mount was designed and built by vertex antennentechnik , duisburg , germany . to protect the telescope from the elements , we use a retractable shelter with seven steel trusses covered with a pvc fabric , which is manufactured by american spaceframe fabricator inc .
the individual parabolic reflectors are built in carbon fiber in order to minimize their weight .
the choice of 0.6 m ( @xmath25 m , to be precise ) for the reflector diameter was driven by the desire to cover angular sizescales from @xmath22@xmath23arcmin in order to extend and overlap with existing cmb data .
the choice of 1.2@xmath0 m for the diameter of the second - generation reflectors was to increase the collecting area and to improve the spatial dynamic range and angular resolution for sze observations while utilizing the entire platform .
the reflectors were designed by asiaa , and manufactured by cotech , taichung , taiwan .
the design includes baffles to shield against crosstalk between individual elements , and gortex covers to shield against direct solar irradiation .
the height of the baffles is approximately 30% of the diameter of the reflectors so that the secondary mirrors are well shielded .
the supports for the secondary are attached to the baffles which are attached to the edges of the primary . replicating an accurate carbon fiber surface against a steel mold
was not simple , and some hand polishing was required .
depositing the actual reflecting surface and a protective coating was also not simple especially for the 1.2@xmath0 m size .
nevertheless , a collaborative effort between asiaa and cotech was successful in delivering the reflectors with surface accuracy better than 50 microns .
this is 2% of the operating wavelength of 3@xmath0 mm , and would allow efficient operation up to 1 mm .
laboratory measurements and outdoor beam pattern measurements showed that the performances of the reflectors met the specifications ( for details see koch et al .
2006 ) .      _
the sixth design criterion was to choose heterodyne receiver systems , which operate between 86 and 102@xmath0ghz .
_ this choice of frequency covers an excellent region of atmospheric transparency .
the design , construction , and integration of the receiver systems were performed by asiaa staff led by ming - tang chen .
dual polarization capabilities are provided by waveguide orthomode transducers , which follow the circular corrugated feedhorns and the circular - to - linear polarizers . each polarization
is then fed into a jpl monolithic - microwave - integrated - circuit ( mmic ) inp hemt low noise amplifier ( lna ) cooled to 15@xmath0k with standard cti22 refrigerators .
these lnas performed well with measured noise temperatures of 3550@xmath0k across a 20@xmath0ghz bandwidth ( weinreb et al .
we used subharmonically pumped mixers ( shm ) using the same mmic technology .
the lo / if systems are designed and built by tah - hsiung chu and his group at the ntu ee - department .
the phase - locked lo signal at 21@xmath0ghz is doubled to 42@xmath0ghz and then phase switched with walsh functions before being combined with the sky signal at the shm .
variable attenuators and amplifiers , before and after mixing , control the if levels before correlation . slope equalizers , phase stabilized cables , and adjustable delays are used to further adjust the if signals .
we measure the effective receiver temperatures across the @xmath22@xmath26ghz if window to be @xmath27@xmath28k .
figure 3 shows the amiba fully loaded with all the receivers and correlator modules .
more details on the receiver system are given by chen et al .
( 2008 ) .      _
the seventh design criterion was to choose a wide band correlator with analog technology .
_ this is a joint development effort between asiaa and atnf , with chao - te li , derek kubo , and warwick wilson leading the effort ( li et al .
a digital correlator would have required sampling at too high a rate to be practical .
the amiba correlator uses balanced diode mixers to multiply the signals from each pair of antennas .
after application of the 4 different lags , the cross - correlated signal is then digitized and read out with integrated circuits designed by tzi - dar chiueh of the ntu ee department .
the readout ics demodulate the phase - switched signals and accumulate the counts for specified times .
the choice of a small number of lags is to reduce bandwidth smearing effects while maintaining the maximum bandwidth response through the cascade of electronics .
the bandwidth smearing effect is described by the product of the fractional bandwidth and the displacement of the source from the field center .
a @xmath29 bandwidth will keep the distortion below @xmath30 ( @xmath31 ) at the fwhm of the primary beam of the @xmath32 m
( @xmath33 m ) reflectors .
because of the small number of lags to cover a large bandwidth , the conversion to complex visibility is strongly affected by gain variations over the passband , differential delays between lags , and non - linear phase response for each lag .
the lag - to - visibility transformation is calibrated with a noise source .
this is discussed further by lin et al .
( 2008 ) .
the hexapod drive system is more complex than the conventional azimuth - elevation drive systems used in radio astronomy .
the difficulty lies in the necessity of driving all six hexapod jacks without over - extension or collisions .
vertex antennentechnik provided the control software .
the performance of the drive system was checked via pointing and tracking tests ( koch et al .
pointing of the amiba utilized an optical telescope mounted on the platform .
we identified a misalignment of the anchor cone and the mount , tilt of the optical telescope with respect to the mount , flexure in the platform as a function of the platform polarization angle , and local platform deformation .
we also measured the repeatability of pointing on both short and long timescales , and we performed photogrammetry to measure the stability of the platform . pointing was found to be repeatable at the @xmath34 level on a timescale of several hours .
the absolute rms pointing error appears to be about @xmath35 for the platform set at zero polarization angle , and up to 3@xmath36 if the platform is allowed to rotate to different polarization angles .
however , with the implementation of an interpolation table , absolute rms pointing can be reduced to @xmath37 over all sky .
this is less than @xmath29 of the primary beams of the 0.6@xmath0 m reflectors , and is within specifications for the operation of the 7-element compact array ( koch et al .
however , when the 1.2@xmath0 m reflectors are deployed , pointing needs to be improved by a factor of 2 .
the segmented approach instead of a monolithic design for our carbon fiber platform resulted in unforeseen difficulties .
the bolted joints between the six outer segments were not stiff enough . in spite of efforts to strengthen the joints with additional plates and brackets
, the platform still deforms under operational loads to a level beyond our specifications .
a saddle - shaped deformation pattern is present with an amplitude dependent on the hexapod azimuth , elevation and polarization position .
fortunately , the deformation appears repeatable when measured with photogrammetry which means that it can be modeled .
maximum deformation at the edge of the platform under simulated full loading is 0.38@xmath0 mm .
the @xmath0maximum deformation in the inner 3@xmath0 m of the platform is more modest at 0.120@xmath0 mm , and within specifications for the operation of the 7-element compact array .
the deformation is due principally to the platform not being stiff enough , and also because the stresses from the drive system are transmitted to the platform in spite of a steel interface ring ( koch et al .
the immediate ramification is that the individual primary beams of the interferometer elements will be mis - pointed relative to each other .
this is equivalent to the primary mirror of a single dish telescope being deformed under gravity or atmospheric distortion .
the equivalent adaptive optics approach for an interferometer involves adjusting the gains and phases of the individual elements either through a look - up table or via self calibration if the signals are strong enough .
the advantage of an interferometer is that we can adjust the signals during the postprocessing phase before adding or multiplying the signals from different baselines . for our current configuration , such correction schemes for
the platform deformation do not need to be implemented .
calibration of the amiba system was led by kai - yang lin and chao - te li .
the receiver temperatures were measured by the standard hot / cold load method to be 55 to @xmath28k .
the contribution from the sky and ground pickup to the total system temperatures is about @xmath38k .
since the correlator only has 4 lags , only 2 frequency channels can be extracted .
the response of the correlator is therefore not easy to measure , and we use planets to calibrate the correlator response periodically .
we find the complex gain to be stable to 5% in gain and 0.1 rad in phase on time scales of a few hours .    using the detected fluxes for saturn and jupiter , the overall system efficiency was estimated to be @xmath39@xmath40 for each of the baselines of the interferometer .
this efficiency accounts for all losses due to illumination , blockage , spillover , alignment of the reflectors with respect to the radio axis of the platform , deformation of the platform , and pointing , as well as the narrower correlator response ( lin et al .
2008 ) . the losses in front of the receiver account for the significant contributions from the sky and ground pickup to the system temperatures .
in two - patch differencing observations a sensitivity of about @xmath41mjy is achieved in 1@xmath0hour of on - source integration under good sky conditions ( chen et al .
2008 ; lin et al .
2008 ) .
since the sze signals and the cmb anisotropy are quite faint , the behavior of the noise is very important .
this analysis has been led by hiroaki nishioka .
we conducted tests of blank sky data as well as data where the inputs are terminated by absorbers .
statistical analysis of different samples of such data showed no significant differences .
time variable signals from the correlator due to electronics and ground pickups can be seen in the various datasets . however , a power spectrum analysis demonstrated white noise behavior for frequencies between 10@xmath42 and @xmath43hz .
hence two - patch differencing of the sky signal at 600  sec intervals is adequate to remove the slowly varying contaminations .    by cross
correlating the outputs from the different lags of the correlator , we find only a weak correlation at less than 10% . by dividing the datasets into smaller samples , we verify that the sample variance is consistent with gaussian random noise which would integrate down over time .
we also applied the kolmogorov - smirnov ( k - s ) test to our cluster data , and we find that 90% of the data are consistent with gaussian behavior at the 5% significance level .
the data which fail the test at 5% level are often found to be associated with hardware problems , and we use this criterion to discard questionable samples .
more details on the gaussianity tests are given in nishioka et al .
( 2008 ) .
since the sze signals are faint , they can be confused with the cmb anisotropies and foreground emission structures .
the analysis of possible contaminations has been led by guo - chin liu . examining the wmap data
, we can see that cmb structures on the scale of a degree are suppressed by our two - patch observing technique , while contamination from galactic emission is typically fainter than the cmb emission by an order of magnitude . from wmap data
, we also estimate possible contamination from primary cmb anisotropies , by analyzing 500 simulated cmb fields , for the sizescales which we sample .
we find contamination at the level of @xmath44 to @xmath45mjy for the sampled multipoles of @xmath46 to @xmath47 ( at @xmath48ghz ) .
we also estimated potential contamination from foreground point sources .
the best way to correct for this issue is to have high angular resolution interferometer data at the observed frequencies .
lacking such data , we extrapolate the contributions from known point sources in low frequency radio surveys .
since spectral indices are not known accurately to @xmath49 mm , a statistical approach is used to estimate the contamination in each cluster from known point sources . in all clusters ,
a net positive contribution of the point sources was found in our main - trail / lead differencing amiba observations , indicating that there are more radio sources towards clusters than in the background ( as expected ; see coble et al .
the corrections for contamination are described more fully in liu et al .
( 2008 ) .
first light with the 7-element amiba was achieved in september 2006 .
the array of 0.6@xmath0 m reflectors was in the close - packed configuration on the platform ( figures 2 and 3 ) .
the image of jupiter shown in figure 4 is the first end - to - end test of the system based on drift - scan data , including the pipeline for data analysis .
the image was constructed from drift scans at four equally - divided platform polarization angles to provide better @xmath24-coverage .
only the transit data of about 0.23@xmath0s exposure time from each drift scan were used , yielding an overall noise level in the dirty map of about 3@xmath0jy , much lower than the jupiter flux density of 844@xmath0jy .
successful images were obtained for saturn , venus , and the crab nebula , all with fluxes on the order of 200@xmath0jy . in particular , the measured angular size of the crab nebula
is consistent with existing optical and radio images . in december 2006 , using saturn as the calibrator , uranus was imaged with a flux density of @xmath50jy ( at @xmath51 confidence level ; in 16 @xmath43s exposures ) , consistent with the expected level of @xmath52jy .      at the achieved sensitivity ,
the close - packed 7-element configuration of amiba is suitable for imaging galaxy clusters via the sze .
during 2007 , we imaged 6 clusters , chosen to be relatively nearby ( a median redshift of @xmath53 ) , already detected in sze , and resolvable with amiba .
target selection and data analysis are described in wu et al .
( 2008b ) . to detect the much fainter sze signals from clusters of galaxies , we need to integrate longer .
at this level of sensitivity , emission from local terrestrial sources such as buildings and the ground will enter through the system either through the sidelobes of the antenna response or via scattering into the primary beam .
the interferometer has the advantage over a bolometer array in that the cross correlation will suppress extended low level emission which enters the system through the sidelobe response of the reflectors .
however , structures on angular scales corresponding to the interferometer baselines will still persist .
the design of the baffles around the edges of the reflectors helps further suppression , but we did not exercise the option of installing ground shields . instead we employ a two- or three - patch observing technique where the target source is preceded and/or followed by tracking over the same range of azimuth and elevation just traversed by the source .
this allows us to subtract and cancel the terrestrial ground and sky emission .
this is the amiba equivalent of the position - switching technique used in single - dish radio astronomy , where we track specifically the same part of the atmosphere . as shown in figure 5 , and in wu et al .
( 2008b ) , this procedure is quite successful and we have mapped the sze decrement towards a number of clusters . these detections serve to demonstrate the scientific potentials of amiba .
amiba will be sensitive to structures as large as @xmath54 with resolution as high as @xmath55 so that sze structures can be resolved . even at the current angular resolution of @xmath56 fwhm , an extended
, elliptical shape can be seen for some of the clusters .
these shapes can be compared to the results from weak gravitational lensing or x - ray imaging .
the sze and x - ray images are sensitive to the distribution of the hot intracluster gas , while the weak lensing images probe the distribution of all the gravitational mass including the dominant `` invisible '' dark matter .
furthermore , the x - ray emission is sensitive to the square of the electron density while the sze measures the thermal electron pressure projected along the line - of - sight , linearly tracing the electron column density .
hence , the complementary x - ray , the sze , and the weak lensing data , probe progressively further out in the mass distribution .    for four of the massive clusters in our sample ( a1689 , a2142 , a2261 , a2390 ) , for which high - quality subaru weak lensing data are available , we have compared the sze results with the subaru weak lensing measurements ( umetsu et al .  2008 ) . in the case of abell cluster a2390 ,
the elliptical shape seen in sze is consistent with the shape of the dark matter distribution as deduced from weak lensing data .
for this comparison , a2142 is of particular interest since this is our brightest , most - nearby ( resolvable ) sze cluster at @xmath57 , known as a merging cluster with two x - ray _ cold fronts _ , which are sharp edges in x - ray dense cores ( markevitch et al .
2000 ) . in figure 6
we compare for a2142 our amiba sze map with the projected mass distribution ( _ white contours _ ) as deduced from our weak lensing analysis . here
the mass map is smoothed to a resolution of @xmath55 fwhm .
our amiba sze map shows an elliptical structure extending in the northwest ( nw ) - southeast ( sw ) direction , similar to x - ray and weak lensing distribution shapes ( figure 6 ; see also figure 10 of okabe & umetsu 2008 ) .
the almost @xmath54 angular extent of the cluster sze justifies the use of observing with antennas as small as @xmath32 m , so that the observations are sensitive to large scale structures .
further , relative positions ( between pixels ) are important in such maps of extended structures , where the interferometer has the inherent advantage of relying on the phase to determine relative positions , while single dish studies might be affected more strongly by pointing of the telescopes .    in the cluster a2142
, a mass subclump located @xmath58 ( @xmath59kpc@xmath60 at the cluster redshift ) to the nw of the cluster center has been detected from the subaru weak lensing analysis of okabe & umetsu ( 2008 ) and umetsu et al .
this nw mass subclump , lying @xmath61 ahead of the nw- edge of the x - ray dense core , has no x - ray counterpart in chandra and xmm - newton observations , but is associated with a slight excess of cluster sequence galaxies ( see figure 10 of okabe & umetsu 2008 ) .
our amiba map exhibits slight excess sze signals in the nw region ( figure 6 ) at @xmath62 significance levels .
we note this slight excess sze appears extended over a couple of synthesized beams , although the per - beam significance level is marginal .
good consistency found between the sze and weak lensing maps is encouraging , and may suggest that the nw excess sze is a pressure increase in the icm associated with the moving nw substructure .
clearly further improvements in both sensitivity and resolution of the sze mapping ( see  [ subsec : expansion ] ) are needed to better constrain the merger geometry and physical properties of the merging substructure .
this however demonstrates the potential of sze observations as a powerful tool for measuring the distribution of icm in cluster outskirts where the x - ray emission measure ( @xmath63 ) is less sensitive .
this also demonstrates the potential and reliability of amiba , and the power of multiwavelength cluster analysis for probing the distribution of mass and baryons in clusters .    finally , for our small sample with a mean virial mass of @xmath64 , the hot gas mass fraction is about @xmath65 ( umetsu et al .
2008 ) . as compared to the cosmic baryon fraction @xmath66 of @xmath67 , as deduced from the wmap 5-year data ( dunkley et al .
2008 ) , possibly @xmath68 of the baryons are missing from the hot cluster environment ( umetsu et al .  2008 ) .
this missing cluster baryon fraction is partially made up by observed stellar and cold gas fractions of @xmath69 several @xmath70 in our x - ray temperature range ( @xmath71kev ) .      by comparing the sze and x - ray imaging results
, we can deduce a value for the hubble constant @xmath72 .
since the sze is induced by inverse - compton scattering of the cmb photons by hot electrons in the cluster environment , it is sensitive to the total column of electrons .
the x - ray emission is due to the bremstrahlung process between protons and electrons , and is therefore sensitive to the square of the electron density . by combining these two relations ,
the angular diameter distance can be deduced .
however this derivation , and the subsequent measurement of @xmath72 depends on the assumed geometry of the emission region .
we have derived a value of @xmath72 from our sample of sze results combined with x - ray data .
the uncertainty is high because of the small size of our sample , and the bias from the assumed geometry .
however , the deduced value is consistent with other sze constraints , and is an indicator of integrity of the amiba data so far .
better limits will come with a much larger sample to be obtained in the future .
while the initial 7-element @xmath32 m reflectors have been commissioned , we are proceeding with the expansion of the amiba to its 13-element configuration . in this configuration
, we will upgrade from @xmath32 m to @xmath33 m antennas .
this will increase the collecting area by a factor of @xmath73 , and the speed of the interferometer by a factor of almost 60 in single pointed observations .
we will place the 13 elements over the platform to generate the longest possible baselines , which will result in angular resolutions up to @xmath55 .
the correlator is also being expanded in order to handle the larger number of cross correlations .    in this second phase of amiba operations
, the first science target will be to measure the angular power spectrum of cmb temperature anisotropies to the higher multipole numbers in order to examine the shape at and beyond the third acoustic peak ( @xmath74 ; nolta et al .
2008 ) , where data taken with the 7-element close - pack 0.6@xmath0 m configuration will be combined with upcoming 13-element 1.2@xmath0 m data .
accurate measurements of the angular power spectrum of cmb temperature anisotropies ( park et al .
2003 ) through @xmath75 up to @xmath76 will allow us to see secondary effects such as the sze ( lin et al .
2004 ) and possible cosmic string structures ( wu 2004 ) .
the second science target will be to resolve cluster sze structures on the sky in order to compare with dark matter structures as deduced from weak gravitational lensing studies ( umetsu & futamase 2000 ; broadhurst et al .
2005 , 2008 ; umetsu & broadhurst 2008 ; okabe & umetsu 2008 )
. the third science target will be to survey for the distribution of galaxy clusters via the sze ( zhang et al .
2002 ; umetsu et al .
2004 ) . to obtain redshifts of cluster candidates , optical follow up
observations will be conducted with ground - based telescopes .    at the time of the publication of this paper ,
the expansion is already in progress .
we anticipate first operations during 2009 .
in this introductory paper , we have presented the design and construction of the yuan - tseh lee amiba project .
the telescope is now in scientific operation , and the first science results are the detection and imaging of six galaxy clusters via the sunyaev zeldovich effect .
our companion papers elaborate on various aspects of our results .
chen et al .  ( 2008 ) describe the technical design of the instruments in detail .
wu et al .  ( 2008b ) describe the observations and data analysis for the first year of amiba sze data .
lin et al .  ( 2008 ) examine the system performance , and discuss the astrophysical properties of planetary calibrators .
koch et al .  ( 2008a ) describe the performance of the hexapod mount and the carbon fiber platform .
nishioka et al .  (
2008 ) examine the integrity of the amiba data and its statistical properties .
liu et al .  (
2008 ) address the issues of foreground and primary - cmb contamination in amiba sze observations .
umetsu et al .
( 2008 ) discuss the combination of amiba sze and subaru weak lensing data in order to probe the structure of dark matter and to derive the cluster gas mass fractions .
koch et al .  ( 2008b ) derive the value of the hubble constant @xmath72 from amiba and x - ray data on the clusters .
huang et al .  ( 2008 ) derive cluster scaling relations between the sze and x - ray observables .
finally , molnar et al .  (
2008 ) discuss the potential of amiba data in its 13-element configuration to constrain the intra - cluster gas distributions .
we thank the ministry of education , the national science council , and the academia sinica for their support of this project .
we thank the smithsonian astrophysical observatory for hosting the amiba project staff at the sma hilo base facility .
we thank the noaa for locating the amiba project on their site on mauna loa .
we thank the hawaiian people for allowing astronomers to work on their mountains in order to study the universe .
bennett , c. et al .  2003 ,
apjs , 148 , 1 birkinshaw , m. 1999 , phys rep , 310 , 97 , t. , takada , m. , umetsu , k. , kong , x. , arimoto , n. , chiba , m. , & futamase , t. 2005 , , 619 , l143 , t. , umetsu , k. , medezinski , e. , oguri , m. , & rephaeli , y. 2008 , apj , 685 , l9 carlstrom , j.d . , holder , g.p . , and reese , e.d .
2002 , araa . , 40 , 643 chen , m .-
2008 , , submitted coble , k. , bonamente , m. , carlstrom , j. e. , dawson , k. , hasler , n. , holzapfel , w. , joy , m. , la roque , s. , marrone , d. p. , & reese , e. d. 2007 , , 134 , 897 dunkley et al .
2008 , submitted to apjs , preprint ( arxiv:0803.0586 ) halverson , n.w .
2008 , , submitted ( arxiv : 0807.4208 ) ho , p.t.p . , moran , j.m .
, & lo , k.y . 2004 , , 616 , l1 ho , p.t.p . et al .  2008 , modern physics letter a , 23 , 1243 huang , c .- w .
l. et al .
2008 , in preparation ( apj companion paper ) koch , p. et al .
2008a , , submitted koch , p. et al .
2008b , , in preparation koch , p. et al .
2006 in _ proceedings of the european conference on antennas and propagation _ : eucap 2006 ( esa sp-626 ) , eds .
: h. lacoste & l. ouwehand . published on cdrom . , p.668.1 kogut , a. et al .  2003 , apjs , 148 , 161 komatsu , e. et al .
2008 , apjs submitted ( arxiv : astro - ph/08030547 ) kovac , j.m .
, et al .
2002 , nature , 420 , 772 leitch , e.m .
et al .  2002a , , 568 , l28 leitch , e.m .
et al .  2002b , nature , 420 , 763 li , c .- t .
2004 in _ millimeter and submillimeter detectors for astronomy ii _ , ed .
j. zmuidzinas , w. holland , and s. withington , spie , bellingham , wa liang , h. 2002 , in _
amiba 2001 : high - z clusters , missing baryons , and cmb polarization , asp conference proceedings , vol .
chen et al .
( astronomical society of pacific ) liao , y .- w .
2008 ,  ( lett . ) , submitted ( this issue ) lin , k .- y . ,
p . , tseng , y .- h . , lin , l. , & chiueh , t. 2004 , , 608 , l1 lin , k .- y .
2008 , , submitted liu , g .- c .
2008 , , in preparation lo , k.y .
2001 , in _ aip conference proceedings _ , 586 , 172 markevitch , m. et al .  2000 , , 541 , 542 mason , b.s . , et al
2003 , , 591 , 540 mather , j.c .
et al .  1990 , , 354 , l37 molnar , s. et al .
2008 , , in preparation mroczkowski , t.et al .
2008 , apj , submitted ( arxiv:0809.5077 ) nishioka , h. et al .
2008 , , in press ( arxiv:0811.1675 ) nolta , m.  r. et al .
2008 , apjs submitted ( arxiv : astro - ph/08030593 ) okabe , n. & umetsu , k. 2008 , , 60 , 345 padin , s. et al .
2002 , , 114 , 83 page , l. et al .
2007 , apjs , 170 , 335 park , c .- g . , ng . , k .- w . ,
park , c. , liu , g .- c . , & umetsu , k. 2003 , , 589 , 67 pearson , t.j .
, et al .
2003 , , 591 , 556 penzias , a.  a. & wilson , r.  w. 1965 , , 142 , 419 pryke , c. et al .
2008 , apj submitted ( arxiv : astro - ph/08051944 ) readhead , a.  c.  s. et al .
2004 , science , 306 , 836 rephaeli , y. 1995 , araa , 33 , 541 scaife , a.  m.  m. et al .
2008 , , 385 , 809 smoot , g.f .
1992 , , 396 , l1 spergel , d.n .
et al .  2003 , apjs , 148 , 175 spergel , d.n .
et al .  2007 ,
apjs , 170 , 377 sunyaev , r.a . & zeldovich , y.b .
1970 , comments astrophys
. space phys . , 2 , 66 sunyaev , r.a . &
zeldovich , y.b .
1972 , comments astrophys .
space phys . , 4 , 173 umetsu , k. & futamase , t. 2000 , , 539 , l5 umetsu , k. , chiueh , t. , lin , k .- y .
, wu , j .-
m . , & tsheng , y .- h .
2004 , monder physics letters a , 19 , 1027 umetsu , k. & broadhurst , t. 2008 , apj , 684 , 177 umetsu , k. et al .
2008 , , in press ( arxiv:0810.0969 ) weinreb , s. , lai , r. , erickson , n. , gaier , t. , wielgus , j. 1999 , _ microwave symposium digest _ , 1999 ieee mtt - s international volume 1 , 101 white , m. , carlstrom , j.e . ,
dragovan , m. , & holzapfel , w.l .
1999 , apj , 514 , 12 wu , j.  h.  p. 2004
, modern physics letters a , 19 , 1019 wu , j.  h.  p.
et al .  2008a , modern physics letters a , 23 , 1675 wu , j.  h.  p. et al .  2008b , , in press ( arxiv:0810.1015 ) zhang , p. , pen ,
, wang , b. 2002 , , 577 , 555    ll interferometer elements & 7 0.6-meter ( 13 1.2-meter ) , f/2.0 cassegrain + telescope mount & hexapod + telescope backup structure & 6-meter carbon fiber platform + primary reflector & monolithic carbon fiber + surface accuracy & 30 microns rms + secondary reflector & carbon fiber , fixed + array configuration & rings at 0.6@xmath0 m spacings + available baselines & 0.6 - 6.0 meters + operating frequencies & 94@xmath0ghz + maximum angular resolution & @xmath4 + primary beam field of view & 23@xmath77 ( 11@xmath77 ) + receiver band & 86102@xmath0ghz + number of receivers & 7 ( 13 ) dual polarization mmic ip hemt + correlator & 4-lag analog 2 @xmath78 21 ( 78 ) baselines + point source sensitivity & 63 ( 8) mjy in 1 hour on - source integration +      m antennas ( _ small solid circles _ ) on a 6@xmath0 m single platform ( _ outer solid circle _ ) .
_ right _ : sensitivity of the 7-element amiba as a function of spherical harmonic multipole @xmath79 at two frequency channels @xmath80 ( _ solid _ ) and 98 ( _ dashed _ ) ghz , shown by window functions for three baselines with different lengths ( see eq . [ 14 ] of white et al .
[ fig2],width=529 ]    m antennas co - mounted on a 6@xmath0 m platform . shown in the right panel
are the receiver packages mounted on the platform together with various electronics such as the correlator and lo / if systems .
the reflectors and receivers can be deployed at various locations on the platform in order to achieve different projected baselines .
[ fig3],title="fig:",width=264 ] m antennas co - mounted on a 6@xmath0 m platform . shown in the right panel
are the receiver packages mounted on the platform together with various electronics such as the correlator and lo / if systems .
the reflectors and receivers can be deployed at various locations on the platform in order to achieve different projected baselines .
[ fig3],title="fig:",width=264 ]        ) , revealing a strong sze decrement of about @xmath81 in the cluster center ( wu et al .
the field size shown is @xmath82 , corresponding to @xmath83mpc@xmath60 at the cluster redshift .
the black circle indicates the size of the amiba field - of - view ( @xmath84 fwhm ) , and the black filled circle at the bottom - right corner shows the size of the amiba synthesized beam ( @xmath56 fwhm ) .
the residual rms noise level in the cleaned map is @xmath85mjy ( wu et al .
overlaid are the contours ( _ white _ ) of the projected mass distribution reconstructed from subaru weak lensing data ( see umetsu et al .
2008 ; okabe & umetsu 2008 ) .
the contours are spaced from the @xmath51 noise level at intervals of @xmath86 .
shown at the bottom - left corner is the @xmath55 fwhm of the gaussian smoothing kernel ( _ white filled circle _ ) used for mass reconstruction .
the weak lensing map shows a mass subclump in the northwest region located about @xmath87 northwest from the cluster center .
a slight excess of the sze signal , extending in the northwest direction , is seen in the northwest region at @xmath62 levels .
[ fig6],width=491 ]


in the quest for new coating materials with outperforming corrosion and wear resistance properties , chromium carbide thin films are of high interest in applications for electrical contacts and decorative purposes  @xcite .
metal carbide films produced by physical vapor deposition ( pvd ) , show different microstructures , ranging from epitaxial single - crystal materials to nanocomposites with crystallites of the carbide in an amorphous carbon / amorphous carbide matrix  @xcite .
completely amorphous metal carbide coatings are also found , in both binary and ternary systems such as cr - c , fe - c , ni - c , and w - fe - c  @xcite .
it is an important methodological challenge to theoretically model such complex non - crystalline environments , in order to understand their structure and bonding properties .    in the case of _
crystalline _ materials , there is a multitude of successful combined experimental and theoretical studies in the literature , see for example refs .  @xcite . on the contrary ,
it is far from straightforward to model _ amorphous _ structures for studying materials properties .
the lack of long - range translational symmetry provides a formidable challenge for obtaining suitable model structures of the systems .
this is also the case experimentally , as the samples become _ x - ray amorphous _ with no usable diffraction peaks to determine the structure . furthermore , in the computation of xanes and exafs spectra , the absence of symmetry means that the excitation response of each atom is unique and can be very distinct . as spectral measurements are effective samplings of individual excitations of atoms in different local chemical environments , an average over a large number of sites must be considered in a theoretical approach .
while detailed explanation for the formation of non - crystalline phases during pvd in certain carbide systems remains unclear , it involves a significant quenching for conditions with neutral atoms and substrate temperatures below the melting point of the carbide , giving rise to limited surface mobility where atoms only find locally favorable sites .
it is known that amorphous films often occur in metallic carbide systems where crystalline phases exhibit structural units that have carbon atoms in a mixture of both octahedral and prismatic sites  @xcite .
previous experimental studies of amorphous cr@xmath0c@xmath1 films  @xcite suggest that the films contain two components ; an amorphous metal - rich carbide phase , and an amorphous carbon - rich phase . with increasing total carbon content ,
the amorphous carbon - rich contribution is believed to grow and incorporate the additional carbon , leaving the metal - rich phase more or less constant in carbon content  @xcite . from c
@xmath3 x - ray photoelectron spectroscopy ( xps ) measurements on the same samples as in this work , it was determined , based on quantitative estimations of carbide and a - c contributions in the spectra , that the amorphous metal carbide phase contained approximately 20 - 30 at.% of carbon .
this is in the same compositional range as e.g.  the crystalline cr@xmath7c@xmath8 and cr@xmath5c@xmath9 phases .
on the other hand , the amorphous carbon phase contained approximately pure carbon with a high amount of @xmath10-hybridized graphene - like carbon for large amounts of total carbon contents  @xcite .
however , to identify the detailed structure and the bonding of the carbide phase in these complex non - equilibrium materials  @xcite , more detailed experimental and theoretical studies are needed , in particular focusing on the surrounding of the metal atoms .    in this work ,
we investigate and compare experiments for nanocomposite cr@xmath0c@xmath1 thin films with theory for _ completely amorphous _ systems as a function of carbon concentration , together with _ ideal crystal _ structures .
x - ray absorption near - edge structure ( xanes ) and extended x - ray absorption fine structure ( exafs ) spectroscopies are ideal tools for characterizing the structure and chemical bonding in crystalline materials  @xcite . since these techniques are element specific , we anticipate that they can also be applied as efficient probes of the local environment in more complex non - crystalline materials .
in particular , we investigate the cr @xmath2 near - edge ( @xmath3 ) and exafs signals of the crc@xmath11 carbide phase in amorphous cr@xmath0c@xmath1 thin films by experiment and _ ab initio _ calculations , including metal cr for comparison . in order to model _ completely _ amorphous systems , we apply the so - called _ stochastic quenching _
( sq ) method  @xcite .
all structural models were obtained using density functional theory ( dft )  @xcite utilizing the projector augmented wave ( paw )  @xcite method as implemented in the vienna ab - initio simulation package ( vasp )  @xcite .
the generalized gradient approximation ( gga )  @xcite for electronic exchange and correlation effects was used in these simulations .
the atoms are first placed randomly in a supercell and then relaxed using dft until the force on each atom is negligible  @xcite . in this way a large set of structures is obtained , which together represent a specific concentration , at a comparatively low computational cost .
finally , to produce the theoretical x - ray absorption spectra , an average is taken over spectra computed at each cr - atom site in the amorphous structure . in order to consider possible microstructures and phases within an amorphous matrix
, we also consider crystalline crc@xmath1 structures that exist experimentally : _ _ cf__116-cr@xmath7c@xmath8 , _ _
op__16-cr@xmath9c , _ _
hp__20-cr@xmath5c@xmath9 , _ _
hp__80-cr@xmath5c@xmath9 , _ _
op__40-cr@xmath5c@xmath9 , _ _
op__20-cr@xmath9c@xmath12 , as well as the cubic b1 ( octahedral ) and hypothetical hexagonal hp2 ( trigonal prismatic ) type of crc .
all spectra are calculated by employing the computationally efficient real - space green s function code feff9  @xcite . in this way it is possible to produce the many individual spectra which are needed to represent amorphous structures .
different magnetic structures , ferromagnetic and antiferromagnetic were considered as initial configurations , although , with the exception of b1 crc , magnetic moments relaxed to zero or almost zero during the electronic structure calculation .
the structural models used for simulations of the spectra of amorphous cr@xmath13c@xmath14 , cr@xmath15c@xmath15 , and cr@xmath14c@xmath13 were obtained using the stochastic quenching approach  @xcite in combination with dft .
cubic supercells containing 250 atoms ( cr@xmath15c@xmath15 ) or 252 atoms ( cr@xmath13c@xmath14 and cr@xmath14c@xmath13 ) were used , allowing for the selected concentrations to be modeled accurately .
the stochastic quenching procedure was carried out in a similar procedure as in ref .
the atomic positions were initiated with a random distribution in the simulation box , avoiding , for numerical stability concerns in the electronic structure calculations , pair distributions shorter than 0.8  , which is about 0.6 times the calculated c - c dimer distance in vacuum .
this random structure is then quenched by gradually relaxing the internal atomic positions inside the supercell , until they reach locally stable positions with forces on each ion being smaller or in the order of 0.02 ev /  .
for each composition , five or six volumes , with 19 - 25 sq structures per volume , were used to obtain the equilibrium volume with a cubic polynomial fit .
a similar number of sq structures were then used to obtain average pair distribution functions g(r ) for each bond type and to simulate spectroscopic properties .
the calculations were initiated with individual cr atoms having either spin up or spin down magnetic polarization with equal probability and a spatially random distribution in order to simulate a paramagnetic configuration .
some cr atoms , in particular in the compositions cr@xmath14c@xmath13 , retained a magnetic moment after the sq procedure with a mean magnitude of about 0.8 @xmath16 with minimum and maximum magnitudes 0.1 and 2.2 @xmath16 , respectively . for the composition
cr@xmath15c@xmath15 the corresponding mean value is about 0.3 @xmath16 with a few of the cr moments being around 2 @xmath16 . for the cr@xmath13c@xmath14 composition
most of the local moments disappear with a mean magnitude of around 0.02 @xmath16 .
overall , the local magnetic moments are clearly less robust as compared to similar simulations for the related amorphous crn - system  @xcite .
a gaussian broadening of the calculated pair distances of 0.2    was used in the simulations .
figure 1 shows one of the supercells ( sq1 ) generated by the procedure of the sq method in order to model amorphous cr@xmath15c@xmath15 .
large blue spheres correspond to cr and smaller brown spheres are c atoms . in order to obtain the theoretical x - ray absorption spectrum for an amorphous system ,
an average of all the spectra for the different sq supercells and cr sites was performed .
the cr - atoms are represented by large blue spheres , while c - atoms are shown as smaller brown spheres .
the supercell was drawn using vesta  @xcite.,width=340 ]      in order to produce theoretical xanes and exafs spectra for the structures described above , we employ the single particle green s function technique within dft  @xcite .
all calculations were performed utilizing the computationally efficient _ ab initio _ real - space green s function ( rsgf ) approach of the feff9 code  @xcite .
for example , in the case of cr@xmath15c@xmath15 the resulting cr @xmath2-edge spectrum is an average over 22 sq structures , each in turn an average over 125 cr sites , amounting to 2750 spectra in total to represent the amorphous system .
all the potentials are calculated self - consistently and the default hedin - lundquist exchange correlation potential is used . in more detail , to obtain accurate near - edge structures considering xanes , we utilized the real - space full multiple scattering ( rsfms ) approach within feff9 . here
, the potentials were calculated for the central scattering cr atom and for the other cr and c atoms .
to take the final state effect into account , a fully screened core - hole was used for the central scattering atom , i.e.  applying the final - state rule , or core - hole approximation .
clusters with more than 550 atoms were employed in order to converge the spectral features in the different systems .
for a brief overview of different methodologies for calculating xanes see ref .
@xcite and references within . for exafs , which considers the energy region far above the near - edge ,
it is sufficient to use multiple scattering paths with no approximation for the core - hole  @xcite .
finally , the theoretical spectra were aligned with experiment by applying a rigid shift .
amorphous chromium carbide ( cr@xmath0c@xmath1 ) samples with different carbon content were made using direct current ( dc ) magnetron sputtering as described elsewhere  @xcite .
it was also found  @xcite that the samples contain about 2 - 5 at.% oxygen present .
the xanes and exafs spectra were measured at the undulator beamline i811 on the max ii ring of the max iv laboratory , lund university , sweden  @xcite .
the energy resolution at the cr @xmath3 edge of the beamline monochromator was better than 0.5 ev .
the x - ray absorption spectra were recorded in reflection mode by collecting the fluorescence yield using a pips detector from -150 ev below to 1000 ev above the cr @xmath3 absorption edge ( 5989 ev ) with 0.5 ev energy steps . to avoid self - absorption effects in the sample and bragg peaks from the substrate , the incidence angle on the sample
was varied in 0.25 degrees steps in a maximum range of @xmath173 degrees from normal incidence with a step motor .    based on the xanes fitting results ( described in the next section ) , two cr - c and two cr - cr scattering paths obtained from feff9  @xcite were included in the exafs fitting procedure with _
_ hp__20-cr@xmath18c@xmath19 as a model system using the viper software package  @xcite .
the reference energy @xmath20 was set to the pre - peak of each cr @xmath3 x - ray absorption spectrum as determined from the first peak of the derivative of each spectrum relative to the pure bcc - cr reference sample . in the fitting procedure , @xmath21
was set to 5988.6 ev ( 25% c ) , + 0.1 ev ( 42 and 47% , as compared to the 25% c sample ) , + 0.4 ev ( 67% c ) and + 0.6 ev ( 85% c ) , respectively .
the exafs functions @xmath22 were extracted from the raw absorption data after removing remaining glitches originating from substrate diffraction peaks , subsequent background subtraction , and summarizing and averaging of 15 absorption spectra . the bond distances ( r ) ,
number of neighbors ( n ) , debye - waller factors ( @xmath23 , representing the amount of disorder ) and @xmath24=reduced @xmath25 as the squared area of the residual , were determined by fitting the back - fourier - transform signal between @xmath26=0 - 12 @xmath27 obtained from the fourier - transform within @xmath28=0.7 - 3  of the first coordination shell using a hanning window function )  @xcite and a global electron reduction factor of @xmath29=0.8 .
the calculated pairs of cr - c and cr - cr 1st and 2nd scattering paths in the model were obtained using the feff9 code  @xcite .
in this section , we start by considering theoretical xanes and exafs cr @xmath2-edge spectra for the amorphous and crystal chromium carbide structures , with a comparison to measurements carried out for the thin film systems and pure cr metal .
an analysis of the xanes data was done by fitting the measured spectra with the calculations using a superposition of an amorphous phase and a smaller contribution of different crystalline phases . using the results of the xanes fitting ,
an analysis of the exafs data was carried out for the experimental data using a standard fitting procedure for crystalline materials , in order to obtain reduced ( radial ) pair distribution functions ( rdfs ) , that are compared to theoretical rdf results . finally , the trend in the structural parameters of bond lengths and coordination numbers extracted from the exafs analysis are compared and discussed with theoretical values .
-edge , a ) xanes and b ) exafs , 125 cr - sites ( gray lines ) , are exemplified for the structure sq1 in figure 1 , together with the spectra for cr@xmath15c@xmath15 over all 22 sq structures . the average of the spectra are shown as black lines.,width=340 ]    in figure 2 , the large differences of the theoretical a ) xanes and b ) exafs spectra in the case of cr@xmath15c@xmath15 are shown .
the upper parts of the two panels a ) and b ) give the cr @xmath2-edge spectra calculated over the 125 different cr - sites ( gray lines ) for the sq1 structure illustrated in fig .  1 .
the averaged spectra are shown as solid black lines .
the lower parts of the panels correspond to the averaged results obtained for each of the 22 sq structures ( gray lines ) .
finally , the resulting xanes and exafs spectra over the whole amorphous cr@xmath15c@xmath15 structure are given by the solid black lines in the lower part of each panel .
as observed in fig .  2 , the x - ray absorption fine structures are highly sensitive to the atom local environment in the amorphous matrix .
similar behavior was found at the other concentrations .
our result in fig .
2 contains an important message to theoreticians working with spectroscopical simulations of disordered materials , and in particular amorphous substances .
it is clear that in order to give even a qualitatively reasonable description of the cr xanes or exafs spectra , averaging over hundreds of atomic sites is needed .
however , for a quantitative simulation of the average signal , which is what is measured experimentally , in particular for the exafs case it might not be enough to consider just one amorphous supercell . it can be seen in the lower curve of fig .  2 that the computed average signals for some of the sq supercells deviates considerably from the average line of all the samples .
thus , it might not be reliable to simulate exafs spectra with just one supercell , say generated with molecular dynamics and containing one- or two hundred atoms , even if the average spectra of all atoms in the supercell is considered .
-edge , a ) xanes and b ) exafs , are compared for amorphous cr@xmath13c@xmath14 ( i ) , cr@xmath15c@xmath15 ( ii ) , cr@xmath14c@xmath13 ( iii ) and with the different crystal structures : bcc - cr ( 1 ) , _ _
cf__116-cr@xmath7c@xmath8 ( 2 ) , _ _ op__40-cr@xmath5c@xmath9 ( 3 ) , _ _ hp__80-cr@xmath5c@xmath9 ( 4 ) , _
_ hp__20-cr@xmath5c@xmath9 ( 5 ) , _ _ op__20-cr@xmath9c@xmath12 ( 6 ) , _
_ op__16-cr@xmath9c ( 7 ) , phases and cubic b1 [ nacl ] ( 8) and hexagonal hp2 [ wc ] ( 9 ) type crc.,width=340 ]    in figure 3 , a comparison of theoretical a ) xanes and b ) exafs spectra is made between amorphous cr@xmath0c@xmath1 and crystalline crc@xmath11 structures . in the upper parts of each panel , results for three different amorphous cr@xmath0c@xmath1 structures with concentrations @xmath30 ( blue dashed lines ) , @xmath31 ( green line ) and @xmath32 ( red dotted line ) are shown . in the lower part of each panel , the corresponding spectra for the crystals
are presented : bcc - cr , _ _
cf__116-cr@xmath7c@xmath8 , _ _
op__16-cr@xmath9c , _ _
hp__20-cr@xmath5c@xmath9 , _ _
hp__80-cr@xmath5c@xmath9 , _ _
op__40-cr@xmath5c@xmath9 , _ _
op__20-cr@xmath9c@xmath12 , structures from experiment and cubic b1 ( octahedral nacl ) and hypothetical hexagonal _ _ hp__2 ( trigonal prismatic wc ) type of crc .
xanes spectra of amorphous cr@xmath0c@xmath1 and crystal pure cr are shown .
the spectra were normalized below the edge and 100 ev above the absorption edge .
the inset shows the experimental intensities at the pre - edge peak . at the bottom
, the first derivative of the measured spectra are shown .
the inset is a close - up of the first peak ( marked with arrow ) .
the edge energies @xmath20 and the chemical shifts were determined by curve fitting of the first peak of the derivative to 5989.0 ev ( pure cr ) , 5988.6 ev ( 25% c ) , + 0.1 ev ( 42 and 47% , as compared with 25% ) , + 0.4 ev ( 67 % ) and + 0.6 ev ( 85 % ) .,width=340 ]    in figure 4 , a direct comparison is made between the cr @xmath2-edge xanes spectra from experimental measurements for different compositions , @xmath33=0.25 , 0.42 , 0.47 , 0.67 and 0.85 ( top ) , and the structures obtained by our theoretical modeling @xmath33=0.33 , 0.5 and 0.67 ( middle ) , including results for pure crystal cr .
the theoretical concentrations are chosen to capture the overall trend in the system with the variation of carbon concentration . at the bottom
, the first derivatives of the experimental spectra are shown .
the chemical shift to higher energy is a signature of higher ionization of the chemical state of the absorbing cr atoms as the carbon content increses . for the exafs spectra
, theory reproduces the general shape of the measurements . as in the case of xanes ,
the amorphous materials yield softer spectral features compared to crystalline materials .
xanes spectra at 25@xmath34 c ( red circles ) with theory for amorphous cr@xmath13c@xmath14 ( green line ) and _ _ hp__20-cr@xmath5c@xmath9 ( blue line ) .
sum ( black line ) shows the best fit , compare with table 1.,width=340 ]    @xmath35    typically , it is most fruitful to analyze the measured exafs data in the form of reduced pair distribution functions , which provide information about the average cr - cr and cr - c bond lengths , which in turn can be compared with corresponding theoretical pair distribution functions . in order to evaluate the probability of different compositions
, we have fitted each measured xanes spectrum with the theoretical fully amorphous 33 @xmath34 c structure together with each one of the crystal structure candidates as listed in table 1 .
the curve fitting analysis was done using the spectral analysis by curve fitting macro package ( spancf ) software  @xcite package . in figure 5 ,
an example is shown for the fit ( black line ) of the 25 @xmath34 c sample ( red dots ) with amorphous cr@xmath36c@xmath37 ( green line ) and _ _ hp__20-cr@xmath5c@xmath9 ( blue line ) .
the background was taken into account both by an integrated and a linear background ( dashed lines ) and the residual is shown by a thin black curve .
the square of the residual was used to compare the quality of different fittings . from the fitting analysis of the xanes data in table 1 ,
the results show that the contribution of the amorphous cr@xmath0c@xmath1 phase increases with carbon content from 40% to 90% . in all simulations
the cr@xmath5c@xmath9 , cr@xmath7c@xmath8 and cr@xmath9c@xmath12 are the most likely crystalline - like coordination while wc , pure cr , cr@xmath19c and b1 are unlikely .    by combining experimental and theoretical xanes and exafs data at the cr @xmath2-edge in figs .  2 - 5
we make several interesting observations .
both in the xanes ( fig .  4 ) and the exafs spectra , we note that the main features observed in experiment are also captured in the calculated spectra . concerning the xanes data ,
it is known that the energy position and the shape of the main peak and the pre - edge depends on the chemical state of the absorbing atom .
the pre - peak is due to a transition of a @xmath3 electron to hybridized @xmath38 and @xmath39 ( cr @xmath40 - c @xmath41 ) states . on the contrary ,
the main peak is a pure cr @xmath3 @xmath42 @xmath43 dipole transition .
the first derivative of the xanes in the bottom of fig .
4 highlights the differences between the samples .
the intensity of the pre - peak is a signature of the amount of @xmath44-@xmath45 hybridization  @xcite , which here decreases as a function of carbon content .
generally , the chemical state of the absorbing cr atom is different in cr metal , cr@xmath9c@xmath12 and cr@xmath12o@xmath9  @xcite and the spectra shift to higher energies with increasing average oxidation states that are higher than in cr metal .
note that this could also be affected by the presence of small amounts of oxygen in the samples , up to 5 at.%  @xcite , though it is difficult to precisely estimate the effect .
although the shape of the peaks and intensities are similar in all the studied cr@xmath0c@xmath1 samples ( fig .
4 ) there is a noticeable chemical energy shift towards higher energy with increasing c content .
it is also possible to distinguish these differences in a close study of the spectral fine structure .
for instance , the samples containing 42 and 47% c show sharper resolved pre - edge peaks ( fig .  4 ) . in the case of the calculated spectra for the amorphous systems ( fig .
2a ) , the individual spectrum of each and every atom is unique with a detailed fine structure of the pre - edge peak and the main line .
the crystallite systems ( fig .
3a ) present very distinct and sharp spectral features both at the main peak and the pre - edge .
however , the spread of the different crystal spectra is smaller than the rich variation found within an amorphous structure . for the amorphous cr@xmath0c@xmath1 materials ,
the averaged spectra over thousands of cr atoms show more smooth features as compared with the crystallites , in agreement with measurements .
still , the overall shape of the spectra with pre - edge and main peak is similar .
these spectral features are also seen for the pure cr metal , similar to previous results  @xcite .
c@xmath1 films are shown for the different c concentrations , with comparison to crystal pure cr .
the fit ( black curve ) is compared to the back - fourier - transform of the first coordination shell within @xmath28=0.7 - 3  and @xmath26=0 - 12 @xmath27 .
, width=340 ]    figure 6 shows experimental exafs data in @xmath46-space of bcc - cr metal and the five amorphous nanocomposite cr@xmath0c@xmath1 films obtained after normalizing and background removal using the viper software package  @xcite . for bcc - cr
, the main sharp oscillations occur in the 4 - 11 @xmath27 region . on the contrary , for the amorphous cr@xmath0c@xmath1 films ,
the main oscillations occur in the region @xmath26=2 - 7 @xmath27 , whereas for higher @xmath26-values the oscillations are damped out in increasing noise from about @xmath26=12 @xmath27 .
this behavior is characteristic of backscattering from light c atoms  @xcite .
c@xmath1 extracted from the fourier transform of the exafs s(q ) data at the cr @xmath3 absorption edge .
the vertical dotted lines indicate the relative positions of the cr - c and cr - cr peaks but the g(r ) functions represent raw data that have not been phase shifted by @xmath47 + 0.5 ev ( compare with table 2).,width=347 ]    figure 7 shows the radial distribution functions , obtained directly by fourier transforming the measured exafs structure factors in fig .  6 by the standard exafs procedure  @xcite .
the quantitative analysis was made using @xmath48 weight as is normal for conductive samples with a reduced @xmath26-range between 0 and 12 @xmath27 .
the main peaks are dominated by the short - range order of the cr - c and cr - cr orbital overlaps observed at @xmath471.5  and @xmath472.2  , in the first coordination shell , corresponding to @xmath26 @xmath49 2 - 6 @xmath27 in fig .
6 . above 2.5
 , corresponding to @xmath26 @xmath4 6 @xmath27 , a smaller peak is observed at 3 - 3.5  due to the second coordination shell .
note that the radius - values on the abscissa shown in fig .
7 represent raw data and are shown non - phase shifted to absolute values .
@xmath50    table 2 shows the results of the exafs fitting using scattering paths from the most probable structure ( _ _ hp__20-cr@xmath5c@xmath9 with a large unit cell with lattice parameters @xmath51=@xmath52=6.96  and @xmath53=4.45  ) from the xanes results in table 1 . for amorphous materials ,
it is not possible to choose a perfect model structure due to the inherent disorder , however , selecting a simpler ( and less probable ) model such as b1-crc results in less accurate fitting with larger residuals .
the phase shift ( @xmath47 + 0.5  ) is included in the fitting model for the bond lengths for easier comparison with the theoretical values .
the details of the fitting procedure are described in the experimental section iii .
the first line in table 2 presents the fitting results of the reference sample of pure bcc - cr with known coordination numbers for the first ( 8) and second ( 6 ) scattering paths .
these nearest neighbor scattering paths represent the bond lengths space diagonal from the corner of the unit cell to the center ( 2.49  ) and the edge length equal to the lattice parameter ( 2.88  ) . as expected , the amorphous samples on the following lines exhibit larger debye - waller factors than crystal bcc - cr due to a greater amount of disorder . for comparison of the coordination numbers ,
the debye - waller factors were kept fixed for the amorphous samples .
generally , we observe a trend of decreasing cr - cr coordination numbers ( n@xmath54+n@xmath55 ) from 9.5 , 7.6 , 6.7 , 5.9 to 4.6 with increasing carbon content .
with exception of the 25% c sample , we also observe that the cr - cr and cr - c radii slightly decrease as a function of carbon content .
for the cr - c coordination numbers , the decreasing trend is less obvious although as it first increases from 1.7 to 4.1 as the carbon content increases from 25% c to 42% c , probably associated with the formation of domain structures with interfaces .
thereafter , the coordination number decreases to 3.6 at 47% c , followed by a small increase to 3.7 and 3.8 as a function of carbon content . the dashed curves in fig .
7 show the individual scattering path components resulting from fitting the pairs of cr - c ( 1,2 ) and cr - cr ( 3,4 ) contributions in the first coordination shell obtained from feff9 of the _ _ h__p20-cr@xmath5c@xmath9 structure .    as observed in fig . 7
( and table 2 ) , the cr - cr peaks in the first coordination shell are not experimentally resolved from the cr - c peaks for the lowest carbon content ( 25 at% ) . as will be discussed below
, the most likely explanation is that the difference between the cr - c and cr - cr bonds becomes smaller at low c - content due to carbon vacancies in the structure .
the cr - cr bond distance is shorter in the carbide compared to the more open reference bcc - cr structure . as the carbon content increases to 42 at% , the cr - cr peak shifts to larger radius values due to the expansion of the structure .
at the same time , the cr - c peak becomes clearly distinguished at lower radius . for carbon contents above 42 at%
the cr - c and cr - cr peaks stays approximately at the same values as indicated by the vertical dotted lines .
this signifies that the carbide structure remains the same , independent of carbon content . for the second coordination shell corresponding to larger cr - cr and cr - c distances between 3.7 - 4.2  , it is interesting to notice that the relative intensity is rather independent of carbon content while the shape changes .
figure 8 shows theoretical ( radial ) pair distribution functions g(r ) obtained directly from the sq structures for amorphous cr@xmath14c@xmath13 , cr@xmath15c@xmath15 , and cr@xmath13c@xmath14 . from top to bottom , the graphs display g(r ) for cr - cr , cr - c and c - c , respectively .
the inset is zoomed in on the first cr - cr peak , illustrating the small but distinct shift to larger cr - cr distances upon increasing carbon content .
the pair distribution function g(r ) is related to the reduced pair distribution function g(r ) by g(r)=g(r)/4@xmath56@xmath57 + 1 , where @xmath57 is the average number density of the material .
the ( radial ) pair - distribution function is normalized so that g(r ) @xmath42 1 when r @xmath42 @xmath58 @xcite .
note that while the theoretical g(r ) is split into separate parts , g(r ) obtained from the measurements includes all the cr bond lengths , in particular , the cr - cr and cr - c peaks in the first coordination shell .
figure 9 shows theoretical cr - cr , cr - c and c - c g(r ) obtained from the optimized geometries of crystalline crc@xmath11 phases : _ _ cf__116-cr@xmath7c@xmath8 , _ _
op__16-cr@xmath9c , _ _
hp__20-cr@xmath5c@xmath9 , _ _
hp__80-cr@xmath5c@xmath9 , _ _
op__40-cr@xmath5c@xmath9 , _ _
op__20-cr@xmath9c@xmath12 , the cubic b1 ( octahedral ) and hp2 ( trigonal prismatic , wc ) type crc and bcc - cr .
these are the same phases whose xanes and exafs spectra are shown previously in fig .
3 . calculated average bond lengths and nearest neighbors are summarized in table 3 for the amorphous systems and crystalline phases , also including c - c bonds .    as shown by the exafs results in table 3
, we generally observe that the crystalline bond lengths are similar to both the theoretical fully amorphous structures and measurements ( table 2 ) , and the coordination numbers are closer to the measurements . considering the bond lengths , cr@xmath5c@xmath9 phases give the closest fit with the xanes measurements ( table 1 ) . from the fitting of the xanes measurements to the theoretical results , we find that the dominant contribution is amorphous crc@xmath1 carbide with a coordination that is most similar to cr@xmath5c@xmath6 carbide in all the samples , ranging from @xmath47 30 at% .
this is followed by less possible contributions of coordination similar to cr@xmath7c@xmath8 , cr@xmath9c@xmath12 and cr@xmath9c .
note that these results are in a general good agreement with the previous experimental findings by maury _
et al . _
@xcite and nygren _
et al . _
@xcite .
@xmath59    in table 3 , a very small increase of both the cr - c and cr - cr bond lengths as a function of c content is predicted in the theoretical purely amorphous structures .
this is consistent with the observed increase in crystal carbide lattice parameters when they are present in nanocrystal form in a carbon matrix  @xcite .
this trend is possibly due to the fact that the carbon - rich environment induces charge - transfer that reduces the electron concentration around the cr atoms with a following reduction of the strength of cr - cr bonding .
in addition , our calculations of pure amorphous materials predict a reduction of the coordination numbers of cr and an increase of coordination numbers of c as a function of carbon content .
the predicted reduction of the coordination number of cr is consistent with our measurements .
however , an increase in the coordination number of c is only experimentally observed between the 25% and 42% samples while at higher c contents , a weak decrease is observed .
this could be due to the fact that the size of the amorphous crc@xmath1 carbide and amorphous carbon domains changes as a function of total carbon content . as observed in transmission electron microscopy ( tem )  @xcite , domain sizes are less than 1 nm for x=0.25 , about 3 nm for x=0.47 , and about 2 nm for the higher carbon contents .
thus , the situation in amorphous nanocomposites is more complex than in pure amorphous materials as there is also a contribution from the interface grain boundaries that changes as a function of total carbon content . furthermore , for the crc@xmath1 carbide phase , there are at least two different competing effects that can give rise to a change in bond length : i ) at low c - content , the cr - rich phase has vacancies that cause an increase or an decrease of the lattice constant that has a maximum for under - stoichiometric carbide .
ii ) more carbon atoms filling into the vacancies of the structure gives rise to a larger lattice parameter @xmath51 , see e.g. , ref .  @xcite .    in a previous experimental study on the same samples  @xcite ,
@xmath33=0.25 was observed to have the most amorphous character , in high resolution transmission electron microscopy ( hrtem ) .
a possible interpretation is that the low carbon concentration sample can be described with a cr@xmath5c@xmath6 coordination with @xmath60 -0.5 , prohibited from forming the optimal structure by carbon vacancies .
the fitting with theoretical xanes spectra ( table 1 ) does not indicate a change in structure or a possible concentration of bcc - cr , while having the lowest amorphous character as compared with the other samples .
we find that in general , exafs modeling can give useful information about the local structure in amorphous materials but largely depends on the chosen model .
therefore , it is important to use _ ab initio _ theory as a complement .    while the present work focus on new theoretical and experimental results for the cr @xmath2-edge , the corresponding results for the _ carbon _ @xmath2-edge can provide valuable complementary information .
therefore , it is also of interest to compare the theoretical c - c bond for the amorphous structures and crystallites ( fig .  8 , 9 and table 3 ) with measurements . generally , the calculated cr - c and cr - cr bond lengths were found to be rather similar between the crystallites and amorphous structures . however , the theoretical c - c bond distances are strikingly different in amorphous carbides ( @xmath471.5  ) compared to crystalline carbide materials ( @xmath473  ) , which are closer to the second c - c coordination shell . as discussed in a previous work  @xcite
, these results qualitatively agree well with the relative distribution of carbon in the carbidic , crc@xmath1 phase and the amorphous carbon phase as reported in ref .
@xcite . for comparison ,
the c - c bond lengths are @xmath471.5  in a - c , graphite and diamond  @xcite .
the experimentally observed c - c bond lengths ( @xmath472  )  @xcite are larger than the theoretical , but significantly shorter than the bonds for the crc crystallites in table 3 , which are similar to previous calculations  @xcite .
note that in this study , the cr - c bond length ( @xmath47 2.0  ) was extracted from exafs modeling , while in ref .  @xcite the bond length ( @xmath47 2.2  ) was only estimated directly from the rdf obtained from electron energy loss spectroscopy with lower energy resolution .
with the combination of experiment and theory for xanes and exafs , we have shown that it is possible to obtain detailed information about the local structure in amorphous nanocomposite materials . using the computationally efficient _ ab initio _
stochastic quenching method , we modeled fully amorphous structures , taking into account local cr magnetic moments in the derivation . experimentally , amorphous structures require advanced fitting procedures where the result is very sensitive to the structure of the model .
amorphous two - phase cr@xmath0c@xmath1 nanocomposite films were investigated , consisting of a chromium - rich amorphous carbide phase and an amorphous carbon - rich phase .
we find that the short - range coordination in the cr - rich amorphous phase is similar to that of a hexagonal cr@xmath5c@xmath6 structure . at total carbon contents below 30 at% ,
the cr - rich phase is clearly under - stoichiometric with carbon vacancies that causes shorter cr - cr bond lengths and directly depend on the packing and orbital overlap with the cr - c bonds .
furthermore , the cr - cr bond length in the amorphous carbide has shorter bond lengths than the more open structure of bcc - cr metal ( space diagonal=2.49  , lattice parameter @xmath51=2.88  ) that indicates stronger bonds .
above 30 at% total carbon content , the structure of the amorphous cr@xmath5c@xmath6-like domains remains , while the excessive carbons assemble in the amorphous carbon phase .
the coordination numbers in the first coordination shell of cr decreases as a function of total carbon content from 9.5 to 4.6 , while the coordination number of c only weakly increases .
this is associated with a smaller number and size of the cr - rich domains causing a larger fraction of interface bonding that lacks cr - nearest neighbors on one side .
theoretically , we also predict that the c - c bond lengths are very sensitive to the local structure , with considerably shorter bond lengths in amorphous materials ( @xmath471.5 ) compared to various crystals ( @xmath472.8 - 3.2 ) . generally , it is important to use theory as an assessment to investigate the trends in complex amorphous materials .
we would like to thank the staff at max - lab for experimental support and u.  jansson and m.  andersson for providing the samples .
this work was supported by the swedish research council ( vr ) linnaeus grant lili - nfm , the funcase project supported swedish strategic research foundation ( ssf ) .
acknowledge financial support from vr grant no .
621 - 2011 - 4426 , the swedish government strategic research area in materials science on functional materials at linkping university ( faculty grant sfo - mat - liu no 2009 00971 ) , knut and alice wallenbergs foundation project strong field physics and new states of matter 2014 - 2019 ( cotxs ) .
b.a .  would like to thank e.  holmstrm and r.  lizrraga for support with the sq method and acknowledges financial support by the swedish research council ( vr ) through the young researcher grant no .
621 - 2011 - 4417 and the international career grant no .  330 - 2014 - 6336 and marie sklodowska curie actions , cofund , project inca 600398 .
the calculations were performed using supercomputer resources provided by the swedish national infrastructure for computing ( snic ) at the national supercomputer centre ( nsc ) and center for parallel computing ( pdc ) .
cavaleiro , a. ; de hosson , j. t. m. _ nanostructured coatings _ ; springer : new york , 2007 .
cheng , y. q. ; ma , e. atomic - level structure and structure - property relationship in metallic glasses .
_ prog . mat .
science _ * 2011 * , _ 56 _ , 379 .
furlan , a. ; lu , j. ; hultman , l. ; jansson , u. ; magnusson , m. crystallization characteristics and chemical bonding properties of nickel carbide thin film nanocomposites .
_ j. phys . :
matt . _ * 2014 * , _ 26 _ , 415501 .
magnuson , m. ; wilhelmsson , o. ; mattesini , m. ; li , s. ; ahuja , r. ; eriksson , o. ; hgberg , h. ; hultman , l. ; jansson , u. anisotropy in the electronic structure of v2gec investigated by soft x - ray emission spectroscopy and first - principles theory .
b _ * 2008 * , _ 78 _ , 035117 .
olovsson , w. ; tanaka , i. ; mizoguchi , t. ; radtke , g. ; puschnig , p. ; ambrosch - draxl , c. al l2,3 edge x - ray absorption spectra in iii - v semiconductors : many - body perturbation theory in comparison with experiment .
b _ * 2011 * , _ 83 _ , 195206 .
olovsson , w. ; weinhardt , l. ; fuchs , o. ; tanaka , i. ; puschnig , p. ; umbach , e. ; heske , c. ; draxl , c. the be k - edge in beryllium oxide and chalcogenides : soft x - ray absorption spectra from first - principles theory and experiment .
_ j. phys . :
matt . _ * 2013 * , _ 25 _ , 315501 .
kdas , k. ; andersson , m. ; holmstrm , e. ; wende , h. ; karis , o. ; urbonaite , s. ; butorin , s. m. ; nikitenko , s. ; kvashnina , k. o. ; jansson , u. ; eriksson , o. structural properties of amorphous metal carbides : theory and experiment .
_ acta materialia _ , 2012 , * 2012 * , _ 60 _ , 4720 - 4728
.    holmstrm , e. ; bock , n. ; peery , t. b. ; lizrraga , r. ; delorenzi - venneri , g. ; chisolm , e. d. ; wallace , d. c. ab initio method for locating characteristic potential - energy minima of liquids .
e _ , 2009 , * 2009 * , _ 80 _ , 051111 .
holmstrm , e. ; bock , n. ; peery , t. ; chisolm , e. ; lizrraga , r. ; delorenzi - venneri , g. ; wallace , d. structure discovery for metallic glasses using stochastic quenching .
_ phys . rev .
b _ * 2010 * , _ 82 _ , 024203 .                      lindmaa , a. ; lizrraga , r. ; holmstrm , e. ; abrikosov , i. a. ; alling , b. exchange interactions in paramagnetic amorphous and disordered crystalline crn - based systems .
b _ * 2013 * , _ 88 _ , 054414 .
pantelouris , a. ; modrow , h. ; pantelouris , m. ; hormes , j. ; reinen , d. the influence of coordination geometry and valency on the k - edge absorption near edge spectra of selected chromium compounds .
phys . _ * 2004 * , _ 300 _ , 13 - 22 .
arcon , i. ; mirtic , b. ; kodre , a. determination of valence states of chromium in calcium chromates by using x - ray absorption near - edge structure ( xanes ) spectroscopy .
ceram . _ * 1998 * , _ 81 _ , 222 - 224 .
maury , f. ; oquab , d. ; manse , j. c. ; morancho , r. structural characterization of chromium carbide coatings deposited at low temperature by low pressure chemical vapour decomposition using dicumene chromium .
coat . tech . _ * 1990 * , _ 41 _ , 51 - 61 .
nygren , k ; andersson , m. ; hgstrm , j. ; fredriksson , w. ; edstrm , k. ; nyholm , l. ; jansson , u. influence of deposition temperature and amorphous carbon on microstructure and oxidation resistance of magnetron sputtered nanocomposite cr - c films .
sci . _ * 2014 * , _ 305 _ , 143 - 153 .
darmstadt , h. ; roy , c. ; kaliaguine , s. ; xu , g. ; auger , m. ; tuel , a. ; ramaswamy , v. solid state 13c - nmr spectroscopy and xrd studies of commercial and pyrolytic carbon blacks . _ carbon _ * 2000 * , _ 38 _ , 1279 - 1287 .
li , y. ; gao , j. ; li , q.l . ; peng , m. f. ; sun , x. h. ; li , y. y. ; yuan , g. ; wen , w. j. the electronic , mechanical properties and theoretical hardness of chromium carbides by first - principles calculations . _ j. alloys comp .
_ * 2011 * , _ 509 _ , 5242 - 5249 .