![Page 1: Computational Photonics: Modeling Ultrashort Pulse Propagation …math.arizona.edu/~nrw/Spring_School_2005/moloney.pdf · White Light Continuum The Physics • White light continuum](https://reader034.vdocument.in/reader034/viewer/2022042802/5f397991417cc81eee013cf9/html5/thumbnails/1.jpg)
Computational Photonics: Modeling Ultrashort Pulse
Propagation and Iteraction with Nanophotonics Structures.
J.V. MoloneyArizona Center for Mathematical Sciences
andOptical Sciences Center
Support: AFOSR F49620-02-1-0380; F49620-03-1-0194FA9550-04-1-0213; FA9550-04-1-0355
NSF DMS0335101
![Page 2: Computational Photonics: Modeling Ultrashort Pulse Propagation …math.arizona.edu/~nrw/Spring_School_2005/moloney.pdf · White Light Continuum The Physics • White light continuum](https://reader034.vdocument.in/reader034/viewer/2022042802/5f397991417cc81eee013cf9/html5/thumbnails/2.jpg)
ACMS Research Team
• Moysey Brio (ACMS/Mathematics)• Miro Kolesik (ACMS/OSC)• Joerg Hader (ACMS/OSC)• Armis Zakharian (ACMS/OSC)• Arash Mafi (ACMS/OSC)• Walter Hoyer (ACMS/Mathematics)• Jens Foerstner (ACMS/Mathematics)• Colm Dineen (ACMS/Mathematics)
• Patrick Kano(Mathematics)• Krishna Mohan Gundu (ACMS/OSC)• Marcelo Matus (ACMS/ECE)• Yong Xie (ACMS/OSC)• Hongbo Li (ACMS/OSC)
Collaborators:
Stephan Koch (Marburg/OSC) + many more from OSC!!
![Page 3: Computational Photonics: Modeling Ultrashort Pulse Propagation …math.arizona.edu/~nrw/Spring_School_2005/moloney.pdf · White Light Continuum The Physics • White light continuum](https://reader034.vdocument.in/reader034/viewer/2022042802/5f397991417cc81eee013cf9/html5/thumbnails/3.jpg)
ACMS Photonics Supercomputing Laboratory
Dual Opteron Workstations
Infiniband
2 TB GFS Home1 TB XFS Home
46 CPU AMD Opteron Cluster 32 CPU Itanium 2 Altix 3700
2 TB PVFS Scratch Directory
Graphic Head
Dual Graphic Headto Prism Graphics
Gigabit Network
![Page 4: Computational Photonics: Modeling Ultrashort Pulse Propagation …math.arizona.edu/~nrw/Spring_School_2005/moloney.pdf · White Light Continuum The Physics • White light continuum](https://reader034.vdocument.in/reader034/viewer/2022042802/5f397991417cc81eee013cf9/html5/thumbnails/4.jpg)
Talk Outline• Femtosecond Atmospheric Light Strings
• Beyond the Nonlinear Schrödinger Equation
• Role of nanostructures in high power SCL design
• Computational Nanophotonics: Issues and Challenges
• Nanophotonics and Plasmonics Applications
• Summary
![Page 5: Computational Photonics: Modeling Ultrashort Pulse Propagation …math.arizona.edu/~nrw/Spring_School_2005/moloney.pdf · White Light Continuum The Physics • White light continuum](https://reader034.vdocument.in/reader034/viewer/2022042802/5f397991417cc81eee013cf9/html5/thumbnails/5.jpg)
Maxwell’s Equations
;
0, 0
B DE Ht t
D B
→ →→ →∂ ∂
∇× = − ∇× =∂ ∂
∇• = ∇• =r r
0 0,B H D E Pµ ε→ → → →
= = +r
Constitutive Relations:
Semiconductor passive/active media – rigorous microscopic many-body modelOther material systems - phenomenological models (Lorenz, Drüde, 2-level, .)
Challenges - reliable physical models for ultra-broadband excitations- numerical algorithms to implement them
Link between E.M and material excitation
![Page 6: Computational Photonics: Modeling Ultrashort Pulse Propagation …math.arizona.edu/~nrw/Spring_School_2005/moloney.pdf · White Light Continuum The Physics • White light continuum](https://reader034.vdocument.in/reader034/viewer/2022042802/5f397991417cc81eee013cf9/html5/thumbnails/6.jpg)
Atmospheric Femtosecond Probe
White Light Continuum
The Physics
• White light continuum spectroscopic probe
• Optical breakdown creates narrow plasma filaments - RF emission/lightning control.
Wöste et al., Laser und Optoelektronik, Vol. 29, p51 1997.
New Scientist, February 19, 2000Arizo
na C
enter
for M
athem
atica
l Scie
nces
![Page 7: Computational Photonics: Modeling Ultrashort Pulse Propagation …math.arizona.edu/~nrw/Spring_School_2005/moloney.pdf · White Light Continuum The Physics • White light continuum](https://reader034.vdocument.in/reader034/viewer/2022042802/5f397991417cc81eee013cf9/html5/thumbnails/7.jpg)
Light String Physics
• Nonlinear self-focusing in air – Pth=7 GW
• Extreme self-phase modulation - remote white light supercontinuum spectroscopic source.
• Dilute plasma channel generation – THz source,remote LIBS spectroscopy.
• Light string diameters below inner turbulence scale – obscurant penetration
![Page 8: Computational Photonics: Modeling Ultrashort Pulse Propagation …math.arizona.edu/~nrw/Spring_School_2005/moloney.pdf · White Light Continuum The Physics • White light continuum](https://reader034.vdocument.in/reader034/viewer/2022042802/5f397991417cc81eee013cf9/html5/thumbnails/8.jpg)
Nonlinear Spatial Replenishment- Low power, radial symmetry
- Nonlinear spatial replenishment (Mlejnek et al. Optics Letters 1998)
![Page 9: Computational Photonics: Modeling Ultrashort Pulse Propagation …math.arizona.edu/~nrw/Spring_School_2005/moloney.pdf · White Light Continuum The Physics • White light continuum](https://reader034.vdocument.in/reader034/viewer/2022042802/5f397991417cc81eee013cf9/html5/thumbnails/9.jpg)
Extended 3D NLS Model
( ) ( )
∫∞−
− −+−
+−−+∂∂
−∇=∂∂
tR
NN
RT
AdtttAtRknifAA
AiAAknfitAkiA
ki
zA
'2''2
22)(
222
2''2
|)(|)(||2
12
||122
β
ρωτσ
22)(
22
||||1 ρω
βρσρ aN
AAEnt
NN
gb−+=
∂∂
h
Diffraction GVD Kerr Nonlinearity Plasma Absorption/refraction
Multi-photon absorption Delayed Raman Response
- Plasma Drude Model
Avalanche generation Multi-photon generation Plasma recombination
Arizo
na C
enter
for M
athem
atica
l Scie
nces
![Page 10: Computational Photonics: Modeling Ultrashort Pulse Propagation …math.arizona.edu/~nrw/Spring_School_2005/moloney.pdf · White Light Continuum The Physics • White light continuum](https://reader034.vdocument.in/reader034/viewer/2022042802/5f397991417cc81eee013cf9/html5/thumbnails/10.jpg)
Rigorous Carrier Resolved 3D Propagator
Arizo
na C
enter
for M
athem
atica
l Scie
nces • Captures all relevant physics of vector Maxwell in absence
of backward generated field.
• Measured absorption/refractive index spectra over relevant bandwidth can be input directly to code
• Includes extreme self-focusing down to the wavelength of light in the material.
![Page 11: Computational Photonics: Modeling Ultrashort Pulse Propagation …math.arizona.edu/~nrw/Spring_School_2005/moloney.pdf · White Light Continuum The Physics • White light continuum](https://reader034.vdocument.in/reader034/viewer/2022042802/5f397991417cc81eee013cf9/html5/thumbnails/11.jpg)
Unidirectional Pulse Propagation Equation Unidirectional Pulse Propagation Equation
Equation written in spectral domain.Not a PDE in real-space representation
Linear propagation,Contains space-time focusing “terms”
Nonlinear response couplingContribution from self-steepening
div .E - related termNonzero due to gradients.
Nonlinear polarization of medium calculatedfrom material equation (Kerr effect,plasma,...)
Arizo
na C
enter
for M
athem
atica
l Scie
nces
This equation is exact as long as D is full field
To close the system of equations for numerical simulations, we approximate:
Nonlinear Dispersion unknown!
![Page 12: Computational Photonics: Modeling Ultrashort Pulse Propagation …math.arizona.edu/~nrw/Spring_School_2005/moloney.pdf · White Light Continuum The Physics • White light continuum](https://reader034.vdocument.in/reader034/viewer/2022042802/5f397991417cc81eee013cf9/html5/thumbnails/12.jpg)
Shock regularization in femtosecond pulse propagation25-fs, 8GW, 775nm light pulse propagating in air (animation records 60cm propagation):
•shock forms after 0.5m in the trailing edge of the pulse
•steep wave-form generates higher-frequency field components
Electric induction(fully resolved by algorithm)
On-axis intensity
![Page 13: Computational Photonics: Modeling Ultrashort Pulse Propagation …math.arizona.edu/~nrw/Spring_School_2005/moloney.pdf · White Light Continuum The Physics • White light continuum](https://reader034.vdocument.in/reader034/viewer/2022042802/5f397991417cc81eee013cf9/html5/thumbnails/13.jpg)
Shock regularization in femtosecond pulse propagationdetail: “shock” portion of the pulse
• Makes possible 3D carrier-resolved propagation over tens of meters!
• Describes fs pulse propagation with focusing to λ• Other “envelope models” derived in seamless fashion
![Page 14: Computational Photonics: Modeling Ultrashort Pulse Propagation …math.arizona.edu/~nrw/Spring_School_2005/moloney.pdf · White Light Continuum The Physics • White light continuum](https://reader034.vdocument.in/reader034/viewer/2022042802/5f397991417cc81eee013cf9/html5/thumbnails/14.jpg)
band gap/stop band
Periodic Hermitian Eigenproblems in 1d
ε1 ε2 ε1 ε2 ε1 ε2 ε1 ε2 ε1 ε2 ε1 ε2
ε(x) = ε(x+a)a
k is periodic:k + 2π/a equivalent to k
k
ω
0 π/a–π/a
irreducible Brillouin zone
Bragg Grating
![Page 15: Computational Photonics: Modeling Ultrashort Pulse Propagation …math.arizona.edu/~nrw/Spring_School_2005/moloney.pdf · White Light Continuum The Physics • White light continuum](https://reader034.vdocument.in/reader034/viewer/2022042802/5f397991417cc81eee013cf9/html5/thumbnails/15.jpg)
Slowing Down Light- engineering material properties!
cons
erve
d fr
eque
ncyω
conserved wavevector k
photonic band gap
band diagram (dispersion relation)
dω/dk → 0: slow light(e.g. bandedge lasing)
backwards slope:negative refraction
strong curvature:super-prisms, …
(+ negative refraction)
![Page 16: Computational Photonics: Modeling Ultrashort Pulse Propagation …math.arizona.edu/~nrw/Spring_School_2005/moloney.pdf · White Light Continuum The Physics • White light continuum](https://reader034.vdocument.in/reader034/viewer/2022042802/5f397991417cc81eee013cf9/html5/thumbnails/16.jpg)
Optically-Pumped Semiconductor Laser
λ/4 Mirror Stack
InGaAs MQW(14-16)
External Mirror Stack
Heat sinks
AR Coating
Output
MM Pump
∼ 5-6 µm!
Very efficient heat extraction!
Highly efficient incoherent diode bar to coherent output converter!
Optional heatspreader
- OSC Expts. M. Fallahi, Li Fan, Marc Schillgalies
![Page 17: Computational Photonics: Modeling Ultrashort Pulse Propagation …math.arizona.edu/~nrw/Spring_School_2005/moloney.pdf · White Light Continuum The Physics • White light continuum](https://reader034.vdocument.in/reader034/viewer/2022042802/5f397991417cc81eee013cf9/html5/thumbnails/17.jpg)
High Power OPSL Chips
US Patent # 6298076
2 mm
![Page 18: Computational Photonics: Modeling Ultrashort Pulse Propagation …math.arizona.edu/~nrw/Spring_School_2005/moloney.pdf · White Light Continuum The Physics • White light continuum](https://reader034.vdocument.in/reader034/viewer/2022042802/5f397991417cc81eee013cf9/html5/thumbnails/18.jpg)
Record High Brightness OPSL
TEM00 460 nm (20 °C)
TEM00 488 nm (10 °C)
MM 920 nm (20 °C)
MM 980 nm (10 °C)
0 10 20 30 40 50 60 70 800
5
10
15
20
25
30
35
40
Net Pump Power (W)
Out
put P
ower
(W)
Can one achieve > 1kW from a single chip?
- Pump spot 500 – 600 µm.
![Page 19: Computational Photonics: Modeling Ultrashort Pulse Propagation …math.arizona.edu/~nrw/Spring_School_2005/moloney.pdf · White Light Continuum The Physics • White light continuum](https://reader034.vdocument.in/reader034/viewer/2022042802/5f397991417cc81eee013cf9/html5/thumbnails/19.jpg)
Semiconductor Sub-CavityActive Mirror
Nanostructured Sub-cavity
- A. Zakharian, J. Hader(ACMS)
Active Mirror Reflectance
![Page 20: Computational Photonics: Modeling Ultrashort Pulse Propagation …math.arizona.edu/~nrw/Spring_School_2005/moloney.pdf · White Light Continuum The Physics • White light continuum](https://reader034.vdocument.in/reader034/viewer/2022042802/5f397991417cc81eee013cf9/html5/thumbnails/20.jpg)
VECSEL/OPSL Modeling
Optical field spectral domainpropagator in the external cavity
Optical field bidirectional Bi_BPMpropagator and carrier density rate eqnin active semiconductor sub-cavity
Thermal transport from active region through to heat sink
- A. Zakharian, M. Kolesik
![Page 21: Computational Photonics: Modeling Ultrashort Pulse Propagation …math.arizona.edu/~nrw/Spring_School_2005/moloney.pdf · White Light Continuum The Physics • White light continuum](https://reader034.vdocument.in/reader034/viewer/2022042802/5f397991417cc81eee013cf9/html5/thumbnails/21.jpg)
High Brightness OPSL Design Cycle
Sub-cavity optimization for pump/signalOptimization of MQW Semiconductor Epi Structure
PL-spectra(wafe diagnostic)
Gain Spectra(running laserdiagnositc)
3D Thermal Analysis of Optically-Pumped VCSELDevice growth and processing
Send for growth
- collaboration involving UofA (ACMS/OSC), AFRL- M. Fallahi, Li Fan, Marc Schillgalies, S.W. Koch, W. Stolz, T. Nelson, R. Bedford
![Page 22: Computational Photonics: Modeling Ultrashort Pulse Propagation …math.arizona.edu/~nrw/Spring_School_2005/moloney.pdf · White Light Continuum The Physics • White light continuum](https://reader034.vdocument.in/reader034/viewer/2022042802/5f397991417cc81eee013cf9/html5/thumbnails/22.jpg)
Sub-cavity Optimization for Bandedge Lasing
Nontrivial nonlinear optimization problem for pump absorption and signal gain!
Conventional sub-cavity design
Band-edge lasing sub-cavity design Original reflectance spectrum
Sharper and enhancedBGE reflectance spectrum
Reflectance as seen by external cavity
![Page 23: Computational Photonics: Modeling Ultrashort Pulse Propagation …math.arizona.edu/~nrw/Spring_School_2005/moloney.pdf · White Light Continuum The Physics • White light continuum](https://reader034.vdocument.in/reader034/viewer/2022042802/5f397991417cc81eee013cf9/html5/thumbnails/23.jpg)
Actively Cooled kW OPSL
Active Structure+ 60°
Active Structureblowup
Transparent singlecrystal diamond heatspreader
CVD diamond heat spreader
• Can dissipate 1 kW of heat with active cooling ⇒ 1kW signal power!
Copper Mount
![Page 24: Computational Photonics: Modeling Ultrashort Pulse Propagation …math.arizona.edu/~nrw/Spring_School_2005/moloney.pdf · White Light Continuum The Physics • White light continuum](https://reader034.vdocument.in/reader034/viewer/2022042802/5f397991417cc81eee013cf9/html5/thumbnails/24.jpg)
MIT Photonics MicropolisSimulation Goal: Resolve fine nanoscale featuresin an otherwise featureless landscape
3D + time Maxwell solvers for time domain analysis
![Page 25: Computational Photonics: Modeling Ultrashort Pulse Propagation …math.arizona.edu/~nrw/Spring_School_2005/moloney.pdf · White Light Continuum The Physics • White light continuum](https://reader034.vdocument.in/reader034/viewer/2022042802/5f397991417cc81eee013cf9/html5/thumbnails/25.jpg)
Vector Maxwell Simulators
Maxwell’s Equations:
0,0
;
=•∇=•∇∂∂
=×∇∂∂
−=×∇
→→
→→
BDtDH
tBE
rr
Constitutive relations:→→→→
+== PEDHB 00 , εµr
FDTD Method Yee scheme
Simulation size grows as N3 making memory and CPU prohibitive!
![Page 26: Computational Photonics: Modeling Ultrashort Pulse Propagation …math.arizona.edu/~nrw/Spring_School_2005/moloney.pdf · White Light Continuum The Physics • White light continuum](https://reader034.vdocument.in/reader034/viewer/2022042802/5f397991417cc81eee013cf9/html5/thumbnails/26.jpg)
Layout of Computational Domain
PML Layer
Micro-ring add-drop filter
![Page 27: Computational Photonics: Modeling Ultrashort Pulse Propagation …math.arizona.edu/~nrw/Spring_School_2005/moloney.pdf · White Light Continuum The Physics • White light continuum](https://reader034.vdocument.in/reader034/viewer/2022042802/5f397991417cc81eee013cf9/html5/thumbnails/27.jpg)
Example of AMR FDTD discretization for TM-mode in two space dimensions
AMR refines the computational domain locally using nested rectangular grid patches. A standard FDTD update is applied to each patch.
At the coarse/fine grid interfacesthe solution is interpolated. Dashed lines denote boundaries of the ghost cells around the fine region. Arrows show a sample interpolation from coarse to fine values of the electric field.
, yx z z
y xz
HH E Et y t x
H HEt x y
∂∂ ∂ ∂= − =
∂ ∂ ∂ ∂∂ ∂∂
= −∂ ∂ ∂
( )
( )
1/ 2 1/ 2, 1/ 2, , 1/ 2, , 1/ 2, 1/ 2 , 1/ 2, 1/ 2
1/ 2 1/ 2, 1/ 2, , 1/ 2, , 1/ 2, 1/ 2 , 1/ 2, 1/ 2
n n n nx i j x i j z i j z i j
n n n ny i j y i j z i j z i j
tH H E EytH H E Ex
+ −+ + + + + −
+ −+ + + + − +
∆− = − −
∆∆
− = −∆
- A.R.Zakharian, M. Brio
![Page 28: Computational Photonics: Modeling Ultrashort Pulse Propagation …math.arizona.edu/~nrw/Spring_School_2005/moloney.pdf · White Light Continuum The Physics • White light continuum](https://reader034.vdocument.in/reader034/viewer/2022042802/5f397991417cc81eee013cf9/html5/thumbnails/28.jpg)
Nested Grids on a 3D PBG Structure
- confined defect mode of a 3D PBG lattice
- A. Zakharian, C. Dineen
![Page 29: Computational Photonics: Modeling Ultrashort Pulse Propagation …math.arizona.edu/~nrw/Spring_School_2005/moloney.pdf · White Light Continuum The Physics • White light continuum](https://reader034.vdocument.in/reader034/viewer/2022042802/5f397991417cc81eee013cf9/html5/thumbnails/29.jpg)
Resolving a QD in a PBG Lattice-Jens Foerstner - experiment Gibbs et al.
230 240 250 260 270 280 290 300
emis
sion
frequency [THz]
QD Wavefunction
Coupling QD to High Q Defect Mode
Energy level splitting
4-20 nm
AMR Mesh
![Page 30: Computational Photonics: Modeling Ultrashort Pulse Propagation …math.arizona.edu/~nrw/Spring_School_2005/moloney.pdf · White Light Continuum The Physics • White light continuum](https://reader034.vdocument.in/reader034/viewer/2022042802/5f397991417cc81eee013cf9/html5/thumbnails/30.jpg)
PBG+QD Confined/Radiating ModesBare PBG Mode from Top Coupled QD+PBG Energy Transfer
Radiating Pattern
![Page 31: Computational Photonics: Modeling Ultrashort Pulse Propagation …math.arizona.edu/~nrw/Spring_School_2005/moloney.pdf · White Light Continuum The Physics • White light continuum](https://reader034.vdocument.in/reader034/viewer/2022042802/5f397991417cc81eee013cf9/html5/thumbnails/31.jpg)
z Nlevels Runtime ratio
2nm 2 2.3 1nm 3 8.1
Nested 3-level AMR Resolves SP on 20nm Metal Film
![Page 32: Computational Photonics: Modeling Ultrashort Pulse Propagation …math.arizona.edu/~nrw/Spring_School_2005/moloney.pdf · White Light Continuum The Physics • White light continuum](https://reader034.vdocument.in/reader034/viewer/2022042802/5f397991417cc81eee013cf9/html5/thumbnails/32.jpg)
Scattering from a metal Nanosphere – 3 nested AMR Levels
![Page 33: Computational Photonics: Modeling Ultrashort Pulse Propagation …math.arizona.edu/~nrw/Spring_School_2005/moloney.pdf · White Light Continuum The Physics • White light continuum](https://reader034.vdocument.in/reader034/viewer/2022042802/5f397991417cc81eee013cf9/html5/thumbnails/33.jpg)
Nested AMR Mesh for a Plasmon Waveguide
- shows 3 levels of mesh refinement
- A.R. Zakharian, C. Dineen
![Page 34: Computational Photonics: Modeling Ultrashort Pulse Propagation …math.arizona.edu/~nrw/Spring_School_2005/moloney.pdf · White Light Continuum The Physics • White light continuum](https://reader034.vdocument.in/reader034/viewer/2022042802/5f397991417cc81eee013cf9/html5/thumbnails/34.jpg)
Summary
Why Computational Photonics?
• Provides a fast track to device design, dramatically reducing wasteful material re-growth, packaging etc
• Interactive simulation laboratory will make many future experiments redundant – computer laboratory on a bench
• Powerful tool to discover new linear and nonlinear optical phenomena on micro- and nano-scales
-Many basic research challenges remain in physics and algorithmdevelopment
![Page 35: Computational Photonics: Modeling Ultrashort Pulse Propagation …math.arizona.edu/~nrw/Spring_School_2005/moloney.pdf · White Light Continuum The Physics • White light continuum](https://reader034.vdocument.in/reader034/viewer/2022042802/5f397991417cc81eee013cf9/html5/thumbnails/35.jpg)
Lateral ASE and Lasing
4.1 kW/cm2 28.1 kW/cm216.4 kW/cm2
scattering from sample edge
Dark lines indicate carrier depletion from lateral lasing mode
Imaging Photoluminescence
• UA, AFRL and Coherent Inc. Collaboration
• Expt data courtesy of Coherent and Robert Bedford(AFRL)