Abstract: The lattice dynamics of low-frequency rigid-plane modes in metallic (graphene multilayers, GML) and in insulating (hexagonal boron-nitride multilayers, BNML) layered crystals is investigated. The frequencies of shearing and compression (stretching) modes depend on the layer number N and are presented in the form of fan diagrams. The results for GML and BNML are very similar. In both cases, only the interactions (van der Waals and Coulomb) between nearest-neighbor planes are effective, while the interactions between more distant planes are screened. A comparison with recent Raman scattering results on low-frequency shear modes in GML [Tan et al., Nat. Mater., in press, doi: 10.1038/nmat3245, (2012)] is made. Relations with the low-lying rigid-plane phonon dispersions in the bulk materials are established. Master curves, which connect the fan diagram frequencies for any given N, are derived. Static and dynamic thermal correlation functions for rigid-layer shear and compression modes are calculated. The results might be of use for the interpretation of friction force experiments on multilayer crystals.

Abstract: Using highly efficient simulations of the tight-binding Bogoliubov-de-Gennes model, we solved self-consistently for the pair correlation and the Josephson current in a superconducting-bilayer graphene-superconducting Josephson junction. Different doping levels for the non-superconducting link are considered in the short- and long-junction regimes. Self-consistent results for the pair correlation and superconducting current resemble those reported previously for single-layer graphene except at the Dirac point, where remarkable differences in the proximity effect are found, as well as a suppression of the superconducting current in the long-junction regime. Inversion symmetry is broken by considering a potential difference between the layers and we found that the supercurrent can be switched if the junction length is larger than the Fermi length.

Abstract: Monte Carlo (MC) methods have a long-standing history as partners of molecular dynamics (MD) to simulate the evolution of materials at the atomic scale. Among these techniques, the uniform-acceptance force-bias Monte Carlo (UFMC) method [ G. Dereli Mol. Simul. 8 351 (1992)] has recently attracted attention [ M. Timonova et al. Phys. Rev. B 81 144107 (2010)] thanks to its apparent capacity of being able to simulate physical processes in a reduced number of iterations compared to classical MD methods. The origin of this efficiency remains, however, unclear. In this work we derive a UFMC method starting from basic thermodynamic principles, which leads to an intuitive and unambiguous formalism. The approach includes a statistically relevant time step per Monte Carlo iteration, showing a significant speed-up compared to MD simulations. This time-stamped force-bias Monte Carlo (tfMC) formalism is tested on both simple one-dimensional and three-dimensional systems. Both test-cases give excellent results in agreement with analytical solutions and literature reports. The inclusion of a time scale, the simplicity of the method, and the enhancement of the time step compared to classical MD methods make this method very appealing for studying the dynamics of many-particle systems.

Abstract: We consider two possible mechanisms for single-vortex fluctuative entry/exit through the surface barrier in ultrathin superconducting disk-shaped nanoislands made of Pb and consisting of just a few monoatomic layers, which can be fabricated using modern techniques. We estimate tunneling probabilities and establish criteria for the crossover between these two mechanisms depending on magnetic field and system sizes. For the case of vortex entry, quantum tunneling dominates on the major part of the temperature/flux phase diagram. For the case of vortex exit, thermal activation turns out to be more probable. This nontrivial result is due to the subtle balance between the barrier height and width, which determine rates of the thermal activation and quantum tunneling, respectively.

Abstract: We study the structure and phonon spectrum of a two-dimensional bilayer system of classical charged dipoles oriented perpendicular to the plane of the layers for equal density in each layer. This system can be tuned through six different crystalline phases by changing the interlayer separation or the charge and/or dipole moment of the particle. The presence of the charge on the dipole particles is responsible for the nucleation of five staggered phases and a disordered phase which are not found in the magnetic dipole bilayer system. These extra phases are a consequence of the competition between the repulsive Coulomb and the attractive dipole interlayer interaction. We present the phase diagram and determine the order of the phase transitions. The phonon spectrum of the system was calculated within the harmonic approximation, and a nonmonotonic behavior of the phonon spectrum is found as a function of the effective strength of the interparticle interaction. The stability of the different phases is determined.

Abstract: We report on the dynamical properties of interacting colloids confined to one dimension and subjected to external periodic energy landscapes. We particularly focus on the influence of hydrodynamic interactions on the mean-square displacement. Using Brownian dynamics simulations, we study colloidal systems with two types of repulsive interparticle interactions, namely, Yukawa and superparamagnetic potentials. We find that in the homogeneous case, hydrodynamic interactions lead to an enhancement of the particle mobility and the mean-square displacement at long times scales as t(alpha), with alpha = 1/2 + epsilon and epsilon being a small correction. This correction, however, becomes much more important in the presence of an external field, which breaks the homogeneity of the particle distribution along the line and, therefore, promotes a richer dynamical scenario due to the hydrodynamical coupling among particles. We provide here the complete dynamical scenario in terms of the external potential parameters: amplitude and commensurability.

Abstract: We show in the case of a superconducting Nb ladder that a mesoscopic superconductor typically exhibits magnetoresistance oscillations whose amplitude and temperature dependence are different from those stemming from the Little-Parks effect. We demonstrate that these large resistance oscillations (as well as the monotonic background on which they are superimposed) are due to current-excited moving vortices, where the applied current in competition with the oscillating Meissner currents imposes or removes the barriers for vortex motion in an increasing magnetic field. Because of the ever present current in transport measurements, this effect should be considered in parallel with the Little-Parks effect in low-critical temperature (T-c) samples, as well as with recently proposed thermal activation of dissipative vortex-antivortex pairs in high-T-c samples.

Abstract: A mode converter for electron vortex beams is described. Numerical simulations, confirmed by experiment, show that the converter transforms a vortex beam with a topological charge m = +/- 1 into beams closely resembling Hermite-Gaussian HG(10) and HG(01) modes. The converter can be used as a mode discriminator or filter for electron vortex beams. Combining the converter with a phase plate turns a plane wave into modes with topological charge m = +/- 1. This combination serves as a generator of electron vortex beams of high brilliance.

Abstract: 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.

Abstract: Vortex matter in mesoscopic superconductors is known to be strongly affected by the geometry of the sample. Here we show that in nanoscale superconductors with coherence length comparable to the Fermi wavelength the shape resonances of the order parameter results in an additional contribution to the quantum topological confinement-leading to unconventional vortex configurations. Our Bogoliubov-de Gennes calculations in a square geometry reveal a plethora of asymmetric, giant multivortex, and vortex-antivortex structures, stable over a wide range of parameters and which are very different from those predicted by the Ginzburg-Landau theory. These unconventional states are relevant for high-T-c nanograins, confined Bose-Einstein condensates, and graphene flakes with proximity-induced superconductivity.

Abstract: The appearance of massless Dirac fermions in graphene requires two equivalent carbon sublattices of trigonal shape. While the generation of an effective mass and a band gap at the Dirac point remains an unresolved problem for freestanding extended graphene, it is well established by breaking translational symmetry by confinement and by breaking sublattice symmetry by interaction with a substrate. One of the strongest sublattice-symmetry-breaking interactions with predicted and measured band gaps ranging from 400 meV to more than 3 eV has been attributed to the interfaces of graphene with Ni and Co, which are also promising spin-filter interfaces. Here, we apply angle-resolved photoemission to epitaxial graphene on Ni (111) and Co(0001) to show the presence of intact Dirac cones 2.8 eV below the Fermi level. Our results challenge the common belief that the breaking of sublattice symmetry by a substrate and the opening of the band gap at the Dirac energy are in a straightforward relation. A simple effective model of a biased bilayer structure composed of graphene and a sublattice-symmetry-broken layer, corroborated by density-functional-theory calculations, demonstrates the general validity of our conclusions.

Abstract: A numerical model for NxOy remediation in humid air plasma produced with a dielectric barrier discharge at atmospheric pressure is presented. Special emphasis is given to NO2 and N2O reduction with the decrease of O2 content in the feedstock gas. A detailed reaction mechanism including electronic and ionic processes, as well as the contribution of radicals and excited atomic/molecular species is proposed. The temporal evolution of the densities of NO, NO2 and N2O species, and some other by-products, is analyzed, and the major pathways for the NxOy remediation are discussed for one pulse. Subsequently, simulations are presented for a multi-pulses case, where three O2 contents are tested for optimization of the remediation process. It is found that when the gas mixture O2/N2/H2O/NOx has no initial O2 content, the best NOx and N2O remediation is achieved.

Abstract: Modeling results are presented to gain a better insight in the properties of a SiH4/O2/Ar inductively coupled plasma (ICP) and how it interacts with a silicon substrate (wafer), as applied in the microelectronics industry for the fabrication of electronic devices. The SiH4/O2/Ar ICP is used for the filling of microtrenches with isolating material (SiO2), as applied in shallow trench isolation (STI). In this article, a detailed reaction set that describes the plasma chemistry of SiH4/O2/Ar discharges as well as surface processes, such as sputtering, oxidation, and deposition, is presented. Results are presented on the plasma properties during the plasma enhanced chemical vapor deposition process (PECVD) for different gas ratios, as well as on the shape of the filled trenches and the surface compositions of the deposited layers. For the operating conditions under study it is found that the most important species accounting for deposition are SiH2, SiH3O, SiH3 and SiH2O, while SiH+2, SiH+3, O+2 and Ar+ are the dominant species for sputtering of the surface. By diluting the precursor gas (SiH4) in the mixture, the deposition rate versus sputtering rate can be controlled for a desired trench filling process. From the calculation results it is clear that a high deposition rate will result in undesired void formation during the trench filling, while a small deposition rate will result in undesired trench bottom and mask damage by sputtering. By varying the SiH4/O2 ratio, the chemical composition of the deposited layer will be influenced. However, even at the highest SiH4/O2 ratio investigated (i.e., 3.2:1; low oxygen content), the bulk deposited layer consists mainly of SiO2, suggesting that low-volatile silane species deposit first and subsequently become oxidized instead of being oxidized first in the plasma before deposition. Finally, it was found that the top surface of the deposited layer contained less oxygen due to preferential sputtering of O atoms, making the top layer more Si-rich. However, this effect is negligible at a SiH4/O2 ratio of 2:1 or lower.

Abstract: Small stationary diesel engines, like in generator sets, have limited emission control measures and are therefore responsible for 44% of the particulate matter (PM) emissions in the United States. The diesel exhaust composition depends on operating conditions of the combustion engine. Furthermore, the measurements are influenced by the used sampling method. This study examines the effect of engine loading and exhaust gas dilution on the composition of small-scale power generators. These generators are used in different operating conditions than road-transport vehicles, resulting in different emission characteristics. Experimental data were obtained for gaseous volatile organic compounds (VOC) and PM mass concentration, elemental composition and nitrate content. The exhaust composition depends on load condition because of its effect on fuel consumption, engine wear and combustion temperature. Higher load conditions result in lower PM concentration and sharper edged particles with larger aerodynamic diameters. A positive correlation with load condition was found for K, Ca, Sr, Mn, Cu, Zn and Pb adsorbed on PM, elements that originate from lubricating oil or engine corrosion. The nitrate concentration decreases at higher load conditions, due to enhanced nitrate dissociation to gaseous NO at higher engine temperatures. Dilution on the other hand decreases PM and nitrate concentration and increases gaseous VOC and adsorbed metal content. In conclusion, these data show that operating and sampling conditions have a major effect on the exhaust gas composition of small-scale diesel generators. Therefore, care must be taken when designing new experiments or comparing literature results.