#PAGE_PARAMS# #ADS_HEAD_SCRIPTS# #MICRODATA#

Generalized nonlinear Schrödinger equations describing the Second Harmonic Generation of femtosecond pulse, containing a few cycles, and their integrals of motion


Authors: Vyacheslav A. Trofimov aff001;  Svetlana Stepanenko aff001;  Alexander Razgulin aff001
Authors place of work: Faculty of Computational Mathematics and Cybernetics, Lomonosov Moscow State University, Moscow, Russia aff001
Published in the journal: PLoS ONE 14(12)
Category: Research Article
doi: https://doi.org/10.1371/journal.pone.0226119

Summary

An interaction of laser pulse, containing a few cycles, with substance is a modern problem, attracting attention of many researches. The frequency conversion is a key problem for a generation of such pulses in various ranges of frequencies. Adequate description of such pulse interaction with a medium is based on a slowly evolving wave approximation (SEWA), which has been proposed earlier for a description of propagation of the laser pulse, containing a few cycles, in a medium with cubic nonlinear response. Despite widely applicability of the frequency conversion for various nonlinear optics problems solutions, SEWA has not been applied and developed for a theoretical investigation of the frequency doubling process until present time. In this study the set of generalized nonlinear Schrödinger equations describing a second harmonic generation of the super-short femtosecond pulse is derived. The equations set contains terms, describing the pulses self-steepening, and the second order dispersion (SOD) of the pulse, a diffraction of the beam as well as mixed derivatives. We propose the transform of the equations set to a type, which does not contain both the mixed derivatives and time derivatives of the nonlinear terms. This transform allows us to derive the integrals of motion of the problem: energy, spectral invariants and Hamiltonian. We show the existence of two specific frequencies (singularities in the Fourier space) inherent to the problem. They may cause an appearance of non-physical absolute instability of the problem solution if the spectral invariants are not taken into account. Moreover, we claim that the energy preservation at the laser pulses propagation may not occur if these invariants do not preserve. Developed conservation laws, in particular, have to be used for developing of the conservative finite-difference schemes, preserving the conservation laws difference analogues, and for developing of adequate theory of the modulation instability of the laser pulses, containing a few cycles.

Keywords:

Modulation – Lasers – Laser beams – Optics – Electric field – Conservation of energy – Diffraction – Computer modeling

Introduction

A well-known process of second harmonic generation (SHG) by the laser beam was the first nonlinear optical effect, which was observed by P. Franken et al. [1] in 1961. One year later, J. Armstrong, N. Bloembergen et al. [2] published the fundamental paper, dealing with the optical frequency conversion, in which the coupled nonlinear equations were developed and many explicit solutions of these equations were obtained in the framework of the plane wave approximation. During passed decades the SHG was widely observed: in plasma optics [35] and nonlinear optics, including high-harmonic generation [612], and in a medium containing nanoparticles [13, 14], and in semiconductor [15]. In [16] the SHG at the boundary between dielectric media is analyzed. Obviously, this list of papers can be easily increased.

Many researchers pay an attention to an investigation of various factors, which restrict the frequency doubling efficiency. For example, in [17] the effect of a group-velocity dispersion was considered. In [18] the authors have considered an influence of interplaying between two types polarizations of the pulse at the fundamental frequency on the frequency doubling efficiency. With increasing of the laser pulse intensity it is necessary to take into account a cubic nonlinear response of a medium under the SHG. An influence of the self-modulation as well as cross-modulation of the laser pulse on the frequency doubling efficiency was investigated experimentally and theoretically (see, for example [2, 1928]). It should be stressed that the most general analytical theory of this process has been developed recently in [28], taking into account the time-dependent pulse shape and without using the basic wave energy non-depletion approximation on the base of using the integrals of motion (conservation laws) of the problem. We notice that the authors of [22] have obtained the most (but not all) solutions of this problem in the framework of long pulse duration approximation. In [29] an effective way for increasing of SHG efficiency in bulk medium was demonstrated by using of the incident beam with a ring (tubular) profile. This leads to re-profiling of the beams due to their diffraction and, therefore, to decreasing the beam phase distortions.

New opportunities appeared at SHG under the big phase-mismatching, so-called cascading SHG [3033]. This process allows us to achieve an effective cubic nonlinear response in a medium with a quadratic nonlinear response. This effect was firstly predicted by Yu. N. Karamzin and A. P. Sukhorukov [34] in 1974. Later, it was used for a compression of the soliton with femtosecond duration [30, 3450], and for the beam self-focusing [5153], and for using the cascading SHG for suppression of the optical pulse intensity fluctuations, occurring in a sequence of pulses with randomly variations of their maximal intensities [5457]. Such sequence of the laser pulses is produced by a laser, operating at the free generation mode. The total duration of the pulses sequences is equal to milliseconds. Therefore, the cascading SHG allows realizing the laser generation mode that is similar to the Kerr-locked mode applying for the femtosecond laser system.

There are a lot of other schemes of SHG, using the modern substances (see, for example, [5884]) and at a present time the problem is still actual and important for various applications. Among these applications, we stress using of SHG for a measurement of pulse parameters and of parameters of a medium as well as for visualization of various fast processes, occurring in the substances. The most famous approaches of parameters measuring are SPIDER (spectral interferometry for direct electric-field reconstruction) [85] and FROG (frequency-resolved optical gating) [86, 87]. Also, in [88], a technique, based on combination of FROG spectra measurements was proposed to completely characterize the amplitude and phase of an ultrashort pulse in space and time. Mix FROG method was implemented in [89]. Thus, the frequency doubling of the optical pulse, containing a few cycles, is of practical interest until present time. However, in this case it is necessary to take into account the self-steepening of the laser pulses on both fundamental and doubled frequencies. According to our knowledge, the first analysis of this problem was made numerically in [90] at the group velocity synchronism and phase matching, without taking into account the beam diffraction and the second order dispersion influence on the pulses interaction. For writing the corresponding equations, the authors have applied an approach, which was developed first by N. Tzoar and M. Jain [91] under the investigation of the single pulse propagation in a medium with a cubic nonlinear response. Experimental observation of the self-steepening influence on the pulse spectrum deformation was demonstrated in [92]. Then, this approach was successfully used for describing the processes of a laser pulse propagation in the optical fiber [93]. The corresponding nonlinear Schrödinger equation (NLSE) contains the first time derivative of a nonlinear medium response. This equation was derived using the slowly varying envelope approximation (SVEA) and named as generalized NLSE.

Taking into account a laser beam diffraction, T. Brabec and F. Krausz proposed in 1997 a new approach, so-called slowly evolving wave approximation (SEWA). Derived equation was also named as the generalized nonlinear Schrödinger equation (GNLSE) [94]. Let us note that the SEWA requires not only an envelope of a wave packet, but also the phase changing must slowly vary as the pulse covers a distance, which is equal to the wavelength. The main feature of GNLSE consists in its containing of mixed derivatives on time and spatial coordinates, and of the second-order time derivative of a medium nonlinear response (which describes the pulse self-steepening).

It should be emphasized that a differential operator of the GNLSE, written in [91, 93], coincides with a differential operator of the GNLSE, derived in [94] if the optical beam diffraction does not take into account. They differ only by the factor two at the coefficient characterizing the self-steepening term of the equation. Obviously, these equations written in the dimensionless units look the same.

Finding of conservation laws attracts an attention of various authors (see, for example, recently published papers [9599]). In papers [100, 101] we derived the conservation laws for the frequency doubling process if the self-steepening of the pulses occurs. However, until present time, the set of GNLSEs describing a SHG of super-short femtosecond pulse with taking into account a diffraction of the laser beam (it means that this process is described in the framework of the SEWA) as well as their conservation laws were not derived. We do this in the present paper.

As is well-known, the invariants are very useful for computer simulation of the laser pulse propagation because the equations describing various problems of nonlinear optics are nonlinear ones and as a rule they do not allow us to find their analytical solution. As a consequence, the finite-difference time-domain (FDTD) method [102, 103] for computer simulation of a propagation of the optical pulse, containing a few cycles, is usually employed. For example, in [104] the authors develop the time-transformation method based on FDTD. An alternative method, which permits applying fast numerical algorithms, is a generalized source method (GSM) [105]. Linear GSM was adapted to describe the SHG in diffraction gratings, containing non-centrosymmetric materials [106]. Obviously, there are other papers, in which a computer simulation is applied for an investigation of the SHG problem. As a rule, a split-step method was used at computer simulation.

Another approach for computer simulation is based on developing of the conservative finite-difference schemes [107]. They allow preserving the difference analogues of the corresponding conservation laws. Obviously, it is necessary firstly to write these conservation laws. With this aim, we propose a new transform for the derived GNLSEs, which reduces the equations to the form containing neither mixed derivatives on time and spatial coordinates nor time derivatives of the nonlinear response of a medium. Using a new set of GNLSEs, we obtain the energy invariant, spectral invariants and Hamiltonian. We add some specific requirements to the problem statement (excluding of the singularities in frequency space) to avoid a development of non-physical instability of the SHG process at taking into account the pulses self-steepening. These requirements, together with the spectral invariants, are also important for a validity of the energy conservation law.

Generalized nonlinear Schrödinger equations derivation

Derivation of the equations for SHG process in the framework of SEWA

Derivation of the equations describing the optical frequency doubling in the framework of SEWA for an electric field strength E(z, x, t) starts from the well-known wave equation:


written above in 2D case, for example. The laser pulse propagates along z coordinate, t is a time, c is the light velocity in a vacuum, x is a transverse coordinate to the optical pulse propagation direction, an operator ∇ ⊥ 2 = ∂ x 2 is the transversal Laplace operator. We choose only x-coordinate as a transverse one for simplicity. This does not restrict our consideration. The function D(z, x, t) describes an electric field induction and is written in the form:

the function P(z, x, t) is a medium polarization, containing its linear and nonlinear responses:

For linear non-instantaneous medium response Plin(z, x, t) the following representation


is widely used, which yields in the relation for a dielectric permittivity:

here χ ˜ ( 1 ) ( ω ) is a linear electric susceptibility of a medium at the frequency ω.

Function Pnl(z, x, t) describes a nonlinear response of a medium. If we consider a medium with the quadratic nonlinearity, then this function is defined as Pnl = χ(2)E2, and χ(2) is a quadratic electric susceptibility.

The next step for the wave equation reduction consists in representation of an electric field strength E(z, x, t) in a following manner:


where E1(z, x, t), E2(z, x, t) are the electric field strengths of waves, propagating with carrier frequency ωo, and with doubled frequency 2ωo and with wave-numbers k1, k2, correspondingly. In (6) we introduce the slowly varying envelopes for both wave packets and denote them as A1(z, x, t) and A2(z, x, t). Symbols c.c. denote a complex conjugation. The linear polarization of a medium at these frequencies can be written as

The corresponding nonlinear polarization of a medium at chosen frequencies is written as follows


Let us remind that the following relations for wave-numbers are valid:


where nω, n2ω are the medium refractive indexes at the frequencies ωo and 2ωo, correspondingly.

Taking into account the expressions (6) and (8) and (1) transforms to the following form:


where Dj(z, x, t), j = 1, 2 is a linear part of the electric field induction at the frequencies ωo and 2ωo, correspondingly:

Let us multiply (9) by exp(− ik1z) or exp(− ik2z). Then, we integrate each of these expressions over the corresponding wave period, taking into account the orthogonal property of the sine and cosine functions. Then, (9) can be reduced to the following equations:


Let us note the Fourier transform property:


Therefore, (11) can be written as:


Taking into account the dispersion relation:


one can write (13) as follows


Normally, we expand k(ω) in a series near the frequency ωo in (15):


and near the frequency 2ωo in (16):

here βn and sm are the n–th or m–th order derivatives of the wave number at the frequencies ωo and 2ωo, correspondingly.

Let us neglect the terms βn, sm, n, m = 3, 4, … referring to the third order dispersion and other higher terms in (17) and (18). Also, we neglect the term (ωωo)4. Thus, (17) and (18) take the forms:



After Inverse Fourier transform, we obtain the following equations set for describing the frequency doubling process:



Above Δk = k2 − 2k1 is the phase mismatch. (21) and (22) are called as the set of generalized nonlinear Schrödinger equations (GNLSEs).

Coordinate system transform

First type of coordinate system transform

In a coordinate system moving with the first wave packet


(21) and (22) could be re-written as follows:


Let us raise to the second power the expressions in brackets and then neglect the terms with the second order derivative on ξ − coordinate:



By grouping the terms in (26) and (27) and dividing them by 2ikj, j = 1, 2, we obtain:



The set of Eqs (28) and (29) describes the SHG process for the pulse with super-short duration, which is about a few femtoseconds.

Further reduction of these equations is a consequence of the following assumptions.

  • The group velocities of pulses and their phase velocities differ insignificantly:



  • Difference between inverse values of the phase velocities of waves is many times less than unity:


As a result, (28) written with respect to the wave for basic frequency ωo, reduces to


and (29) takes the following form:

Second type of coordinate system transform

In a coordinate system moving with overage velocity


the set of Eqs (21) and (22) could be re-written as follows:


If the assumptions (30)(32) are valid, then (36), with respect to basic frequency ωo, reduces to


and (37), with respect to doubled frequency 2ωo. takes the following form:

Below we use the First type of Coordinate system transform.

Dimensionless variables and problem statement

Let us introduce dimensionless variables:



Above the parameter τp is a duration of the incident pulse for basic frequency ωo, Io is its maximal intensity, the parameter a is a beam radius.

In new variables, the GNLSE set (33) and (34) takes the following form:


which is considered in the following domain

with the boundary conditions (BCs) for the functions Aj(ξ, x, t), j = 1, 2:

or

and with the initial conditions

Above γ is a parameter, which is inversely proportional to the doubled frequency 2ωo and the pulse duration τp:


The parameter ν describes group-velocity mismatch (GVM) normalized by pulse duration τp:


Parameters D21, D22 are equal to a ratio between the diffraction length of the beam and a dispersion length of the pulse for each wave:



They characterize the second order dispersion of a wave packet. Parameter D is equal to unity for chosen normalization of spatial coordinate along which the laser pulse propagates. Parameter α describes the nonlinear coupling of waves and is defined in the following form


which is written using the following assumption:

which does not restrict our consideration. Let us note that if the equality (51) is not valid, then one can transform the equations with respect to new variables, in which we obtain only single coefficient (it characterizes the nonlinear coupling of waves for both equations, see [21]).

Usually, at theoretical analysis of the problem (42), the exponentially decaying functions Aj(ξ, x, t), j = 1, 2 are used as the initial laser pulse distributions:


and instead of the BCs (43) and (44) the set of equations (42) is solved in the following unbounded (increasing of |x|, |t| and z) domain

Below, for definiteness, we use BCs (43). In addition, for convenience, below we use the previous notation for a longitudinal coordinate: z instead of ξ.

Equations transform

For transforming the equations set (42) to the form, which does not contain a time derivative of the nonlinear response of medium and mixed derivatives on time and spatial coordinate, we use the following equalities:


Also, due to a finiteness of the initial distributions of complex amplitudes and boundedness of the laser pulse propagation distance, let us state the following additional conditions:


Using the transforms (53) and condition (54), GNLSE set (42) can be written as follows (see detailed derivation in Appendix [A]):



It should be stressed that these equations will be also used for writing the conservation laws.

For further reduction of the equations, let us introduce the functions:


or

Keeping in mind the relations (53) and (58), the set of Eqs (55) and (56) can be written as follows



Further, let us introduce the new functions Ej(z, x, t), j = 1, 2 in accordance with a rule:


and the functions E ˜ j ( z , x , t ), j = 1, 2 as

Then, (59) and (60) can be written in the following manner:



These equations should be solved together with the relaxation equations with respect to the functions Ej(z, x, t), E ˜ j ( z , x , t ), Aj(z, x, t):



Taking into account the BCs (43) or (43′), and initial conditions (45) for the functions Aj, j = 1, 2, we obtain the following BCs with respect to the functions Ej, E ˜ j on time coordinate:



As consequence of BCs (43), we obtain the following conditions:



if we solve the problem with BCs (43′).

Obviously, the BCs for these functions on x-coordinate are


The initial conditions for the functions Ej, E ˜ j, j = 1, 2 are



All conditions written above will be used at the construction of the conservation laws.

Problem invariants

Using new variables, some conservation laws can be derived.

The conservation law of energy

Theorem 1. The problem (42) with BCs (43) and initial conditions (45) possesses the conservation law of energy


Proof. Let us suppose that the functions Aj satisfy the conditions (43) and (54). Then, we multiply (55) and (56) by A j *, and the equations, conjugated to (55) and (56),—by Aj. Then, we sum the obtained equations and integrate the resulted equation with respect to x, t coordinates in the domain Ωo. Thus, we obtain the following equation in Appendix [B]:


For the validity of the invariant (73) one has to show that the last integrals in (74) are equal to zero:


In fact, from (55) and (56) at the time moment t = Lt and from the BCs (54), it follows that


This yields to the relations


It means that the Fourier harmonic amplitudes at the frequencies 1 γ and 1 2 γ vary with x-coordinate growing. This property occurs even if we consider the laser pulse linear propagation (in (42) the parameter α equals zero). Obviously, the physical reason of such harmonic amplitude evolution is absent. Therefore, the functions C1j(z), C2j(z) have to be equal zero. If C1j(z) will be equal to nonzero constant then the energy of the laser beam will increase with a propagation distance. Obviously, a physical reason of this is absent for the case under consideration. Consequently, the amplitudes of Fourier harmonics at the frequencies 1 γ and 1 2 γ must be equal zero:


Thus, the invariant (73) is valid. In a case of the laser pulse nonlinear interaction (α ≠ 0), the same conditions take place.

Corollary 1. The Fourier harmonics at the frequencies 1 γ and 1 2 γ must be absent in the incident pulses spectra the energy conservation law preservation for the set of equations (42):


Let us note that if the parameter γ tends to zero, then the set of equations (42) reduces to the set of the NLSEs and the energy invariant (73) is valid. If the spectral harmonics 1 γ or 1 2 γ are absent in the incident pulse spectrum, then these spectral harmonics will be absent in the pulse spectrum computed in any section of a medium. This statement validity is stated by the spectral invariants, which are formulated below.

Spectral invariants

Spectral invariants describe an evolution of spectral harmonic amplitudes at the frequencies 1 γ and 1 2 γ along z-coordinate and show that the amplitudes should not increase during the laser pulses interaction. For writing of the Spectral invariants, let us introduce the following additional conditions:


Actually, in laser physics the propagation distance is boundedness and the domain under consideration can be chosen in such a way that the additional conditions are valid in this domain. We note that these conditions are important only for the invariant obtaining and they are not mandatory for the problem statement. The problem could be solved without additional conditions. Obviously, instead of BCs (80) one can use BCs (52). We stress that the laser sources generate the finite beam distribution in spatial coordinates. Therefore, the statements mentioned above are valid with respect to the domain choice.

Theorem 2. The problem (42) with BCs (43) and additional conditions (78) and (80), and initial conditions (45) possesses the spectral invariants



Proof. Let us consider (63) at a time moment t = Lt and integrate them with respect to x-coordinate:



Taking into account the relaxation Eqs (64) and (65) at time moment t = Lt



and additional conditions (78) we obtain:

It is easy to see that the third terms in (83) and (84) are equal to zero:


The last integral (88) equals zero due to BCs (43). Thus, we obtain the following problems:



with the initial conditions


Therefore, the invariants (81) and (82) can be obtained as the solutions of the problems (89)(92).

Hamiltonian of the problem

Let us use the following substitution for the function A2(z, x, t):


Then, the functions E2(z, x, t), E ˜ 2 ( z , x , t ) can be replaced as follows:



Using the substitution (93)(95), the set of equations (42) can be written in the following form:


which does not contain the terms with exp(±iΔkz). Consequently, with respect to the substitutions (94) and (95), the set of equations (96) takes the form:


Theorem 3. The problem (42) with BCs (43) and additional conditions (78), and initial conditions (45) possesses a Hamiltonian


Proof. Let us multiply (63) by ∂ A * ∂ z and the equation, conjugated to it,—by ∂ A ∂ z. Then, we sum the obtained equations and integrate the resulting expression with respect to x, t coordinates in the domain Ωo under consideration. Then, we write an equation in the form of the sum of integrals:


Analysis of the integrals is made in Appendix [C]. On its base, one can obtain the Hamiltonian of the problem under consideration.

Corollary 2. Using the inverse transforms for (93)(95) the Hamiltonian takes the following form:


Discussion

A few words about applicable range of the laser pulse parameters and a medium at which it is necessary to use the GNLSEs for a description of the process under consideration. It should be noticed that the first experiment, which demonstrated self-steepening of one pulse, was made in paper [91] for the pulse with picosecond duration propagating in optical fiber about 5 km long. The pulse, propagated this distance, possessed non-symmetrical spectrum. Therefore, a key role for application of GNLSEs plays a relation between the dispersion length of the pulse and the laser pulse propagation distance. However, the pulse self-steepening appears at short length (about a few centimeters) of a medium if the pulse duration is about 10-50 fs.

Another case of using GNLSE for the laser pulse propagation description corresponds to falling of incident laser pulse with sharp intensity distribution on a medium because an influence of a time derivative of the nonlinear response of the medium enhances many times in comparison with the Gaussian pulse propagation. However, contrast temporal nonlinear response can be induced by the laser pulse if an optical bistability occurs. As is well-known, in this case the explosive changing of the characteristics of a medium occurs. As a result, the temporal structure with strong gradient appears. Therefore, a time derivative of the nonlinear response of a medium increases its influence on the phase modulation.

A detail analysis of SHG describing in the framework of the GNLSEs has to be made using computer simulation. Nevertheless, one can suppose some characteristics features of the frequency conversion in such conditions. First, the group velocity of each of wave packets will depend on complex amplitude of another wave packet. Second, obviously that the pulse spectra will distort and become non-symmetrical. Third, the nonlinear length, that defines the distance at which the energy of base wave transfers to the energy of the wave with doubled frequency, will differ for front of the pulse and its trailing part. Fourth, because of the presence of mixed derivatives, the optical beam diffraction will depend on time moment of propagating pulse.

Conclusions

In the framework of the SEWA approximation, we derived the GNLSEs, describing the SHG in a medium with a quadratic nonlinear response for the pulses, containing a few cycles. The main feature of the equations set concludes in a presence of the second order time derivative of the nonlinear response of a medium (dispersion of a medium nonlinear response) as well as the mixed derivatives on time and spatial coordinate.

We proposed the equations transform, which reduced these equations to the other ones, containing neither mixed derivatives of complex amplitudes nor time derivatives of the nonlinear response. These equations are more convenient for the computer simulation and for theoretical analysis of the frequency doubling process of the optical pulses, containing a few cycles.

Based on this transform we derived some conservation laws (invariants) for the SHG problem. We showed an existence of two specific frequencies (singularities in the Fourier space) inherent to the problem. They may cause an appearance of non-physical absolute instability of the problem solution if the spectral invariants are not taken into account. It should be mentioned that a presence of such singularities was also discussed in [32] at analysis of the modulation instability occurring in χ(2)-medium. We showed that the energy conservation law is valid at certain conditions on the incident pulse spectra: the spectra must not contain non-zero spectral amplitudes at two specific frequencies. Their zero-value amplitudes at the pulse propagation are provided by the spectral invariants. We derived also the Hamiltonian (the third invariant) of the problem.

All invariants mentioned above should be taken into account at least for developing of the conservative finite-difference schemes at computer simulation of the problem under consideration.

Appendix

A Equation transform derivation

Taking into account the transforms (53), one can write


– for the first equation of the set (42), and

– for the second equation of the set (42).

Then, multiplying (101) by exp ( - i t 2 γ ) ( 2 i γ ) - 1, and (102) by exp ( - i t γ ) ( i γ ) - 1, and integrating them with respect to t-coordinate, we obtain:



correspondingly.

Due to a boundedness of the initial distribution and boundedness of the laser pulse propagation distance, let us state the following additional type of the conditions:


Then, (103) and (104) become the following ones:



which contain the first derivative of the quadratic nonlinearity. Multiplying (106) by exp ( i t 2 γ ), and (107) by exp ( i t γ ), we obtain the following equations:


Let us stress that these equations are used also for writing of conservation laws.

B Energy invariant derivation

Using (55) and (56) we obtain the equations for the functions Aj and A j *, j = 1, 2, which can be re-written as a sum of integrals:


Below we discuss transforms of the integrals mentioned above. The first one contains the exact differential of functions A1, A2 intensities. Therefore, it is easy to see that the integral I1 transforms to the following expression:


Also, the integral I2 contains the time derivative of function A2 intensity:


which equals zero, due to BCs (43). As a consequence of the integration by parts, the equalities

take place under the conditions (43) validity.

Let us show, that the integral of the nonlinear terms I5 equals zero. The terms in integral I5 without the parameter γ cancel each other and integral I5 can be written as


For the integrals I41, I42, taking into account the functions Pj, j = 1, 2 definitions (57) and the BC (43), we obtain:


Thus, (110) transforms to the following equation:


C Hamiltonian derivation

As mentioned in (99), we obtain a sum of the integrals. Let us transform them. First, taking into account (64) and (78) and integrating of the I1j by parts we obtain:


The integrals I2j are transformed into the following expression:


It is easy to see that the integrals J21j and J22j are transformed into the following integrals:



Thus, the integrals I2j take the form:


The integrals I3j can be written in the following way:


where

and

Let us note that multiplying equation, conjugated to (63), by ∂ E ˜ j ∂ z, and taking into account the BCs (43) and additional conditions (78), we obtain, that the second term in (123) equals zero. Thus,


The integrals J33j can be written as follows


Also, the integral I4, by using (64) and integrating by parts, reduces to the following one:


The integrals in I5 can be written as


Using the condition (78) and integrating by parts, one can transform the integrals I51, I52 to the form


and

Obviously, the integral I6 can be written as follows:


Substituting the expressions (116), (120), (121), (126), (127) and (130) into (99), we obtain the invariant (98).


Zdroje

1. Franken PA, Hill AE, Peters CW, Weinreich G. Generation of optical harmonics. Phys. Rev. Lett. 1961; 7(4): 118–120. doi: 10.1103/PhysRevLett.7.118

2. Armstrong J, Bloembergen N, Ducuing J, Pershan P. Interactions between Light Waves in a Nonlinear Dielectric. Phys. Rev. 1962; 127(6): 1918–39. doi: 10.1103/PhysRev.127.1918

3. Linde D, Schulz H, Engers T, Schüler H. Second Harmonic Generation in Plasmas Produced by Intense Femtosecond Laser Pulses. IEEE Jour. of Quant. Electr. 1992; 28(10): 2388–97. doi: 10.1109/3.159545

4. Streeter M, Foster PS, Cameron FH, Borghesi M, Brenner C, Carroll DC, et. al. Relativistic plasma surfaces as an efficient second harmonic generator. New J. Phys. 2011; 13: 023041. doi: 10.1088/1367-2630/13/2/023041

5. Singh M, Gupta DN, Suk H. Efficient second- and third-harmonic radiation generation from relativistic laser-plasma interactions. Phys. of Plasm. 2015; 22: 063303. doi: 10.1063/1.4922435

6. Kim S, Jin J, Kim YJ, Park IY, Kim Y, Kim SW. High-harmonic generation by resonant plasmon field enhancement. Nature. 2008; 453: 757–760. doi: 10.1038/nature07012 18528390

7. Steingrube D, Schulz E, Binhammer T, Gaarde MB, Couairon A, Morgner U, Kovacev M. High-order harmonic generation directly from a filament. New J. Phys. 2011; 13: 043022. doi: 10.1088/1367-2630/13/4/043022

8. Zheng J, Qiu E, Lin Q. High harmonic generation with sub-cycle pulses. J. Opt. 2011; 13: 075206. doi: 10.1088/2040-8978/13/7/075206

9. Lucchini M, Calegari F, Kim K, Sansone G, Nisoli M. Nonadiabatic quantum path analysis of the high-order harmonic generation in a highly ionized medium. New J. Phys. 2012; 14(2): 023025.

10. Kemper A, Moritz B, Freericks JK, Devereaux TP. Theoretical description of high-order harmonic generation in solids. New J. Phys. 2013; 15: 023003. doi: 10.1088/1367-2630/15/2/023003

11. Vampa G, Hammond TJ, Thiré N, Schmidt BE, Légaé F, McDonald CR, Brabec T, Corkum PB. Linking high harmonics from gases and solids. Nature. 2015; 522: 462–464. doi: 10.1038/nature14517 26108855

12. Neyra E, Videla F, Ciappina MF, Pérez-Hernández JA, Roso L, Lewenstein M, Torchia GA. High-order harmonic generation driven by inhomogeneous plasmonics fields spatially bounded: influence on the cut-off law. J. Opt. 2018; 20(3): 034002. doi: 10.1088/2040-8986/aaa6f7

13. Zdanowicz M, Kujala S, Husu H, Kauranen M. Effective medium multipolar tensor analysis of second-harmonic generation from metal nanoparticles. New J. Phys. 2011; 13: 023025. doi: 10.1088/1367-2630/13/2/023025

14. Pavlyukh Y, Berakdar J, Hübner W. Semi-classical approximation for second-harmonic generation in nanoparticles. New J. Phys. 2012; 14: 093044. doi: 10.1088/1367-2630/14/9/093044

15. Kautek W, Sorg N, Krüger J. Femtosecond pulse laser second harmonic generation on semiconductor electrodes. Electrochimica Acta. 1994; 39(8/9): 1245–49. doi: 10.1016/0013-4686(94)E0043-Y

16. Mlejnek M, Wright E, Moloney J, Bloembergen N. Second Harmonic Generation of Femtosecond Pulses at the Boundary of a Nonlinear Dielectric. Phys. Rev. Lett. 1999; 83(15): 2934–37. doi: 10.1103/PhysRevLett.83.2934

17. Sidick E, Knoesen A, Dienes A. Ultrashort-pulse second-harmonic generation. I. Transform-limited fundamental pulses. J. Opt. Soc. Am. B. 1995; 12(9): 1704–12. doi: 10.1364/JOSAB.12.001704

18. Kim D-W, Xiao G-Y, Ma G-B. Temporal properties of the second-harmonic generation of a short pulse. Appl. Opt. 1997; 36(27): 6788–93. doi: 10.1364/ao.36.006788 18259545

19. Telegin LS, Chirkin AS. Interaction in frequency doubling of ultrashort laser pulses. Sov. J. Quantum Electron. 1982; 12: 1358. doi: 10.1070/QE1982v012n10ABEH006061

20. Razumikhina TB, Telegin LS, Kholodnykh AI, Chirkin AS. Three-frequency interactions of high-intensity light waves in media with quadratic and cubic nonlinearities. Sov. J. Quantum Electron. 1984; 14(10): 1358–63. doi: 10.1070/QE1984v014n10ABEH006408

21. Ditmire T, Rubenchik AM, Eimerl D, Perry MD. Effects of cubic nonlinearity on frequency doubling of high-power laser pulses. J. Opt. Soc. Am. B. 1996; 13(4): 649–652. doi: 10.1364/JOSAB.13.000649

22. Choe W, Banerjee PP, Caimi FC. Second-harmonic generation in an optical medium with second- and third-order nonlinear susceptibilities. J. Opt. Soc. Am. B. 1991; 8(5): 1013–22. doi: 10.1364/JOSAB.8.001013

23. Komissarova MV, Sukhorukov AP, Trofimov VA. Self-compression of the fundamental and second harmonic pulses in the media with quadratic and cubic nonlinearities. Bulletin of the Russian Academy of Sciences. Physics supplement. Physics of Vibration. 1993; 57: 189–192.

24. Lysak TM, Trofimov VA. The bistable mode of second harmonic generation by femtosecond pulses. Technical Physics. 2001; 46: 1401–06. doi: 10.1134/1.1418503

25. Lysak TM, Trofimov VA. Bistability and uniqueness of solutions in the problem of second harmonic generation of femtosecond pulses. Computational Mathematics and Mathematical Physics. 2001; 41: 1214–26.

26. Lin R, Gao Y. Observation of the modulation instability and frequency-doubling in self-defocusing crystal. Phys. Lett. A. 2011; 375: 3228–31. doi: 10.1016/j.physleta.2011.07.007

27. Kasumova RJ, Safarova GA, Ahmadova AR, Kerimova NV. Influence of self- and cross-phase modulations on an optical frequency doubling process for metamaterials. Appl. Opt. 2018; 57(25): 7385–90. doi: 10.1364/AO.57.007385 30182960

28. Trofimov VA, Kharitonov DM, Fedotov MV. Theory of SHG in a medium with combined nonlinear response. J. Opt. Soc. Am. B. 2018; 35(12): 3069–87. doi: 10.1364/JOSAB.35.003069

29. Trofimov VA, Trofimov VV. High effective SHG of femtosecond pulse with ring profile of beam in bulk medium with cubic nonlinear response. Proceedings of SPIE. 2007; 66100: 66100R. doi: 10.1117/12.740023

30. Ashihara S, Nishina J, Shimura T, Kuroda K. Soliton compression of femtosecond pulses in quadratic media. J. Opt. Soc. Am. B. 2002; 19(10) 2505–10. doi: 10.1364/JOSAB.19.002505

31. Liu X, Qian L, Wise F. High-energy pulse compression by use of negative phase shifts produced by the cascade χ(2): χ(2) nonlinearity. Opt. Lett. 1999; 24(23): 1777–79. doi: 10.1364/ol.24.001777 18079931

32. Wyller J, Królikowski WZ, Bang O, Petersen DE, Rasmussen JJ. Modulational instability in the nonlocal χ(2)-model. Physica D. 2007; 227: 8–25. doi: 10.1016/j.physd.2007.01.002

33. Wang J, Ma Z, Li Y, Lu D, Guo Q, Hu W. Stable quadratic solitons consisting of fundamental waves and oscillatory second harmonics subject to boundary confinement. Phys. Rev. A. 2015; 91: 033801. doi: 10.1103/PhysRevA.91.033801

34. Karamzin YN, Sukhorukov AP. Nonlinear interaction of diffracted light beams in a medium with quadratic nonlinearity: mutual focusing of beams and limitation on the efficiency of optical frequency converters. JETP Lett. 1974; 20: 339.

35. Torruellas WE, Wang Z, Hagan DJ, VanStryland EW, Stegeman GI, Torner L, et. al. Observation of two-dimensional spatial solitary waves in a quadratic medium. Phys. Rev. Lett. 1995; 74: 5036. doi: 10.1103/PhysRevLett.74.5036 10058667

36. Schiek R, Baek Y, Stegeman GI. One-dimensional spatial solitary waves due to cascaded second-order nonlinearities in planar waveguides. Phys. Rev. E. 1996; 53: 1138. doi: 10.1103/PhysRevE.53.1138

37. Fuerst RA, Baboiu DM, Lawrence B, Torruellas WE, Stegeman GI, Trillo S, et.al. Spatial modulational instability and multisolitonlike generation in a quadratically nonlinear optical medium. Phys. Rev. Lett. 1997; 78: 2756. doi: 10.1103/PhysRevLett.78.2756

38. Costantini B, De Angelis C, Barthelemy A, Bourliaguet B, Kermene V. Collisions between type II two-dimensional quadratic solitons. Opt. Lett. 1998; 23: 424–426. doi: 10.1364/ol.23.000424 18084532

39. Di Trapani P, Caironi D, Valiulis G, Dubietis A, Danielius R, Piskarskas A. Observation of temporal solitons in second-harmonic generation with tilted pulses. Phys. Rev. Lett. 1998; 81: 570. doi: 10.1103/PhysRevLett.81.570

40. Liu X, Qian LJ, Wise FW. Generation of optical spatiotemporal solitons. Phys. Rev. Lett. 1999; 82: 23. doi: 10.1103/PhysRevLett.82.4631

41. Kim DH, Kang JU, Khurgin JB. Cascaded Raman self-frequency shifted soliton generation in an Er/Yb-doped fiber amplifier. Appl. Phys. Lett. 2002; 81: 2695. doi: 10.1063/1.1512823

42. Kharenko DS, Bednyakova AE, Podivilov EV, Fedoruk MP, Apolonski A, Babin SA. Cascaded generation of coherent Raman dissipative solitons. Opt. Lett. 2016; 41: 175. doi: 10.1364/OL.41.000175 26696187

43. Buryak AV, Di Trapani P, Skryabin DV, Trillo S. Optical solitons due to quadratic nonlinearities: from basic physics to futuristic applications. Phys. Rep. 2002; 370: 63. doi: 10.1016/S0370-1573(02)00196-5

44. Cheng Z, Fu HY, Li Q. Cascaded photonic crystal fibers for three stage non-integer order soliton compression. Opt. Comm. and Net. 2017; 1–2.

45. Bache M, Wise FW. Type-I cascaded quadratic soliton compression in lithium niobate: Compressing femtosecond pulses from high-power fiber lasers. Phys. Rev. A. 2010; 81: 053815. doi: 10.1103/PhysRevA.81.053815

46. Zeng X, Ashihara S, Shimura T, Kuroda K. Adiabatic femtosecond pulse compression and control by using quadratic cascading nonlinearity. Nonlinear Optics: Technologies and Applications.—International Society for Optics and Photonics. 2008; 6839: 68390B.

47. Li Q, Kutz J, Wai P. Cascaded higher-order soliton for non-adiabatic pulse compression. J. Opt. Soc. B. 2010; 27: 2180. doi: 10.1364/JOSAB.27.002180

48. Bache M, Królikowski WZ, Moses J, Wise FW. Limits to compression with cascaded quadratic soliton compressors. Opt. Expr. 2008; 16: 3273. doi: 10.1364/OE.16.003273

49. Moses J, Wise FW. Soliton compression in quadratic media: high-energy few-cycle pulses with a frequency-doubling crystal. Opt. Lett. 2006; 31: 1881. doi: 10.1364/ol.31.001881 16729102

50. Šuminas R, Tamošauskas G, Valiulis G, Dubietis A. Spatiotemporal light bullets and supercontinuum generation in β-BBO crystal with competing quadratic and cubic nonlinearities. Opt. Lett. 2016; 41: 2097. doi: 10.1364/OL.41.002097 27128083

51. Conti C, Trillo S, Di Trapani P, Kilius J, Bramati A, Minardi S, et. al. Effective lensing effects in parametric frequency conversion. J. Opt. Soc. Am. B. 2002; 19(4): 852–859. doi: 10.1364/JOSAB.19.000852

52. Di Trapani P, Bramati A, Minardi S, Chinaglia W, Trillo S, Conti C, et. al. Focusing versus defocusing nonlinearities in self-trapping due to parametric frequency conversion. Phys. Rev. Lett. 2001; 87: 183902. doi: 10.1103/PhysRevLett.87.183902

53. Trofimov V, Lysak T. Catastrophic self-focusing of axially symmetric laser beams due to cascading SHG. Proceedings of SPIE. 2011; 7822: 78220E–78220E-11. doi: 10.1117/12.891265

54. Lysak TM, Trofimov VA. Achieving high-efficiency second harmonic generation in a sequence of laser pulses with random peak intensity. Part I. Efficient generation in optical fibers. Comp. Math. and Modeling. 2008; 19: 333–342. doi: 10.1007/s10598-008-9012-z

55. Lysak TM, Trofimov VA. Achieving high-efficiency second harmonic generation in a sequence of laser pulses with random peak intensity Part II. Suppression of intensity fluctuations in a quadratic-nonlinearity medium. Comp. Math. and Modeling. 2009; 20: 1–25. doi: 10.1007/s10598-009-9015-4

56. Lysak TM, Trofimov VA. Achieving high-efficiency second harmonic generation in a sequence of laser pulses with random peak intensity Part III. Propagation of pulses in a bulk medium. Comp. Math. and Modeling. 2009; 20: 101–112. doi: 10.1007/s10598-009-9030-5

57. Trofimov VA, Lysak TM. Highly efficient SHG of a sequence of laser pulses with a random peak intensity and duration. Opt. Spectrosc. 2009; 107: 399–406. doi: 10.1134/S0030400X0909015X

58. Huttunen M, Mäkitalo J, Bautista G, Kauranen M. Multipolar second-harmonic emission with focused Gaussian beams. New J. Phys. 2012; 14: 113005.

59. O’Donnell K, Torre R. Characterization of the second-harmonic response of a silver-air interface. New J. Phys. 2005; 7: 154.

60. Rodrigo S, Laliena V, Martin-Moreno L. Second-harmonic generation from metallic arrays of rectangular holes. J. Opt. Soc. Am. B. 2015; 32: 15–25. doi: 10.1364/JOSAB.32.000015

61. Luo M, Liu Q. Extraordinary enhancement of second harmonic generation in a periodically patterned distributed Bragg reflector. J. Opt. Soc. Am. B. 2015; 32: 1193–1201. doi: 10.1364/JOSAB.32.001193

62. Yudovich S, Shwartz S. Second-harmonic generation of focused ultrashort x-ray pulses. J. Opt. Soc. Am. B. 2015; 32: 1894–1900. doi: 10.1364/JOSAB.32.001894

63. Butet J, Gallinet B, Thyagarajan K, Martin OJF. Second-harmonic generation from periodic arrays of arbitrary shape plasmonic nanostructures: a surface integral approach. J. Opt. Soc. Am. B. 2013; 30: 2970–79. doi: 10.1364/JOSAB.30.002970

64. Kolmychek IA, Krutyanskiy VL, Murzina TV, Sapozhnikov MV, Karashtin EA, Rogov VV, et. al. First and second order in magnetization effects in optical second-harmonic generation from a trilayer magnetic structure. J. Opt. Soc. Am. B. 2015; 32: 331–338. doi: 10.1364/JOSAB.32.000331

65. Samim M, Krouglov S, Barzda V. Double Stokes Mueller polarimetry of second-harmonic generation in ordered molecular structures. J. Opt. Soc. Am. B. 2015; 32: 451–461. doi: 10.1364/JOSAB.32.000451

66. Arjmand A, Abolghasem P, Han J, Helmy AS. Interface modes for monolithic nonlinear photonics. J. Opt. Soc. Am. B. 2015; 32: 577–587. doi: 10.1364/JOSAB.32.000577

67. Hardhienata H, Alejo-Molina A, Reitböck C, Prylepa A, Stifter D, Hingerl K. Bulk dipolar contribution to second-harmonic generation in zincblende. J. Opt. Soc. Am. B. 2016; 33: 195–201. doi: 10.1364/JOSAB.33.000195

68. Zhang S, Zhang X. Strong second-harmonic generation from bilayer-graphene embedded in one-dimensional photonic crystals. J. Opt. Soc. Am. B. 2016; 33: 452–460. doi: 10.1364/JOSAB.33.000452

69. Sabouri SG, Khorsandi A. Thermal dephasing compensation in high-power and high-repetition-rate second-harmonic generation using spillover loss. J. Opt. Soc. Am. B. 2016; 33: 1640–48. doi: 10.1364/JOSAB.33.001640

70. Tang D, Wang J, Zhou B, Xie G, Ma J, Yuan P, et. al. Temperature-insensitive second-harmonic generation based on noncollinear phase matching in a lithium triborate crystal. J. Opt. Soc. Am. B. 2017; 34: 1659–68. doi: 10.1364/JOSAB.34.001659

71. Guo S, Ge Y, Han Y, He J, Wang J. Investigation of optical inhomogeneity of MgO:PPLN crystals for frequency doubling of 1560 nm laser. Optics Communications. 2014; 326: 114–120. doi: 10.1016/j.optcom.2014.04.008

72. Yuan J.-H, Zhang Y, Mo H, Chen N, Zhang Z. The second-harmonic generation susceptibility in semiparabolic quantum wells with applied electric field. Optics Communications. 2015; 356: 405–410. doi: 10.1016/j.optcom.2015.08.030

73. Kanseri B, Bouillard M, Tualle-Brouri R. Efficient frequency doubling of femtosecond pulses with BIBO in an external synchronized cavity. Optics Communications. 2016; 380: 148–153. doi: 10.1016/j.optcom.2016.05.067

74. Kato K, Umemura N, Petrov V. Sellmeier and thermo-optic dispersion formulas for CdGa2S4 and their application to the nonlinear optics of Hg1−xCdxGa2S4. Optics Communications. 2017; 386: 49–52. doi: 10.1016/j.optcom.2016.10.054

75. Zhang Y, Hyodo M, Okada-Shudo Y, Zhu Y, Wang X. Characteristics of pulse width for an enhanced second harmonic generation. Optics Communications. 2017; 387: 241–244. doi: 10.1016/j.optcom.2016.11.058

76. Cai L, Wang Y, Hu H. Efficient second harmonic generation in χ(2) profile reconfigured lithium niobate thin film. Optics Communications. 2017; 387: 405–408. doi: 10.1016/j.optcom.2016.10.064

77. Leo N, Meier D, Becker P, Bohatý L, Fiebig M. Magnetically driven second-harmonic generation with phase matching in MnWO4. Optics Express. 2015; 23: 27700. doi: 10.1364/OE.23.027700 26480432

78. Zeng J, Li J, Li H, Dai Q, Tie S, Lan S. Effects of substrates on the nonlinear optical responses of two-dimensional materials. Optics Express. 2015; 23: 31817. doi: 10.1364/OE.23.031817 26698973

79. Kim S, Qi M. Broadband second-harmonic phase-matching in dispersion engineered slot waveguides. Optics Express. 2016; 24(2): 773. doi: 10.1364/OE.24.000773 26832462

80. Stegeman GI, Schiek R, Fang H, Malendevich R, Jankovic L, Torner L, et. al. Beam evolution in quadratically nonlinear one-dimensional media: LiNbO3 slab waveguides. Laser Phys. 2003; 13: 137–147.

81. Su W, Qian L, Luo H, Fu X, Zhu H, Wang T, et. al. Induced group-velocity dispersion in phase-mismatched second-harmonic generation. J. Opt. Soc. Am. B. 2006; 23(1): 51–55. doi: 10.1364/JOSAB.23.000051

82. Bache M, Bang O, Zhou BB, Moses J, Wise FW. Optical Cherenkov radiation by cascaded nonlinear interaction: an efficient source of few-cycle energetic near- to mid-IR pulses. Opt. Expr. 2011; 19: 22557–62. doi: 10.1364/OE.19.022557

83. Cai Y, Xu S, Zeng X, Zou D, Li J. High-efficiency intracavity second-harmonic enhancement for a few-cycle laser pulse train. J. Opt. 2012; 14: 105202.

84. Ota Y, Watanabe K, Iwamoto S, Arakawa Y. Measuring the second-order coherence of a nanolaser by intracavity frequency doubling. Phys. Rev. A. 2014; 89: 023824. doi: 10.1103/PhysRevA.89.023824

85. Iaconis C, Walmsley IA. Self-Referencing Spectral Interferometry for Measuring Ultrashort Optical Pulses. IEEE J. Quantum Electron. 1999; 35: 501–509. doi: 10.1109/3.753654

86. Kane DJ, Trebino R. Characterization of arbitrary femtosecond pulses using frequency-resolved optical gating. IEEE J. Quantum Electron. 1993; 29: 571–579. doi: 10.1109/3.199311

87. Baltuska A, Pshenichnikov M, Wiersma D. Second-harmonic generation frequency-resolved optical gating in the single-cycle regime. IEEE J. Quant. El. 1999; 35(4): 459–478. doi: 10.1109/3.753651

88. Bragheri F, Faccio D, Bonaretti F, Lotti A, Clerici M, Jedrkiewicz O, et. al. Complete retrieval of the field of ultrashort optical pulses using the angle-frequency spectrum. Opt. Lett. 2008; 33(24): 2952–54. doi: 10.1364/ol.33.002952 19079503

89. Hause A, Kraft S, Rohrmann P, Mitschke F. Reliable multiple-pulse reconstruction from second-harmonic-generation frequency-resolved optical gating spectrograms. J. Opt. Soc. Am. B. 2015; 32: 868–877. doi: 10.1364/JOSAB.32.000868

90. Pirogova I, Sukhorukov A. Effect of nonlinear-wave coupling dispersion on the frequency doubling of subpicosecond light-pulses. Opt. and Spect. 1985; 59(3): 694–696.

91. Tzoar N, Jain M. Self-phase modulation in long-geometry optical waveguides. Phys. Rev. A. 1981; 23(3): 1266–70. doi: 10.1103/PhysRevA.23.1266

92. Anderson D, Lisak M. Nonlinear asymmetric self-phase modulation and self-steepening of pulses in long optical waveguides. Phys. Rev. A. 1983; 27(3): 1393–98. doi: 10.1103/PhysRevA.27.1393

93. Agrawal G. Nonlinear Fiber Optics. (4th ed. Academic Press). 2007.

94. Brabec T, Krausz F. Nonlinear optical pulse propagation in the single-cycle regime. Phys. Rev. Lett. 1997; 78(17): 3282–85. doi: 10.1103/PhysRevLett.78.3282

95. Dong M-J, Tian S-F, Yan X-W, Zhang T-T. Nonlocal symmetries, conservation laws and interaction solutions for the classical Boussinesq–Burgers equation. Nonlinear Dynamics. 2019; 95(1): 273–291. doi: 10.1007/s11071-018-4563-9

96. Yan X-W, Tian S-F, Dong M-J, Wang X-B, Zhang T-T. Nonlocal symmetries, conservation laws and interaction solutions of the generalised dispersive modified Benjamin–Bona–Mahony equation. Zeitschrift für Naturforschung A. 2018; 73(5): 399–405. doi: 10.1515/zna-2017-0436

97. Wang X-B, Tian S-F, Qin C-Y, Zhang T-T. Lie symmetry analysis, conservation laws and analytical solutions of a time-fractional generalized KdV-type equation. J. Nonlinear Math. Phys. 2017; 24(4): 516–530. doi: 10.1080/14029251.2017.1375688

98. Peng W-Q, Tian S-F, Zhang T-T. Dynamics of breather waves and higher-order rogue waves in a coupled nonlinear Schrödinger equation. Europhysics Letters. 2018; 123(5): 50005.

99. Wang X-B, Tian S-F, Qin C-Y, Zhang T-T. Lie symmetry analysis, analytical solutions, and conservation laws of the generalised Whitham–Broer–Kaup–Like equations. Zeitschrift für Naturforschung A. 2017; 72(3): 269–279. doi: 10.1515/zna-2016-0389

100. Varentsova SA, Trofimov VA. Lagrangian of the process of generation of the second harmonic of femtosecond light pulses. Diff. Eq. 1998; 34(7): 997–999.

101. Varentsova SA, Trofimov VA. Invariants in the process of second harmonic generation by femtosecond light pulses. Moscow University Comput. Math. and Cybern. 1998; 4: 45–47.

102. Hovhannisyan D, Stepanyan K, Avagyan R. Computational modelling of second-harmonic generation by a femtosecond laser pulse of a few optical cycles. J. Mod. Opt. 2005; 52(1): 97–107. doi: 10.1080/09500340410001703832

103. Szarvas T, Kis Z. Numerical simulation of nonlinear second harmonic wave generation by the finite difference frequency domain method. J. Opt. Soc. Am. B. 2018; 35(4): 731–740. doi: 10.1364/JOSAB.35.000731

104. Xiao Y, Maywar DN, Agrawal GP. Propagation of few-cycle pulses in nonlinear Kerr media: harmonic generation. Opt. Lett. 2013; 38(5): 724–726. doi: 10.1364/OL.38.000724 23455278

105. Shcherbakov AA, Tishchenko AV. New fast and memory-sparing method for rigorous electromagnetic analysis of 2D periodic dielectric structures. J. Quant. Spectrosc. Radiat. Transfer. 2012; 113: 158–171. doi: 10.1016/j.jqsrt.2011.09.019

106. Weismann M, Gallagher D, Panoiu N. Nonlinear generalized source method for modeling second-harmonic generation in diffraction gratings. J. Opt. Soc. Am. B. 2015; 32: 523–533. doi: 10.1364/JOSAB.32.000523

107. Karamzin YN, Sukhorukov AP, Trofimov VA. Mathematical modeling in nonlinear optics. Moscow: Publishing house of Moscow University [in Russian]; 1989.


Článek vyšel v časopise

PLOS One


2019 Číslo 12
Nejčtenější tento týden
Nejčtenější v tomto čísle
Kurzy Podcasty Doporučená témata Časopisy
Přihlášení
Zapomenuté heslo

Zadejte e-mailovou adresu, se kterou jste vytvářel(a) účet, budou Vám na ni zaslány informace k nastavení nového hesla.

Přihlášení

Nemáte účet?  Registrujte se

#ADS_BOTTOM_SCRIPTS#