Computation of rare transitions in the barotropic quasi-geostrophic equations

We investigate the theoretical and numerical computation of rare transitions in simple geophysical turbulent models. We consider the barotropic quasi-geostrophic and two-dimensional Navier–Stokes equations in regimes where bistability between two coexisting large-scale attractors exist. By means of large deviations and instanton theory with the use of an Onsager–Machlup path integral formalism for the transition probability, we show how one can directly compute the most probable transition path between two coexisting attractors analytically in an equilibrium (Langevin) framework and numerically otherWe adapt a class of numerical optimization algorithms known as minimum action methods to simple geophysical turbulent models. We show that by numerically minimizing an appropriate action functional in a large deviation limit, one can predict the most likely transition path for a rare transition between two states. By considering examples where theoretical predictions can be made, we show that the minimum action method successfully predicts the most likely transition path. Finally, we discuss the application and extension of such numerical optimization schemes to the computation of rare transitions observed in direct numerical simulations and experiments and to other, more complex, turbulent systems.

A comparison of variational and Markov chain Monte Carlo methods for inference in partially observed stochastic dynamic systems

In recent work we have developed a novel variational inference method for partially observed systems governed by stochastic differential equations. In this paper we provide a comparison of the Variational Gaussian Process Smoother with an exact solution computed using a Hybrid Monte Carlo approach to path sampling, applied to a stochastic double well potential model. It is demonstrated that the variational smoother provides us a very accurate estimate of mean path while conditional variance is slightly underestimated. We conclude with some remarks as to the advantages and disadvantages of the variational smoother. © 2008 Springer Science + Business Media LLC.

Gaussian process approximations of stochastic differential equations

Stochastic differential equations arise naturally in a range of contexts, from financial to environmental modeling. Current solution methods are limited in their representation of the posterior process in the presence of data. In this work, we present a novel Gaussian process approximation to the posterior measure over paths for a general class of stochastic differential equations in the presence of observations. The method is applied to two simple problems: the Ornstein-Uhlenbeck process, of which the exact solution is known and can be compared to, and the double-well system, for which standard approaches such as the ensemble Kalman smoother fail to provide a satisfactory result. Experiments show that our variational approximation is viable and that the results are very promising as the variational approximate solution outperforms standard Gaussian process regression for non-Gaussian Markov processes.

Dissipative nonlinear structures in fiber optics

Optical fiber materials exhibit a nonlinear response to strong electric fields, such as those of optical signals confined within the small fiber core. Fiber nonlinearity is an essential component in the design of the next generation of advanced optical communication systems, but its use is often avoided by engineers because of its intractability. The application of nonlinear technologies in fiber optics offers new opportunities for the design of photonic systems and devices. In this chapter, we make an overview of recent progress in mathematical theory and practical applications of temporal dissipative solitons and self-similar nonlinear structures in optical fiber systems. The design of all-optical high-speed signal processing devices, based on nonlinear dissipative structures, is discussed.

A variational method and approximations of a Cauchy problem for elliptic equations

A Cauchy problem for general elliptic second-order linear partial differential equations in which the Dirichlet data in H½(?1 ? ?3) is assumed available on a larger part of the boundary ? of the bounded domain O than the boundary portion ?1 on which the Neumann data is prescribed, is investigated using a conjugate gradient method. We obtain an approximation to the solution of the Cauchy problem by minimizing a certain discrete functional and interpolating using the finite diference or boundary element method. The minimization involves solving equations obtained by discretising mixed boundary value problems for the same operator and its adjoint. It is proved that the solution of the discretised optimization problem converges to the continuous one, as the mesh size tends to zero. Numerical results are presented and discussed.

Dissipative nonlinear structures in fiber optics

Optical fiber materials exhibit a nonlinear response to strong electric fields, such as those of optical signals confined within the small fiber core. Fiber nonlinearity is an essential component in the design of the next generation of advanced optical communication systems, but its use is often avoided by engineers because of its intractability. The application of nonlinear technologies in fiber optics offers new opportunities for the design of photonic systems and devices. In this chapter, we make an overview of recent progress in mathematical theory and practical applications of temporal dissipative solitons and self-similar nonlinear structures in optical fiber systems. The design of all-optical high-speed signal processing devices, based on nonlinear dissipative structures, is discussed.

Population lensing vs thermal lensing in Yb:YAG rods and disks

The induced lenses in the Yb:YAG rods and disks end-pumped by a Gaussian beam were analyzed both analytically and numerically. The thermally assisted mechanisms of the lens formation were considered to include: the conventional volume thermal index changes ("dn/dT"), the bulging of end faces, the photoelastic effect, and the bending (for a disk). The heat conduction equations (with an axial heat flux for a disk and a radial heat flux for a rod), and quasi-static thermoelastic equations (in the plane-stress approximation with free boundary conditions) were solved to find the thermal lens power. The population rate equation with saturation (by amplified spontaneous emission or an external wave) was examined to find the electronic lens power in the active elements.

Ion irradiation effects in nanocrystalline TiN coatings

Nitride materials and coatings have attracted extensive research interests for various applications in advanced nuclear reactors due to their unique combination of physical properties, including high temperature stability, excellent corrosion resistance, superior mechanical property and good thermal conductivity. In this paper, the ion irradiation effects in nanocrystalline TiN coatings as a function of grain size are reported. TiN thin films (thickness of 100 nm) with various grain sizes (8-100 nm) were prepared on Si substrates by a pulsed laser deposition technique. All the samples were irradiated with He ions to high fluences at room temperature. Transmission electron microscopy (TEM) and high resolution TEM on the ion-irradiated samples show that damage accumulation in the TiN films reduces as the grain size reduces. Electrical resistivity of the ion-irradiated films increases slightly compared with the as-deposited ones. These observations demonstrate a good radiation-tolerance property of nanocrystalline TiN films. © 2007 Elsevier B.V. All rights reserved.

Modeling of the spectrum in a random distributed feedback fiber laser within the power balance modes

The simplest model for a description of the random distributed feedback (RDFB) Raman fiber laser is a power balance model describing the evolution of the intensities of the waves over the fiber length. The model predicts well the power performances of the RDFB fiber laser including the generation threshold, the output power and pump and generation wave intensity distributions along the fiber. In the present work, we extend the power balance model and modify equations in such a way that they describe now frequency dependent spectral power density instead of integral over the spectrum intensities. We calculate the generation spectrum by using the depleted pump wave longitudinal distribution derived from the conventional power balance model. We found the spectral balance model to be sufficient to account for the spectral narrowing in the RDFB laser above the threshold of the generation. © 2014 SPIE.

Parabolic similaritons in optical fibres

Recent developments in nonlinear optics have brought to the fore of intensive research an interesting class of pulses with a parabolic intensity profile and a linear instantaneous frequency shift or chirp. Parabolic pulses propagate in optical fibres with normal group-velocity dispersion in a self-similar manner, holding certain relations (scaling) between pulse power, duration and chirp parameter, and can tolerate strong nonlinearity without distortion or wave breaking. These solutions, which have been dubbed similaritons, were demonstrated theoretically and experimentally in fibre amplifiers in 2000. Similaritons in fibre amplifiers are, along with solitons in passive fibres, the most well-known classes of nonlinear attractors for pulse propagation in optical fibre, so they take on major fundamental importance. The unique properties of parabolic similaritons have stimulated numerous applications in nonlinear optics, ranging from ultrashort high-power pulse generation to highly coherent continuum sources and to optical nonlinear processing of telecommunication signals. In this work, we review the physics underlying the generation of parabolic similaritons as well as recent results obtained in a wide range of experimental configurations.

Impact of the carrier relaxation paths on two-state operation in quantum dot lasers

We study InGaAs QD laser operating simultaneously at ground (GS) and excited (ES) states under 30ns pulsed-pumping and distinguish three regimes of operation depending on the pump current and the carrier relaxation pathways. An increased current leads to an increase in ES intensity and to a decrease in GS intensity (or saturation) for low pump range, as typical for the cascade-like pathway. Both the GS and ES intensities are steadily increased for high current ranges, which prove the dominance of the direct capture pathway. The relaxation oscillations are not pronounced for these ranges. For the mediate currents, the interplay between the both pathways leads to the damped large amplitude relaxation oscillations with significant deviation of the relaxation oscillation frequency from the initial value during the pulse.

Optical wave turbulence and wave condensation in a nonlinear optical experiment

We present theory, numerical simulations and experimental observations of a 1D optical wave system. We show that this system is of a dual cascade type, namely, the energy cascading directly to small scales, and the photons or wave action cascading to large scales. In the optical context the inverse cascade is particularly interesting because it means the condensation of photons. We show that the cascades are induced by a six-wave resonant interaction process described by weak turbulence theory. We show that by starting with weakly nonlinear randomized waves as an initial condition, there exists an inverse cascade of photons towards the lowest wavenumbers. During the cascade nonlinearity becomes strong at low wavenumbers and, due to the focusing nature of the nonlinearity, it leads to modulational instability resulting in the formation of solitons. Further interaction of the solitons among themselves and with incoherent waves leads to the final condensate state dominated by a single strong soliton. In addition, we show the existence of the direct energy cascade numerically and that it agrees with the wave turbulence prediction.

Bifurcating trajectory of non-diffractive electromagnetic airy pulse

The explicit expression for spatial-temporal Airy pulse is derived from the Maxwell's equations in paraxial approximation. The trajectory of the pulse in the time-space coordinates is analysed. The existence of a bifurcation point that separates regions with qualitatively different features of the pulse propagation is demonstrated. At this point the velocity of the pulse becomes infinite and the orientation of it changes to the opposite.