Dynamical Behaviours of Matter Wave Positons in Bose-Einstein Condensates

1. Introduction

Theoretical investigations into the nonlinear collective excitations of matter waves have emerged as a highly intriguing and pertinent field, especially in light of experimental observations of Bose-Einstein condensation (BEC) in vapors of alkali-metal atoms [1,2]. Among the captivating manifestations of localized waves in atomic matter, solitons hold particular interest. The concept of solitons was initially introduced to describe nonlinear solitary waves that exhibit remarkable properties, such as non-dispersive behavior, preserving their localized form and speeds both during propagation and after collisions [3,4,5,6]. These inherent advantages of solitons have sparked significant interest in the study of nonlinear systems across various fields of physics, particularly in high-rate telecommunications involving optical fibers, fluid dynamics, capillary waves, hydrodynamics, plasma physics and so on [4,7,8,9]. In addition to this, BECs exhibit the emergence of Faraday and resonant density waves when subjected to harmonic driving [10]. The characteristics of density waves in dipolar condensates at absolute zero temperature using both mean-field variational and full numerical approaches has been investigated. The breaking of symmetry resulting from the anisotropy of the dipole-dipole interaction was found to be a crucial factor in this phenomenon.

From an experimental perspective, precise control over the existence of matter waves in a BEC system can be achieved by effectively manipulating the nonlinear atom interactions through the Feshbach resonance technique, as well as by varying the external trap potential [11,12,13,14]. This flexibility permits us to consider that the coefficient of nonlinearity and the external potential terms in the Gross-Pitaevskii (GP) equation/generalized nonlinear Schrödinger (NLS) equation can vary as functions of both time and/or space. Consequently, investigating the distinctive features of matter waves (solitons, breathers and rogue waves) becomes highly intriguing, given their spatial and temporal localization, particularly in the context of BEC experiments. Motivated by these achievements, extensive research has been dedicated to investigating localized matter waves within quasi-one-dimensional BECs [14,15,16,17,18,19,20,21]. Furthermore, studies focusing on the varible coefficient NLS equation have unveiled the potential to manipulate and enhance these localized density profiles through the utilization of inhomogeneity parameters [22,23,24,25]. The analysis of soliton propagation in optical and the condensed matter systems with $\mathcal{PT}$-symmetry, particularly in inhomogeneous setups, has gained significant attention. Consequently, considerable efforts have been devoted to showcasing the existence of stable bright solitons, dark solitons, and vortices within the NLS equation featuring $\mathcal{PT}$-symmetric potentials [26,27,28]. In the case of weakly interacting toroidal BECs, the occurrence of rotational fluxons (commonly known as Josephson vortices) is linked to the spontaneous disruption of the rotational symmetry within the tunneling superflows [29]. To explore the impact of controllable symmetry breaking on the resulting state of merged counter-propagating superflows, a weakly dissipative mean-field model was employed. In line with this research trajectory, our aim is to construct an intriguing type of localized solution known as positons within the GP equation. We further endeavor to explore the effects of time-dependent modulation of nonlinearity parameters on the characteristics of positon profiles.

Positons, unlike exponentially decaying soliton solutions, are weakly localized nonlinear waves that hold significant importance in the field of nonlinear physics [30,31,32,33]. These solutions are obtained by constraining degenerate eigenvalues within the widely recognized N-soliton algorithm. For positon solutions, the corresponding eigenvalue in the spectral problem is positive and lies within the continuous spectrum. It has been observed that when two positons collide, they retain their individual identities, whereas the soliton remains unchanged following a collision with a positon. However, the positon experiences an influence from the carrier wave and envelope, resulting in a finite phase shift [34,35]. Notably, Matveev’s positon solution to the Korteweg-de Vries (KdV) equation exhibited a spectral singularity [32]. Building on this pioneering work, positon solutions have been derived for other nonlinear evolution equations successfully [36,37,38,39,40]. Recent efforts by Cen et al. introduced the concept of smooth positons or degenerate soliton solutions by allowing the spectral parameter to take complex values [41,42], thereby eliminating the singularity in the KdV equation. Following these advancements, endeavors have been made to construct smooth positon solutions for various equations including the focusing mKdV equation [43], complex mKdV equation [46], derivative NLS equation [44,45], NLS-Maxwell-Bloch equation [47], higher-order Chen-Lee-Liu equation [49] and Gerdjikov-Ivanov equation [48]. More recently, smooth positons and breather positons have been derived for the generalized NLS equation with higher-order nonlinearity, along with higher-order solutions for an extended NLS equation featuring cubic and quartic nonlinearity [50,51]. Inspired by these advancements in the field of positons, our research aims to construct positon solutions within the GP equation, incorporating time-varying nonlinearity and trap potential.

The crucial step in this attempt involves utilizing the similarity transformation on a meticulously chosen ansatz solution. This transformation effectively converts the variable parameter GP equation into the conventional NLS equation with constant coefficients. By implementing this transformation, the modified variables allow us to derive new solutions for the considered equation by expressing the known positon solutions in the altered coordinate system. By leveraging the combination of known positons of typical NLS equation solutions with similarity transformation functions, one can derive novel (non-autonomous positon) solutions for the GP equation. The integrability requirements, which establish the relationship between variable parameters (modulated nonlinearity and trap potential), and the proposed ansatz solution, serve as the underlying considerations for this procedure. However, despite the associated costs, this approach holds significant value as it not only reveals new analytical solutions for the GP equation but also empowers users to control the outcomes by judiciously selecting appropriate nonlinearity strengths and trap potentials.

Motivated by the experimental feasibility of studying BECs, our research focuses on exploring the dynamical characteristics of positons. To achieve this objective, we construct second- and third-order matter wave positon solutions for the one-dimensional GP equation, considering a variable nonlinearity parameter and an external trap potential. The construction of these solutions involves transforming the time-modulated GP equation into a ccNLS equation using a similarity transformation. We establish that the trapping potential and nonlinearity modulated parameter must satisfy a constraint for the considered equation to be integrable and yield the desired solutions. By leveraging the known smooth positon solutions (second- and third-order) of the constant coefficient NLS equation, we present matter wave smooth positon solutions of the GP equation. We investigate the deformation of positon density profiles with respect to three different forms of variable nonlinearity parameters, namely (i) kink-like nonlinearity $R\left(t\right)=R_0+R_{1}tanh(R_{2}t+R_3)$, (ii) localized or sech-type nonlinearity $R\left(t\right)=R_0+R_1\mathrm{sech}(R_{2}t+R_3)$, and (iii) periodic nonlineaity $R\left(t\right)=R_0+R_{1}sin(R_{2}t+R_3)$, where $R_0$, $R_1$, $R_2$ and $R_3$ are arbitrary parameters. Our findings reveal that a range of nonlinear physical phenomena, including stretching, curving, annihilation, breathing, oscillating, enhancement, and suppression, are manifested in the underlying matter wave positon density profiles. When considering a kink-like modulated nonlinearity, the position density profiles of second- and third-order smooth matter wave positions experience stretching, while their amplitudes can either be enhanced or suppressed. It is important to note that these profiles vanish during different time intervals, with disappearance occurring for $t<0$ when the parameter $R_2$ assumes positive values, and for $t>0$ when $R_2$ takes negative values. In the case of a localized or sech-type modulated nonlinearity, the positon density profiles become compressed and curved within the background density of the condensate. For periodic modulated nonlinearity, positons exhibit a periodic behavior, and adjusting the strengths of nonlinearity leads to an increase in their periodicity, as observed in our analysis. This observation provides valuable insights for experimentalists analyzing novel density profiles in BECs.

We structure our work as follows. In Section 2, we take into account the GP equation with time-modulated nonlinearity and trap potentials. The second and third-order smooth positon solutions are deduced for this equation using the similarity transformation. The integrable requirement between the modulated nonlinearity and trap potential is obtained while employing the integrable technique to the considered equation. In Section 3, by suitably choosing the different forms of variable nonlinearity parameter, we explore the various dynamical characteristics in the density of matter wave positon profiles. Finally, in Section 4, we conclude our observations.

2. BEC Model and Similarity transformation

The behavior of a BEC confined within an external potential can be effectively characterized using the renowned NLS equation derived from mean field theory, commonly referred to as the GP equation. In the specific scenario of a cigar-shaped trapping potential, where simplicity and physical significance coincide, the radial degree of freedom in the three-dimensional GP equation can be eliminated through integration, leading to the derivation of a dimensionless quasi-one-dimensional equation [1,2,5,14,17]

$$\begin{array}{c}\hfill i\frac{\partial \psi }{\partial t}+\frac{1}{2}\frac{{\partial }^2\psi }{\partial x^2}+R\left(t\right){\left|\psi \right|}^2\psi +\frac{1}{2}{\lambda }^2\left(t\right)x^2\psi =0,\end{array}$$

(1)

where $\psi (x,t)$ is the condensate wave function. In Equation (1), t and x are time and spatial coordinates that are expressed in units ${\omega }_{\perp }^{-1}$ and $a_{\perp }=\sqrt{\frac{\hslash }{\left(m{\omega }_{\perp }\right)}}$, m is the atomic mass, respectively. The atom-atom interaction term (scattering length between atoms) denoted by the representation, $R\left(t\right)$, is tuneable using the Feshbach resonance (FR). In a series of exquisite experiments using sodium and rubidium condensates, Feshbach resonances were investigated. They have also been employed in a variety of significant experimental studies, such as the creation of bright and dark matter-wave solitons, among others. The time-modulated trap potential is described by ${\displaystyle {\lambda }^2\left(t\right)=\frac{{\omega }^2\left(t\right)}{{\omega }_{\perp }}}$, where $\omega$ and ${\omega }_{\perp }$ are the trap frequency in the axial direction and the radial trap frequency, respectively. In the study of trapped BECs, the trap frequency along the elongated axis, denoted as $\omega$, has been intentionally selected to vary with time, t, in order to investigate the characteristics of BECs within the trap. As a result, both the coefficient of nonlinearity (R) and the potential parameter ($\lambda$) can exhibit time dependence. By appropriately choosing these two time-dependent parameters, the GP equation (1) can effectively capture the dynamics and manipulation of BECs. These parameters serve as powerful tools for controlling and manipulating localized matter waves in BECs, achieved through the adjustment of external magnetic fields and the optically controlled interactions using techniques such as the FR method [13,17,20].

To study the matter wave positons in (1), we adopt the similarity transformation mentioned below to map the time-modulated GP equation (1) to the ccNLS equation [5,17,26,28]:

$$\begin{array}{c}\hfill \psi (x,t)=s\left(t\right)\varphi (\eta ,\tau )exp\left[i\theta \right(x,t\left)\right].\end{array}$$

(2)

In Equation (2), the unknown functions, namely $s\left(t\right),\eta (x,t),\tau \left(t\right)$ and $\theta (x,t)$ are the amplitude, similarity spatial variable, the dimensionless time and the phase factor, respectively which are to be computed. Upon involving the substitution of (2) into (1), we obtain the following set of partial differential equations that are related to the unknown functions, such as

$$\begin{array}{cc}\hfill {\eta }_{xx}& =0,\hfill \\ \hfill {\eta }_t+{\eta }_x{\theta }_x& =0,\hfill \\ \hfill {\eta }_x^2-R\left(t\right)s^2\left(t\right)& =0,\hfill \\ \hfill s^{\prime }\left(t\right)+s\left(t\right)\frac{1}{2}{\theta }_{xx}& =0,\hfill \\ \hfill {\tau }^{\prime }\left(t\right)-R\left(t\right)s^2\left(t\right)& =0,\hfill \\ \hfill {\theta }^{\prime }\left(t\right)+\frac{1}{2}{\theta }_x^2-\frac{{\lambda }^2\left(t\right)}{2}{x}^2& =0.\hfill \end{array}$$

(3)

The explicit expressions of the unknown functions can be acquired by solving the aforementioned set of equations, and they take the form

$$\begin{array}{cc}\hfill s\left(t\right)=& s_0\sqrt{R\left(t\right)},\hfill \\ \hfill \eta (x,t)=& s_{0}R\left(t\right)x-b{s}_0^3\int R^2\left(t\right)dt,\hfill \\ \hfill \theta (x,t)=& -\frac{R{\left(t\right)}_t}{2R\left(t\right)}{x}^2+b{s}_0^{2}R\left(t\right)x-\frac{1}{2}{b}^2{s}_0^4\int R^2\left(t\right)dt,\hfill \\ \hfill \tau \left(t\right)=& \frac{1}{2}{s}_0^2\int R^2\left(t\right)dt.\hfill \end{array}$$

(4)

where b and $s_0$ are arbitrary constants. Additionally, we have found an integrability condition that imposes a connection between the time-modulated nonlinearity and the trap potential parameter, as shown by: [5,17,28]

$$\frac{d}{dt}\left(\frac{R_t}{R}\right)-{\left(\frac{R_t}{R}\right)}^2+{\lambda }^2\left(t\right)=0.$$

(5)

Using this equation (5), one can find the $\lambda \left(t\right)$ by fixing the $R\left(t\right)$ and vice versa. In this work, we consider physically intriguing function $R\left(t\right)$ and determine the $\lambda \left(t\right)$ by the expression, that is ${\displaystyle \lambda \left(t\right)=\frac{\sqrt{R^{\prime }{\left(t\right)}^2+R^{\prime }\left(t\right)-R\left(t\right)R^{\prime \prime }\left(t\right)}}{R\left(t\right)}}$.

In Equation (2), the function $\varphi (\eta ,\tau )$ is found to fulfill the ccNLS equation as

$$\begin{array}{c}\hfill i\frac{\partial \varphi }{\partial \tau }+\frac{1}{2}\frac{{\partial }^2\varphi }{\partial {\eta }^2}+{\left|\varphi \right|}^2\varphi =0.\end{array}$$

(6)

The equation under consideration (6) exhibits a wide range of localized solutions, including solitons, breathers, rogue waves, and their corresponding profiles. In this study, we focus on the positon solutions of the NLS equation to investigate matter wave positons in quasi-one-dimensional BECs.

By selecting an appropriate functional form for the time-modulated nonlinearity function $R\left(t\right)$ while ensuring the satisfaction of condition (5), we can derive matter wave positon solutions for the GP equation (1) in the following form:

$$\begin{array}{cc}\hfill \psi (x,t)=& s_0\sqrt{R\left(t\right)}\varphi (\eta ,\tau )exp\left[i\left(-\frac{1}{2}{b}^2{s}_0^4\int R{\left(t\right)}^{2}dt+b{s}_0^{2}xR\left(t\right)-\frac{x^2{R}^{\prime }\left(t\right)}{2R\left(t\right)}\right)\right],\hfill \end{array}$$

(7)

where $\varphi (\eta ,\tau )$ is the solution of the ccNLS equation (6). The solution (7) has the potential to generate a multitude of novel positon structures that could be experimentally realized. To summarize the current progress, one can generate several solutions (positons) for the GP model (1) by first obtaining solutions for the ccNLS equation (6) while satisfying the mentioned relationships. An intriguing advantage and potential perspective of the similarity transformation is worth emphasizing, as it allows the extension of this approach to models featuring variable nonlinearity and external trap potential coefficients dependent on both longitudinal and spatial coordinates. By appropriately tailoring and imposing constraints, the resulting dynamics of physical systems can be attainable. Consequently, in Sec. III, we construct positon solutions for the ccNLS equation (6) to analyze matter wave positons in quasi-one-dimensional BECs.

3. Deformation of matter wave positons in BECs

In this section, we intend to investigate the dynamical characteristics of matter wave positons using the solution (7) by considering three different forms of time-varying nonlinearity parameter and the associated trap potentials. For this investigation, in the following, we first derive the second-order and third-order smooth positon solution of the ccNLS equation (6).

3.1. Second-order smooth matter wave positons

Now, we explore the various novel density profiles of second-order smooth matter wave positons in BECs. For this objective, we first derive the second-order smooth positon solution for the ccNLS Equation (6) as [50]

$$\begin{array}{c}\hfill \varphi \left[2\right]=\frac{P_{11}}{Q_{11}},\end{array}$$

(8)

where

$$\begin{array}{cc}\hfill P_{11}=& 4\left({\alpha }_1-{\alpha }_1^{*}\right)e^{2i{\alpha }_1\left(2{\alpha }_1\tau +\eta \right)}\left(-{\alpha }_1\left(\eta -4{\alpha }_1^{*}\tau \right)+{\alpha }_1^{*}\eta -4{\alpha }_1^2\tau -i\right)\hfill \\ \hfill & +4\left({\alpha }_1^{*}-{\alpha }_1\right)e^{2i{\alpha }_1^{*}\left(2{\alpha }_1^{*}\tau +\eta \right)}\left(-{\alpha }_1\left(4{\alpha }_1^{*}\tau +\eta \right)+{\alpha }_1^{*}\eta +4\left({\alpha }_1^{*}\right){}^2\tau +i\right),\hfill \\ \hfill Q_{11}=& -2e^{\left(2i\left({\alpha }_1+{\alpha }_1^{*}\right)\eta +4i\left({\alpha }_1^2+{\alpha }_1^{*2}\right)\tau \right)}\left(-1+2\left({\alpha }_1-{\alpha }_1^{*2}\right){\eta }^2+32{\alpha }_1{\alpha }_1^{*}{\left({\alpha }_1-{\alpha }_1^{*}\right)}^2{\tau }^2\right.\hfill \\ & \left.+8\eta \tau \left({\alpha }_1+{\alpha }_1^{*}\right)\left({\alpha }_1-{\alpha }_1^{*2}\right)\right)+e^{4i{\alpha }_1\left(2{\alpha }_1\tau +\eta \right)}+e^{4i{\alpha }_1^{*}\left(2{\alpha }_1^{*}\tau +\eta \right)},\hfill \end{array}$$

(9)

where ${\alpha }_1$ is the eigenvalue of the spectral parameter, ${\alpha }_1^{*}$ is the complex conjugate of ${\alpha }_1$, $\eta$ and $\tau$ are provided in Equation (4). Substituting this solution into (7) along with the suitable form of $R\left(t\right)$, we obtain the second-order smooth matter wave positon solution of (1).

Now, utilizing the above mentioned solution, we move to investigate its dynamical evolutions through three different forms of time-modulated nonlinearity parameter, such as (i) kink-like nonlinearity $R\left(t\right)=R_0+R_{1}tanh(R_{2}t+R_3)$, (ii) localized or sech-type nonlinearity $R\left(t\right)=R_0+R_1\mathrm{sech}(R_{2}t+R_3)$, and (iii) periodic nonlineaity $R\left(t\right)=R_0+R_{1}sin(R_{2}t+R_3)$, where $R_0$, $R_1$, $R_2$ and $R_3$ are arbitrary parameters. In the following, we provide a comprehensive demonstration of the impact of modulated nonlinearity parameters on the positon density profiles.

To begin, we take the kink-like nonlinearity parameter, that is $R\left(t\right)=R_0+R_{1}tanh(R_{2}t+R_3)$ to reveal the novel features in BECs. Substituting this nonlinearity term in the generalized solution (7), we find

$$\begin{array}{cc}\hfill \psi (x,t)=& s_0\sqrt{R_{1}tanh\left(R_{2}t+R_3\right)+R_0}\varphi (\eta ,\tau )exp\left[-\frac{i}{4R_2}\left(b^2{s}_0^4\left(-2R_1^2\right.\right.\right.\hfill \\ & \times tanh\left(R_{2}t+R_3\right)+\left(R_0-R_1\right){}^{2}log\left(tanh\left(R_{2}t+R_3\right)+1\right)-\left(R_0+R_1\right){}^2\hfill \\ & \left.\times log\left(1-tanh\left(R_{2}t+R_3\right)\right)\right)-4b{s}_0^{2}x\left(R_{1}tanh\left(R_{2}t+R_3\right)+R_0\right)\hfill \\ & \left.\left.+\frac{2R_1{R}_2{x}^2{\sec \mathrm{h}}^2\left(R_{2}t+R_3\right)}{R_{1}tanh\left(R_{2}t+R_3\right)+R_0}\right)\right],\hfill \end{array}$$

(10)

where $\varphi (\eta ,\tau )$ is the second-order smooth positon solution of the ccNLS equation (6). Utilizing this solution (10), we examine a thorough analysis of the positon density profiles, investigating their diverse characteristics as we vary the strength of the time-modulated nonlinearity parameter.

Figures 1a-f illustrates the qualitative profile of a second-order smooth matter positon in BECs corresponding to a kink-like nonlinearity modulated parameter $R\left(t\right)=R_0+R_{1}tanh(R_{2}t+R_3)$. By selecting specific parameter values, $R_0=1.05$, $R_1=0.01$, $R_2=1.05$, and $R_3=0.5$, we obtain a well-localized second-order smooth positon profile associated with the eigenvalue ${\alpha }_1=0.2+0.5i$, as depicted in Figure 1a. Notably, the orientation of the smooth positon density profile changes when varying the eigenvalue associated with the solution. For instance, by modifying the eigenvalue to ${\alpha }_1=0.3+0.6i$, Figure 1b clearly demonstrates a shift in the orientation of the positon profile accompanied by an enhancement in its amplitude. Furthermore, we investigate the effects of adjusting the strengths of the nonlinearity parameters, specifically $R_0$, $R_1$, and $R_2$. For instance, when $R_0$ is changed to $1.85$, the positon profile’s orientation relocates while experiencing a slight increase in amplitude, as depicted in Figure 1c. Additionally, Figure 1d illustrates that increasing the value of $R_1$ to $0.5$ leads to the collapse of the condensate profile on one side of the positon profile ($t<0$). Conversely, when $R_1$ takes a negative value, such as $R_1=-0.5$, the reverse phenomenon occurs, resulting in the disappearance of the density profile in the corresponding plane ($t>0$) as represented in Figure 1e. Moreover, increasing the value of $R_2$ to $3.5$ yields a well-localized positon profile that exhibits a compression within the condensate density background. Additionally, as depicted in Figure 1f, the width of the crest of the positon profile widens over time.

Next, we consider the localized-type or sech-type nonlinearity, namely $R\left(t\right)=R_0+R_1\mathrm{sech}(R_{2}t+R_3)$ for investigating the distortion of positon profiles in the condensate density background. Plugging this form of $R\left(t\right)$ in (7), we obtain the matter positon solution in the form

$$\begin{array}{cc}\hfill \psi (x,t)=& s_0\sqrt{R_1\sec \mathrm{h}\left(R_{2}t+R_3\right)+R_0}\varphi (\xi ,\tau )exp\left[\frac{i}{2R_2}\left(-b^2{s}_0^4\left(R_2{R}_0^{2}t\right.\right.\right.\hfill \\ & \left.+R_1^{2}tanh\left(R_{2}t+R_3\right)+2R_1{R}_0{tan}^{-1}\left(sinh\left(R_{2}t+R_3\right)\right)\right)\hfill \\ & \left.\left.+b{s}_0^{2}x\left(R_1\sec \mathrm{h}\left(R_{2}t+R_3\right)+R_0\right)+\frac{R_1{R}_2{x}^{2}tanh\left(R_{2}t+R_3\right)}{2\left(R_{0}cosh\left(R_{2}t+R_3\right)+R_1\right)}\right)\right],\hfill \end{array}$$

(11)

where $\varphi (\eta ,\tau )$ is the positon solution of the ccNLS equation (6).

Figure 2a-f represents the density profile of a second-order smooth matter positon with a modulated nonlinearity function given by $R\left(t\right)=R_0+R_1\mathrm{Sec}\mathrm{h}(R_{2}t+R_3)$. In this study, we consider specific parameter values: $R_0=1.05$, $R_1=0.01$, $R_2=1.05$, and $R_3=0.5$, aiming to investigate the intriguing properties exhibited by the positon density profiles. With these initial parameter values, we obtain the second-order smooth positon profile as shown in Figure 2a. By increasing $R_0$ to $1.85$, we observe a stretching of the positon and an enhancement in its amplitude, as depicted in Figure 2b. Similarly, when we tune the parameter $R_1$ to $0.85$, a curvature appears in the condensate profile, as shown in Figure 2c. Notably, when $R_1$ takes a negative value, we observe the formation of two peaks at the center ($t=0$), accompanied by a suppression in amplitude, as demonstrated in Figure 2d. Additionally, in the case where $R_0=1.85$ and $R_1=1.75$, we observe a gradual increase in amplitude, further stretching of the positon, and the formation of a curved profile, as shown in Figure 2e. At $R_0=2.5$ and $1.5$, a compressed density profile with an identical amplitude is obtained, as depicted in Figure 2f.

Finally, we investigate the influence of a periodically modulated nonlinearity given by $R\left(t\right)=R_0+R_1\mathrm{sech}(R_{2}t+R_3)$ on the positon density profiles. By inserting this $R\left(t\right)$ in (7), we obtain the explicit form of matter wave positon solution given by

$$\begin{array}{cc}\hfill \psi (x,t)=& s_0\sqrt{R_{1}sin\left(R_{2}t+R_3\right)+R_0}\varphi (\xi ,\tau )exp\left(-\frac{i}{8R_2}\left(b^2{s}_0^4\left(4R_0^2\left(R_{2}t+R_3\right)\right.\right.\right.\hfill \\ & \left.+R_1^2\left(2R_{2}t-sin\left(2\left(R_{2}t+R_3\right)\right)+2R_3\right)-8R_1{R}_{0}cos\left(R_{2}t+R_3\right)\right)\hfill \\ & \left.\left.-8b{s}_0^{2}x\left(R_{1}sin\left(R_{2}t+R_3\right)+R_0\right)+\frac{4R_1{R}_2{x}^{2}cos\left(R_{2}t+R_3\right)}{R_{1}sin\left(R_{2}t+R_3\right)+R_0}\right)\right),\hfill \end{array}$$

(12)

where $\varphi (\eta ,\tau )$ is the positon solution of the ccNLS equation (6).

The second-order smooth matter wave positon density profiles are displayed in Figure 3 for periodic nonlinearity modulated function $R\left(t\right)=R_0+R_{1}sin{R}_{2}t+R_3$ in the context of BECs. By appropriately tuning the parameters $R_0$, $R_1$, and $R_2$, we are able to achieve periodic positon solutions. In Figure 3a, using initial parameter values of $R_0=1.5$, $R_1=0.05$, $R_2=1.25$, and $R_3=0.5$, we observe a periodic behavior in the positon condensate profile. Subsequently, when the nonlinearity strength $R_0$ is increased to $2.25$, the oscillation of the positon becomes more elongated and the amplitude is raised, as shown in Figure 3b. Moreover, as we increase the value of $R_1$ to $0.55$, the periodicity of the positon profile increases along with an increase in amplitude which can be seen in 3c. This trend is further amplified when the value of $R_1$ is increased, as depicted in Figure 3d. When the parameters $R_2=2.75$ and $R_3=5.25$ are chosen, the oscillation in the positon profile increases and retains a similar orientation, as evident in Figure 3e and Figure 3f, respectively.

3.2. Third-order smooth matter wave positons

In the previous sub-section, we investigated the modifications of second-order smooth matter wave positon profiles by varying the distributed coefficients (nonlinearity function) in the GP equation (1). In this sub-section, we delve into the deformations of third-order smooth matter wave positons in condensates as we manipulate the strength of the modulated nonlinearity parameter. To explore these intriguing characteristics, we consider the third-order smooth positon solution for the ccNLS equation (6), which takes the form of [50]

$$\begin{array}{c}\hfill \varphi \left[3\right]=\frac{P_2}{Q_2},\end{array}$$

(13)

where

$$\begin{array}{cc}\hfill P_2=& e^{-i\left(2\tau \left({\alpha }_1^2+{\alpha }_1^{*2}\right)+\eta ({\alpha }_1+{\alpha }_1^{*})\right)}\left(P_{21}{e}^{2i({\alpha }_1-{\alpha }_1^{*})(\eta +2\tau ({\alpha }_1+{\alpha }_1^{*}))}\right.\hfill \\ & \left.+P_{22}{e}^{-2i({\alpha }_1-{\alpha }_1^{*})(\eta +2\tau ({\alpha }_1+{\alpha }_1^{*}))}+P_{23}\right),\hfill \end{array}$$

(14)

with

$$\begin{array}{cc}\hfill P_{21}=& -4({\alpha }_1-{\alpha }_1^{*})\left(-3+2\eta (3i+8\tau {\alpha }_1\left({\alpha }_1-{\alpha }_1^{*}\right)+2{\eta }^2{({\alpha }_1-{\alpha }_1^{*})}^2\right.\hfill \\ & \left.+32{\tau }^2{\alpha }_1^2{({\alpha }_1-{\alpha }_1^{*})}^2-4i\tau \left({\alpha }_1-{\alpha }_1^{*}\right)\left(-7{\alpha }_1+{\alpha }_1^{*}\right)\right),\hfill \\ \hfill P_{22}=& -4({\alpha }_1-{\alpha }_1^{*})\left(-3+4i\tau ({\alpha }_1-7{\alpha }_1^{*})({\alpha }_1-{\alpha }_1^{*})+2{\eta }^2{({\alpha }_1-{\alpha }_1^{*})}^2\right.\hfill \\ & \left.+32{\tau }^2{({\alpha }_1-{\alpha }_1^{*})}^2{\alpha }_1^{*2}-2\eta ({\alpha }_1-{\alpha }_1^{*})(3i+8\tau {\alpha }_1^{*}(-{\alpha }_1+{\alpha }_1^{*})\right),\hfill \\ \hfill P_{23}=& 8({\alpha }_1-{\alpha }_1^{*})\left(3+4{\eta }^2{({\alpha }_1-{\alpha }_1^{*})}^2-2x^4{({\alpha }_1-{\alpha }_1^{*})}^4+{\alpha }_1^2-512{\tau }^4{\alpha }_1^2{({\alpha }_1-{\alpha }_1^{*})}^4{\alpha }_1^{*2}\right.\hfill \\ & +8i\eta {\alpha }_1^{*3}-4{\eta }^2{\alpha }_1^{*4}-8\eta {\alpha }_1^3(-i+\eta {\alpha }_1^{*})-2{\alpha }_1{\alpha }_1^{*}(-3+4i\eta {\alpha }_1^{*}+4x^2{\alpha }_1^{*2})\hfill \\ & +16\tau {({\alpha }_1-{\alpha }_1^{*})}^2(-i+i{\eta }^2{({\alpha }_1-{\alpha }_1^{*})}^2+\eta ({\alpha }_1+{\alpha }_1^{*})-{\eta }^3{({\alpha }_1-{\alpha }_1^{*})}^2({\alpha }_1+{\alpha }_1^{*}))\hfill \\ & +8{\tau }^3{({\alpha }_1-{\alpha }_1^{*})}^3({\alpha }_1+{\alpha }_1^{*})(-i{\alpha }_1^{*}+{\alpha }_1(-i+4\eta {\alpha }_1^{*}))\hfill \\ & \left.+8{\tau }^2{({\alpha }_1-{\alpha }_1^{*})}^2({\alpha }_1^2-8i\eta {\alpha }_1^2{\alpha }_1^{*}+{\eta }^2(1+24{\alpha }_1^2{\alpha }_1^{*2}))\right),\hfill \end{array}$$

(15)

and

$$\begin{array}{cc}\hfill Q_2=& 4e^{-3i\left({\alpha }_1-{\alpha }_1^{*}\right)+(\eta +2\tau ({\alpha }_1+{\alpha }_1^{*}))}\left(1+Q_{21}{e}^{4i({\alpha }_1-{\alpha }_1^{*})(\eta +2\tau ({\alpha }_1+{\alpha }_1^{*}))}\right.\hfill \\ & \left.+Q_{22}{e}^{2i({\alpha }_1-{\alpha }_1^{*})(\eta +2\tau ({\alpha }_1+{\alpha }_1^{*}))}+e^{6i({\alpha }_1-{\alpha }_1^{*})(\eta +2\tau ({\alpha }_1+{\alpha }_1^{*}))}\right),\hfill \end{array}$$

(16)

where

$$\begin{array}{cc}\hfill Q_{21}=& 3+1024{\tau }^4{\alpha }_1^2{({\alpha }_1-{\alpha }_1^{*})}^4{\alpha }_1^{*2}+48{\tau }^2{({\alpha }_1-{\alpha }_1^{*})}^2({\alpha }_1^2-6{\alpha }_1{\alpha }_1^{*}+{\alpha }_1^{*2})\hfill \\ \hfill & +4{\eta }^4{({\alpha }_1-{\alpha }_1^{*})}^4-128i{\tau }^3{({\alpha }_1-{\alpha }_1^{*})}^3({\alpha }_1+{\alpha }_1^{*})({\alpha }_1^2-4{\alpha }_1{\alpha }_1^{*}+{\alpha }_1^{*2})\hfill \\ \hfill & +8{\eta }^3{({\alpha }_1-{\alpha }_1^{*})}^3(i+4\tau ({\alpha }_1^2-{\alpha }_1^{*2}))+4{\eta }^2{({\alpha }_1-{\alpha }_1^{*})}^2\left(-3+12i\tau ({\alpha }_1^2-{\alpha }_1^{*2})\right.\hfill \\ & \left.+16{\tau }^2{({\alpha }_1-{\alpha }_1^{*})}^2({\alpha }_1^2+4{\alpha }_1{\alpha }_1^{*}+{\alpha }_1^{*2}\right)-16\eta \tau {({\alpha }_1-{\alpha }_1^{*})}^2\left(3{\alpha }_1^{*}-32{\tau }^2{\alpha }_1^4{\alpha }_1^{*}\right.\hfill \\ & \left.+32{\tau }^2{\alpha }_1^3{\alpha }_1^{*2}+8\tau {\alpha }_1^2{\alpha }_1^{*}(-3i+4\tau {\alpha }_1^{*2})+{\alpha }_1(3+24i\tau {\alpha }_1^{*2}-32{\tau }^2{\alpha }_1^{*4})\right),\hfill \end{array}$$

$$\begin{array}{cc}\hfill Q_{22}=& 3+1024{\tau }^4{\alpha }_1^2{({\alpha }_1-{\alpha }_1^{*})}^4{\alpha }_1^{*2}+48{\tau }^2{({\alpha }_1-{\alpha }_1^{*})}^2({\alpha }_1^2-6{\alpha }_1{\alpha }_1^{*}+{\alpha }_1^{*2})\hfill \\ \hfill & +4{\eta }^4{({\alpha }_1-{\alpha }_1^{*})}^4+128i{\tau }^3{({\alpha }_1-{\alpha }_1^{*})}^3({\alpha }_1+{\alpha }_1^{*})({\alpha }_1^2-4{\alpha }_1{\alpha }_1^{*}+{\alpha }_1^{*2})\hfill \\ \hfill & +8{\eta }^3{({\alpha }_1-{\alpha }_1^{*})}^3(i+4\tau ({\alpha }_1^2-{\alpha }_1^{*2}))+4{\eta }^2{({\alpha }_1-{\alpha }_1^{*})}^2\left(-3-12i\tau ({\alpha }_1^2-{\alpha }_1^{*2})\right.\hfill \\ & \left.+16{\tau }^2{({\alpha }_1-{\alpha }_1^{*})}^2({\alpha }_1^2+4{\alpha }_1{\alpha }_1^{*}+{\alpha }_1^{*2}\right)-16\eta \tau {({\alpha }_1-{\alpha }_1^{*})}^2\left(3{\alpha }_1^{*}-32{\tau }^2{\alpha }_1^4{\alpha }_1^{*}\right.\hfill \\ & \left.+32{\tau }^2{\alpha }_1^3{\alpha }_1^{*2}+8\tau {\alpha }_1^2{\alpha }_1^{*}(3i+4\tau {\alpha }_1^{*2})+{\alpha }_1(3+768i{\tau }^3{\alpha }_1^{*6})\right).\hfill \end{array}$$

By substituting the third-order smooth positon solution of the ccNLS Equation (6), along with the appropriate form of the modulated parameter $R\left(t\right)$, into equation (7), we delve into the analysis of the underlying dynamics of GP equation (1).

Figures 4a-f present the condensate profiles of third-order smooth positons with a kink-like nonlinearity modulated function $R\left(t\right)=R_0+R_{1}tanh(R_{2}t+R_3)$. By considering initial parameters $R_0=1.25$, $R_1=0.01$, $R_2=1.05$, and $R_3=0.5$, an appropriate third-order positon profile is formed within the condensate density background, as shown in Figure 4a. Notably, a decrease in $R_0$ leads to a sudden rise in one of the subcrests, as demonstrated in Figure 4b for $R_0=0.85$. Furthermore, Figure 4c illustrates an increase in amplitude when tuning the value of $R_0$ to $1.75$. On the other hand, in Figure 4d, when the parameter $R_1$ is enhanced to $1.25$, the positon tends to disappear within the condensate profile for $t<0$. Interestingly, a reverse scenario occurs when $R_1=-0.65$, resulting in a higher peak in one of the profiles and an increased amplitude, as depicted in Figure 4e. Similarly, an increase in $R_2$ to $3.75$ yields a compressed three-positon profile with the same amplitude as Figure 4a, as shown in Figure 4f.

In Figures 5a-f, we show density profile of third-order smooth matter positons for sech-like nonlinearity modulated function $R\left(t\right)=R_0+R_1\mathrm{sech}(R_{2}t+R_3)$ for the GP model (1). Taking initial parameters $R_0=1.25$, $R_1=0.01$, $R_2=1.05$, $R_3=0.5$, we obtain the conventional smooth three-positon profile, as shown in Figure 5a. By adjusting $R_0$ to $1.85$, the positon undergoes stretching, resulting in an enhanced amplitude, as depicted in Figure 5b. Similarly, when $R_1$ is modified to $1.5$, the positon exhibits curvature in one of its wave crest, as seen in Figure 5c. Figure 5d illustrates a decrease in the soliton amplitude when $R_1$ is set to $0.35$. However, no significant changes are observed in the smooth three-positon profile when altering the value of R2 compared to Figure 5a. This observation is displayed in Figure 5e. Finally, for $R_0=1.85$ and $R_1=0.8$, the three-positon becomes compressed, accompanied by a decrease in its amplitude, as shown in Figure 5f.

Finally, we investigate the density profile of the third-order smooth matter positon in one-component BECs utilizing a periodic nonlinearity modulated function, denoted as $R\left(t\right)=R_0+R_{1}sin(R_{2}t+R_3)$. Setting the initial nonlinearity strengths as $R_0=1.5$, $R_1=0.05$, $R_2=0.85$, and $R_3=0.5$, we aim to obtain an appropriate periodic third-order smooth positon density profile, as illustrated in Figure 6a. Especially, a remarkable change occurs in the profile when the parameter $R_0$ is selected as $0.85$, resulting in the formation of a sharp ascent on one side of the profile, as depicted in Figure 6b. Similarly, when we increase $R_0$ to $1.75$, the profile returns to its original position, as displayed in Figure 6c. Furthermore, in Figure 6d, we observe that the periodicity of the positon density profiles increases when $R_1$ is set to $0.5$. Moreover, raising the value of $R_2$ to $2.25$ leads to an increase in the oscillation amplitude of the positon, as demonstrated in Figure 6e. Finally, when the value of $R_2$ is further increased to $5.25$, we observe no change in the positon amplitude, but a variation in the oscillation behavior, which is clearly visualized in Figure 6f.

4. Conclusion

In our study, we have derived the second- and third-order smooth matter wave positon solutions of the GP equation. These solutions capture the dynamics of one-component BECs subjected to time-modulated nonlinearity (represented by the effective scattering lengths) and external harmonic trap potentials. Through a similarity transformation technique, we have mapped the time-modulated GP equation onto the ccNLS equation, ensuring an integrability condition between the nonlinearity coefficient and the external trap potential. We have investigated three distinct forms of modulated nonlinearities: (i) kink-like, (ii) localized or sech-like, and (iii) periodic. By varying the parameters associated with the nonlinearity strength, we have observed various nonlinear phenomena in the positon density profiles. These phenomena include stretching, curving, oscillating, breathing, collapsing, amplification, and suppression. In the case of a kink-like modulated nonlinearity, the positon density profiles (represented by the second- and third-order smooth matter wave positons) undergo stretching, while their amplitudes can be enhanced or suppressed. It is noteworthy that these profiles vanish for different time intervals, with disappearance occurring for $t<0$ and $t>0$ when the parameter $R_2$ takes positive and negative values, respectively. For the localized or sech-type modulated nonlinearity, the density profiles of positons become compressed and curved within the condensate density background. In the case of periodic modulated nonlinearity, positons exhibit a periodic nature, and we have observed an increase in periodicity as the nonlinearity strengths are adjusted. Our findings contribute to a deeper understanding of the behavior of matter wave positons in BECs under different types of modulated nonlinearities. These results shed light on the intricate interplay between nonlinearity, external trapping potentials, and the corresponding effects on the density profiles of positons. The theoretical findings presented in this study, along with previous research in the literature, offer a valuable groundwork for experimental researchers to explore and validate the deformation of solitons/positons in $\mathcal{PT}$-symmetric systems with spatiotemporal modulation. These investigations can be extended to various fields, such as Bose-Einstein condensates and nonlinear optics, that are currently of great interest. Additionally, as a future direction, this theoretical study can be readily expanded to examine higher-order solitons, breathers, and rogue waves. It can also encompass the exploration of combined spatial and longitudinally varying trap potentials, nonlinear effects, and novel forms of $\mathcal{PT}$-symmetric potentials, potentially leading to the discovery of new applications.