Abstract: Results of dielectric and resonance frequency (f(r)) measurements below room temperature are presented for Pb(ZrxTi1-x)O-3, x = 0.515 and 0.520. It is shown that the temperature coefficient of f(r) changes sign from negative to positive around 210 and 265 K for x = 0.520 and 200 and 260 K for x = 0.515. Anomalies in the real part of the dielectric constant (epsilon') are observed around the same temperatures at which the temperature coefficient of f(r) changes sign because of the electrostrictive coupling between the elastic and dielectric responses. Low-temperature powder x-ray-diffraction (XRD) data, however, reveal only one transition from the tetragonal to monoclinic phase similar to that reported by Noheda et al. [Phys. Rev. B, 61, 8687 (2000)]. Electron-diffraction data, on the other hand, reveal yet another structural transition at lower temperatures corresponding to the second anomaly in the epsilon' vs T and f(r) vs T curves. This second transition is shown to be a cell-doubling transition not observed by Noheda et al. in their XRD studies. The observation of superlattice reflections raises doubts about the correctness of the Cm space group proposed by Noheda et al. for the monoclinic phase of Pb(ZrxTi(1-x))O-3 below the second transition temperature.

Abstract: Quasiperiodic pinning arrays, as recently demonstrated theoretically and experimentally using a fivefold Penrose tiling, can lead to a significant enhancement of the critical current Ic as compared to traditional regular pinning arrays. However, while regular arrays showed only a sharp peak in Ic(Φ) at the matching flux Φ1 and quasiperiodic arrays provided a much broader maximum at Φ<Φ1, both types of pinning arrays turned out to be inefficient for fluxes larger than Φ1. We demonstrate theoretically and experimentally the enhancement of Ic(Φ) for Φ>Φ1 by using non-Penrose quasiperiodic pinning arrays. This result is based on a qualitatively different mechanism of flux pinning by quasiperiodic pinning arrays and could be potentially useful for applications in superconducting microelectronic devices operating in a broad range of magnetic fields.

Abstract: The effect of strain in graphene is usually modeled by a pseudomagnetic vector potential which is, however, derived in the limit of small strain. In realistic cases deviations are expected in view of graphene's very high strain tolerance, which can be up to 25%. Here we investigate the pseudomagnetic field generated by a Gaussian bump and we show that it exhibits significant differences with numerical tight-binding results. Furthermore, we calculate the electronic states in the strained region for a hexagon shaped flake with armchair edges. We find that the sixfold symmetry of the wave functions inside the Gaussian bump is directly related to the different effects of strain along the fundamental directions of graphene: zigzag and armchair. Low energy electrons are strongly confined in the armchair directions and are localized on the carbon atoms of a single sublattice.

Abstract: When subjected to the interfacially induced Dzyaloshinskii-Moriya interaction, the ground state in thin ferromagnetic films with high perpendicular anisotropy is cycloidal. The period of this cycloidal state depends on the strength of the Dzyaloshinskii-Moriya interaction. In this work, we have studied the effect of confinement on the magnetic ground state and excited states, and we determined the phase diagram of thin strips and thin square platelets by means of micromagnetic calculations. We show that multiple cycloidal states with different periods can be stable in laterally confined films, where the period of the cycloids does not depend solely on the Dzyaloshinskii-Moriya interaction strength but also on the dimensions of the film. The more complex states comprising skyrmions are also found to be stable, though with higher energy.

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: Using highly efficient GPU-based simulations of the tight-binding Bogoliubov-de Gennes equations we solve self-consistently for the pair correlation in rhombohedral (ABC) and Bernal (ABA) multilayer graphene by considering a finite intrinsic s-wave pairing potential. We find that the two different stacking configurations have opposite bulk/surface behavior for the order parameter. Surface superconductivity is robust for ABC stacked multilayer graphene even at very low pairing potentials for which the bulk order parameter vanishes, in agreement with a recent analytical approach. In contrast, for Bernal stacked multilayer graphene, we find that the order parameter is always suppressed at the surface and that there exists a critical value for the pairing potential below which no superconducting order is achieved. We considered different doping scenarios and find that homogeneous doping strongly suppresses surface superconductivity while nonhomogeneous field-induced doping has a much weaker effect on the superconducting order parameter. For multilayer structures with hybrid stacking (ABC and ABA) we find that when the thickness of each region is small (few layers), high-temperature surface superconductivity survives throughout the bulk due to the proximity effect between ABC/ABA interfaces where the order parameter is enhanced. DOI: 10.1103/PhysRevB.87.134509

Abstract: The ground-state properties of a paramagnetic Mott insulator at half-filling are investigated in the presence of an external electric field using the inhomogeneous Gutzwiller approximation for a single-band Hubbard model in a slab geometry. We find that the metal-insulator transition is shifted toward higher Hubbard repulsions by applying an electric field perpendicular to the slab. The main reason is the accumulation of charges near the surface. The spatial distribution of site-dependent quasiparticle weight shows that it is maximal in a few layers beneath the surface, while the central sites where the field is screened have a very low quasiparticle weight. Our results show that above a critical-field value, states near the surface will be metallic, while the bulk quasiparticle weight is extremely suppressed but never vanishing, even for large Hubbard repulsions above the bulk zero-field critical value. Below the critical-field value, our results hint toward an insulating state in which the electric field is totally screened and the slab is again at half-filling.

Abstract: Recent experiments showed that nonuniform strain can be produced by depositing graphene over pillars. We employed atomistic calculations to study the nonuniform strain and the induced pseudomagnetic field in graphene on top of nanopillars. By decreasing the distance between the nanopillars a complex distribution for the pseudomagnetic field can be generated. Furthermore, we performed tight-binding calculations of the local density of states (LDOS) by using the relaxed graphene configuration obtained from atomistic calculations. We find that the quasiparticle LDOS are strongly modified near the pillars, both at low energies showing sublattice polarization and at high energies showing shifts of the van Hove singularity. Our study shows that changing the specific pattern of the nanopillars allows us to create a desired shape of the pseudomagnetic field profile while the LDOS maps provide an input for experimental verification by scanning tunneling microscopy.

Abstract: Atomistic simulations are used to study the bending of rectangular graphene nanoribbons subjected to axial stress both for free boundary and supported boundary conditions. The shapes of the deformations of the buckled graphene nanoribbons, for small values of the stress, are sine waves where the number of nodal lines depend on the longitudinal size of the system and the applied boundary condition. The buckling strain for the supported boundary condition is found to be independent of the longitudinal size and estimated to be 0.86%. From a calculation of the free energy at finite temperature we find that the equilibrium projected two-dimensional area of the graphene nanoribbon is less than the area of a flat sheet. At the optimum length the boundary strain for the supported boundary condition is 0.48%.