Controlling Matter-Wave Smooth Positons in Bose–Einstein Condensates

: In this investigation, we explore the existence and intriguing features of matter-wave smooth positons in a non-autonomous one-dimensional Bose–Einstein condensate (BEC) system with attractive interatomic interactions. We focus on the Gross–Pitaevskii (GP) equation/nonlinear Schrödinger-type equation with time-modulated nonlinearity and trap potential, which govern nonlinear wave propagation in the BEC. Our approach involves constructing second-and third-order matter-wave smooth positons using a similarity transformation technique. We also identify the constraints on the time-modulated system parameters that give rise to these nonlinear localized proﬁles. This study considers 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 strengths, we observe a rich variety of captivating behaviors in the matter-wave smooth positon proﬁles. These behaviors include stretching, curving, oscillating, breathing, collapsing, ampliﬁcation, and suppression. Our comprehensive studies shed light on the intricate density proﬁle of matter-wave smooth positons in BECs, providing valuable insights into their controllable behavior and characteristics in the presence of time-modulated nonlinearity and trap potential effects.


Introduction
Theoretical investigations into the nonlinear collective excitations of matter waves have emerged as a highly intriguing and pertinent field, especially in light of the 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 inspire particular interest. The concept of a soliton 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][7][8][9]. 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,6,[10][11][12]. On the other hand, BECs emit Faraday and resonant density waves when subjected to harmonic driving [13]. The characteristics of density waves in dipolar condensates at a temperature of absolute zero have been investigated using both mean-field variational and full numerical approaches [13]. 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 nonlinear atom interactions using the Feshbach resonance (FR) technique as well as by varying the external trap potential [14][15][16][17]. 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-onedimensional BECs [17][18][19][20][21][22][23][24]. Furthermore, studies focusing on the variable coefficient NLS equation have unveiled the possibility of manipulating and enhancing these localized profiles through the utilization of inhomogeneity parameters [5][6][7][25][26][27][28][29][30][31]. Moreover, investigations have explored the identification of soliton solutions in the generalized NLS equation with spatially modulated parameters and an external potential [32][33][34][35]. Notably, in [32], the authors identified alternate solitons and determined their stability regions in both two-dimensional and one-dimensional BEC models. These investigations were conducted in the presence of optical lattices created by laser beams illuminating the condensates. The analysis of soliton propagation in optical and condensed matter systems with parity-time (P T ) symmetry, particularly in inhomogeneous setups, has received significant attention. Accordingly, considerable efforts have been devoted to showcasing the existence of stable bright solitons, dark solitons, and vortices within the NLS equation featuring P T -symmetric potentials [36][37][38]. 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 [39]. 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 smooth positons within the GP equation. We further endeavor to explore the effects of the time-dependent modulation of nonlinearity parameters on the characteristics of smooth positon profiles.
Positons, unlike exponentially decaying soliton solutions, are weakly localized nonlinear waves that hold significant importance in the field of nonlinear physics [40][41][42][43]. 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 a soliton remains unchanged following a collision with a positon. However, the positon is influenced by the carrier wave and envelope, resulting in a finite phase shift [44,45]. Notably, Matveev's positon solution to the Korteweg-de Vries (KdV) equation exhibited a spectral singularity [42]. Building on this pioneering work, positon solutions have been successfully derived for other nonlinear evolution equations; one may refer to [46][47][48] for more information. Recent efforts by Cen et al. introduced the concept of smooth positons or degenerate soliton solutions by allowing for the spectral parameter to take complex values [49,50], thereby eliminating the singularity in the KdV equation. Following these advancements, endeavors have been made to construct smooth positon solutions for various nonlinear evolution equations, including the focusing mKdV equation [51], the complex mKdV equation [52], the derivative NLS equation [53,54], the NLS-Maxwell-Bloch equation [55], the higher-order Chen-Lee-Liu equation [56], and the Gerdjikov-Ivanov equation [57]. 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 [58,59]. Inspired by these advancements in the field of positons, our research aims to construct smooth positon solutions within the GP equation, incorporating time-varying nonlinearity and trap potentials.
The crucial step in this attempt involves utilizing a similarity transformation [5][6][7][8]20] on a meticulously chosen ansatz solution. This transformation effectively converts the GP equation with time-varying coefficients 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 smooth positon solutions in the altered coordinate system [6,8,20,23]. By leveraging the combination of known positon solutions of typical NLS equations 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 potentials), 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 outcomes by judiciously selecting appropriate nonlinearity strengths and trap potentials [20,25,26].
Motivated by the experimental feasibility of studying BECs, our research focuses on exploring the characteristics of matter-wave smooth 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 constant coefficient NLS (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 ccNLS equation, we present matter-wave smooth positon solutions of the GP equation. We investigate the controllable behavior of positon density profiles with respect to three different forms of variable nonlinearity parameters, namely, (i) kink-like nonlinearity R(t) = R 0 + R 1 tanh (R 2 t + R 3 ), (ii) localized or sech-type nonlinearity R(t) = R 0 + R 1 sech (R 2 t + R 3 ), and (iii) periodic nonlinearity R(t) = R 0 + R 1 sin (R 2 t + R 3 ), for which 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 matter-wave positon density profiles of a second-and third-order experience stretching, while their amplitudes can either be enhanced or suppressed. It is important to note that these profiles vanish at 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 have structured our work as follows. In Section 2, we account for the GP equation with time-modulated nonlinearity and trap potentials. The second and third-order smooth positon solutions are deduced for this equation using a similarity transformation. The integrability requirement between the modulated nonlinearity and trap potential is obtained while applying the integrable technique to the considered equation. In Section 3, by suitably choosing the different forms of the variable nonlinearity parameter, we explore the various characteristics in the density of matter-wave smooth positon profiles. Finally, in Section 4, we provide conclusions based on our observations.

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, which is 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,8,17,20] i ∂ψ ∂t where ψ(x, t) is the condensate wave function and t and x represent dimensionless temporal (propagation direction) and spatial (transverse) coordinates, respectively. The atom-atom interaction term (i.e., the scattering length between atoms) denoted by the representation, R(t), can be tuned using the FR technique. In a series of exquisite experiments using sodium and rubidium condensates, FRs were investigated [14,15]. They have also been employed in a variety of significant experimental studies, such as the creation of bright and dark matterwave solitons, among others. The time-modulated trap potential is denoted by λ 2 (t). In the study of trapped BECs, the trap potential along the elongated axis has been intentionally selected to vary with time, t, in order to investigate the characteristics of the BECs within the trap. As a result, both the coefficient of nonlinearity (R) and the potential parameter (λ) can exhibit time dependence. By appropriately selecting 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, which can be achieved through the adjustment of external magnetic fields and optically controlled interactions using techniques such as the FR method [16,20,23]. It is worth noting that such an analysis can be extended to study the influence of spacedependent parameters on solitons/positons in the system under consideration (1) or in the context of the generalized NLS equation with the presence of optical lattices. A noteworthy and unique characteristic of solitons in the presence of nonlinear lattices is the existence of a finite threshold value of the soliton norm, which is necessary for the lattices' formation [32]. This property does not occur in the absence of nonlinear lattices. In nonlinear lattices, solitons are created directly as opposed to bifurcating from Bloch modes as seen in linear lattices. This distinction highlights the unique nature and behavior of solitons in the presence of spatially modulated nonlinearity [32,33]. As a consequence of these observations, substantial progress have been made with regard to investigating the localized solutions in the generalized NLS equation with spatially modulated parameters and the external potential [32][33][34][35]. This opens up promising avenues for future investigations with respect to understanding the behavior and dynamics of solitons and positons in spatially modulated systems.
In this work, to study the matter-wave smooth positons in (1), we adopt the similarity transformation mentioned below to map the time-modulated GP Equation (1) to the ccNLS equation [5][6][7][8]20,38]: In Equation (2), the unknown functions, namely, s(t), η(x, t), τ(t), and θ(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: The explicit expressions of the unknown functions can be acquired by solving the aforementioned set of equations, and they take the following form 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 in [8,20,38].
Using Equation (5), one can find the λ(t) by fixing the R(t), and vice versa. In this work, we consider the physically intriguing function R(t) and determine the λ(t) via the following expression: In Equation (2), the function φ(η, τ) was found to fulfill the ccNLS equation: i ∂φ ∂τ The equation under consideration (6) exhibits a wide range of localized solutions, including solitons, breathers, rogue waves, and their corresponding profiles [6,20,43]. In this study, we focus on the smooth 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(t) while ensuring the satisfaction of condition (5), we can derive matter-wave positon solutions for the GP Equation (1) in the following form: where φ(η, τ) is the smooth positon solution of the ccNLS Equation (6). The solution (7) has the potential to generate a multitude of novel positon structures that can 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 of and potential new perspective regarding the similarity transformation is worth emphasizing, as it allows for 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 Section 3, we construct smooth positon solutions for the ccNLS Equation (6) to analyze matter-wave positons in quasi-one-dimensional BECs.

Characteristics of Matter-Wave Smooth Positons in BECs
In this section, we investigate the characteristics of matter-wave positons using solution (7) by considering three different forms of the time-varying nonlinearity parameter and the associated trap potentials. For this investigation, in the following, we first derive the secondorder and third-order smooth positon solution for the ccNLS Equation (6).

Second-Order Matter-Wave Smooth Positons
Now, we explore the various novel density profiles of second-order matter-wave smooth positons in BECs. For this objective, we first derive the second-order smooth positon solution for the ccNLS Equation (6) as provided by where where α 1 is the eigenvalue of the spectral parameter, α * 1 is the complex conjugate of α 1 , and η and τ are provided in Equation (4). The above-mentioned second-order smooth positon solution (8) can be obtained from the material presented in [58] by considering ν = 0. Substituting this solution into (7) along with the suitable form of R(t), we obtain the second-order matter-wave smooth positon solution of (1). Now, utilizing the above mentioned solution, we move to investigate its controllable behaviors through three different forms of a time-modulated nonlinearity parameter, namely, (i) kink-like nonlinearity R(t) = R 0 + R 1 tanh (R 2 t + R 3 ), (ii) localized or sech-type nonlinearity R(t) = R 0 + R 1 sech (R 2 t + R 3 ), and (iii) periodic nonlineaity R(t) = R 0 + R 1 sin (R 2 t + R 3 ), where R 0 , R 1 , R 2 , and R 3 are arbitrary parameters. Moreover, it is important to mention that this investigation can be extended to study the impact of spatially modulated coefficients [32][33][34][35] on the matter-wave smooth positons of (1), which will provide fruitful directions for future investigations. In the following, we provide a comprehensive demonstration of the impact of time-modulated nonlinearity parameters on the positon density profiles.
To begin, we employ the kink-like nonlinearity parameter, that is, R(t) = R 0 + R 1 tanh (R 2 t + R 3 ), to reveal the novel features in BECs. Substituting this nonlinearity term into the generalized solution (7) yields where φ(η, τ) is the second-order smooth positon solution of the ccNLS Equation provided in Equation (8). Utilizing this solution (10), we conduct a thorough analysis of the positon density profiles, investigating their diverse characteristics as we vary the strength of the time-modulated nonlinearity parameter.
Figure 1a-f illustrate the qualitative profiles of a second-order matter-wave smooth positon in BECs corresponding to a kink-like nonlinearity modulated parameter R(t) = R 0 + R 1 tanh (R 2 t + R 3 ). By selecting specific parameter values, i.e., 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 α 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, when modifying the eigenvalue to α 1 = 0.3 + 0.6i, Figure 1b clearly demonstrates a shift in the orientation of the positon profile accompanied by an enhancement of its amplitude. Furthermore, we investigate the effects of adjusting the strengths of the nonlinearity parameters, specifically R 0 , R 1 , and R 2 . For example, 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 assumes 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 localized-type or sech-type nonlinearity, namely, R(t) = R 0 + R 1 sech (R 2 t + R 3 ), in order to investigate the distortion of positon profiles in the condensate density background. By inserting this form of R(t) into (7), we obtain the matter positon solution in the following form where φ(η, τ) is the positon solution of the ccNLS Equation, which is provided in Equation (8). Figure 2a-f represent the density profiles of a second-order matter-wave smooth positon with a modulated nonlinearity function yielded by R(t) = R 0 + R 1 sech (R 2 t + R 3 ). In this study, we consider specific parameter values, i.e., 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 assumes a negative value, we observe the formation of two peaks at the center (t = 0), which is 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, the further stretching of the positon, and the formation of a curved profile, as shown in Figure 2e. At R 0 = 2.5 and R 2 = 1.5, a compressed density profile with an identical amplitude is obtained, as depicted in Figure 2f.
(e) (f) Figure 1. Density profile of second-order matter-wave smooth positons for (1) with the kink-like nonlinearity modulated function R(t) = R 0 + R 1 tanh (R 2 t + R 3 ). The parameters are (a) R 0 = 1.05, In (b-f), the mentioned parameters are varied, while the remaining parameters are maintained as they are in (a). The other parameters are α 1 Finally, we investigate the influence of a periodically modulated nonlinearity given by R(t) = R 0 + R 1 sin (R 2 t + R 3 ) on the positon density profiles. By inserting this R(t) into (7), we obtain the explicit form of a matter-wave positon solution yielded in the following expression: where φ(η, τ) is the positon solution of the ccNLS Equation given in (8).
The second-order matter-wave smooth positon density profiles are displayed in Figure 3 for a periodic nonlinearity-modulated function R(t) = 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 can 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 can 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 Figure 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 5.25 are chosen, the oscillation in the positon profile increases and retains a similar orientation, as evident in Figure 3e,f, respectively.
(e) (f) Figure 3. Density profile of second-order matter-wave smooth positons for (1) with the periodic nonlinearity modulated function R(t) = R 0 + R 1 sin (R 2 t + R 3 ). The parameters are (a) R 0 = 1.5, In (b-f), the mentioned parameters are varied, while the remaining parameters are maintained as they are in (a). The other parameters are the same as in Figure 1.

Third-Order Matter-Wave Smooth Positons
In the previous sub-section, we investigated the modifications of second-order matterwave smooth positon profiles by varying the distributed coefficients (the nonlinearity function) in the GP Equation (1). In this sub-section, we delve into the deformations of third-order matter-wave smooth 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 following form: where P 2 =e −i(2τ(α 2 1 +α * 2 1 )+η(α1+α * 1 )) P 21 e 2i(α 1 −α * 1 )(η+2τ(α 1 +α * 1 )) with where It is important to note that the third-order smooth positon solution can be deduced from the research in [58] by setting ν = 0. By substituting this third-order smooth positon solution of the ccNLS Equation (13), along with the appropriate form of the modulated parameter R(t), into Equation (7), we delve into an analysis of the underlying characteristics of the GP Equation (1).
Figure 4a-f present the condensate profiles of third-order smooth positons with a kink-like nonlinearity-modulated function R(t) = R 0 + R 1 tanh (R 2 t + R 3 ). By considering the 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 changing the value of R 0 to 1.75. On the other hand, in Figure 4d, when the parameter R 1 is increased to 1.25, the positon tends to disappear within the condensate profile when 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 Figure 5a-f, we show the density profiles of third-order matter-wave smooth positons for the sech-like nonlinearity-modulated function R(t) = R 0 + R 1 sech (R 2 t + R 3 ) for the GP model (1). Employing the following initial parameters, namely, R 0 = 1.25, R 1 = 0.01, R 2 = 1.05, and 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 crests, 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 R 2 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 profile becomes compressed, which is accompanied by a decrease in its amplitude, as shown in Figure 5f.
(e) (f) Figure 5. Density profile of third-order matter-wave smooth positons for (1) with nonlinearitymodulated function R(t) = R 0 + R 1 sech (R 2 t + R 3 ). The parameters are (a) R 0 = 1.25, R 1 = 0.01, In (b-f), the mentioned parameters are varied, while the remaining parameters are maintained as they are in (a). The other parameters are the same as in Figure 4.
Finally, we investigate the density profiles of the third-order matter-wave smooth positon in one-component BECs utilizing a periodic nonlinearity-modulated function, denoted as R(t) = 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. Notably, in the results, 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 there is a variation in the oscillation behavior, which is clearly depicted in Figure 6f.
(e) (f) Figure 6. Density profile of third-order matter-wave smooth positons for (1) with the nonlinearitymodulated function R(t) = R 0 + R 1 sin (R 2 t + R 3 ). The parameters are (a) R 0 = 1.5, R 1 = 0.05, In (b-f), the mentioned parameters are varied, while the remaining parameters are maintained as they are in (a). The other parameters are the same as in Figure 4

Conclusions
In our study, we have derived the second-and third-order matter-wave smooth positon solutions of the GP equation. These solutions capture the features 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 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 matter-wave smooth 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 localized or sech-type modulated nonlinearity, the density profiles of the positons become compressed and curved within the condensate density background. In the case of periodically modulated nonlinearity, the positons exhibit a periodic nature, and we observed an increase in periodicity as the nonlinearity strengths were 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 those presented in previous research, offer a valuable groundwork for experimental researchers to explore and validate the deformation of solitons/positons in P T -symmetric systems with spatiotemporal modulation. These investigations can be extended to various fields, such as BECs and nonlinear optics, which 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 P T -symmetric potentials, potentially leading to the discovery of new applications.