Abstract: We performed a detailed investigation of the oxidation of the technologically relevant GaAs(001)-beta 2(2x4) surface via density functional calculations. The purpose is to gain insights on the atomistic mechanisms and local bondings that underlie the degradation of the surface properties once exposed to oxygen. The study comprises the adsorption of single O atoms, through the sampling of several adsorption sites, and the subsequent formation of the O adsorbate at increasing coverage by taking into account multiple-atom adsorption. Based on the evaluation of the energetics and the structural properties of the atomistic models generated, the results here reported delineate a consistent picture of the initial stage of the surface oxidation: (i) at low coverage, in the limit of single O insertions, oxygen is incorporated on the surface forming a twofold-bridging Ga-O-As bond; (ii) at increasing coverage, as multiple O atoms are involved, this is accompanied by the formation of a threefold-coordinated bond (with two Ga and one As atoms); (iii) the latter has important implications regarding the electronic properties of the adsorbate since this O bonding may result in the formation of As dangling bonds. Moreover, a clear trend of increased energy gain for the incorporation of neighboring O atoms compared to single O insertions indicates that the formation of oxide clusters is favored over a regime of uniform oxidation. Our findings provide a detailed description of the O bonding and stress the importance of modeling the adsorption of multiple O atoms for an accurate description of the surface oxidation.

Abstract: A data-acquisition and data-processing method is proposed that aims at minimizing the effect of retardation on the Kramers-Kronig analysis of valence electron energy-loss spectra. This method is applied to diamond, which, due to its high dielectric constant, is a material that shows strong retardation effects and thus is a challenging material to be studied by valence electron energy-loss spectroscopy. The results obtained show a significant improvement but still show small discrepancies with respect to optical data, which are most likely due to the residual retardation contributions and the fact that nonzero momentum transfers are measured.

Abstract: Using atomistic simulations we investigate the morphological properties of graphene deposited on top of a nanostructured substrate. Sinusoidally corrugated surfaces, steps, elongated trenches, one-dimensional and cubic barriers, spherical bubbles, Gaussian bumps, and Gaussian depressions are considered as support structures for graphene. The graphene-substrate interaction is governed by van der Waals forces and the profile of the graphene layer is determined by minimizing the energy using molecular dynamics simulations. Based on the obtained optimum configurations, we found that (i) for graphene placed over sinusoidally corrugated substrates with corrugation wavelengths longer than 2 nm, the graphene sheet follows the substrate pattern while for supported graphene it is always suspended across the peaks of the substrate, (ii) the conformation of graphene to the substrate topography is enhanced when increasing the energy parameter in the van der Waals model, (iii) the adhesion of graphene into the trenches depends on the width of the trench and on the graphene's orientation, i. e., in contrast to a small-width (3 nm) nanoribbon with armchair edges, the one with zigzag edges follows the substrate profile, (iv) atomic-scale graphene follows a Gaussian bump substrate but not the substrate with a Gaussian depression, and (v) the adhesion energy due to van der Waals interaction varies in the range [0.1-0.4] J/m(2).

Abstract: The interaction between electrons and plasmons in trilayer graphene is investigated within the Overhauser approach resulting in the “plasmaron” quasiparticle. This interaction is cast into a field theoretical problem, and its effect on the energy spectrum is calculated using improved Wigner-Brillouin perturbation theory. The plasmaron spectrum is shifted with respect to the bare electron spectrum by ΔE(k)∼150−200meV for ABC stacked trilayer graphene and for ABA trilayer graphene by ΔE(k)∼30−150 meV[ ΔE(k) ∼1 −5meV] for the hyperbolic (linear) part of the spectrum. The shift in general increases with the electron concentration and electron momentum. The dispersion of plasmarons is more pronounced in ABC stacked than in ABA stacked trilayer graphene, because of the different energy band structure and their different plasmon dispersion.

Abstract: The split diffuse maxima around the {110} and {100} positions in the diffraction pattern of short-range-ordered Cu1-xPdx alloys (x=0.10...0.60) are attributed to small atomic clusters, being part of the underlying fee lattice. By analyzing the reciprocal space geometry, our cluster method identifies two prominent cluster types: the tetrahedron of nearest neighbors and a linear three-points cluster along the [110] directions. Since both cluster types contain different information on the same nearest-neighbor correlations, local anisotropy has to be assumed. It is shown that the three interatomic pair interactions within these basic clusters are sufficient to generate the spot splitting in the diffraction pattern. A ground-state analysis with these interactions reproduces the results of the anisotropic next-nearest-neighbor Ising model.

Abstract: Vortex loops are generated by the inhomogeneous stray field of a magnetic dipole on top of a current-carrying mesoscopic superconductor. Cutting and recombination processes unfold under the applied drive, resulting in periodic voltage oscillations across the sample. We show that a direct and detectable consequence of the cutting and recombination of these vortex loops in the present setup is the onset of vortices at surfaces where they were absent prior to the application of the external current. The nonlinear dynamics of vortex loops is studied within the time-dependent Ginzburg-Landau theory to describe the profound three-dimensional features of their time evolution.

Abstract: We study magnetic flux interacting with arrays of pinning sites (APSs) placed on vertices of hyperbolic tesselations (HTs). We show that, due to the gradient in the density of pinning sites, HT APSs are capable of trapping vortices for a broad range of applied magnetic fluxes. Thus, the penetration of magnetic field in HT APSs is essentially different from the usual scenario predicted by the Bean model. We demonstrate that, due to the enhanced asymmetry of the surface barrier for vortex entry and exit, this HT APS could be used as a “capacitor” to store magnetic flux.

Abstract: We study the influence of the inelastic relaxation time τ̃E of the quasiparticle distribution function f(E) on the phase slip process in quasi-one-dimensional superconducting rings at a temperature close to the critical temperature Tc. We find that the initial time of growth of the order parameter |Δ| in the phase slip core after the phase slip is a nonmonotonic function of τ̃E which has a maximum at τ̃E≃τ̃GL=πℏ/8kB(Tc−T) and has a tendency to saturate for large τ̃E⪢τ̃GL. The effective heating of the electron subsystem due to the increase in |Δ| in the phase slip center together with the above effect result in a nonmonotonic dependence of the number of subsequent phase slips on τ̃E in rings of relatively large radius (in which each phase slip reduces the current density to a small fraction of its initial value). During the phase slip process the order parameter distribution has two peaks near the phase slip core due to the diffusion of the nonequilibrium quasiparticles from that region.

Abstract: Quantum confinement of the perpendicular motion of electrons in single-crystalline metallic superconducting nanofilms splits the conduction band into a series of single-electron subbands. A distinctive feature of such a nanoscale multi-band superconductor is that the energetic position of each subband can vary significantly with changing nanofilm thickness, substrate material, protective cover and other details of the fabrication process. It can occur that the bottom of one of the available subbands is situated in the vicinity of the Fermi level. We demonstrate that the character of the superconducting pairing in such a subband changes dramatically and exhibits a clear molecule-like trend, which is very similar to the well-known crossover from the Bardeen-Cooper-Schrieffer regime to Bose-Einstein condensation (BCS-BEC) observed in trapped ultracold fermions. For Pb nanofilms with thicknesses of 4 and 5 monolayers (MLs) this will lead to a spectacular scenario: up to half of all the Cooper pairs nearly collapse, shrinking in the lateral size (parallel to the nanofilm) down to a few nanometers. As a result, the superconducting condensate will be a coherent mixture of almost molecule-like fermionic pairs with ordinary, extended Cooper pairs.

Abstract: Nonlinear dynamics of Josephson vortices (fluxons) in artificial stacks of superconducting-normal-superconducting Josephson junctions under simultaneously applied time-periodic ac and constant biasing dc currents is studied using the time dependent Ginzburg-Landau formalism with a Lawrence-Doniach extension. At zero external magnetic field and dc biasing current the resistive state of the system is characterized by periodic nucleation and annihilation of fluxon-antifluxon pairs, relative positions of which are determined by the state of neighboring junctions. Due to the mutual repulsive interaction, fluxons in different junctions move out of phase. Their collective motion can be synchronized by adding a small ac component to the biasing dc current. Coherent motion of fluxons is observed for a broad frequency range of the applied drive. In the coherent state the maximal output voltage, which is proportional to the number of junctions in the stack, is observed near the characteristic frequency of the system determined by the crossing of the fluxons across the sample. However, in this frequency range the dynamically synchronized state has an alternative-a less ordered state with smaller amplitude of the output voltage. Collective behavior of the junctions is strongly affected by the sloped sidewalls of the stack. Synchronization is observed only for weakly trapezoidal cross sections, whereas irregular motion of fluxons is observed for larger slopes of the sample edge.

Abstract: Using atomistic simulations we investigate the thermodynamical properties of a single atomic layer of hexagonal boron nitride (h-BN). The thermal induced ripples, heat capacity, and thermal lattice expansion of large scale h-BN sheets are determined and compared to those found for graphene (GE) for temperatures up to 1000 K. By analyzing the mean-square height fluctuations < h(2)> and the height-height correlation function H(q) we found that the h-BN sheet is a less stiff material as compared to graphene. The bending rigidity of h-BN (i) is about 16% smaller than the one of GE at room temperature (300 K), and (ii) increases with temperature as in GE. The difference in stiffness between h-BN and GE results in unequal responses to external uniaxial and shear stress and different buckling transitions. In contrast to a GE sheet, the buckling transition of a h-BN sheet depends strongly on the direction of the applied compression. The molar heat capacity, thermal-expansion coefficient, and Gruneisen parameter are estimated to be 25.2 J mol(-1) K-1, 7.2 x 10(-6) K-1, and 0.89, respectively.

Abstract: Crystal structure and magnetic properties of the layered vanadium phosphate PbZnVO(PO4)2 are studied using x-ray powder diffraction, magnetization and specific-heat measurements, as well as band-structure calculations. The compound resembles AA′VO(PO4)2 vanadium phosphates and fits to the extended frustrated square-lattice model with the couplings J1, J1′ between nearest neighbors and J2, J2′ between next-nearest neighbors. The temperature dependence of the magnetization yields estimates of averaged nearest-neighbor and next-nearest-neighbor couplings, J̅ 1≃−5.2 K and J̅ 2≃10.0 K, respectively. The effective frustration ratio α=J̅ 2/J̅ 1 amounts to −1.9 and suggests columnar antiferromagnetic ordering in PbZnVO(PO4)2. Specific-heat data support the estimates of J̅ 1 and J̅ 2 and indicate a likely magnetic ordering transition at 3.9 K. However, the averaged couplings underestimate the saturation field, thus pointing to the spatial anisotropy of the nearest-neighbor interactions. Band-structure calculations confirm the identification of ferromagnetic J1, J1′ and antiferromagnetic J2, J2′ in PbZnVO(PO4)2 and yield (J1′−J1)≃1.1 K in excellent agreement with the experimental value of 1.1 K, deduced from the difference between the expected and experimentally measured saturation fields. Based on the comparison of layered vanadium phosphates with different metal cations, we show that a moderate spatial anisotropy of the frustrated square lattice has minor influence on the thermodynamic properties of the model. We discuss relevant geometrical parameters, controlling the exchange interactions in these compounds and propose a strategy for further design of strongly frustrated square-lattice materials.