Nematic dispersive shock waves from nonlocal to local

The structure of optical dispersive shock waves in nematic liquid crystals is investigated as the power of the optical beam is varied, with six regimes identified, which complements previous work pertinent to low power beams only. It is found that the dispersive shock wave structure depends critically on the input beam power. In addition, it is known that nematic dispersive shock waves are resonant and the structure of this resonant is also critically dependent on the beam power. Whitham modulation theory is used to find solutions for the six regimes with the existence intervals for each identified. These dispersive shock wave solutions are compared with full numerical solutions of the nematic equations and excellent agreement is found.


Introduction
Nematic liquid crystals form an ideal medium to study nonlinear optics due to their "huge" nonlinearity, which is orders of magnitude larger than that of optical fibres, so that nonlinear effects can be observed over millimetre distances [1,2,3,4]. In particular, the refractive index of nematic liquid crystals increases with optical intensity, so that they form a focussing medium. When a light beam propagates through a nematic liquid crystal, the electric field of the electromagnetic wave induces dipoles in the nematic molecules, which then rotate, changing the refractive index. In addition, nematic liquid crystals have a non-local response to an optical beam in that the elastic response of the nematic extends far beyond the optical forcing [4]. An optical beam propagating through a nematic medium can then form its own waveguide, resulting in a self-guided beam, an optical solitary wave, termed a "nematicon" [2,3,5,6], which was first experimentally generated and observed in 2000 [6]. Since this first observation, nematicons, and related solitary-type waves, such as optical vortices, have become a theme of intense experimental and theoretical research effort, driven both by interest in the nonlinear optics of nematic liquid crystals and also by their potential applications in optical devices [7,8,9,10]; see Refs. [2,3,5,11] for general reviews on the nonlinear optics of nematic liquid crystals.
Solitary waves are generic wave forms for nonlinear dispersive wave equations [12], first observed and studied in the context of water waves [12,13] and fluid dynamics [12], but are widespread in nature arising, e.g., in solid mechanics [14], biology [15], ecology [16], and the above mentioned context of nonlinear optics [2,17,18], for instance. One of the appealing features of solitary waves, in addition to their widespread occurrence in nature, is that they are localized waves with steady profiles, which makes them easier to study theoretically. In addition to this, many generic nonlinear dispersive wave equations, such as the Korteweg-de Vries (KdV), nonlinear Schrödinger (NLS) and Sine-Gordon equations, are completely integrable systems via the Inverse Scattering Transform method [12,18]. Thus, a general initial condition for these equations will form a finite number of solitary waves, plus dispersive radiation. In addition, solitary wave solutions of integrable nonlinear dispersive wave equations, solitons, interact "elastically", i.e., they emerge unscathed out of the interaction without any change in their form; hence, due to this particle-like behaviour, solitary waves are termed solitons for such equations.
Another generic wave form supported by nonlinear dispersive wave equations are the dispersive shock waves (DSWs), also termed undular bores; these structures are as widespread in nature as solitary waves, with well known examples being tidal bores and tsunamis [19]. In contrast to a solitary wave, a DSW is a non-steady wave form which continuously expands. A DSW is a dispersive regularization of a discontinuity and is a modulated periodic wavetrain with solitary waves at one edge and linear, dispersive waves at the other-see [19] for a general review of DSWs. Since DSWs are non-steady waveforms, their study is more difficult than that for solitary waves. The development of DSW solutions of nonlinear dispersive wave equations relies chiefly on Whitham modulation theory [12,20,21,22], which is a version of the asymptotic method of multiple scales that is used to analyse slowly varying periodic wavetrains. Whitham modulation equations are a system of partial differential equations which govern the parameters of a slowly varying wavetrain, such as its amplitude, wavenumber, frequency and mean height. If this system is hyperbolic, then the underlying wavetrain is modulationally stable, while if it is elliptic, the wavetrain is unstable [12]. A major achievement of Whitham modulation theory was the development of the modulation equations for the KdV equation [12,21]. These modulation equations form a hyperbolic system, so that the cnoidal wave solution of the KdV equation is modulationally stable. It was subsequently realized that a simple wave solution of the KdV modulation equations is a DSW [23], even though the initial condition is a step, which is not slowly varying. This DSW solution is in excellent agreement with numerical solutions of the KdV equation [24]. The key to the determination of the simple wave DSW solution is the ability to set the modulation equations in Riemann invariant form. If the nonlinear dispersive wave equation governing the DSW is integrable, then its Whitham modulation equations can be automatically set in Riemann invariant form [25], so that the DSW solution can easily be found.
As mentioned above, the standard DSW form, termed of KdV type [19], is a modulated periodic wave with solitary waves at one edge and linear dispersive waves at the other. A non-standard DSW type is a resonant DSW [26,27], for which the waves of the DSW are in resonance with (linear) dispersive waves, resulting in a resonant wavetrain being emitted from the DSW. Resonant DSWs also occur for the KdV equation with next higher-order dispersion, i.e., fifth-order dispersion, namely for the Kawahara equation [28], and the NLS equation with next order, third-order dispersion [29,30,31,32]. If the emitted resonant wavetrain is of large enough amplitude, the KdV-type DSW structure can be destroyed; this results in the so-called traveling dispersive shock wave (TDSW) regime [26,27], consisting of a resonant wavetrain with negative polarity solitary wave, which is the remnant of the DSW, linking this to the level behind [26]-see Figure 2(d) below for an example of such a TDSW. A nematic liquid crystal is a focusing medium and, thus, optical waves are modulationally unstable; as a result, an optical DSW is not supported. However, the addition of azo dyes to the nematic medium changes its response so that it becomes defocusing [33]; in this case, nematics can support DSWs [34,35,36]. A nematic DSW is an example of a resonant DSW [34,35,36]. In these works, the nematic DSW was studied in the highly nonlocal limit, for which the nematic elastic response extends far beyond the light beam, with the nematic DSW generated by a step jump in the optical intensity. While the nematic equations are of NLStype [2], in the highly nonlocal limit the nematic bore is of KdV-type and is well described by the DSW solution of the KdV equation. The nematic DSW structure is highly dependent on the size of the jump of the optical electric field intensity generating it, with six distinct DSW types identified [36].
As stated above, in the highly nonlocal limit the nematic DSW is of KdV-type with the DSW having positive polarity. However, in the limit of weak nonlocality, the nematic equations reduce to the NLS equation [2,4], and the nematic DSW becomes the NLS DSW, which is non-resonant. The degree of nonlocality of the optical response of a nematic is inversely proportional to the power of the optical beam, with the response being highly nonlocal for lower power beams, transitioning to local as the beam power increases [2,4], as will be detailed in Section 2. In this work, the evolution of the nematic DSW structure as the degree of nonlocality ranges from highly nonlocal (low-power beams) [34,35,36] to local (high-power beams) will be studied. As the nonlocality decreases, the changes in the DSW structure from those previously found [36] in the limit of high nonlocality to the standard NLS DSW [38] will be identified and the solutions for these will be derived. It is found that there exist two additional DSW regimes over those for large nonlocality, including the NLS DSW for zero nonlocality. The new DSW regime is a transition between the KdV DSW behavior for large nonlocality and the NLS DSW behavior for very small nonlocality. In this regime, the DSW structure consists of a resonant wavetrain headed by a partial DSW which takes the solution to the initial level ahead, similar to the resonant DSW for the KdV equation with Figure 1: Sketch of nematic cell. A coherent light beam (yellow region) whose electric field E is polarised in the Y direction propagates in the Z direction through a cell filled with a dye-doped nematic liquid crystal. Thin film electrodes (black) are deposited on the upper and lower cell walls (grey). An external low frequency voltage bias V LF creates an electric field E LF to pre-tilt the molecules at an angle θ 0 to Z. The nematic molecules which are located at the boundaries are held tightly by the virtue of the anchoring films. The far right inset (black) dashed box exhibits the angular rotation of a nematic molecule with respect to the propagation direction Z in the absence (θ 0 ) and presence (θ 0 + φ) of the optical beam.
fifth-order dispersion [27,39]. As the nonlocality decreases, the optical power increases, the resonant wavetrain contracts with the leading partial DSW expanding and becoming a full NLS DSW. The analytical solutions for the various DSW types will be compared with full numerical solutions of the nematic equations.

Nematic Equations
Let us consider the propagation of a linearly, extra-ordinarily polarized, coherent light beam of wavenumber k 0 , wavelength λ 0 = 2π/k 0 , through a planar cell filled with nematic liquid crystals. The optical beam is assumed to propagate down the cell along in the Z direction, with its electric field E polarised in the Y direction. The coordinate X then completes the coordinate system. Nematic liquid crystals are a uniaxial medium consisting of elongated molecules, with the long axis termed the molecular director. The refractive index of the medium is n for optical beams polarized along the molecular director and n ⊥ for fields polarized orthogonal to the director. A fundamental property of nematic liquid crystals is that the feature the so-called Freédericksz threshold, whereby a minimum optical power is needed to rotate the nematic molecules and thus change the refractive index of the medium [1]. However, high optical powers lead to heating of the nematic medium, which can cause the loss of the nematic state if the temperature change is high enough [1,6]. One method to overcome this is to pre-tilt the nematic molecules at an angle θ 0 with respect to the Z-direction upon the application of an external static electric field E LF , so that milli-Watt power beams can rotate the nematic molecules [6]. Let us denote the optically induced rotation of the nematic by φ, so that in the presence of an optical beam the total angle of the nematic director to the Z direction is θ = θ 0 + φ. This configuration of the nematic cell is illustrated in Figure  1. The dimensional equations governing the propagation of the optical beam in the nematic cell are then of the following form: for the electric field of the beam and for the nematic response [2,3,4,5]. Here, the extraordinary refractive index of the nematic is: In the above equations, ∆ǫ = n 2 − n 2 ⊥ is the optical anisotropy, ∆ǫ LF is the low-frequency dielectric anisotropy and ǫ 0 is the electrical permittivity of free space. In addition, the constant K is the elastic medium constant in the one constant approximation for which the elastic constants of bend, twist and splay are taken equal [1,2].
The nematic equations (1) and (2) are highly nonlinear and difficult to analyse. However, for milli-Watt power beams the optical induced rotation φ is small, |φ| ≪ θ 0 , so that these equations can be expanded in Taylor series around |φ|. In addition, these equations can be put in dimensionless form using typical scales L Z down the cell and W transverse to the down cell direction, as well as a typical scale A b for the electric field of the optical beam, so that Here, (x, y, z) is the non-dimensional coordinate system and u is the non-dimensional electric field of the optical beam. The electric field scale is obtained by assuming that the input optical beam is a Gaussian beam of power P b , amplitude A b and width W b , so that Substituting these into the nematic equations (1) and (2), and expanding in Taylor series for small |φ|, we find [5,40] that suitable scalings are and the resulting non-dimensional equations read Here, the dimensionless elasticity and pre-tilting parameters, ν and q, are given by Typical experimental beam parameter values are power P b = 2mW and half-width W b = 1.5µm, with a wavelength λ 0 = 2π/k 0 = 1.064µm in the near infrared [2,5]. For the liquid crystal E7, a typical elastic constant is K = 1.2 × 10 −11 N . These parameter values give the elasticity parameter ν = O(100), as found in other studies [5,41,42]. This high value of ν means that the nematic is operating in the highly nonlocal regime, in that the elastic response of the nematic extends far beyond the waist of the optical beam [2,3,4]. However, ν is inversely proportional to the beam power P b . Note that for ν = 0 the nematic equations (7) and (8) reduce to the standard NLS equation This is the local response limit for the nematic. Note that in (2 + 1)-dimensions, beams governed by this equation are unstable and can show catastrophic collapse above a critical power [17]. It is known that a nonlocal response, ν large, stabilizes (2+ 1)-dimensional optical beams [2,3,4,5]. This is because the nematic response equation (2) is elliptic and so its solution depends on u in the entire domain, the origin of the physical concept of nonlocality. Hence, by adjusting the beam power P b , the response of the medium can be adjusted from nonlocal to local, as long as the induced heating does not destroy the nematic phase at high power.
The nematic system (7) and (8) is a focusing NLS-type system, that is the refractive index in the dimensional equation (1) increases with beam intensity |u| 2 . Since focusing NLS equations do not possess (stable) DSW solutions, the equation needs to be defocusing; in such a case, the refractive index decreases with beam intensity, and DSW solutions do exist [19]. The nematic medium can feature a defocusing response through the addition of azodyes [33]. The change in the nematic response due to the addition of the azo-dye is physically complicated, with the "order parameter" change being opposite to that without the presence of the dye. A simple model of this response change is to modify the electric field equation (7) from focusing to defocusing. In addition, the analysis of DSWs is simplest in (1 + 1)dimensions as then there are no geometric spreading effects. With these assumptions and simplifications, the nematic equations become The same system of equations also describes optical beam propagation in thermal optical media for which the refractive index depends on the temperature of the medium [43]. Such thermal optical media typically have a defocusing response. For these defocusing nematic equations, a suitable initial condition which will generate a DSW is the intensity jump initial condition The DSW solution of the defocusing nematic equations (11) and (12) has been studied in the highly nonlocal limit ν ≫ 1 in previous work [34,35,36]. As stated above, in the local limit ν ≪ 1 this system reduces to a perturbed defocusing NLS equation and becomes the defocusing NLS equation for ν = 0. The DSW solution of the defocusing NLS equation is known [38], so that the perturbed local DSW solution can be found using perturbed Whitham modulation theory [44]. Before studying the behaviour of the nematic DSW as ν varies from large to small, some previously derived results [34,35,36] will be briefly summarized. The analytical DSW solutions derived in this Section will be compared with full numerical solutions of the nematic equations (11) and (12) with the initial condition (13). The electric field equation (11) will solved using the pseudo-spectral method of Fornberg and Whitham [24], as extended [45,46] to improve the stability for high wavenumbers through the use of an integrating factor. The x derivatives are calculated using the Fast Fourier Transform (FFT) and the solution is advanced in z in Fourier space employing the fourth-order Runge-Kutta method, as detailed in previous work [36,37]. The step initial condition (13) was smoothed using the hyperbolic tangent function, as detailed in [36,37]. The director equation (12) was also solved using the FFT, with does not have a singularity at zero wavenumber due to the 2qθ term. Figure 2 displays a summary of the nematic DSW types as the nonlocality ν varies from large ν, corresponding to a highly nonlocal response for low beam power, to small ν, pertinent to a local medium response for high beam power. The terminology for the DSW regimes will be taken from previous work on radiating DSWs [26,36]. In particular, the DSW regimes from nonlocal to local response will now be detailed, for u − = 1 and u + = 0.8, with the existence intervals for the various DSW types given in Table 1.
• PDSW (perturbed DSW): This regime is illustrated in Figure 2(a). The nematic DSW is essentially a KdV DSW governed by the Kawahara equation (44) and its solution can be found as a perturbed KdV DSW [36,47]. This DSW regime is non-resonant and there is no resonant wavetrain attached to the leading edge of the DSW.
• RDSW (radiating DSW): As the nonlocality ν decreases, the DSW becomes resonant with a resonant wavetrain attached to the leading edge of the DSW, as illustrated in Figure 2(b). As all individual waves in the DSW are resonant, resonant waves are emitted from the DSW, which results in the DSW not being rank ordered [36].
• CDSW (crossover DSW): As the nonlocality decreases and the amplitude of the resonant wavetrain grows, the DSW becomes unstable with a total loss of the rank ordering of the waves of the DSW [26,36]. This regime is illustrated in Figure 2(c).
• TDSW (travelling DSW): As the amplitude of the resonant wavetrain grows, the shedding of conserved quantities into resonant radiation eventually destroys the DSW, leaving a high amplitude resonant wavetrain with a negative polarity solitary wave linking this wavetrain to the intermediate level [26] as seen in Figure 2(d). While there is a solitary wave linking the resonant wavetrain to the intermediate level, this linking can be conveniently treated as a Whitham shock [48], a shock wave in the Whitham modulation equation variables. The resonant wavetrain is brought down to the level u + ahead by a partial DSW [36]. This partial DSW has linear dispersive waves at its leading edge, but has a finite wavelength wave at its trailing edge, that is, it is not bounded by solitary waves at the trailing edge [39,49], as for a standard DSW.
• RNLS DSW (radiating NLS DSW): Further decrease in the nonlocality results in the amplitude of the linking solitary wave becoming negligible, so that the wave form consists of a (stable) resonant wavetrain headed by a partial DSW which brings the wavetrain down to the level u + ahead, as seen in Figure 2(e). This DSW regime does not occur in the high nonlocality limit as it is a "bridge" to the local NLS DSW for ν = 0.
• NLS DSW: As the nonlocality ν decreases to ν = 0 and the nematic equations (11) and (12) reduce to the NLS equation, the resonant wavetrain length contracts and the leading partial DSW evolves to a full DSW, with linear dispersive waves at the leading edge and solitary waves at the trailing edge. The resonant wavetrain then disappears and the leading DSW bringing the wavetrain down to the level u + ahead attaches to the intermediate level with |u| = u i , as in Figure 2(f). The resulting DSW is essentially an NLS DSW, which completes the transition from the KdV-type DSW for high nonlocality, that is for low power beams, to a NLS DSW for low nonlocality, that is for high power beams.
In addition to these six DSW types, when u − − u + is large enough as u + → 0 there is an additional DSW type which is a sub-case of the TDSW regime, the vacuum DSW (VDSW) Table 1: Regime boundaries for type classifications of Figure 2. Here u − = 1.0, u + = 0.8 and q = 2.
[38], for which the electric field u of the resonant radiation vanishes at a point [35,36]. As the solution for this DSW type has been previously analysed [36] and it is a sub-case of the TDSW regime, it will not be considered here. In addition, the analytical work of this manuscript is based on u − − u + being small, which is not valid in the VDSW regime. Table 1 shows that the transition regimes of TDSW and RNLS DSW which evolve the DSW from the KdV-type DSW for high nonlocality to the NLS-type DSW for weak nonlocality exist for very restricted ranges of ν, so that this transition is rapid. Over most of the range of ν the DSW is of KdV-type, one of the PDSW, RDSW and CDSW types. As noted, the first four DSW regimes also occur for the Kahawara equation and the nematic equations in the lower power, high nonlocality limit.
The standard method to analyse DSW solutions of nonlinear dispersive wave equations of the defocusing NLS-type is to transform the equation into the hydrodynamic form using the Madelung transformation [19] u where the real functions ρ and ψ denote the density and phase of the field u, while v is the fluid velocity. Upon substituting, it is found that the nematic equations (11) and (12) become The nematic equations are characterized by the linear dispersion relation [34,35] ω = kv + √ρ k νk 2 + 2q for waves around the mean levelρ for ρ andv for v. In the short wave and high nonlocality limit, νk 2 ≫ 1, this dispersion relation can be approximated by In the opposite local limit with ν small, the dispersion relation can be expanded as As expected, at the leading-order, O(1), this dispersion relation is the same as that for the NLS equation [19,50]. These dispersion relations are needed for the determination of the resonant wavetrain generated by the nematic DSW by which the linear phase velocity is matched to the velocity of the DSW front. As seen from Figure 2, the solution outside of the DSW and the resonant wavetrain is non-dispersive. On neglecting dispersion, the nematic equations (15)-(17) become the shallow water equations [12] with ρ playing the role of fluid depth and v the (horizontal) fluid velocity. In Riemann invariant form, the dispersionless nematic equations read: The initial level behind u − is linked to the intermediate shelf by a simple wave on the characteristic C − . This simple wave solution has been derived previously [34,35,36] and is of the form: We can see from the above calculations that the phase gradient on the intermediate level In the small jump limit |u − − u + | ≪ 1 the nematic equations (11) and (12) can be reduced -in the high nonlocality regime (ν ≫ 1) under consideration-to a KdV equation with fifthorder dispersion [35,51]. This will be justified below upon employing a multiscale expansion method.

Derivation of the extended KdV and Kawahara equations
We seek solutions of Eqs. (14)- (16) in the form of the following asymptotic expansions in the formal small parameter ε ≡ √ u i − u + (with 0 < ε ≪ 1) where ρ + = u 2 + , and the unknown functions ρ j , V j and φ j (j = 1, 2, 3, . . .) depend on the stretched variables Here, U will be treated as an unknown velocity, which will be determined self-consistently. Substituting the expansions (26)-(28) into Eqs. (14)- (16), and using the stretched coordinates (29), we obtain a set of equations at the different orders in ε. In particular, at the leading order, we derive the following linear equations: and where subscripts denote partial derivatives. The compatibility of the above equations suggests that the squared velocity U 2 is given by U 2 = 2u 2 + /q. Next, Eqs. (14)- (16) yield a set of nonlinear equations, namely: and The compatibility condition at this order can be found upon eliminating the fields ρ 2 , V 2 and φ 2 upon using Eqs. (32)- (33) and the definition of the velocity U . This yields the following KdV equation: To the next order of approximation, we obtain: and It is now possible to follow the procedure used at the previous order and eliminate the fields ρ 3 , V 3 and φ 3 from Eqs. (35)- (37). Indeed, solving Eq. (35) for φ 3 , Eq. (37) for ρ 3ξ and substituting into Eq. (36) eliminates every term with index 3 (recall U 2 = 2u 2 + /q). Furthermore, employing the equations obtained at the previous orders, we can express the fields φ 1,2 and V 1,2 in terms of the amplitudes ρ 1 and ρ 2 , which yields To this end, we multiply Eq. (38) by ε 2 and add it to the KdV equation Eq. (34). Then, introducing the combined amplitude function we solve for ρ 1 = P − ερ 2 and substitute the result into the above equation (38). We hence obtain the nonlinear evolution equation for the field P (ξ, η) The coefficients b j (j = 1, 2, 3, 4) appearing in Eq. (42) are given by Then, we seek an asymptotic expansion in the optical beam intensity |u| as and use the relation |u| = √ ρ. This asymptotically gives P = 2u + Q. The reductive nonlinear equation (40) can now be written in terms of the field Q(ξ, η) as The coefficients c j (j = 1, 2, 3, 4) appearing in Eq. (42) are given by Notice that Eq. (42) is the so-called extended KdV equation (eKdV), which can model the evolution of steeper waves, with shorter wavelengths, than those governed by the KdV equation. As such, the eKdV equation has been used to describe solitary waves in plasmas [52] and shallow water waves [53] in the presence of higher order effects. We note that the coefficient of the third derivative dispersive term changes sign when Hence, in the high nonlocality, low power, limit, such that ν > q 2 /(4u 2 + ), the coefficient of the third derivative in eKdV equation is positive, so that its DSW (and solitary wave) solutions have positive polarity, with solitary waves at its leading edge and linear dispersive waves at its trailing edge. On the other hand, in the local limit, ν < q 2 /(4u 2 + ), the coefficient of the third derivative is negative and the DSW has negative polarity, with linear dispersive waves at its leading edge and solitary waves at its trailing edge, so that it resembles the standard NLS DSW [38]. The nematic DSW then undergoes a change of form from a KdV-type DSW to an NLS-type DSW as ν decreases at the value of the nonlocality parameter given by (43).
At this point, it is useful to make the following remarks. First, in the highly nonlocal limit, the dominant higher-order coefficient is the one of the fifth-order dispersion term, namely c 4 ∝ ν 2 . Thus, in this limit, the eKdV equation (40) may be approximated by the Kawahara equation as was done in previous work [35,36].

Nonlocal to local nematic DSWs
The form and details of the DSW solution of the defocusing nematic equations (11) and (12) will be found as the nonlocality ν changes from O(100) to 0, that is from the nematic having a highly nonlocal response to a local response as the optical power increases. For high nonlocality, the DSW is of KdV type, with the leading edge of the DSW consisting of solitary waves of elevation [34,35,36]. The reason for this can be seen from the eKdV equation (42) as for ν > q 2 /(4u 2 + ) the coefficient of the third derivative is positive, so that the DSW is of KdV type and has positive polarity. For 0 ≤ ν < q 2 /(4u 2 + ) the sign of the third derivative is negative and the DSW has negative polarity, with solitary waves at the trailing edge and linear dispersive waves at the leading edge, as for the NLS DSW [38]. A nematic DSW of NLS type is illustrated in Figure 2(e) for ν = 0.5. Indeed, for ν = 0, the nematic equations (11) and (12) reduce to the standard NLS equation on substituting for φ from (12) into (11). The DSW solution of the NLS equation is well known [38] as the NLS equation is integrable and so the solution is completely determined. Hence, the solution for this regime will not be considered here.
The existence regions for the various nematic DSW types as the nonlocality parameter ν varies (optical power varies), found from full numerical solutions of the nematic equations (11) and (12) are shown in Figure 3 as the initial level ahead u + varies. The range 0.3 ≤ u + ≤ 0.9 was chosen as this encompassed all six of the DSW types studied here. In addition, most of the theoretical expressions for the boundaries between these regions and the solutions with each region were based on u − − u + small, for example the boundary (43), which is based on the eKdV equation (42). It can be seen that over most of the (u + , ν) domain the nematic DSW is of CDSW or TDSW type, so that it is typically unstable. In addition, it is deduced that the nematic DSW is of NLS-type only for small values of the nonlocality parameter ν of 2 and below. The nematic DSW is then nonlocal, except for high enough optical powers for which ν is small. The nonlocality parameter ν is given by (9). For the nematic liquid crystal 4-(trans-4-n-hexylcyclohexyl)-isothiocyanato-benzene (6CHBT) the parameter values are K ∼ 10 −11 N , n = 1.6335 and n ⊥ = 1.4967 [54]. Let us take the pre-tilt angle θ 0 to be π/4 so that the nematic response is maximised [2]. A typical beam wavelength is 1064nm and a typical half width W b is 1.5µm [54]. With these parameter values, it is found that ν = 2 when the beam power P b is 288mW , far in excess of typical beam powers of a few milliWatts to a few tens of milliWatts [2,54]. Such a large optical power can result in the nematic medium being heated enough so that its temperature goes above the critical temperature, 43 o C for 6CHBT [54], so that it undergoes a phase change out of the nematic state. In this regard, it should be noted that experimental nematic cells are small, of the order of 1mm in the down cell propagation direction of the beam and 100µm × 10mm in cross-section. The thin cross-section is the direction in which the pre-tilting electric field is applied, which results in a stable and uniform molecular pre-tilt. We then deduce that for experimental beam powers, the nematic bore will be in the nonlocal response regime with ν large, that is low optical power.

PDSW and RDSW regimes
Typical PDSW and RDSW solutions are illustrated in Figure 2(a) and (b). In the PDSW regime the DSW is not in resonance with linear diffractive radiation, so that the DSW is of KdV type, see the boundary (43), and is a perturbed KdV-type DSW. In the PDSW regime the DSW is in resonance with diffractive radiation, so that it consists of a KdV-type DSW with resonant radiation propagating ahead of it. This resonant radiation is not large enough, however, as to destroy the KdV-type DSW structure, as in the CDSW regime, see Figure 2(c). As the DSW in the PDSW and RDSW regimes are perturbed KdV DSWs [36], the solutions in these two regimes can be found using the perturbed KdV DSW solution of [47]. In this previous work, the general eKdV equation, a particular case of which is (42), was asymptotically transformed to the KdV equation, whose known DSW solution [19,23] was then used to find the asymptotic DSW solution of the original eKdV equation. This asymptotic DSW solution can be used here based on the eKdV equation reduction (42) of the full nematic equations in the limit u − − u + small. The work [47] then gives the following PDSW and RDSW nematic DSW solutions. The amplitude of the DSW is its wavenumber is and its mean level is Here, K(m) and E(m) are complete elliptic integrals of the first and second kinds of modulus squared m, respectively. The modulus squared m is a parameter in these amplitude, wavelength and mean height expressions. It is determined in terms of the simple wave (similarity) variable x/z by This expression for x/z derives from the characteristic of the KdV modulation equations on which the simple wave DSW solution occurs [23,47,24]. The coefficients C j (j = 1, 2, 3, 4) in the above solutions are connected to c j (j = 1, 2, 3, 4) through the relations At the leading, solitary wave edge of the DSW m → 1 and at the trailing, harmonic wave edge of the DSW m → 0. It can then be found from the characteristics (48) that the DSW lies in the range where s i and s + are the harmonic and solitary wave edge velocities of the DSW, respectively.
Comparisons between the lead solitary wave amplitude a + as given by the asymptotic DSW solution (45), with m = 1, and numerical solutions are given in Figure 4 as the nonlocality parameter ν varies. The existence regions for the PDSW and RDSW types depend on both u + and ν, see Figure 3, so that the comparison curves for each u + were stopped at the boundary between the RDSW and CDSW regimes. Figure 4 shows the lead solitary wave amplitude for the full eKdV equation (42) and the Kawahara equation, which is (42) with c 1 = c 2 = c 3 = 0. Previous work on the nematic bore [35,36] was based on the Kawahara equation, that is only the higher order fifth derivative was included in the asymptotic eKdV equation, so the lead wave amplitude based on this equation is given in the Figure to determine the effect of the extra higher order terms in the full eKdV equation (42). A key observation is that the height of the lead wave of the DSW depends very weakly on the strength of the nonlocality, with little variation even down to ν = O(10) from the high nonlocality amplitudes with ν = O(100). In some sense, the nematic DSW is then nonlocal down to small values of ν in the PDSW and RDSW regimes, which was also deduced above from Figure 3. It can also be seen that the inclusion of the extra higher order terms in the eKdV equation (42) over the Kawahara equation improves the agreement with numerical solutions on the whole, especially as the nonlocality parameter ν decreases, but the effect of these extra terms is small, with the Kawahara equation giving good agreement over the whole range of ν and for all values of u + , except near the RDSW/CDSW borderline at u + = 0.7. This is expected as the weak dependence on the nonlocality parameter ν means that the DSW is nonlocal, so that ν can be taken as large. The dominant higher order term in the eKdV equation is ε 2 c 4 Q ξξξξξ in this limit, as noted in (44). As u + = 0.7 is approached the resonant radiation shed by the DSW is of relatively large amplitude as in this limit the RDSW/CDSW boundary is approached. This results in oscillations in the lead wave amplitude as the resonant radiation moves through the DSW and is shed. In these cases, the numerical amplitude shown in Figure 4 was calculated as an average in z over the last few amplitude oscillations in the numerical solution. It is noted that except for u + = 0.75 and 0.7 the amplitude grows as ν decreases. This change in behaviour is due to the DSW changing form as it transitions from the RDSW to the CDSW regime, for which the resonant radiation has a major effect on the DSW with its amplitude decreasing markedly due to the large amount of mass being shed as resonant radiation [36].
The full linear dispersion relation (18) for the nematic equations (11) and (12) is when the mean 2ρ/q is added. This mean term arises on integrating v = ψ x to obtain the dispersion relation for ψ from that for v [35]. In the limit νk 2 large, this dispersion relation becomes (19) again with the carrier wave phase shift 2ρ/q added to obtain the dispersion relation for ψ. In the RDSW regime the resonant wavetrain has low amplitude, see Figure 2(b). In addition, this wavetrain has high frequency relative to the DSW, again see Figure 2 giving so that on settingρ = √ u + andv + = 0 as the resonant wavetrain propagates on the level ahead.
This gives the wavenumber of the resonant wavetrain based on this criterion. The resonant wavetrain then exists if s + ≥ 2u + / √ q, so that the borderline between the PDSW and RDSW regimes is s + = 2u + / √ q. Previous work [35,36] has shown that this theoretical borderline is in excellent agreement with numerical solutions in the high nonlocality limit ν large. For fixed u + , as ν decreases a PDSW changes to an RDSW, then to a CDSW, see Figure 3. For instance, for u + = 0.8, the DSW changes from PDSW to RDSW at ν = 88, then to CDSW at ν = 34. However, the resonance condition (53), or (54), gives that the DSW changes from PDSW to RDSW at ν = 3, which is the TDSW regime according to Table 1. This resonance condition is based on the limit νk 2 ≫ 1, but even if the full dispersion relation (51) is used for the resonance condition (53) the predicted PDSW/RDSW borderline is ν = 2.38, which is still far from the numerical value and close to that for νk 2 ≫ 1. The resonance condition (53) is based on resonance between the lead wave of the DSW and diffractive radiation. However, as pointed out previously [36], a DSW is modulated periodic wave so that all waves of the DSW can resonant with diffractive radiation, not just the lead wave, as seen in Figure 7 for a PDSW. Internal resonance will be discussed in detail in Section 5. The phase velocity of a component wave of the DSW is (72). Equating this bore component phase velocity with the nematic diffractive radiation phase velocity determined from the dispersion relation (18) determines the internal resonance. However, even using this internal resonance does not result in a borderline between the PDSW and RDSW regimes in any reasonable accord with numerical solutions, see Figure 3. A resonant wavetrain will then exist if the internally resonant waves can propagate out of the DSW, that is their group velocity is greater than the velocity of the lead solitary wave of the DSW. However, even this condition does not give the correct boundary between the PDSW and RDSW regimes as the nonlocality parameter ν decreases. The issue of internal resonance and its relation to the existence of the PDSW and RDSW regimes merits further study. In this regard, the recent work [55] on the interaction of linear wavepackets and DSWs could be relevant.

CDSW regime
For fixed a nonlocality parameter ν, as u + decreases, the jump height u − − u + increases, the nematic DSW changes from RDSW to CDSW form, see Figures 3(b) and (c). The reason for this is that as the amplitude of the resonant wavetrain grows, it takes more conserved quantities from the DSW, so that its amplitude is reduced. In addition, the DSW becomes unstable in the CDSW regime, as for the Kahawara equation DSW [26], noting that in the small jump height limit the nematic equations reduce to the eKdV equation (42), which becomes the Kahawara equation (44) in the limit of large nonlocality ν. Figure 2(c) shows a typical nematic DSW in the CDSW regime. It can be seen that the DSW has changed from a modulated wavetrain with a monotonically decreasing amplitude from front to rear to a disordered wavetrain with an essentially uniform amplitude on average, except at its rear. This structure is in agreement with unstable DSW structure for the focussing NLS equation [56]. The unstable DSW can then be approximated by a train of equal amplitude solitary waves, which has been found to give good results for DSW solutions [57], particularly unstable DSWs, and for the particular case of the nematic CDSW in the high nonlocality, low optical power, regime [36]. The amplitude of the solitary waves of the CDSW is determined from mass and energy conservation equations for the underlying nonlinear dispersive wave equation [57]. If we set for simplicity, the eKdV equation (42) has the mass conservation equation The derivation of the energy conservation for the eKdV equation (42) is not as straightforward. Multiplying the eKdV equation by Q and integrating by parts gives ∂ ∂η The final term on the right had side of this equation cannot be expressed as a perfect derivative. However, it can be approximately expressed in this form on noting that ε is small, so that at first order the eKdV equation (42) is the KdV equation We then have at leading order that This expression may now be used to eliminate the term in QQ ξξξ in equation (57) to give the final energy conservation equation which is accurate to O(ε 2 ).
To obtain an approximation to the nematic CDSW, let us assume that at position η the CDSW consists of N equal solitary waves of amplitudeã s and widthw s [57], where we shall use tildes to denote scaled variables in the moving and stretched coordinates (ξ, η). It is also assumed that the CDSW sheds a uniform downstream resonant wavetrain of (scaled) amplitudeã r . Then as ξ → −∞, Q → 1 and as ξ → ∞, Q →ã r cos k r ξ −ω r η , since |u| = u + + ε 2 Q with ε 2 = u i − u + . As the CDSW is approximated by a train of solitary waves, the solitary wave solution of the eKdV equation (42) is also needed. While there is no known exact solitary wave solution of this equation, there is an asymptotic solution for ε ≪ 1 [58]. To use this solution, the eKdV equation (42) needs to be rescaled to conform with the eKdV scaling of [58]. Performing this, we find that the asymptotic solitary wave solution of the eKdV equation (42) is where The rescaled coefficients c i , i = 1, . . . , 4, of the eKdV equation (42), denoted by C i , i = 1, . . . , 4, are given by The final quantity to determine is the amplitude of a r of the resonant radiation, which is related to the scaled amplitudeã r by a r = ε 2ã r = (u i − u + )ã r . This resonant radiation was determined as a WKB solution of the nematic equations (11) and (12) by linearising about the mean level u + of the resonant radiation [35]. This WKB solution gives the amplitude of the resonant radiation as Here, s + is the unscaled velocity of the CDSW, which is [36] s + = 2 The resonant radiation wavenumber k r is determined by the resonance condition (54) and the group velocity of the resonant radiation is given by the k derivative of the short wave dispersion relation (52). The resonant radiation is a solution of the nematic equations in the limit νk 2 ≫ 1, so that the appropriate group velocity for it is that from the dispersion relation (52), not the linearised KdV group velocity of the eKdV equation (42) [36]. Figure 5 displays comparisons of the nematic DSW amplitude a s and resonant wave amplitude a r in the CDSW regime as given by (68) and (70) and numerical solutions. A typical CDSW is shown in Figure 2(c), with details of the actual CDSW of this figure shown in Figure 6(a). It can be seen that the lead waves of the DSW have an approximately uniform amplitude, with a rapid decrease of the amplitude towards the trailing edge of the CDSW, as also illustrated in Figure 6(b). This solution, and that of Figure 6(a), are typical structures for an unstable DSW [56]. The numerical DSW amplitude for the comparisons of this figure was calculated as an average over the approximately uniform waves at the leading edge, which is the same assumption on which the equal amplitude approximation used to calculate the solitary wave amplitude a s was based. Figure 5(a) shows comparisons for the amplitude a r of the resonant wavetrain leading the CDSW, see Figure 2(c). It can be seen that there is excellent agreement between the theoretical amplitude with the numerical amplitude for the larger values of the level ahead u + , with the increase of a r as the nonlocality parameter ν decreases being correctly given. This agreement is much improved through the inclusion of the all the higher order terms in the eKdV equation (42) than that of previous work [36] based on the Kawahara equation, for which c 1 = c 2 = c 3 = 0, as given by the green dashed line in the figure. As u + decreases and the TDSW regime is approached, the agreement between theory and numerical solutions decreases. This is shown particularly in the final comparison of Figure 5(a) for u + = 0.5, which is near the TDSW boundary, see Figure 3. The reason for this decreasing agreement is that as the TDSW regime is approached the number of waves in the CDSW decreases so that only one lead wave is left, see Figure 6(c). The approximation that an average can be taken over an equal amplitude wavetrain then breaks down. Figure 5(b) displays similar comparisons for the amplitude a s of the nematic CDSW with numerical solutions. It can be seen that the comparison for the DSW amplitude is similar to that for the resonant wave amplitude. It should be noted that different values of u + have been chosen for the DSW amplitude comparisons for the sake of clarity. The inclusion of all the higher order terms in the eKdV equation (42) results in a significant improvement in the agreement with numerical solutions over that based on the Kawahara equation with c 1 = c 2 = c 3 = 0 when the level ahead u + is away from the TDSW/CDSW boundary of Figure 3. As u + approaches the TDSW/CDSW boundary, the DSW amplitude as given by the eKdV equation differs significantly from the numerical amplitude. The reason for this is that discussed in the previous paragraph for the resonant wave amplitude, the fact that the CDSW ceases to be a train of equal amplitude solitary waves, but reduces to a few waves.
The next thing that we discuss is the analytical borderline between the nematic CDSW and RDSW regimes. This borderline can be found by determining when the resonant amplitude (70) approaches zero or approaches a minimum as a function of the nonlocality parameter ν. The determination of this borderline is similar to that of previous work [36] in which the resonant amplitude was found as a function of the initial state u + , rather than ν. This approach gave a borderline for u + ≤ 0.73 for large ν, with the CDSW regime not existing above this value of u + for large ν. Above this limit, the resonance condition (53) ceases to work as a function of ν and an alternative method needs to be found to determine the borderline. To determine the borderline in this case we exploit the fact that the (interior) structure of the nematic DSW is resonant, as evidenced in Figure 7, not only the leading, solitary wave edge. A DSW is an extended modulated periodic wavetrain, so that all its component waves can resonant, not just the leading edge. To verify this internal resonance,  (11) and (12) in PDSW regime at z = 1500 for initial condition (13)  the phase velocity of the modulated cnoidal waves forming the DSW needs to be matched with the nematic linear phase velocity on the local mean level of the DSW. By way of illustration, we equate the DSW phase velocity [47] where C 1 , C 3 and C 4 are given in (49), with the nematic phase velocity from the full linear dispersion relation (51) on the DSW background (47) and solved for internal resonant wavenumbers, which are always positive and real, as the modulus m varies from near zero to near one. As the DSW parameter expressions (47) and (72) are only valid for a well-ranked DSW (stable DSW), such as the PDSW and the RDSW, and the CDSW is an ill-ranked DSW (unstable DSW), then a borderline exists when these wave parameter expressions result in imaginary internal resonant wavenumbers, which are unphysical, at a borderline value of ν. The mean flowv in the dispersion relation (51) was determined from the extended KdV reduction of the nematic equations of Section 4. Substituting the O(ε 2 ) mean flow v 1 given by (31) into the mean flow perturbation expansion (27) for v and then averaging gives This then completes the determination of the CDSW/RDSW borderline. A comparison between the theoretical and numerical borderlines between the RDSW and CDSW regimes is shown in Figure 3. It can be seen that the theoretical borderline is in excellent agreement with the numerical borderline for large values of the nonlocality parameter ν down to around ν = 50, with poorer agreement towards the local limit with ν small for which the DSW is changing form from KdV type to NLS type. Figure 2(d) displays a typical DSW in the TDSW regime. There is (almost) a constant amplitude resonant wavetrain which at its trailing edge is connected to the intermediate level u i . At its leading edge there is a modulated wavetrain which takes u down to the level u + ahead. This wavetrain leading the resonant wave is a partial DSW [27,36,49]. A partial DSW differs from a standard DSW in that it connects a uniform state to a periodic wavetrain, unlike a standard DSW which links two different levels. While there is a negative polarity solitary wave connecting the resonant wavetrain to the intermediate level [26], this connection can be approximated by a Whitham shock [36,48], a shock, a jump, in the modulation parameters of the wavetrain, wavelength, frequency, amplitude and mean level, of the Whitham modulation equations for the modulated periodic wavetrain [12]. A Whitham shock is determined from the Whitham modulation equations. As noted above, the nematic DSW is in the nonlocal regime, so the appropriate Whitham modulation equations are those for ν large. The Whitham modulation equations in the highly nonlocal limit ν ≫ 1 (low optical power) have previously been determined [36], and so these modulation equations will just be quoted here. These modulation equations determine the mean levelρ of ρ, the amplitude a and the wavenumber k of the Stokes' wave solution of the nematic equations (15)- (17). As there is no known general periodic wave solution of the nematic equations, the highly nonlocal Whitham modulation equations are derived based on the weakly nonlinear Stokes wavetrain for the nematic equations [36]. It can be seen from Figure 2(d) that the resonant wavetrain has small amplitude, so the weakly nonlinear limit is appropriate. In the highly nonlocal limit ν ≫ 1, the Stokes' wave solution of the nematic equations (15)- (17) is
where the uniform phase is ϕ = kx − ωz and the over bar,ρ andv, denotes the mean value of a wave parameter [36]. The amplitude a of the (Stokes) wave is assumed to be small. The work of [36] gives that at O(a) the nematic Stokes coefficients and the O(a 2 ) correction to the dispersion relation are The coefficient ω 1 is set to zero, ω 1 = 0, to eliminate secular terms. The weakly nonlinear Whitham modulation equations can then be derived by averaging conservation laws deduced from Nöther's Theorem [59]. The Lagrangian for the hydrodynamic form of the nematic equations (15)- (17) is Applying Nöther's Theorem we have that translation invariance with respect to the phase ψ gives the mass conservation law Substituting the Stokes expansions (74)-(76) into these conservation laws and averaging by integrating in ϕ from 0 to 2π [12] yields the modulation equations, truncated to O(1/ν), for the (slowly varying) amplitude a, wavenumber k and meansρ andv of the modulated Stokes wave [12]. The modulation equation (87) is that for optical power conservation, equation (88) is momentum conservation and (89) is energy conservation. The modulation equation (86) is the equation for conservation of waves, k z + ω x = 0, on noting that the x derivative term is just the dispersion relation (77). In this regard, it should be noted that the dispersion relation for the Stokes' wave from which the modulation equations are calculated has the carrier waves' phase shift term 2ρ/q added [35,36], as explained above,  [26]. As discussed above, this link between the resonant wavetrain and the intermediate level can be approximated by a Whitham shock, a jump in the modulation equation variables [36,48], so that the Whitham shock links the resonant wavetrain with the level u i behind, in a similar manner to a gas dynamic shock wave links two compressible flow states [12]. Let us denote the amplitude, wavenumber, mean level and mean phase gradient of the resonant wavetrain by a r , k r ,ρ r and v r , respectively. Matching the Whitham shock velocity U s to the Stokes' wave velocity (90), as these are co-propagating, gives Ahead of the Whitham shock there is the resonant wavetrain and behind it is a flat shelf, the intermediate level, which is a wavetrain of zero amplitude. The mass, momentum and energy conservation equations (87)-(89) then give the jump conditions across the Whitham shock. Together with the resonance condition (91), these form four equations for the five unknowns U s , a r , k r ,ρ r andv r , noting that ρ i and v i are given by (24) and (25), respectively. The final equation is obtained by assuming that the Riemann invariant R − (22) is conserved through the resonant wavetrain and its lead partial DSW, which is valid for a full DSW [19,50]. This then determines the mean of the resonant phase gradientv r [36]. This Riemann invariant condition gives The above nematic Whitham modulation equation jump conditions (92)-(94) with (95) can be solved numerically for U s , a r , k r andρ r using Newton's method. The full details for this numerical solution of the Whitham shock jump conditions can be found in [36]. Figure  8 shows comparisons for the Whitham shock velocity from the nonlocal to local limits, the optical power increasing, as given by the jump conditions and full numerical solutions of the nematic equations (11) and (12). The values of the level ahead u + were chosen to lie in the TDSW regime, see Figure 3. It can be seen that there is excellent agreement between the theoretical and numerical solutions from high nonlocality, ν large, down to ν = O(10). As for the lead solitary wave amplitude comparison of Figure 4 there is little change in the Whitham shock velocity as the nonlocality parameter ν varies by an order of magnitude, with only a small increase in the velocity. There is a small, increasing deviation between the theoretical and numerical values towards ν = 20. This is due to the onset on the VDSW regime for which u vanishes at a point, a vacuum point [35,36]. Once the vacuum point is reached, |u| cannot decrease further, so that the Whitham shock jump conditions need to be modified [36]. This will not be pursued further here.
For a fixed level ahead u + , as the nonlocality parameter ν decreases (optical power increases), the nematic DSW evolves from CDSW to TDSW type, see Figure 3. The borderline between the CDSW and TDSW regimes can be determined from the Whitham shock jump conditions (92)-(94) and the resonance condition (91) based on the following condition. For a fixed nonlocality parameter ν, as the level ahead u + increases, it is found that the Whitham shock velocity becomes greater than the linear group velocity of the resonant wavetrain. This is unphysical as this would mean that the resonant wavetrain could not form. Figure 3 shows this theoretical bound between the CDSW and TDSW regimes as a red line. It can be seen that the agreement with numerical solutions is excellent over the entire range of ν, even for jump heights u − − u + which are not small. As ν decreases, it would be expected that the high nonlocality modulation equations (86)-(89) cease to be applicable. The Whitham modulation equations for the nematic equations in the local limit, ν small, were calculated based on the equivalent of the Stokes' wave expansions (74)-(77) and (78)-(81), expanding in ν rather than 1/ν. However, these were found not to give solutions in good agreement with numerical solutions. This is expected as Figures 3 and  4 show that the nematic DSW is highly nonlocal down to very small values of ν, which are unphysical, as discussed at the end of Section 3.

RNLS and NLS DSW Regimes
The KdV approximation (44) gives that the DSW changes from KdV to NLS type at ν = q 2 /(4u 2 + ) due to the change in sign of the coefficient of the third derivative, as discussed above, and as shown in Figure 3. It should be noted that numerical solutions do not show a distinct change in DSW type, but a transition from a KdV-type DSW to an NLS-type DSW, as seen in Figures 2(d) and (e) and as shown by the two regimes RNLS and NLS type in Figure 3. The TDSW regime is characterised by a negative polarity solitary wave connecting the resonant wavetrain to the intermediate level, as in Figure 2(d). As the nonlocality ν decreases, the beam power in increases, the height of this solitary wave decreases, resulting in the RNLS regime which consists of a Whitham shock connecting a resonant wavetrain to the intermediate level u i . Ahead of the resonant wavetrain is a partial NLS-type DSW which connects to the level ahead u + , see Figure 2(e). As ν decreases further, the resonant wavetrain shrinks and the partial DSW becomes a full NLS DSW with solitary waves at its trailing edge and linear waves at its leading edge, resulting in an NLS DSW for sufficiently small ν, see Figure 3(b). In the limit the NLS DSW alone links the intermediate level to the level ahead. As ν decreases in the RNLS regime, the waves at the trailing edge of the partial DSW evolve from weakly nonlinear Stokes waves to fully nonlinear periodic waves, then solitary waves in the NLS DSW regime.
The intermediate RNLS state, illustrated in Figure 2(e), consists of a resonant wavetrain with the height of the negative polarity solitary wave at the Whitham shock negligible. To compare the KdV-NLS DSW boundary ν = q 2 /(4u 2 + ) with numerical solutions, the choice of the height of this solitary wave being 5 × 10 −3 above u i was chosen for the onset of the RNLS regime in numerical solutions. It can be seen that there is good agreement for this regime boundary for u + close to u − , but there is increasing disagreement as u + decreases. This is expected as the reductive nematic eKdV equation (42) was derived under the assumption that |u − − u + | is small.

Conclusions
The structure of the nematic DSW (dispersive shock wave) solution of the defocussing nematic equations governing the propagation of an optical beam through a cell filled with nematic liquid crystals has been investigated using a combination of numerical solutions of the equations governing the beam, consisting of an NLS-type equation for the optical beam and an elliptic medium response equation, and solutions of the governing nematic equations using Whitham modulation theory and/or asymptotic solutions. In contrast to previous work [34,35,36] the evolution of the DSW structure was studied as the power of the optical beam varied, from the experimental low power for which the nematic response is nonlocal to high power for which the response is local. As the beam power varies, it was found that the nematic DSW transitions between six regimes, four of which were studied in previous work [34,35,36]. The two NLS-type DSWs do not exist in the low power regime studied in this previous work. However, the experimental verification of these high power DSW types is questionable as the powers required for their existence are unrealistic due to the possible excessive medium heating the high beam powers would cause. Excellent agreement was found between numerical solutions and analytical solutions for the four physically relevant DSW regimes, the PDSW, RDSW, CDSW and TDSW regimes displayed in Figure 2(a)-(d). In particular, the analytical theory gives good agreement for the boundaries between the existence regions for five of the DSW types as the nonlocality parameter ν varies, the exception being the boundary between the PDSW and RDSW regimes.
It has been found that the details of the nematic DSW, for instance its lead wave amplitude and velocity and the amplitude of the associated resonant radiation, are well approximated by the nematic equations in the high nonlocality limit, the nonlocality parameter ν large, as the DSW transitions to the local limit. This holds for the PDSW, RDSW, CDSW and TDSW regimes with the nonlocality parameter ν ranging from O(100) to O(1). As the analysis of the nematic DSW is much easier in the high nonlocality limit based on asymptotic analyses with ν ≫ 1, this is an important result for future analysis of the nematic DSW in its various regimes and over its nonlocality range.
There are still a number of issues which could be addressed by the future work. An outstanding issue is the correct determination of resonance between the nematic DSW and diffractive radiation. In contrast to the high nonlocality limit with ν large [34,35,36], as ν decreases from the high to the low nonlocality limit, the beam power increases, the resonance condition used in previous work in the RDSW and CDSW regimes that the velocity of the lead solitary wave of the DSW matches the linear phase velocity of the resonant waves does not agree with numerical solutions. The theoretical transition from the PDSW regime for high nonlocality to the RDSW regime due to the onset of resonance as ν decreases occurs for ν = O(1), while numerical solutions show the transition for ν = O(100) to O(10). As noted in previous work [36], and in contrast with other work on resonant DSWs, not only is the lead wave of the DSW in resonance with diffractive radiation, but the whole modulated periodic wave which forms the DSW is in resonance. However, even this observation does not yield the correct resonance condition for the RDSW regime as the nonlocality parameter ν decreases below the highly nonlocal limit. The correct resonance condition between the DSW and diffractive radiation requires further study. This should be an important general issue for all resonant DSWs beyond the specific application to nematic liquid crystals.