Abstract: We apply neutron diffraction, high-resolution synchrotron x-ray diffraction, magnetization measurements, electronic structure calculations, and quantum Monte-Carlo simulations to unravel the structure and magnetism of (CuCl)LaTa2O7. Despite the pseudo-tetragonal crystallographic unit cell, this compound features an orthorhombic superstructure, similar to the Nb-containing (CuX)LaNb2O7 with X = Cl and Br. The spin lattice entails dimers formed by the antiferromagnetic fourth-neighbor coupling J(4), as well as a large number of nonequivalent interdimer couplings quantified by an effective exchange parameter J(eff). In (CuCl)LaTa2O7, the interdimer couplings are sufficiently strong to induce the long-range magnetic order with the Neel temperature T-N similar or equal to 7 K and the ordered magnetic moment of 0.53 mu(B), as measured with neutron diffraction. This magnetic behavior can be accounted for by J(eff)/J(4) similar or equal to 1.6 and J(4) similar or equal to 16 K. We further propose a general magnetic phase diagram for the (CuCl)LaNb2O7-type compounds, and explain the transition from the gapped spin-singlet (dimer) ground state in (CuCl)LaNb2O7 to the long-range antiferromagnetic order in (CuCl)LaTa2O7 and (CuBr)LaNb2O7 by an increase in the magnitude of the interdimer couplings J(eff)/J(4), with the (CuCl)LaM2O7 (M = Nb, Ta) compounds lying on different sides of the quantum critical point that separates the singlet and long-range-ordered magnetic ground states.

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.

Abstract: Carefully converged calculations are performed for the band gap of ZnO within many-body perturbation theory (G0W0 approximation). The results obtained using four different well-established plasmon-pole models are compared with those of explicit calculations without such models (the contour-deformation approach). This comparison shows that, surprisingly, plasmon-pole models depending on the f-sum rule gives less precise results. In particular, it confirms that the band gap of ZnO is underestimated in the G0W0 approach as compared to experiment, contrary to the recent claim of Shih et al. [ Phys. Rev. Lett. 105 146401 (2010)].

Abstract: The giant isotope effect on the ferro- and antiferroelectric transition temperature upon deuteration of hydrogen-bonded systems is well known experimentally since various decades. Yet, theoretically only recently a microscopic understanding of this effect has been achieved which, specifically, took into account the geometry of the O ... H ... O bond. The implications of this modeling are multiple as numerous hydrogen-bonded organic systems show the same effects as ferro- and antiferroelectrics, i.e., cooperative proton tunneling at a well-defined temperature.

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: The thermodynamic properties of the electron gas in multilayer graphene depend strongly on the number of layers and the type of stacking. Here we analyze how those properties change when we vary the number of layers for rhombohedral stacked multilayer graphene and compare our results with those from a conventional two-dimensional electron gas. We show that the highly degenerate zero-energy Landau level which is partly filled with electrons and partly with holes has a strong influence on the values of the different thermodynamic quantities.

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: 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: We calculate the spectrum of states that are localized at the interface between two regions of opposite bias in trilayer graphene (TLG). These potential profiles, also known as potential kinks, have been predicted to support two different branches of localized states for the case of bilayer graphene, and show similarities to the surface states of topological insulators. On the other hand, we found that ABC stacked TLG exhibits three different unidimensional branches of states in each valley that are confined to the kink interface. They have the property E(k(y)) = -E(-k(y)) when belonging to the same valley and E-K(k(y)) = -E-K' (-k(y)). A kink-antikink potential profile opens a gap in the spectrum of these one-dimensional states.

Abstract: Low-temperature magnetoconductance measurements were made in the vicinity of the charge neutrality point (CNP). Two origins for the fluctuations were identified close to the CNP. At very low magnetic fields there exist only mesoscopic magnetoconductance quantum interference features which develop rapidly as a function of density. At slightly higher fields (>0.5 T), close to the CNP, additional fluctuations track the quantum Hall (QH) sequence expected for monolayer graphene. These additional features are attributed to effects of locally charging individual QH localized states. These effects reveal a precursor to the quantum Hall effect since, unlike previous transport observations of QH dot charging effects, they occur in the absence of quantum Hall plateaus or Shubnikov-de Haas oscillations. From our transport data we are able to extract parameters that characterize the inhomogeneities in our device.