Nonlinear Dynamics of Nd:YAG lasers: Hopf bifurcation, Multistability and Chaotic Synchronization

In this thesis we have presented some aspects of the nonlinear dynamics of Nd:YAG lasers including synchronization, Hopf bifurcation, chaos control and delay induced multistability.We have chosen diode pumped Nd:YAG laser with intracavity KTP crystal operating with two mode and three mode output as our model system.Different types of orientation for the laser cavity modes were considered to carry out the studies. For laser operating with two mode output we have chosen the modes as having parallel polarization and perpendicular polarization. For laser having three mode output, we have chosen them as two modes polarized parallel to each other while the third mode polarized orthogonal to them.

Numerical study of reverse period doubling route from chaos to stability in a two-mode intracavity doubled Nd-YAG laser

We have numerically studied the behavior of a two-mode Nd-YAG laser with an intracavity KTP crystal. It is found that when the parameter, which is a measure of the relative orientations of the KTP crystal with respect to the Nd-YAG crystal, is varied continuously, the output intensity fluctuations change from chaotic to stable behavior through a sequence of reverse period doubling bifurcations. The graph of the intensity in the X-polarized mode against that in the Y-polarized mode shows a complex pattern in the chaotic regime. The Lyapunov exponent is calculated for the chaotic and periodic regions.

Studies on some nonlinear schrodinger type equations describing pulse propagation through optical fibers

The discovery of the soliton is considered to be one of the most significant events of the twentieth century. The term soliton refers to special kinds of waves that can propagate undistorted over long distances and remain unaffected even after collision with each other. Solitons have been studied extensively in many fields of physics. In the context of optical fibers, solitons are not only of fundamental interest but also have potential applications in the field of optical fiber communications. This thesis is devoted to the theoretical study of soliton pulse propagation through single mode optical fibers.

Femtosecond imaging-mode laser-induced breakdown spectroscopy

Many nonlinear optical microscopy techniques based on the high-intensity nonlinear phenomena were developed recent years. A new technique based on the minimal-invasive in-situ analysis of the specific bound elements in biological samples is described in the present work. The imaging-mode Laser-Induced Breakdown Spectroscopy (LIBS) is proposed as a combination of LIBS, femtosecond laser material processing and microscopy. The Calcium distribution in the peripheral cell wall of the sunflower seedling (Helianthus Annuus L.) stem is studied as a first application of the imaging-mode LIBS. At first, several nonlinear optical microscopy techniques are overviewed. The spatial resolution of the imaging-mode LIBS microscope is discussed basing on the Point-Spread Function (PSF) concept. The primary processes of the Laser-Induced Breakdown (LIB) are overviewed. We consider ionization, breakdown, plasma formation and ablation processes. Water with defined Calcium salt concentration is used as a model of the biological object in the preliminary experiments. The transient LIB spectra are measured and analysed for both nanosecond and femtosecond laser excitation. The experiment on the local Calcium concentration measurements in the peripheral cell wall of the sunflower seedling stem employing nanosecond LIBS shows, that nanosecond laser is not a suitable excitation source for the biological applications. In case of the nanosecond laser the ablation craters have random shape and depth over 20 µm. The analysis of the femtosecond laser ablation craters shows the reproducible circle form. At 3.5 µJ laser pulse energy the diameter of the crater is 4 µm and depth 140 nm for single laser pulse, which results in 1 femtoliter analytical volume. The experimental result of the 2 dimensional and surface sectioning of the bound Calcium concentrations is presented in the work.

Variational calculation of static and dynamic vibrational nonlinear optical properties

The vibrational configuration interaction method used to obtain static vibrational (hyper)polarizabilities is extended to dynamic nonlinear optical properties in the infinite optical frequency approximation. Illustrative calculations are carried out on H2 O and N H3. The former molecule is weakly anharmonic while the latter contains a strongly anharmonic umbrella mode. The effect on vibrational (hyper)polarizabilities due to various truncations of the potential energy and property surfaces involved in the calculation are examined

The counter-propagating Rossby-wave perspective on baroclinic instability. Part IV: Nonlinear life cycles

Pairs of counter-propagating Rossby waves (CRWs) can be used to describe baroclinic instability in linearized primitive-equation dynamics, employing simple propagation and interaction mechanisms at only two locations in the meridional plane—the CRW ‘home-bases’. Here, it is shown how some CRW properties are remarkably robust as a growing baroclinic wave develops nonlinearly. For example, the phase difference between upper-level and lower-level waves in potential-vorticity contours, defined initially at the home-bases of the CRWs, remains almost constant throughout baroclinic wave life cycles, despite the occurrence of frontogenesis and Rossby-wave breaking. As the lower wave saturates nonlinearly the whole baroclinic wave changes phase speed from that of the normal mode to that of the self-induced phase speed of the upper CRW. On zonal jets without surface meridional shear, this must always act to slow the baroclinic wave. The direction of wave breaking when a basic state has surface meridional shear can be anticipated because the displacement structures of CRWs tend to be coherent along surfaces of constant basic-state angular velocity, U. This results in up-gradient horizontal momentum fluxes for baroclinically growing disturbances. The momentum flux acts to shift the jet meridionally in the direction of the increasing surface U, so that the upper CRW breaks in the same direction as occurred at low levels

Stochastic resonance in a nonlinear model of a rotating, stratified shear flow, with a simple stochastic inertia-gravity wave parameterization

We report on a numerical study of the impact of short, fast inertia-gravity waves on the large-scale, slowly-evolving flow with which they co-exist. A nonlinear quasi-geostrophic numerical model of a stratified shear flow is used to simulate, at reasonably high resolution, the evolution of a large-scale mode which grows due to baroclinic instability and equilibrates at finite amplitude. Ageostrophic inertia-gravity modes are filtered out of the model by construction, but their effects on the balanced flow are incorporated using a simple stochastic parameterization of the potential vorticity anomalies which they induce. The model simulates a rotating, two-layer annulus laboratory experiment, in which we recently observed systematic inertia-gravity wave generation by an evolving, large-scale flow. We find that the impact of the small-amplitude stochastic contribution to the potential vorticity tendency, on the model balanced flow, is generally small, as expected. In certain circumstances, however, the parameterized fast waves can exert a dominant influence. In a flow which is baroclinically-unstable to a range of zonal wavenumbers, and in which there is a close match between the growth rates of the multiple modes, the stochastic waves can strongly affect wavenumber selection. This is illustrated by a flow in which the parameterized fast modes dramatically re-partition the probability-density function for equilibrated large-scale zonal wavenumber. In a second case study, the stochastic perturbations are shown to force spontaneous wavenumber transitions in the large-scale flow, which do not occur in their absence. These phenomena are due to a stochastic resonance effect. They add to the evidence that deterministic parameterizations in general circulation models, of subgrid-scale processes such as gravity wave drag, cannot always adequately capture the full details of the nonlinear interaction.

Empirical mode decomposition: a novel technique for the study of tremor time series

Tremor is a clinical feature characterized by oscillations of a part of the body. The detection and study of tremor is an important step in investigations seeking to explain underlying control strategies of the central nervous system under natural (or physiological) and pathological conditions. It is well established that tremorous activity is composed of deterministic and stochastic components. For this reason, the use of digital signal processing techniques (DSP) which take into account the nonlinearity and nonstationarity of such signals may bring new information into the signal analysis which is often obscured by traditional linear techniques (e.g. Fourier analysis). In this context, this paper introduces the application of the empirical mode decomposition (EMD) and Hilbert spectrum (HS), which are relatively new DSP techniques for the analysis of nonlinear and nonstationary time-series, for the study of tremor. Our results, obtained from the analysis of experimental signals collected from 31 patients with different neurological conditions, showed that the EMD could automatically decompose acquired signals into basic components, called intrinsic mode functions (IMFs), representing tremorous and voluntary activity. The identification of a physical meaning for IMFs in the context of tremor analysis suggests an alternative and new way of detecting tremorous activity. These results may be relevant for those applications requiring automatic detection of tremor. Furthermore, the energy of IMFs was visualized as a function of time and frequency by means of the HS. This analysis showed that the variation of energy of tremorous and voluntary activity could be distinguished and characterized on the HS. Such results may be relevant for those applications aiming to identify neurological disorders. In general, both the HS and EMD demonstrated to be very useful to perform objective analysis of any kind of tremor and can therefore be potentially used to perform functional assessment.

State estimation using the particle filter with mode tracking

A particle filter is a data assimilation scheme that employs a fully nonlinear, non-Gaussian analysis step. Unfortunately as the size of the state grows the number of ensemble members required for the particle filter to converge to the true solution increases exponentially. To overcome this Vaswani [Vaswani N. 2008. IEEE Trans Signal Process 56:4583–97] proposed a new method known as mode tracking to improve the efficiency of the particle filter. When mode tracking, the state is split into two subspaces. One subspace is forecast using the particle filter, the other is treated so that its values are set equal to the mode of the marginal pdf. There are many ways to split the state. One hypothesis is that the best results should be obtained from the particle filter with mode tracking when we mode track the maximum number of unimodal dimensions. The aim of this paper is to test this hypothesis using the three dimensional stochastic Lorenz equations with direct observations. It is found that mode tracking the maximum number of unimodal dimensions does not always provide the best result. The best choice of states to mode track depends on the number of particles used and the accuracy and frequency of the observations.

Nonlinear stability of Euler flows in two-dimensional periodic domains

The non-quadratic conservation laws of the two-dimensional Euler equations are used to show that the gravest modes in a doubly-periodic domain with aspect ratio L = 1 are stable up to translations (or structurally stable) for finite-amplitude disturbances. This extends a previous result based on conservation of energy and enstrophy alone. When L 1, a saturation bound is established for the mode with wavenumber |k| = L −1 (the next-gravest mode), which is linearly unstable. The method is applied to prove nonlinear structural stability of planetary wave two on a rotating sphere.

A new understanding of El Niño's impact over East Asia: dominance of the ENSO combination mode

Previous studies have shown that the Indo-Pacific atmospheric response to ENSO comprises two dominant modes of variability: a meridionally quasi-symmetric response (independent from the annual cycle) and an anti-symmetric response (arising from the nonlinear atmospheric interaction between ENSO variability and the annual cycle), referred to as the combination mode (C-Mode). This study demonstrates that the direct El Niño signal over the tropics is confined to the equatorial region and has no significant impact on the atmospheric response over East Asia. The El Niño-associated equatorial anomalies can be expanded towards off-equatorial regions by the C-Mode through ENSO’s interaction with the annual cycle. The C-Mode is the prime driver for the development of an anomalous low-level anticyclone over the western North Pacific (WNP) during the El Niño decay phase, which usually transports more moisture to East Asia and thereby causes more precipitation over southern China. We use an Atmospheric General Circulation Model that well reproduces the WNP anticyclonic anomalies when both El Niño sea surface temperature (SST) anomalies as well as the SST annual cycle are prescribed as boundary conditions. However, no significant WNP anticyclonic circulation anomaly appears during the El Niño decay phase when excluding the SST annual cycle. Our analyses of observational data and model experiments suggest that the annual cycle plays a key role in the East Asian climate anomalies associated with El Niño through their nonlinear atmospheric interaction. Hence, a realistic simulation of the annual cycle is crucial in order to correctly capture the ENSO-associated climate anomalies over East Asia.

Blowout bifurcation and spatial mode excitation in the bubbling transition to turbulence

The transition to turbulence (spatio-temporal chaos) in a wide class of spatially extended dynamical system is due to the loss of transversal stability of a chaotic attractor lying on a homogeneous manifold (in the Fourier phase space of the system) causing spatial mode excitation Since the latter manifests as intermittent spikes this has been called a bubbling transition We present numerical evidences that this transition occurs due to the so called blowout bifurcation whereby the attractor as a whole loses transversal stability and becomes a chaotic saddle We used a nonlinear three-wave interacting model with spatial diffusion as an example of this transition (C) 2010 Elsevier B V All rights reserved

Bubbling transition to spatial mode excitation in an extended dynamical system

We investigated the transition to spatio-temporal chaos in spatially extended nonlinear dynamical systems possessing an invariant subspace with a low-dimensional attractor. When the latter is chaotic and the subspace is transversely stable we have a spatially homogeneous state only. The onset of spatio-temporal chaos, i.e. the excitation of spatially inhomogeneous modes, occur through the loss of transversal stability of some unstable periodic orbit embedded in the chaotic attractor lying in the invariant subspace. This is a bubbling transition, since there is a switching between spatially homogeneous and nonhomogeneous states with statistical properties of on-off intermittency. Hence the onset of spatio-temporal chaos depends critically both on the existence of a chaotic attractor in the invariant subspace and its being transversely stable or unstable. (C) 2008 Elsevier B.V. All rights reserved.

Nonlinear refractive index of RECOB (RE = Gd and La) crystals

The Z-scan technique is employed to obtain the nonlinear refractive index (n (2)) of the Ca(4)REO(BO(3))(3) (RECOB, where RE = Gd and La) single crystals using 30 fs laser pulses centered at 780 nm for the two orthogonal orientations determined by the optical axes (X and Z) relative to the direction of propagation of the laser beam (k//Y// crystallographic b-axis). The large values of n (2) indicate that both GdCOB and LaCOB are potential hosts for Yb:RECOB lasers operating in the Kerr-lens mode locking (KLM) regime.

Time-resolved study electronic and thermal contributions to the nonlinear refractive index of Nd3+: SBN laser crystals

The propagation of an optical beam through dielectric media induces changes in the refractive index, An, which causes self-focusing or self-defocusing. In the particular case of ion-doped solids, there are thermal and non-thermal lens effects, where the latter is due to the polarizability difference, Delta alpha, between the excited and ground states, the so-called population lens (PL) effect. PL is a pure electronic contribution to the nonlinearity, while the thermal lens (TL) effect is caused by the conversion of part of the absorbed energy into heat. In time-resolved measurements such as Z-scan and TL transient experiments, it is not easy to separate these two contributions to nonlinear refractive index because they usually have similar response times. In this work, we performed time-resolved measurements using both Z-scan and mode mismatched TL in order to discriminate thermal and electronic contributions to the laser-induced refractive index change of the Nd3+-doped Strontium Barium Niobate (SrxBa1-xNb2O6) laser crystal. Combining numerical simulations with experimental results we could successfully distinguish between the two contributions to An. (C) 2007 Elsevier B.V. All rights reserved.