Maximizing the Bandwidth from Supercontinuum Generation in Photonic Crystal Chalcogenide Fibers

Similar documents
Use of Wijsman's Theorem for the Ratio of Maximal Invariant Densities in Signal Detection Applications

Report Documentation Page

SW06 Shallow Water Acoustics Experiment Data Analysis

Modulation Instability of Spatially-Incoherent Light Beams and Pattern Formation in Incoherent Wave Systems

3. DATES COVERED (From - To) Final Report 09/28/ /31/2015

Report Documentation Page

A report (dated September 20, 2011) on. scientific research carried out under Grant: FA

Z-scan Measurement of Upconversion in Er:YAG

Final Report for AOARD grant FA Measurement of the third-order nonlinear susceptibility of graphene and its derivatives

Ocean Acoustics Turbulence Study

Closed-form and Numerical Reverberation and Propagation: Inclusion of Convergence Effects

Scattering of Internal Gravity Waves at Finite Topography

Diagonal Representation of Certain Matrices

P. Kestener and A. Arneodo. Laboratoire de Physique Ecole Normale Supérieure de Lyon 46, allée d Italie Lyon cedex 07, FRANCE

Broadband matched-field source localization in the East China Sea*

Estimation of Vertical Distributions of Water Vapor and Aerosols from Spaceborne Observations of Scattered Sunlight

Thermo-Kinetic Model of Burning for Polymeric Materials

Development and Application of Acoustic Metamaterials with Locally Resonant Microstructures

System Reliability Simulation and Optimization by Component Reliability Allocation

REGENERATION OF SPENT ADSORBENTS USING ADVANCED OXIDATION (PREPRINT)

Analysis Comparison between CFD and FEA of an Idealized Concept V- Hull Floor Configuration in Two Dimensions. Dr. Bijan Khatib-Shahidi & Rob E.

Crowd Behavior Modeling in COMBAT XXI

FRACTAL CONCEPTS AND THE ANALYSIS OF ATMOSPHERIC PROCESSES

CRS Report for Congress

7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 8. PERFORMING ORGANIZATION REPORT NUMBER

INFRARED SPECTRAL MEASUREMENTS OF SHUTTLE ENGINE FIRINGS

DIRECTIONAL WAVE SPECTRA USING NORMAL SPREADING FUNCTION

Estimation of Vertical Distributions of Water Vapor from Spaceborne Observations of Scattered Sunlight

Report Documentation Page

HIGH-POWER SOLID-STATE LASER: LETHALITY TESTING AND MODELING

LAGRANGIAN MEASUREMENTS OF EDDY CHARACTERISTICS IN THE CALIFORNIA CURRENT

Attribution Concepts for Sub-meter Resolution Ground Physics Models

REPORT DOCUMENTATION PAGE

Quantitation and Ratio Determination of Uranium Isotopes in Water and Soil Using Inductively Coupled Plasma Mass Spectrometry (ICP-MS)

High-Fidelity Computational Simulation of Nonlinear Fluid- Structure Interaction Problems

SENSORS FOR MEASURING THE VOLUME SCATTERING FUNCTION OF OCEANIC WATERS

An Examination of 3D Environmental Variability on Broadband Acoustic Propagation Near the Mid-Atlantic Bight

REPORT DOCUMENTATION PAGE. Theoretical Study on Nano-Catalyst Burn Rate. Yoshiyuki Kawazoe (Tohoku Univ) N/A AOARD UNIT APO AP

International Conference on Information Sciences, Machinery, Materials and Energy (ICISMME 2015)

High Resolution Surface Characterization from Marine Radar Measurements

Determining the Stratification of Exchange Flows in Sea Straits

Metrology Experiment for Engineering Students: Platinum Resistance Temperature Detector

Wavelet Spectral Finite Elements for Wave Propagation in Composite Plates

Shallow Water Fluctuations and Communications

Dynamics of Droplet-Droplet and Droplet-Film Collision. C. K. Law Princeton University

USMC Enlisted Endstrength Model

Grant Number: N IP To compare obtained theoretical results with NPAL experimental data.

Reverse Ion Acceleration by Laser-Matter Interaction

Predicting Tropical Cyclone Formation and Structure Change

Real-Time Environmental Information Network and Analysis System (REINAS)

Coastal Engineering Technical Note

Hurricane Wave Topography and Directional Wave Spectra in Near Real-Time

Sensitivity of West Florida Shelf Simulations to Initial and Boundary Conditions Provided by HYCOM Data-Assimilative Ocean Hindcasts

Volume 6 Water Surface Profiles

uniform distribution theory

Marginal Sea - Open Ocean Exchange

PIPS 3.0. Pamela G. Posey NRL Code 7322 Stennis Space Center, MS Phone: Fax:

Super-Parameterization of Boundary Layer Roll Vortices in Tropical Cyclone Models

Super-Parameterization of Boundary Layer Roll Vortices in Tropical Cyclone Models

Analysis of South China Sea Shelf and Basin Acoustic Transmission Data

Improved Parameterizations Of Nonlinear Four Wave Interactions For Application In Operational Wave Prediction Models

Comparative Analysis of Flood Routing Methods

Understanding Near-Surface and In-cloud Turbulent Fluxes in the Coastal Stratocumulus-topped Boundary Layers

MODELING SOLAR UV/EUV IRRADIANCE AND IONOSPHERE VARIATIONS WITH MT. WILSON MAGNETIC INDICIES

VLBA IMAGING OF SOURCES AT 24 AND 43 GHZ

Rogue Wave Statistics and Dynamics Using Large-Scale Direct Simulations

INTERACTION AND REMOTE SENSING OF SURFACE WAVES AND TURBULENCE

Range-Dependent Acoustic Propagation in Shallow Water with Elastic Bottom Effects

Analysis of Infrared Measurements of Microbreaking and Whitecaps

Flocculation, Optics and Turbulence in the Community Sediment Transport Model System: Application of OASIS Results

Direct Numerical Simulation of Aeolian Tones

Predictive Model for Archaeological Resources. Marine Corps Base Quantico, Virginia John Haynes Jesse Bellavance

THE EULER FUNCTION OF FIBONACCI AND LUCAS NUMBERS AND FACTORIALS

Periodic Magnetoresistance Oscillations in Side-Gated Quantum Dots

Swash Zone Dynamics: Modeling and Data Analysis

Optimizing Robotic Team Performance with Probabilistic Model Checking

Coastal Mixing and Optics

Catalytic Oxidation of CW Agents Using H O in 2 2 Ionic Liquids

Improvements in Modeling Radiant Emission from the Interaction Between Spacecraft Emanations and the Residual Atmosphere in LEO

Extension of the BLT Equation to Incorporate Electromagnetic Field Propagation

Weak Turbulence in Ocean Waves

Assimilation of Synthetic-Aperture Radar Data into Navy Wave Prediction Models

Long-Range Underwater Sound Propagation: Environmental Variability, Signal Stability and Signal Coherence

Internal Tide Generation in the Indonesian Seas

Sums of the Thue Morse sequence over arithmetic progressions

Computer Simulation of Sand Ripple Growth and Migration.

Satellite Observations of Surface Fronts, Currents and Winds in the Northeast South China Sea

Regional Stratification and Shear of the Various Streams Feeding the Philippine Straits ESR Component

Coupled Ocean-Atmosphere Modeling of the Coastal Zone

Models of Marginal Seas Partially Enclosed by Islands

Design for Lifecycle Cost using Time-Dependent Reliability

DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited.

Wavelets and Affine Distributions A Time-Frequency Perspective

Advanced Numerical Methods for NWP Models

Evolution of Tropical Cyclone Characteristics

Molecular Characterization and Proposed Taxonomic Placement of the Biosimulant 'BG'

Mine Burial Studies with a Large Oscillating Water-Sediment Tunnel (LOWST)

Generation and Propagation of Internal Solitary Waves on the Continental Shelf and Slope

Exact Solution of a Constrained. Optimization Problem in Thermoelectric Cooling

The Mechanics of Bubble Growth and Rise in Sediments

Transcription:

Maximizing the Bandwidth from Supercontinuum Generation in Photonic Crystal Chalcogenide Fibers Curtis R. Menyuk based on the PhD dissertation of: Dr. Jonathan Hu now at Princeton University 1

Report Documentation Page Form Approved OMB No. 0704-0188 Public reporting burden for the collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson Davis Highway, Suite 1204, Arlington VA 22202-4302. Respondents should be aware that notwithstanding any other provision of law, no person shall be subject to a penalty for failing to comply with a collection of information if it does not display a currently valid OMB control number. 1. REPORT DATE SEP 2010 2. REPORT TYPE N/A 3. DATES COVERED - 4. TITLE AND SUBTITLE Maximizing the Bandwidth from Supercontinuum Generation in Photonic Crystal Chalcogenide Fibers 5a. CONTRACT NUMBER 5b. GRANT NUMBER 5c. PROGRAM ELEMENT NUMBER 6. AUTHOR(S) 5d. PROJECT NUMBER 5e. TASK NUMBER 5f. WORK UNIT NUMBER 7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) NATO; UMBC 8. PERFORMING ORGANIZATION REPORT NUMBER 9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 10. SPONSOR/MONITOR S ACRONYM(S) 12. DISTRIBUTION/AVAILABILITY STATEMENT Approved for public release, distribution unlimited 11. SPONSOR/MONITOR S REPORT NUMBER(S) 13. SUPPLEMENTARY NOTES See also ADA564694. Mid-Infrared Fiber Lasers (Les fibres laser infrarouge moyen). RTO-MP-SET-171 14. ABSTRACT 15. SUBJECT TERMS 16. SECURITY CLASSIFICATION OF: 17. LIMITATION OF ABSTRACT SAR a. REPORT unclassified b. ABSTRACT unclassified c. THIS PAGE unclassified 18. NUMBER OF PAGES 39 19a. NAME OF RESPONSIBLE PERSON Standard Form 298 (Rev. 8-98) Prescribed by ANSI Std Z39-18

Maximizing the Bandwidth from Supercontinuum Generation in Photonic Crystal Chalcogenide Fibers With: Dr. L. Brandon Shaw, J. S. Sanghera, and I. D. Aggarwal at the Naval Research Laboratory 2

Project Goal GOAL: GOAL: To make a broadband (2 10 µm) mid-ir source 10 0 Power (a.u.) A mid-ir Light bulb Bandwidth 10 60 2 Wavelength (µm) 8 3

Why mid-ir sources? Many important materials radiate or absorb in this range 1.55 µ m 9.5 µ m 4 Spectral response of ammonia And it is not alone!

Why chalcogenide? Attenuation in silica grows rapidly beyond 2.5 µ m Attenuation in the chalcogenides remains small beyond 10 µm 5 Source: Oxford Electronics www.oxford-electronics.com

What is chalcogenide? glass is based on chalcogens mixed with As losses ~ 0.1 1 db/m Kerr nonlinearity = 1000X silica fiber CW peak power = 50 125 kw/cm 2 pulse peak power = 1 2 GW/cm 2 6

Photonic crystal fiber (PCF) Solid-core PCF Photonic bandgap fiber (PBGF) 7 solid core holey cladding forms effective low-index material Air core periodic cladding forms photonic band gap

Photonic crystal fiber (PCF) Solid-core PCF We focus on solid-core PCFs to make use of the nonlinearity solid core 8 holey cladding forms effective low-index material

Supercontinuum generation Supercontinuum generation Kerr nonlinearity Raman effect Dispersion It is a complicated, incoherent process! 9

Supercontinuum generation Supercontinuum generation Kerr nonlinearity Raman effect Dispersion Supercontinuum generation using photonic crystal fiber (PCF) 1 Wide single-mode region Enhanced nonlinearity Tailored dispersion 10 1 Dudley, et al., Rev. Mod. Phys. 78, 1135 (2006).

Supercontinuum generation Supercontinuum generation in chalcogenide fibers is not the same as in silica fibers! WHY? Different material properties There are no good sources beyond 2.5 3.0 µm Our design goal is to increase the maximum wavelength of the spectrum as rapidly as possible! 11 1 Dudley, et al., Rev. Mod. Phys. 78, 1135 (2006).

Supercontinuum generation Supercontinuum generation in chalcogenide fibers is not the same as in silica fibers! A key finding: supercontinuum generation proceeds in two stages Stage 1: four-wave mixing Stage 2: soliton self-frequency shift Each stage should be as large as possible! 12 1 Dudley, et al., Rev. Mod. Phys. 78, 1135 (2006).

Prior work Delmonte, et al. 1 show experimental supercontinuum generation from 0.9 to 2.5 µm using a tellurite fiber with a wagon-wheel structure. Price, et al. 2 theoretically demonstrate supercontinuum generation from 2 to 4 µm using a bismuth glass fiber with a wagon-wheel structure. Shaw, et al. 3 show experimental supercontinuum generation from 2.1 to 3.2 µm in a As 2 Se 3 based chalcogenide PCF with one ring of air holes. 13 1 Delmonte, et al., CLEO, CTuA4 (2006) 2 Price, et al., J. Sel. Topics Quantum Electron. 13, 738 (2007). 3 Shaw, et al., Adv. Solid State Photonics TuC5 (2005)

Model validation Shaw, et al. 1 show experimental supercontinuum generation from 2.1 to 3.2 µm in a As 2 Se 3 based chalcogenide PCF with one ring of air holes. We use the Shaw, et al. results to validate our model 14 1 Shaw, et al., Adv. Solid State Photonics, TuC5 (2005).

Design criteria Supercontinuum generation is a complicated process BUT there are general design criteria that work well 1. Design the fiber so that it is single-mode increases the effective nonlinearity 2. Ensure that four-wave mixing is phase-matched with the largest possible Stokes wavelength Rapidly moves energy to a large wavelength 3. Make the second zero dispersion wavelength as large as possible Allows the soliton self-frequency shift to go to long wavelengths 15

A specific example Fixed fiber and pulse features As 2 Se 3 fiber Five-ring hexagonal structure A pump wavelength of 2.5 µ m Fiber parameters to vary: Air-hold diameter (d ) Pitch (Λ) Pulse parameters to vary: Peak power Pulse duration Λ 16

A specific example Needed fiber quantities (experimentally determined) Kerr coefficient Raman gain Material dispersion Needed fiber quantities (calculated) Total Raman response calculated once Total dispersion calculated for each set of fiber parameters Λ 17

Generalized nonlinear Schrödinger equation (GNLS) In principle: We can optimize by solving the GNLS for a broad set of fiber and pulse parameters A( z, t) z iift ( + Ω) ( ) Ω '( ) %(, Ω) {[ β ω ] } 0 β ω0 β ω0 A z i t 2 = iγ 1 A( z, t) R( t ') A( z, t t ') dt ' + ω0 t A(z,t) : Electric field envelope β : Propagation constant = n /( ca ) : Kerr coefficient γ ω 2 0 eff R( t) = (1 f ) δ ( t) + f h ( t) R R R 18 Kerr effect Raman effect

Generalized nonlinear Schrödinger equation (GNLS) In practice: We use our design criteria to reduce the labor In any case: We must solve the GNLS for a broad enough parameter set to verify the design criteria 19

Third-order susceptibility Imaginary part 1 1 Real part Silica N N As 2 Se 3 0 0 0 1000 0 Frequency (cm 1 ) Frequency (cm 1 ) 1000 The real part is obtained by a Hilbert transform of the imaginary part (Kramers-Kronig Kronig relation) 20

Raman response function 5 h R (t) (a.u.) Silica As 2 Se 3-5 0 1500 t (fs) Chalcogenide fiber has a longer response time than silica fiber 21

Raman gain and Raman response function R( t) = (1 f ) δ ( t) + f h ( t) R R R Kerr effect Raman effect h R (t) [ ] Ω = f 0.1 g( ) (2 ω / c) n f Im FFT( h ( t)) p 2 R R R Raman gain Pump frequency Raman fraction of the nonlinearity 22 1 Stolen, et al., J. Opt. Soc. Am. B 6, 1159 (1989). 2 Slusher, et al., J. Opt. Soc. Am. B 21, 1146 (2004).

Fiber geometry Experiment 1 Simulation D d 23 Preform PCF 1 Shaw, et al., Adv. Solid State Photonics, TuC5 (2005) Λ D = 10 µm d/λ = 0.8

Dispersion 50 D (ps/nm-km) -100 2 6 Wavelength (µm) Material Fiber d Λ D = 10 µm d/λ = 0.8 24

Supercontinuum generation 10 Power (a.u.) 10 0-8 2 Wavelength (µm) 3.4 Experiment 1 Simulation L = 1 m FWHM = 100 fs Energy = 100 pj λ p = 2.5 µm Measured nonlinear response can completely account for the supercontinuum generation 25 1 Shaw, et al., Adv. Solid State Photonics, TuC5 (2005)

Single-mode analysis Step-index fiber Solid-core PCF 2π 2 2 V = a n co n cl Fundamental space-filling mode λ n FSM 26

Fundamental space-filling mode 3 n FSM 1 0 λ / Λ 10 d / Λ 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 27

Endlessly single-mode region λ /Λ 10 1 10 10 1 Endlessly Single-mode region 2 0 Single-mode region Multi-mode region n=2.8 n=2.4 n=1.45 d / Λ 0.9 28

Endlessly single-mode region What we learned: d / Λ = 0.4 When, the fiber is single mode AND We have the best mode confinement. d / Λ = 0.4 We set from this point on. 29

Four-wave mixing (FWM ) 2ω p ω s ω a Phase-matching condition ( n ω + n ω 2 n ω ) / c + 2(1 f ) γ P = 0 s s a a p p R p 30

FWM wavelength P = 0 kw P = 1 kw P = 0.1 kw Λ = 3 µm P = 0.1 kw Λ = 4 µ m Λ = 3 µ m Λ = 2 µ m At P = 0.1 kw, Λ = 3 µ m gives a large Stokes wavelength 31

5Bandwidth (µm) Bandwidth as a function of pitch 0 2 Λ (µm) 4 32

Output spectrum 10 0 Power (a.u.) 10 60 Bandwidth 20 db 1 Wavelength (µm) 8 L = 0.1 m FWHM = 500 fs d/λ= 0.4 Λ = 3 µm 33

34 Spectrogram t (ps) Distance (m) Wavelength (µm)

5Bandwidth ( µm) Bandwidth as a function of input pulse width 35 0 100 FWHM (fs) 1000

6Bandwidth (µm) Bandwidth as a function of input peak power 0 0 Peak power (kw) 1.5 36

Application of this approach to other fibers Weiblen, et al. As 2 S 3 fiber Presented at CLEO 2010 Stay tuned for more! 37

Conclusions We have developed a design approach that allows us to maximize the supercontinuum bandwidth in chalcogenide fibers We showed that a bandwidth of 4 µm can be generated using an As 2 Se 3 PCF with d/λ = 0.4 and Λ = 3 µm at a pump wavelength of 2.5 µm 38 This same approach can be applied to a wide variety of chalcogenide fibers.

39 Thank you!