1 Introduction
As a field of research, non equilibrium physics is simultaneously very new and very old. Non equilibrium physics is inherently interdisciplinary. It has grown into a major enterprise. Non equilibrium systems that form patterns or exhibit complex behaviors are especially sensitive to small perturbations; thus their detailed behavior may be extremely challenging to predict or control, and as such is therefore a very interesting field of research. Pattern formation is ubiquitous in nature (Cross & Hohenberg Reference Cross and Hohenberg1993); it occurs for instance in systems undergoing convection, in crystals developing in supercooled or supersaturated environments and even in growing biological systems. In gas discharge systems, patterns often can be observed in the glow emitted by the plasma, for example, in direct current discharges with one electrode made of a semiconductor material (Astrov et al. Reference Astrov, Müller, Ammelt and Purwins1998; Strümpel, Astrov & Purwins Reference Strümpel, Astrov and Purwins2002) and in dielectric barrier discharges (Ammelt, Schweng & Purwins Reference Ammelt, Schweng and Purwins1993; Breazeal, Flynn & Gwinn Reference Breazeal, Flynn and Gwinn1995). Lane formation (Sarma, Baruah & Ganesh Reference Sarma, Baruah and Ganesh2020) is another representative example of non-equilibrium phase transition. When two kinds of particles are driven in opposite directions, the system exhibits a self-organization from a uniformly mixed state into strongly ordered anisotropic patterns. This driven segregation phenomenon, first found in computer simulations (Dzubiella, Hoffmann & Löwen Reference Dzubiella, Hoffmann and Löwen2002; Sarma et al. Reference Sarma, Baruah and Ganesh2020) has been observed in laboratory experiments, such as mixtures of oppositely charged colloids (Leunissen et al. Reference Leunissen, Christova, Hynninen, Royall, Campbell, Imhof, Dijkstra, Roij and van Blaaderen2005a) or dusty plasmas in the presence of an external electric field (Sütterlin et al. Reference Sütterlin, Wysocki, Ivlev, Räth, Thomas, Rubin-Zuzic, Goedheer, Fortov, Lipaev and Molotkov2009). For a fundamental understanding of such non-equilibrium phenomena requires the study of non-equilibrium phase transition phenomena from ‘first principle’ approach, such as Molecular Dynamics (MD) simulation, which integrates 
$N$ particle system without any ansatz or assumptions (Sarma et al. Reference Sarma, Baruah and Ganesh2020).
Phase transitions are common in equilibrium systems, and the related analogies have represented an important contribution to the understanding of non-equilibrium processes. While equilibrium bulk phase transitions are now well understood both by computer simulations (Allen et al. Reference Allen, Allen, Tildesley, Allen and Tildesley1989) and by statistical theories (Hansen & McDonald Reference Hansen and McDonald1986; Löwen Reference Löwen1994), non-equilibrium situations may induce a much richer scenario of phase transformations. Until now, very extensive theoretical and simulation studies have been performed on such transitions in driven diffusive systems (Katz, Lebowitz & Spohn Reference Katz, Lebowitz and Spohn1983; Marro & Dickman Reference Marro and Dickman1999), particularly, models of identical particles have been studied, which couple with a different sign to an external uniform field (Aertsens & Naudts Reference Aertsens and Naudts1991), and the particles then form stripes perpendicular to the field direction. Later, this transition has been put forward as a concept of panic theory applied to pedestrian zones (Helbing, Farkas & Vicsek Reference Helbing, Farkas and Vicsek2000). Similar situations with two groups of particles moving in the other direction have been extensively studied in the context of lane and pattern formations not only of pedestrians (Ikeda & Kim Reference Ikeda and Kim2017; Feliciani, Murakami & Nishinari Reference Feliciani, Murakami and Nishinari2018), but also of various physical particles, such as charged colloids (Vissers et al. Reference Vissers, Wysocki, Rex, Löwen, Royall, Imhof and van Blaaderen2011b; Vissers, van Blaaderen & Imhof Reference Vissers, van Blaaderen and Imhof2011a; Tarama, Egelhaaf & Löwen Reference Tarama, Egelhaaf and Löwen2019), granular systems (Aranson & Tsimring Reference Aranson and Tsimring2006), migrating macroions (Netz Reference Netz2003), microswimmers (Kogler & Klapp Reference Kogler and Klapp2015) and plasmas (Sarma et al. Reference Sarma, Baruah and Ganesh2020). However, the mechanism of phase transition is much less clear for non-equilibrium conditions. Additionally, the phase stability of pair plasmas in the strongly coupled regime has attracted special attention. In realization of the fact, Swati Baruah et al. (Baruah, Ganesh & Avinash Reference Baruah, Ganesh and Avinash2015) carried out investigations for addressing phase transitional study in strongly coupled two-dimensional (2-D) pair-ion plasma (PIP) system with a soft core. Furthermore, our recent work on a self-organization effect (Sarma et al. Reference Sarma, Baruah and Ganesh2020) where the ‘plasma’ is structured as a set of lanes that alternatively contain positive and negative ions provides an opportunity for the in-depth fundamental understanding of non-equilibrium lane formation phenomena in the naturally occurring PIP systems.
PIP is a novel state of multi-ion plasma. This kind of plasma is made up of opposite charges and equal masses. The PIPs represent an interesting state of matter with unique thermodynamic properties drastically different from ordinary electron–ion plasmas (Iwamoto Reference Iwamoto1993). The PIP plays a significant role in plasma physics, for example, numerous astrophysical environments such as the pulsar magnetosphere, active galactic nuclei, neutron stars, etc. (Piran Reference Piran2005), where intense energies create electrons and positrons through pair production and annihilation. In the last few years, a lot of research work has been performed to explore the linear as well as nonlinear properties of PIPs so as to understand such types of astrophysical plasmas (Sahu et al. Reference Sahu, Pal, Poria and Roychoudhury2015). In semiconductors, non-thermal pair plasmas are ubiquitous in the form of electron and ion holes (Piran Reference Piran2005). The PIP is expected to be used for the synthesis of dimers directly from carbon allotropes, as well as for nanotechnology (Oohara & Hatakeyama Reference Oohara and Hatakeyama2003a). The physics of the PIP became more interesting when its application extended to the terrestrial laboratory. Generally, attention has been focused largely on the relativistic and weakly relativistic regime (Iwamoto Reference Iwamoto1993; Zank & Greaves Reference Zank and Greaves1995; Asseo Reference Asseo2003) of PIP, since such plasmas are thought to be generated naturally under certain astrophysical conditions (Asseo Reference Asseo2003). In many of the cases above, magnetic fields – static as well as dynamic fields – are known to play an important role (Pilch et al. Reference Pilch, Piel, Trottenberg and Koepke2007; Pilch, Reichstein & Piel Reference Pilch, Reichstein and Piel2008; Baruah & Das Reference Baruah and Das2010; Ott & Bonitz Reference Ott and Bonitz2011). In PIPs, pair annihilation can take place (Zank & Greaves Reference Zank and Greaves1995). This process is analogous to recombination in electron–ion plasmas. Another interesting fact is that in such plasmas the nonlinear phenomena such as solitary wave structures in stationary plasmas and in the beam plasma system, cross-field transport, and the phenomena of holes and clumps are different from electron–ion plasmas (Berman, Tetreault & Dupree Reference Berman, Tetreault and Dupree1985). Many of the phenomena found in conventional electron–ion plasmas exist in modified form, while others, notably the whistler wave, the lower hybrid wave and Faraday rotation, are absent (Zank & Greaves Reference Zank and Greaves1995).
A number of experimental approaches have been proposed for studying PIP in the laboratory, and a few of them are now being actively pursued. Oohara & Hatakeyama (Reference Oohara and Hatakeyama2003b, Reference Oohara and Hatakeyama2007) reported that fullerene PIP 
$C^{\pm }_{60}$ can be produced experimentally. Further, electron-free, positive and negative-ion plasmas have been achieved in the afterglow of pulsed-power 
$Cl_{2}$ discharges (Kanakasabapathy & Overzet Reference Kanakasabapathy and Overzet1998), and some control on the amount of background electron population has also been achieved. Remaining very low density background electrons tend to play the role of ‘shielding’. Recently, in another experimental work, PIP with positive and negative Hydrogen ions 
$(H^{+}+H^{-})$ has been successfully created (Oohara & Hatakeyama Reference Oohara and Hatakeyama2007). In this work, it is assumed that pair-ions interact with a potential which is of Yukawa form, where the shielding of the PIP particles is believed to come from the ‘background’ electrons, as observed in the 
$C_{60}^{\pm }$ experiment (Oohara, Date & Hatakeyama Reference Oohara, Date and Hatakeyama2005).
It is well proven that an external magnetic field influences the dynamics of the charged particles as well as their collective behaviour in plasma systems. Over the last decade, a series of research focusing on the effect of the externally applied magnetic field over the particle dynamics has been reported. For example, dusty plasma offers an excellent medium for studying various phenomena at single particle and fluid levels with remarkable temporal and spatial resolutions in the presence of an external magnetic field. An experimental study of fine particles in magnetized plasma was done by Sato et al. (Reference Sato, Uchida, Kaneko, Shimizu and Iizuka2001). They found that fine particles rotate in the azimuthal direction on the horizontal plane in presence of a magnetic field applied in the vertical direction. In their experiment, the rotation is induced under the condition that the particles are strongly coupled with very small interparticle distance. Konopka et al. (Reference Konopka, Samsonov, Ivlev, Goree, Steinberg and Morfill2000) studied the effect of a vertical magnetic field on the particle cloud in the sheath of a radio frequency He plasma at various discharge conditions. They observed that the cloud rotates under the influence of the magnetic field. Uchida, Konopka & Morfill (Reference Uchida, Konopka and Morfill2004) on the other hand, studied the wave dispersion relation in a 2-D strongly coupled plasma crystal by using both theoretical analysis and MD simulation, taking into account a constant magnetic field parallel to the crystal normal. Very recently, a dusty plasma device capable of operating with magnetized electrons, ions and charged dust particles was designed at Auburn University, USA, under the direction of Thomas et al. (Reference Thomas, Lynch, Konopka, Menati, Williams, Merlino and Rosenberg2019), to study of highly magnetized plasmas and dusty plasmas, where, experiments can be performed at magnetic field up to 4T. Their study shows the formation of new types of self- and imposed-ordered structures that form in both plasma and among the microparticles.
The PIP in a magnetic field is unique in that both species are equally and strongly magnetized (Zank & Greaves Reference Zank and Greaves1995). Additionally, the behaviour of such a kind of plasma system in the presence of the external magnetic field in 2-D is largely unexplored. Therefore, here, we focus on a 2-D PIP system and explore the lane formation dynamics, using computer simulations, under the influence of an external magnetic field. Here, the entire work is performed for 
$\kappa \sim 0$ (
$\kappa =10^{-4}$), where, 
$\kappa ={a}/{\lambda _{D}}$ is the screening parameter resulting from the shielding dynamics of ‘background’ charged particles (Sarma et al. Reference Sarma, Baruah and Ganesh2020), 
$\lambda _{D}$ is the Debye length and ‘
$a$’ is the Wigner–Seitz radius; for the rest of the presentation this is the value of 
$\kappa$ used. In this study, the particle motion is governed by Langevin dynamics (LD). In our previous work (Sarma et al. Reference Sarma, Baruah and Ganesh2020), lane formation dynamics of driven 2-D PIPs were investigated in the underdamped case, where we studied the behaviour of the system in the presence of both constant and time-varying external electric field taking into consideration the finite screening case (
$\kappa \neq 0$). There, comparative study performed between the results in the overdamped and underdamped limit brought new and interesting information as far as the self-organization effect is concerned. Here, we focus on some aspects of the lane formation dynamics of a PIP system in a magnetic field in the presence of both constant and oscillatory electric fields in zero 
$\kappa$ limit. First we perform our study by applying the external electric field (both constant and oscillatory) in the 
$Y$-direction, considering the external magnetic field direction is out of the simulation chamber, and obtain the critical value of the strength of the external electric field above which lanes are formed distinctly in the presence of the external magnetic field. The resulting critical electric field strength is then compared with the one in the absence of the external magnetic field and it is observed that the presence of the external magnetic field retards the phase transition process and leads to a disintegration of lane formation. Through this present study, we expect that our results may provide an opportunity for an in-depth fundamental understanding of non-equilibrium lane formation phenomenon in the naturally occurring PIP systems as well as for their relevance to technological applications. Additionally, due to the appearance of electric field (
$E \times B$) drift in our system, the self-organized lanes start to move along a direction perpendicular to both the directions of the applied electric and magnetic fields. It is also observed that the self-organization is correlated to a strong decrease of the transverse diffusion coefficient in the absence of the external magnetic field. However, due to the appearance of the electric field drift, the value of the long-time self-diffusion coefficient increases along with the electric field strength after getting into the self-organized phase. On the other hand, from a Coulomb coupling parameter varying order parameter plot we observed the existence of a lane formation-disintegration phase of the PIP system both in the presence and absence of the external magnetic field, taking into consideration of the constant electric field case, which reveals that, in our system, a range of critical coupling parameters exists in which only the system shows a self-organized phase.
Using Langevin Dynamics (LD) simulation technique, we further demonstrate the existence of a non-stationary lane-like self-organization behaviour, with oscillation between self-organized and random states for the small frequency values, when submitted to a time-varying electric field and,as a result of the guiding centre drift, lanes are found to move along a direction perpendicular to both the electric and magnetic field directions. It is also seen that, for each electric field there exists a critical frequency above which the system cannot reach the self-organized state and the presence of the external magnetic field significantly affects the critical frequency value. However, in second part of the paper we discuss the results obtained from the application of the oblique external electric field (both constant and oscillatory field) to the PIP system in the presence of the external magnetic field. Our motivation for studying a PIP system comes from the fact that the magnetic field may affect the lane formation dynamics of these particle species differently. On the other hand, phase transitions in conventional and model complex fluids are of significant scientific and technological interest. The phase transition mechanism in a driven system reveals a wealth of novel instabilities and pattern formations induced by non-equilibrium conditions (Sarma et al. Reference Sarma, Baruah and Ganesh2020). However, a complete understanding of the key parameters which determine the existence and order of non-equilibrium phases still remains to be explored fully in a PIP system and is an active area of research. Additionally, the behaviour of such a kind of plasma system in the presence of the external magnetic field in two dimensions is largely unexplored. Therefore, for this study we have chosen a pair-ion system to explore the lane formation dynamics, using computer simulations, under the influence of an external magnetic field.
The rest of the paper is organized as follows: In § 2, we define the model used. The simulation technique, along with the lane formation order parameter and long-time self-diffusion coefficient, are discussed in § 3. In § 4, the results are presented and discussed in detail. Finally we finish with conclusions in § 5.
2 The model
We simulate the Langevin molecular dynamics for a 2-D PIP system of 
$N$ charged particles in a regular square of length 
$X=Y=L$ keeping the total density 
$\rho$ and the temperature constant. The periodic boundary condition is used in both the 
$X$ and 
$Y$-directions. Here, half of the PIP particles carry positive charge and the other half carry a negative charge of the same magnitude 
$|Q|$ and of equal masses subjected to a perpendicular external magnetic field 
${\boldsymbol {B}}=B\hat {e}_{z}$ (
$\hat {e}_{z}$ is the unit vector along 
$Z$-direction). The number density is 
$\rho = N/L^{2}$, where 
$L$ is the dimension of the system.
The particle trajectories are governed by the LD
\begin{equation} m_{i}\frac{\textrm{d}^{2}\boldsymbol{r}_{i}}{\textrm{d}t^{2}}={-}\gamma \frac{\textrm{d}\boldsymbol{r}_{i}}{\textrm{d}t}-\sum_{i < j}\boldsymbol{\nabla} U(r_{ij})+\boldsymbol{F}_{\textrm{ext}}+\boldsymbol{F}^{(R)}_{i},\end{equation}
where 
$i=1,\ldots ,N$, 
$m_{i}$, 
$r_{i}$ are, respectively, the number of PIP particles, the mass and position of the 
$i\textrm {th}$ particle, 
$-\sum _{i < j}\boldsymbol {\nabla } U(r_{ij})$ is the force comprised of interparticle interactions. As is well known, in Langevin model, each PIP particle experiences a Brownian motion from the background interactions which results in dissipation of PIP particles via 
$-\gamma ({\textrm {d}\boldsymbol {r}_{i}}/{\textrm {d}t})$, here 
$\gamma$ is the friction constant. The force, comprised of the external force fields, is
\begin{equation} \boldsymbol{F}_{\textrm{ext}}=\boldsymbol{F}_{E}+\boldsymbol{F}_{B}, \end{equation}
were 
$\boldsymbol {F}_{E}$ is the space-independent external electric force acting on the particles pointing in the 
$Y$-direction, and is modelled as
\begin{equation} \boldsymbol{F}_{E}(t)=\hat{y}E_{0}\cos(\omega t)\end{equation}
where 
$E_{0}$ is the strength or amplitude of the external electric field, 
$\omega$ is the frequency of the external electric field (with 
$\omega \equiv 0$ leading to the constant-field case) and 
$\hat {y}$ is the unit vector along the 
$Y$-direction. Temporally oscillatory electric fields are often found both in experiments (either controlled form or in noise form) and in other naturally occurring PIPs, such as in active galactic nuclei (Begelman, Blandford & Rees Reference Begelman, Blandford and Rees1984), in stars (Michel Reference Michel1991), in pulsar magnetospheres (Paul N. Arendt & Eilek Reference Paul N. Arendt and Eilek2002), in gamma ray fire balls (Piran Reference Piran1999), in solar flares (Tandberg-Hanssen & Emslie Reference Tandberg-Hanssen and Emslie1988), etc. As is well known, any bounded function 
$f(t)$ may be represented as aFourier sum of frequencies 
$\omega$ and their appropriate weights. Thus, performing numerical studies with a cosine or sine function of 
$\omega$ (see (2.3)) should be able to give generic insight for either a controlled or noise form of external electric field. Additionally, as the occurrence of lane formation appears to be very general, this non-equilibrium transition could be experimentally verified in quite different systems ( 2001) as well as observable in real systems, e.g. driven by an oscillatory electric field. Hence, in our case, for a fundamental understanding of the lane formation phenomenon occurs in realistic system (here, in a naturally occurring PIP system), the electric field is modelled as mentioned in (2.3). The force due to the external magnetic field applied along the 
$Z$-direction (i.e. out of the simulation box) is
\begin{equation} \boldsymbol{F}_{B}=\frac{Q_{A,B}B}{m_{i}c}\dot{r}_{i}\times\hat{e}_{z}, \end{equation}
where 
$c$ is the speed of light. The force 
$\boldsymbol {F}^{(R)}_{i}$ describes random kicks of the background acting on the 
$i$th PIP particle by the ‘bath’. These kicks are mimicked by Gaussian random numbers with zero mean, 
$\overline {\boldsymbol {F}^{(R)}_{i}}=0$, and variance
\begin{equation} \overline{(\boldsymbol{F}^{(R)}_{i})_{\alpha}(t)(\boldsymbol{F}^{(R)}_{j})_{\eta}(t^{\prime})}=2k_{B}T\gamma\delta_{\alpha \eta}\delta_{ij}\delta(t-t^{\prime}). \end{equation}
The subscripts 
$\alpha$ and 
$\eta$ stand for two Cartesian components, 
$i$ and 
$j$ stands for 
$i\textrm {th}$ and 
$j\textrm {th}$ particles, 
$t$ and 
$t^{\prime }$ for two distinct time, where the first and second 
$\delta$ are Kronecker and the third 
$\delta$ stand for Dirac delta function and 
$k_{B}T$ is the thermal energy of the bath.
The interaction potential reads as
\begin{equation} U(r_{ij})=\left(\frac{Q_{A}Q_{B}}{4{\rm \pi}\epsilon_{0}r_{ij}}\right) \exp\left({-\kappa\left(\frac{r_{ij}}{a}-1\right)}\right),\end{equation}
where, 
$r_{ij}$ is the distance between 
$i{\textrm {th}}$ and 
$j{\textrm {th}}$ pair, 
$Q_{A, B}$ is the charge of particle 
$A$ 
$(B)$, such that 
$Q_{A}Q_{B}=-1$ and 
$Q_{A}Q_{A} = Q_{B}Q_{B} = +1$. This is a valid model for charge-stabilized suspensions in three (Löwen & Kramposthuber Reference Löwen and Kramposthuber1993) and two (Lowen Reference Lowen1992) spatial dimensions. It has also been shown experimentally that the potential as mentioned above ((2.6), i.e. Yukawa potential) is a suitable effective interaction for this kind of system (Leunissen et al. Reference Leunissen, Christova, Hynninen, Royall, Campbell, Imhof, Dijkstra, van Roij and van Blaaderen2005b; Royall et al. Reference Royall, Leunissen, Hynninen, Dijkstra and van Blaaderen2006). Here, the length, time, energy, density and electric field are normalized, respectively, by 
$a$, 
$\omega _{p}^{-1}$, 
$k_BT$, 
$a^{2}$ and 
${Q}/{4{\rm \pi} \epsilon _{0}a^{2}}$. In these units, the strength of the interaction potential is defined as 
$\varGamma$, where 
$\varGamma ={Q_{A,B}^{2}}/{4{\rm \pi} \epsilon _{0}ak_{B}T}$; here, 
$k_{B}T$ denotes the thermal energy. In this study, the strength of the magnetic field 
$B$ is given by 
$\beta ={\omega _{c}}/{\omega _{p}}\propto B$, i.e. the ratio of the cyclotron frequency 
$\omega _{c}={Q_{A,B}B}/{m_{i}c}$ (here, 
$c$ represents speed of light) and the plasma frequency 
$\omega _{p}=[{2Q_{A,B}^{2}}/{4{\rm \pi} \epsilon _{0}a^{3}m_{i}}]^{1/2}$. Here ‘
$a$’ is the Wigner–Seitz radius, 
$\epsilon _{0}$ is the permittivity of free space and 
$Q_{A,B}$ is the charge of particle 
$A(B)$. Here, simulations are performed with 
$\beta$ values varying from 0.001 to 1.0. In real units, 
$\beta =1.0$ is equivalent to magnetic field strength 
$B=4.6388$ T for the PIP system as created in the 
$C^{\pm }_{60}$ experiment (Oohara et al. Reference Oohara, Date and Hatakeyama2005). As mentioned in the introductory section above, already experimental works (Thomas et al. Reference Thomas, Lynch, Konopka, Menati, Williams, Merlino and Rosenberg2019) are going on where highly magnetic fields with strength 
$4$ to 5 T are used; therefore, the range of magnetic field strength that is used in this study can be considered as reasonable enough to describe this problem.
For this study, the stochastic Langevin normalized equations for the particle trajectories 
$\boldsymbol {r}_{i}(t)$ 
$(i = 1,\ldots , N)$ read as
\begin{equation} \frac{\textrm{d}^{2}\boldsymbol{r}_{i}}{\textrm{d}t^{2}}={-}\gamma\frac{\textrm{d}\boldsymbol{r}_{i}}{\textrm{d}t} +\sum_{i < j}\frac{\varGamma}{r_{ij}}\left(\kappa+\frac{1}{r_{ij}}\right) \exp[-\kappa(r_{ij}-1)]+\boldsymbol{F}_{\textrm{ext}}(t)+\boldsymbol{F}^{(R)}_{i}(t). \end{equation} In this work, we study two cases in the presence of the external magnetic field: in first case, the external constant or oscillatory electric force 
$\boldsymbol {F}_{E}$ is pointing in the 
$Y$-direction as discussed in (2.3) above; in second case, the external electric field 
$E_{0}\cos (\omega t)$ intersects the 
$X$–
$Y$ plane at an oblique angle 
$\theta$. The force due to the oblique electric field reads as
\begin{equation} \boldsymbol{F}_{E}=\boldsymbol{F}_{\textrm{EOblique}}=F_{\textrm{Ex}}\hat{e}_{x}+F_{\textrm{Ey}}\hat{e}_{y}. \end{equation}The corresponding components of the oblique electric field are modelled as
\begin{equation} \left.\begin{gathered} F_{\textrm{Ex}}=E_0\cos(\omega t)\cos(\theta),\\ F_{\textrm{Ey}}=E_0\cos(\omega t)\sin(\theta). \end{gathered}\right\} \end{equation}3 Simulation procedure
Extensive LD simulations using the OpenMP parallel program are carried out to analyse the effect of the external magnetic field (
$\beta \neq 0$) on the non-equilibrium phase transition which is associated with lane formation parallel to the external field of the 2-D PIP. The simulations are performed using 
$N=500$ particles in a square simulation cell of side 
$L$. In the process of deciding the appropriate number of particles 
$N$, we have performed a series of runs from 
$N=250$ to 
$2000$ (
$N=250$, 500, 600, 700, 1000 and 2000) at fixed density 
$\rho$ and found that for particles 
$N \geqslant 500$, the finite size effects tend to become negligible. Hence 
$N=500$ has been chosen for this study (Sarma et al. Reference Sarma, Baruah and Ganesh2020). After equilibration with the desired temperature 
$\varGamma ^{-1}$, using the Langevin equation, the force due to the external electric field 
$\boldsymbol {F_{E}}$ is applied to the particles. Here, the external magnetic field is applied from the beginning of the simulation until the entire study has been performed. The lane or structure order parameter (defined later) is evaluated for various values of the applied electric field strength 
$E_0$ in the presence of the external magnetic field. As shown earlier (Sarma et al. Reference Sarma, Baruah and Ganesh2020), here also it is found that beyond a critical value of 
$E_0$, the order parameter tends to unity. It is also observed that the critical electric field strength value obtained from our simulations is independent of the system size, which suggests that 
$N=500$ is enough to describe this problem and to minimize the finite size effect.
To begin with, each particle is assigned an initial velocity which is random in direction, such that the average kinetic energy corresponds to the chosen temperature. The Langevin equations of motion including the external electric and magnetic fields are solved using the velocity Verlet method with a sufficiently small time step 
$\Delta t$. The simulation is done in the following steps.
For both the constant and oscillatory field cases including the oblique force field cases, with and without the external magnetic field (for a typical run, say, at 
$\varGamma =2.5$, 
$\rho =0.5$) typically 
$8\times 10^{5} \Delta t$ steps are simulated which corresponds to a simulation time of 
$t\omega _{\textrm {pd}}=2400$ with 
$\Delta t=3\times 10^{-3}$. After equilibration with the thermal energy using the Langevin equation, the external force field is applied to the particles at a time 
$t=2\times 10^{5}\Delta t$ steps. The data collection is obtained for the last 
$2\times 10^{5}\Delta t$ steps after the system has reached a far-from-equilibrium quasi-steady-state.
3.1 Order parameter
To characterize the degree of lane formation, a convenient order parameter is proposed for evaluating the number of particles of the same species along the driven direction (the 
$Y$-direction here). An order parameter in effect counts the percentage of particles in a lane-like configuration. It helps to distinguish two different phases (or orders) (Sarma et al. Reference Sarma, Baruah and Ganesh2020). For example, in one phase, the value of the order parameter may be zero and in another phase it may attain a non-zero value. Following Ikeda & Kim (Reference Ikeda and Kim2017) and Sarma et al. (Reference Sarma, Baruah and Ganesh2020), we define the following as the order parameter of lane formation. The simulation system is divided into 
$n_{\textrm {div}}$ numbers of identical rectangular strips along the driven direction (i.e. along 
$\hat {y}$). Thus, the width of each strip along the 
$x$-direction is given by 
$l_{\textrm {div}} = L/n_{\textrm {div}}$. The order parameter for the lane formation, say, in the 
$k\textrm {th}$ strip (
$\phi _{k}$) is determined by the difference in numbers between the positively and negatively driven particles, which is normalized by the total particle number in the 
$k$-th strip as
\begin{equation} \phi_{k}=\frac{n^{+}_{k}-n^{-}_{k}}{n^{+}_{k}+n^{-}_{k}}, \end{equation}
where 
$n^{\pm }_{k}$ represents the number of ‘
$\pm$’ species in the 
$k\textrm {th}$ rectangular strip. Thus, the range of 
$\phi _{k}$ is from 
$-1$ to 1, and it becomes positive if the positive particles become dominant in the strip, and negative if the negative particles become dominant. The order parameter 
$\phi$ for the total system is then defined as
\begin{equation} \phi=\left\langle\frac{1}{n_{\textrm{div}}}\sum^{n_{\textrm{div}}}_{k=1}\mid\phi_{k}\mid\right\rangle, \end{equation}
where the bracket denotes a time average. The order parameter 
$\phi$ begins at around 
$0$ at the initial state since oppositely charged particles are found inside the system, and finally reaches 
$1$ when the system exhibits the lane formation in the steady state. The former corresponds to the randomly mixed state, whereas the latter corresponds to the lane formation structure by the particles of the same species along the driven direction. We note that our definition of 
$\phi$ is essentially the same as those introduced previously (Dzubiella et al. Reference Dzubiella, Hoffmann and Löwen2002; Ikeda, Wada & Hayakawa Reference Ikeda, Wada and Hayakawa2012; Sarma et al. Reference Sarma, Baruah and Ganesh2020).
3.2 Long-time self-diffusion coefficient (
$D_{\textrm {xL}}$)
The long-time self-diffusion coefficient 
$D_{\textrm {xL}}$ perpendicular to the external field direction is defined as (Sarma et al. Reference Sarma, Baruah and Ganesh2020)
\begin{equation} D_{\textrm{xL}} = \lim_{t \to \infty} \frac{1}{2t} \left\langle\sum_{i=1}^{N}[(r_{i}(t)-r_{i}(0))\cdot \hat x]^2\right\rangle, \end{equation}
where 
$r_{i}(t)$ is the trajectory of the 
$i\textrm {th}$ particle at an instant 
$t$. Here, 
$\langle \ldots \rangle$ is the equilibrium average for a system unperturbed by the external force.
The calculation of the long-time self-diffusion coefficient 
$D_{\textrm {xL}}$ perpendicular to the applied field direction will help us to distinguish between the mixed orderless state and lane state (Sarma et al. Reference Sarma, Baruah and Ganesh2020). The details of the working of the long-time self-diffusion coefficient 
$D_{\textrm {xL}}$ in this context is discussed in the § 4.
4 Results and discussions
In this paper, we perform our study to see the effect of an external magnetic field on lane formation dynamics in PIP using both constant and oscillatory external electric fields using the Langevin equation of motion (2.1). This work will primarily focus on finding out the parameter regime of existence of lane structures in a PIP system in the zero screening limit in the presence of the external magnetic field and details of it is discussed in § 4.1. From the comparative study between the results obtained from the presence and absence of the external magnetic field it is observed that a higher value of the electric field strength (constant field case) will be needed to enter into the lane formation regime in the presence of the external magnetic field, however, comparatively at a lower value of the electric field frequency (known as the critical frequency for oscillatory field case), the system exhibits a transition back to the disordered state in the presence of the external magnetic field. The paper will then discuss how this lane formation mechanism will be significant in the presence of the oblique external electric field both in the presence and absence of the external magnetic field in § 4.2.
4.1 Lane formation in the presence of the external magnetic field ($\beta \neq 0$)
4.1.1 Constant field
In this subsection, we report on comprehensive simulation studies of lane formation in a PIP system, that are carried out in the presence of a constant external electric field. Here, we focus the influence of the external magnetic field on the lane dynamics. The magnetic field is applied from the beginning of the simulation in every case being studied.
Three simulation snapshots in the steady state for increasing external electric forces and with 
$\beta =1.0$ are shown in figure 1. Figure 1 shows the simulation snapshots associated with a situation with an external electric field below the critical field point (
$E_{c}$), at the critical field point and above the critical field strength, in the presence of the external magnetic field (
$\beta =1.0$). As shown earlier (Sarma et al. Reference Sarma, Baruah and Ganesh2020) for the finite 
$\kappa$ and 
$\beta =0$ case, in this study one also clearly sees lane formation parallel to the external electric field in the presence of the external magnetic field. These figures indicate that large 
$E_{0}$ leads to lane formation. With increasing 
$E_{0}$, the collision between oppositely driven particles was enhanced. This effect resulted in a bundle of particles in the same direction and generated the interface between oppositely driven particles; this is the main mechanism for the lane formation. From figure 1, it is seen that the thickness of the lanes perpendicular to the driven direction involve several particle layers for the external field strength 
$E_{0}>E_{c}$ and is comparable to the size of the simulation box. This implies that lane formation is a macrophase separation. However, below the critical strength needed to generate lane formation, some very weak patterns are generated with a width of one particle. Moreover, it is also observed that due to the presence of both uniform electric and magnetic fields the PIP particles experience not only the usual circular Larmor gyration, but also a drift of the guiding centre (
$E\times B$ or electric field drift), and as a result of which the lanes seem to move along the 
$X$-direction of the simulation chamber. As the magnitude of this electric field drift of the guiding centre is generally defined as 
$v_{E}=E_{0}/B$, therefore, from our simulation results it is observed that with increasing electric field strength 
$E_{0}$ the movement of the lanes towards the 
$X$-direction increases. In order to characterize the degree of lane formation in the presence of an external magnetic field (
$\beta =1.0$), the order parameter (
$\phi$) is calculated for various values of the electric field strength (
$E_{0}$) ranging from 1.0 to 200, taking density value (
$\rho =1.0$) as shown in figure 2, both in the presence (
$\beta =1.0$) and absence (
$\beta =0.0$) of the external magnetic field. In real units, the above-mentioned range of 
$E_{0}$ examined here are equivalent to 
$E_{0} \approx$ 14.405 to 2881.133 volts per metre, for the PIP system as created in the 
$C^{\pm }_{60}$ experiment (Oohara et al. Reference Oohara, Date and Hatakeyama2005). It is observed that, for both the cases, 
$\phi$ is small for a smaller value of field strength, and it grows when a critical electric field strength (
$E_{c}$) is approached. As observed, with turning on the external field strength 
$E_{0}$, PIP particles only follow the external field by eluding each other, so that the diffusion perpendicular to the external field has to increase with 
$E_{0}$. This effect grows until the critical field strength is reached and the system begins to form lanes. The critical value of the electric field strength 
$E_{c}$ is determined at the electric field strength value where the generation of lane formation has just started, and the corresponding value of the order parameter is recorded as 
$\phi =0.5$ for the simulation runs with increasing field strength. It is noted that this particular value is also considered as the threshold value to characterize the critical lane value by other authors (Dzubiella et al. Reference Dzubiella, Hoffmann and Löwen2002; Rex & Löwen Reference Rex and Löwen2008; Klymko, Geissler & Whitelam Reference Klymko, Geissler and Whitelam2016; Ikeda & Kim Reference Ikeda and Kim2017; Sarma et al. Reference Sarma, Baruah and Ganesh2020). Thus, the value 
$\phi =0.5$ plays a significant role in identifying the phase transition point. The critical field strength 
$E_{c}$ thus obtained in the presence of the external magnetic field (
$\beta =1.0$) is 
$E_{c}= 70$. However, the location of the phase transition in the absence of the external magnetic field (
$\beta =0.0$) is found at electric field strength value 
$E_{c}=60$, which indicates that the presence of the external magnetic field retards the phase transition process, as can be verified from figure 2.
Typical simulation snapshots of the 2-D PIP system with constant external electric field (
$\omega =0$), taking 
$\beta =1.0$, 
$\varGamma =2.5$, 
$\rho =1.0$. Panel (a) shows disordered state with electric field 
$E_{0} < E_{c}$ (
$\phi =0.12 < 0.5$), panel (b) shows a less disordered state with electric field 
$E_{0}=E_{c}$ (
$\phi =0.5$) and panel (c), shows ordered lane formation with 
$E_{0} > E_{c}$ (
$\phi =0.85 > 0.5$).
Order parameter 
$\phi$ versus electric field strength 
$E_{0}$ plot for constant electric field case (
$\omega =0$) both in the presence of the external magnetic field (
$\beta =1.0$) and in the absence of the magnetic field (
$\beta =0.0$), taking other parameters as 
$\varGamma =2.5$, 
$\rho =1.0$.
Here, the long-time self-diffusion coefficient 
$D_{\textrm {xL}}$ is studied along 
$\hat {x}$, i.e. in the direction perpendicular to both the external electric field (along 
$\hat {y}$) and magnetic field (along 
$\hat {z}$). It measures the mean square particle displacement along 
$\hat {x}$ during a characteristic time interval. Figure 3 shows a comparative plot for the perpendicular long-time self-diffusion coefficient (
$D_{\textrm {xL}}$), both in the presence and absence of the magnetic field, and the other parameters used are 
$\varGamma =2.5$, 
$\rho =1.0$. The figure clearly signifies the effect of the magnetic field on the diffusion coefficient. In the zero 
$\beta$ case, it is observed that the diffusion increases with 
$E_{0}$ and reaches the maximum of 
$D_{\textrm {xL}}$, the corresponding 
$E_{c}$ is found as 
$50$ for the 
$\beta =0$ case, where the system begins to form the lanes as indicated by the value of the order parameter 
$\phi =0.5$. Now the lateral transport decreases strongly once the particles are in a lane, which reduces the perpendicular diffusion again. As a result, drastic decrease of 
$D_{\textrm {xL}}$ with 
$E_{0}$ takes place beyond 
$E_{c}$, which strongly indicates the location of the lane formation. However, as indicated in figure 3, unlike the zero 
$\beta$ case, for finite 
$\beta$ case the value of 
$D_{\textrm {xL}}$ goes on increasing along with the electric field strength 
$E_{0}$, beyond 
$E_{c}$. This increase in transport perpendicular to the both electric and magnetic field strength indicates the presence of the 
$E\times B$ drift effect. This may be explained as follows: as the PIP particles tend to move parallel in the direction of the electric field (
$y$-direction) with a certain velocity, they experience a (
$E \times B$) drift velocity of magnitude 
$v_{E}={E_{0}}/{\beta }$ (in MD units) due to the presence of both electric field and magnetic field, and as this drift velocity is directly proportional to the electric field strength 
$E_{0}$ and inversely proportional to 
$\beta$, therefore, for a fixed value of 
$\beta$ with increasing 
$E_{0}$, drift velocity increases rapidly, which in-turn increases the magnitude of the long-time self-diffusion coefficient as shown in figure 3. However, by modifying the calculation of the self-diffusion coefficient and removing the effect of the 
$E \times B$ drift velocity (which is not done here), it may be possible to get a similar type of 
$D_{\textrm {xL}}$ versus 
$E_{0}$ plot as observed for the 
$\beta =0$ case.
Perpendicular diffusion coefficient (
$D_{\textrm {xL}}$) versus 
$E_{0}$ plot for constant external electric field case (
$\omega =0$), both in presence (
$\beta =1.0$) and absence (
$\beta =0$) of the external magnetic field, taking 
$\varGamma =2.5$, 
$\rho =1.0$.
The phase diagram obtained from a simulation with a constant (
$\omega =0$) external force is shown in figure 4. The location of phase transition is identified as the ‘critical field strength 
$E_{c}$’ which is obtained by setting 
$\phi =0.5$ for a set of runs with increasing 
$E_{0}$. This exercise is performed for different values of the density 
$\rho =0.5$, 
$1.0$ and external magnetic field strength (
$\beta$). An interesting phenomenon observed in this study is that the critical value of electric field strength 
$E_{c}$ increases with the application of the magnetic field (
$\beta$), which indicates that the strength of the external magnetic field also plays an important role in lane formation dynamics. The higher the magnetic field strength, the higher the strength of external electric field required to reach the transition point. From the figure it is also observed that a higher value of the electric field strength (
$E_{c}$) is required by the system to enter into the lane formation state with an increasing value of the number density of the system in the presence of the magnetic field.
Critical electric field 
$E_{c}$ versus 
$\beta$ plot for constant electric field case (
$\omega =0$), taking 
$\rho =0.5$, 
$1.0$, 
$\varGamma =2.5$.
A study related to the dependence of the Coulomb coupling parameter (
$\varGamma$) on lane formation dynamics is investigated in the presence of the external magnetic field 
$\beta \neq 0$ in the constant electric field case. Data for the order parameter 
$\phi$ which decides the location of phase transition are plotted in figure 5 in the presence of the different values of the external magnetic field (
$\beta =0.06, 1.0$). They are shown versus the Coulomb coupling parameter (
$\varGamma$). The other parameters used are 
$E_{0}=150$, 
$\omega =0$ and 
$\rho =1.0$. The plots clearly show that in the beginning for smaller 
$\varGamma$ values, 
$\phi$ is small, it grows and reaches 
$\phi =0.5$ with increasing 
$\varGamma$, from where the formation of lanes start in the PIP system and that value of 
$\varGamma$ is denoted as 
$\varGamma _{\textrm {lc}}$, the lower limit of critical 
$\varGamma$, which is found as 
$1.9$ for 
$\beta =1.0$, and for 
$\beta =0.06$, 
$\varGamma _{\textrm {lc}}$ is found at 
$1.8$. Further increasing the 
$\varGamma$ value beyond 
$\varGamma _{\textrm {lc}}$, the 
$\phi$ value grows and gradually reaches unity. This increase in 
$\phi$ values gradually reduces and reaches 
$\phi =0.5$ again for a 
$\varGamma$ value on the higher side, this 
$\varGamma$ is indicated as the higher limit of critical 
$\varGamma$, denoted by 
$\varGamma _{\textrm {hc}}$. However, for further increase in the 
$\varGamma$ value beyond 
$\varGamma _{\textrm {hc}}$, the 
$\phi$ value again decreases further. The region between 
$\varGamma _{\textrm {lc}}$ and 
$\varGamma _{\textrm {hc}}$ indicates the lane formation regime (where, 
$0.5 \geqslant \phi \leqslant 1.0$) and this whole trend is observed for both the cases. From our study it is observed that the lane phase extends to a certain 
$\varGamma$ range, for 
$\beta =1.0$, from 
$\varGamma _{\textrm {lc}}=1.9$ (lower limit of critical 
$\varGamma$) to 
$\varGamma _{\textrm {hc}}=14.2$ (indicates as higher limit of critical 
$\varGamma$), on the other hand for 
$\beta =0.06$ it extends up to 
$\varGamma _{\textrm {hc}}=14.1$. Further increase of the 
$\varGamma$ value yields 
$\phi$ values that decrease gradually as shown in figure 5, indicating the decreasing order of the system. However, at 
$\varGamma =14.0$ near phase transition boundary (see, figure 6a) a sudden rise of the order parameter value (
$\phi =0.6$) at 
$t\omega _{\textrm {pd}}=2400$ is observed indicating a lane ordered state (this can also be understood from the time variation of the 
$\phi$ plot (figure 7a), along with figure 7b for 
$\beta =1.0$, the corresponding times when the instantaneous positions of particles are represented), which happens possibly due to the reason that, near phase boundary there is some discontinuity in free energy takes place. Later, at some points beyond 
$\varGamma _{\textrm {hc}}$, the system again transits to a completely disordered phase (see 
$\phi$ versus 
$t\omega _{\textrm {pd}}$ plots along with the instantaneous particles snapshots in figures 9, 10, respectively, for 
$\varGamma =14.1$, 
$15$). In addition to this, at 
$\varGamma =14.2$ for 
$\beta =1.0$, the system shows a very peculiar behaviour, in which it again suddenly transfers to the lane state with 
$\phi =0.82$ (see figure 11; the physics is not understood), however, beyond that point the system again transits to a disordered state with 
$\varGamma =15$ (see figure 10). The behaviour of the 2-D PIP system during the 
$\varGamma$ scan in the presence of the different values of the external magnetic field (
$\beta =0.06,1.0$) may be recognized as the lane formation-disintegration phase. In figure 8, we represent a sequence of snapshots that illustrate this evolution for 
$\beta =1.0$.
Phase diagram for constant electric field case (
$\omega =0$) in the presence of the external magnetic field; location of the phase diagram is estimated via the behaviour of the order parameter for various values of the 
$\varGamma$ in the presence of the different magnetic field (
$\beta =0.06, 1.0$). Other parameters are 
$E_{0}=150$, 
$\rho =1.0$, 
$\varGamma _{\textrm {lc}}$ and 
$\varGamma _{\textrm {hc}}$, respectively, and indicate lower value of critical limit and higher value of critical limit of 
$\varGamma$. For 
$\beta =0.06$, 
$1.0$, 
$\varGamma _{\textrm {lc}}$ values obtained from simulation are 
$1.8$, 
$1.9$, respectively, and 
$\varGamma _{\textrm {hc}}$ values are 
$14.1$, 
$14.2$, respectively.
Phase diagram indicating ordered lane state and disordered no-lane state. (a) Phase diagram for constant electric field case (
$\omega =0$) both in the presence (
$\beta =1.0$) and absence (
$\beta =0$) of the external magnetic field; location of the phase diagram is estimated via the behaviour of the order parameter for various values of the 
$\varGamma$. Other parameters are 
$E_{0}=150$, 
$\rho =1.0$, 
$\varGamma _{\textrm {lc}}$ and 
$\varGamma _{\textrm {hc}}$, respectively, which indicate lower value of critical limit and higher value of critical limit of 
$\varGamma$. For 
$\beta =0$, 
$\varGamma _{\textrm {lc}}$ and 
$\varGamma _{\textrm {hc}}$ are found as 
$1.9$, 
$14.1$, respectively, and for 
$\beta =1.0$ 
$\varGamma _{\textrm {lc}}=1.9$ and 
$\varGamma _{\textrm {hc}}=14.2$ are found. (b) Phase diagram for constant electric field case (
$\omega =0$) in the presence of the external magnetic field (
$\beta =0.06$) for both 
$\kappa =0.0001$, 
$4.0$; location of the phase diagram is estimated via the behaviour of the order parameter for various values of the 
$\varGamma$. Other parameters are 
$E_{0}=150$, 
$\rho =1.0$, simulated values of 
$\varGamma _{\textrm {lc}}=1.1$ and 
$\varGamma _{\textrm {hc}}=3.5$, respectively, for 
$\kappa =4.0$, and for 
$\kappa =0.0001$ 
$\varGamma _{\textrm {lc}}=1.8$ and 
$\varGamma _{\textrm {hc}}=14.1$ are observed. Here 
$\varGamma _{\textrm {lc}}$ indicates the lower value of critical limit and 
$\varGamma _{\textrm {hv}}$ indicates the higher value of critical limit of 
$\varGamma$.
Lane formation dynamics both in the presence and absence of the external magnetic field for constant electric field cases. (a) Order parameter 
$\phi$ versus 
$t\omega _{\textrm {pd}}$ plot for constant field case, both in the presence (
$\beta =1.0$) and absence (
$\beta =0.0$) of the external magnetic field. Here, electric field is applied at 
$t\omega _{\textrm {pd}}=600$ time step. The vertical dashed lines indicate the time at which the simulation snapshots are collected and plotted, as shown in figure 7(b). The other parameters used are 
$E_{0}=150$, 
$\varGamma =14.0$ and 
$\rho =1.0$. (b) Typical simulation snapshots taken at 
$t\omega _{\textrm {pd}}=300, 600, 2400$, showing lane formation dynamics in the constant electric field case, both in the presence (
$\beta =1.0$) and absence (
$\beta =0.0$) of the external magnetic field. The other parameters used are 
$E_{0}=150$, 
$\varGamma =14.0$ (near phase transition location) and 
$\rho =1.0$.
Instantaneous positions of particles for different values of 
$\varGamma$ for constant electric field (
$\omega =0$) case in the presence of the external magnetic field 
$\beta =1.0$, the other parameters used are 
$E_{0}=150$, 
$\rho =1.0$.
Lane formation dynamics both in the presence and absence of the external magnetic field for the constant electric field case. (a) Order parameter 
$\phi$ versus 
$t\omega _{\textrm {pd}}$ plot for the constant field case, both in the presence (
$\beta =1.0$) and absence (
$\beta =0.0$) of the external magnetic field. Here, electric field is applied at 
$t\omega _{\textrm {pd}}=600$ time step. The vertical dashed lines indicate the time at which the simulation snapshots are collected and plotted as shown in figure 9(b). The other parameters used are 
$E_{0}=150$, 
$\varGamma =14.1$ and 
$\rho =1.0$. (b) Typical simulation snapshots taken at 
$t\omega _{\textrm {pd}}=300$, 600, 2400, showing lane formation dynamics in the constant electric field case, both in the presence (
$\beta =1.0$) and absence (
$\beta =0.0$) of the external magnetic field. The other parameters used are 
$E_{0}=150$, 
$\varGamma =14.1$ and 
$\rho =1.0$.
Lane formation dynamics both in the presence and absence of the external magnetic field for constant electric field case. (a) Order parameter 
$\phi$ versus 
$t\omega _{\textrm {pd}}$ plot for constant field case, both in the presence (
$\beta =1.0$) and absence (
$\beta =0.0$) of the external magnetic field. Here, electric field is applied at 
$t\omega _{\textrm {pd}}=600$ time step. The vertical dashed lines indicate the time at which the simulation snapshots are collected and plotted as shown in figure 10(b). The other parameters used are 
$E_{0}=150$, 
$\varGamma =15.0$ and 
$\rho =1.0$. (b) Typical simulation snapshots taken at 
$t\omega _{\textrm {pd}}=300$, 600, 2400, showing lane formation dynamics in the constant electric field case, both in the presence (
$\beta =1.0$) and absence (
$\beta =0.0$) of the external magnetic field. The other parameters used are 
$E_{0}=150$, 
$\varGamma =15.0$ and 
$\rho =1.0$.
Lane formation dynamics both in the presence and absence of the external magnetic field for the constant electric field case. (a) Order parameter 
$\phi$ versus 
$t\omega _{\textrm {pd}}$ plot for the constant field case, both in the presence (
$\beta =1.0$) and absence (
$\beta =0.0$) of the external magnetic field. Here, electric field is applied at 
$t\omega _{\textrm {pd}}=600$ time step. The vertical dash lines indicate the time at which the simulation snapshots are collected and plotted as shown in figure 11(b) below. The other parameters used are 
$E_{0}=150$, 
$\varGamma =14.2$ (phase transition location, see figure 6a) and 
$\rho =1.0$. (b) Typical simulation snapshots taken at 
$t\omega _{\textrm {pd}}=300$, 600, 2400, showing lane formation dynamics in the constant electric field case, both in the presence (
$\beta =1.0$) and absence (
$\beta =0.0$) of the external magnetic field. The other parameters used are 
$E_{0}=150$, 
$\varGamma =14.2$ and 
$\rho =1.0$.
The lane formation–disintegration phase has also been recognized in our PIP system in the absence of the external magnetic field as shown in figure 6(a), where the 
$\varGamma$ variation of the order parameter (
$\phi$) values in the absence of the external magnetic field (
$\beta =0.0$) are plotted with that of the results obtained in the presence of the external magnetic field. For the 
$\beta =0$ case, 
$\varGamma _{\textrm {lc}}$ and 
$\varGamma _{\textrm {hc}}$ are observed as 
$1.9$, 
$14.1$, respectively, and for 
$\beta =1.0$, 
$\varGamma _{\textrm {lc}}$ and 
$\varGamma _{\textrm {hc}}$ read as 
$1.9$ and 
$14.2$, respectively. The lane formation–disintegration phase diagram is also plotted for the finite 
$\kappa$ case (
$\kappa =4.0$), as shown in figure 6(b), in the presence of the external magnetic field strength 
$\beta =0.06$ and compared the same with that of the zero 
$\kappa$ case in order to see the effect of screening on the lane disintegration process. The finite 
$\kappa$ case is studied only for comparison, as the potential used in this study is not accurate for the finite 
$\kappa$ case. From our study it is observed that for the finite 
$\kappa$ case, the lane formation area shrinks to a smaller region as compared with that of the zero 
$\kappa$ case, which indicates that the screening parameter has a significant effect on the lane formation–disintegration process. The lane area extends from 
$\varGamma _{\textrm {lc}}=1.1$ to 
$\varGamma _{\textrm {hc}}=3.5$ for 
$\kappa =4.0$, and for 
$\kappa =0.0001$, it ranges from 
$\varGamma _{\textrm {lc}}=1.8$ to 
$\varGamma _{\textrm {hc}}=14.1$. Additionally, from our study it is also been observed that in the absence of the external magnetic field (
$\beta =0$), at 
$\varGamma =14.0$ and 
$\varGamma =14.1$, as shown in figures 7(a) and 9(a), respectively, the order parameter value suddenly raises beyond 
$\phi =0.5$ at 
$t\omega _{\textrm {pd}}=1300, 1000$ time steps, respectively, although, the external electric field is applied at 
$t\omega _{\textrm {pd}}=600$ time step. However, this behaviour of the order parameter plot in the 
$\beta =0$ case is not clearly understood. In addition to that, for both 
$\beta =0$ and 
$\beta \neq 0$ cases, the physical reason behind this lane disintegration process has yet to be understood.
As shown in figures 12(a) and 12(b), a comparative study on time variation of the order parameter (
$\phi$) is performed between the 
$\beta =0$ case and the 
$\beta \neq 0$ case (taking different values of the magnetic field strength 
$\beta =0.001$, 0.005, 0.01, 0.06) taking 
$\varGamma =2.5$ (near the lower limit of the critical 
$\varGamma$ value) and 
$\varGamma =5.0$ (far from the 
$\varGamma _{\textrm {lc}}$), respectively. Here 
$\varGamma =2.5$ taken for figure 12(a) is just above the lower value of the critical 
$\varGamma$ limit (i.e. 
$\varGamma _{\textrm {lc}}=1.9$) both in the presence and absence of the external magnetic field (see, figures 5 and 6a) and 
$\varGamma =5.0$ used for figure 12(b) is far beyond the critical lower limit, (see figures 5 and 6a). For all the cases the other input parameters are taken as 
$E_{0}=150$, 
$\omega =0.0$ and 
$\rho =1.0$. Here, the magnetic field is applied from the beginning of the simulation, i.e. at 
$t\omega _{\textrm {pd}}=0$ time steps, and electric field at 
$t\omega _{\textrm {pd}}=600$ time steps. From the 
$\phi$ versus 
$t\omega _{\textrm {pd}}$ plot as shown in figure 12(b), taking 
$\varGamma =5.0$, an appearance of an oscillatory pattern of the order parameter plots is observed in the presence of the external magnetic field and it increases with an increase in 
$\beta$ value. From this plot it is also observed that the frequency of the corresponding oscillation of the order parameter plot increases with an increase in the 
$\beta$ value. This result is supported by the frequency (
$f$) of oscillation of the order parameter 
$\phi$ versus 
$\beta$ plot as shown in figure 12(c). This is seen for the whole range of the magnetic field strength under study for 
$\varGamma =5.0$. Using a linear fit method, a relation may be obtained between 
$f$ and 
$\beta$ as 
$f=0.0076973 + 0.75581 \beta$. However, this oscillatory pattern of the time varying order parameter plot is not visible for lower 
$\varGamma$ value (i.e. 
$\varGamma =2.5$, see figure 12a) even with increasing 
$\beta$ value. Additionally, in the lower 
$\varGamma ({=}2.5)$ case, the order parameter varies in the similar way with time both in the presence and absence of the magnetic field, which is not visible in the presence of a higher 
$\varGamma ({=}5.0)$ value.
(a) Variation of order parameter (
$\phi$) versus 
$t\omega _{\textrm {pd}}$ for varying magnetic field strength 
$\beta =0.001, 0.005, 0.01, 0.06$ compared with that of 
$\beta =0$ case taking 
$\varGamma =2.5$. The other parameters used are 
$E_{0}=150$, 
$\rho =1.0$, 
$\omega =0.0$. (b) Variation of order parameter (
$\phi$) versus 
$t\omega _{\textrm {pd}}$ for varying magnetic field strength 
$\beta =0.001$, 0.005, 0.01, 0.06 compared with that of 
$\beta =0$ case taking 
$\varGamma =5.0$. The figures in the second column are zoomed plots from 
$t\omega _{\textrm {pd}}=600$ to 
$1000$. The other parameters used are 
$E_{0}=150$, 
$\rho =1.0$, 
$\omega =0.0$. (c) The 
$\beta$ dependence of the frequency of oscillation (
$f$) of the order parameter (
$\phi$); the dashed black line is the linear fitting and the blue solid line locates 
$f = \beta$. The other parameters used are 
$\varGamma =5.0$, 
$E_{0}=150$, 
$\rho =1.0$, 
$\omega =0.0$.
4.1.2 Oscillatory field
In this subsection, the lane formation dynamics of 2-D PIP in the presence of an external oscillating electric field (
$\omega \neq 0$) by simulating Langevin molecular dynamics is presented. In particular, we have examined the effect of external magnetic field (
$\beta \neq 0$) on lane formation. Here, the external magnetic field is applied from the beginning of the simulation, and external electric field is applied from 
$t\omega _{\textrm {pd}}=600$ time step in the case under study.
In this subsection, in figure 13(a), the simulation snapshots associated with a situation are plotted at different simulation times for an external electric field having sinusoidal oscillation and of strength 
$E_{0}=150$, frequency 
$\omega =0.001$, 
$\varGamma =2.6$, 
$\rho =0.5$ and in the presence of the external magnetic field with 
$\beta =1.0$. The initial snapshot at 
$t\omega _{\textrm {pd}}=300$ with 
$\phi =0.037$ corresponds a disordered phase. This snapshot is taken just before the external electric field is switched on at 
$t\omega _{\textrm {pd}}=600$. The second snapshot is taken at 
$\omega =910$ having 
$\phi =0.71$, when the system reaches to a lane formation state. The next snapshot is taken at 
$t\omega _{\textrm {pd}}=1634$ with 
$\phi =0.125$ showing the disordered state and at 
$t\omega _{\textrm {pd}}=2242$, the system again reaches to ordered lane state with 
$\phi =0.791$, as shown in the last snapshot of figure 13(a). This whole process of spontaneous formation and breaking of lane structures is observed to repeat during our extensive long-time simulations in the presence of external magnetic field in the lower electric frequency range, like in our previous study in the absence of the external magnetic field
$^{6}$. This result is also supported by the time variation of the lane order parameter plot as shown in figure 13(b). From figure 13(b), it is also observed that for the frequency 
$\omega =0.001$, the order parameter (
$\phi$) value decreases beyond the critical value of 
$\phi =0.5$ at some points of simulation, which indicates that the system has a tendency to transit from an ordered lane state to a disordered mixed phase. The lane structure (with 
$\phi \approx 1$) formed at some special moments breaks up spontaneously, but the resulting disorganized state (
$\phi \approx 0$) is unstable, and the lane formation immediately proceeds again, resulting in the spikes of 
$\phi$ in figure 13(b). This phase transition observed here in figure 13(a,b) in the presence of the oscillatory external electric field and uniform external magnetic field can be understood in the following way. In the segregated mixture, the PIP particles are moving collectively with the external electric field. Due to the oscillatory external electric field the velocity of the particles in field direction changes sign but roughly has the modulus of the ideal drift velocity. At the interfaces between two lanes, there is an additional friction due to the oppositely moving particles of the other type. Hence, the drift velocity near an interface changes sign periodically according to the shaking external field. Now, for the stability of two lanes at their interface, the field frequency (
$\omega$) has to be small enough in order to provide a sufficiently long time period in which the two lanes can slide against each other avoiding a mixing of different particle species, and lanes are not destroyed. If this time is getting small, diffusion perpendicular to the field direction will dominate and destroy the sharp interface, and this can be observed from figure 13(c) where time variation of 
$\phi$ is plotted for larger, intermediate and lower frequency values. From figure 13(c), it is also observed that 
$\phi$ oscillation tends to disappear for the intermediate frequency value (
$\omega =0.7$) before the plasma loses its self-organization for the larger frequency value (
$\omega =7.0$). For intermediate frequency values (
$\omega =0.7$), the particles get comparatively less chance for interaction with different kind of particles, this restrains effectively the random motion of the particles, as a result, the particles are not likely to get entangled into another lane and lanes are not destroyed. Moreover, like the constant electric field with finite 
$\beta$ case, here, in the oscillatory electric field with 
$\beta =1.0$ case also the lanes seem to move along the 
$X$-direction of the simulation chamber, which indicates the (
$E\times B$) drift effect.
Lane formation dynamics in the presence of the external magnetic field (with 
$\beta =1.0$) for oscillatory electric field case (
$\omega \neq 0$). (a) Typical simulation snapshots taken at 
$t\omega _{\textrm {pd}}=300$, 910, 1634, 2242, showing spontaneous lane formation and deformation dynamics in the oscillatory electric field case (
$\omega =0.001$), in the presence of external magnetic field (
$\beta =1.0$). The other parameters used are 
$E_{0}=150$, 
$\varGamma =2.5$, 
$\rho =0.5$. (b) Time variation of order parameter (
$\phi$) plot taking the external oscillatory electric field frequency 
$\omega =0.001$ in the presence of external magnetic field (
$\beta =1.0$). The vertical dashed lines indicate the time at which instantaneous position of the PIP particles (see, figure 13c) are recorded. The other parameters used are 
$E_{0}=150$, 
$\varGamma =2.5$, 
$\rho =0.5$. (c) Time variation of order parameter (
$\phi$) plot for 
$\omega =0.001$, 0.7, 3.5, 7.0 in the presence of external magnetic field (
$\beta \neq 0$). The other parameters used are 
$E_{0}=150$, 
$\varGamma =2.5$, 
$\rho =0.5$ and 
$\beta =1.0$.
As shown in figure 13(c), for the oscillatory external electric field case, from the order parameter 
$\phi$ versus 
$t\omega _{\textrm {pd}}$ plot for various frequency values, the critical frequency (
$\omega = \omega _c$) upon which a transition of the system back to the disordered state occurs can be understood, if an oscillating field with amplitude 
$E_{0} > E_{c}$ is present. Like our previous study (Sarma et al. Reference Sarma, Baruah and Ganesh2020) in the absence of the external magnetic field, here, in the presence of the external magnetic field also, it is evident from figure 13(c) that lane formation with the oscillating sinusoidal electric field depends on the value of 
$\omega$. If an oscillating field with amplitude 
$E_{0} > E_{c}$ is present and the frequency 
$\omega$ is increased, the order parameter 
$\phi$ decreases with increasing 
$\omega$. For low frequencies the system remains in the lane state, while above a critical frequency (
$\omega _{c}$), the system goes back to the disordered state. In the present case, 
$\omega _c=3.5$ where 
$\phi \approx 0.5$. Above 
$\omega _c=3.5$ the system remains in the disordered state with 
$\phi <0.5$ as shown in figure 13(c). With this concept, the phase diagram for an underdamped PIP system with an oscillating sinusoidal external electric field and in the presence of an external uniform magnetic field taking density 
$\rho =1.0$ is plotted and is shown in figure 14(a). The same is plotted in figure 14(b) for different values of densities (
$\rho =0.5, 1.0$) taking into consideration both the presence (
$\beta =1.0$) and absence (
$\beta =0.0$) of the external magnetic field. In both the cases (figure 14a,b), the critical external oscillating field frequency 
$\omega _{c}$ is plotted versus 
$E_{0}$. The trend of the 
$\omega _{c}$ increases with increasing 
$E_{0}$, as can be seen in both figures 14(a) and 14(b). From figure 14(b), it is also observed that the trend of 
$\omega _{c}$ increases with increasing density 
$\rho$. In figure 14(b), the 
$\omega _{c}$ values obtained in the absence of the external magnetic field (
$\beta =0.0$), for density values 
$\rho =0.5, 1.0$, are plotted along with the results obtained in the presence of the external magnetic field 
$\beta \neq 0.0$, and clearly it is observed that upon application of magnetic field, a lower value of critical field frequency is required to transit the system from the ordered lane state to the disordered state.
Non-equilibrium phase diagram for 2-D PIP system with oscillatory external electric field in the application of the external magnetic field. (a) Non-equilibrium phase diagram with oscillatory external electric field (
$\omega \neq 0$) in the presence of the external magnetic field (
$\beta =1.0$) showing lane zone and disordered zone. The other parameters used are 
$\varGamma =2.5$, 
$\rho =0.5$. (b) Non-equilibrium phase diagram with oscillatory external electric field (
$\omega \neq 0$), both in the presence (
$\beta =1.0$) and absence (
$\beta =0$) of the external magnetic field for different 
$\rho$ values. The other parameters used are 
$\varGamma =2.5$, 
$E_{0}=2.5$.
In figure 15, time variation of the order parameter 
$\phi$ is plotted for different values of magnetic field, ranging from 
$\beta =0.0$ to 
$\beta =1.0$ (although few are plotted here) taking input parameters as 
$E_{0}=150$, 
$\omega =0.001$, 
$\varGamma =2.5$ and 
$\rho =1.0$. The figure shows that the varying magnetic field strength does not bring any significant change in the trend of order parameter 
$\phi$ with time for a lower 
$\varGamma =2.5$ value. This can also be understood from comparative time variation of order parameter plots both in the presence and absence of the external magnetic field, as shown in figure 16(a). Here, the input parameters are the same as used in the case of figure 15. In figure 16(b), a similar type of time variation of order parameter data are plotted, but, taking 
$\varGamma =5.0$, and it is observed that with an increase in 
$\beta$ value, an oscillatory pattern of the order parameter plots appears and the frequency of corresponding oscillation increases with an increase in 
$\beta$ value. This is seen for the whole range of magnetic field strength under study for 
$\varGamma =5.0$.
Time variation of order parameter (
$\phi$) plot for various values of the external magnetic field strength 
$\beta =0$, 0.005, 0.01, 0.06, 0.1, 0.5, 1.0, in the presence of the oscillatory (
$\omega \neq 0$) external electric field having strength 
$E_{0}=150$. The other parameters used are 
$\varGamma =2.5$, 
$\rho =1.0$.
Comparative time variation of order parameter (
$\phi$) plot for two different values of Coulomb coupling parameter (
$\varGamma =2.5, 5.0$), both in the absence (
$\beta =0$) and presence (
$\beta \neq 0$) of the external magnetic field. (a) Time variation of order parameter (
$\phi$) plot for the external magnetic field strength values 
$\beta =0$, 0.005, 0.01, 0.06, for a PIP system in the presence of the oscillatory external electric field having frequency (
$\omega =0.001$). The Coulomb coupling parameter value used here is 
$\varGamma =2.5$; the other parameters used are 
$E_{0}=150$, 
$\rho =1.0$. (b) Time variation of order parameter (
$\phi$) plot for the external magnetic field strength values 
$\beta =0$, 0.005, 0.01, 0.06, for a PIP system in the presence of the oscillatory external electric field having frequency (
$\omega =0.001$). The figures in the second column are zoomed plots from 
$t\omega _{\textrm {pd}}=600$ to 
$1000$. The plots are showing an oscillatory pattern of 
$\phi$, which increases with increasing 
$\beta$, for clarity, zoomed plots are mentioned on the right-hand side. The parameters used are 
$\varGamma =5.0$, 
$E_{0}=150$, 
$\rho =1.0$.
For an oscillating external electric field case, to study the effect of external magnetic field on critical frequencies (
$\omega _{c}$), data for the 
$\omega _{c}$ are also plotted versus 
$\beta$, as shown in figure 17. The study is performed by taking a range of 
$\beta$ values starting from 
$\beta =0.001, 0.005$ to 
$\beta =1.0$, along with considering the case where external magnetic field is absent (
$\beta =0.0$), taking 
$E_{0}=150$, 
$\varGamma =2.5$ and 
$\rho =1.0$. The critical frequency 
$\omega _{c}$ is obtained by setting 
$\phi =0.5$, for the set of runs with increasing frequency 
$\omega$, corresponding to all the values of 
$\beta$. As can be identified from the figure, the 
$\omega _{c}$ value corresponding to 
$\beta =0.0$ is 
$7.7$, however, for 
$\beta =0.001, 0.005$, the corresponding 
$\omega _{c}$ values are obtained as 
$w_{c}=7.6, 7.0$, respectively. It is observed that for a further increase in 
$\beta$ values to 
$\beta =0.007$, the 
$\omega _{c}$ value decreases further to 
$6.5$. However, from 
$\beta =0.01$, it starts to increase gradually and for the rest of the higher 
$\beta$ values it maintains the same pattern and finally, as can be seen, for 
$\beta =1.0$, 
$\omega _{c}$ is obtained as 
$8.1$. Here, one sees that a critical value of the electric field frequency below which self-organization takes place also exists and the critical frequency is an increasing function of higher 
$\beta$ values.
4.2 Lane formation in the presence of oblique electric field
In this section, we have presented the lane formation dynamics of PIP by choosing the external electric force to be operative in an oblique way by making an angle 
$\theta$ with the 
$X$–
$Y$ plane of the simulation box. The study is performed both in the presence and absence of the external magnetic field. In particular, we have examined the effect of the external magnetic field on lane formation. The effect of system asymmetry on lane formation has also been studied here.
Variation of the critical value of the electric field frequency (
$\omega _{c}$) with 
$\beta$ plot. Input parameters used are 
$E_{0}=150$, 
$\varGamma =2.5$, 
$\rho =1.0$, 
$\Delta t=0.003$.
Instantaneous position of particles in the application of constant (
$\omega =0$) oblique electric field taking 
$\theta =30^{\circ }, 45^{\circ }, 60^{\circ }, 90^{\circ }$ (
$\theta$ is the angle made by the electric field with the 
$x$-direction), external electric force 
$\boldsymbol {F}_{E}$ acting on the particles pointing in the 
$Y$-direction; other parameters used are 
$E_{0}=170$, 
$\varGamma =2.5$, 
$\rho =1.0$. This study is performed both in the presence (
$\beta =1$) and absence (
$\beta =0$) of the external magnetic field.
4.2.1 Symmetric system
In this subsection, the whole study has been performed by using a regular square shaped simulation chamber (
$X=Y=L$, a symmetric system). Here, first, we will discuss the results obtained from a study performed for a constant electric field (
$\omega =0$) case which is applied obliquely for various values of 
$\theta$, in the presence (
$\beta =1.0$) of the external magnetic field. Second, the results obtained from the obliquely applied oscillatory electric field case (
$\omega \neq 0, \theta \neq 0$), in the presence of the external magnetic field have also been discussed. For both the above mentioned cases, a comparative study on lane formation dynamics is performed both in the presence and absence of the external magnetic field and results are mentioned below.
In figure 18, the instantaneous positions of the particles are plotted at different values of 
$\theta$, ranging from 
$\theta =30^{\circ }, 45^{\circ }$ to 
$\theta =225^{\circ }$ taking constant electric field case with 
$\beta =1.0$. The other parameters used are 
$E_{0}=170$, 
$\omega =0$, 
$\varGamma =2.5$ and 
$\rho =1.0$. In this case, the forcing due to the oblique electric field is presented throughout the simulation duration after 
$t\omega _{\textrm {pd}}=600$ time steps. In figure 18, from the snapshot obtained for 
$\beta =1.0$, 
$\theta =30^{\circ }$, one can clearly observe lane structures where the direction of the lanes coincides with the 
$X$-axis of the simulation box. For 
$\theta =45^{\circ }$ with 
$\beta =1.0$, some steady lane states are observed with tilted bands, see the snapshot in figure 18. These tilted lanes seen in the simulations might be artefacts of the external electric field used obliquely (
$\theta =45^{\circ }$) in the simulation. With further increasing the oblique angle to 
$\theta =60^{\circ }$, lane formation occurs where the direction of the lanes coincides with the 
$Y$-direction of the simulation chamber. Lastly with 
$\theta =90^{\circ }$, the typical lane structures are formed parallel to the drive direction, see the corresponding snapshot in figure 18. To verify whether the tilted bands that are observed for 
$\theta =45^{\circ }$ (figure 18) will be reformed, several simulations are performed for different 
$\theta$ values, and for 
$\theta =135^{\circ }$ and 
$225^{\circ }$, reappearance of the similar kind of tilted bands are observed for both in the presence and absence of the external magnetic field, the related particle snapshots for the constant electric field case are plotted, as shown in figure 19.
Instantaneous position of particles in the application of constant (
$\omega =0$) oblique electric field taking 
$\theta =135^{\circ }, 225^{\circ }$ (
$\theta$ is the angle made by the electric field with 
$x$-direction), external electric force 
$\boldsymbol {F}_{E}$ acting on the particles pointing in the 
$Y$-direction; other parameters used are 
$E_{0}=170$, 
$\varGamma =2.5$, 
$\rho =1.0$. This study is performed both in the presence (
$\beta =1$) and absence (
$\beta =0$) of the external magnetic field.
From our study similar types of pattern formations are also observed in the case of the oscillatory electric field having frequency of oscillation 
$\omega = 2.0$, applied obliquely making the same magnitude of the oblique angles (
$\theta =30^{\circ }, 45^{\circ }, 60^{\circ }, 90^{\circ }$) with the 
$X$–
$Y$ plane of the simulation box as used in the constant field case, in the presence of the external magnetic field. The related particle snapshots are plotted as shown in figure 20. The other parameters used for this study are 
$E_{0}=170$, 
$\varGamma =2.5$ and 
$\rho =1.0$. From our study, a similar trend of the lane structures is also observed in the absence of the external magnetic field (
$\beta =0$), as it is obtained in the presence of both the obliquely applied constant (see figure 18 for the 
$\omega =0$ case) and oscillatory electric fields (see figure 19 for the 
$\omega \neq 0$ case). The only difference observed between the results obtained from the oblique force field cases in the presence and absence of the external magnetic field is that, in the presence of the external magnetic field (
$\beta =1.0$), the movement of the lanes is observed in the plane (
$X$–
$Y$) perpendicular to the direction of the applied external magnetic field (
$Z$-direction), which is not observed in the absence of the magnetic field. This movement of the lanes in the 
$X$–
$Y$ plane in the presence of the external magnetic field indicates the presence of the (
$E\times B$) drift in the system.
Instantaneous position of particles in the application of oscillatory (
$\omega \neq 0$) oblique electric field taking 
$\theta =30^{\circ }, 45^{\circ }, 60^{\circ }, 90^{\circ }$ (
$\theta$ is the angle made by the electric field with 
$x$-direction), external electric force 
$\boldsymbol {F}_{E}$ acting on the particles pointing in the 
$Y$-direction; other parameters used are 
$E_{0}=170$, 
$\omega =2.0$, 
$\varGamma =2.5$, 
$\rho =1.0$. This study is performed both in the presence (
$\beta =1$) and absence (
$\beta =0$) of the external magnetic field.
4.2.2 Asymmetric system
In this subsection, the whole study has been performed by using an asymmetric simulation chamber, 
$X=20Y$, to observe the effect of the aspect ratio on lane formation dynamics. The aspect ratio is a fundamental quantity that denotes the proportional ratio between the individual system sizes. We find that the lane formation does not depend on either 
$X$ or 
$Y$ separately (
$X$, 
$Y$ represent system lengths), rather, the aspect ratio 
${X}/{Y}$ appears to be important. In this study, both constant (
$\omega =0$) and oscillatory (
$\omega \neq 0$) electric fields are used, which are applied obliquely for various values of 
$\theta$. To see the effect of external magnetic field, here a comparative studies has also been performed both in the presence and absence of the magnetic field.
In figure 21, instantaneous particle positions are plotted for 
$\theta =30^{\circ }$, 
$45^{\circ }$, 
$60^{\circ }$ and 
$90^{\circ }$, keeping all the other input parameters values as same as used in the case of the symmetric system in the presence of the constant electric field (see, figure 18). This present study reveals that lane formation is observed for 
$\theta =30^{\circ }$ only, using an asymmetric simulation chamber. There is no lane formation occurring for higher values of 
$\theta$, i.e. 
$\theta > 30^{\circ }$. However, this dependence on aspect ratio is not understood theoretically. A similar lane formation trend is also observed both in the presence and absence of the external magnetic field. However, for the obliquely applied oscillatory (
$\omega =2.0$) electric field case, the lane formation observed inside the simulation chamber for the oblique angle 
$\theta =60^{\circ }$ only, the lanes are aligned with the 
$X$-direction of the simulation chamber as shown in figure 22. The other parameters used are 
$E_{0}=170$, 
$\varGamma =2.5$ and 
$\rho =1.0$. A similar lane formation trend is also observed in the absence of the external magnetic field by taking into consideration an obliquely applied oscillatory electric field, as shown in figure 22.
Instantaneous position of particles in the application of constant (
$\omega =0$) oblique electric field taking 
$\theta =30^{\circ }, 45^{\circ }, 60^{\circ }, 90^{\circ }$ (
$\theta$ is the angle made by the electric field with 
$x$-direction), external electric force 
$\boldsymbol {F}_{E}$ acting on the particles pointing in the 
$Y$-direction, an asymmetric system (
$X=20Y$) is used here; the other parameters are 
$E_{0}=170$, 
$\omega =0.0$, 
$\varGamma =2.5$, 
$\rho =1.0$. This study is performed both in the presence (
$\beta =1$) and absence (
$\beta =0$) of the external magnetic field.
Instantaneous position of particles in the application of oscillatory (
$\omega \neq 0$) oblique electric field taking 
$\theta =30^{\circ }, 45^{\circ }, 60^{\circ }, 90^{\circ }$ (
$\theta$ is the angle made by the electric field with 
$x$-direction), external electric force 
$\boldsymbol {F}_{E}$ acting on the particles pointing in the 
$Y$-direction, an asymmetric system (
$X=20Y$) is used here; the other parameters are 
$E_{0}=170$, 
$\omega =2.0$, 
$\varGamma =2.5$, 
$\rho =1.0$. This study is performed both in the presence (
$\beta =1$) and absence (
$\beta =0$) of the external magnetic field.
5 Conclusions
We have carried out the LD simulations to study lane formation dynamics in a 2-D PIP medium in the presence of forcing due to an external electric field. More specifically, the authors focused on the effect of an external magnetic field on the self-organization effect in the zero 
$\kappa$ (screening parameter) case, where, the ‘plasma’ is structured as a set of lanes that alternatively contain positive and negative ions. The study is based on a plasma model, where, the ion–ion interaction is described with a screened electrostatic potential, and the ion–neutral background interaction is described using an overall friction force characterized by a damping factor and a zero-average stochastic collisional term that enables the description of the diffusion effect. This model is used to describe the dynamics of a set of 
$N$ ions submitted to a stationary or time-varying electric field in a 2-D domain in the presence of an external magnetic field. This effect is quantitatively characterized in terms of an order parameter 
$\phi$. Although the dynamics of lane formation has already been reported in our earlier work (Sarma et al. Reference Sarma, Baruah and Ganesh2020), where self-organization was investigated only in the case of an external electric field in finite 
$\kappa$ limit. In the present work, we consider a much broader situation where the applied external magnetic field is significant and the obliquely applied external electric field is also taken into consideration along with the constant and oscillatory fields.
The paper brings new and interesting information as far as this self-organization effect is concerned in the presence of an external magnetic field. In particular, we finish with the following list of conclusions.
(i) A critical value of the electric field above which self-organization takes place exists in the presence of the external magnetic field (
$\beta \neq 0$).(ii) This critical value of the electric field strength is larger than the one obtained in the absence of the external magnetic field. The higher the magnetic field strength, the higher the strength of external electric field required to reach the transition point.
(iii) The critical value of field strength increases with an increase of the density (
$\rho$) of the system in the presence of the external magnetic field.(iv) Due to the presence of both uniform electric and magnetic fields, the PIP particles experience a drift of the guiding centre (
$E\times B$ drift) and as a result of that, the self-organized lanes start to move along a direction perpendicular to both the directions of the applied electric and magnetic fields with a velocity $v_{E}=E_{0}/\beta$ (in MD units). It is also observed that the speed of these lane movements increases with an increase of the applied electric field strength in the presence of the uniform external magnetic field.(v) The self-organization is correlated to a strong decrease of the transverse diffusion coefficient (
$D_{\textrm {xL}}$) in the zero $\beta$ case (absence of the magnetic field). Unlike the zero $\beta$ case, for the finite $\beta$ case, the value of $D_{\textrm {xL}}$ goes on increasing along with the electric field strength $E_{0}$, which signifies the presence of the $E\times B$ drift in the system.(vi) The phase diagram obtained from the
$\varGamma$ scan of the order parameter ($\phi$) values, for a PIP system with constant external electric field both in the presence and absence of the external magnetic field shows a lane formation–disintegration phase. In all the cases, for the lower $\varGamma$ value side, a critical Coulomb coupling parameter ($\varGamma _{\textrm {lc}}$) appears below which the PIP system is in a disordered phase with $\phi < 0.5$. The lane formation area extends for a wide rage of $\varGamma$ values starting from $\varGamma _{\textrm {lc}}$ before to disappear to a disordered phase above a critical Coulomb coupling parameter ($\varGamma _{\textrm {hc}}$) value in higher side of $\varGamma$. However, a further increase of the $\varGamma$ value lanes reappeared again in some particular point in the presence of an external magnetic field only.(vii) Our investigation demonstrates that due to the presence of both electric and magnetic fields, the PIP particles experience (
$E \times B$) drift which occurs in the form of fluctuations in the order parameter plot in a strongly coupled regime for both the presence of constant and oscillatory external electric fields. The frequency of these fluctuations vary with varying magnitude of the magnetic field strength.(viii) The PIP shows a non-stationary lane-like self-organization behaviour, with oscillation between self-organized and random states for the small frequency values, when submitted to a time-varying electric field in the presence of external magnetic field also. The lanes are found to move along a direction perpendicular to the both
$E$ and $B$ field directions as a result of the guiding centre drift.(ix) For each electric field, there exists a critical frequency
$\omega _{c}$, above which the system cannot reach the self-organized state.(x) The critical frequency is an increasing function of the electric field in the presence of an external magnetic field.
(xi) Upon application of a magnetic field, a lower value of critical field frequency is required to transit the system from an ordered lane state to a disordered state.
(xii) The stability study of lanes identifies the intermediate frequency value, in the presence of the external magnetic field, where
$\phi$ oscillation tends to disappear before the plasma loses its self-organization (for larger frequency). This intermediate frequency value in the absence of the external magnetic field is much higher than that in the presence of an external magnetic field.
Finally, the study of lane dynamics has been performed using an oblique (oblique angle 
$\theta \neq 0$) electric force on a PIP system both in the presence and absence of the external magnetic field for a range of oblique angles 
$\theta = 30^{\circ }\longleftrightarrow 90^{\circ }$. This application of oblique force reveals the following.
(i) Variation of the oblique angle
$\theta$ plays a significant role for achieving the lane state. For a small oblique angle ($\theta =30^{\circ }$), strikingly, the system first develops structures along the $X$-axis. For intermediate oblique angles, i.e. for $\theta =45^{\circ }, 60^{\circ }$, respectively, segregated stripes of PIP particles tilted or oriented perpendicular to the $X$-direction of the simulation chamber. On the other hand, for $\theta =90^{\circ }$, lanes are formed parallel to the direction of the electric field.(ii) The only difference observed in this study in the presence of the external magnetic field is that, in the presence of the
$B$ field, movement of the lanes appears in a direction perpendicular to both the directions of the $E$ and $B$ fields. This effect is due to the presence of the electric field drift, which is not observed in the absence of the $B$ field.(iii) Lastly, the system aspect ratio
${X}/{Y}$ plays a major role on the lane formation dynamics in the presence of an oblique electric field.
A study of self-organized lane structures for PIP system using LD simulation in the presence of the external magnetic field is reported here. Our simulation results may help to identify and extend the parameter range and precise location of lane formation phase in a 2-D PIP system in the presence of the external magnetic field.
Acknowledgements
This research work is supported by the Board of Research in Nuclear Sciences (BRNS), DAE, project sanctioned no. 39/14/25/2016-BRNS/34428, date 20/01/2017; revalidation no. 39/14/25/2016-BRNS/34014, date 03/04/2017. S.B. and U.S. would like to acknowledge the Institute for Plasma Research (IPR), Bhat, Gandhinagar, for allowing us to use the HPC cluster at IPR.
Editor Roger Blandford thanks the referees for their advice in evaluating this article.
Declaration of interests
The authors report no conflict of interest.
