Positive-Negative Birefringence in Multiferroic Layered Metasurfaces

We uncover and identify the regime for a magnetically and ferroelectrically controllable negative refraction of light traversing multiferroic, oxide-based metastructure consisting of alternating nanoscopic ferroelectric (SrTiO$_2$) and ferromagnetic (Y$_3$Fe$_2$(FeO$_4$)$_3$, YIG) layers. We perform analytical and numerical simulations based on discretized, coupled equations for the self-consistent Maxwell/ferroelectric/ferromagnetic dynamics and obtain a biquadratic relation for the refractive index. Various scenarios of ordinary and negative refraction in different frequency ranges are analyzed and quantified by simple analytical formula that are confirmed by full-fledge numerical simulations. Electromagnetic-waves injected at the edges of the sample are propagated exactly numerically. We discovered that for particular GHz frequencies, waves with different polarizations are characterized by different signs of the refractive index giving rise to novel types of phenomena such as a positive-negative birefringence effect, and magnetically controlled light trapping and accelerations

Electromagnetically controlled multiferroic thermal diode

We propose an electromagnetically tunable thermal diode based on a two phase multiferroics composite. Analytical and full numerical calculations for prototypical heterojunction composed of Iron on Barium titanate in the tetragonal phase demonstrate a strong heat rectification effect that can be controlled externally by a moderate electric field. This finding is of an importance for thermally based information processing and sensing and can also be integrated in (spin)electronic circuits for heat management and recycling.Comment: Accepted in Phys. Rev.

Nonlinear Band Gap Transmission in Optical Waveguide Arrays

The effect of nonlinear transmission in coupled optical waveguide arrays is theoretically investigated via numerical simulations on the corresponding model equations. The realistic experimental setup is suggested injecting the beam in a single boundary waveguide, linear refractive index of which ($n_0$) is larger than one ($n$) of other identical waveguides in the array. Particularly, the effect holds if $\omega(n_0-n)/c>2Q$, where $Q$ is a linear coupling constant between array waveguides, $\omega$ is a carrier wave frequency and $c$ is a light velocity. Making numerical experiments in case of discrete nonlinear Schr\"odinger equation it is shown that the energy transfers from the boundary waveguide to the waveguide array above certain threshold intensity of the injected beam. This effect is explained by means of the creation and propagation of gap solitons in full analogy with the similar phenomenon of nonlinear supratransmission [F. Geniet, J. Leon, PRL, {\bf 89}, 134102, (2002)] in case of discrete sine-Gordon lattice.Comment: 4 pages, 6 figures. Phys. Rev. Lett. (in press

Bistability in sine-Gordon: the ideal switch

The sine-Gordon equation, used as the representative nonlinear wave equation, presents a bistable behavior resulting from nonlinearity and generating hysteresis properties. We show that the process can be understood in a comprehensive analytical formulation and that it is a generic property of nonlinear systems possessing a natural band gap. The approach allows to discover that sine-Gordon can work as an it ideal switch by reaching a transmissive regime with vanishing driving amplitude.Comment: Phys. Rev. E, (to be published, May 2005

Driving light pulses with light in two-level media

A two-level medium, described by the Maxwell-Bloch (MB) system, is engraved by establishing a standing cavity wave with a linearly polarized electromagnetic field that drives the medium on both ends. A light pulse, polarized along the other direction, then scatters the medium and couples to the cavity standing wave by means of the population inversion density variations. We demonstrate that control of the applied amplitudes of the grating field allows to stop the light pulse and to make it move backward (eventually to drive it freely). A simplified limit model of the MB system with variable boundary driving is obtained as a discrete nonlinear Schroedinger equation with tunable external potential. It reproduces qualitatively the dynamics of the driven light pulse

Delocalization and spreading in a nonlinear Stark ladder

We study the evolution of a wave packet in a nonlinear Schr\"odinger lattice equation subject to a dc bias. In the absence of nonlinearity all normal modes are spatially localized giving rise to a Stark ladder with an equidistant eigenvalue spectrum and Bloch oscillations. Nonlinearity induces frequency shifts and mode-mode interactions and destroys localization. With increasing strength of nonlinearity we observe: (I) localization as a transient, with subsequent subdiffusion (weak mode-mode interactions); (II) immediate subdiffusion (strong mode-mode interactions); (III) single site trapping as a transient, with subsequent explosive spreading, followed by subdiffusion. For single mode excitations and weak nonlinearities stability intervals are predicted and observed upon variation of the dc bias strength, which affect the short and long time dynamics.Comment: 4 pages, 5 figure

Bistable light detectors with nonlinear waveguide arrays

Bistability induced by nonlinear Kerr effect in arrays of coupled waveguides is studied and shown to be a means to conceive light detectors that switch under excitation by a weak signal. The detector is obtained by coupling two single 1D waveguide to an array of coupled waveguides with adjusted indices and coupling. The process is understood by analytical description in the conservative and continuous case and illustrated by numerical simulations of the model with attenuation.Comment: Phys. Rev. Lett., v.94, (2005, to be published

Supersonic Discrete Kink-Solitons and Sinusoidal Patterns with "Magic" wavenumber in Anharmonic Lattices

The sharp pulse method is applied to Fermi-Pasta-Ulam (FPU) and Lennard-Jones (LJ) anharmonic lattices. Numerical simulations reveal the presence of high energy strongly localized ``discrete'' kink-solitons (DK), which move with supersonic velocities that are proportional to kink amplitudes. For small amplitudes, the DK's of the FPU lattice reduce to the well-known ``continuous'' kink-soliton solutions of the modified Korteweg-de Vries equation. For high amplitudes, we obtain a consistent description of these DK's in terms of approximate solutions of the lattice equations that are obtained by restricting to a bounded support in space exact solutions with sinusoidal pattern characterized by the ``magic'' wavenumber $k=2\pi/3$. Relative displacement patterns, velocity versus amplitude, dispersion relation and exponential tails found in numerical simulations are shown to agree very well with analytical predictions, for both FPU and LJ lattices.Comment: Europhysics Letters (in print

Driven Macroscopic Quantum Tunneling of Ultracold Atoms in Engineered Optical Lattices

Coherent macroscopic tunneling of a Bose-Einstein condensate between two parts of an optical lattice separated by an energy barrier is theoretically investigated. We show that by a pulsewise change of the barrier height, it is possible to switch between tunneling regime and a self-trapped state of the condensate. This property of the system is explained by effectively reducing the dynamics to the nonlinear problem of a particle moving in a double square well potential. The analysis is made for both attractive and repulsive interatomic forces, and it highlights the experimental relevance of our findings

Gap soliton dynamics in an optical lattice as a parametrically driven pendulum

A long wavelength optical lattice is generated in a two-level medium by low-frequency contrapropagating beams. Then a short wave length gap soliton generated by evanescent boundary instability (supratransmission) undergoes a dynamics shown to obey the Newton equation of the parametrically driven pendulum, hence presenting extremely rich, possibly chaotic, dynamical behavior. The theory is sustained by numerical simulations and provides an efficient tool to study soliton trajectories