vikram dharodi | Institute for plasma research (original) (raw)

Papers by vikram dharodi

Research paper thumbnail of Signatures of an energetic charge bunch moving in a plasma

A charge bunch moving in a plasma can excite a variety of linear and nonlinear waves in the form ... more A charge bunch moving in a plasma can excite a variety of linear and nonlinear waves in the form of trailing wakes, fore-wake shocks and precursor solitons. These structures can further interact with the background plasma to create secondary effects that can serve as signatures of the passage of the charge bunch. Using particle-in-cell simulations we investigate in detail the dynamics of a plasma system that is being traversed by an energetic charged ion bunch. Using two different shapes of the charge source, namely, an idealized one dimensional line source and a two dimensional thin rectangular source we examine the differences in the nature of the excited wave structures and their consequent impact on the background plasma. Our simulations reveal interesting features such as the dependence of the precursor speeds on the total charge of the ion bunch, local particle trapping, and energization of the trapped ions in various regions along the traversal path leading to the formation of energetic ion beam-lets. The collective excitations and the signatures in the ambient plasma could prove useful in practical applications such as in ion beam heating of plasmas. They can also help in analysing the trajectories of charged objects like space debris orbiting in the ionosphere.

Research paper thumbnail of Kelvin-Helmholtz instability in a strongly coupled dusty plasma medium

Research paper thumbnail of Turbulence in two dimensional visco - elastic medium

arXiv: Plasma Physics, 2014

The properties of decaying turbulence is studied with the help of a Generalized Hydrodynamic (GHD... more The properties of decaying turbulence is studied with the help of a Generalized Hydrodynamic (GHD) fluid model in the context of two dimensional visco - elastic medium such as a strongly coupled dusty plasma system. For the incompressible case considered here however, the observations are valid for a wider class of visco - elastic systems not necessarily associated with plasmas only. Our observations show that an initial spectrum that is confined in a limited domain of wave numbers becomes broad, even when the Reynold’s number is much less than the critical value required for the onset of turbulence in Newtonian fluids. This is a signature of elastic turbulence where Weissenberg’s number also plays a role in the onset of turbulence. This has been reported in several experiments. It is also shown that the existence of memory relaxation time parameter and the transverse shear wave inhibit the normal process (for 2-D systems) of inverse spectral cascade in this case. A detailed simulat...

Research paper thumbnail of A numerical study of gravity-driven instability in strongly coupled dusty plasma. Part 1. Rayleigh–Taylor instability and buoyancy-driven instability

Journal of Plasma Physics, 2021

Rayleigh–Taylor (RT) and buoyancy-driven (BD) instabilities are driven by gravity in a fluid syst... more Rayleigh–Taylor (RT) and buoyancy-driven (BD) instabilities are driven by gravity in a fluid system with inhomogeneous density. The paper investigates these instabilities for a strongly coupled dusty plasma medium. This medium has been represented here in the framework of the generalized hydrodynamics (GHD) fluid model which treats it as a viscoelastic medium. The incompressible limit of the GHD model is considered here. The RT instability is explored both for gradual and sharp density gradients stratified against gravity. The BD instability is discussed by studying the evolution of a rising bubble (a localized low-density region) and a falling droplet (a localized high-density region) in the presence of gravity. Since both the rising bubble and falling droplet have symmetry in spatial distribution, we observe that a falling droplet process is equivalent to a rising bubble. We also find that both the gravity-driven instabilities get suppressed with increasing coupling strength of th...

Research paper thumbnail of Rotating vortices in two-dimensional inhomogeneous strongly coupled dusty plasmas: Shear and spiral density waves

Dusty plasma experiments can be performed quite easily in a strong coupling regime. In our previo... more Dusty plasma experiments can be performed quite easily in a strong coupling regime. In our previous work [V. S. Dharodi, S. K. Tiwari, and A. Das, Physics of Plasmas 21, 073705 (2014)]PHPAEN1070-664X10.1063/1.4888882, we numerically explored such plasmas with constant density and observed the transverse shear (TS) waves from the rotating vortex. Laboratory dusty plasmas are good examples of homogeneous plasmas; however, heterogeneity (e.g., density, temperature, and charge) may be due to the existence of voids, different domains with different orientations, presence of external forces like magnetic and/or electric, size or charge imbalance, etc. Here, we examine how the density heterogeneity in dusty plasmas responds to the circularly rotating vortex monopoles, specifically, smooth and sharp cutoff. For this purpose, we have carried out a series of two-dimensional fluid simulations in the framework of the incompressible generalized hydrodynamics fluid model. The rotating vortices ar...

Research paper thumbnail of Transport and mixing in strongly coupled dusty plasma medium

Research paper thumbnail of Elastic Turbulence in strongly coupled dusty plasma medium

Research paper thumbnail of Functional form of distribution and correlation functions of strongly coupled plasmas

Research paper thumbnail of Sculpted ultracold neutral plasmas

Physical Review E

Ultracold neutral plasma (UNP) experiments allow for careful control of plasma properties across ... more Ultracold neutral plasma (UNP) experiments allow for careful control of plasma properties across Coulomb coupling regimes. Here, we examine how UNPs can be used to study heterogeneous, nonequilibrium phenomena, including nonlinear waves, transport, hydrodynamics, kinetics, stopping power, and instabilities. Through a series of molecular dynamics simulations, we have explored UNPs formed with spatially modulated ionizing radiation. We have developed a computational model for such sculpted UNPs that includes an ionic screened Coulomb interaction with a spatiotemporal screening length, and Langevin-based spatial ion-electron and ion-neutral collisions. We have also developed a hydrodynamics model and have extracted its field quantities (density, flow velocity, and temperature) from the molecular dynamics simulation data, allowing us to investigate kinetics by examining moment ratios and phase-space dynamics; we find that it is possible to create UNPs that vary from nearly perfect fluids (Euler limit) to highly kinetic plasmas. We have examined plasmas in three geometries: a solid rod, a hollow rod, and a gapped slab; we have studied basic properties of these plasmas, including the spatial Coulomb coupling parameter. By varying the initial conditions, we find that we can design experimental plasmas that would allow the exploration of a wide range of phenomena, including shock and blast waves, stopping power, two-stream instabilities, and much more. Using an evaporative cooling geometry, our results suggest that much larger Coulomb couplings can be achieved, possibly in excess of 10.

Research paper thumbnail of Sub- and super-luminar propagation of structures satisfying Poynting-like theorem for incompressible generalized hydrodynamic fluid model depicting strongly coupled dusty plasma medium

Physics of Plasmas, 2016

The strongly coupled dusty plasma has often been modelled by the Generalized Hydrodynamic (GHD) m... more The strongly coupled dusty plasma has often been modelled by the Generalized Hydrodynamic (GHD) model used for representing visco-elastic fluid systems. The incompressible limit of the model which supports transverse shear wave mode is studied in detail. In particular dipole structures are observed to emit transverse shear waves in both the limits of sub and super luminar propagation, where the structures move slower and faster than the phase velocity of the shear waves, respectively. In the sub-luminar limit the dipole gets engulfed within the shear waves emitted by itself, which then backreacts on it and ultimately the identity of the structure is lost. However, in the super-luminar limit the emission appears like a wake from the tail region of the dipole. The dipole, however, keeps propagating forward with little damping but minimal distortion in its form. A Poynting like conservation law with radiative, convective and dissipative terms being responsible for the evolution of W , which is similar to 'enstrophy' like quantity in normal hydrodynamic fluid systems, has also been constructed for the incompressible GHD equations. The conservation law is shown to be satisfied in all the cases of evolution and collision amidst the nonlinear structures to a great accuracy. It is shown that monopole structures which do not move at all but merely radiate shear waves, the radiative term and dissipative losses solely contribute to the evolution of W. The, dipolar structures, on the other hand, propagate in the medium and hence convection also plays an important role in the evolution of W .

Research paper thumbnail of Turbulence in two dimensional visco - elastic medium

The properties of decaying turbulence is studied with the help of a Generalized Hydrodynamic (GHD... more The properties of decaying turbulence is studied with the help of a Generalized Hydrodynamic (GHD) fluid model in the context of two dimensional visco - elastic medium such as a strongly coupled dusty plasma system. For the incompressible case considered here however, the observations are valid for a wider class of visco - elastic systems not necessarily associated with plasmas only. Our observations show that an initial spectrum that is confined in a limited domain of wave numbers becomes broad, even when the Reynold's number is much less than the critical value required for the onset of turbulence in Newtonian fluids. This is a signature of elastic turbulence where Weissenberg's number also plays a role in the onset of turbulence. This has been reported in several experiments. It is also shown that the existence of memory relaxation time parameter and the transverse shear wave inhibit the normal process (for 2-D systems) of inverse spectral cascade in this case. A detailed simulation study has been carried out for the understanding of this inhibition.

Research paper thumbnail of Elastic Turbulence in strongly coupled dusty plasma medium

Research paper thumbnail of Sanat 10

Research paper thumbnail of Collective dynamics in strongly coupled dusty plasma medium

Journal of Plasma Physics, 2014

A simplified description of dynamical response of strongly coupled medium is desirable in many co... more A simplified description of dynamical response of strongly coupled medium is desirable in many contexts of physics. The dusty plasma medium can play an important role in this regard due to its uniqueness, as its dynamical response typically falls within the perceptible grasp of human senses. Furthermore, even at room temperature and normal densities it can be easily prepared to be in a strongly coupled regime. A simplified phenomenological fluid model based on the visco - elastic behaviour of the medium is often invoked to represent the collective dynamical response of a strongly coupled dusty plasma medium. The manuscript reviews the role of this particular Generalized Hydrodynamic (GHD) fluid model in capturing the collective properties exhibited by the medium. In addition the paper also provides new insights on the collective behaviour predicted by the model for the medium, in terms of coherent structures, instabilities, transport and mixing properties.

Research paper thumbnail of Kelvin-Helmholtz instability in dusty plasma medium: Fluid and particle approach

Journal of Plasma Physics, 2014

The Kelvin-Helmholtz (KH) instability is studied in a two dimensional strongly coupled dusty plas... more The Kelvin-Helmholtz (KH) instability is studied in a two dimensional strongly coupled dusty plasma medium using a fluid approach as well as through a molecular dynamic (MD) simulation. For the fluid description the generalized hydrodynamic (GH) model which treats the strongly coupled dusty plasma as a visco-elastic fluid is adopted. For the MD studies the ensemble of particles are assumed to interact through a Yukawa potential. Both the approaches predict a stabilization of the KH growth rate with an increase in the strong coupling parameter. The present study also delineates the temporal evolution and the interaction of transverse shear waves with the collective dynamics of the dusty plasma medium within the framework of both these approaches.

Research paper thumbnail of Turbulence in strongly coupled dusty plasmas using generalized hydrodynamic description

Physics of Plasmas, 2015

The properties of decaying turbulence is studied with the help of a Generalized Hydrodynamic (GHD... more The properties of decaying turbulence is studied with the help of a Generalized Hydrodynamic (GHD) fluid model in the context of two dimensional visco-elastic medium such as a strongly coupled dusty plasma system. For the incompressible case considered here however, the observations are valid for a wider class of visco-elastic systems not necessarily associated with plasmas only. Our observations show that an initial spectrum that is confined in a limited domain of wave numbers becomes broad, even when the Reynold's number is much less than the critical value required for the onset of turbulence in Newtonian fluids. This is a signature of elastic turbulence where Weissenberg's number also plays a role in the onset of turbulence. This has been reported in several experiments. It is also shown that the existence of memory relaxation time parameter and the transverse shear wave inhibit the normal process (for 2-D systems) of inverse spectral cascade in this case. A detailed simulation study has been carried out for the understanding of this inhibition.

Research paper thumbnail of Visco-elastic fluid simulations of coherent structures in strongly coupled dusty plasma medium

Research paper thumbnail of Evolution of sheared flow structure in visco-elastic fluids

A visco-elastic fluid is represented frequently by the Generalized Hydrodynamic (GHD) model. Unli... more A visco-elastic fluid is represented frequently by the Generalized Hydrodynamic (GHD) model. Unlike Navier Stokes fluid such a system supports transverse shear waves. Thus, while a sheared flow structure in a Navier Stokes fluid system is marred by the development of the Kelvin-Helmholtz (KH) instability alone, in the case of visco-elastic GHD fluid, a shear flow pattern can be expected to be effected by both, namely the KH instability and the emission of transverse shear waves. The manuscript explores this feature in detail. For this purpose (a) the modifications on the KH instability mode due to strong correlations in linear as well as nonlinear regime and (b) the conditions under which either the KH mode and/or the transverse shear wave emission dominates have been identified in the manuscript.

Research paper thumbnail of Signatures of an energetic charge bunch moving in a plasma

A charge bunch moving in a plasma can excite a variety of linear and nonlinear waves in the form ... more A charge bunch moving in a plasma can excite a variety of linear and nonlinear waves in the form of trailing wakes, fore-wake shocks and precursor solitons. These structures can further interact with the background plasma to create secondary effects that can serve as signatures of the passage of the charge bunch. Using particle-in-cell simulations we investigate in detail the dynamics of a plasma system that is being traversed by an energetic charged ion bunch. Using two different shapes of the charge source, namely, an idealized one dimensional line source and a two dimensional thin rectangular source we examine the differences in the nature of the excited wave structures and their consequent impact on the background plasma. Our simulations reveal interesting features such as the dependence of the precursor speeds on the total charge of the ion bunch, local particle trapping, and energization of the trapped ions in various regions along the traversal path leading to the formation of energetic ion beam-lets. The collective excitations and the signatures in the ambient plasma could prove useful in practical applications such as in ion beam heating of plasmas. They can also help in analysing the trajectories of charged objects like space debris orbiting in the ionosphere.

Research paper thumbnail of Kelvin-Helmholtz instability in a strongly coupled dusty plasma medium

Research paper thumbnail of Turbulence in two dimensional visco - elastic medium

arXiv: Plasma Physics, 2014

The properties of decaying turbulence is studied with the help of a Generalized Hydrodynamic (GHD... more The properties of decaying turbulence is studied with the help of a Generalized Hydrodynamic (GHD) fluid model in the context of two dimensional visco - elastic medium such as a strongly coupled dusty plasma system. For the incompressible case considered here however, the observations are valid for a wider class of visco - elastic systems not necessarily associated with plasmas only. Our observations show that an initial spectrum that is confined in a limited domain of wave numbers becomes broad, even when the Reynold’s number is much less than the critical value required for the onset of turbulence in Newtonian fluids. This is a signature of elastic turbulence where Weissenberg’s number also plays a role in the onset of turbulence. This has been reported in several experiments. It is also shown that the existence of memory relaxation time parameter and the transverse shear wave inhibit the normal process (for 2-D systems) of inverse spectral cascade in this case. A detailed simulat...

Research paper thumbnail of A numerical study of gravity-driven instability in strongly coupled dusty plasma. Part 1. Rayleigh–Taylor instability and buoyancy-driven instability

Journal of Plasma Physics, 2021

Rayleigh–Taylor (RT) and buoyancy-driven (BD) instabilities are driven by gravity in a fluid syst... more Rayleigh–Taylor (RT) and buoyancy-driven (BD) instabilities are driven by gravity in a fluid system with inhomogeneous density. The paper investigates these instabilities for a strongly coupled dusty plasma medium. This medium has been represented here in the framework of the generalized hydrodynamics (GHD) fluid model which treats it as a viscoelastic medium. The incompressible limit of the GHD model is considered here. The RT instability is explored both for gradual and sharp density gradients stratified against gravity. The BD instability is discussed by studying the evolution of a rising bubble (a localized low-density region) and a falling droplet (a localized high-density region) in the presence of gravity. Since both the rising bubble and falling droplet have symmetry in spatial distribution, we observe that a falling droplet process is equivalent to a rising bubble. We also find that both the gravity-driven instabilities get suppressed with increasing coupling strength of th...

Research paper thumbnail of Rotating vortices in two-dimensional inhomogeneous strongly coupled dusty plasmas: Shear and spiral density waves

Dusty plasma experiments can be performed quite easily in a strong coupling regime. In our previo... more Dusty plasma experiments can be performed quite easily in a strong coupling regime. In our previous work [V. S. Dharodi, S. K. Tiwari, and A. Das, Physics of Plasmas 21, 073705 (2014)]PHPAEN1070-664X10.1063/1.4888882, we numerically explored such plasmas with constant density and observed the transverse shear (TS) waves from the rotating vortex. Laboratory dusty plasmas are good examples of homogeneous plasmas; however, heterogeneity (e.g., density, temperature, and charge) may be due to the existence of voids, different domains with different orientations, presence of external forces like magnetic and/or electric, size or charge imbalance, etc. Here, we examine how the density heterogeneity in dusty plasmas responds to the circularly rotating vortex monopoles, specifically, smooth and sharp cutoff. For this purpose, we have carried out a series of two-dimensional fluid simulations in the framework of the incompressible generalized hydrodynamics fluid model. The rotating vortices ar...

Research paper thumbnail of Transport and mixing in strongly coupled dusty plasma medium

Research paper thumbnail of Elastic Turbulence in strongly coupled dusty plasma medium

Research paper thumbnail of Functional form of distribution and correlation functions of strongly coupled plasmas

Research paper thumbnail of Sculpted ultracold neutral plasmas

Physical Review E

Ultracold neutral plasma (UNP) experiments allow for careful control of plasma properties across ... more Ultracold neutral plasma (UNP) experiments allow for careful control of plasma properties across Coulomb coupling regimes. Here, we examine how UNPs can be used to study heterogeneous, nonequilibrium phenomena, including nonlinear waves, transport, hydrodynamics, kinetics, stopping power, and instabilities. Through a series of molecular dynamics simulations, we have explored UNPs formed with spatially modulated ionizing radiation. We have developed a computational model for such sculpted UNPs that includes an ionic screened Coulomb interaction with a spatiotemporal screening length, and Langevin-based spatial ion-electron and ion-neutral collisions. We have also developed a hydrodynamics model and have extracted its field quantities (density, flow velocity, and temperature) from the molecular dynamics simulation data, allowing us to investigate kinetics by examining moment ratios and phase-space dynamics; we find that it is possible to create UNPs that vary from nearly perfect fluids (Euler limit) to highly kinetic plasmas. We have examined plasmas in three geometries: a solid rod, a hollow rod, and a gapped slab; we have studied basic properties of these plasmas, including the spatial Coulomb coupling parameter. By varying the initial conditions, we find that we can design experimental plasmas that would allow the exploration of a wide range of phenomena, including shock and blast waves, stopping power, two-stream instabilities, and much more. Using an evaporative cooling geometry, our results suggest that much larger Coulomb couplings can be achieved, possibly in excess of 10.

Research paper thumbnail of Sub- and super-luminar propagation of structures satisfying Poynting-like theorem for incompressible generalized hydrodynamic fluid model depicting strongly coupled dusty plasma medium

Physics of Plasmas, 2016

The strongly coupled dusty plasma has often been modelled by the Generalized Hydrodynamic (GHD) m... more The strongly coupled dusty plasma has often been modelled by the Generalized Hydrodynamic (GHD) model used for representing visco-elastic fluid systems. The incompressible limit of the model which supports transverse shear wave mode is studied in detail. In particular dipole structures are observed to emit transverse shear waves in both the limits of sub and super luminar propagation, where the structures move slower and faster than the phase velocity of the shear waves, respectively. In the sub-luminar limit the dipole gets engulfed within the shear waves emitted by itself, which then backreacts on it and ultimately the identity of the structure is lost. However, in the super-luminar limit the emission appears like a wake from the tail region of the dipole. The dipole, however, keeps propagating forward with little damping but minimal distortion in its form. A Poynting like conservation law with radiative, convective and dissipative terms being responsible for the evolution of W , which is similar to 'enstrophy' like quantity in normal hydrodynamic fluid systems, has also been constructed for the incompressible GHD equations. The conservation law is shown to be satisfied in all the cases of evolution and collision amidst the nonlinear structures to a great accuracy. It is shown that monopole structures which do not move at all but merely radiate shear waves, the radiative term and dissipative losses solely contribute to the evolution of W. The, dipolar structures, on the other hand, propagate in the medium and hence convection also plays an important role in the evolution of W .

Research paper thumbnail of Turbulence in two dimensional visco - elastic medium

The properties of decaying turbulence is studied with the help of a Generalized Hydrodynamic (GHD... more The properties of decaying turbulence is studied with the help of a Generalized Hydrodynamic (GHD) fluid model in the context of two dimensional visco - elastic medium such as a strongly coupled dusty plasma system. For the incompressible case considered here however, the observations are valid for a wider class of visco - elastic systems not necessarily associated with plasmas only. Our observations show that an initial spectrum that is confined in a limited domain of wave numbers becomes broad, even when the Reynold's number is much less than the critical value required for the onset of turbulence in Newtonian fluids. This is a signature of elastic turbulence where Weissenberg's number also plays a role in the onset of turbulence. This has been reported in several experiments. It is also shown that the existence of memory relaxation time parameter and the transverse shear wave inhibit the normal process (for 2-D systems) of inverse spectral cascade in this case. A detailed simulation study has been carried out for the understanding of this inhibition.

Research paper thumbnail of Elastic Turbulence in strongly coupled dusty plasma medium

Research paper thumbnail of Sanat 10

Research paper thumbnail of Collective dynamics in strongly coupled dusty plasma medium

Journal of Plasma Physics, 2014

A simplified description of dynamical response of strongly coupled medium is desirable in many co... more A simplified description of dynamical response of strongly coupled medium is desirable in many contexts of physics. The dusty plasma medium can play an important role in this regard due to its uniqueness, as its dynamical response typically falls within the perceptible grasp of human senses. Furthermore, even at room temperature and normal densities it can be easily prepared to be in a strongly coupled regime. A simplified phenomenological fluid model based on the visco - elastic behaviour of the medium is often invoked to represent the collective dynamical response of a strongly coupled dusty plasma medium. The manuscript reviews the role of this particular Generalized Hydrodynamic (GHD) fluid model in capturing the collective properties exhibited by the medium. In addition the paper also provides new insights on the collective behaviour predicted by the model for the medium, in terms of coherent structures, instabilities, transport and mixing properties.

Research paper thumbnail of Kelvin-Helmholtz instability in dusty plasma medium: Fluid and particle approach

Journal of Plasma Physics, 2014

The Kelvin-Helmholtz (KH) instability is studied in a two dimensional strongly coupled dusty plas... more The Kelvin-Helmholtz (KH) instability is studied in a two dimensional strongly coupled dusty plasma medium using a fluid approach as well as through a molecular dynamic (MD) simulation. For the fluid description the generalized hydrodynamic (GH) model which treats the strongly coupled dusty plasma as a visco-elastic fluid is adopted. For the MD studies the ensemble of particles are assumed to interact through a Yukawa potential. Both the approaches predict a stabilization of the KH growth rate with an increase in the strong coupling parameter. The present study also delineates the temporal evolution and the interaction of transverse shear waves with the collective dynamics of the dusty plasma medium within the framework of both these approaches.

Research paper thumbnail of Turbulence in strongly coupled dusty plasmas using generalized hydrodynamic description

Physics of Plasmas, 2015

The properties of decaying turbulence is studied with the help of a Generalized Hydrodynamic (GHD... more The properties of decaying turbulence is studied with the help of a Generalized Hydrodynamic (GHD) fluid model in the context of two dimensional visco-elastic medium such as a strongly coupled dusty plasma system. For the incompressible case considered here however, the observations are valid for a wider class of visco-elastic systems not necessarily associated with plasmas only. Our observations show that an initial spectrum that is confined in a limited domain of wave numbers becomes broad, even when the Reynold's number is much less than the critical value required for the onset of turbulence in Newtonian fluids. This is a signature of elastic turbulence where Weissenberg's number also plays a role in the onset of turbulence. This has been reported in several experiments. It is also shown that the existence of memory relaxation time parameter and the transverse shear wave inhibit the normal process (for 2-D systems) of inverse spectral cascade in this case. A detailed simulation study has been carried out for the understanding of this inhibition.

Research paper thumbnail of Visco-elastic fluid simulations of coherent structures in strongly coupled dusty plasma medium

Research paper thumbnail of Evolution of sheared flow structure in visco-elastic fluids

A visco-elastic fluid is represented frequently by the Generalized Hydrodynamic (GHD) model. Unli... more A visco-elastic fluid is represented frequently by the Generalized Hydrodynamic (GHD) model. Unlike Navier Stokes fluid such a system supports transverse shear waves. Thus, while a sheared flow structure in a Navier Stokes fluid system is marred by the development of the Kelvin-Helmholtz (KH) instability alone, in the case of visco-elastic GHD fluid, a shear flow pattern can be expected to be effected by both, namely the KH instability and the emission of transverse shear waves. The manuscript explores this feature in detail. For this purpose (a) the modifications on the KH instability mode due to strong correlations in linear as well as nonlinear regime and (b) the conditions under which either the KH mode and/or the transverse shear wave emission dominates have been identified in the manuscript.