A $\phi^6$ soliton with a long-range tail

TL;DR

This paper introduces an analytically solvable sextic potential model with non-parity symmetric solitons exhibiting unique long-range tail behavior, and analyzes their spectral properties with potential for further theoretical exploration.

Contribution

The paper presents a new solvable $^6$ model with explicit Lambert W function solutions, revealing long-range tails and spectral characteristics, expanding understanding of non-parity symmetric solitons.

Findings

Solvable model with Lambert W function solutions

Solitons with power-law and exponential asymptotics

Analytical spectrum of boson and fermion fields

Abstract

We propose an analytically solvable sextic potential model with non-trivial soliton solutions connecting the trivial vacua. The model does not respect parity symmetry, and like $\phi^4$ theory has two minima. The soliton solutions and the consequent results are obtained in terms of the Lambert W function, i.e., the inverse function of $f(W) = We^W$. They have power-law asymptotics at one spatial infinity and exponential asymptotics at the other. We compare the solution with the kink of $\phi^4$ theory, which preserves the parity symmetry and has exponential asymptotics at both spatial infinities. Moreover, we study the full spectrum (bound and continuum states) of boson and fermion fields in the presence of the proposed soliton. We consider two types of coupling for the boson-soliton interaction and Yukawa coupling for the fermion-soliton interaction. Most results are derived analytically. This property renders the model a fertile ground for further study, including parity breaking related phenomena and long-range soliton-soliton interactions.

Equations

$$\mathcal{L}=\frac{1}{2}{\partial }_{\mu }\varphi {\partial }^{\mu }\varphi -V(\varphi ),$$

$$V(\varphi )=\frac{{\lambda }^2}{2}{\left(1-{\varphi }^2\right)}^2{\left(1-\varphi \right)}^2.$$

$${\partial }_x\varphi -\left(1-{\varphi }^2\right)\left(1-\varphi \right)=0,$$

$$\frac{1}{4}\left[\log \left(\frac{\varphi +1}{\varphi -1}\right)-\frac{2}{\varphi -1}\right]=x+C,$$

$${\varphi }_s(x)=1-\frac{2}{1+W\left[e^{1+4x}\right]}.$$

$$\{\begin{array}{ll}{\varphi }_s(x)\to -1+2e^{1+4x}& x\to -\infty ,\\ {\varphi }_s(x)\to 1-1/(2x)& x\to \infty .\end{array}$$

$$M_{cl}$$

$$={\int }_{-\infty }^{\infty }\left[\frac{1}{2}{\left(\frac{\mathrm{d}\varphi }{\mathrm{d}x}\right)}^2+V(x)\right]\mathrm{d}x={\int }_{-1}^1(1-{\varphi }^2{)}^2(1-\varphi {)}^2\mathrm{d}\varphi =\frac{4}{3}(\lambda ),$$

$$\displaystyle\left[-\partial_{x}^{2}+U[\phi(x)]\Bigr{|}_{\phi_{s}(x)}\right]\eta_{n}(x)=\omega_{n}^{2}\eta_{n}(x)\,,$$

$$\displaystyle U[\phi(x)]\Bigr{|}_{\phi_{s}(x)}$$

$$=\frac{16\left(1-8W\left[e^{1+4x}\right]+6{\left(W\left[e^{1+4x}\right]\right)}^2\right)}{{\left(1+W\left[e^{1+4x}\right]\right)}^4},$$

$${\eta }_0={\partial }_x\varphi =\frac{8W\left[e^{1+4x}\right]}{{\left(W\left[e^{1+4x}\right]+1\right)}^3}.$$

$${\frac{d^{2}V}{d{\varphi }^2}|}_{\varphi =1}=0,{\frac{d^{2}V}{d{\varphi }^2}|}_{\varphi =-1}=16({\lambda }^2),$$

$$\mathcal{L}=\frac{1}{2}{\partial }_{\mu }\varphi {\partial }^{\mu }\varphi -V(\varphi )+\frac{1}{2}{\partial }_{\mu }\chi {\partial }^{\mu }\chi +\frac{1}{2}{m}^2{\chi }^2-g\varphi \chi ,$$

$$\left(\square -m^2\right)\chi =-g\varphi .$$

$$\left({\partial }_x^2+k^2\right){\chi }_s=g\varphi ,$$

$${\chi }_s(x)$$

$$-\frac{g}{\kappa }\int \sinh \left(\kappa u\right)\left(1-\frac{2}{W\left[e^{1+4(x-u)}\right]+1}\right)\text{d}u=-\frac{g}{{\kappa }^2}+\frac{2g}{\kappa }\int \frac{\sinh \left(\kappa u\right)}{1+W\left[e^{1+4(x-u)}\right]}\text{d}u.$$

$$-\frac{g}{{\kappa }^2}-\frac{g}{2\kappa }\int \sinh \left[-\frac{\kappa }{4}\left(\ln (w)+w-1-4x\right)\right]\frac{1}{w}\text{d}w$$

$$=$$

$$-$$

$${\chi }_s(x)$$

$$-$$

$$B=\frac{g}{16n}{\left(-n\right)}^{-n}{e}^{-n}\Gamma \left(n\right).$$

$${\chi }_s(x)=A{e}^{ikx}+B{e}^{-ikx}+\frac{g}{k}{\int }^x\sin [k(x-y)]\varphi (y)\text{d}y,$$

$$-\frac{g}{2k}\int \sin \left[-\frac{k}{4}\left(\ln (w)+w-1-4x\right)\right]\frac{1}{w}\text{d}w$$

$$=$$

$${\chi }_s(x)$$

$$+\frac{g}{2k}\mathrm{Im}\left[\Gamma (-ik/4,ikW[e^{1+4x}]/4)e^{i[1+4x+\ln (k/4)+i\pi /2]k/4}\right].$$

$${\chi }_s(x\to +\infty )$$

$${\chi }_s(x\to -\infty )$$

$$+\frac{g}{2k}\text{Im}\left[{\left(\frac{k}{4}\right)}^{ik/4}{e}^{(\pi +i)k/4}\Gamma \left(-i\frac{k}{4}\right)e^{ikx}\right].$$

$$\mathcal{L}=\frac{1}{2}{\partial }_{\mu }\varphi {\partial }^{\mu }\varphi -V(\varphi )+\frac{1}{2}{\partial }_{\mu }\chi {\partial }^{\mu }\chi +\frac{1}{2}{m}^2{\chi }^2+g\varphi {\chi }^2,$$

$$\left(\square -m^2\right)\chi -2g\varphi \chi =0.$$

Full text

1 Introduction

There is a small class of nonlinear differential equations with soliton solutions. A soliton is a stable solution with localized energy density. It arises as a consequence of the interaction between nonlinearity and dispersion when a nonlinear sharpening term counterbalances the dispersive term. The competition between these two contributions shapes the structure of the soliton and provides its stability. In the language of topology, the soliton configuration has an associated conserved topological charge or winding number, which protects it against decay into a trivial configuration. Solitons are fascinating due to their mathematical properties. However, their usefulness extends far beyond that, touching multiple areas of science. In particular, they are subject of research in diverse areas of physics, including high energy physics, nonlinear optics and condensed matter physics [1, 2, 3, 4, 5]. Amongst the most known solitons are skyrmions and domain walls in magnetic materials [6, 7, 8], vortices in superconductors and fluids [9, 10, 11, 12] as well as magnetic monopoles, Q-balls, cosmic strings and instantons in high-energy physics [13, 14, 15, 16, 17, 18, 19, 20, 21]. Besides the theoretical applications of solitons, they play an increasingly important role in technology, e.g., in communications [22, 23, 24].

Since the solitons are not isolated objects in most physical systems, their interaction with other fields has been subject to intense research in the literature. Boson and Dirac fields interacting with a soliton are known to affect or even create many intriguing phenomena including vacuum polarization and Casimir effect [25, 26], superconductivity and Bose-Einstein condensation [27, 28], localization of fermions in the braneworld scenarios [29], charge and fermion number fractionalization [30] as well as conducting polymers [31]. Massless Dirac fermions behave as the quasiparticles in materials such as graphene and topological insulators [32, 33].

Exactly solvable models are considered indispensable tools to explore the physics of a system and the symmetries behind it. In this paper, we introduce a parity breaking model with an analytical soliton solution. The potential includes powers up to sixth order in the scalar field $\phi$, where odd powers of $\phi$ exist alongside the even, causing a parity asymmetry. In [34], the authors considered a massless Dirac field interacting with a skyrmion-like planar defect in a system that does not respect the parity symmetry. They studied the fermion bound spectrum as well as the scattering of fermions from the localized topological structure and found a closed form for the scattering cross-section for small fermion-skyrmion coupling. Parity or inversion symmetry breaking models with topological solutions are of importance in many areas of physics, for example in the context of superconductivity [35, 36, 37, 38, 39], fractional quantum Hall effect [40], mesoscopic electron transport [41], current of abnormal parity [42], heavy-ion collisions [43], nonlinear Schrödinger equations [44] and hydrodynamics [45].

In this paper, we consider a parity-breaking model with two minima where the soliton solutions connect the two vacua in a nonsymmetric form. Unlike the kink of $\phi^{4}$ theory, they have a power-law tail at one side and exponential asymptotics at the other. This behavior can be found most frequently in models where the potential has higher than sextic power in the scalar field [46]. These types of solitons are especially interesting in the context of the soliton-soliton interactions (see, e.g., [47, 48, 49, 50]). Although this is not the main focus of this work, we will comment on it when we find it relevant. The goal here is to find the soliton solutions and stability equation analytically as well as to study the interaction of the soliton with boson and fermion fields. We consider two types of interactions with the boson field and Yukawa interaction for the fermion field. The boson bound and scattering states, as well as the fermion zero mode, are expressed in closed analytical forms. However, the massive fermion bound states and energy spectrum are solved numerically. Most analytical results are expressed in terms of the Lambert W function. In Sec. 2, we introduce the model, find the corresponding topological solution and analyze the small oscillations of the soliton. In Sec. 3, we study the interaction of boson and fermion fields with the soliton of our model. Finally, in Sec. 4, we summarize the results of the current work. The appendices provide the details of the calculations.

2 Model

We propose the theory described by the following Lagragian

[TABLE]

where the potential term is given by

[TABLE]

The potential presents two minima, $\phi_{0}=\pm 1$, which allows one to obtain solitonic solutions interpolating between them. Although the potential is sixth-order, it has no parity symmetry since odd powers of $\phi$ are also included, unlike the classical $\phi^{6}$ theory. An equivalent potential could be considered by mapping $\phi\to-\phi$, resulting in an interchange of the roles of the kink and antikink solutions. The coupling $\lambda$ has mass dimension one, defining a natural mass scale in the system, which we use to rescale all the parameters. Nevertheless, when deemed relevant, we explicit the mass dimension as a function of $\lambda$.

Although Lagragian (1) yields a second order equation of motion, thanks to the BPS condition one can obtain an equivalent first order equation

[TABLE]

where the field $\phi$ is static. Integrating the above equation we find

[TABLE]

where $C$ is the integration constant. We choose the center of the soliton at $\phi(0)=0$, implying $C=\frac{1}{4}(2-i\pi)$. This can be solved in terms of Lambert $W$ function111For the properties of Lambert $W$ function check, e.g., [51].. The details of this calculation are provided in Appendix A. The solution is

[TABLE]

The corresponding antikink solution can be obtained by mapping $\phi\to-\phi$. Figure 1a shows the kink profiles for $\phi^{4}$ and our models. Notice that in the kink profile of our model the parity is explicitly broken.

At large $x$ the behavior of the kink is as follows

[TABLE]

The above asymptotic behavior means that the kink at large $x$ has a long-range power-law fall-off in contrast with the opposite tail, $x\to-\infty$, with exponential asymptotics.

Using the BPS condition, it is straightforward to calculate the energy of the soliton configuration, the so-called classical mass of the soliton,

[TABLE]

where the energy density ${\mathcal{E}}(x)$ is shown in Fig. 1b for our model and $\phi^{4}$ kink. Interestingly, despite the difference in the energy density of the two models, the resulting mass is the same.

Having the profile of the soliton, it is relevant to analyze the small fluctuations of the boson field described by the linear stability equation

[TABLE]

with the stability potential

[TABLE]

where $\eta_{n}$’s are the normal modes of the fluctuations around the static solution. Due to the translational symmetry of the system there exists a zero mode, $\omega_{0}=0$. It is possible to show that it is as follows

[TABLE]

Knowing that

[TABLE]

and requiring $\omega_{n}^{2}$ to be non-negative, implies that the zero mode is the only discrete mode, i.e., there are no bound oscilation modes of the soliton apart from translation. Figure 2 shows the stability potential $U(x)$ (panel (a)) as well as the zero mode $\eta_{0}$ (panel (b)). As one can see, the potential presents different limits at $x\to\pm\infty$. Due to this fact, only waves whose energy exceeds $U(-\infty)$ are permitted when travelling from the left. In contrast, incoming waves from the right are allowed for lower energies starting from 0, the value of $U(+\infty)$. In this case, oscillations with an energy smaller than $U(-\infty)$ are totally reflected by the potential barrier.

Until now, we have been concerned with the soliton solutions in isolation. In what follows, we consider the interaction of the soliton of our model with other fields, including boson and fermion fields. We analyze two different types of couplings responsible for the soliton-boson interaction, and a Yukawa coupling between the soliton and the Dirac field. In all three cases, we consider the soliton a background field.

3 Interaction with a scalar field

3.1 Model I

First, let us consider the interaction of a real massive scalar field $\chi$ with the soliton of our model in the following form

[TABLE]

where $m$ is the mass of the field $\chi$ and $g$ is the scalar-soliton coupling constant. This interaction yields a non-homogeneous Klein-Gordon equation

[TABLE]

Separating the time dependence as $\chi=\chi_{s}e^{-iEt}$, we find the equation

[TABLE]

where we define $k^{2}\equiv E^{2}-m^{2}$. First consider the bound states, for which we have $k^{2}<0$. The solution to the equation of motion, eq. (14), is

[TABLE]

introducing $k^{2}\equiv-\kappa^{2}$. The first two terms come from the solution of the homogeneous equation and the last one is a particular solution. Focusing only on the integral term in the above solution and performing the change of variables $u=x-y$ results in

[TABLE]

Two more changes of variables, $v=e^{1+4(x-u)}$ followed by $w=W[v]$, allow us to rewrite the integral in a form that can be directly solved

[TABLE]

Therefore, the general solution takes the form

[TABLE]

To obtain a real $\chi_{s}$ one has to impose the restriction $\kappa=4n$ where $n$ is an integer number. This means that $E^{2}=m^{2}-16n^{2}$ and also $n<m/4(\lambda)$ following the fact that $E^{2}$ is non-negative. Looking at the limit of $\chi_{s}$ when $x\to+\infty$ it is easy to see that, for the solution to be finite, $A$ should be null. At this limit, the last term in the above expression vanishes and therefore the solution converges to $-g/\kappa^{2}$. It remains to determine the value of $B$ which can be found by requiring the solution to be finite when $x\to-\infty$. Doing so, it can be shown that

[TABLE]

The detailed calculations are provided in Appendix B. Figure 3 shows the bound states for three values of $n$, 1, 2 and 3. As it can be seen, at the limits $x\to\pm\infty$ the solution converges to $\mp\frac{g}{\kappa^{2}}=\mp\frac{g}{16n^{2}}$. This result is expected considering eq. (14) when $\phi(x\to\pm\infty)=\pm 1$.

Notice that, interstingly, bound states are also solitons which means that the original soliton can trap another boson field in the form of a soliton configuration.

Now, let us look at the case $k^{2}>0$ which corresponds to the scattering states. The general solution for this equation is in the form

[TABLE]

where again the first two terms come from the solution of the homogeneous equation and the last one is a particular solution. Following the same series of change of variables the integral in the above expression changes to

[TABLE]

where $c.c.$ stands for the complex conjugate. After some simplifications, the solution (20) takes the form

[TABLE]

To verify the result, one can look at the limits $x\to\pm\infty$. At the limit $x\to+\infty$, the last term in the above solution tends to zero and we recover the expected result using eq. (14) when $\phi\to 1$

[TABLE]

The same goes for the limit $x\to-\infty$ where the eq. (3.1) tends to

[TABLE]

The last term can be removed through a redefinition of the coefficients $A$ and $B$, which gives the expected result using eq. (14) when $\phi\to-1$.

3.2 Model II

Now we introduce a different type of coupling between the soliton field $\phi$ and the scalar field $\chi$. Consider the following Lagragian

[TABLE]

where the coupling between the fields is analogous to a Yukawa interaction. This interaction yields the equation of motion

[TABLE]

Considering $\chi=\chi_{s}e^{-iEt}$ and rearranging the terms we arrive at

[TABLE]

Replacing the solitonic solution of our model in the above equation leads to

[TABLE]

which has the formal structure of the Schrödinger equation with energy equal to $(k^{2}+2g)$. Figure 4 shows the form of the potential term in the above Schrödinger-like equation. In [52], the author solved a similar equation. To map our system to the quantum mechanical system solved in the aforementioned paper, we need first to consider the change of variables $1+4x\to-y$ which results in

[TABLE]

Now, the map between their system and ours is given by $2m/\hbar^{2}\to 1$, $V_{0}\to g/4$, $E\to(k^{2}+2g)/16$ and $\sigma\to 1$. As a result, the solution to our system is in the following form

[TABLE]

with $z=W\left[e^{1+4x}\right]$, $\delta^{\pm}=\frac{1}{2}\sqrt{k^{2}\pm 2g}$, $a=\left(\delta^{+}+\delta^{-}\right)^{2}/(4\delta^{+})$,

[TABLE]

where $C_{1}$ and $C_{2}$ are constants and ${}_{1}F_{1}$ and $U$ are the Kummer and Tricomi confluent hypergeometric functions, respectively. The system does not have any bound state, which is easy to recognize from the form of the potential (see Fig. 4). The scattering states from the right and the left are shown in Fig. 5a. Besides that, in the same figure, one can see the scattering from the right where the energy is beneath the threshold required to surpass the barrier. In this case, the wave is totally reflected. Figure 5b shows the reflection coefficient as a function of the momentum for the waves from the right and left. In the case of incoming waves from the right, the reflection coefficient is 1 for momenta associated with energies below the barrier, as expected. Also, for waves coming from both directions, the reflection coefficient drops to zero at high energies since the wave does not see the barrier.

4 Interaction with a fermion field

Fermions can be coupled to the soliton in various ways. We introduce a fermion field $\psi$ coupled to the soliton through a Yukawa coupling in the following form

[TABLE]

where $g$ is a coupling constant. The resulting equation of motion in the background of the soliton reads

[TABLE]

Opening the spinor field $\psi$ in components as $\psi=e^{-iEt}\begin{pmatrix}\psi_{1}\\ \psi_{2}\end{pmatrix}$ one can find the pair of equations

[TABLE]

where the representation for the Dirac matrices is taken as $\gamma^{0}=\sigma_{1}$, $\gamma^{1}=i\sigma_{3}$ and $\gamma^{5}=\sigma_{2}$. For the case of $\phi^{4}$ model, a zero energy bound state or zero mode is known to exist, which is also the case for our model. The zero mode is given by

[TABLE]

where $\mathcal{N}$ is the normalization constant. Since one of the components is null the soliton never receives backreaction from this state and the solution is exact [53]. Performing the above integration we can obtain an explicit solution to the state

[TABLE]

Performing the change of variables $y=\exp(1+4x)$ this becomes

[TABLE]

using the property of Lambert $W$ function

[TABLE]

Therefore, the wavefunction becomes

[TABLE]

with the normalization constant

[TABLE]

Details of the calculation above are supplied in appendix C. In Fig. 6, we show the fermionic zero mode for two different values of the coupling. The resulting parity asymmetry in the zero mode of our model is visible, especially when the fermion-soliton coupling $g$ increases.

Besides that, we solve the equations of motion in (33) for nonzero bound states numerically where the result for the upper and lower components of the first and second fermionic bound states is presented in Fig. 7. A close inspection of the figure reveals that the states do not respect parity. Moreover, we plot the bound and threshold energies as a function of the bound state number as well as the fermion-soliton coupling $g$ in Fig. 8, where the system is solved numerically. It is not difficult to show that the system has energy-reflection symmetry, which is given by $\gamma^{1}$ in our model. In Fig. 8 the symmetry manifests itself by the symmetric form of the spectrum around $E=0$ line. For very small values of $g$, the only discrete mode is the zero mode. However, gradually increasing $g$ from zero supports more and more bound states. Besides the bound states, one can explore the scattering ones considering energies above the threshold in the equation of motion (33). We show the upper and lower components of the fermionic scattering states for the scattering from both directions in Fig. 9. Again, it is easy to observe that the states do not respect parity symmetry.

5 Conclusion

In this work, we have designed a parity-breaking solitonic model where the potential is up to sixth order in the scalar field $\phi$, with two minima. The soliton solutions connecting the two minima in a nonsymmetric form, with one long-range power-law tail and one exponential asymptotics, has been solved in terms of the Lambert W function. Although the system lacks $Z_{2}$ symmetry, changing $\phi\to-\phi$ only swaps the role of the soliton and antisoliton solutions. We have found the soliton mass, which is equal to the one for the kink of $\phi^{4}$ theory, despite a very different energy density. Studying the linear stability equation for the small perturbations around the static soliton solutions, we have concluded that the only discrete mode is the zero mode associated with the translational invariance, in contrast with the parity-symmetric $\phi^{4}$ model. Besides that, we have studied the interaction of the boson and fermion fields with the soliton considering two different types of interaction terms for the bosonic one and the Yukawa interaction for the fermionic one. The first interaction we have examined has led to a non-homogeneous Klein-Gordon equation with interesting results. For example, we have shown that the boson bound state also acquires the form of a defect, which means that the soliton in our model traps the bosonic field in a kink configuration. Considering the second interaction term, a Yukawa-like interaction, we have shown that one can write the equation of motion in the form of a Schrödinger equation. With a change of variables and mapping the parameters with the results obtained in [52], we have found the bound and continuum states analytically. We have also studied the scattering of the waves from the left and right as well as the reflection coefficient, knowing the barrier shape potential term. We have shown that the reflection coefficient is unity for waves with energies beneath the barrier and goes to zero at high-energy, as expected. We have solved the system analytically for both types of interactions. It is not common to find systems that can be fully solved analytically, and this makes the model more valuable for follow-up studies and applications. Moreover, in both cases, we have verified that the results match the expectations in the limiting cases where $\phi(x\to\pm\infty)\to\pm 1$. Finally, the interaction of the fermion field with the soliton has been considered. In this case, we have been able to find the normalized fermion zero mode analytically. We have also obtained the nonzero bound energy spectrum as a function of bound state number as well as the fermion-soliton coupling $g$ numerically. The system has energy reflection symmetry given by $\gamma^{1}$ resulting in a symmetric bound energy spectrum. We have shown that, for very small values of the coupling, the only discrete mode is the zero mode, with a growing number of bound states appearing as we gradually increase the coupling. Finally, the scattering oscillating modes for the waves coming from the left and the right have been shown.

In future work, we plan to apply the model proposed here to address the soliton-soliton long-range interactions, taking advantage of the analytical properties of the model.

Acknowledgments

AA would like to thank CAPES for financial support under the PNPD fellowship. AM acknowledges the financial support from CNPq (process number 305893/2017-3), CAPES and Universidade Federal de Pernambuco Edital Qualis A. AA and AM are thankful to Carlos Batista for fruitful discussions.

Appendix A Calculation of the soliton solution

Performing the change of variables $\chi=\phi-1$ in eq. 4 we obtain

[TABLE]

which leads to

[TABLE]

Recalling that the Lambert $W$ function is defined as the inverse function of

[TABLE]

the previous equation can be written as

[TABLE]

Now solving for $\chi$ and reintroducing $\phi$ we have the final result

[TABLE]

Appendix B Integration constants of the bound states in model I

We start with the general form of the solution for model I, eq. (3.1),

[TABLE]

where

[TABLE]

Let us first look at the limit $x\to+\infty$. At this limit, $f_{n}(x)$ becomes

[TABLE]

Therefore,

[TABLE]

and similarly

[TABLE]

Since the $A$ term diverges in this limit and there is no other term to compensate it, the constant $A$ should be set to zero.

Now considering the limit $x\to-\infty$ we can determine the remaining constant $B$. Using the expansion of the Lambert function for small arguments we have

[TABLE]

In the limit $x\to-\infty$, we can ignore the second term and replace $e^{n\,e^{1+4x}}$ with 1, which gives

[TABLE]

We need to deal with $f_{-n}(x)$ differently. For this function, we have

[TABLE]

where $\gamma(s,z)$ is the lower incomplete gamma function. Therefore, we obtain

[TABLE]

The first term diverges at $x\to-\infty$ and should be cancelled by the $B$ term in the full solution. As a result, eq. (46) becomes

[TABLE]

Appendix C Normalization of the fermionic zero mode

In eq. (38), we can find the normalization factor in the following way

[TABLE]

Choosing the transformation $y=e^{1+4x}$, it results in

[TABLE]

Now let’s consider the change of variables $w=W(y)$ (notice that $y=we^{w}$, by the very definition of the Lambert $W$ function). Therefore,

[TABLE]

which leads to

[TABLE]