Abstract: Large-scale atomistic simulations using the reactive force field approach are implemented to investigate the thermomechanical properties of fluorinated graphene (FG). A set of parameters for the reactive force field potential optimized to reproduce key quantum mechanical properties of relevant carbon-fluorine cluster systems are presented. Molecular dynamics simulations are used to investigate the thermal rippling behavior of FG and its mechanical properties and compare them with graphene, graphane and a sheet of boron nitride. The mean square value of the height fluctuations < h(2)> and the height-height correlation function H(q) for different system sizes and temperatures show that FG is an unrippled system in contrast to the thermal rippling behavior of graphene. The effective Young's modulus of a flake of fluorinated graphene is obtained to be 273 N/m and 250 N/m for a flake of FG under uniaxial strain along armchair and zigzag directions, respectively. DOI: 10.1103/PhysRevB.87.104114

Abstract: Within the Ginzburg-Landau theory, we study the vortex structures in three-dimensional anisotropic mesoscopic superconductors in the presence of a uniform magnetic field. Anisotropy is included through varied Tc in different layers of the sample and leads to distinct differences in the vortex states and their free energy. Several unconventional states are found, some comprising vortex clusters or exhibiting asymmetry. In a tilted magnetic field, we found second-order transitions between different vortex states, although vortex entry is generally a first-order transition in mesoscopic samples. In multilayered samples the kinked vortex strings are formed owing to the competing interactions of vortices with Meissner currents and the weak-link boundaries. The length and deformation of vortex fragments are determined solely by the inclination and strength of applied magnetic field, and this lock-in does not depend on the degree of anisotropy between the superconducting layers.

Abstract: The technique of beam optics is applied to the description of the wave function of Dirac electrons. This approach is illustrated by considering electron transmission through simple nonhomogeneous structures, such as flat and bent p-n junctions and superlattices. We found that a convex p-n junction compresses the beam waist, while a concave interface widens it without loosing its focusing properties. At a flat p-n junction the waist of the transmitted Gaussian beam can be narrowed or widened, depending on the angle of incidence. A general condition is derived for the occurrence of beam collimation in a superlattice which is less stringent than previous discussed.

Abstract: We investigate the effect of Rashba and intrinsic spin-orbit couplings on the electronic properties and spin configurations of Dirac fermions confined in: (i) a flat graphene sheet, (ii) a graphene wire with p-n-p structure, and (iii) a superlattice of graphene wires. The interplay between the spin-orbit interaction mechanisms breaks the electron-hole symmetry and the spin configuration induced by Rashba spin-orbit coupling lacks inversion symmetry in k space. We show that the Rashba spin-orbit interaction doubles the Fabry-Perot resonant modes in the transmission spectrum of a graphene wire and opens new channels for the electron transmission. Moreover, it leads to the appearance of spin split extra Dirac cones in the energy spectrum of a graphene superlattice. It is shown that the spin of the electrons and holes confined in a flat graphene sheet is always perpendicular to their motion while this is not the case for the other nanostructures.

Abstract: Full ab initio techniques are applied to study the electronic and dynamical properties of free standing, hydrogen-passivated Si/Ge core-shell nanowires oriented along the [110] direction. All studied wires exhibit a direct band gap and are found to be structurally stable. The different contributions of the core and shell atoms to the phonon spectra are identified. The acoustic phonon velocities and the frequencies of some typical optical modes are compared with those of pure Si and Ge nanowires. These depend either on the concentration or on the type of core material. Optical modes are hardened and longitudinal acoustic velocities are softened with decreasing wire diameter.

Abstract: Si nanowires with a ⟨111⟩ orientation, synthesized by vapor-liquid-solid process with low silane partial pressure reactant and gold as the catalyst, are known to exhibit sawtooth facets containing gold adsorbates. We report herein the study of the nanowire morphology by means of transmission electron microscopy and scanning tunneling microscopy. The nanowires consist of faceted sidewalls. The number of the sidewalls changes from 12 to 6 along the growth axis, giving rise to nanowires with an irregular hexagonal cross section at their base. The sidewalls are covered with Au-rich clusters. Their facets also exhibit atomic structures that reveal the presence of gold, resulting from the diffusion of gold during the growth. Based on these observations, the tapering of the nanowire is found to be related to two contributions: the reduction in the catalyst particle size during the growth and lateral overgrowth from the direct incorporation of Si species onto the nanowire sidewalls. Because the rearrangement of atoms at surfaces and interfaces might affect the growth kinetics, the trigonal symmetry as well as the higher lateral growth rate on the widest sidewalls are explained from the existence of an interfacial atomic structure with two inequivalent parts in the unit cell. Finally, spectroscopic measurements were performed on the major facets and revealed a metallic behavior at 77 K.

Abstract: Rashba-type spin-orbit interaction (SOI) controlled band structure of a two-dimensional superlattice allows for the modulation of the conductance of finite size devices by changing the strength of the SOI. We consider rectangular arrays and find that the temperature dependence of the conductance disappears for high temperatures, but the strength of the SOI still affects the conductance at these temperatures. The modulation effect can be seen even in the presence of strong dephasing, which can be important for practical applications.

Abstract: An interface modification of (LaCa)MnO3-BaTiO3 superlattices was found to massively influence magnetic and magnetotransport properties. Moreover it determines the crystal structure of the manganite layers, changing it from orthorhombic (Pnma) for the conventional superlattice (cSL), to rhombohedral (R3̅ c) for the modified one (mSL). While the cSL shows extremely nonlinear ac transport, the mSL is an electrically homogeneous material. The observations go beyond an oversimplified picture of dead interface layers and evidence the importance of electronic correlations at perovskite interfaces.

Abstract: A unified theory of phonon dispersions and piezoelectricity in bulk and multilayers of hexagonal boron nitride (h-BN) is derived. The dynamical matrix is calculated on the basis of an empirical force constant model of intralayer valence and interlayer van der Waals interactions. Coulomb interactions are calculated by Ewalds method, adapted for the three-dimensional (3D) and the multilayer case. The deformation of the ionic charge distribution with long-wave lattice displacements is taken into account. Special attention is devoted to the nonanalytic long-range Coulomb contribution to the dynamical matrix which is different for the 3D crystal and the multilayer case. Consequently there is a splitting of the transverse optical (TO) and longitudinal optical (LO) phonon branches of E1u symmetry and a discontinuity of the A2u branch at the Γ point in 3D h-BN. No such splitting and discontinuity at Γ are present in multilayer crystals with a finite number N of layers. There a diverging bundle of N overbending optical phonon branches emerges from Γ. Borns long-wave theory is applied and extended for the study of piezoelectricity in layered crystals. While 3D h-BN and h-BN multilayers with an even number of layers (symmetry D6h) are not piezoelectric, multilayers with an uneven number of Nu layers (symmetry D3h) are piezoelectric; the piezoelectric coefficient e1,11 is inversely proportional to Nu.

Abstract: Electrons in graphene, in addition to their spin, have two pseudospin degrees of freedom: sublattice and valley pseudospin. Valleytronics uses the valley degree of freedom as a carrier of information similarly to the way spintronics uses electron spin. We show how a double-barrier structure consisting of electric and vector potentials can be used to filter massive Dirac electrons based on their valley index. We study the resonant transmission through a finite number of barriers and we obtain the energy spectrum of a superlattice consisting of electric and vector potentials. When a mass term is included, the energy bands and energy gaps at the K and K′ points are different and they can be tuned by changing the potential.

Abstract: The effect of different stacking order of graphene multilayers on the electric field induced band gap is investigated. We considered a positively charged top and a negatively charged back gate in order to independently tune the band gap and the Fermi energy of three and four layer graphene systems. A tight-binding approach within a self-consistent Hartree approximation is used to calculate the induced charges on the different graphene layers. We found that the gap for trilayer graphene with the ABC stacking is much larger than the corresponding gap for the ABA trilayer. Also we predict that for four layers of graphene the energy gap strongly depends on the choice of stacking, and we found that the gap for the different types of stacking is much larger as compared to the case of Bernal stacking. Trigonal warping changes the size of the induced electronic gap by approximately 30% for intermediate and large values of the induced electron density.

Abstract: The ground state energy and the extent of the wavefunction of a negatively charged donor (D − ) located near a semiconductormetal or a semiconductordielectric interface are obtained. We apply the effective mass approximation and use a variational two-electron wavefunction that takes into account the influence of all image charges that arise due to the presence of the interface, as well as the correlation between the two electrons bound to the donor. For a semiconductormetal interface, the D − binding energy is enhanced for donor positions d > 1.5aB (aB is the effective Bohr radius) due to the additional attraction of the electrons with their images. When the donor approaches the interface (i.e. d < 1.5aB) the D − binding energy drops and eventually it becomes unbound. For a semiconductordielectric (or a semiconductorvacuum) interface the D − binding energy is reduced for any donor position as compared to the bulk case and the system becomes rapidly unbound when the donor approaches the interface.