Nonlinear optics and wavelength translation via cavity-optomechanics

The field of cavity-optomechanics explores the interaction of light with sound in an ever increasing array of devices. This interaction allows the mechanical system to be both sensed and controlled by the optical system, opening up a wide variety of experiments including the cooling of the mechanical resonator to its quantum mechanical ground state and the squeezing of the optical field upon interaction with the mechanical resonator, to name two.

In this work we explore two very different systems with different types of optomechanical coupling. The first system consists of two microdisk optical resonators stacked on top of each other and separated by a very small slot. The interaction of the disks causes their optical resonance frequencies to be extremely sensitive to the gap between the disks. By careful control of the gap between the disks, the optomechanical coupling can be made to be quadratic to first order which is uncommon in optomechanical systems. With this quadratic coupling the light field is now sensitive to the energy of the mechanical resonator and can directly control the potential energy trapping the mechanical motion. This ability to directly control the spring constant without modifying the energy of the mechanical system, unlike in linear optomechanical coupling, is explored.

Next, the bulk of this thesis deals with a high mechanical frequency optomechanical crystal which is used to coherently convert photons between different frequencies. This is accomplished via the engineered linear optomechanical coupling in these devices. Both classical and quantum systems utilize the interaction of light and matter across a wide range of energies. These systems are often not naturally compatible with one another and require a means of converting photons of dissimilar wavelengths to combine and exploit their different strengths. Here we theoretically propose and experimentally demonstrate coherent wavelength conversion of optical photons using photon-phonon translation in a cavity-optomechanical system. For an engineered silicon optomechanical crystal nanocavity supporting a 4 GHz localized phonon mode, optical signals in a 1.5 MHz bandwidth are coherently converted over a 11.2 THz frequency span between one cavity mode at wavelength 1460 nm and a second cavity mode at 1545 nm with a 93% internal (2% external) peak efficiency. The thermal and quantum limiting noise involved in the conversion process is also analyzed and, in terms of an equivalent photon number signal level, are found to correspond to an internal noise level of only 6 and 4 times 10x^-3 quanta, respectively.

We begin by developing the requisite theoretical background to describe the system. A significant amount of time is then spent describing the fabrication of these silicon nanobeams, with an emphasis on understanding the specifics and motivation. The experimental demonstration of wavelength conversion is then described and analyzed. It is determined that the method of getting photons into the cavity and collected from the cavity is a fundamental limiting factor in the overall efficiency. Finally, a new coupling scheme is designed, fabricated, and tested that provides a means of coupling greater than 90% of photons into and out of the cavity, addressing one of the largest obstacles with the initial wavelength conversion experiment.

Acceleration Sensing, Feedback Cooling, and Nonlinear Dynamics with Nanoscale Cavity-Optomechanical Devices

Light has long been used for the precise measurement of moving bodies, but the burgeoning field of optomechanics is concerned with the interaction of light and matter in a regime where the typically weak radiation pressure force of light is able to push back on the moving object. This field began with the realization in the late 1960's that the momentum imparted by a recoiling photon on a mirror would place fundamental limits on the smallest measurable displacement of that mirror. This coupling between the frequency of light and the motion of a mechanical object does much more than simply add noise, however. It has been used to cool objects to their quantum ground state, demonstrate electromagnetically-induced-transparency, and modify the damping and spring constant of the resonator. Amazingly, these radiation pressure effects have now been demonstrated in systems ranging 18 orders of magnitude in mass (kg to fg).

In this work we will focus on three diverse experiments in three different optomechanical devices which span the fields of inertial sensors, closed-loop feedback, and nonlinear dynamics. The mechanical elements presented cover 6 orders of magnitude in mass (ng to fg), but they all employ nano-scale photonic crystals to trap light and resonantly enhance the light-matter interaction. In the first experiment we take advantage of the sub-femtometer displacement resolution of our photonic crystals to demonstrate a sensitive chip-scale optical accelerometer with a kHz-frequency mechanical resonator. This sensor has a noise density of approximately 10 micro-g/rt-Hz over a useable bandwidth of approximately 20 kHz and we demonstrate at least 50 dB of linear dynamic sensor range. We also discuss methods to further improve performance of this device by a factor of 10.

In the second experiment, we used a closed-loop measurement and feedback system to damp and cool a room-temperature MHz-frequency mechanical oscillator from a phonon occupation of 6.5 million down to just 66. At the time of the experiment, this represented a world-record result for the laser cooling of a macroscopic mechanical element without the aid of cryogenic pre-cooling. Furthermore, this closed-loop damping yields a high-resolution force sensor with a practical bandwidth of 200 kHZ and the method has applications to other optomechanical sensors.

The final experiment contains results from a GHz-frequency mechanical resonator in a regime where the nonlinearity of the radiation-pressure interaction dominates the system dynamics. In this device we show self-oscillations of the mechanical element that are driven by multi-photon-phonon scattering. Control of the system allows us to initialize the mechanical oscillator into a stable high-amplitude attractor which would otherwise be inaccessible. To provide context, we begin this work by first presenting an intuitive overview of optomechanical systems and then providing an extended discussion of the principles underlying the design and fabrication of our optomechanical devices.

Part I - The radiation of elastic waves from a spherical cavity in a half space. Part II - Precision determination of focal depths and epicenters of earthquakes.

Part I: The dynamic response of an elastic half space to an explosion in a buried spherical cavity is investigated by two methods. The first is implicit, and the final expressions for the displacements at the free surface are given as a series of spherical wave functions whose coefficients are solutions of an infinite set of linear equations. The second method is based on Schwarz's technique to solve boundary value problems, and leads to an iterative solution, starting with the known expression for the point source in a half space as first term. The iterative series is transformed into a system of two integral equations, and into an equivalent set of linear equations. In this way, a dual interpretation of the physical phenomena is achieved. The systems are treated numerically and the Rayleigh wave part of the displacements is given in the frequency domain. Several comparisons with simpler cases are analyzed to show the effect of the cavity radius-depth ratio on the spectra of the displacements.

Part II: A high speed, large capacity, hypocenter location program has been written for an IBM 7094 computer. Important modifications to the standard method of least squares have been incorporated in it. Among them are a new way to obtain the depth of shocks from the normal equations, and the computation of variable travel times for the local shocks in order to account automatically for crustal variations. The multiregional travel times, largely based upon the investigations of the United States Geological Survey, are confronted with actual traverses to test their validity.

It is shown that several crustal phases provide control enough to obtain good solutions in depth for nuclear explosions, though not all the recording stations are in the region where crustal corrections are considered. The use of the European travel times, to locate the French nuclear explosion of May 1962 in the Sahara, proved to be more adequate than previous work.

A simpler program, with manual crustal corrections, is used to process the Kern County series of aftershocks, and a clearer picture of tectonic mechanism of the White Wolf fault is obtained.

Shocks in the California region are processed automatically and statistical frequency-depth and energy depth curves are discussed in relation to the tectonics of the area.

Minimum drag profiles in infinite cavity flows

The problem considered is that of minimizing the drag of a symmetric plate in infinite cavity flow under the constraints of fixed arclength and fixed chord. The flow is assumed to be steady, irrotational, and incompressible. The effects of gravity and viscosity are ignored.

Using complex variables, expressions for the drag, arclength, and chord, are derived in terms of two hodograph variables, Γ (the logarithm of the speed) and β (the flow angle), and two real parameters, a magnification factor and a parameter which determines how much of the plate is a free-streamline.

Two methods are employed for optimization:

(1) The parameter method. Γ and β are expanded in finite orthogonal series of N terms. Optimization is performed with respect to the N coefficients in these series and the magnification and free-streamline parameters. This method is carried out for the case N = 1 and minimum drag profiles and drag coefficients are found for all values of the ratio of arclength to chord.

(2) The variational method. A variational calculus method for minimizing integral functionals of a function and its finite Hilbert transform is introduced, This method is applied to functionals of quadratic form and a necessary condition for the existence of a minimum solution is derived. The variational method is applied to the minimum drag problem and a nonlinear integral equation is derived but not solved.

A singular perturbation method of calculating the behavior of supercavitating hydrofoils with rounded noses

A simple, direct and accurate method to predict the pressure distribution on supercavitating hydrofoils with rounded noses is presented. The thickness of body and cavity is assumed to be small. The method adopted in the present work is that of singular perturbation theory. Far from the leading edge linearized free streamline theory is applied. Near the leading edge, however, where singularities of the linearized theory occur, a non-linear local solution is employed. The two unknown parameters which characterize this local solution are determined by a matching procedure. A uniformly valid solution is then constructed with the aid of the singular perturbation approach.

The present work is divided into two parts. In Part I isolated supercavitating hydrofoils of arbitrary profile shape with parabolic noses are investigated by the present method and its results are compared with the new computational results made with Wu and Wang's exact "functional iterative" method. The agreement is very good. In Part II this method is applied to a linear cascade of such hydrofoils with elliptic noses. A number of cases are worked out over a range of cascade parameters from which a good idea of the behavior of this type of important flow configuration is obtained.

Some of the computational aspects of Wu and Wang's functional iterative method heretofore not successfully applied to this type of problem are described in an appendix.

Dynamic analysis of V-folded cavity for TEM00 operation of end-pumped solid-state laser

Based on graphic analysis design method of optical resonator, a simple design expression of V-folded cavity of end-pumped solid-state lasers with TEM00 operation is described, which satisfies two criterias of the resonator design. We give numerical simulation of spot size as a function of thermal focal length using this design approach whose advantages are validated experimentally.

A novel method to produce small droplets from large nozzles.

This work presents a new method to generate droplets with diameters significantly smaller than the nozzle from which they emerge. The electrical waveform used to produce the jetting consists of a single square negative pulse. The negative edge of the pressure wave pulls the meniscus in, overturning the surface in such a way that a cavity is created. This cavity is then forced to collapse under the action of the positive edge of the pressure wave. This violent collapse produces a thin jet that eventually breaks up and produces droplets. Four droplet generator prototypes that demonstrate the capabilities of this novel mechanism are described. It is also shown that the proposed mechanism extends the existing limits of the commonly accepted inkjet operating regime.

Finite element analysis of compressor disc heat transfer using the transient conduction method with sparse boundary collocation points

Surface temperature measurements from two discs of a gas turbine compressor rig are used as boundary conditions for the transient conduction solution (inverse heat transfer analysis). The disc geometry is complex, and so the finite element method is used. There are often large radial temperature gradients on the discs, and the equations are therefore solved taking into account the dependence of thermal conductivity on temperature. The solution technique also makes use of a multigrid algorithm to reduce the solution time. This is particularly important since a large amount of data must be analyzed to obtain correlations of the heat transfer. The finite element grid is also used for a network analysis to calculate the radiant heat transfer in the cavity formed between the two compressor discs. The work discussed here proved particularly challenging as the disc temperatures were only measured at four different radial locations. Four methods of surface temperature interpolation are examined, together with their effect on the local heat fluxes. It is found that the choice of interpolation method depends on the available number of data points. Bessel interpolation gives the best results for four data points, whereas cubic splines are preferred when there are considerably more data points. The results from the analysis of the compressor rig data show that the heat transfer near the disc inner radius appears to be influenced by the central throughflow. However, for larger radii, the heat transfer from the discs and peripheral shroud is found to be consistent with that of a buoyancy-induced flow.

Design of novel three port optical gates scheme for the integration of large optical cavity electroabsorption modulators and evanescently-coupled photodiodes

This paper presents a novel scheme to monolithically integrate an evanescently-coupled uni-travelling carrier photodiode with a planar short multimode waveguide structure and a large optical cavity electroabsorption modulator based on a multimode waveguide structure. By simulation, both electroabsorption modulator and photodiode show excellent optical performances. The device can be fabricated with conventional photolithography, reactive ion etching, and chemical wet etching.

Room-temperature photonic crystal laser in H3 cavity on InGaAsP/InP slab

We fabricate and investigate two-dimensional photonic crystal H3 microcavities in an InGaAsP slab. The lasing action at room temperature is observed. The lasering threshold is 7mW under the pulsed pump of 0.75% duty cycle. The Q factor and the lasing mode characteristics are simulated by three-dimensional finite difference time domain method. The simulation result matches well with the experiment.

Influence of fused interface on the optical and thermal characteristics of vertical cavity lasers

From the effective absorption coefficient of bonded interface and the relationship of interface to reflectivity at cavity mode for double bonded vertical cavity laser, it can be seen that bonded interfaces should be positioned at the null of standing wave distribution, and the thickness of interface should be less than 20 nm. Using the finite elements method, the temperature contour map of laser can be calculated. Results showed that the influence of thin interface to thermal characteristics of VCSELS is slight, while thick interface will lead to temperature increase of active region. SEM images demonstrate that hydrophobic bonding is suitable for the fabrication of the device, while hydrophilic bonding interface is unfavorable to optical and thermal properties of devices with interface thickness larger than 40 nm.

Narrow line-width resonant cavity enhanced photodetectors operating at 1.55 mu m

A method for fabrication of long-wavelength narrow line-width InGaAs resonant cavity enhanced (RCE) photodetectors in a silicon substrate operating at the wavelength range of 1.3-1.6 mu m has been developed. A full width at half maximum (FWHM) of 0.7 nm and a peak responsivity of 0. 16 A/W at the resonance wavelength of 1.55 mu m have been accomplished by using a thick InP layer as part of the resonant cavity. The effects of roughness and tilt of the InP layer surface, and its free carrier absorption, as well as the thickness deviation of the mirror pair on the resonance wavelength shift and the peak quantum efficiency of the RCE photodetectors are analyzed in detail, and approaches for minimizing them toward superior performance are suggested. (C) 2007 Elsevier B.V. All rights reserved.

Slow wave effect of 2-D photonic crystal coupled cavity array

We designed a two-dimensional coupled photonic crystal resonator array with hexagonal lattice. The calculation by plane-wave-expansion method shows that the dispersion curve of coupled cavity modes in the bandgap are much flattened in all directions in the reciprocal space. We simulated the transmission spectra of transverse electric (TE) wave along the Gamma K direction. Compared with the PC single cavity structure, the transmission ratio of the coupled cavity array increases more than three orders of magnitude, while the group velocity decreases to below 1/10, reaching 0.007c. The slow wave effect has potential application in the field of miniaturized tunable optical delay components and low-threshold photonic crystal lasers.

Reduced divergence angle of photonic crystal vertical-cavity surface-emitting laser

The reduced divergence angle of the photonic crystal vertical-cavity surface-emitting laser (PC-VCSEL) was investigated in both theory and experiment. The photonic crystal waveguide possessed the weakly guiding waveguide characteristic, which accounted for the reduction of the divergence angle. The three-dimensional finite-difference time-domain method was used to simulate the designed PC-VCSEL, and a calculated divergence angle of 5.2 degrees was obtained. The measured divergence angles of our fabricated PC-VCSEL were between 5.1 degrees and 5.5 degrees over the entire drive current range, consistent with the numerical results. This is the lowest divergence angle of the fabricated PC-VCSEL ever reported.

Single mode vertical-cavity surface emitting laser with surface plasmon nanostructure

The vertical-cavity surface-emitting laser(VCSEL) has proved to be a low cost light source with attractive properties such as surface emission, circular and low divergence output beam, and simple integration in two-dimensional array. Many new applications such as in spectroscopy, optical storage, short distance fiber optic interconnects, and in longer distance communication, are continuously arising. Many of these applications require stable and single-mode high output power. Several methods that affect the transverse guiding and/or introduce mode selective loss or gain have been developed. In this study, a method for improving the single mode output power by using metal surface plasmons nanostructure is proposed. Theoretical calculation shows that the outpout power is improved about 50% compared to the result of standard VCSELs.