Effective nonlinear optical properties of shape distributed composite media

The effective linear and nonlinear optical properties of metal/dielectric composite media, in which ellipsoidal metal inclusions are distributed in shape, are investigated. The shape distribution function P(L _x, L _y) is assumed to be 2Δ^-2θ(L _x - 1/3 + Δ/3)θ(L _y - 1/3 + Δ/3)θ(2/3 + Δ/3 - L _x - L _y), where θ( ^{. . .} ) is the Heaviside function, Δ is the shape variance and L_i are the depolarization factors of the ellipsoidal inclusions along i-symmetric axes (i = x, y). Within the spectral representation, we adopt Maxwell-Garnett type approximation to study the effect of shape variance Δ on the effective nonlinear optical properties. Numerical results show that both the effective linear optical absorption α ∼ ωIm() and the modulus of the effective third-order optical nonlinearity enhancement |χ⁽³⁾ _e|/χ⁽³⁾ ₁ exhibit the nonmonotonic behavior with Δ. Moreover, with increasing Δ, the optical absorption and the nonlinearity enhancement bands become broad, accompanied with the decrease of their peaks. The adjustment of Δ from 0 to 1 allows us to examine the crossover behavior from no separation to large separation between optical absorption and nonlinearity enhancement peaks. As Δ → 0, i.e., the ellipsoidal shape deviates slightly from the spherical one, the dependence of |χ⁽³⁾ _e|/χ⁽³⁾ ₁ on Δ becomes strong first and then weak with increasing the imaginary part of inclusions' dielectric constant. In the dilute limit, the exact formula for the effective optical nonlinearity is derived, and the present approximation characterizes the exact results better than old mean field one does. Received 10 December 2002 Published online 4 June 2003 RID="a" ID="a"e-mail: lgaophys@pub.sz.jsinfo.net     

Topological and dynamical excitations in a classical 2D easy-axis Heisenberg model

The properties of dynamical solitons (magnon droplets) in the classical, two-dimensional anisotropic Heisenberg model with easy-axis exchange anisotropy are studied. The solution of the Landau-Lifshitz equation in the continuum limit for the soliton with topological charge q = 1 is obtained numerically using a shooting method. We analized a wide range of the anisotropy parameter and our results are in good agreement with results obtained from spin dynamics simulations. The dependence of an internal precession frequency of the soliton on both the anisotropy parameter and the radius of the soliton is also investigated. Finally, the limits of applicability of the continuum approach are discussed. Received 22 August 2000     

The appearance of gap solitons in a nonlinear Schrödinger lattice

We study the appearance of discrete gap solitons in a nonlinear Schrödinger model with a periodic on-site potential that possesses a gap evacuated of plane-wave solutions in the linear limit. For finite lattices supporting an anti-phase (q=π/2) gap edge phonon as an anharmonic standing wave in the nonlinear regime, gap solitons are numerically found to emerge via pitchfork bifurcations from the gap edge. Analytically, modulational instabilities between pairs of bifurcation points on this “nonlinear gap boundary” are found in terms of critical gap widths, turning to zero in the infinite-size limit, which are associated with the birth of the localized soliton as well as discrete multisolitons in the gap. Such tunable instabilities can be of relevance in exciting soliton states in modulated arrays of nonlinear optical waveguides or Bose-Einstein condensates in periodic potentials. For lattices whose gap edge phonon only asymptotically approaches the anti-phase solution, the nonlinear gap boundary splits in a bifurcation scenario leading to the birth of the discrete gap soliton as a continuable orbit to the gap edge in the linear limit. The instability-induced dynamics of the localized soliton in the gap regime is found to thermalize according to the Gibbsian equilibrium distribution, while the spontaneous formation of persisting intrinsically localized modes (discrete breathers) from the extended out-gap soliton reveals a phase transition of the solution.     

Long-range effects on superdiffusive algebraic solitons in anharmonic chains

Studies on thermal diffusion of lattice solitons in Fermi-Pasta-Ulam (FPU)-like lattices were recently generalized to the case of dispersive long-range interactions (LRI) of the Kac-Baker form. The variance of the soliton position shows a stronger than linear time-dependence (superdiffusion) as found earlier for lattice solitons on FPU chains with nearest-neighbour interactions (NNI). Since the superdiffusion seems to be generic for nontopological solitons, we want to illuminate the role of the soliton shape on the superdiffusive mechanism. Therefore, we concentrate on an FPU-like lattice with a certain class of power-law long-range interactions where the solitons have algebraic tails instead of the exponential tails in the case of FPU-type interactions (with or without Kac-Baker LRI).Despite of structurally similar Langevin equations which hold for the soliton position and width of the two soliton types, the algebraic solitons reach the superdiffusive long-time limit with a characteristic t^3/2 time-dependence much faster than exponential solitons. The soliton shape determines the diffusion constant in the long-time limit that is approximately a factor of π smaller for algebraic solitons. Our results appear to be generic for nonlinear excitaitons in FPU-chains, because the same superdiffusive time-dependence was also observed in simulations with discrete breathers.     

Solitary waves in the Madelung's fluid: Connection between the nonlinear Schrödinger equation and the Korteweg-de Vries equation

An investigation to deepen the connection between the family of nonlinear Schr?dinger equations and the one of Korteweg-de Vries equations is carried out within the context of the Madelung's fluid picture. In particular, under suitable hypothesis for the current velocity, it is proven that the cubic nonlinear Schr?dinger equation, whose solution is a complex wave function, can be put in correspondence with the standard Korteweg-de Vries equation, is such a way that the soliton solutions of the latter are the squared modulus of the envelope soliton solution of the former. Under suitable physical hypothesis for the current velocity, this correspondence allows us to find envelope soliton solutions of the cubic nonlinear Schr?dinger equation, starting from the soliton solutions of the associated Korteweg-de Vries equation. In particular, in the case of constant current velocities, the solitary waves have the amplitude independent of the envelope velocity (which coincides with the constant current velocity). They are bright or dark envelope solitons and have a phase linearly depending both on space and on time coordinates. In the case of an arbitrarily large stationary-profile perturbation of the current velocity, envelope solitons are grey or dark and they relate the velocity u₀ with the amplitude; in fact, they exist for a limited range of velocities and have a phase nonlinearly depending on the combined variable x-u_{0 s} (s being a time-like variable). This novel method in solving the nonlinear Schr?dinger equation starting from the Korteweg-de Vries equation give new insights and represents an alternative key of reading of the dark/grey envelope solitons based on the fluid language. Moreover, a comparison between the solutions found in the present paper and the ones already known in literature is also presented. Received 20 February 2002 and Received in final form 22 April 2002 Published online 6 June 2002     

Dichromatic squeezing generation

We propose a novel scheme for the joint generation of two squeezed beams at arbitrary frequencies ω ₁ and ω ₂. The scheme consists of two successive steps, both involving nonlinear interactions in χ⁽²⁾ crystals. The dynamics of the setup is analyzed both quantum mechanically and classically within the parametric approximation. An experimental implementation involving the fundamental and the harmonics of a Nd:YAG laser pulse, and β-BaB ₂ O ₄ nonlinear crystals is suggested. Received 17 May 2000 and Received in final form 9 October 2000     

Finite-size scaling of the level compressibility at the Anderson transition

We compute the number level variance Σ ₂ and the level compressibility χ from high precision data for the Anderson model of localization and show that they can be used in order to estimate the critical properties at the metal-insulator transition by means of finite-size scaling. With N, W, and L denoting, respectively, linear system size, disorder strength, and the average number of levels in units of the mean level spacing, we find that both χ(N, W) and the integrated Σ ₂ obey finite-size scaling. The high precision data was obtained for an anisotropic three-dimensional Anderson model with disorder given by a box distribution of width W/2. We compute the critical exponent as ν≈ 1.45±0.12 and the critical disorder as W _c≈ 8.59±0.05 in agreement with previous transfer-matrix studies in the anisotropic model. Furthermore, we find χ≈ 0.28±0.06 at the metal-insulator transition in very close agreement with previous results. Received 1st November 2001 and Received in final form 8 March 2002 Published online 6 June 2002     

Short-range plasma model for intermediate spectral statistics

We propose a plasma model for spectral statistics displaying level repulsion without long-range spectral rigidity, i.e. statistics intermediate between random matrix and Poisson statistics similar to the ones found numerically at the critical point of the Anderson metal-insulator transition in disordered systems and in certain dynamical systems. The model emerges from Dysons one-dimensional gas corresponding to the eigenvalue distribution of the classical random matrix ensembles by restricting the logarithmic pair interaction to a finite number k of nearest neighbors. We calculate analytically the spacing distributions and the two-level statistics. In particular we show that the number variance has the asymptotic form Σ²(L) ∼χL for large L and the nearest-neighbor distribution decreases exponentially when s→∞, P(s) ∼ exp(- Λs) with Λ = 1/χ = kβ + 1, where β is the inverse temperature of the gas (β = 1, 2 and 4 for the orthogonal, unitary and symplectic symmetry class respectively). In the simplest case of k = β = 1, the model leads to the so-called Semi-Poisson statistics characterized by particular simple correlation functions e.g. P(s) = 4s exp(- 2s). Furthermore we investigate the spectral statistics of several pseudointegrable quantum billiards numerically and compare them to the Semi-Poisson statistics. Received 13 September 2000     

Nonlinear frequency conversion in 2D χ (2) photonic crystals and novel nonlinear double-circle construction

The analytic solution to the wave equation for small-signal sum-frequency process is derived in 2D χ ⁽²⁾ photonic crystals with use of the Green function method. It is predicted that the sum-frequency electrical field at quasi-phase matching (QPM) resonance is proportional to the angle-dependent effective crystal length. This implies that multiple wavelength QPM frequency conversion with controllable intensity output can be realized in a single 2D χ ⁽²⁾ photonic crystal. It is revealed that efficient frequency conversion requires both the QPM and the proper structure matching. A novel double-circle construction, different from the conventional Ewald construction, is presented to reflect important QPM processes. It is also shown that the QPM resonance tuning of second-harmonic generation can operate over the whole transparent wavelength range of crystals. Received 19 April 2001     

The lifetime of the soliton in the improved Davydov model at the biological temperature 300 K for protein molecules

We study the effects of quantum fluctuations and thermal perturbations on the lifetime of the soliton in the improved Davydov model proposed by us with two-quanta and with an added interaction. By using quantum perturbation theory, we compute the soliton lifetime for a wide ranges of parameter values relevant for protein molecules. The lifetime of the new soliton at the biological temperature 300 K is of the order of 10^-10 second or τ/τ≥ 500 for parameters appropriate to α-helical protein molecules. This shows clearly that the new soliton in the improved model is a viable mechanism for the bio-energy transport in the α-helix region of proteins. Received 7 January 1999 and Received in final form 16 August 2000     

Weierstrass \wp -lumps

We present our results of a numerical investigation of the behaviour of a system of two solitons in the (2+1) dimensional CP¹ model on a torus. Defined by the elliptic function of Weierstrass, and working in the Skyrme version of the model, the soliton lumps exhibit splitting, scattering at right angles and motion reversal in the various configurations considered. The work is restricted to systems with no initial velocity. Received 24 April 2001     

Vibrational wave packet dynamics in the silver tetramer probed by NeNePo femtosecond pump-probe spectroscopy

We present results on the ultrafast dynamics of mass-selected neutral Ag₄ clusters using NeNePo (negative ion - neutral - positive ion) femtosecond pump-probe spectroscopy. One-color pump-probe spectra of the Ag₄ ^-/Ag₄/Ag₄ ⁺ system measured at 385 nm and an internal cluster temperature of 20 K display a complex beat structure over more than 60 ps. The oscillatory structure is attributed to vibrational wave packet dynamics in an excited “dark" state of neutral Ag₄. A dominant 740 fs wave packet period as well as wave packet dephasing and rephasing are observed in the spectra. Fourier analysis of the spectra yields a group of frequencies centered around 45 cm^-1 and an anharmonicity χ _e2νχ _eχ _e of 2.65 cm^-1 for the active vibrational mode. Received 30 November 2000     

Neutron halos in hypernuclei

Properties of single- Λ and double- Λ hypernuclei for even-N Ca isotopes ranging from the proton dripline to the neutron dripline are studied using the relativistic continuum Hartree-Bogoliubov theory with a zero-range pairing interaction. Compared with ordinary nuclei, the addition of one or two Λ-hyperons lowers the Fermi level. The predicted neutron dripline nuclei are, respectively, ⁷⁵ _ΛCa and ⁷⁶ _2ΛCa, as the additional attractive force provided by the Λ-N interaction shifts nuclei from outside to inside the dripline. Therefore, the last bound hypernuclei have two more neutrons than the corresponding ordinary nuclei. Based on the analysis of two-neutron separation energies, neutron single-particle energy levels, the contribution of continuum and nucleon density distribution, giant halo phenomena due to the pairing correlation, and the contribution from the continuum are suggested to exist in Ca hypernuclei similar to those that appear in ordinary Ca isotopes. Received: 21 October 2002 / Accepted: 11 January 2003 / Published online: 8 April 2003     

Spatial solitons in photonic lattices with large-scale defects

We address the existence,stability and propagation dynamics of solitons supported by large-scale defects surrounded by the harmonic photonic lattices imprinted in the defocusing saturable nonlinear medium.Several families of soliton solutions,including flat-topped,dipole-like,and multipole-like solitons,can be supported by the defected lattices with different heights of defects.The width of existence domain of solitons is determined solely by the saturable parameter.The existence domains of various types of solitons can be shifted by the variations of defect size,lattice depth and soliton order.Solitons in the model are stable in a wide parameter window,provided that the propagation constant exceeds a critical value,which is in sharp contrast to the case where the soliton trains is supported by periodic lattices imprinted in defocusing saturable nonlinear medium.We also find stable solitons in the semi-infinite gap which rarely occur in the defocusing media.     

Analytic doubly periodic wave patterns for the integrable discrete nonlinear Schrödinger (Ablowitz–Ladik) model

We derive two new solutions in terms of elliptic functions, one for the dark and one for the bright soliton regime, for the semi-discrete cubic nonlinear Schrödinger equation of Ablowitz and Ladik. When considered in the complex plane, these two solutions are identical. In the continuum limit, they reduce to known elliptic function solutions. In the long wave limit, the dark one reduces to the collision of two discrete dark solitons, and the bright one to a discrete breather.     

Cascaded χ(2) nonlinearity in QND measurement

 QND measurement schemes often use Kerr nonlinearity to couple the intensity fluctuations of the signal beam to the phase fluctuations of the probe beam by means of the cross-phase modulation effect. Such schemes use materials with intrinsic χ⁽³⁾nonlinearity. However, it has already been shown that Kerr-effect-like correlation between two light waves may be achieved through cascaded χ⁽²⁾: χ⁽²⁾ processes. The value of the nonlinearity induced by cascading χ⁽²⁾ can, in many cases, be much higher than fibre χ⁽³⁾. The authors point out the possibility of using the cascaded χ⁽²⁾ nonlinearity instead of intrinsic χ⁽³⁾ nonlinearity in the scheme which resembles the well-known photon number QND measurement scheme. The non-degenerate second-harmonic generation (SHG) with a large phase mismatch is considered. The harmonic wave influence on the measurements is shown. High QND measurement characteristics can be achieved with the use of an appropriate set of parameters. Received: 6 April 1996     

Slow relaxation experiments in disordered charge and spin density waves: collective dynamics of randomly distributed solitons

We show that the dynamics of disordered charge density waves (CDWs) and spin density waves (SDWs) is a collective phenomenon. The very low temperature specific heat relaxation experiments are characterized by: (i) “interrupted” ageing (meaning that there is a maximal relaxation time); and (ii) a broad power-law spectrum of relaxation times which is the signature of a collective phenomenon. We propose a random energy model that can reproduce these two observations and from which it is possible to obtain an estimate of the glass cross-over temperature (typically T _g≃ 100-200 mK). The broad relaxation time spectrum can also be obtained from the solutions of two microscopic models involving randomly distributed solitons. The collective behavior is similar to domain growth dynamics in the presence of disorder and can be described by the dynamical renormalization group that was proposed recently for the one dimensional random field Ising model [D.S. Fisher, P. Le Doussal, C. Monthus, Phys. Rev. Lett. 80, 3539 (1998)]. The typical relaxation time scales like ∼τexp(T _g/T). The glass cross-over temperature T_g related to correlations among solitons is equal to the average energy barrier and scales like T _g∼ 2xξΔ. x is the concentration of defects, ξ the correlation length of the CDW or SDW and Δ the charge or spin gap. Received 12 December 2001     

Self-trapped spatially localized states in combined linear-nonlinear periodic potentials

We analyze the existence and stability of two kinds of self-trapped spatially localized gap modes,gap solitons and truncated nonlinear Bloch waves,in one-and two-dimensional optical or matter-wave media with self-focusing nonlinearity,supported by a combination of linear and nonlinear periodic lattice potentials.The former is found to be stable once placed inside a single well of the nonlinear lattice,it is unstable otherwise.Contrary to the case with constant self-focusing nonlinearity,where the latter solution is always unstable,here,we demonstrate that it nevertheless can be stabilized by the nonlinear lattice since the model under consideration combines the unique properties of both the linear and nonlinear lattices.The practical possibilities for experimental realization of the predicted solutions are also discussed.     

Nonlinear surface magneto-optics of ferromagnetic Ni/Cu(001) from first principles

We present first-principles calculations of the nonlinear optical (NLO) response of a Ni/Cu(001) bilayer. The calculations are based on the full potential linearized augmented plane wave (FLAPW) method with the additional implementation of spin–orbit coupling (SOC). On the basis of this set of eigenstates the magneto-optical transition-matrix elements are evaluated. Using the surface-sheet model the optical reflection properties are determined for the cases of the magnetization vector perpendicular to the surface (polar magneto-optical configuration (MOC)) and for the in-plane magnetization (longitudinal MOC). The nonlinear optical susceptibility tensor elements χ⁽²⁾ _ijk for different magnetization directions as well as the spectral dependence of χ⁽²⁾ _ijk, the resulting intensities, and Kerr angles are presented for the Ni/Cu(001) bilayer. The results show that the magnetic tensor elements of the χ⁽²⁾ _ijk tensor are smaller than the nonmagnetic ones by only one order of magnitude, confirming the important role of magnetic properties in the NLO response. Received: 16 October 2001 / Revised version: 8 March 2002 / Published online: 29 May 2002     

二阶孤子间相互作用引起孤子衰变

对于二阶孤子脉冲对,通过数值求解非线性Schrdinger方程,得到了二阶孤子在光纤中传输时的衰变规律.由于脉冲之间的相互作用,原来的二阶孤子衰变为两个分离的基孤子.二阶孤子对同位相时,相互作用为“吸引”,衰变后峰值小的基孤子位于内侧,二阶孤子对反位相时,相互作用为“排斥”,衰变后峰值小的基孤子位于外侧.两个二阶孤子脉冲间隔越小,反位相时,在相同传输距离处,衰变产生的两个基孤子的分离程度越大.