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: We consider an interface between two strong time-reversal invariant topological insulators having surface states with opposite spin chirality, or equivalently, opposite mirror Chern number. We show that such an interface supports gapless modes that are protected by mirror symmetry. The interface states are investigated with a continuum model for the Bi2Se3 class of topological insulators that takes into account terms up to third order in the crystal momentum, which ensures that the model has the correct symmetry. The model parameters are obtained from ab initio calculations. Finally, we consider the effect of rotational mismatch at the interface, which breaks the mirror symmetry and opens a gap in the interface spectrum.

Abstract: We consider a junction between two topological insulators and calculate the properties of the interface states with an effective low-energy Hamiltonian for topological insulators with a single cone on the surface. This system bears a close resemblance to bilayer graphene, as both result from the hybridization of Dirac cones. We find gapless interface states not only when the helicity directions of the topological surface states are oppositely oriented, but they can also exist if they are equally oriented. Furthermore, we find that the existence of the interface states can be understood from the closing of the bulk gap when the helicity changes orientation. Recently superluminal tachyonic excitations were also claimed to exist at the interface between topological insulators. However, here we show that these interface states do not exist. DOI: 10.1103/PhysRevB.87.115113

Abstract: We conduct ballistic transport calculations to study the electronic transport properties of diazatriptycene molecule which can be self-assembled on metallic surfaces with uniform coverage and upright orientation of the functional head group. Due to its structural asymmetry, the molecule shows a clear current rectification, where the level of the rectification depends on the nature of the head group. For example, current rectification can be increased by more than a factor of 2 by anchoring the molecules to the electrode by CN functional group or introducing insulating CH2 group between the thiol end group and the adjacent phenyl ring. Our findings show the possibility of creating self-assembled monolayer of DAT molecules with controlled electronic transport properties through functionalization of the head group.

Abstract: We compare three first-principles methods of calculating the Curie temperature in two-dimensional (2D) ferromagnetic materials (FM), modeled using the Heisenberg model, and propose a simple formula for estimating the Curie temperature with high accuracy that works for all common 2D lattice types. First, we study the effect of exchange anisotropy on the Curie temperature calculated using the Monte Carlo (MC), the Green's function, and the renormalized spin-wave (RNSW) methods. We find that the Green's function method overestimates the Curie temperature in high-anisotropy regimes compared to the MC method, whereas the RNSW method underestimates the Curie temperature compared to the MC and the Green's function methods. Next, we propose a closed-form formula for calculating the Curie temperature of 2D FMs, which provides an estimate of the Curie temperature that is greatly improved over the mean-field expression for magnetic material screening. We apply the closed-form formula to predict the Curie temperature 2D magnets screened from the C2DB database and discover several high Curie temperature FMs, with Fe2F2 and MoI2 emerging as the most promising 2D ferromagnets. Finally, by comparing to experimental results for CrI3, CrCl3, and CrBr3, we conclude that for small effective anisotropies, the Green's-function-based equations are preferable, while for larger anisotropies, MC-based results are more predictive.

Abstract: We compare the magnetic properties of Co cluster assembled films with different degrees of oxidation. Clusters with grain size (2.3 +/- 0.7) nm are produced in a laser vaporization cluster source and soft-landed in ultrahigh vacuum conditions, forming highly porous nanogranular films. After exposure to air for different periods of time, the Co clusters oxidize and the sample may be considered as a thin antiferromagnetic Co oxide matrix containing ferromagnetic Co clusters. Magnetization measurements were performed in a temperature range from 300 down to 5 K, at applied magnetic fields up to 30 kOe. The exchange bias value at 5 K for the strongly oxidized sample is 4.8 kOe against the value of 0.75 kOe for the less oxidized sample. The mean values of the thicknesses of the Co oxide layers are estimated to be 0.6 and 0.3 nm for the more and less oxidized sample, respectively. We propose a method of measuring the exchange bias inducing temperature, i.e., the temperature at which exchange anisotropy is established. We determined the mean inducing temperatures for both samples, which are 55 and 25 K, respectively, for the more and less oxidized samples. Both temperatures are well below the bulk CoO Neel temperature of 292 K. A low value of the inducing temperature of the Co oxide layer is a consequence of its subnanometer thickness, while a large exchange bias value is a consequence of different dimensionality of Co clusters and Co oxide matrix.

Abstract: We compare classical versus quantum electron transport in recently investigated magnetic focusing devices (Bhandari et al 2016 Nano Lett. 16 1690) exposed to the perturbing potential of a scanning gate microscope (SGM). Using the Landauer-Buttiker formalism for a multi-terminal device, we calculate resistance maps that are obtained as the SGM tip is scanned over the sample. There are three unique regimes in which the scanning tip can operate (focusing, repelling, and mixed regime) which are investigated. Tip interacts mostly with electrons with cyclotron trajectories passing directly underneath it, leaving a trail of modified current density behind it. Other (indirect) trajectories become relevant when the tip is placed near the edges of the sample, and current is scattered between the tip and the edge. We point out that, in contrast to SGM experiments on gapped semiconductors, the STM tip can induce a pn junction in graphene, which improves contrast and resolution in SGM. We also discuss possible explanations for spatial asymmetry of experimentally measured resistance maps, and connect it with specific configurations of the measuring probes.

Abstract: We calculated the temperature dependence of the Debye-Waller factors for a variety of group IV, III-V and II-VI semiconductors from 0.1 to 1000 K. The approach used to fit the temperature dependence is described and resulting fit parameters are tabulated for each material. The Debye-Waller factors are deduced from generalized phonon densities of states which were derived from first principles using the WIEN2k and the ABINIT codes.

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: We calculate the quasiparticle energy spectrum of SnO2 within the GW approximation, properly taking into account the contribution of core levels to the energy corrections. The calculated fundamental gap is of 3.85 eV. We propose that the difference with respect to the experimental optical gap (3.6 eV) is due to excitonic effects in the latter. We further consider the effect applied on uniaxial pressure along the c-axis. Compared to GW, the effect of pressure on the quasiparticle energies and band gap is underestimated by the local-density approximation. The quasiparticle effective masses, however, appear to be well described by the latter.

Abstract: We calculate the excitonic spectrum of few-layer black phosphorus by direct diagonalization of the effective mass Hamiltonian in the presence of an applied in-plane electric field. The strong attractive interaction between electrons and holes in this system allows one to investigate the Stark effect up to very high ionizing fields, including also the excited states. Our results show that the band anisotropy in black phosphorus becomes evident in the direction-dependent field-induced polarizability of the exciton.

Abstract: We calculate the Curie temperature of layered ferromagnets, chromium tri-iodide (CrI3), chromium tri-bromide (CrBr3), chromium germanium tri-telluride (CrGeTe3), and the Ned temperature of a layered antiferromagnet iron di-chloride (FeCl2), using first-principles density functional theory calculations and Monte Carlo simulations. We develop a computational method to model the magnetic interactions in layered magnetic materials and calculate their critical temperature. We provide a unified method to obtain the magnetic exchange parameters (J) for an effective Heisenberg Hamiltonian from first principles, taking into account both the magnetic ansiotropy as well as the out-of-plane interactions. We obtain the magnetic phase change behavior, in particular the critical temperature, from the susceptibility and the specific-heat, calculated using the three-dimensional Monte Carlo (METROPOLIS) algorithm. The calculated Curie temperatures for ferromagnetic materials (CrI3, CrBr3, and CrGeTe3), match well with experimental values. We show that the interlayer interaction in bulk CrI3 with R (3) over bar stacking is significantly stronger than the C2/m stacking, in line with experimental observations. We show that the strong interlayer interaction in R (3) over bar CrI3 results in a competition between the in-plane and the out-of-plane magnetic easy axes. Finally, we calculate the Ned temperature of FeCl2 to be 47 +/- 8 K and show that the magnetic phase transition in FeCl2 occurs in two steps with a high-temperature intralayer ferromagnetic phase transition and a low-temperature interlayer antiferromagnetic phase transition.

Abstract: We calculate bubble nucleation rates in a Lennard-Jones fluid through explicit molecular dynamics simulations. Our approach-based on a recent free energy method (dubbed reweighted Jarzynski sampling), transition state theory, and a simple recrossing correction-allows us to probe a fairly wide range of rates in several superheated and cavitation regimes in a consistent manner. Rate predictions from this approach bridge disparate independent literature studies on the same model system. As such, we find that rate predictions based on classical nucleation theory, direct brute force molecular dynamics simulations, and seeding are consistent with our approach and one another. Published rates derived from forward flux sampling simulations are, however, found to be outliers. This study serves two purposes: First, we validate the reliability of common modeling techniques and extrapolation approaches on a paradigmatic problem in materials science and chemical physics. Second, we further test our highly generic recipe for rate calculations, and establish its applicability to nucleation processes.

Abstract: We are reporting an unexpected metal insulator transition at the ferromagnetic phase-transition temperature for thin films of La0.9Sr0.1MnO3 (< 50 nm), grown on a (100) face of SrTiO3 substrate. For the thicker films (> 50 nm), similar to the single crystal, no such transition is observed below T-C. Additionally, we observe the suppression of the features associated with charge or orbital ordering in intentionally La-deficient thin films of La0.88Sr0.1MnO3 (< 75 nm). In thin films, transmission electron microscopy reveals a compressive strain due to the epitaxial growth, that is, lattice parameters adopt those of the cubic lattice of SrTiO3. As the film thickness increases, coherent microtwinning is observed in the films and the films relax to a orthorhombic structure. (C) 2000 American Institute of Physics. [S0003-6951(00)00402-2].

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.