Dynamics and uncertainty for maximally entangled bipartite system constrained on a helicoid

The classical and quantum dynamics of particles constrained on a right helicoid is discussed via Dirac approach. We show how the uncertainty in measurement for observables of maximally entangled system is affected by the total number of constrained particles σ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma $$\end{document}, external magnetic field B→\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\vec {\mathcal {B}}$$\end{document}; as well as by geometric parameters like the pitchρ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rho $$\end{document} and the radial position of particles. In doing so we also highlight numeric bounds on the external field strengths to tune both intraparticle and bipartite entanglement and we remark that the bipartite entanglement is more robust to changes in the fields than the intraparticle entanglement, in this framework. We also highlight specific parameter regimes which lead the uncertainty (in measurement) to achieve respective parameter independence and, for a particular subset of commutation relations, the system remains confined in the quantum regime even in the limit σ→∞\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma \rightarrow \infty $$\end{document}. It is observed that the uncertainties are strongly influenced by the geometric parameters e.g., ρ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rho $$\end{document}, and the strength of bipartite as well as intraparticle entanglement might be controllable through ρ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rho $$\end{document}. The energy equation for this setup is obtained and the additional terms are discussed which arise due to quantum correlations, orbit–orbit interaction and the normal Zeeman effect, which leads to the splitting of the energy level into 11 non-degenerate levels. Finally we comment that, a linkage of this phenomenology with Aharonov–Bohm like effect might be possible by strictly confining B→\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\vec {\mathcal {B}}$$\end{document} along the central axis of the helicoid.


Introduction
From a classical physics perspective it is possible to specify all the degrees of freedom of a system simultaneously with arbitrary precision, but uncertainty relation is fundamental in quantum mechanics and is one of the chief features that underlines a difference between the classical and quantum physics. In the seminal EPR paper it was shown that if two systems are correlated then the uncertainty of measurement between a conjugate pair of observables of these systems reduces [1]. The extended uncertainty relation (by Robertson) sets a limit on the precision of simultaneous measurement of any pair of observables in terms of a lower bound, which gives a quantitative measure of the uncertainty of measurement [2].
The uncertainty relations in more recent literature are usually expressed in the well known entropic form e.g., see Ref. [3], in which the lower bounds quantifying uncertainty of measurement reduce if the two parties are maximally entangled, implying that quantum entanglement has a significant effect on the product of dispersion. In [4] it was reported that the lower bound deviates from the traditional Heisenberg Uncertainty Relation (HUR) in such a way that the precision of simultaneous measurement increases for maximally entangled bipartite and tripartite system.
One can also consider various types of entanglement where quantum correlations can exist between distinct degrees of freedom of a single particle i.e., intraparticle entanglement [5], or two particles (bipartite entanglement), as well as various background geometries.
Dirac approach for constrained Hamiltonian dynamics provides a modified symplectic structure to obtain the constrained dynamics of the system [6]. To study the quantum dynamics, one must work in the Hamiltonian formalism. However, for singular Lagrangians i.e., one can not express allẋ α in terms of p α and the Hamiltonian dynamics can not be obtained by the usual procedures. Dirac proposed a generic technique to handle singular Lagrangians in which the constrained dynamics can be obtained by Dirac brackets instead of Poisson brackets [6,7]. Dirac brackets can be used to assign geometric meaning to Hamiltonian formulation in the presence of second class constraints. A similar study in this direction was done in Ref. [8], and a different approach to find geometric duality for quantum equations in the pilot-wave limit was presented in Refs. [9][10][11]. Studies on Dirac approach for particle constrained on a helicoid and for a free particle on S 3 were carried out in [12,13]. Dirac quantization on curved spaces leads to extra terms in the energy equation which can be linked with the curvature of space or shift in the energy spectrum leading to novel physics. It is therefore also of interest to discuss such energy terms for constrained quantum correlated particles on curved spaces [14], which was done in [15] for a maximally entangled system constrained on S 1 . In view of the usefulness of generic constrained quantum systems for modern avenues in quantum research, e.g., quantum matter, it is instructive to explore some fundamental aspects of such systems e.g., the effect of entanglement on product of dispersion for quantum systems constrained on certain geometries. We address this question in this paper for a maximally entangled bipartite system on a right helicoid, which is a minimal surface. We will develop uncertainty relations and study the effect of various parameters on the product of dispersion. For this, we choose a pure bipartite state [16,17].
and compute the lower bounds using Robertson relation to quantify the uncertainty of measurement for various observables acting on this bipartite state. Note that throughout this paper, the statements Uncertainty of (simultaneous) measurement/product of dispersion/precision of simultaneous measurement will be used interchangeably, as appropriated by context. The paper is organized as: In Sect. 2, we present the constrained classical dynamics on a right helicoid for σ = 1 and σ = 2. In Sect. 3, we quantize the classical system via Dirac quantization rule for particles constituting a maximally entangled bipartite system. In Sect. 4, we discuss the results of the effect of σ , A, ρ and radial position on the uncertainty of measurement. In Sect. 5, we find the energy equation and address the significance of additional terms arising in the equation. Section 6 is the summary of results.
2 Constrained classical dynamics on a right helicoid 2.1 Single particle case (σ = 1) We consider a free particle (σ = 1) constrained on a right helicoid. The right helicoid is parameterized by the following equation [18] where η is the radial distance from z-axis (hereafter called x c -axis); ϑ and ρ is the angle of the twist and pitch of the helicoid, respectively. The Lagrangian of this system can be written as Since the particle is constrained on the helicoid, this constraint is given by [12] f In order to get the full dynamics of the system, we must incorporate this constraint in Eq. (5) into the Lagrangian in Eq. (4).
where ω is the Lagrange multiplier which serves as a degree of freedom. We have assumed the particle to be constrained on helicoid under the effect of external magnetic field and no additional force acts on the particle (i.e., V = 0). The Hamiltonian of the system is As the Hamiltonian in Eq. (7) is not independent of velocities, therefore we need to follow the Dirac's approach to obtain the dynamics of the system. In such cases one needs to consider the relations which prevent the Hamiltonian to be independent ofη. These relations which emerge directly from the Lagrangian are the primary constraints ζ λ ≈ 0.
One can obtain a series of second class constraints from such relations via the consistency condition [6]: where H T is the total Hamiltonian which is obtained by multiplying a function of time (κ), with the primary constraint and including this term into the Hamiltonian Eq. (7).
whereω is absorbed into κ. The constraint emerging directly from the Lagrangian is given by p ω = ∂L ∂ω = 0. We will call this the primary constraint and denote it as ζ 1 = p ω ≈ 0. The secondary constraints are obtained by means of the consistency conditionζ = [ζ, H T ] ≈ 0. These constraints are listed below, In such cases, the dynamics of the system can not be obtained by the usual Poisson brackets but by Dirac brackets. For any two physical quantities U and V the Dirac bracket is given by where, μν = ( −1 ) μν is the inverse matrix of second class constraints with μν = [ζ μ , ζ ν ] [6]. In case of singular Lagrangians, we need to include the constraints into the Lagrangian, Hamiltonian as well as the brackets to obtain the correct constrained dynamics of the system. The inverse matrix of second class constraints is given by Using Dirac's approach, the classical dynamics is found to be The corresponding quantum dynamics of a free particle constrained on a helicoid can be obtained by the Dirac quantization rule: 2.2 Two particle case (σ = 2) We consider two particle case under the influence of scalar and vector potential. The Lagrangian is given by where α is the particle index such that α = 1, . . . , σ anḋ The particles are assumed to keep their motion on the surface of the helicoid under the action of a uniform magnetic field B directed along x c -axis. Applying this constraint on the motion of these particles leads to The Hamiltonian is: Using V = 0, and also This is the primary constraint since it is derived directly from the Lagrangian. The total Hamiltonian is thus given by: whereω is absorbed in κ.
If a constraint is initially zero, it must be zero for all times.
The resulting equation is taken as another constraint. Following the consistency condition, one obtains a series of constraints as follows: Using ζ 1 = ζ 2 = 0, the total Hamilton (24) is To find the Dirac brackets for this system, we follow (13), where ζ represents the constraint evaluated above, and μν is a matrix such that μν = ( −1 ) μν and μν = [ζ μ , ζ ν ].
The inverse matrix is computed and is given below with, The classical dynamics obtained via the Dirac's approach is given in the Appendix. Following Dirac's approach we achieved four sets of equations (listed in Appendix) which we have labelled as Type-I, Type-II, Type-III and Type-IV equations.
Type-I This set of equations reflects how the particle's motion deviates from a linear trajectory along a particular axis. According to classical mechanics, the trajectory of a particle in a particular direction is encapsulated in the Poisson brackets. Type-I Dirac brackets differ from Poisson brackets due to the 47) in Appendix for details) originating from the constraints. The particle's dynamics is deviated from their original axis by these additional factors which include constituents (ρ, x a , x b ) underpinning the fact that underlying geometry is a Helicoid. Type-II All canonical brackets of type-II vanish in classical mechanics; on the contrary, if a particle is moving along x a axis then p b can generate its translation along x b axis due to the constraints incorporated into the dynamics of the system. Therefore, these brackets are non-zero in this framework.
Type-III This is the set of equations where Dirac brackets coincide with the standard canonical brackets of classical mechanics. This is because none of the ζ μ is function of p only, thus these brackets are not affected by the constraints. Type-IV These equations contain Dirac brackets between various momentum components of both particles. These equations contain 2 and ρ define the radius and pitch of the helicoid. Thus, these brackets specify the trajectory of particles on the surface of right helicoid.

Quantum dynamics
We consider bipartitioning the system in such a way that τ (σ ) = τ (α) + τ (β) (i.e., the system with σ number of particles is partitioned into two subsystems containing α and β number of particles) The bipartitioning of the classical system discussed above treats two particles as two parties (A and B) which constitute a global system, whose Hilbert space can be written as The basis sets spanning H A and H B are taken to be {ψ α } σ α=1 and {φ β } σ β=1 respectively. We assume these subsystems to be maximally entangled having unit concurrence. The maximally entangled pure state with respect to this partitioning is given by an equation of the following form where C = 2 | εδ − ξυ |= 1 for a maximally entangled state [19]. The state of a single free (unentangled) particle on helicoid is defined as: (η, ϑ) = exp(ilϑ)ψ(η), where l is angular momentum [20]. The maximally entangled pure bipartite state for two particles constrained on a helicoid is given by an equation of the following form [15], To quantize the Dirac brackets we use Dirac quantization rule. The commutators in Eq. (67) contain x c -component of angular momentum. The operator is defined as [7]: It is important to note that the operators acting on an entangled state lead to a reduced uncertainty of measurement in view of the single free (unentangled) particle discussed in Sect. 2.1. This deviation in uncertainty leads to an enhanced precision in simultaneous measurement of observables. The Robertson relation serves as a measure of this deviation by computing the numerical lower bounds for a pure bipartite state of the form given in [2] and in Eq. (37). The lower bounds are computed under the assumption that ψ(η) is a constant function. These lower bounds are affected by a series of parameters including: the number of particles (σ ), external fields, pitch of the helicoid (ρ) and the radial position of particles on the helicoid. The lower bounds are computed numerically using Robertson relation and a few of these are given in the Appendix. These lower bounds will be plotted as a function of some parameters in the following section.

Uncertainty of measurement as a function of σ
The plots in Figs. 1 and 2 present the product of dispersion as a function of number of particles constrained on the right helicoid. In Figs. 1 and 2, the lower bounds of Type-I and Type-IV equations are evaluated at x a , ρ = 10, e = 1 and the angular position (in radians) of the particles (constrained on the helicoid) is: Figure 1 represents uncertainty of measurement of Eqs. (45)-(47) as a function of σ . These plots essentially tell us how the lower bounds of the subsystem of two particles are affected by varying the total number of particles in the global system. The precision of simultaneous measurement of position and momentum decreases for both particles but this trend is not as pronounced for x c and x a component of particle 2 as in the case of x b component. However, for particle 1, x a and x b manifest a sharp deviation with increasing σ . The lower bounds of Fig. 1 become independent of σ as σ → ∞, i.e., for a very large number of particles, the system achieves σ -independence. Moreover, the uncertainty increases with increasing σ (for observables of the same particle) and never touches zero even in the limit σ → ∞, implying that the system remains in the quantum regime for an arbitrarily large number of identical constrained entangled particles, a similar conclusion was given for a different setup in [21]. Figure 2 illustrates that the uncertainty decreases exponentially for all commutation relations implying that the system tends to exhibit a classical-like behaviour for σ → ∞. It is evident from Figs. 1 and 2 that the product of dispersion between the observables of the same particle increases and that between the observables of different particles decreases in view of the single particle case. This indicates that the entanglement between the two parties increases (in accordance with Ref. [3] which suggests that increased entanglement leads to smaller uncertainty), however no intraparticle entanglement exists between position and momentum observables of either particle.

Uncertainty of measurement as a function of external field
The left and right panel plots in Fig. 3 present the lower bounds of Eqs. (67) and (73) as a function of A, respectively.
In the left panel plot, red (squares) and green (circles) curves represent the lower bound of Eq. (67) for particle 1 and 2, respectively, and in the right panel, red and green curves represent the lower bound of Eq. (73) i.e., [ p a 1 , p b 2 ] and [ p a 2 , p b 1 ], respectively. The (dashed) black line in both plots corresponds to the single free particle case discussed in Sect. 2.1. The left plot depicts that the lower bounds quantifying uncertainty of simultaneous measurement between various momentum degrees of freedom of each particle (i.e., [ p a 1 , p b 1 ] and [ p a 2 , p b 2 ]) initially decrease as a function of A and hit minimum at A 1 = −0.4 (for subsystem 1) and A 2 = 2.3 (for subsystem 2) which, for our choice of parameters correspond to the maximum intraparticle entanglement that can be achieved for subsystem 1 and 2, respectively. The inset reveals that in view of the single particle case the intraparticle entanglement for both subsystems exists for lower bounds below 2 × 10 −3 ; the intraparticle entanglement is zero beyond this cutoff. Region-I (−0.8 ≤ A ≤ 0) and region-II (−0.2 ≤ A ≤ 4.7) identified in the inset of the left plot correspond to the entanglement regions for subsystem 1 and subsystem 2, respectively. The range of vector potential that favors intraparticle entanglement in the total system (subsystem 1 and 2) is thus −0.8 ≤ A ≤ 4.7.
Lower bounds quantifying the uncertainty of measurement between momentum degrees of freedom of distinct particles, as presented in the right panel plot of Fig. 3, also illustrates an initially decreasing trend for negative values of A. These bounds are zero for red (squares) and green (circles) curves at A 12 = 0 and A 21 = 0, respectively. The bipartite entanglement exists (for both red and green curves) for It is important to note, from both plots in Fig. 3, that at A = 0, the intraparticle entanglement for subsystem 1 is zero whereas the bipartite entanglement is maximum which suggests that, for a physical setup like ours, the intraparticle entanglement can only be achieved in the presence of external field B. However, the bipartite entanglement can be achieved and appropriately tuned (in the absence of exter-nal magnetic field) for maximum strength in regions of the parameter space respecting the corresponding cutoff values, as discussed above.
It is possible to calculate the minimum/maximum strength of the external field B that enforces intraparticle and bipartite entanglement in the system, from vector potential. Since the external field is assumed to be uniform and directed along x caxis, the vector potential corresponding to this field rotates about x c -axis and is: A = nB /2, where n is the number of twists of the helicoid. We choose a right helicoid with n = 2, the minimum and maximum strength of external magnetic field required for intraparticle entanglement is B i l = 0 and B i u = 0.47, respectively. The minimum/maximum strength of the field required for bipartite entanglement is B bi l = 0.12 and B bi u = 0.51, respectively. Thus it suggests that the bipartite entanglement is more robust against small fluctuations in B as compared to the intraparticle entanglement.  Fig. 4 show the uncertainty as a function of ρ for particle 1 and 2, respectively. The red (solid), green (circles) and magenta (dashed) curves correspond to the lower bounds of [x a , p a ], [x b , p b ] and [x c , p c ], respectively. The left plot in Fig. 4 infers that the uncertainty of measurement in the case of [x a , p a ], [x b , p b ] (i.e., red and green curves) drops sharply for ρ ≤ 10 and stabilizes for higher ρ > 10; a phenomenon which is reversed for the case of [x c , p c ] (magenta) curve where it rises sharply for small ρ regime. The plot in the right panel of Fig. 4 also portrays a similar scenario but with a much smaller lower bound for green (circles) ([x b , p b ]) that tend to zero for ρ > 10 approximately, and a higher bound for red (solid) curve ([x a , p a ]). Figure 4 illustrates that uncertainty of measurement is not stable against variations in the pitch for ρ < 10; with the param-eter values utilized. It however achieves robustness in parameter regimes with large ρ. Thus, the lower bounds (and hence uncertainty of measurement) can achieve ρ-independence for particles constrained on helicoidal geometries with large ρ. The left plot in Fig. 5 depicts that the lower bounds increase with increasing ρ, and reach a somewhat stable point for regimes in parameter ranges ρ > 10. Besides, there is no intraparticle entanglement between momentum degrees of freedom (of each particle) in the region 1 ≤ ρ ≤ 3. For ρ > 3, intraparticle entanglement between momentum degrees of freedom for subsystem 2 can be achieved. The inset in right panel plot tells that initially the bounds increase and then decrease as a function of ρ, for both [ p a 1 , p b 2 ] and [ p a 2 , p b 1 ]. It is interesting to note (from right panel plot) that both cases lie below the black (dashed) curve which corresponds to single particle bounds. Thus, one can conclude from Fig. 5 that bipartite entanglement between the , only the plot is flipped by π 2 for observables depicting the location of particles at right angles to x a -axis. The parameter regions in bright yellow colour in Fig. 6 correspond to relatively higher uncertainty ranges, and it is easy to see that the top left and middle panel plots hint at an increased uncertainty away from the origin for both cases. Increasing and decreasing (x a σ , x b σ ) imply that the particles move towards the edges or towards the central axis (x c -axis) of the helicoid, respectively. The top right plot for the x c -component represents the opposite feature with the uncertainty increasing with decrease in (x a σ , x b σ ) and attaining the maximum value of 0.25 at the origin (0, 0). The uncertainty of measurement increases with the decrease in the pitch ρ and it is evident from the left and middle panel plots in the bottom row of Fig. 6 that the uncertainty of measurement increases for significantly large parameter regimes. Although, the lower bounds reduce with increase in ρ, no reduction in uncertainty is obtained in view of the singe particle case and no intraparticle entanglement exists between position and momentum observables. Fig. 7 give the lower bounds of Eq. (73), i.e., ( p a 1 ) 2 ( p b 2 ) 2 and ( p a 2 ) 2 ( p b 1 ) 2 respectively, as a function of (x a 1 , x a 2 ). The left plot shows a small faint peak and one relatively large and broad peak leading to the global maxima value 0.093403. In the right plot there are two peaks located at (−6.28, −0.01) and (6.28, −0.01). Figure 7 also exhibits the sensitive dependence of uncertainty on the radial coordinates; the uncertainty decreases quickly as one deviates from the peaks.

Energy equation
The Dirac quantization leads to some additional terms in the energy equation which can be of great physical significance, in connection to curved spaces [7,14,15] and to the external fields. In this section, we obtain the energy equation for our setup and study the effect of various parameters i.e., B, A etc., on these terms. By utilizing Eq. (38), the total Hamiltonian in Eq. (29) can be written as 2 is the radius of right helicoid (assuming the full helicoid to perfectly fit inside a cylinder of radius ≈ 2 ). The solution of the Schrödinger equation for maximally entangled bipartite state in Eq. (37) leads to the following energy equation, where,˜ = − 1,h = 1 and Q 1 aψ = e il(ϑ 1 +ϑ 2 +ϑ 1 ϑ 2 ) indicates the variation in state and QE represents the additional energy terms associated with Q . Keeping in view the definition of Q, the indices in the above equation can be interpreted as: the index a, b, c, n in the subscript refers to the derivative of ψ(φ) with respect to x a , x b , x c , η respectively, and the subscript (ab) refers to the second order derivative with respect to x a and x b . The superscript represents the particle number and the index ψ(φ) in the subscript refers to the derivative of ψ(φ). The first term in the energy equation has the form E = L 2 2mr 2 , with r → ( 2 α + ρ 2 ); 2 α and ρ being the radius and pitch of the helicoid respectively. Hence this term characterizes the energy of particles (constrained to retain their motion on a right helicoid) as a function of angular position ϑ α . One can see that: E = E(1, 2) = E 1 + E 2 i.e., this is the case of interacting particles, where it is not possible to specify each particle independently of the other (particle(s)) making up the system. The particles are maximally entangled and the energy of each particle is affected by the simultaneous angular position ϑ α of the other particle. It is also to be noted that the energy equation contains additional terms in view of the equation for particles under the effect of a vector potential (but free from any additional constraints on its trajectory), as given by the following equation The energy equation for free particle f (x a , x b , x c ) = V = 0 andh = 1 takes the following form The additional energy term (second term) in Eq. (40) has non-constant terms, which involve derivatives of the function characterizing the radial distance of particles from the central axis of the helicoid. Also note that, since we choose the particles with l = 5, the interaction between the external field B and the field generated by the orbital motion of the particles splits the energy level into 11 different energy level (i.e., normal Zeeman effect leads to 11 split levels with slightly different energies). In the absence of spin-spin interaction, the total angular momentum vector (obtained as a result of orbit-orbit coupling of the two particles) precess about the magnetic field and the projection of orbital angular momentum vector of these particles along x c -axis is given by the following equation The total x c -component of angular momentum is given by From the above discussion one can see that, the additional terms arising in Eq. (40) represent the energies of split levels which are a consequence of a perturbation collectively caused by quantum correlations, orbit-orbit interaction and the normal Zeeman effect. Lastly it is interesting to note that, if the external field B is strictly confined to only the region along the central axis of the helicoid (along x c -axis), then it is very likely to achieve a situation where no field penetrates the surface of the helicoid which may lead to a possible linkage of this geometric scheme to the Aharonov-Bohm effect. The energy equation Eq. (40) gives the dependence of energy spectrum on the vector potential instead of B and hence on magnetic flux through A = 2π , in accordance with Aharonov-Bohm effect.

Summary
In this work, the uncertainty of simultaneous measurement for various observables is quantified by computing their respective lower bounds using Robertson relation [2]. The classical dynamics of two constrained particle system under the influence of an external magnetic field is handled via Dirac formalism, the classical system is then quantized by Dirac quantization rule assuming the two particles constitute a maximally entangled bipartite system. We have shown how these uncertainties, in a constrained maximally entangled system on a helicoid, are effected under various geometric and non-geometric system parameters.
• Effect of σ : It was shown that the uncertainty between various observables of the two particles reduces as a function of σ . However, no reduction in uncertainty is possible for position and momentum observables of a single particle i.e., there is no intraparticle entanglement between these degrees of freedom. In view of a quantum to classical transition perspective, our results reveal that for Type-I equations the quantum regime persists for a system with a significantly large number of identical entangled particles constrained on helical geometries. Besides, the lower bounds of Type-I equations achieve σ -independence at the expense of a very large number of constrained particles i.e., in the limit σ → ∞. • Effect of B: It was shown that the external field B plays an important role in preserving quantum correlations and suitable strengths of A were revealed for tuning intraparticle as well as bipartite entanglement. Numerical bounds on the strength of B required for intraparticle and bipartite entanglement were given. We also show that for this system, it is possible to bring the uncertainty between momentum observables (of the different particles) to zero at A = 0, which corresponds to maximal bipartite entanglement. However, intraparticle entanglement can only be tuned in the presence of non-zero external fields.
• Effect of ρ: We showed that the product of dispersion (between position and momentum observables of the same particle) is susceptible to variations in the pitch for ρ ≤ 10, and achieves ρ-independence for large values of ρ. Moreover, It is observed that the bipartite entanglement exists for all ρ and is tunable by adjusting ρ. • Effect of radial position of particles: It has also been highlighted that the lower bounds are sensitive to changes in radial coordinates. The lower bounds are also affected by the pitch ρ, in such a way that increasing ρ reduces the lower bounds between position and momentum observables of the same particles. However, this reduction does not imply existence of intraparticle entanglement in view of the single particle case.
It has been observed that the bipartite entanglement is more robust to fluctuations in the fields as compared with the intraparticle entanglement. The energy equation is obtained and few additional energy terms arise in comparison to Eq. (42). These new energy terms arise due to coupling between the angular momentum of two particles, quantum correlations and the normal Zeeman effect which leads to the splitting of energy level into 11 levels. Finally, we remark that a suitable linkage might be possible between this geometric setup and the Aharonov-Bohm effect under the A dependence of energy, which might be an interesting avenue to explore further.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: The work is purely mathematical and all the techniques are standard. The simulation reported is also simple equation plotting which is straightforward. Thus there is no data required to be deposited.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 . SCOAP 3 supports the goals of the International Year of Basic Sciences for Sustainable Development.