To the use of Sellmeier formula

Similar documents
Ray Optics Theory and Mode Theory. Dr. Mohammad Faisal Dept. of EEE, BUET

REGRESSION (Physics 1210 Notes, Partial Modified Appendix A)

Free Space Optical Wireless Communications under Turbulence Channel Effect

Topic 9: Sampling Distributions of Estimators

Statisticians use the word population to refer the total number of (potential) observations under consideration

The axial dispersion model for tubular reactors at steady state can be described by the following equations: dc dz R n cn = 0 (1) (2) 1 d 2 c.

TR/46 OCTOBER THE ZEROS OF PARTIAL SUMS OF A MACLAURIN EXPANSION A. TALBOT

Lecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting

Stability of fractional positive nonlinear systems

Numerical Methods in Fourier Series Applications

Expectation and Variance of a random variable

Lecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting

4. Partial Sums and the Central Limit Theorem

Optimally Sparse SVMs

Comparison Study of Series Approximation. and Convergence between Chebyshev. and Legendre Series

CEE 522 Autumn Uncertainty Concepts for Geotechnical Engineering

Linear Regression Models

Topic 9: Sampling Distributions of Estimators

OPTICAL CHARACTERISTICS OF BIOPOLYMER FILMS FROM PECTIN AND GELATIN

Castiel, Supernatural, Season 6, Episode 18

a b c d e f g h Supplementary Information

Kinetics of Complex Reactions

On a Smarandache problem concerning the prime gaps

A Cobb - Douglas Function Based Index. for Human Development in Egypt

is also known as the general term of the sequence

Section 19. Dispersing Prisms

Section 19. Dispersing Prisms


62. Power series Definition 16. (Power series) Given a sequence {c n }, the series. c n x n = c 0 + c 1 x + c 2 x 2 + c 3 x 3 +

Commutativity in Permutation Groups

t distribution [34] : used to test a mean against an hypothesized value (H 0 : µ = µ 0 ) or the difference

New Version of the Rayleigh Schrödinger Perturbation Theory: Examples

Synopsis of Euler s paper. E Memoire sur la plus grande equation des planetes. (Memoir on the Maximum value of an Equation of the Planets)

INF-GEO Solutions, Geometrical Optics, Part 1

The Method of Least Squares. To understand least squares fitting of data.

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity

Optimization Methods: Linear Programming Applications Assignment Problem 1. Module 4 Lecture Notes 3. Assignment Problem

Riesz-Fischer Sequences and Lower Frame Bounds

A sequence of numbers is a function whose domain is the positive integers. We can see that the sequence

17 Phonons and conduction electrons in solids (Hiroshi Matsuoka)

Sequences A sequence of numbers is a function whose domain is the positive integers. We can see that the sequence

Analysis of Experimental Data

Topic 9: Sampling Distributions of Estimators

FIR Filter Design: Part II

On an Application of Bayesian Estimation

Analytic Theory of Probabilities

Introduction to Laser Diffraction

Exam II Covers. STA 291 Lecture 19. Exam II Next Tuesday 5-7pm Memorial Hall (Same place as exam I) Makeup Exam 7:15pm 9:15pm Location CB 234

Confidence intervals summary Conservative and approximate confidence intervals for a binomial p Examples. MATH1005 Statistics. Lecture 24. M.

The target reliability and design working life

The McClelland approximation and the distribution of -electron molecular orbital energy levels

11 Correlation and Regression

Section 1.1. Calculus: Areas And Tangents. Difference Equations to Differential Equations

ERT 318 UNIT OPERATIONS

Provläsningsexemplar / Preview TECHNICAL REPORT INTERNATIONAL SPECIAL COMMITTEE ON RADIO INTERFERENCE

TRACEABILITY SYSTEM OF ROCKWELL HARDNESS C SCALE IN JAPAN

MODELS FOR FITTING REFRACTIVE INDEX n vs. λ

Cork Institute of Technology Bachelor of Science (Honours) in Applied Physics and Instrumentation-Award - (NFQ Level 8)

Infinite Sequences and Series

Dr. Clemens Kroll. Abstract

Linear Regression Demystified

Confidence Interval for Standard Deviation of Normal Distribution with Known Coefficients of Variation

Mechanical Efficiency of Planetary Gear Trains: An Estimate

Discrete-Time Systems, LTI Systems, and Discrete-Time Convolution

Sequences, Mathematical Induction, and Recursion. CSE 2353 Discrete Computational Structures Spring 2018

Direction of Arrival Estimation Method in Underdetermined Condition Zhang Youzhi a, Li Weibo b, Wang Hanli c

Signals & Systems Chapter3

The Growth of Functions. Theoretical Supplement

Vector Quantization: a Limiting Case of EM

WHAT IS THE PROBABILITY FUNCTION FOR LARGE TSUNAMI WAVES? ABSTRACT

Hoggatt and King [lo] defined a complete sequence of natural numbers

Curve Sketching Handout #5 Topic Interpretation Rational Functions

Stability Analysis of the Euler Discretization for SIR Epidemic Model

Chapter 4. Fourier Series

Effect of Temperature Variations on Strain Response of Polymer Bragg Grating Optical Fibers

DETERMINATION OF MECHANICAL PROPERTIES OF A NON- UNIFORM BEAM USING THE MEASUREMENT OF THE EXCITED LONGITUDINAL ELASTIC VIBRATIONS.

Weighted Approximation by Videnskii and Lupas Operators

Recursive Algorithms. Recurrences. Recursive Algorithms Analysis

Section 5.5. Infinite Series: The Ratio Test

Reliability and Queueing

Using An Accelerating Method With The Trapezoidal And Mid-Point Rules To Evaluate The Double Integrals With Continuous Integrands Numerically

ECE 8527: Introduction to Machine Learning and Pattern Recognition Midterm # 1. Vaishali Amin Fall, 2015

THE KALMAN FILTER RAUL ROJAS

Fizeau s Experiment with Moving Water. New Explanation. Gennady Sokolov, Vitali Sokolov

For example suppose we divide the interval [0,2] into 5 equal subintervals of length

Quantum Annealing for Heisenberg Spin Chains

A statistical method to determine sample size to estimate characteristic value of soil parameters

multiplies all measures of center and the standard deviation and range by k, while the variance is multiplied by k 2.

6.867 Machine learning, lecture 7 (Jaakkola) 1

ANALYSIS OF EXPERIMENTAL ERRORS

Basics of Probability Theory (for Theory of Computation courses)

Sample Size Estimation in the Proportional Hazards Model for K-sample or Regression Settings Scott S. Emerson, M.D., Ph.D.

Finally, we show how to determine the moments of an impulse response based on the example of the dispersion model.

Central Limit Theorem the Meaning and the Usage

Symmetric Division Deg Energy of a Graph

Disjoint set (Union-Find)

x a x a Lecture 2 Series (See Chapter 1 in Boas)

Sequences of Definite Integrals, Factorials and Double Factorials

MDIV. Multiple divisor functions

Independence number of graphs with a prescribed number of cliques

Transcription:

To the use of Sellmeier formula by Volkmar Brücker Seior Experte Service (SES) Bo ad HfT Leipzig, Germay Abstract Based o dispersio of pure silica we proposed a geeral Sellmeier formula for various dopats like GeO 2, P 2 O 5, B 2 O 3 ad F suitable for calculatio of dispersio i glass fibers. 1 Itroductio Already i 1871 Wolfgag Sellmeier [1] foud that the refractive idex depedece o wavelegth is very similar i may trasparet ad o-trasparet media i visible ad ear ifrared rage. It ca be described by the well-kow Sellmeier formula ( ) 1 M i1 ² Ai ² ² i Eq. 1 For M= 3 Eq. 1 was amed 3-term Sellmeier equatio. That meas measurig the refractive idex of the medium at least at 6 differet wavelegths we are able to calculate 6 Sellmeier costats A 1, A 2, A 3, λ 1, λ 2 ad λ 3 ad to approximate the dispersio curve. Today Sellmeier formula is widely used i optical sciece ad optical idustry to describe ad characterize the dispersio i glasses ad crystals. Ad may compaies which offer optical materials deliver their products together with the correspodig Sellmeier coefficiets, e.g. Schott [2]. As a example we ca use Sellmeier costats of the most popular glass BK7 [2] to describe the dispersio curve. I practical cases it is ofte ecessary to kow dispersio (λ) roughly that s why we tried to develop a more geeral Sellmeier formula which ca be used for may materials like a rule of thumb. 2 Glasses for fiber optics For fibers with core ad claddig(s) we eed special glasses, which have very low losses of about 0.2 db/km (at 1550 m) combied with well-kow refractive idices. To prepare special fibers ad/or calculate their trasmissio properties the refractive idexes of core ( core ) ad claddig ( clad ) as well as the correspodig dispersio properties should be flexible ad wellkow. I fiber optics, pure silica (SiO 2 ) with its specific depedece of refractive idex o wavelegth serves as the basic material. 2.1 Dopats of glasses For a log time [3] it is kow that the refractive idex ca be chaged by additives (dopats). To icrease the refractive idex, oe has to add, e.g., GeO 2, TiO 2 or P 2 O 5, to decrease it, e.g., B 2 O 3 or fluorie (F). Followig [3] we developed empiric formulas to describe this depedece o cocetratio d (i Mol%) of dopats: for GeO 2 Ge (d) = 7 + 0.0016. d Eq. 2 for TiO 2 Ti (d) = 7 + 0.0051. d for P 2 O 5 P (d) = 7 + 0.0009. d for B 2 O 3 B (d) = 7-0.0008. d+0.00003. d 2 Eq. 3 Eq. 4 Eq. 5

for F F (d) = 7-0.004. d+0.00015. d 2 Eq. 6 Results are depicted i Fig. 1. As referece idex we used the refractive idex i SiO 2 at λ = 633 m. 1.50 TiO 2 GeO 2 P 2 O 5 B 2 O 3 F 1.40 0 5 10 15 % Fig. 1 Chages of refractive idex i glass by additives To describe the depedece of refractive idex o wavelegths for differet dopats as fuctio of dopats cocetratio we used data from [4] ad obtaied the Sellmeier parameters give i Tab. 1. Tab. 1 Sellmeier coefficiets of glasses with differet dopats [4] Material Name i Fig. 2 A 1 λ 1 A 2 λ 2 A 3 λ 3 SiO 2 Si 0.696750 0.069066 0.408218 0.115662 0.890815 9.900559 13,5% GeO 2 + 86,5% SiO 2 Ge 0.711040 0.064270 0.451885 0.129408 0.704048 9.425478 9,1% P 2 O 5 + 90,9% SiO 2 P 0.695790 0.061568 0.452497 0.119921 0.712513 8.656641 13,3% B 2 O 3 + 86,7% SiO 2 B 0.690618 0.061900 0.401996 0.123662 0.898817 9.098960 1% F + 99% SiO 2 F 0.691116 0.068227 0.399166 0.116460 0.890423 9.993707 Oe ca see from i Fig. 2 that the dispersio shapes for differet dopats are very similar ad oly shifted by a offset Δ. It should be oted, however, that other dopats, e.g. silica with 16.9% Na 2 O ad 32.5% B 2 O 3 lead to dispersio curves with other shapes. This ca be advatageously used for suitable core-claddig compositios to maage dispersio i the fiber. Ge P Si 1.44 F B 1 1.5 Fig. 2 Dispersio i glasses with differet dopats

2.2 Variatio of cocetratio of GeO 2 dopats There exist a huge umber of papers which provide the Sellmeier coefficiets ad the correspodig dispersio for various udoped ad doped glasses. May authors studied the dispersio i silica (SiO 2 ) doped with GeO 2. I Tab. 2 we give the data of [5] = source 1, ad [6] = source 2. Tab. 2 Sellmeier coefficiets of silica for various GeO 2 cocetratios GeO 2 (mole%) source Sellmeier costat A 1 A 2 A 3 λ 1 λ 2 λ 3 0.0 1 0.6961663 0.4079426 0.8974794 0.0684043 0.1162414 9.896161 0.0 2 0.6965325 0.4083099 0.8968766 0.0660932 0.1181101 9.896160 3.1 1 0.7028554 0.4146307 0.8974540 0.0727723 0.1143085 9.896161 3.5 1 0.7042038 0.4160032 0.9074049 0.0514415 0.1291600 9.896156 4.1 1 0.6867178 0.4348151 0.8965658 0.0726752 0.1151435 10.002398 5.8 1 0.7088876 0.4206803 0.8956551 0.0609053 0.1254514 9.896162 6.3 2 0.7083925 0.4203993 0.8663412 0.0853842 0.1024839 9.896175 7.0 1 0.6869829 0.4447950 0.7907351 0.0780876 0.1155184 10.436628 7.9 1 0.7136824 0.4254807 0.8964226 0.0617167 0.1270814 9.896161 8.7 2 0.7133103 0.4250904 0.863198 0.0831439 0.1079664 9.896131 11.2 2 0.7186243 0.4301997 0.8543265 0.0634539 0.1277683 9.896181 13.5 1 0.73454395 0.4271083 0.8210340 0.0869769 0.1119519 10.48654 15.0 2 0.7249180 0.4381220 0.8221368 0.0871572 0.1078145 9.896197 19.3 2 0.7347008 0.4461191 0.8081698 0.0764679 0.1246081 9.896203 I Fig. 3 dispersio curves are give separately for two sources of Sellmeier coefficiets. Icreasig refractive idex correspods to icreasig cocetratio of dopats betwee 0 19.3% GeO 2 (Fig. 3a) ad 0 13.5% GeO 2 (Fig. 3b). As wavelegth rage we cosidered the rage of fiber optical widows (0.9 1.7 µm). 19.3% 13.5% 0% 0% 1.44 1 1.2 1.4 1.6 1.44 1 1.2 1.4 1.6 (a) (b) Fig. 3 Dispersio i silica for differet GeO 2 cocetratio from sources 2 (a) ad 1 (b) As oe ca see i Fig. 3 early all dispersio curves have the same shape. The oly differece is a offset i refractive idex Δ. For example we calculated the differece Δ betwee maximum

(19.3 % i Fig. 3a, 13.5% i Fig. 3b) ad miimum (0%) cocetratios of dopats GeO 2 (Fig. 4). Δ 0.028 (a) 0.026 0.024 0.022 (b) 0.02 1 1.2 1.4 1.6 Fig. 4 Maximum differece Δ at GeO 2 cocetratios correspodig to Fig. 3a ad b Furthermore; at all wavelegths we foud a liear depedece of refractive idex o GeO 2 cocetratio calculated by Sellmeier formula. As a example, this depedece is illustrated for 1300 m i Fig. 5 separately for refractive idexes from source 1 ad source 2. A slight differece was oly foud i the slope (0.00157 /% for source 1 ad 0.00142/% for source 2). source 1 source 2 1.44 0 5 10 15 20 %GeO 2 Fig. 5 Depedece of refractive idex o cocetratio of dopats at λ = 1300 m Havig this behavior i mid we suggest a approximatio of these dispersio curves by a geeral formula. It starts from dispersio of pure silica, ad we have to add the offset i agreemet with dopats cocetratio. To this ed we have to select a referece wavelegth, where our choice is 1300 m. That meas our geeral formula is based o Sellmeier equatio of pure silica: 3 0.6961663 0.068404 ² Si ( ) 1Ai with A i 0.4079426 ad i 0.1162414 ² Eq. 7 i1 i ² 0.8974794 9.896161 The we have to add oly a offset depedig o the cocetratio of dopats (d Ge ) usig the slope of curve i Fig. 5 (we used the curve correspodig to source 2). For GeO 2 dopats we get the followig equatio: Ge (λ)= Si (λ) + 1.4145. 10-3. d Eq. 8

Calculatios usig this equatio are i very good agreemet with data described i the literature (see). It works perfectly i the wavelegth rage used i optical commuicatio techologies like DWDM (1500 1570 m). Here the maximum error is less tha 10-4. calculated (15%) source 2 for (15%) calculated (10%) 1 1.5 Fig. 6 Dispersio curves for d Ge = 10% ad 15% GeO 2 calculated by usig Eq. 8 ad curve correspodig to source 2 for d Ge = 15% For practical use it is ofte useful to use derivatives of refractive idex, i.e. group idex d g ad group velocity dispersio (GVD) resultig i the material dispersio d d g d² parameter D Mat = = -. For pure silica it is depicted i Fig. 7. Note that d c c d² D Mat (λ) is exactly the same for all dopats ad cocetratios discussed here. 1.49 g g (15% GeO 2 ) g (10% GeO 2 ) D Mat (ps/m. km) 0 g (silica) 1 1.2 1.4 1.6 50 Fig. 7 Group idex g i doped ad udoped silica ad dispersio parameter D Mat i pore silica

2.3 Variatio of cocetratios of other dopats Ufortuately, util ow here are ot available sufficiet experimetal data to make similar calculatios for other dopats (P 2 O 5, B 2 O 3 ad F). We ca use oly make use of dopats cocetratios give i Tab. 1. Agai we selected 1300 m as referece wavelegth. I this case the geeral formula is based agai o dispersio of silica ad we have to add a offset depedig o cocetratio of dopats. The slope ca be obtaied from Fig. 2 at 1300 m. Thus we get the followig formulas: for P 2 O 5 P (λ) = Si (λ) + 1.652. 10-3. d P Eq. 9 for B 2 O 3 B (λ) = Si (λ) 3.760. 10-4. d B Eq. 10 for F F (λ) = Si (λ) 4.665. 10-3. d F Eq. 11 At λ = 633 m equatios 8 11 are i rather good agreemet with results of Fig. 1 correspodig empirical equatios 2 ad 4 6. Note that group idex ad group velocity dispersio behavior is i geeral the same as depicted i Fig. 7.. 2.4 Use of geeralized formulas Ofte there is aother task: To create ad produce special (maily sigle mode) fibers oe has to select material for core ad claddig(s) to achieve a certai umerical aperture 2 core 2 clad NA (applicatio-specific fibers). For example, we assume NA = 0.14 is required i stadard sigle-mode fibers SMF-28 of Corig Glass compay. Usig equatios 8 11 we are able to calculate for this fiber clad at give values of core ad NA, or core at give values of clad ad NA, respectively. For calculatio of ecessary cocetratio of dopats to realize a certai refractive idex 1 we developed the followig formulas (ote that 1 should be larger tha Si for GeO 2 ad P 2 O 5 ad less tha Si for B 2 O 3 ad F): for GeO 2 : for P 2 O 5 : for B 2 O 3 : for F: 1( ) Si( ) dge ( ) 3 1.414510 1( ) Si( ) dp ( ) -3 1.65210 d B Si( ) 1( ) ( ) -4 3.76010 Si( ) 1( ) df( ) -3 4.66510 Eq. 12 Eq. 13 Eq. 14 Eq. 15 Examples for calculatios usig equatios 12 14 are give i Tab. 3. Here we used NA = 0.14 ad λ = 1550 m.

Tab. 3 Calculatios of core ad claddig refractive idexes Core Claddig core dopats cocetratio equatio cladd dopats coce- equ. (%) tratio (%) 14 GeO 2 4.8 1.4446 pure SiO 2 - Eq. 7 E q. 12 60 GeO 2 8.5 1.4493 GeO 2 3.7 E q. 12 Eq. 12 14 P 2 O 5 4.1 1.4446 pure SiO 2 - Eq. 7 E q. 13 60 GeO 2 8.5 1.4493 P 2 O 5 2.8 E q. 12 Eq. 13 1.4446 pure SiO 2 - Eq. 7 14 B 2 O 3 18.0 Eq. 14 1.4446 pure SiO 2 - Eq. 7 14 F 1.5 Eq. 15 This algorithm ca be used for preparatio of applicatio-specific glass fibers. At the same time oe ca use the correspodig Sellmeier formulas to describe material, waveguide ad chromatic dispersio [7]. All calculatios have bee performed by the program MathCad. Correspodig programs may be dowloaded from the same olie-plus website at Vieweg+Teuber publisher. Refereces [1] Sellmeier W.: Zur Erklärug der aborme Farbefolge im Spektrum eiiger Substaze, Aale der Physik. 143, 1871, S. 272 282 [2] Schott data sheets: www.uqgoptics.com/pdf/schott%20bk7.pdf [3] Hultzsch, H. (Editor): Optische Telekommuikatiossysteme. Gelsekirche: Damm Verlag 1996. [4] Kerste, R. Th.: Eiführug i die Optische Nachrichtetechik. Berli Heidelberg New York: Spriger Verlag 1983. [5] Huyh T. L.: Dispersio i photoic systems, Techical report MECSE-10-2004, Dept. of Electrical ad computer systems egieerig, Moash Uiversity, Clayto, Australia 2004 [6] Brücker V.: Elemets of Optical Networkig. Spriger Fachmedie Wiesbade GmbH: Vieweg+Teuber Verlag 2011 [7] Brücker V., Bylia M.S. et al.: Material, waveguide ad chromatic dispersio i differet sigle-mode fibers (i preparatio)

http://www.spriger.com/978-3-8348-1302-2