Abstract: Density-functional tight-binding and classical molecular dynamics simulations are used to investigate the structural deformations and melting of planar carbon nanoclusters C-N with N = 2-55. The minimum-energy configurations for different clusters are used as starting configurations for the study of the temperature effects on the bond breaking and rotation in carbon lines (N < 6), carbon rings (5 < N < 19), and graphene nanoflakes. The larger the rings (graphene nanoflakes) the higher the transition temperature (melting point) with ring-to-line (perfect-to-defective) transition structures. The melting point was obtained by using the bond energy, the Lindemann criteria, and the specific heat. We found that hydrogen-passivated graphene nanoflakes (CNHM) have a larger melting temperature with a much smaller dependence on size. The edges in the graphene nanoflakes exhibit several different metastable configurations (isomers) during heating before melting occurs. DOI: 10.1103/PhysRevB.87.134103

Abstract: It is well known that the formation of discrete electron levels strongly influences the pairing in metallic nanograins. Here, we focus on another effect of quantum confinement in superconducting grains that was not studied previously, i.e., spatially nonuniform pairing. This effect is very significant when single-electron levels form bunches and/or a kind of shell structure. We find that, in highly symmetric grains, the order parameter can exhibit variations with position by an order of magnitude. Nonuniform pairing is closely related to a quantum-confinement-induced modification of the pairing-interaction matrix elements and size-dependent pinning of the chemical potential to groups of degenerate or nearly degenerate levels. For illustrative purposes, we consider spherical metallic nanograins and also rectangular shapes. We show that the relevant matrix elements are, as a rule, enhanced in the presence of quantum confinement, which favors spatial variations of the order parameter, compensating the corresponding energy cost. The size-dependent pinning of the chemical potential further increases the spatial variation of the pair condensate. The role of nonuniform pairing is smaller in less symmetric confining geometries and/or in the presence of disorder. However, it always remains of importance when the energy spacing between discrete electron levels δ is approaching the scale of the bulk gap ΔB, i.e., δ>0.10.2 ΔB.

Abstract: We numerically investigate the mixed pairing symmetries inmesoscopic superconducting loops in the presence of spin correlations by solving the Bogoliubov-de Gennes equations self-consistently. The spatial variations of the superconducting order parameters and the spontaneous magnetization are determined by the band structure. When the threaded magnetic flux turns on, the charge and spin currents both emerge and depict periodic evolution. In the case of a mesoscopic loop with dominant triplet p(x) +/- ip(y)-wave symmetry, a slight change of the chemical potential may lead to novel flux-dependent evolution patterns of the ground-state energy and the magnetization. The spin-polarized currents show pronounced quantum oscillations with fractional periods due to the appearance of energy jumps in flux, accompanied with a steplike feature of the enhanced spin current. Particularly, at some appropriate flux, the peaks of the zero-energy local density of states clearly indicate the occurrence of the odd-frequency pairing. In the case of a superconducting loop with dominant singlet d(x2-y2)-wave symmetry, the spatial profiles of the zero-energy local density of states and the magnetization show spin-dependent features on different sample diagonals. Moreover, the evolution of the flux-induced spin current always exhibits an hc/e periodicity.

Abstract: The electronic tunneling properties of the two stable forms of trilayer graphene (TLG), rhombohedral ABC and Bernal ABA, are examined for p-n and p-n-p junctions as realized by using a single gate (SG) or a double gate (DG). For the rhombohedral form, due to the chirality of the electrons, the Klein paradox is found at normal incidence for SG devices, while at high-energy interband scattering between additional propagation modes can occur. The electrons in Bernal ABA TLG can have a monolayer- or bilayer-like character when incident on a SG device. Using a DG, however, both propagation modes will couple by breaking the mirror symmetry of the system, which induces intermode scattering and resonances that depend on the width of the DG p-n-p junction. For ABC TLG the DG opens up a band gap which suppresses Klein tunneling. The DG induces also an unexpected asymmetry in the tunneling angle for single-valley electrons.

Abstract: Nanoindentation of bilayer graphene is studied using molecular-dynamics simulations. We compared our simulation results with those from elasticity theory as based on the nonlinear Föppl-Hencky equations with rigid boundary condition. The force-deflection values of bilayer graphene are compared to those of monolayer graphene. Youngs modulus of bilayer graphene is estimated to be 0.8 TPa which is close to the value for graphite. Moreover, an almost flat bilayer membrane at low temperature under central load has a 14% smaller Youngs modulus as compared to the one at room temperature.

Abstract: The electric-field response of a one-dimensional ring of interacting fermions, where the interactions are described by the extended Hubbard model, is investigated. By using an accurate real-time propagation scheme based on the Chebyshev expansion of the evolution operator, we uncover various nonlinear regimes for a range of interaction parameters that allows modeling of metallic and insulating (either charge density wave or spin density wave insulators) rings. The metallic regime appears at the phase boundary between the two insulating phases and provides the opportunity to describe either weakly or strongly correlated metals. We find that the fidelity susceptibility of the ground state as a function of magnetic flux piercing the ring provides a very good measure of the short-time response. Even completely different interacting regimes behave in a similar manner at short time scales as long as the ground-state fidelity susceptibility is the same. Depending on the strength of the electric field we find various types of responses: persistent currents in the insulating phase, a dissipative regime, or damped Bloch-like oscillations with varying frequencies or even irregular in nature. Furthermore, we also consider the dimerization of the ring and describe the response of a correlated band insulator. In this case the distribution of the energy levels is more clustered and the Bloch-like oscillations become even more irregular.

Abstract: Using a numerical tight-binding approach based on the Chebyshev–Bogoliubov–de Gennes method we describe Josephson junctions made of multilayer graphene contacted by top superconducting gates. Both Bernal (ABA) and rhombohedral (ABC) stacking are considered and we find that the type of stacking has a strong effect on the proximity effect and the supercurrent flow. For both cases the pair amplitude shows a polarization between dimer and nondimer atoms, being more pronounced for rhombohedral stacking. Even though the proximity effect in nondimer sites is enhanced when compared to single-layer graphene, we find that the supercurrent is suppressed. The spatial distribution of the supercurrent shows that for Bernal stacking the current flows only in the topmost layers while for rhombohedral stacking the current flows throughout the whole structure.

Abstract: We study in detail the optoelectronic properties of graphene. Considering the electron interactions with photons and phonons, we employ the mass- and energy-balance equations to self-consistently evaluate the photoinduced carrier densities, the optical conductance, and the transmission coefficient in the presence of a linearly polarized radiation field. We demonstrate that the photoinduced carrier densities increase around the electron-photon-phonon resonant transition. They depend strongly on the radiation intensity and frequency, temperature, and dark carrier density. For short-wavelength radiation (L<3 μm), we obtain the universal optical conductance σ0=e2/(4ℏ). Importantly, there exists an optical-absorption window in the radiation wavelength range 4100 μm, which is induced by different transition energies required for interband and intraband optical absorption. The position and width of this window depend sensitively on the temperature and the carrier density of the system. These theoretical results are in line with recent experimental findings and indicate that graphene exhibits important features not only in the visible regime but also in the midinfrared bandwidth.

Abstract: A theoretical model is presented for the description of the metal-insulator transition which accompanies the structural phase transition at T approximate to 50 K in polymerized KC60. The model involves orientational charge density waves (along the C-60 polymer chains) which were introduced previously for a description of the structural phase transition. A satisfactory qualitative and quantitative understanding is obtained when the three-dimensionality of the crystal and the presence of the K+ counterions is properly taken into account.