Unveiling the Electronic Structure of Pseudotetragonal WO3 Thin Films

WO3 is a 5d compound that undergoes several structural transitions in its bulk form. Its versatility is well-documented, with a wide range of applications, such as flexopiezoelectricity, electrochromism, gating-induced phase transitions, and its ability to improve the performance of Li-based batteries. The synthesis of WO3 thin films holds promise in stabilizing electronic phases for practical applications. However, despite its potential, the electronic structure of this material remains experimentally unexplored. Furthermore, its thermal instability limits its use in certain technological devices. Here, we employ tensile strain to stabilize WO3 thin films, which we call the pseudotetragonal phase, and investigate its electronic structure using a combination of photoelectron spectroscopy and density functional theory calculations. This study reveals the Fermiology of the system, notably identifying significant energy splittings between different orbital manifolds arising from atomic distortions. These splittings, along with the system's thermal stability, offer a potential avenue for controlling inter- and intraband scattering for electronic applications.

Free surfaces recast superconductivity in few-monolayer MgB2: Combined first-principles and ARPES demonstration

Two-dimensional materials are known to harbour properties very different from those of their bulk counterparts. Recent years have seen the rise of atomically thin superconductors, with a caveat that superconductivity is strongly depleted unless enhanced by specific substrates, intercalants or adatoms. Surprisingly, the role in superconductivity of electronic states originating from simple free surfaces of two-dimensional materials has remained elusive to date. Here, based on first-principles calculations, anisotropic Eliashberg theory, and angle-resolved photoemission spectroscopy (ARPES), we show that surface states in few-monolayer MgB2 make a major contribution to the superconducting gap spectrum and density of states, clearly distinct from the widely known, bulk-like σ- and π-gaps. As a proof of principle, we predict and measure the gap opening on the magnesium-based surface band up to a critical temperature as high as ~30 K for merely six monolayers thick MgB2. These findings establish free surfaces as an unavoidable ingredient in understanding and further tailoring of superconductivity in atomically thin materials.

A novel explanation for the increased conductivity in annealed Al-doped ZnO: an insight into migration of aluminum and displacement of zinc

A combined experimental and first-principles study is performed to study the origin of conductivity in ZnO:Al nanoparticles synthesized under controlled conditions via a reflux route using benzylamine as a solvent. The experimental characterization of the samples by Raman, nuclear magnetic resonance (NMR) and conductivity measurements indicates that upon annealing in nitrogen, the Al atoms at interstitial positions migrate to the substitutional positions, creating at the same time Zn interstitials. We provide evidence for the fact that the formed complex of Al_Zn and Zn_i corresponds to the origin of the Knight shifted peak (KS) we observe in ²⁷Al NMR. As far as we know, the role of this complex has not been discussed in the literature to date. However, our first-principles calculations show that such a complex is indeed energetically favoured over the isolated Al interstitial positions. In our calculations we also address the charge state of the Al interstitials. Further, Zn interstitials can migrate from Al_Zn and possibly also form Zn clusters, leading to the observed increased conductivity.

Structural and electronic properties of defects at grain boundaries in CuInSe2

We report on a first-principles study of the structural and electronic properties of a Σ3 (112) grain boundary model in CuInSe₂. The study focuses on a coherent, stoichiometry preserving, cation-Se terminated grain boundary, addressing the properties of the grain boundary as such, as well as the effect of well known defects in CuInSe₂. We show that in spite of its apparent simplicity, such a grain boundary exhibits a very rich phenomenology, providing an explanation for several of the experimentally observed properties of grain boundaries in CuInSe₂ thin films. In particular, we show that the combined effect of Cu vacancies and cation antisites can result in the observed Cu depletion with no In enrichment at the grain boundaries. Furthermore, Cu vacancies are unlikely to produce a hole barrier at the grain boundaries, but Na may indeed have such an effect. We find that Na-on-Cu defects will tend to form abundantly at the grain boundaries, and can provide a mechanism for the carrier depletion and/or type inversion experimentally reported.

Ab initio study of hydrogenic effective mass impurities in Si nanowires

The effect of B and P dopants on the band structure of Si nanowires is studied using electronic structure calculations based on density functional theory. At low concentrations a dispersionless band is formed, clearly distinguishable from the valence and conduction bands. Although this band is evidently induced by the dopant impurity, it turns out to have purely Si character. These results can be rigorously analyzed in the framework of effective mass theory. In the process we resolve some common misconceptions about the physics of hydrogenic shallow impurities, which can be more clearly elucidated in the case of nanowires than would be possible for bulk Si. We also show the importance of correctly describing the effect of dielectric confinement, which is not included in traditional electronic structure calculations, by comparing the obtained results with those of G₀W₀ calculations.

The work function of few-layer graphene

A theoretical and experimental study of the work function of few-layer graphene is reported. The influence of the number of layers on the work function is investigated in the presence of a substrate, a molecular dipole layer, and combinations of the two. The work function of few-layer graphene is almost independent of the number of layers with only a difference between monolayer and multilayer graphene of about 60 meV. In the presence of a charge-donating substrate the charge distribution is found to decay exponentially away from the substrate and this is directly reflected in the work function of few-layer graphene. A dipole layer changes the work function only when placed in between the substrate and few-layer graphene through a change of the charge transfer between the two.

Extension of the basis set of linearized augmented plane wave (LAPW) method by using supplemented tight binding basis functions

In order to increase the accuracy of the linearized augmented plane wave (LAPW) method, we present a new approach where the plane wave basis function is augmented by two different atomic radial components constructed at two different linearization energies corresponding to two different electron bands (or energy windows). We demonstrate that this case can be reduced to the standard treatment within the LAPW paradigm where the usual basis set is enriched by the basis functions of the tight binding type, which go to zero with zero derivative at the sphere boundary. We show that the task is closely related with the problem of extended core states which is currently solved by applying the LAPW method with local orbitals (LAPW+LO). In comparison with LAPW+LO, the number of supplemented basis functions in our approach is doubled, which opens up a new channel for the extension of the LAPW and LAPW+LO basis sets. The appearance of new supplemented basis functions absent in the LAPW+LO treatment is closely related with the existence of the u̇l-component in the canonical LAPW method. We discuss properties of additional tight binding basis functions and apply the extended basis set for computation of electron energy bands of lanthanum (face and body centered structures) and hexagonal close packed lattice of cadmium. We demonstrate that the new treatment gives lower total energies in comparison with both canonical LAPW and LAPW+LO, with the energy difference more pronounced for intermediate and poor LAPW basis sets.

Giant magnetic anisotropy in doped single layer molybdenum disulfide and fluorographene

Stable monolayer materials based on existing, well known and stable two-dimensional crystal fluorographene and molybdenum disulfide are predicted to exhibit a huge magnetocrystalline anisotropy when functionalized with adsorbed transition metal atoms at vacant sides. Ab initio calculations within the density-functional theory formalism were performed to investigate the adsorption of the transitional metals in a single S (or F) vacancy of monolayer molybdenum disulfide (or fluorographene). We found strong bonding of the transitional metal atoms to the vacant sites with binding energies ranging from 2.5 to 5.2 eV. Our calculations revealed that these systems with adsorbed metal atoms exhibit a magnetic anisotropy, specifically the structures including Os and Ir show a giant magnetocrystalline anisotropy energy of 31-101 meV. Our results demonstrate the possibility of obtaining stable monolayer materials with huge magnetocrystalline anisotropy based on preexisting, well known and stable two-dimensional crystals: fluorographene and molybdenum disulfide. We believe that the results obtained here are useful not only for deeper understanding of the origin of magnetocrystalline anisotropy but also for the design of monolayer optoelectronic devices with novel functionalities.

Comment on "Generalized exclusion processes: Transport coefficients"

In a recent paper, Arita et al. [Phys. Rev. E 90, 052108 (2014)PLEEE81539-375510.1103/PhysRevE.90.052108] consider the transport properties of a class of generalized exclusion processes. Analytical expressions for the transport-diffusion coefficient are derived by ignoring correlations. It is claimed that these expressions become exact in the hydrodynamic limit. In this Comment, we point out that (i) the influence of correlations upon the diffusion does not vanish in the hydrodynamic limit, and (ii) the expressions for the self- and transport diffusion derived by Arita et al. are special cases of results derived in Becker et al. [Phys. Rev. Lett. 111, 110601 (2013)PRLTAO0031-900710.1103/PhysRevLett.111.110601].

Native point defects in CuIn(1-x)Ga(x)Se2: hybrid density functional calculations predict the origin of p- and n-type conductivity

We have performed a first-principles study of the p- and n-type conductivity in CuIn(1-x)Ga(x)Se2 due to native point defects, based on the HSE06 hybrid functional. Band alignment shows that the band gap becomes larger with x due to the increasing conduction band minimum, rendering it hard to establish n-type conductivity in CuGaSe2. From the defect formation energies, we find that In/GaCu is a shallow donor, while V(Cu), V(In/Ga) and Cu(In/Ga) act as shallow acceptors. Using the total charge neutrality of ionized defects and intrinsic charge carriers to determine the Fermi level, we show that under In-rich growth conditions InCu causes strongly n-type conductivity in CuInSe2. Under increasingly In-poor growth conditions, the conductivity type in CuInSe2 alters to p-type and compensation of the acceptors by In(Cu) reduces, as also observed in photoluminescence experiments. In CuGaSe2, the native acceptors pin the Fermi level far away from the conduction band minimum, thus inhibiting n-type conductivity. On the other hand, CuGaSe2 shows strong p-type conductivity under a wide range of Ga-poor growth conditions. Maximal p-type conductivity in CuIn(1-x)Ga(x)Se2 is reached under In/Ga-poor growth conditions, in agreement with charge concentration measurements on samples with In/Ga-poor stoichiometry, and is primarily due to the dominant acceptor Cu(In/Ga).

The role of the VZn-NO-H complex in the p-type conductivity in ZnO

Past research efforts aiming at obtaining stable p-type ZnO have been based on complexes involving nitrogen doping. A recent experiment by (J. G. Reynolds et al., Appl. Phys. Lett., 2013, 102, 152114) demonstrated a significant (∼10(18) cm(-3)) p-type behavior in N-doped ZnO films after appropriate annealing. The p-type conductivity was attributed to a VZn-NO-H shallow acceptor complex, formed by a Zn vacancy (VZn), N substituting O (NO), and H interstitial (Hi). We present here a first-principles hybrid functional study of this complex compared to the one without hydrogen. Our results confirm that the VZn-NO-H complex acts as an acceptor in ZnO. We find that H plays an important role, because it lowers the formation energy of the complex with respect to VZn-NO, a complex known to exhibit (unstable) p-type behavior. However, this additional H atom also occupies the hole level at the origin of the shallow behavior of VZn-NO, leaving only two states empty higher in the band gap and making the VZn-NO-H complex a deep acceptor. Therefore, we conclude that the cause of the observed p-type conductivity in experiment is not the presence of the VZn-NO-H complex, but probably the formation of the VZn-NO complex during the annealing process.

Diffusion of interacting particles in discrete geometries: Equilibrium and dynamical properties

We expand on a recent study of a lattice model of interacting particles [Phys. Rev. Lett. 111, 110601 (2013)PRLTAO0031-900710.1103/PhysRevLett.111.110601]. The adsorption isotherm and equilibrium fluctuations in particle number are discussed as a function of the interaction. Their behavior is similar to that of interacting particles in porous materials. Different expressions for the particle jump rates are derived from transition-state theory. Which expression should be used depends on the strength of the interparticle interactions. Analytical expressions for the self- and transport diffusion are derived when correlations, caused by memory effects in the environment, are neglected. The diffusive behavior is studied numerically with kinetic Monte Carlo (kMC) simulations, which reproduces the diffusion including correlations. The effect of correlations is studied by comparing the analytical expressions with the kMC simulations. It is found that the Maxwell-Stefan diffusion can exceed the self-diffusion. To our knowledge, this is the first time this is observed. The diffusive behavior in one-dimensional and higher-dimensional systems is qualitatively the same, with the effect of correlations decreasing for increasing dimension. The length dependence of both the self- and transport diffusion is studied for one-dimensional systems. For long lengths the self-diffusion shows a 1/L dependence. Finally, we discuss when agreement with experiments and simulations can be expected. The assumption that particles in different cavities do not interact is expected to hold quantitatively at low and medium particle concentrations if the particles are not strongly interacting.

Electronic and optical properties of core-shell nanowires in a magnetic field

The electronic and optical properties of zincblende nanowires are investigated in the presence of a uniform magnetic field directed along the [001] growth direction within the k · p method. We focus our numerical study on core-shell nanowires consisting of the III-V materials GaAs, Al(x)Ga(1-x)As and (Al(y)Ga(1-y))₀.₅₁In₀.₄₉P. Nanowires with electrons confined in the core exhibit a Fock-Darwin-like spectrum, whereas nanowires with electrons confined in the shell show Aharonov-Bohm oscillations. Thus, by properly choosing the core and the shell materials of the nanowire, the optical properties in a magnetic field can be tuned in very different ways.

Perovskite transparent conducting oxides: an ab initio study

We present an ab initio study of the electronic structure and of the formation energies of various point defects in BaSnO3 and SrGeO3. We show that La and Y impurities substituting Ba or Sr are shallow donors with a preferred 1 + charge state. These defects have a low formation energy within all the suitable equilibrium growth conditions considered. Oxygen vacancies behave as shallow donors as well, preferring the 2 + charge state. Their formation energies, however, are higher in most growth conditions, indicating a limited contribution to conductivity. The calculated electron effective mass in BaSnO3, with a value of 0.21 m(e), and the very high mobility reported recently in La-doped BaSnO3 single-crystals, suggest that remarkably low scattering rates can be achieved in the latter. In the case of SrGeO3, our results point to carrier density and mobility values in the low range for typical polycrystalline TCOs, in line with experiment.

Diffusion of interacting particles in discrete geometries

We evaluate the self-diffusion and transport diffusion of interacting particles in a discrete geometry consisting of a linear chain of cavities, with interactions within a cavity described by a free-energy function. Exact analytical expressions are obtained in the absence of correlations, showing that the self-diffusion can exceed the transport diffusion if the free-energy function is concave. The effect of correlations is elucidated by comparison with numerical results. Quantitative agreement is obtained with recent experimental data for diffusion in a nanoporous zeolitic imidazolate framework material, ZIF-8.

Quasiparticle band structure of rocksalt-CdO determined using maximally localized Wannier functions

CdO in the rocksalt structure is an indirect band gap semiconductor. Thus, in order to determine its band gap one needs to calculate the complete band structure. However, in practice, the exact evaluation of the quasiparticle band structure for the large number of k-points which constitute the different symmetry lines in the Brillouin zone can be an extremely demanding task compared to the standard density functional theory (DFT) calculation. In this paper we report the full quasiparticle band structure of CdO using a plane-wave pseudopotential approach. In order to reduce the computational effort and time, we make use of maximally localized Wannier functions (MLWFs). The MLWFs offer a highly accurate method for interpolation of the DFT or GW band structure from a coarse k-point mesh in the irreducible Brillouin zone, resulting in a much reduced computational effort. The present paper discusses the technical details of the scheme along with the results obtained for the quasiparticle band gap and the electron effective mass.

Spin-orbit-mediated manipulation of heavy-hole spin qubits in gated semiconductor nanodevices

A novel spintronic nanodevice is proposed that is able to manipulate the single heavy-hole spin state in a coherent manner. It can act as a single quantum logic gate. The heavy-hole spin transformations are realized by transporting the hole around closed loops defined by metal gates deposited on top of the nanodevice. The device exploits Dresselhaus spin-orbit interaction, which translates the spatial motion of the hole into a rotation of the spin. The proposed quantum gate operates on subnanosecond time scales and requires only the application of a weak static voltage which allows for addressing heavy-hole spin qubits individually. Our results are supported by quantum mechanical time-dependent calculations within the four-band Luttinger-Kohn model.

Atomic scale dynamics of ultrasmall germanium clusters

Starting from the gas phase, small clusters can be produced and deposited with huge flexibility with regard to composition, materials choice and cluster size. Despite many advances in experimental characterization, a detailed morphology of such clusters is still lacking. Here we present an atomic scale observation as well as the dynamical behaviour of ultrasmall germanium clusters. Using quantitative scanning transmission electron microscopy in combination with ab initio calculations, we are able to characterize the transition between different equilibrium geometries of a germanium cluster consisting of less than 25 atoms. Seven-membered rings, trigonal prisms and some smaller subunits are identified as possible building blocks that stabilize the structure.

Electronic structure of transparent oxides with the Tran-Blaha modified Becke-Johnson potential

We present electronic band structures of transparent oxides calculated using the Tran-Blaha modified Becke-Johnson (TB-mBJ) potential. We studied the basic n-type conducting binary oxides In(2)O(3), ZnO, CdO and SnO(2) along with the p-type conducting ternary oxides delafossite CuXO(2) (X=Al, Ga, In) and spinel ZnX(2)O(4) (X=Co, Rh, Ir). The results are presented for calculated band gaps and effective electron masses. We discuss the improvements in the band gap determination using TB-mBJ compared to the standard generalized gradient approximation (GGA) in density functional theory (DFT) and also compare the electronic band structure with available results from the quasiparticle GW method. It is shown that the calculated band gaps compare well with the experimental and GW results, although the electron effective mass is generally overestimated.

Hole subbands in freestanding nanowires: six-band versus eight-band k·p modelling

The electronic structure of GaAs, InAs and InSb nanowires is studied using the six-band and the eight-band k·p models. The effect of the different Luttinger-like parameters (in the eight-band model) on the hole band structure is investigated. Although GaAs nanostructures are often treated within a six-band model because of the large bandgap, it is shown that an eight-band model is necessary for a correct description of its hole spectrum. The camel-back structure usually found in the six-band model is not always present in the eight-band model. This camel-back structure depends on the interaction between light and heavy holes, especially the ones with opposite spin. The latter effect is less pronounced in an eight-band model, but could be very sensitive to the Kane inter-band energy (E(P)) value.

Influence of an ellipsoid on the angular order in a two-dimensional cluster

The influence of an ellipsoid on the angular order of two-dimensional classical clusters is investigated through Brownian dynamics simulations. We found the following: (1) The presence of an ellipsoid does not influence the start of the angular melting, but reduces the rate at which the inner rings can rotate with respect to each other. (2) Even a small eccentricity of the ellipsoid leads to a stabilization of the angular order of the system. (3) Depending on the position of the ellipsoid in the cluster, a reentrant behavior in the angular order is observed before full radial melting of the cluster sets in. (4) The ellipsoid can lead to a two-step angular melting process: First, the rotation of the inner rings with respect to each other is hindered by the ellipsoid, but on further increasing the kinetic energy of the system, the ellipsoid just starts to behave as a spherical particle with different mobility. The effect of an ellipsoid on the molten system does not depend crucially on the interparticle interaction, but a softer parabolic confinement reduces the angular stabilization.

The quasiparticle band structure of zincblende and rocksalt ZnO

We present the quasiparticle band structure of ZnO in its zincblende (ZB) and rocksalt (RS) phases at the Γ point, calculated within the GW approximation. The effect of the p-d hybridization on the quasiparticle corrections to the band gap is discussed. We compare three systems, ZB-ZnO which shows strong p-d hybridization and has a direct band gap, RS-ZnO which is also hybridized but includes inversion symmetry and therefore has an indirect band gap, and ZB-ZnS which shows a weaker hybridization due to a change of the chemical species from oxygen to sulfur. The quasiparticle corrections are calculated with different numbers of valence electrons in the Zn pseudopotential. We find that the Zn(20+) pseudopotential is essential for the adequate treatment of the exchange interaction in the self-energy. The calculated GW band gaps are 2.47 eV and 4.27 eV respectively, for the ZB and RS phases. The ZB-ZnO band gap is underestimated compared to the experimental value of 3.27 by ∼ 0.8 eV. The RS-ZnO band gap compares well with the experimental value of 4.5 eV. The underestimation for ZB-ZnO is correlated with the strong p-d hybridization. The GW band gap for ZnS is 3.57 eV, compared to the experimental value of 3.8 eV.

Ground state configurations and melting of two-dimensional non-uniformly charged classical clusters

We consider classical two-dimensional (2D) Coulomb clusters consisting of two species containing five particles with charge q(1) and five with charge q(2), respectively. Using Monte Carlo and molecular dynamics (MD) simulations, we investigated the ground state configurations as well as radial and angular displacements of particles as a function of temperature and their dependence on the ratio q = q(2)/q(1). We found new configurations and a new multi-step melting behavior for q sufficiently different from the uniform charge limit q = 1.

Coherent three-level mixing in an electronic quantum dot

We observe magnetic-field-induced level mixing and quantum superposition phenomena between three approaching single-particle states in a quantum dot probed via the ground state of an adjacent quantum dot by single-electron resonant tunneling. The mixing is attributed to anisotropy and anharmonicity in realistic dot confining potentials. The pronounced anticrossing and transfer of strengths (both enhancement and suppression) between resonances can be understood with a simple coherent level mixing model. Superposition can lead to the formation of a dark state by complete cancellation of an otherwise strong resonance, an effect resembling coherent population trapping in a three-level-system of quantum and atom optics.

Phonon band structure of si nanowires: a stability analysis

We present full ab initio calculations of the phonon band structure of thin Si nanowires oriented along the [110] direction. Using these phonon dispersion relations, we investigate the structural stability of these wires. We found that all studied wires were stable also when doped with either B or P, if the unit cell was taken sufficiently large along the wire axis. The evolution of the phonon dispersion relations and of the sound velocities with respect to the wire diameters is discussed. Softening is observed for acoustic modes and hardening for optical phonon modes with increasing wire diameters.

Flow analyses in the lower airways: patient-specific model and boundary conditions

Computational fluid dynamics (CFD) is increasingly applied in the respiratory domain. The ability to simulate the flow through a bifurcating tubular system has increased the insight into the internal flow dynamics and the particular characteristics of respiratory flows such as secondary motions and inertial effects. The next step in the evolution is to apply the technique to patient-specific cases, in order to provide more information about pathological airways. This study presents a patient-specific approach where both the geometry and the boundary conditions (BC) are based on individual imaging methods using computed tomography (CT). The internal flow distribution of a 73-year-old female suffering from chronic obstructive pulmonary disease (COPD) is assessed. The validation is performed through the comparison of lung ventilation with gamma scintigraphy. The results show that in order to obtain agreement within the accuracy limits of the gamma scintigraphy scan, both the patient-specific geometry and the BC (driving pressure) play a crucial role. A minimal invasive test (CT scan) supplied enough information to perform an accurate CFD analysis. In the end it was possible to capture the pathological features of the respiratory system using the imaging and computational fluid dynamics techniques. This brings the introduction of this new technique in the clinical practice one step closer.

Dynamics of scattering on a classical two-dimensional artificial atom

A classical two-dimensional (2D) model for an artificial atom is used to make a numerical "exact" study of elastic and nonelastic scattering. Interesting differences in the scattering angle distribution between this model and the well-known Rutherford scattering are found in the small energy and/or small impact parameter scattering regime. For scattering off a classical 2D hydrogen atom different phenomena such as ionization, exchange of particles, and inelastic scattering can occur. A scattering regime diagram is constructed as function of the impact parameter (b) and the initial velocity (v) of the incoming particle. In a small regime of the (b,v) space the system exhibits chaos, which is studied in more detail. Analytic expressions for the scattering angle are given in the high impact parameter asymptotic limit.

Correlation between severity of sleep apnea and upper airway morphology based on advanced anatomical and functional imaging

Determination of the apnea hypopnea index (AHI) as a measure of the severity of obstructive sleep apnea/hypopnea syndrome (OSAHS) is a widely accepted methodology. However, the outcome of such a determination depends on the method used, is time consuming and insufficient for prediction of the effect of all treatment modalities. For these reasons more methods for evaluating the severity of OSAHS, based on different imaging modalities, have been looked into and recent studies have shown that anatomical properties determined from three-dimensional (3D) computed tomography (CT) images are good markers for the severity of the OSAHS. Therefore, we correlated anatomical measurements of a 3D reconstruction of the upper airway together with flow simulation results with the severity of OSAHS in order to find a combination of variables to indicate the severity of OSAHS in patients. The AHI of 20 non-selected, consecutive patients has been determined during a polysomnography. All patients also underwent a CT scan from which a 3D model of the upper airway geometry was reconstructed. This 3D model was used to evaluate the anatomical properties of the upper airway in OSAHS patients as well as to perform computational fluid dynamics (CFD) computations to evaluate the airflow and resistance of this upper airway. It has been shown that a combination of the smallest cross-sectional area and the resistance together with the body mass index (BMI) form a set of markers that predict very well the severity of OSAHS in patients within this study. We believe that these markers can be used to evaluate the outcome of an OSAHS treatment.

Inhomogeneous melting in anisotropically confined two-dimensional clusters

Molecular dynamic simulations are performed to investigate the melting process of two-dimensional clusters of classical charged particles trapped in an anisotropic parabolic potential. The confined particles interact through a repulsive potential. We find that the eccentricity of the confinement potential strongly affects the melting pattern of such clusters. Increasing the eccentricity of the confinement potential drives the system through three different melting regimes. Inhomogeneous melting is the typical melting process for anisotropically confined clusters and its appearance in small systems occurs in a distinct form called here internal intershell melting. The latter involves only particles in the center of the cluster while particles on the far left and right of the cluster are still ordered having a much higher melting temperature. Using the Lindemann's criterion the melting temperatures are determined as a function of the different parameters. The internal intershell melting process is found for both long-range (i.e., logarithmic) and short-range (i.e., screened Coulomb) interparticle interaction. Decreasing the range of the interparticle interaction increases the eccentricity of the confinement potential for which internal intershell melting can occur.

Exciton states in cylindrical nanowires

The exciton ground state and excited state energies are calculated for a model system of an infinitely long cylindrical wire. The effective Coulomb potential between the electron and the hole is studied as function of the wire radius. Within the adiabatic approximation, we obtain 'exact' numerical results for the effective exciton potential and the lowest exciton energy levels which are fitted to simple analytical expressions. Furthermore, we investigated the influence of a magnetic field parallel to the nanowire on the effective potential and the exciton energy.

Spectrum of classical two-dimensional Coulomb clusters

The frequency spectrum of a system of classical charged particles interacting through a Coulomb repulsive potential and which are confined in a two-dimensional parabolic trap is studied. It is shown that, apart from the well-known center-of-mass and breathing modes, which are independent of the number of particles in the cluster, there are more "universal" modes whose frequencies depend only slightly on the number of particles. To understand these modes the spectrum of excitations as a function of the number of particles is compared with the spectrum obtained in the hydrodynamic approach. The modes are classified according to their averaged vorticity and it is shown that these "universal" modes have the smallest vorticity and follow the hydrodynamic behavior.

Structure and spectrum of anisotropically confined two-dimensional clusters with logarithmic interaction

We studied the structural and spectral properties of a classical system consisting of a finite number of particles, moving in two dimensions, and interacting through a repulsive logarithmic potential and held together by an anisotropic harmonic potential. Increasing the anisotropy of the confinement potential can drive the system from a two-dimensional (2D) to a one-dimensional (1D) configuration. This change occurs through a sequence of structural transitions of first and second order which are reflected in the normal mode frequencies. Our results of the ground state configurations are compared with recent experiments and we obtained a satisfactory agreement. The transition from the 1D line structure to the 2D structure occurs through a zigzag transition which is of second order. We found analytical expressions for the eigenfrequencies before the zigzag transition, which allowed us to obtain an analytical expression for the anisotropy parameter at which the zigzag transition occurs as a function of the number of particles in the system.

Bubble, stripe, and ring phases in a two-dimensional cluster with competing interactions

A system of classical charged particles confined in a two-dimensional trap and interacting through a competing short-range attraction and long-range repulsion potential is studied. The structure of the system strongly depends on the interaction range of the short-range attraction potential and the total number of particles. Depending on the appropriate choice of parameters for the attractive potential, the particles organize themselves in bubbles, stripe phases, and ringlike configurations, or combinations of both of them. Detailed phase diagrams are presented for systems consisting of N=2 up to N=6 particles. General rules are derived for the different subsequent transitions between those configurations and how the ground state configuration of a large system can be deduced from the one of a small number of particles.

Multistep radial melting in small two-dimensional classical clusters

We report on a molecular dynamics study of small classical two-dimensional clusters with ringlike configurations. We focus on the particles motion at low temperatures before the radial and angular melting sets in. It is shown that in magic number configurations a local radial melting of subshells occur, which is related to the intershell rotation.

Structural phase transitions and unusual melting behavior in a classical two-dimensional Coulomb bound cluster

The melting properties of a cluster with N equally charged particles confined by a Coulomb potential are studied. The system exhibits a structural phase transition before it melts. The melting process is not dominated by the usual thermal hops between stable states. We also show that the symmetry of the ground state configuration is a dominant factor in determining the melting temperature and that more confined particles in the cluster do not necessarily have a higher melting temperature.

Nonlinear screening in large two-dimensional Coulomb clusters

The distortion due to a fixed point impurity with variable charge placed in the center of a classical harmonically confined two-dimensional (2D) large Coulomb cluster is studied. We find that the net topological charge (N(-)-N+ ) of the system is always equal to six independent of the position and charge of the impurity. In comparison with a 2D cluster without impurity charge, only the breathing mode remains unchanged. The screening length is found to be a highly nonlinear function of the impurity charge. For values of the impurity charge smaller than the charge of the other particles, the system has almost the same screening strength. When the impurity charge is larger, the screening length is strongly enhanced. This result can be explained by the competition between the different forces active in the system.

Influence of a defect particle on the structure of a classical two-dimensional cluster

A system of classical charged particles interacting through a Coulomb repulsive potential which are confined in a two-dimensional parabolic trap is studied. We allow one or two particles, called defect particles, to have a different mass and/or charge than the other particles. The structure of the whole system depends on the mass and the charge of the defects and the total number of particles in the system. The ground state configurations are investigated and phase diagrams are constructed, which explain the recent experimental results of Phys. Rev. E 64, 11 603 (2001)]]. We found that several of the experimental configurations are metastable and that replacing the Coulomb interparticle potential by an inversely quadratic one has only a minor effect on the results.

Structure and spectrum of two-dimensional clusters confined in a hard wall potential

The structural and dynamical properties of two-dimensional (2D) clusters of equally charged classical particles, which are confined in an external hard wall potential, are investigated through the Monte Carlo simulation technique. The ground-state configuration is investigated as a function of the interparticle interaction (Coulomb, dipole, logarithmic, and screened Coulomb). The excitation spectrum corresponding to the ground-state configuration of the system is discussed. The eigenmodes are investigated and the corresponding divergence and rotor are calculated, which indicates the "shearlike" and "compressionlike" aspects of the different modes. Both small and large clusters are considered.

Topological defects and nonhomogeneous melting of large two-dimensional Coulomb clusters

The configurational and melting properties of large two-dimensional (2D) clusters of charged classical particles interacting with each other via the Coulomb potential are investigated through the Monte Carlo simulation technique. The particles are confined by a harmonic potential. For a large number of particles in the cluster (N>150), the configuration is determined by two competing effects, namely, the fact that in the center a hexagonal lattice is formed, which is the groundstate for an infinite 2D system, and the confinement that imposes its circular symmetry on the outer edge. As a result, a hexagonal Wigner lattice is formed in the central area while at the border of the cluster the particles are arranged in rings. In the transition region defects appear as dislocations and disclinations at the six corners of the hexagonal-shaped inner domain. Many different arrangements and types of defects are possible as metastable configurations with a slightly higher energy. The particle motion is found to be strongly related to the topological structure. Our results clearly show that the melting of the clusters starts near the geometry induced defects, and that three different melting temperatures can be defined corresponding to the melting of different regions in the cluster.

Transition between ground state and metastable states in classical two-dimensional atoms

Structural and static properties of a classical two-dimensional system consisting of a finite number of charged particles that are laterally confined by a parabolic potential are investigated by Monte Carlo simulations and the Newton optimization technique. This system is the classical analog of the well-known quantum dot problem. The energies and configurations of the ground and all metastable states are obtained. In order to investigate the barriers and the transitions between the ground and all metastable states we first locate the saddle points between them, then by walking downhill from the saddle point to the different minima, we find the path in configurational space from the ground state to the metastable states, from which the geometric properties of the energy landscape are obtained. The sensitivity of the ground-state configuration on the functional form of the interparticle interaction and on the confinement potential is also investigated.

Molecule-type phases and Hund's rule in vertically coupled quantum dots

We study the ground state of two vertically coupled quantum dots as a function of the interdot distance within the spin density functional theory. The tunneling between the dots is included. For small and large interdot distances the atomic phases are recovered. For intermediate distances new molecule-type phases are predicted which can be observed experimentally in the addition energies. The results are interpreted in terms of an effective single particle picture and we find that Hund's rule breaks down for 11 and 12 electrons. The results are summarized in a phase diagram in which spin and isospin blockade regions are also found.