Search

A new combined PIC-MCC approach is developed for accurate and fast simulation of a radio frequency discharge at low gas pressure and high density of plasma. Test calculations of transition between different modes of electron heating in a ccrf discharge in helium and argon show a good agreement with experimental data. We demonstrate high efficiency of the combined PIC-MCC algorithm, especially for the collisionless regime of electron heating., Comment: 6 paged, 8 figures

The transitions between the different vortex states of thin mesoscopic superconducting disks and rings are studied using the non-linear Ginzburg-Landau functional. They are saddle points of the free energy representing the energy barrier which has to be overcome for transition between the different vortex states. In small superconducting disks and rings the saddle point state between two giant vortex states, and in larger systems the saddle point state between a multivortex state and a giant vortex state and between two multivortex states is obtained. The shape and the height of the nucleation barrier is investigated for different disk and ring configurations., Comment: 10 pages, 18 figures

As first pointed out by Bardeen and Ginzburg in the early sixties, the amount of magnetic flux carried by vortices depends on their distance to the sample edge and can be smaller than one flux quantum, f0 = h/2e. In bulk superconductors, this reduction of flux becomes negligible already at submicron distances from the edge but, in thin films, the effect may survive at much larger distances. In the absence of any experimental observation, such flux reduction is perceived to be an exotic or unimportant effect, and it is often assumed that magnetic field enters type-II superconductors in units of f0. Here we report the measurements of the amount of flux associated with the entrance of individual vortices in a superconducting film and show that the flux that they bring in, always differs substantially from f0. We have observed vortices that carry as little as 0.001 f0 as well as seemingly "negative vortices" whose penetration leads to the expulsion of magnetic field. We distinguish two phenomena responsible for non-quantized flux penetration: the finite-size effect and a non-linear screening of magnetic field in the presence of a surface barrier. The latter effect has previously not been considered but is likely to cause non-quantized penetration in many cases., Comment: pdf-file of a slighly longer version than appeared in Nature

The groundstate energy and binding energy of an exciton, confined in a^M quantum disk, are calculated as a function of an external magnetic field. The confinement potential is a hard wall of finite height. The diamagnetic shift is investigated for magnetic fields up to 40$T$. Our results are applied to $In_{y}Al_{1-y}As/Al_{x}Ga_{1-x}As$ self-assembled quantum dots and very good agreement with experiments is obtained. Furthermore, we investigated the influence of the dot size on the diamagnetic shift by changing the disk radius. The exciton excited states are found as a function of the magnetic field. The relative angular momentum is not a quantum number and changes with the magnetic field strength., Comment: 10 pages, 17 figures

Using the non-linear Ginzburg-Landau (GL) theory, we study the magnetic response of different shaped samples in the field-cooled regime (FC). For high external magnetic fluxes, the conventional diamagnetic response under cooling down can be followed by the paramagnetic Meissner effect (PME). A second-order transition from a giant vortex state to a multi-vortex state, with the same vorticity, occurs at the second critical field which leads to the suppression of PME., Comment: 4 pages, 6 figures

Melting of two dimensional (2D) clusters of classical particles is studied using Brownian dynamics and Langevin molecular dynamics simulations. The particles are confined by a circular hard wall or a parabolic external potential and interact through a dipole or a screened Coulomb potential. We found that with decreasing strength of the inter--particle interaction clusters with short-range inter-particle interaction which are confined by a hard wall exhibit a re-entrant behavior in its orientational order., Comment: 4 pages, 8 figures

Phase transitions between different (i.e. giant and multi-vortex) superconducting states and between the superconducting-normal state of mesoscopic disks and rings are studied in the presence of an external magnetic field by solving the two non-linear Ginzburg-Landau equations self-consistently. The flux through a circular disk with a hole in the middle is not quantized., Comment: 8 pages, 10 figures; to appear in Physica C (proceedings of the conference on Vortex matter, Crete (september 1999)

Using a linear analysis, we study the stability of giant-vortex states in very thin disks. The vortex expulsion and penetration fields are obtained for finite thickness disks from a numerical solution of the non-linear Ginzburg-Landau (GL) equations. Using an extension of the London approximation, in which the phase distribution of the order parameter is prescribed and the superconducting electron density is found numerically, we consider the free energy behavior for transitions between different superconducting states., Comment: 8 pages, 5 figures Accepted for publication in Physica C

The superconducting state of a thin superconducting disk with a hole is studied within the non-linear Ginzburg-Landau theory in which the demagnetization effect is accurately taken into account. We find that the flux through the hole is not quantized, the superconducting state is stabilized with increasing size of the hole for fixed radius of the disk, and a transition to a multi-vortex state is found if the disk is sufficiently large. Breaking the circular summetry through a non central location of the hole in the disk enhances the multi-vortex state., Comment: 11 pages, 23 figures (postscript). To appear in Physical Review B, Vol. 61 (2000)

Using an expansion of the order parameter over the eigenfunctions of the linearized first Ginzburg-Landau (GL) equation, we obtain numerically the saddle points of the free energy separating the stable states with different number of vortices. In contrast to known surface and geometrical barrier models, we find that in a wide range of magnetic fields below the penetration field, the saddle point state for flux penetration into a disk does not correspond to a vortex located nearby the sample boundary, but to a region of suppressed superconductivity at the disk edge with no winding of the current, and which is {\it a nucleus} for the following vortex creation. The height of this {\it nucleation barrier}, which determines the time of flux penetration, is calculated for different disk radii and magnetic fields., Comment: Accepted for publication in Physical Review Letters

The stability and melting transition of a single layer and a bilayer crystal consisting of charged particles interacting through a Coulomb or a screened Coulomb potential is studied using the Monte-Carlo technique. A new melting criterion is formulated which we show to be universal for bilayer as well as for single layer crystals in the case of (screened) Coulomb, Lennard--Jones and 1/r^{12} repulsive inter-particle interactions. The melting temperature for the five different lattice structures of the bilayer Wigner crystal is obtained, and a phase diagram is constructed as a function of the interlayer distance. We found the surprising result that the square lattice has a substantial larger melting temperature as compared to the other lattice structures. This is a consequence of the specific topology of the defects which are created with increasing temperature and which have a larger energy as compared to the defects in e.g. a hexagonal lattice., Comment: Accepted for publication in Physical Review B

The nucleation field for surface superconductivity, $H_{c3}$, depends on the geometrical shape of the mesoscopic superconducting sample and is substantially enhanced with decreasing sample size. As an example we studied circular, square, triangular and wedge shaped disks. For the wedge the nucleation field diverges as $H_{c3}/H_{c2}=\sqrt{3}/\alpha$ with decreasing angle ($\alpha$) of the wedge, where $H_{c2}$ is the bulk upper critical field., Comment: 4 pages, 3 figures. Accepted for publication in Phys. Rev. B

The melting transition of the five different lattices of a bilayer crystal is studied using the Monte-Carlo technique. We found the surprising result that the square lattice has a substantial larger melting temperature as compared to the other lattice structures, which is a consequence of the specific topology of the temperature induced defects. A new melting criterion is formulated which we show to be universal for bilayers as well as for single layer crystals., Comment: 4 pages, 5 figures (postscript files). Accepted in Physical Review Letters

Using the non-linear Ginzburg-Landau (GL) eqs. type I superconducting disks of finite radius ($R$) and thickness ($d$) are studied in a perpendicular magnetic field. Depending on $R$ and $d$, first or second order phase transitions are found for the normal to superconducting state. For sufficiently large $R$ several transitions in the superconducting phase are found corresponding to different angular momentum giant vortex states. In increasing magnetic field the superconductor is in its ground state, while in field down sweep it is possible to drive the system into metastable states. We also present a quantitative analysis of the relation between the detector output and the sample magnetization. The latter, and the incorporation of the finite thickness of the disks, are essential in order to obtain quantitative agreement with experiment., Comment: A brief review with new results

Solving numerically the 3D non linear Ginzburg-Landau (GL) equations, we study equilibrium and nonequilibrium phase transitions between different superconducting states of mesoscopic disks which are thinner than the coherence length and the penetration depth. We have found a smooth transition from a multi-vortex superconducting state to a giant vortex state with increasing both the disk thickness and the magnetic field. A vortex phase diagram is obtained which shows, as function of the magnetic field, a re-entrant behavior between the multi-vortex and the giant vortex state., Comment: 5 figures (post script files) included

We investigate the effect of the position of the donor in the quantum well on the energy spectrum and the oscillator strength of the D- system in the presence of a perpendicular magnetic field. As a function of the magnetic field we find that when the D- centers are placed sufficiently off-center they undergo singlet-triplet transitions which are similar to those found in many-electron parabolic quantum dots. The main difference is that the number of such transitions depends on the position of the donor and only a finite number of such singlet-triplet transitions are found as function of the strength of the magnetic field. For sufficiently large magnetic fields the two electron system becomes unbound. For the near center D- system no singlet-triplet and no unbinding of the D- is found with increasing magnetic field. A magnetic field vs. donor position phase diagram is presented that depends on the width of the quantum well., Comment: 16 pages, 17 figures. Accepted for publication in Phys. Rev. B

The magnetization behavior of mesoscopic superconducting disks can show hysteretic behavior which we explain by using the Ginzburg-Landau (GL) theory and properly taking into account the de-magnetization effects due to geometrical form factors. In large disks the Bean-Livingston surface barrier is responsible for the hysteresis. While in small disks a volume barrier is responsible for this hysteresis. It is shown that although the sample magnetization is diamagnetic (negative), the measured magnetization can be positive at certain fields as observed experimentally, which is a consequence of the de-magnetization effects and the experimental set up., Comment: Latex file, 4 ps files

The groundstate configuration and the eigenmodes of two parallel two-dimensional classical atoms are obtained as function of the inter-atomic distance (d). The classical particles are confined by identical harmonic wells and repel each other through a Coulomb potential. As function of d we find several structural transitions which are of first or second order. For first (second) order transitions the first (second) derivative of the energy with respect to d is discontinuous, the radial position of the particles changes discontinuously (continuously) and the frequency of the eigenmodes exhibit a jump (one mode becomes soft, i.e. its frequency becomes zero)., Comment: 4 pages, RevTex, 5 ps figures, to appear in Phys.Rev.Lett

Solutions of Ginzburg-Landau eqns. coupled with three dimensional Maxwell eqns. reveal intriguing magnetic response of small superconducting particles, qualitatively different from the two dimensional approximation but in agreement with recent experiments. Depending on the radius and thickness first or second order transitions are found for the normal to superconducting state. For a sufficient large radius of the disc several transitions in the superconducting phase are obtained which correspond to different angular momentum giant vortex states. The incorporation of the finite thickness in the calculation is crucial in order to obtain agreement with the position and the size of these jumps, and the line shape and magnitude of the magnetization curves., Comment: Some minor corrections and more accurate citation to earlier works. To appear in Phys. Rev. Lett

We present a study of the spectral properties like the energy spectrum, the eigenmodes and density of states of a classical finite system of two-dimensional (2D) charged particles which are confined by a quadratic potential. Using the method of Newton optimization we obtain the ground state and the metastable states. For a given configuration the eigenvectors and eigenfrequencies for the normal modes are obtained using the Householder diagonalization technique for the dynamical matrix whose elements are the second derivative of the potential energy. For small clusters the lowest excitation corresponds to an intershell rotation. The energy barrier for such rotations is calculated. For large clusters the lowest excitation consists of a vortex/anti-vortex pair. The Lindeman melting criterion is used to calculate the order-disorder transition temperature for intershell rotation and intershell diffusion. The value of the transition temperature at which intershell rotation becomes possible depends very much on the configuration of the cluster, i.e. the distribution of the particles between the different shells. Magic numbers are associated to clusters which are most stable against intershell rotation. The specific heat of the cluster is also calculated using, Comment: 27 pages, 2 tables, 13 figures, RevTex-file