Abstract

We consider a particle in a spatial symmetric/asymmetric potential driven by time periodic bichromatic AC fields of ratchet type. The associated time-dependent Schrödinger equation is conveniently tackled with the Floquet theory. We next proceed to investigate the ratchet effect induced by the driver, comparing the symmetric with the asymmetric cases. It turns out that the current in the asymmetric case is stronger than that of the symmetric one. Besides, we also investigate the case where the driver is a delta kicked acting on our spatial potential with more emphasis on its chaotic behaviour. Here we check that the current emerges as the phase space is mixed and that the system with asymmetric spatial potential becomes more chaotic than the symmetric one at low kicking strength.

1. Introduction

The ratchet effect is the possibility of generating directed current from a fluctuating environment in the absence of external unbalanced forces or any perturbation of any sort [1, 2].

This fluctuating environment is usually being induced by noisy dynamics governed by thermal fluctuations. This thermal fluctuation is usually brought about by nonlinear interaction between molecules. The Brownian ratchet has recently been discovered to be such a system that takes advantage of this noise as it converts random fluctuations into directed motion in the absence of external forces [1]. This concept was proposed to understand the functioning of molecular motors (e.g., kinesins) and has also found a wide range of applications in quantum system dynamics like nanoscale devices in which artificial ratchets may serve as new electron pumps, microscopic systems with lattices, experimental realization of quantum ratchets in semiconductor heterostructures, thermal rectifiers in molecular chains, biomotors and directed transport in molecular wires and proteins, molecular switch, particle selectors, transistors [3], laser matter interaction, and laser induced symmetric breaking [4]. This has also been largely applicable in studding the ratchet dynamics of cold atoms in optical lattices.

To observe the appearance of a DC current in a fluctuating environment of a dynamical system, the breaking of symmetries, that is, the spatial and temporal symmetries, is a necessary but not a sufficient condition. It has recently been shown that, in the area of atomic physics, directed transport has been achieved in driven ratchets with only the temporal symmetry being broken by dissipative processes [5]. Also in quantum transport a necessary condition for a nonzero transport is a mixed phase space with coexisting regular and chaotic regions [2]. This means that, for directed chaotic transport to occur, if phase space cannot be decomposed into invariant sets (chaotic sea and regular islands), there cannot be directed transport even with all symmetries broken.

Optical lattices are artificial crystals of light, consisting of hundreds and thousands of optical microtraps (potential wells) created by interfering optical laser beams. These optical lattices act as potential wells to trap cold quantum particles like Bosons and Fermions. Cold atoms in optical lattices provide an excellent experimental demonstration of the phenomenon of dynamical localization [5]. Dynamical localization is simply the quantum representation of classical chaos. Hamiltonian ratchet has recently been successfully realized with cold atoms in optical lattices with two harmonics driving and tunable weak dissipation [4]. Many motivated researchers in the domain of quantum Hamiltonian ratchets have mainly focused on the Hamiltonian kicked rotor model [4] which can be simulated easy. The quantum kicked rotor model is a model that best describes the phenomenon of dynamical localization, tunnelling of electrons, and quantum resonances [3].

2. The Model

Let us start this section with a cloud of atoms equally populating the wells of double-well periodic potential exposed to periodically driven two-color lasers. The cloud of atoms is said to be at low density [4]. Neglecting dissipation and nonlinear interactions among the particles, we can model the situation using the Schrödinger equation: where is rationalized Planck's constant ,   is the wave function, and is the model Hamiltonian given by where is the ratchet potential with being a phase factor. is the kicking strength, is a factor that breaks the symmetry of the system, and is the two-color laser of the form (see [4]), which is periodic that is with . Since the Hamiltonian is periodic in time, in other words, to solve (1), we employ the Floquet method. To start, we will revisit the Floquet formalism. According to the Floquet theorem [6], there exists a solution to (1) of the form where are the Floquet states that are the set of wave function solutions to the Schrödinger equation, are periodic Floquet modes, that is, , and are the quasienergy levels.

If the Floquet modes and quasienergies are known for a particular Hamiltonian , any initial wave function can be decomposed in the Floquet states and the solution is easily obtained for arbitrary : where the coefficients are determined by the initial wave function

Now substituting for into the Schrödinger equation, we have that and we can assign This gives us

This eigenvalue problem can be solved analytically or numerically.

The Floquet states are eigenstates of the Floquet operator where is the quasienergy. If and are solutions, so are and , where stands for frequency. are also called Dressed states.

Now let us take a look at the time evolution operator for Schrödinger equation (1) which satisfies Replacing the Floquet states by their expression in (5), we have but we know that , so we have which shows that the Floquet states are eigenstates of the one-period propagator [6]. Therefore we can then find Floquet modes and quasienergies by numerically integrating and diagonalizing it. This method is useful to find by calculating and diagonalizing .

The Floquet modes at arbitrary time can now be found by propagating to using the wave function propagator which for the Floquet modes yields Inverting the above equation, we have

By using the gauge transformation [4], where is the periodic vector potential. Substituting in (1), we have that Making use of the fact that , the above equation simplifies to Assigning which is our new Hamiltonian, we have that

Integrating the above equation, we have that which is the solution of the Schrödinger equation for the new Hamiltonian. Now assigning we call the propagator which acts on a wave function at an initial time to translate it to a final time as The propagator describes the time evolution of the wave function. forms a family of unitary operators, continues in , and satisfies the Chapman-Kolmogorov chain rule: If our Hamiltonian is time-independent, to get the propagator, we just necessarily need to integrate the equation which gives us with Now for the case where is periodic with a period of   like in our case, the form of the propagator is given above but satisfies

A section of the Floquet mode spectrum of asymmetric potential with values of the quasienergies on the top right-hand side of the graph (see Figure 1).

A section of the Floquet mode spectrum of asymmetric potential under periodic perturbation (see Figure 2).

3. Classical Ratchet

Solving the equation of motion for our Hamiltonian We have that Hamilton’s equation of motion is given by

This gives us

The current can be obtained by numerically solving nonlinear equation (33) above with the Runge Kutta Algorithm of Fourth Order.

For the purpose of simulation, we will restrict ourselves to the case where , , and and we will vary , and for our experiment. The graphs of the above results are displayed as follows.

In Figure 3 when comparing the symmetric (red line) and asymmetric (green line) potentials for very weak perturbation, we see that their ratchets are similar.

In Figure 4 we subjected the symmetric potential to both weak and strongly driving AC fields where we see that with strongly driven AC field the kinetic energy is greater than with weakly driven AC fields.

In Figure 5 we compare the ratchets of both symmetric and asymmetric double-well potentials in the presence of strongly driven AC fields. We realised that the ratchet effect in the presence of spatiotemporal asymmetry with (red) is greater than that with (green) which is greater than the symmetric double-well potential case (blue). Hence we can merely state that the ratchet effect in an asymmetric double-well potential is more than that in the symmetric case. We can affiliate this result to the concept of electron tunnelling which is not our main goal in this work.

In Figure 6, in the presence of a strong driving force with temporal symmetry conserved but spatial symmetry broken, the asymmetric potential with and yields similar ratchet effect but is greater than that with symmetry. Therefore we can conclude that the breaking of either the spatial or the temporal symmetry or both is a necessary condition for the appearance of a DC current.

4. Kicked Ratchet

The quantum Hamiltonian ratchet is based on the kicked rotor model. In this section, we are going to investigate a kicked rotor model for the Hamiltonian. By “kick” we mean a periodic impulse as forcing term [3]. This model is a prototype model for chaos and quantum chaos studies. It describes a particle that is constrained to move in a ring. Hence we are just going to replace our Cartesian coordinate and momentum with angular coordinate and momentum or . The particle is kicked periodically by a homogeneous field. The model is described by the Hamiltonian: where is the Dirac-delta function and is the kicking strength. Its dynamics is described by the standard map

The standard map is an area-preserving map of two generalised variables: position and momentum [7]. This map is periodic with a period of in both variables, so it can be thought of as being existing on a torus. The standard map describes a universal characteristic of area-preserving maps, namely, a phase space composed of regular islands embedded in a chaotic sea [4]. The standard map by varying a control parameter goes from exhibiting regular periodic behaviour to exhibiting unpredictable, chaotic behaviour. What is so unique about this map is that it is Hamiltonian or conservative; that is, the standard map is area-preserving. For strong enough kicks, , the system is chaotic and has a positive maxima Lyapunov exponent [7]. The Poincaré sections at different kicking strengths are displayed below. To look for the stable points, we have that This gives us This gives us the points ,  , and .

Figure 7 is the Poincaré sections for the periodically symmetric potential at a kicking strength of . This potential shows stable points.

In Figure 8, the asymmetric potential is at ; these stable points are replaced by a limiting circle. The trajectories are already showing chaos. This is an example of a pitchfork bifurcation.

In Figure 9, we see that the effect of increasing the kicking strength and the phase space for symmetric potential is now mixed showing the presence of regular islands and chaotic sea which is another necessary condition for nonzero transport. At this kicking strength, the stable point has now been replaced by limiting circles.

In Figure 10, at , the asymmetric potential shows a single stable island embedded in a chaotic sea.

In Figure 11, at , the symmetric potential shows partial full chaos.

In the asymmetric case, full chaos is achieved at (see Figure 12).

5. Conclusion

To conclude, we have studied the ratchet effect of cold atoms in periodically symmetric and asymmetric potentials driven by periodic bichromatic AC fields where it was found that the Floquet method has been preferred because of the periodicity of the time-dependent Schrödinger equation. In the classical limit, we show that the rectification effect of the ratchet current in asymmetric potential was more than that for a symmetric one. This result can be used to explain the concept of DC current enhancement due to incoherent electron tunnelling in an asymmetric potential, while, in the kicked ratchet, however, we found that the mix phase space of asymmetric potential becomes more chaotic at low kicking strength than of the symmetric case which equally explains why there is a great variation in their DC currents.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors would like to extend their sincere thanks to AIMS-NEI network and AIMS-Cameroon in particular for their wonderful financial and material support and also wish to knowledge God for His strength and inspirations.