Physical Chemistry Chemical Physics

Golub, Pavlo;Manzhos, Sergei

doi: 10.1039/C8CP06433Dpmid: 30525136

Abstract: We study the performance of fourth-order gradient expansions of the kinetic energy density (KED) in semi-local kinetic energy functionals depending on the density-dependent variables. The formal fourth-order expansion is convergent for periodic systems and small molecules but does not improve over the second-order expansion (Thomas-Fermi term plus one-ninth of von Weizsäcker term). Linear fitting of the expansion coefficients somewhat improves on the formal expansion. The tuning of the fourth order expansion coefficients allows for better reproducibility of Kohn-Sham kinetic energy density than the tuning of the second-order expansion coefficients alone. The possibility of a much more accurate match with the Kohn-Sham kinetic energy density by using neural networks trained using the terms of the 4th order expansion as density-dependent variables is demonstrated. We obtain ultra-low fitting errors without overfitting. Small single hidden layer neural networks can provide good accuracy in separate KED fits of each compound, while for joint fitting of KEDs of multiple compounds multiple hidden layers were required to achieve good fit quality. The critical issue of data distribution is highlighted. We also show the critical role of pseudopotentials in the performance of the expansion, where in the case of a too rapid decay of the valence density at the nucleus with some pseudopotentials, numeric instabilities arise.

Xia, Yi-qi;Shan, Wen-jie;Tian, Wen-de;Chen, Kang;Ma, Yu-qiang

doi: 10.1039/C8CP05976Dpmid: 30734786

Abstract: Folding and unfolding of biopolymers are often manipulated in experiment by tuning pH, temperature, single-molecule force or shear field. Here we carry out Brownian dynamics simulations to explore the behavior of a single self-attracting chain in the suspension of self-propelling particles (SPPs). As the propelling force increases, globule-stretch (G-S) transition of the chain happens due to the enhanced disturbance from SPPs. Two distinct mechanisms of the transition in the limits of low and high rotational diffusion rates of SPPs have been observed: shear effect at low rate and collision-induced melting at high rate. The G-S and S-G (stretch-globule) curves form hysteresis loop at low rate, while they merge at high rate. Besides, we find two competing effects result in the non-monotonic dependence of the G-S transition on the SPP density at low rate. Our results suggest an alternative approach to manipulating the folding and unfolding of (bio)polymers by utilizing active agents.

Rossi, Kevin;Asara, Gian Giacomo;Baletto, Francesca

doi: 10.1039/C8CP05720Fpmid: 30484450

Abstract: With a focus on platinum nanoparticles of different sizes (diameter of 1-9 nm) and shapes, we sequence their geometrical genome by recording the relative occurrence of all the non equivalent active site, classified according to the number of neighbours in their first and second coordination shell. The occurrence of sites is morphology and size dependent, with significant changes in the relative occurrence up to 9 nm. Our geometrical genome sequencing approach is immediately transferable to address the effects of the morphological polydispersivity in size-selected samples and the influence of temperature, including ionic vibrations and thermal activated processes. The proposed geometrical genome forecasts an enhancement of the catalytic reduction of molecular oxygen on stellated and anisotropic platinum twinned nanoparticles, with their shortest axes of ~2 nm, and an irreversible disruption of the Pt nanocatalyst's structure above 1000 K.

Kolesov, Grigory;Lin, Chungwei;Knyazev, Andrew;Kojima, Keisuke;Katz, Joseph;Akiyama, Koichi;Nakai, Eiji;Kawahara, Hiroyuki

doi: 10.1039/C9CP00122Kpmid: 30839041

Abstract: Band alignment between two materials is of fundamental importance for multitude of applications. However, density functional theory (DFT) either underestimates the bandgap - as is the case with local density approximation (LDA) or generalized gradient approximation (GGA) - or is highly computationally demanding, as is the case with hybrid-functional methods. The latter can become prohibitive in electronic-structure calculations of supercells which describe quantum wells. We propose to apply the DFT$+U$ method, with $U$ for each atomic shell being treated as set of tuning parameters, to automatically fit the bulk bandgap and the lattice constant, and then use thus obtained $U$ parameters in large supercell calculations to determine the band alignment. We apply this procedure to InP/In$_{0.5}$Ga$_{0.5}$As, In$_{0.5}$Ga$_{0.5}$As/In$_{0.5}$Al$_{0.5}$As and InP/In$_{0.5}$Al$_{0.5}$As quantum wells, and obtain good agreement with experimental results. Although this procedure requires some experimental input, it provides both meaningful valence and conduction band offsets while, crucially, lattice relaxation is taken into account. The computational cost of this procedure is comparable to that of LDA. We believe that this is a practical procedure that can be useful for providing accurate estimate of band alignments between more complicated alloys.

Borghetti, P.;Mouchaal, Y.;Dai, Z.;Cabailh, G.;Chenot, S.;Lazzari, R.;Jupille, J.

doi: 10.1039/C6CP08595Dpmid: 28379222

Abstract: Orientation-dependent reactivity and band-bending are evidenced upon Ti deposition (1-10 Å) on the polar ZnO(0001)-Zn and ZnO(000$\bar{1}$)-O surfaces. At the onset of the Ti deposition, a downward band-bending was observed on ZnO(000$\bar{1}$)-O while no change occurred on ZnO(0001)-Zn. Combining this with the photoemission analysis of the Ti 2p core level and Zn L$_3$(L$_2$)M$_{45}$M$_{45}$ Auger transition, it is established that the Ti/ZnO reaction is of the form Ti + 2 ZnO $\rightarrow$ TiO$_2$ + 2 Zn on ZnO(0001)-Zn and Ti + y ZnO $\rightarrow$ TiZn$_x$O$_y$ + (y-x) Zn on ZnO(000$\bar{1}$)-O. Consistently, upon annealing thicker Ti adlayers, the metallic zinc is removed to leave ZnO(0001)-Zn surfaces covered with TiO$_2$-like phase and ZnO(000$\bar{1}$)-O surfaces covered with a defined (Ti, Zn, O) compound. Finally, a difference in the activation temperature between the O-terminated (500 K) and Zn-terminated (700 K) surfaces is observed, which is tentatively explained by different electric fields in the space charge layer at ZnO surfaces.

Cai, Jin;Xie, Yuee;Chang, Po-Yao;Kim, Heung-Sik;Chen, Yuanping

doi: 10.1039/C8CP02810Apmid: 30083674

Abstract: Coexistence of topological elements in a topological metal/semimetal (TM) has gradually attracted attentions. However, the non-topological factors always mess up the Fermi surface and cover interesting topological properties. Here, we find that Ba3Si4 is a "clean" TM in which coexists nodal-chain network, intersecting nodal rings (INRs) and triple points, in the absence of spin-orbit coupling (SOC). Moreover, the nodal rings in the topological phase exhibit diverse types: from type-I, type-II to type-III rings according to band dispersions. All the topological elements are generated by crossings of three energy bands, and thus they are correlated rather than mutual independence. When some structural symmetries are eliminated by an external strain, the topological phase evolves into another phase including Hopf link, one-dimensional nodal chain and new INRs.

de Lima, F. Crasto;Ferreira, G. J.;Miwa, R. H.

doi: 10.1039/C8CP03792Bpmid: 30132483

Abstract: The realization of the Quantum anomalous Hall effect (QAHE) in two dimensional (2D) metal organic frameworks (MOFs), (MC$_4$S$_4$)$_3$ with M = Mn, Fe, Co, Ru and Rh, has been investigated based on a combination of first-principles calculations and tight binding models. Our results for the magnetic anisotropy energy (MAE) reveal that the out-of-plane (in-plane) magnetization is favored for M = Mn, Fe, and Ru (Co, and Rh). Given the structural symmetry of (MC$_4$S$_4$)$_3$, the QAHE takes place only for M = Mn, Fe and Ru. Such a quantum anomalous Hall phase has been confirmed through the calculation of the Chern number, and examining the formation of topologically protected (metallic) edge states. Further electron ($n$-type) doping of the MOFs has been done in order to place the Fermi level within the non-trivial energy gap; where we find that in (RuC$_4$S$_4$)$_3$, in addition to the up-shift of the Fermi level, the MAE energy increases by 40\%. Finally, we show that in MOF/graphene (vdW) interfaces, the Fermi level tunning can be done with an external electric field, which controls the charge transfer at the MOF/graphene interface, giving rise to switchable topologically protected edge currents in MOFs.

Shukla, Vivekanand;Grigoriev, Anton;Jena, Naresh K.;Ahuja, Rajeev

doi: 10.1039/C8CP03815Epmid: 30156222

Abstract: Experimental realization of two-dimensional boron sheets was reported very recently by Feng et. al. using molecular beam epitaxy on silver (111) surface. These boron sheets possess promising electronic and transport properties. We performed the density functional theory (DFT) calculation to see the stability of two $\beta_{12}$ and $\chi$ polymorphs of boron under strain and further studied electronic and transport properties. We verified the directional dependency in electron transport properties in these two boron polymorphs. Here we report tunable anisotropic transport properties of the borophenes. We also investigated current-voltage characteristics in low bias regime after applying strain on these systems to see how this external strain affects the anisotropy of current.

Webster, Lucas;Liang, Liangbo;Yan, Jia-An

doi: 10.1039/C8CP03599Gpmid: 30183024

Abstract: We apply the density-functional theory to study various phases (including non-magnetic (NM), anti-ferromagnetic (AFM), and ferromagnetic (FM)) in monolayer magnetic chromium triiodide (CrI$_3$), a recently fabricated 2D magnetic material. It is found that: (1) the introduction of magnetism in monolayer CrI$_3$ gives rise to metal-to-semiconductor transition; (2) the electronic band topologies as well as the nature of direct and indirect band gaps in either AFM or FM phases exhibit delicate dependence on the magnetic ordering and spin-orbit coupling; and (3)the phonon modes involving Cr atoms are particularly sensitive to the magnetic ordering, highlighting distinct spin-lattice and spin-phonon coupling in this magnet. First-principles simulations of the Raman spectra demonstrate that both frequencies and intensities of the Raman peaks strongly depend on the magnetic ordering. The polarization dependent $A_{1g}$ modes at 77 cm$^{-1}$ and 130 cm$^{-1}$ along with the $E_g$ mode at about 50 cm$^{-1}$ in the FM phase may offer a useful fingerprint to characterize this material. Our results not only provide a detailed guiding map for experimental characterization of CrI$_3$, but also reveal how the evolution of magnetism can be tracked by its lattice dynamics and Raman response.

Schmitteckert, Peter

doi: 10.1039/C8CP03763Apmid: 30371691

Abstract: In this work we discuss the extraction of mean field single particle Hamiltonians from a many body wave function of a fermionic system. It allows us to discuss the result of a many particle wave function in terms of a non-interacting description.