arxiv:astro-ph/ v1 21 Nov 2005

Transcription

1 Astronomy & Astrophysis manusript no. dimeq June 25, 218 (DOI: will be inserted by hand later) Dimensional Analysis applied to Spetrum Handling in Virtual Observatory Context Pedro Osuna, Jesus Salgado arxiv:astro-ph/511616v1 21 Nov 25 ESAVO, European Spae Astronomy Centre, ESA Apartado 5727, E-288 Madrid, Spain Abstrat. The handling of units in an automated way by software systems an be a umbersome proedure when the units are parsed as strings. Software systems parsing units have to take into aount extensive tables of unit names, not always idential between different standards. In the Virtual Observatory ontext, the transfer of metadata in aess protools is speially ritial for the understanding of the ontent of the data (many of whih are legay data). Driven by this issue, we present a way to handle units automatially whih is based in dimensional analysis onsiderations. Although the approah presented here an be valid to any other branh of physis, we onentrate mainly in the appliation to spetrum handling in the VO ontext, and show how the proposed solution has been implemented in the VOSpe, a tool to handle VO-ompatible spetra. Key words. Dimensional Analysis, Spetrum, Virtual Observatory, Unit, CS 1. Introdution Handling of units in spetra an be umbersome when dealing with different Flux energy densities. For instane, the onversion between F λ and is normally done by parsing the unit name strings. Despite several efforts to homogenize unit name strings (See Greisen et al.(24), Taylor (1995), George et al. (1995)), several different standards do exist, foring the parsingstring mehanism to stik to one or several of those standards. On top of this, the adaptation of already existing legay data to one or other standard might be a umbersome -and even sometimes impossible- task. The aforementioned problems led us to try and figure out a way to automatially handle units without the neessity to parse strings. Using dimensional analysis, we have devised an algorithmi way to onvert between dimensionally homogeneous quantities. As a simple example, the string-parsing algorithm needed to onvert between /m 2 /µm and erg/m 2 /s/å (both F λ ) would be substituted for an algorithm to go from 1 1 ML 1 T 3 to 1 7 ML 1 T 3, i.e., dividing by a fator of 1 3. In setion 2 we give a general mathematial formalism of the relevant parts of Dimensional Analysis tehniques thatwillbeneededforthehandlingofthisproblem.inthis setion we will follow with very slight modifiations the exellent book by Szirtes(1997). e also desribe the handling of the dimensional matrix to unveil how to extrat Send offprint requests to: P. Osuna, J. Salgado dimensional relations between different units on the same problem.in setion3weapplythe theoryto the aseoff λ and onversions. And in setion 4, we give a general algorithmi method for the onversion between different unit systems for the ase of spetra together with details about the use of this tehnique in the VOSpe, a tool to handle SSAP (Simple Spetrum Aess Protool) ompatible spetra developed at the European Spae Astronomy Centre (ESAC) of the European Spae Ageny. 2. Dimensional Analysis overview Dimensional analysis helps in the understanding of ertain problems for whih no analyti mathematial formulation exists, or for whih the mathematial formulation is too omplex. In these ases, the dimensional analysis allows to extrat ertain onlusions about the behavior of the system without the need of speifi mathematial formulae relating the different variables of the problem at hand. In physis we deal with quantities whih have ertain dimensions. These are ombinations of a redued number of basi or fundamental dimensions. These fundamental dimensions form a dimensional system. Dimensional systems an range from mono-dimensional(only a fundamental dimension is used to represent any physial quantity) to omni-dimensional (all dimensions are fundamental dimensions)(see ref Szirtes(1997) for examples). The intermediate, and mostly used, is the multidimensional system, in whih a redued set of fundamental dimensions is used. In the examples that will follow, we will adhere to the

2 2 Osuna-Salgado: Dimensional Analysis and VO Spetra most widely used system of MASS-LENGTH-TIME system with the addition of Temperature, Eletri urrent and luminosity to deal with other more omplex problems. Following a long tradition in dimensional analysis (see Maxwell(189)) will all M-L-T these fundamental dimensions. Following this onvention, a physial quantity, e.g., a Fore, would be represented as follows: [F] = [m][a] = MLT 2 where square brakets should be read as dimensions of and we will all the rightmost part the dimensional equation of the quantity under onsideration Dimensional Matrix Let F(V 1,V 2,...,V n ) = be a physial relation among a set of variable quantities V i. The dimensional equations of the different variables will be: [V 1 ] = d α11 1 d α d α1n n [V 2 ] = d α21 1 d α d α2n n... [V n ] = d αm1 1 d αm2 2...d αmn n The matrix formed with the exponents of the dimensions is alled the Dimensional Matrix: α 11 α α 1n α 21 α α 2n D =..... (1). α m1 α m2... α mn From the dimensional matrix, we an now takle the following problem: how an we find dimensional (or dimensionless) produts of the variables at hand?. In other words, how do we find the ε 1,ε 2,...,ε n exponents that solve the following equation: [V ε1 1 V ε2 2...V εn n ] = d q1 1 dq2 2...dqp p (2) where: V 1,...,V n are the variables of the problem, d 1,...,d p are the fundamental dimensions of the problem, q 1,...,q p are the given exponents (sought ombinations of fundamental dimensions; to be set equal to zero for dimensionless produts). The problem is therefore redued to solving for ε j the following system of linear equations: n D ij ε j = q i (3) j=1 It an be shown that (see Szirtes(1997)) the solution to the above is: ε j = n E jk η k (4) k=1 where: 1. E = A 1 B.. A 1 and A is a nonsingular square matrix of size the rank of the dimensional matrix, and B is formed by the rest of olumns inluded in D and not inluded in A. The sizes of the different matries are: E : (n) (n) square matrix I : (n p) (n p) identitymatrix : (n p) (p) nullmatrix A : (p) (p) square matrix B : (p) (n p) matrix A 1 B : (p) (n p) matrix the olumn matrix η k is formed by ε n p arbitrary exponents followed by q p sought dimensional exponents. The first ε n p are arbitrary due to the fat that R D rank(d) = p and therefore out of the n exponents we an only determine p, leaving n-p arbitrary. Therefore, η k is: ε 1 ε 2. η k = q 1 q 2. q p ε n p (5) The problem is thus redued to finding the matrix P result of the produt: P ij = n E ik Z kj (6) k=1 where Eik is as defined earlier and Z kj is the result of superimposing the η k olumns for the different variables: ε 11 ε ε 1p ε 21 ε ε 2p Z = ε (n p)1 ε (n p)2... ε (n p)p q 1 q 1... q 1 (7) q 2 q 2... q q p q p... q p both P and Z are of size nxp. In the partiular ase where we are seeking for dimensionless produts, all the q i will be set to zero (dimensional exponents) as a dimensionless quantity is onsidered to have dimension=1 (for a dissertation on the onveniene to all dimensionless these type of quantities, hek Szirtes(1997)).

3 2.2. Complete set of dimensional produts: Bukingham theorem The so alled Bukingham s theorem reads as follows: The number of independent dimensional produts n π whih an be omposed for a given number of variables and dimensions is: n π = n R D for dimensionlessproduts n π = n R D + 1 for dimensionalproduts where n is the number of variables and R DM is the rank of the dimensional matrix. In this ase, the sizes of the matries at would be: I : (n R D ) (n R D ) identitymatrix : (n R D ) (R D ) nullmatrix A : (R D ) (R D ) squarematrix B : (R D ) (n R D ) matrix E : (n) (n) square matrix A 1 B : (R D ) (n R D ) matrix Osuna-Salgado: Dimensional Analysis and VO Spetra 3 (8) An obvious advantage of finding the dimensionless produts of the problem under investigation is that it redues the number of parameters relevant to the system. In the ase that the preise mathematial formula governing the behavior of the system is known, e.g., the Navier-Stokes equation for a fluid system, then the dimensionless produts give information as well on the relative importane of eah of the terms in the equation, helping in the simplifiation of the equations for ertain ombinations of the dimensionless produts (for a beautiful example regarding the Navier-Stokes equations, see Diez-Rohe(198)). 3. Appliation to Spetral Fluxes Using the dimensional analysis priniples, we will extrat the dimensionless produts for flux densities. These produts will be used later, in setion 4, to ahieve the unit onversion between flux densities in different units F λ versus λ To onstrut the dimensional matrix for F λ we will onsider a spetral energy distribution of the form: F λ (,λ,h) = (9) The dimensional equations of the different members of the previous relation are: [F λ ] = ML 1 T 3, [] = LT 1, [λ] = L, [h] = ML 2 T 1 (1) we now onstrut the following table: F λ λ h M 1 1 L T (11) aording to previous setions, we will have the following set of matries: A = B = A 1 B = E = Z = P = E Z = where following the general pratie(see Szirtes(1997)) we have hosen the matrix A to be omposed of the independent variables in relation (9), and he matrix P gives the exponents(olumns) on the dimensions (rows) that give rise to dimensionless produts. Aording to the theory before, the number of dimensionless produts in this ase would be n R DM, i.e., 4-3=1, hene the fat that three of the olumns are linear ombinations of the remaining one (idential, in this ase) in the P matrix only giving result to one dimensionless produt ombination: P = where rows identify the variables in (1) and olumns orrespond to the exponents of those variables. Therefore, the unique dimensional produt thus for this ase would be: π = F λ λ 5 2 h 3.2. versus λ (12) Following exatly the same proedure as before, we would have: λ h M 1 1 L T and therefore: A = B = (13)

4 4 Osuna-Salgado: Dimensional Analysis and VO Spetra A 1 B = E = Z = P = E Z = And therefore the only dimensionless produt would be: π = λ 3 h (14) Both dimensionless quantities in (12) and (14) are desriptions of the same physial problem. As they are lineal in the dependent variable (F λ and respetively) we are able to onlude that they will be equivalent, i.e., that we an write in this ase: π = π and therefore: λ 3 λ 5 F λ 2 h = h = F λ = λ 2 (15) whih is the physial result expeted when transforming flux densities. The dimensionless produt obtained when repeating the above proedure for and F λ with ν as the independent variable give: 2 π = ν 3 h and π = F λ ν 5 h respetively 3 This working example shows how the dimensional analysis an be used in the handling of spetrum fluxes and unit transformation. Muh more omplex problems an be takled using this approah. For a nie ompilation of literally hundreds of examples, onsult the book by Szirtes(1997). In what follows, we give an algorithm to bring the aforementioned ideas to pratie when designing a lient tool to handle Virtual Observatory spetra Comment on apparent veloity as X-axis spetral oordinate In the paper by Greisen et al.(24), the apparent radial veloity is onsidered as one of the possible spetral oordinates. A whole set of possible transformations between the different spetral oordinates ( x-axes ) is given together with their derivatives. The paper deals spetral oordinate transformations, rather than with flux transformations (on the y-axis ). In this work, we deal with energy densities whih depend, generally, on λ or ν. Conversions in the x-axis to veloity spae do need a entral referene wavelength (frequeny) whih will have to aompany the metadata for the data. One this value is known, the transformation of the x-axis will just onsist of a translationplus a dilation, and the Y-axis will be kept as is. In the ase that the data are oming with veloity in the x-axis, they will have to ontain the referene lambda giving rise to those veloity values, and the proess to onvert the data to wavelength values would be simply inverted. 4. Algorithmi approah and use in VOSpe In the Virtual Observatory ontext, there is a need to have a standard protool to make spetra aessible from different projets in a simple way. The idea is to reate a protool for spetra in the same line as the already standard Simple Image Aess Protool (SIAP), that would allow for the reation of on-the-fly Spetral Energy Distributions (SEDs) from heterogeneous data sets. This protool has been alled Simple Spetrum Aess Protool(SSAP) and basially implements a two-step proess: In the first step, a one searh is done on available servies and the math results are sent, together with metadata, in a VO standard VOTable. In the seond step, the pointers to the real data files (spetra) are alled and data are retrieved. The main problem faed when trying to reatean SED using data from different projets/formats is to ompare data different units. e should be able to transform spetral oordinate units and fluxes to a ommon unit system. Our proposed solution is to speify in the metadata (first SSAP step) for every spetrum, the dimensional equation for the spetrum axes, and use dimensional analysis to extrat the onversion formulae needed to go from one to the other. To prove that this on-the-fly onversion was possible, we developed a tool alled VOSpe, able to request different SSAP servers and produe a ommon SED from different spetra in different formats from different projets. The appliation is already available to the general publi at: and it has been used for the AVO demo whih took plae on Jan 25 & 26, 25 at ESAC. See AVO Demo(25) for referene. At the time of writing this paper, spetra from ISO, IUE, HST (FOS), SDSS, HyperLeda and FUSE projets are already providing SSAP aess to their data. All the spetra from these servies an be superimposed in the same display, and the user an generate on-the-fly SEDs as an be seen in Figure 1. Spetra from projets that don t have SSAP servies, an be loaded loally using the SSAP rapper Creator integrated in VOSpe.

5 Osuna-Salgado: Dimensional Analysis and VO Spetra 5 Fig.1. SED generated for the quasar 3273 using VOSpe, overing from Radio to X-Ray spetral data. From left to right: Radio and sub-milimiter points (loal files loaded data taken from the Geneve University database). Infrared spetrum data from the ISO SSAP servie (LS and PHOT data). Optial and ultraviolet spetrum data from FUSE (SSAP servie) and HST(FOS) (SSAP servie). On the right, X-Ray data from XMM-Newton (proessed loal file) 4.1. Algorithmi approah Thereareseveralwaystodesribe aspetrum flux andthe spetral oordinate inside a 1-D spetrum. A onversion table ould be used to make the transformations aordingly, but it is not easy to use in an automati algorithm, and this would limit the number of possible transformations allowed by the system. In this setion, we will desribe an algorithmi way to approah the units problem from the dimensional point of viewandshowhowthis isused inthe VOSpeappliation. Suppose now, we want to onvert one spetrum point [ defined by the pair [ν(hz), (Jy)] to other ones, e.g. λ(µm),f λ ( m 2 µm ], ) i.e.: λ(µm) = 1 6 L F λ ( m 2 µm ) = 11 ML 1 T 3 First we need to generate the matrix E for the original system, as we saw in setion 3: A unit an be desribed in the following way: [UNIT] = Saling M a L b T (16) wherethesalingfatorisdefinedwith respettoaertain ommon system of units and the exponents a,b, define the unit dimensionally. e will hoose the SI as our base dimensional system of referene. In order to understand how the algorithm works, suppose we are dealing with a spetrum in Jansky (y-axis) and Hertz (x-axis). The units dimensions (and sale fators) turn out to be: ν(hz) = T 1 (Jy) = 1 26 MT 2 ν h M 1 1 L 1 2 T so, in this ase: A = A 1 = and onstruting B from the Flux density units: B = 1 A 1 B = (17) where 1 and 1 26 are the referene saling fators to the SI units system for Hz and Jy respetively. so finally, onstruting the matrix E using the rules desribed in setion 3:

6 6 Osuna-Salgado: Dimensional Analysis and VO Spetra 1 E = Now, to go from one system to the other, we have to generate two different Z matries, one for the spetrum oordinate transformation and the seond one for the flux transformation. As we are looking for ertain target dimensions, the Z matrix will have to inlude the sought dimensional exponents, and therefore, as the final spetrum oordinates are µm (L) we will have: Z 1 = 1 M L T where the first zero is imposing no dependene in the flux, and the rest are the exponents in M, L, T respetively. For the final system of units, the F λ is (ML 1 T 3 ): 1 Z 2 = M L T m 2 µm where, as it was defined in previous setion, the first element is the dependene in the flux and the rest of the elements in the Z are the dimensional equation exponents for M,L and T respetively. If we multiply then the original E matrix with these Z vetors we obtain for the spetral oordinate: P 1 = E Z 1 = 1 1 and for the flux: 1 P 2 = E Z 2 = 1 2 [ ν h ν h That means, respetively: m 2 µm [λ] = [] [ν] ] = [Jy]. [ν]2 [] To finalize the transformation, we must inlude the saling fators of the different units. To ahieve this goal, we note that every time a magnitude is used, the saling must appear, i.e., S ( S (µm) λ(µm) = S (Hz) ν(hz) )F λ( m 2 µm m 2 µm ) = S (Jy) (Jy). (S (Hz)ν(Hz)) 2 here the S elements orrespond to the salings with respet to a ommon system of referene units (in this ase SI). Finally we obtain: λ(µm) = 1. (m/s) S (µm).s (Hz) ν(hz) = 1 6. (m/s) ν(hz) F λ ( m 2 µm ) = S (Jy).S(Hz) 2 (Jy) ν2 (Hz) S ( m 2 ) (m/s) µm = (Jy) ν2 (Hz) (m/s) = 1 36 (Jy) ν2 (Hz) (m/s) These final formulae tell us how to express the values of a [x,y] point in the final units as a funtion of the point in the original ones. The algorithm an be summarized as follows: 1. Construt the matrix A, using the spetral oordinate, and h dimensional equations. 2. Construt the vetor B using the flux density dimensional equation. 3. Invert the matrix A, and onstrut the matrix E as it was desribed in point Construt two Z vetors using spetral oordinate and flux density dimensional equations of the final units. 5. Multiply the matrixewith the twozvetorsto obtain the onversion fators. 6. Finally, use the saling fators to finish the onversion. 5. Conlusions e have shown how to make use of Dimensional Analysis tehniques to handle unit onversion in an automated way for the ase of spetral flux densities. e have proposed the IVOA to inlude the SCALEQ and DIMEQ parameters as part of the Simple Spetrum Aess Protool so that lients an use Dimensional Analysis algorithms to handle units automatially. The approah shown is not only relevant to the Spetral aess within the VO, and an be extended -as

7 mentioned in the introdution- to any physial problem by inluding the relevant dimensions. In this sense,weimagineanyunit asomposedofthree main attributes: 1. A name (and possibly a symbol) 2. A SCALEQ(giving the Saling of the unit with respet to the International System of Units) 3. A DIMEQ (giving the dimensions of the unit) An example of a serialization of this units model ould be the FITS representation, in whih a unit like the Jansky ould be represented in the following way: Osuna-Salgado: Dimensional Analysis and VO Spetra 7 CUNIT = Jansky CSCALEQ = 1E 26 CDIMEQ = MT 2 Certainly, to over all possible physial units, more dimensions than M-L-T would be needed, but the set will always be far smaller than the units they an represent. Aknowledgements. e aknowledge Jose Tomas Diez-Rohe, Matteo Guainazzi and Andy Pollok for their useful omments and disussions. Referenes [1] E.. Greisen, M.R. Calabretta, F.G. Valdes, S.L.Allen: Representations of Spetral Coordinates in FITS [2] Barry N. Taylor: Guide for the Use of the International System of Units (SI), NIST Speial Publiation 811 [3] I.M. George, L. Angelini: Speifiation of Physial Units within OGIP FITS files, OGIP Memo OGIP/93-1 [4] Thomas Szirtes: Applied Dimensional Analysis and Modelling, MGraw Hill 1997 [5] Jose Tomas Diez-Rohe: Fluid Mehanis, ETSIN publiations [6] James Clerk Maxwell: Eletriity and Magnetism [7] AVO Demo 25 European Spae Astronomy Centre (ESAC) (Madrid, Spain) ESA [8] International Virtual Observatory Alliane (IVOA)