Abstract: The autoregressive neural networks are emerging as a powerful computational tool to solve relevant problems in classical and quantum mechanics. One of their appealing functionalities is that, after they have learned a probability distribution from a dataset, they allow exact and efficient sampling of typical system configurations. Here we employ a neural autoregressive distribution estimator (NADE) to boost Markov chain Monte Carlo (MCMC) simulations of a paradigmatic classical model of spin-glass theory, namely, the two-dimensional Edwards-Anderson Hamiltonian. We show that a NADE can be trained to accurately mimic the Boltzmann distribution using unsupervised learning from system configurations generated using standard MCMC algorithms. The trained NADE is then employed as smart proposal distribution for the Metropolis-Hastings algorithm. This allows us to perform efficient MCMC simulations, which provide unbiased results even if the expectation value corresponding to the probability distribution learned by the NADE is not exact. Notably, we implement a sequential tempering procedure, whereby a NADE trained at a higher temperature is iteratively employed as proposal distribution in a MCMC simulation run at a slightly lower temperature. This allows one to efficiently simulate the spin-glass model even in the low-temperature regime, avoiding the divergent correlation times that plague MCMC simulations driven by local-update algorithms. Furthermore, we show that the NADE-driven simulations quickly sample ground-state configurations, paving the way to their future utilization to tackle binary optimization problems.

Abstract: In the quest for miniaturizing magnetic devices, the effects of boundaries and surfaces become increasingly important. Here we show how the recently predicted boundary-induced Dzyaloshinskii-Moriya interaction (DMI) affects the magnetization of ferromagnetic films with a C-infinity v symmetry and a perpendicular magnetic anisotropy. For an otherwise uniformly magnetized film, we find a surface twist when the magnetization in the bulk is canted by an in-plane external field. This twist at the surfaces caused by the boundary-induced DMI differs from the common canting caused by internal DMI observed at the edges of a chiral magnet. Furthermore, we find that the surface twist due to the boundary-induced DMI strongly affects the width of the domain wall at the surfaces. We also find that the skyrmion radius increases in the depth of the film, with the average size of the skyrmion increasing with boundary-induced DMI. This increase suggests that the boundary-induced DMI contributes to the stability of the skyrmion.

Abstract: The demand for higher throughput in the processing of materials with ultra-short pulsed lasers has motivated studies on the use of double pulses (DP). It has been observed in such studies that at relatively high time delays between the two pulses, the ablated volume is lower than that for a single pulse (SP). This has been attributed to the shielding of the second pulse and the re-deposition of the material removed by the first pulse. The investigation of re-deposition in copper with the aid of atomistic simulations is the main objective of this study. Nevertheless, a computational investigation of SP-ablation and experimental measurement of the SP-ablation depths and threshold fluence are also covered. The applied computational apparatus comprises a combination of molecular dynamics with the two-temperature model and the Helmholtz wave equation. The analysis of the simulation results shows that the derived quantities like the SP-ablation threshold fluence and the ratio of DP ablation depth to SP-ablation depth are in agreement with the experimental values. An important finding of this study is that the characteristics of the re-deposition process are highly dependent on the fluence.

Abstract: In this contribution, we report a fundamental study of the factors that set the contact resistivity between metals and highly doped n-type 2D and 3D semiconductors. We investigate the case of n-type doped Si contacted with amorphous TiSi combining first principles calculations with Non-Equilibrium Green functions transport simulations. The evolution of the intrinsic contact resistivity with the doping concentration is found to saturate at similar to 2 x 10(-10) Omega.cm(2) for the case of TiSi and imposes an intrinsic limit to the ultimate contact resistance achievable for n-doped Silamorphous-TiSi (aTiSi). The limit arises from the intrinsic properties of the semiconductors and of the metals such as their electron effective masses and Fermi energies. We illustrate that, in this regime, contacting heavy electron effective mass metals with semiconductor helps reducing the interface intrinsic contact resistivity. This observation seems to hold true regardless of the 3D character of the semiconductor, as illustrated for the case of three 2D semiconducting materials, namely MoS2, ZrS2 and HfS2. (C) The Author(s) 2018. Published by ECS.

Abstract: In this work we compute the torque field present in a ferromagnet in contact with a metallic nanowire when a skyrmion is present. If the nanowire is narrow enough, then the current is carried by a single conduction band. In this regime the classical torque model breaks down and we show that a skyrmion driven by spin transfer torque moves in a different direction than predicted by the classical model. However, the amount of charge current required to move a skyrmion with a certain velocity in the single-band regime is similar to a classical model of torque where it is implicitly assumed current transport by many conduction bands. The single-band regime is more efficient creating spin current from charge current because of the perfect polarization of the single band but is less efficient creating torque from spin current. Nevertheless, it is possible to take profit of the single-band regime to move skyrmions even with no net charge or spin current flowing between the device contacts. We have also been able to recover the classical limit considering an ensemble of only a few electronic states. In this limit we have discovered that electron diffusion needs to be considered even in ballistic nanowires due the effect of the skyrmion structure on the electron current.

Abstract: We present a tight-binding investigation of strained bilayer graphene within linear elasticity theory, focusing on the different environments experienced by the A and B carbon atoms of the different sublattices. We find that the inequivalence of the A and B atoms is enhanced by the application of perpendicular strain epsilon(zz), which provides a physical mechanism for opening a band gap, most effectively obtained when pulling the two graphene layers apart. In addition, perpendicular strain introduces electron-hole asymmetry and can result in linear electronic dispersion near the K point. Our findings suggest experimental means for strain-engineered band gaps in bilayer graphene.

Abstract: We investigate scattering of the topological surface state of a three-dimensional time-reversal invariant topological insulator when graphene is deposited on the topological-insulator surface. Specifically, we consider the (111) surface of a Bi2Se3-like topological insulator. We present a low-energy model for the graphene-topological insulator heterostructure and we calculate the transmission probability at zigzag and armchair edges of the deposited graphene, and the conductance through graphene nanoribbon barriers, and show that its features can be understood from antiresonances in the transmission probability.

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: The electronic structure of the conduction and valence bands of a quantum ring containing a layer inside the ring opening is modeled This structure (nanocup) consists of a GaAs nanodisk (the cup s bottom) and a GaAs nanoring (the cup s rim) which encircles the disk The whole system is embedded in an (Al Ga)As matrix and its shape resembles realistic ring structures grown by the droplet epitaxy technique The conduction-band states in the structure are modeled by the single-band effective-mass theory while the 4-band Luttinger-Kohn model is adopted to compute the valence-band states We analyze how the electronic structure of the nanocup evolves from the one of a quantum ring when the size of either the nanodisk or the nanoring is changed For that purpose (1) the width of the ring (2) the disk radius and (3) the disk height are separately varied For dimensions typical for experimentally realized structures we find that the electron wavefunctions are mainly localized inside the ring even when the thickness of the Inner layer is 90% of the ring thickness These calculations indicate that topological phenomena like the excitonic Aharonov-Bohm effect are negligibly affected by the presence of the layer inside the ring (C) 2010 Elsevier Ltd All rights reserved

Abstract: The water-gas-shift reaction (WGSR) is an important industrial process that can be significantly enhanced at suitable catalyst surfaces. In this work, we investigate the catalytic behavior of metallic Cu(1 0 0) and bimetallic Cu–Au(1 0 0) surfaces. With density functional theory calculations, the variation in the Gibbs free energy (ΔG°), the activation barriers, and the rate constants for the WGSR are calculated. The variation in ΔG° for water dissociation shows that the process is spontaneous up to 520 K on the bimetallic surface and up to 229 K on the Cu(1 0 0) surface. The calculated rate constants for the process also show that the bimetallic surface is much more reactive than the Cu(1 0 0) surface. The calculated pressure–temperature phase diagram for water dissociation shows that the partial pressure of H2O required for water dissociation on the bimetallic surface is substantially lower than that on the Cu(1 0 0) surface at all the studied temperatures. Additionally, the calculations demonstrate that the kinetics of the water-gas-shift reaction is dominated by redox processes on both the surfaces.

Abstract: Plasmon excitations in metallic armchair graphene nanoribbons are investigated using the random phase approximation. An exact analytical expression for the polarization function of Dirac fermions is obtained, valid for arbitrary temperature and doping. We find that at finite temperatures, due to the phase space redistribution among inter-band and intra-band electronic transitions in the conduction and valence bands, the full polarization function becomes independent of temperature and position of the chemical potential. It is shown that for a given width of nanoribbon there exists a single plasmon mode whose energy dispersion is determined by the graphene's fine structure constant. In the case of two Coulomb-coupled nanoribbons, this plasmon splits into in-phase and out-of-phase plasmon modes with splitting energy determined by the inter-ribbon spacing.

Abstract: Chrome yellow refers to a group of synthetic inorganic pigments that became popular as an artists material from the second quarter of the 19th century. The color of the pigment, in which the chromate ion acts as a chromophore, is related to its chemical composition (PbCr1-xSxO4, with 0≤x≤0.8) and crystalline structure (monoclinic/orthorhombic). Their shades range from the yellow-orange to the paler yellow tones with increasing sulfate amount. These pigments show remarkable signs of degradation after limited time periods. Pure PbCrO4 (crocoite in its natural form) has a deep yellow color and is relatively stable, while the co-precipitate with lead sulfate (PbCr1-xSxO4) has a paler shade and seems to degrade faster. This degradation is assumed to be related to the reduction of Cr(VI) to Cr(III). We show that on increasing the sulfur(S)-content in chrome yellow, the band gap increases. Typically, when increasing the band gap, one might assume that a decrease in photo activity is the result. However, the photo activity relative to the Cr content, and thus Cr reduction, of sulfur-rich PbCr1-xSxO4 is found to be much higher compared to the sulfur-poor or non-doped lead chromates. This discrepancy can be explained by the evolution of the crystal and electronic structure as function of the sulfur content: first-principles density functional theory calculations show that both the absorption coefficient and reflection coefficients of the lead chromates change as a result of the sulfate doping in such a way that the generation of electron-hole pairs under illumination relative to the total Cr content increases. These changes in the material properties explain why paler shade yellow colors of this pigment are more prone to discoloration. The electronic structure calculations also demonstrate that lead chromate and its co-precipitates are p-type semiconductors, which explains the observed reduction reaction. As understanding this phenomenon is valuable in the field of cultural heritage, this study is the first joint action of photo-electrochemical measurements and first-principles calculations to approve the higher tendency of sulfur-rich lead chromates to darken.

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 report on the observation of complex superlattices at the surface of the topological insulator Bi2Te3. Scanning tunneling microscopy reveals the existence of two different periodic structures in addition to the Bi2Te3 atomic lattice, which is found to strongly affect the local electronic structure. These three different periodicities are interpreted to result from a single small in-plane rotation of the topmost quintuple layer only. Density functional theory calculations support the observed increase in the DOS near the Fermi level, and exclude the possibility that strain is at the origin of the observed Moire pattern. Exploration of Moire superlattices formed by the quintuple layers of topological insulators holds great potential for further tuning of the properties of topological insulators.

Abstract: Thin superconducting films can exhibit negative magnetoresistance when an in-plane external magnetic field is aligned parallel with the transport current. We explain this effect as due to appearance of parallel vortices in the plain of the film at the first critical magnetic field H-c1 which leads to an enhancement of the superconducting properties and impedes the motion of the current induced perpendicular vortices. Our theoretical results are based on a numerical solution of the time-dependent and stationary 3D Ginzburg-Landau equations.