Statistical mechanics and pressure of composite multimoded weakly nonlinear optical systems.

Statistical mechanics can provide a versatile theoretical framework for investigating the collective dynamics of weakly nonlinear-wave settings that can be utterly complex to describe otherwise. In optics, composite systems arise due to interactions between different frequencies and polarizations. The purpose of this work is to develop a thermodynamic theory that takes into account the synergistic action of multiple components. We find that the type of the nonlinearity involved can have important implications in the thermalization process and, hence, can lead to different thermal equilibrium conditions. Importantly, we derive closed-form expressions for the actual optomechanical pressure that is exerted on the system. In particular, the total optomechanical pressure is the sum of the partial pressures due to each component. Our results can be applied to a variety of weakly nonlinear optical settings such as multimode fibers, bulk waveguides, photonic lattices, and coupled microresonators. We present two specific examples, where two colors interact in a one-waveguide array with either a cubic or quadratic nonlinearity.

Dalton's law of partial optical thermodynamic pressures in highly multimoded nonlinear photonic systems.

We show that in highly multimoded nonlinear photonic systems, the optical thermodynamic pressures emerging from different species of the optical field obey Dalton's law of partial pressures. In multimode settings, the optical thermodynamic pressure is defined as the conjugate to the extensive variable associated with the system's total number of modes and is directly related to the actual electrodynamic radiation forces exerted at the physical boundaries of the system. Here, we extend this notion to photonic configuration supporting different species of the optical field. Under thermal equilibrium conditions, we formally derive an equation that establishes a direct link between the partial thermodynamic pressures and the electrodynamic radiation pressures exerted by each polarization species. Our theoretical framework provides a straightforward approach for quantifying the total radiation pressures through the system's thermodynamic variables without invoking the Maxwell stress tensor formalism. In essence, we show that the total electrodynamic pressure in such arrangements can be obtained in an effortless manner from initial excitation conditions, thus avoiding time-consuming simulations of the utterly complex multimode dynamics. To illustrate the validity of our results, we carry out numerical simulations in multimoded nonlinear optical structures supporting two polarization species and demonstrate excellent agreement with the Maxwell stress tensor method.

Inverted pin beams for robust long-range propagation through atmospheric turbulence.

We introduce a new, to the best of our knowledge, class of optical beams, which feature a spatial profile akin to an "inverted pin." In particular, we asymptotically find that close to the axis, the transverse amplitude profile of such beams takes the form of a Bessel function with a width that gradually increases during propagation. We examine numerically the behavior of such inverted pin beams in turbulent environments as measured via the scintillation index and show that they outperform Gaussian beams (collimated and focused) as well as Bessel beams and regular pin beams, which are all optimized, especially in the moderate and strong fluctuation regimes.

Nature of Optical Thermodynamic Pressure Exerted in Highly Multimoded Nonlinear Systems.

The theory of optical thermodynamics provides a comprehensive framework that enables a self-consistent description of the intricate dynamics of nonlinear multimoded photonic systems. This theory, among others, predicts a pressurelike intensive quantity (p[over ^]) that is conjugate to the system's total number of modes (M)-its corresponding extensive variable. Yet at this point, the nature of this intensive quantity is still nebulous. In this Letter, we elucidate the physical origin of the optical thermodynamic pressure and demonstrate its dual essence. In this context, we rigorously derive an expression that splits p[over ^] into two distinct components, a term that is explicitly tied to the electrodynamic radiation pressure and a second entropic part that is responsible for the entropy change. We utilize this result to establish a formalism that simplifies the quantification of radiation pressure under nonlinear equilibrium conditions, thus eliminating the need for a tedious evaluation of the Maxwell stress tensor. Our theoretical analysis is corroborated by numerical simulations carried out in highly multimoded nonlinear optical structures. These results may provide a novel way in predicting and controlling radiation pressure processes in a variety of nonlinear electromagnetic settings.

Tunable self-similar Bessel-like beams of arbitrary order.

We predict that Bessel-like beams of arbitrary integer order can exhibit a tunable self-similar behavior (that take an invariant form under suitable stretching transformations). Specifically, by engineering the amplitude and the phase on the input plane in real space, we show that it is possible to generate higher-order vortex Bessel-like beams with fully controllable radius of the hollow core and maximum intensity during propagation. In addition, using a similar approach, we show that it is also possible to generate zeroth-order Bessel-like beams with controllable beam width and maximum intensity. Our numerical results are in excellent agreement with our theoretical predictions.

Valley Vortex States and Degeneracy Lifting via Photonic Higher-Band Excitation.

We demonstrate valley-dependent vortex generation in photonic graphene. Without breaking inversion symmetry, the excitation of two valleys leads to the formation of an optical vortex upon Bragg reflection to the third equivalent valley, with its chirality determined by the valley degree of freedom. Vortex-antivortex pairs with valley-dependent topological charge flipping are also observed and corroborated by numerical simulations. Furthermore, we develop a three-band effective Hamiltonian model to describe the dynamics of the coupled valleys and find that the commonly used two-band model is not sufficient to explain the observed vortex degeneracy lifting. Such valley-polarized vortex states arise from high-band excitation without a synthetic-field-induced gap opening. Our results from a photonic setting may provide insight for the study of valley contrasting and Berry-phase-mediated topological phenomena in other systems.

Independent amplitude and trajectory/beam-width control of nonparaxial beams.

We show that it is possible to generate non-paraxial optical beams with pre-engineered trajectories and designed maximum amplitude along these trajectories. The independent control of these two degrees of freedom is made possible by engineering both the amplitude and the phase of the optical wave on the input plane. Furthermore, we come to the elegant conclusion that the beam width depends solely on the local curvature of the trajectory. Thus, we can generate beams with pre-defined amplitude and beam-width by appropriately selecting the local curvature. Our theoretical results are in excellent agreement with numerical simulations. We discuss about methods that can be utilized to experimentally generate such beam. Our work might be useful in applications where precise beam control is important such as particle manipulation, filamentation, and micromachining.

Observation of microscale nonparaxial optical bottle beams.

We predict and experimentally observe three-dimensional microscale nonparaxial optical bottle beams based on the generation of a caustic surface under revolution. Such bottle beams exhibit high contrast between the surrounding surface and the effectively void interior. Via caustic engineering, we can precisely control the functional form of the high-intensity surface to achieve microscale bottle beams with longitudinal and transverse dimensions of the same order of magnitude. Although, in general, the phase profile at the input plane can be computed numerically, we find closed-form expressions for bottle beams with various types of surfaces both in the real and in the Fourier space.

Spatiotemporal diffraction-free pulsed beams in free-space of the Airy and Bessel type.

We investigate the dynamics of spatiotemporal optical waves with one transverse dimension obtained as the intersections of the dispersion cone with a plane. We show that, by appropriate spectral excitations, the three different types of conic sections (elliptic, parabolic, and hyperbolic) can lead to optical waves of the Bessel, Airy, and modified Bessel types, respectively. We find closed form solutions that accurately describe the wave dynamics and unveil their fundamental properties.

Nonlinear imaging in photonic lattices.

We show that nonlinear imaging is possible in periodic waveguide configurations, provided that we use two different segments of nonlinear media with opposite signs of the Kerr nonlinearity with, in general, no other restriction about their magnitudes. The second medium is used to implement effective "reverse propagation." A main ingredient in achieving nonlinear imaging is the control of the sign and the amplitude of the coupling coefficient. We numerically test our results in one- and two-dimensional square arrangements of waveguides.

Nonparaxial abruptly autofocusing beams.

We study nonparaxial autofocusing beams with pre-engineered trajectories. We consider the case of linearly polarized electric optical beams and examine their focusing properties, such as contrast, beam width, and numerical aperture. Such beams are associated with larger intensity contrasts, can focus at smaller distances, and have smaller spot sizes as compared to the paraxial regime.

Modifying the optical path in a nonlinear double-slit experiment.

In this Letter, we study a nonlinear interferometric setup based on diffraction, rather than beam combining. It consists of a nonlinear analog of Young's double-slit experiment where a nonlinear material is placed exactly after one of the slits. The presence of nonlinearity breaks the transverse spatial symmetry of the system and, thus, modifies the optical path. For moderate nonlinearities, this leads to a self-induced shift of the intensity pattern in the transverse plane. A simple theoretical model is developed which is surprisingly accurate in predicting the intensity profile of the main lobes for a wide range of parameters. We discuss possible applications of our model in nonlinear interferometry, for example in measuring the nonlinearities of optical materials.

Composite multi-vortex diffraction-free beams and van-Hove singularities in honeycomb lattices.

We find diffraction-free beams for graphene and MoS2-type honeycomb optical lattices. The resulting composite solutions have the form of multi-vortices, with spinor topological charges (n, n±1). Exact solutions for the spinor components are obtained in the Dirac limit. The effects of the valley degree of freedom and the mass are analyzed. Passing through the van Hove singularity, the topological structure of the solutions is modified. Exactly at the singularity, the diffraction-free beams take the form of strongly localized one-dimensional stripes.

Closed-form expressions for nonparaxial accelerating beams with pre-engineered trajectories.

In this Letter, we propose a general real-space method for the generation of nonparaxial accelerating beams with arbitrary predefined convex trajectories. Our results lead to closed-form expressions for the required phase at the input plane. We present such closed-form results for a variety of caustic curves: beside circular, elliptic, and parabolic, we find for the first time general power-law and exponential trajectories. Furthermore, by changing the initial amplitude, we can design different intensity profiles along the caustic.

Observation of self-accelerating Bessel-like optical beams along arbitrary trajectories.

We experimentally demonstrate self-accelerating Bessel-like optical beams propagating along arbitrary trajectories in free space. With computer-generated holography, such beams are designed to follow different controllable trajectories while their main lobe transverse profiles remain nearly invariant and symmetric. Examples include parabolic, snake-like, hyperbolic, hyperbolic secant, and even three-dimensional spiraling trajectories. The self-healing property of such beams is also demonstrated. This new class of optical beams can be considered as a hybrid between accelerating and nonaccelerating nondiffracting beams that may find a variety of applications.

Caustic design in periodic lattices.

We study curved trajectory dynamics and design in discrete array settings. We find that beams with power law phases produce curved caustics associated with the fold and cusp type catastrophes. A parabolic phase produces a focus that suffers from spherical aberrations. More important, we find that by designing the initial phase or wavefront of the beam we can construct trajectories with pure power law caustics as well as aberration-free focusing of discrete waves.

Linear and nonlinear waves in surface and wedge index potentials.

We study optical beams that are supported at the surface of a medium with a linear index potential and by a piecewise linear wedge-type potential. In the linear limit the modes are described by Airy functions. In the nonlinear regime we find families of solutions that bifurcate from the linear modes and study their stability for both self-focusing and self-defocusing Kerr nonlinearity. The total power of such nonlinear waves is finite without the need for apodization.

Bessel-like optical beams with arbitrary trajectories.

A method is proposed for generating Bessel-like optical beams with arbitrary trajectories in free space. The method involves phase-modulating an optical wavefront so that conical bundles of rays are formed whose apexes write a continuous focal curve with pre-specified shape. These ray cones have circular bases on the input plane; thus their interference results in a Bessel-like transverse field profile that propagates along the specified trajectory with a remarkably invariant main lobe. Such beams can be useful as hybrids between non-accelerating and accelerating optical waves that share diffraction-resisting and self-healing properties.

Reflection and refraction of an Airy beam at a dielectric interface.

Reflection and refraction of a finite-power Airy beam at the interface between two dielectric media are investigated analytically and numerically. The formulation takes into account the paraxial nature of the optical beams to derive convenient field evolution equations in coordinate frames moving along Snell's refraction and reflection axes. Through numerical simulations, the self-accelerating dynamics of the Airy-like refracted and reflected beams are observed. Of special interest are the cases of critical incidence at Brewster and total-internal-reflection (TIR) angles. In the former case, we find that the reflected beam achieves self-healing, despite the severe suppression of a part of its spectrum, while, in the latter case, the beam remains nearly unaffected except for the Goos-Hänchen shift. The self-accelerating quality persists even if the beam is trapped by multiple TIRs inside a dielectric film. The grazing incidence of an Airy beam at the interface between two media with close refractive indices is also investigated, revealing that the interface can act as a filter depending on the beam scale and tilt. We finally consider reverse refraction and perfect imaging of an Airy beam into a left-handed medium.

Airy trajectory engineering in dynamic linear index potentials.

We study the propagation of Airy beams in transversely linear index potentials with a gradient that is dynamically changing along the propagation direction. We find exact solutions in the case of Airy and apodized (Gaussian and exponentially) Airy beams in 1+1 and 2+1 dimensions. More important, we find that the Airy beam can follow any predefined path, in which case the potential gradient is determined as a function of this path.