Abstract: Using the valence force field model of Perebeinos and Tersoff (2009 Phys. Rev. B 79 241409(R)), different energy modes of suspended graphene subjected to tensile or compressive strain are studied. By carrying out Monte Carlo simulations it is found that: (i) only for small strains (vertical bar epsilon vertical bar (sic) 0.02) is the total energy symmetrical in the strain, while it behaves completely differently beyond this threshold; (ii) the important energy contributions in stretching experiments are stretching, angle bending, an out-of-plane term, and a term that provides repulsion against pi-pi misalignment; (iii) in compressing experiments the two latter terms increase rapidly, and beyond the buckling transition stretching and bending energies are found to be constant; (iv) from stretching-compressing simulations we calculated the Young's modulus at room temperature 350 +/- 3.15 N m(-1), which is in good agreement with experimental results (340 +/- 50 N m(-1)) and with ab initio results (322-353) N m(-1); (v) molar heat capacity is estimated to be 24.64 J mol(-1) K-1 which is comparable with the Dulong-Petit value, i. e. 24.94 J mol(-1) K-1, and is almost independent of the strain; (vi) nonlinear scaling properties are obtained from height-height correlations at finite temperature; (vii) the used valence force field model results in a temperature independent bending modulus for graphene, and (viii) the Gruneisen parameter is estimated to be 0.64.

Abstract: We investigate the stability of circular monolayer graphene subjected to a radial load using non-equilibrium molecular dynamics simulations. When monolayer graphene is radially stressed, after some small circular strain (~0.4%) it buckles and bends into a new bowl-like shape. Young's modulus is calculated from the linear relation between stress and strain before the buckling threshold, which is in agreement with experimental results. The prediction of elasticity theory for the buckling threshold of a radially stressed plate is presented and its results are compared to the one of our atomistic simulation. The Jarzynski equality is used to estimate the difference between the free energy of the non-compressed states and the buckled states. From a calculation of the free energy we obtain the optimum radius for which the system feels the minimum boundary stress.

Abstract: We investigate the dynamics of a charged particle moving in a graphene layer and in a two-dimensional electron gas, where it obeys the Dirac and the Schrödinger equations, respectively. The charge carriers are described as Gaussian wavepackets. The dynamics of the wavepackets is studied numerically by solving both quantum-mechanical and relativistic equations of motion. The scattering of such wavepackets by step-like magnetic and potential barriers is analysed for different values of wavepacket energy and width. We find: (1) that the average position of the wavepacket does not coincide with the classical trajectory, and (2) that, for slanted incidence, the path of the centre of mass of the wavepacket does not have to penetrate the barrier during the scattering process. Trembling motion of the charged particle in graphene is observed in the absence of an external magnetic field and can be enhanced by a substrate-induced mass term.

Abstract: We consider a one-dimensional (1D) superlattice (SL) on graphene consisting of very high and very thin (δ-function) magnetic and potential barriers with zero average potential and zero magnetic field. We calculate the energy spectrum analytically, study it in different limiting cases, and determine the condition under which an electron beam incident on an SL is highly collimated along its direction. In the absence of the magnetic SL the collimation is very sensitive to the value of W/Ws and is optimal for W/Ws = 1, where W is the distance between the positive and negative barriers and L = W + Ws is the size of the unit cell. In the presence of only the magnetic SL the collimation decreases and the symmetry of the spectrum around ky is broken for W/W_s\neq 1 . In addition, a gap opens which depends on the strength of the magnetic field. We also investigate the effect of spatially separated potential and magnetic δ-function barriers and predict a better collimation in specific cases.

Abstract: It is well known that, in bulk, the solution of the Bogoliubovde Gennes equations is the same whether or not the HartreeFock term is included. Here the HartreeFock potential is position independent and so gives the same contribution to both the single-electron energies and the Fermi level (the chemical potential). Thus, the single-electron energies measured from the Fermi level (they control the solution) stay the same. This is not the case for nanostructured superconductors, where quantum confinement breaks the translational symmetry and results in a position-dependent HartreeFock potential. In this case its contribution to the single-electron energies depends on the relevant quantum numbers. We numerically solved the Bogoliubovde Gennes equations with the HartreeFock term for a clean superconducting nanocylinder and found a shift of the curve representing the thickness-dependent oscillations of the critical superconducting temperature to larger diameters.

Abstract: We consider classical two-dimensional (2D) Coulomb clusters consisting of two species containing five particles with charge q1 and five with charge q2, 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 = q2/q1. We found new configurations and a new multi-step melting behavior for q sufficiently different from the uniform charge limit q = 1.

Abstract: Quantum mechanical coupling and strain in two vertically arranged InP/InGaP quantum dots is studied as a function of the size of the dots and the spacer thickness. The strain distribution is determined by the continuum mechanical model, while the single-band effective-mass equation and the multiband k (.) p theory are employed to compute the conduction and valence band energy levels, respectively. The exciton states are obtained from an exact diagonalization approach, and we also compute the oscillator strength for recombination. We found that the light holes are confined by strain to the spacer, which is the reason that the hole states exhibit coupling at much larger distances as compared with the electrons. At small d, the doublet structure of the hole energy levels arises as a consequence of the relocation of the light hole from the matrix to the regions located-outside the stack, close to the dot-matrix interface. When d varies, the exciton ground state exhibits numerous anticrossings with other states, which are related to the changing spatial localization of the hole as a function of d. The oscillator strength of the exciton recombination is strongly reduced in a certain range of spacer thicknesses, which effectively turns a bright exciton state into a dark one. This effect is associated with anticrossings between exciton energy levels.

Abstract: Explicit analytic expressions are obtained for the density of states D(E) and Fermi energy E-F of a two-dimensional electron gas in the presence of a weak and periodic unidirectional electric or magnetic modulation and of a uniform perpendicular magnetic field B. The Landau levels broaden into bands and their width, proportional to the modulation strength, oscillates with B and gives rise to Weiss oscillations in D(E), E-F and the transport coefficients. When both electric and magnetic modulations are present the position of the resulting oscillations depends on the ratio delta between the two modulation strengths. When the modulations are out of phase there is no shift in the position of the oscillations when delta varies and for a particular value of delta the oscillations are suppressed.

Abstract: The dynamics of Josephson vortices (fluxons) in artificial stacks of superconducting-normal-superconducting Josephson junctions is investigated using the anisotropic time-dependent Ginzburg-Landau theory in the presence of a square/rectangular array of pinning centers (holes). For small values of the applied drive, fluxons in different junctions move out of phase, forming a periodic triangular lattice. A rectangular lattice of moving fluxons is observed at larger currents, which is in agreement with previous theoretical predictions (Koshelev and Aranson 2000 Phys. Rev. Lett. 85 3938). This 'superradiant' flux-flow state is found to be stable in a wide region of applied current. The stability range of this ordered state is considerably larger than the one obtained for the pinning-free sample. Clear commensurability features are observed in the current-voltage characteristics of the system with pronounced peaks in the critical current at (fractional) matching fields. The effect of density and strength of the pinning centers on the stability of the rectangular fluxon lattice is discussed. Predicted synchronized motion of fluxons in the presence of ordered pinning can be detected experimentally using the rf response of the system, where enhancement of the Shapiro-like steps is expected due to the synchronization.

Abstract: Using a combination of the phenomenological GinzburgLandau theory and micromagnetic simulations, we study properties of a superconducting film with an array of soft magnetic dots on top. An external in-plane magnetic field gradually drives the magnets from an out-of-plane or magnetic vortex state to an in-plane single-domain state, which changes spatially the distribution of the superconducting condensate. If induced by the magnets, the vortexantivortex molecules exhibit rich transitions as a function of the applied in-plane field. At the same time, we show how the magnetic dots act as very effective dynamic pinning centers for vortices in an applied perpendicular magnetic field.