Elucidating Solvatochromic Shifts in Two-Dimensional Photocatalysts by Solving the Bethe-Salpeter Equation Coupled with Implicit Solvation Method.

Many studies have focused on tailoring the photophysical properties of two-dimensional (2D) materials for photocatalytic (PC) or photoelectrochemical (PEC) applications. To understand the optical properties of 2D materials in solution, we established a computational method that combined the Bethe-Salpeter equation (BSE) calculations with our GW-GPE method, allowing for GW/BSE-level calculations with implicit solvation described using the generalized Poisson equation (GPE). We applied this method to MoS2, phosphorene (PP), and g-C3N4 and found that when the solvent dielectric increased, it reduced the exciton binding energy and quasiparticle bandgap, resulting in almost no solvatochromic shift in the excitonic peaks of MoS2 and PP, which is consistent with previous experiments. However, our calculations predicted that the solvent dielectric had a significant impact on the excitonic properties of g-C3N4, exhibiting a large solvatochromic shift. We expect that our GW/BSE-GPE method will offer insights into the design of 2D materials for PC and PEC applications.

Gold(I)-Thiolate Coordination Polymers as Multifunctional Materials: The Case of Au(I)-p-Fluorothiophenolate.

Gold-sulfur interaction has vital importance in nanotechnologies and material chemistry to design functional nanoparticles, self-assembled monolayers, or molecular complexes. In this paper, a mixture of only two basic precursors, such as the chloroauric acid (HAu(III)Cl4) and a thiol molecule (p-fluorothiophenol (p-HSPhF)), are used for the synthesis of gold(I)-thiolate coordination polymers. Under different conditions of synthesis and external stimuli, five different functional materials with different states of [Au(I)(p-SPhF)]n can be afforded. These gold-thiolate compounds are (i) red emissive, flexible, and crystalline fibers; (ii) composite materials made of these red emissive fibers and gold nanoparticles; (iii) amorphous phase; (iv) transparent glass; and (v) amorphous-to-crystalline phase-change material associated with an ON/OFF switch of luminescence. The different functionalities of these materials highlight the great versatility of the gold(I) thiolate coordination polymers with easy synthesis and diverse shaping that may have great potential as sustainable phosphors, smart textiles, sensors, and phase change memories.

Tuning the electronic and optical properties of small organic acenedithiophene molecular crystals for photovoltaic applications: First principles calculations.

Periodic density functional theory was employed to investigate the impact of chemical modifications on the properties of π-conjugated acenedithiophene molecular crystals. Here, we highlight the importance of the ß-methylthionation effect, the position of the sulfur atoms of the thiacycle group and their size, and the number of central benzene rings in the chemical modification strategy. Our results show that the introduction of the methylthio groups at the ß-positions of the thiophene and the additional benzene ring at the center of the BDT crystal structure are a promising strategy to improve the performance of organic semiconductors, as observed experimentally. We found that ß-MT-ADT exhibits large charge carrier mobility, which is in good agreement with the experimental results and comparable to that of rubrene. In addition, the electronic and optical properties of these ambipolar materials suggest promising performances with ß-MT-ADT > ADT >ß-MT-NDT > NDT > BEDT-BDT >ß-MT-BDT > BDT. Moreover, functionalization with thiacycle-fused sulfur atoms of different sizes and numbers improve the properties of BDT but is still less efficient than the methylthionation effect. Overall, our findings suggest a promising molecular modification strategy for possibly high performance ambipolar organic semiconducting materials.

Tuning the 1D-2D dimensionality upon ligand exchange in silver thiolate coordination polymers with photoemission switching.

Silver nanoparticles are known and widely used for their antimicrobial activities. Nevertheless, once they are released into the natural or biological environments, they can become toxic with time, because of the dissolution of some Ag(I) ions that can then react with thiol-based molecules, such as glutathione and/or compete with copper proteins. These assumptions are based on the high affinity of the soft acid Ag(I) and the soft base thiolates and the exchange reactions that are involved in complex physiological media. Here we synthesized and fully characterized two new 2D silver thiolate coordination polymers (CPs) that exhibit a reversible 2D-to-1D structural transformation in the presence of an excess of thiol molecules. This dimensionality change induces also a switch of the yellow emission of the Ag-thiolate CP. This study highlights that these highly stable silver-thiolate CPs, in basic, acidic and oxidant media can undergo a complete dissolution-recrystallization mechanism upon thiol exchange reactions.

Solving the Schrödinger Equation in the Configuration Space with Generative Machine Learning.

The configuration interaction approach provides a conceptually simple and powerful approach to solve the Schrödinger equation for realistic molecules and materials but is characterized by an unfavorable scaling, which strongly limits its practical applicability. Effectively selecting only the configurations that actually contribute to the wave function is a fundamental step toward practical applications. We propose a machine learning approach that iteratively trains a generative model to preferentially generate the important configurations. By considering molecular applications it is shown that convergence to chemical accuracy can be achieved much more rapidly with respect to random sampling or the Monte Carlo configuration interaction method. This work paves the way to a broader use of generative models to solve the electronic structure problem.

First-principles study of the structural and electronic properties of tetragonal ZrOX (X = S, Se, and Te) monolayers and their vdW heterostructures for applications in optoelectronics and photocatalysis.

Using first-principles calculations, we have studied the structural and electronic properties of ZrOX (X = S, Se, and Te) monolayers and their van der Waals heterostructures in the tetragonal structure. Our results show that these monolayers are dynamically stable and are semiconductors with electronic bandgaps ranging from 1.98 to 3.16 eV as obtained with the GW approximation. By computing their band edges, we show that ZrOS and ZrOSe are of interest for water splitting applications. In addition, the van der Waals heterostructures formed by these monolayers show a type I band alignment for ZrOTe/ZrOSe and a type II alignment for the other two heterostructures, making them potential candidates for certain optoelectronic applications involving electron/hole separation.

Ab initio molecular dynamics investigation of the co-adsorption of iodine species with CO and H2O in silver-exchanged chabazite.

In the field of nuclear energy, there is particular interest for the trapping of harmful iodine species (I2 and CH3I) that could be released during a nuclear accident, due to their dangereous impact on the human metabolic processes and the environment. Here, the adsorption of these iodine molecules versus several inhibitory compounds (CO, H2O, CH3Cl and Cl2) in the silver exchanged chabazite zeolite is studied in detail using ab initio molecular dynamics simulations at a realistic temperature and composition. Interestingly, we found that the iodine molecules remain attached to the cations even when the number of water molecules inside the structure is greater than two times the number of cations per cell at T = 413 K. For CO, we found that CH3I is more perturbed than I2 by the presence of this inhibitor. Overall, our results indicate that the silver-exchanged chabazite zeolite is a promising candidate to trap iodine species in the case of a severe nuclear accident.

GW Quasiparticle Energies and Bandgaps of Two-Dimensional Materials Immersed in Water.

Computational simulations have become of major interest to screen potential photocatalysts for optimal band edge positions which straddle the redox potentials. Unfortunately, these methods suffer from a difficulty in resolving the dynamic solvent response on the band edge positions. We have developed a computational method based on the GW approximation coupled with an implicit solvation model that solves a generalized Poisson equation (GPE), that is, GW-GPE. Using GW-GPE, we have investigated the band edge locations of (quasi) 2D materials immersed in water and found a good agreement with experimental data. We identify two contributions of the solvent effect, termed a "polarization-field effect" and an "environmental screening effect", which are found to be highly sensitive to the atomic and charge distribution of the 2D materials. We believe that the GW-GPE scheme can pave the way to predict band edge positions in solvents, enabling design of 2D material-based photocatalysts and energy systems.

Assessing the Accuracy of Machine Learning Thermodynamic Perturbation Theory: Density Functional Theory and Beyond.

Machine learning thermodynamic perturbation theory (MLPT) is a promising approach to compute finite temperature properties when the goal is to compare several different levels of ab initio theory and/or to apply highly expensive computational methods. Indeed, starting from a production molecular dynamics trajectory, this method can estimate properties at one or more target levels of theory from only a small number of additional fixed-geometry calculations, which are used to train a machine learning model. However, as MLPT is based on thermodynamic perturbation theory (TPT), inaccuracies might arise when the starting point trajectory samples a configurational space which has a small overlap with that of the target approximations of interest. By considering case studies of molecules adsorbed in zeolites and several different density functional theory approximations, in this work we assess the accuracy of MLPT for ensemble total energies and enthalpies of adsorption. It is shown that problematic cases can be detected even without knowing reference results and that even in these situations it is possible to recover target level results within chemical accuracy by applying a machine-learning-based Monte Carlo (MLMC) resampling. Finally, on the basis of the ideas developed in this work, we assess and confirm the accuracy of recently published MLPT-based enthalpies of adsorption at the random phase approximation level, whose high computational cost would completely hinder a direct molecular dynamics simulation.

Establishing the accuracy of density functional approaches for the description of noncovalent interactions in ionic liquids.

We test a number of dispersion corrected versatile Generalized Gradient Approximation (GGA) and meta-GGA functionals for their ability to predict the interactions of ionic liquids, and show that most can achieve energies within 1 kcal mol-1 of benchmarks. This compares favorably with an accurate dispersion corrected hybrid, ωB97X-V. Our tests also reveal that PBE (Perdew-Burke-Ernzerhof GGA) calculations using the plane-wave projector augmented wave method and Gaussian Type Orbitals (GTOs) differ by less than 0.6 kJ mol-1 for ionic liquids, despite ions being difficult to evaluate in periodic cells - thus revealing that GTO benchmarks may be used also for plane-wave codes. Finally, the relatively high success of explicit van der Waals density functionals, compared to elemental and ionic dispersion models, suggests that improvements are required for low-cost dispersion correction models of ions.

Structure and Energetics of Dye-Sensitized NiO Interfaces in Water from Ab Initio MD and Large-Scale GW Calculations.

The energy-level alignment across solvated molecule/semiconductor interfaces is a crucial property for the correct functioning of dye-sensitized photoelectrodes, where, following the absorption of solar light, a cascade of interfacial hole/electron transfer processes has to efficiently take place. In light of the difficulty of performing X-ray photoelectron spectroscopy measurements at the molecule/solvent/metal-oxide interface, being able to accurately predict the level alignment by first-principles calculations on realistic structural models would represent an important step toward the optimization of the device. In this respect, dye/NiO surfaces, employed in p-type dye-sensitized solar cells, are undoubtedly challenging for ab initio methods and, also for this reason, much less investigated than the n-type dye/TiO2 counterpart. Here, we consider the C343-sensitized NiO surface in water and combine ab initio molecular dynamics (AIMD) simulations with GW (G0W0) calculations, performed along the MD trajectory to reliably describe the structure and energetics of the interface when explicit solvation and finite temperature effects are accounted for. We show that the differential perturbative correction on the NiO and molecule states obtained at the GW level is mandatory to recover the correct (physical) interfacial energetics, allowing hole transfer from the semiconductor valence band to the highest occupied molecular orbital (HOMO) of the dye. Moreover, the calculated average driving force quantitatively agrees with the experimental estimate.

Assessment and prediction of band edge locations of nitrides using a self-consistent hybrid functional.

Due to their optimal bandgap size and large defect tolerance, nitrides are becoming pivotal materials in several optoelectronic devices, photovoltaics, and photocatalysts. A computational method that can accurately predict their electronic structures is indispensable for exploring new nitride materials. However, the relatively small bandgap of nitrides, which stems from the subtle balance between ionic and covalent bond characteristics, makes conventional density functional theory challenging to achieve satisfactory accuracy. Here, we employed a self-consistent hybrid functional where the Hartree-Fock mixing parameter is self-consistently determined and thus the empiricism of the hybrid functional is effectively removed to calculate the bandgaps of various nitride compounds. By comparing the bandgaps from the self-consistent hybrid functional calculations with the available experimental and high-level GW calculation results, we found that the self-consistent hybrid functional can provide a computationally efficient approach for quantitative predictions of nitride electronic structures with an accuracy level comparable to the GW method. Additionally, we aligned the band edge positions of various nitride compounds using self-consistent hybrid functional calculations, providing material design principles for heterostructures of nitride-based optoelectronic devices. We anticipate the wide use of the self-consistent hybrid functional for accelerating explorations and predictions of new nitride-based functional materials in various photoactive applications.

Hybrid localized graph kernel for machine learning energy-related properties of molecules and solids.

Nowadays, the coupling of electronic structure and machine learning techniques serves as a powerful tool to predict chemical and physical properties of a broad range of systems. With the aim of improving the accuracy of predictions, a large number of representations for molecules and solids for machine learning applications has been developed. In this work we propose a novel descriptor based on the notion of molecular graph. While graphs are largely employed in classification problems in cheminformatics or bioinformatics, they are not often used in regression problem, especially of energy-related properties. Our method is based on a local decomposition of atomic environments and on the hybridization of two kernel functions: a graph kernel contribution that describes the chemical pattern and a Coulomb label contribution that encodes finer details of the local geometry. The accuracy of this new kernel method in energy predictions of molecular and condensed phase systems is demonstrated by considering the popular QM7 and BA10 datasets. These examples show that the hybrid localized graph kernel outperforms traditional approaches such as, for example, the smooth overlap of atomic positions and the Coulomb matrices.

Ab Initio Study of Helicity-Dependent Light-Induced Demagnetization: From the Optical Regime to the Extreme Ultraviolet Regime.

We use ab initio real-time time-dependent density functional theory to investigate the effect of optical and extreme ultraviolet (XUV) circularly polarized femtosecond pulses on the magnetization dynamics of ferromagnetic materials. We demonstrate that the light induces a helicity-dependent reduction of the magnitude of the magnetization. In the XUV regime, where the 3p semicore states are involved, a larger helicity dependence persisting even after the passage of light is exhibited. Finally, we were able to separate the part of the helicity-dependent dynamics due to the absorption from the part due to the inverse Faraday effect. Doing so, we show that the former has, overall, a greater impact on the magnetization than the latter, especially after the pulse and in the XUV regime. This work hints at the yet experimentally unexplored territory of the XUV light-induced helicity-dependent dynamics, which, according to our prediction, could magnify the helicity-dependent dynamics already exhibited in the optical regime.

Lifshitz Transition and Non-Fermi Liquid Behavior in Highly Doped Semimetals.

The classical Fermi liquid theory and Drude model have provided fundamental ways to understand the resistivity of most metals. The violation of the classical theory, known as non-Fermi liquid (NFL) transport, appears in certain metals, including topological semimetals, but quantitative understanding of the NFL behavior has not yet been established. In particular, the determination of the non-quadratic temperature exponent in the resistivity, a sign of NFL behavior, remains a puzzling issue. Here, a physical model to quantitatively explain the Lifshitz transition and NFL behavior in highly doped (a carrier density of ≈1022 cm-3 ) monoclinic Nb2 Se3 is reported. Hall and magnetoresistance measurements, the two-band Drude model, and first-principles calculations demonstrate an apparent chemical potential shift by temperature in monoclinic Nb2 Se3 , which induces a Lifshitz transition and NFL behavior in the material. Accordingly, the non-quadratic temperature exponent in the resistivity can be quantitatively determined by the chemical potential shift under the framework of Fermi liquid theory. This model provides a new experimental insight for nontrivial transport with NFL behavior or sign inversion of Seebeck coefficients in emerging materials.

Adsorption mechanisms of fatty acids on fluorite unraveled by infrared spectroscopy and first-principles calculations.

HYPOTHESIS: The adsorption mechanisms of fatty acids on minerals are largely debated from years, and their understanding is now required to improve flotation processing in the critical context of raw materials. Three wavenumbers have been observed in the literature for the asymmetric stretching vibration of COO- after the adsorption of fatty acids on mineral surfaces. They have been interpreted as different adsorbed forms, such as a precipitate formation, an adsorption of sole or bridged carboxylates, an anion exchange, or adsorbed modes, such as monodentate or bidentate configurations. EXPERIMENTS/THEORY: Diffuse reflectance infrared Fourier transform spectroscopy was combined with ab initio molecular dynamics simulations and simulation of infrared spectra. Fluorite and sodium octanoate - or longer-chain fatty acids - were used as prototypical materials for all the investigations. FINDINGS: At low fatty acids concentration, the asymmetric stretching vibration of COO- peaks at 1560 cm-1 while, at higher concentration, this infrared band converts into a doublet peaking at 1535 and 1575 cm-1. Using simulations, we assign the band at 1560 cm-1 to the adsorption of a carboxylate molecule bridged on a sodium counter-cation and the doublet at 1535 and 1575 cm-1 to the adsorption of the sole carboxylate anion under a monodentate or a bidentate binuclear configuration, respectively. The formation of an adsorbed layer on the mineral surface is initiated by the adsorption of a sodium carboxylate and followed by the adsorption of mixed sole anionic forms. The role of the carboxylate counter-cation is highlighted for the first time, which was totally ignored in the literature beforehand. This particularly opens the path to the development of innovative strategies to enhance the separation contrast between minerals, which is of uttermost importance for the recovery of critical raw materials.

Establishing the accuracy of density functional approaches for the description of noncovalent interactions in biomolecules.

Biomolecules have complex structures, and noncovalent interactions are crucial to determine their conformations and functionalities. It is therefore critical to be able to describe them in an accurate but efficient manner in these systems. In this context density functional theory (DFT) could provide a powerful tool to simulate biological matter either directly for relatively simple systems or coupled with classical simulations like the QM/MM (quantum mechanics/molecular mechanics) approach. Additionally, DFT could play a fundamental role to fit the parameters of classical force fields or to train machine learning potentials to perform large scale molecular dynamics simulations of biological systems. Yet, local or semi-local approximations used in DFT cannot describe van der Waals (vdW) interactions, one of the essential noncovalent interactions in biomolecules, since they lack a proper description of long range correlation effects. However, many efficient and reasonably accurate methods are now available for the description of van der Waals interactions within DFT. In this work, we establish the accuracy of several state-of-the-art vdW-aware functionals by considering 275 biomolecules including interacting DNA and RNA bases, peptides and biological inhibitors and compare our results for the energy with highly accurate wavefunction based calculations. Most methods considered here can achieve close to predictive accuracy. In particular, the non-local vdW-DF2 functional is revealed to be the best performer for biomolecules, while among the vdW-corrected DFT methods, uMBD is also recommended as a less accurate but faster alternative.

Ternary Chalcogenides BaMxTe2 (M = Cu, Ag): Syntheses, Modulated Crystal Structures, Optical Properties, and Electronic Calculations.

Standard solid-state methods produced black crystals of the compounds BaCu0.43(3)Te2 and BaAg0.77(1)Te2 at 1173 K; the crystal structures of each were established using single-crystal X-ray diffraction data. Both crystal structures are modulated. The compound BaCu0.43(3)Te2 crystallizes in the monoclinic superspace group P2(αß1/2)0, having cell dimensions of a = 4.6406(5) Å, b = 4.6596(5) Å, c = 10.362(1) Å, ß = 90.000(9)°, and Z = 2 and an incommensurate vector of q = 0.3499(6)b* + 0.5c*. The compound BaAg0.77(1)Te2 crystallizes in the orthorhombic P21212(α00)000 superspace group with cell dimensions of a = 4.6734(1) Å, b = 4.6468(1) Å, c = 11.1376(3) Å, and Z = 2 and an incommensurate vector of q = 0.364(2)a*. The asymmetric unit of the BaCu0.43(3)Te2 structure comprises eight crystallographically independent sites; that for BaAg0.77(1)Te2 comprises four. In these two structures, each of the M (M = Cu, Ag) atoms is connected to four Te atoms to make two-dimensional layers of [MxTe4/4]n- that are separated by layers of Ba atoms and square nets of Te. A Raman spectroscopic study at 298(2) K on a pelletized polycrystalline sample of BaAg0.8Te2 shows the presence of Ag-Te (83, 116, and 139 cm-1) and Ba-Te vibrations (667 and 732 cm-1). A UV-vis-NIR spectroscopic study on a powdered sample of BaAg0.8Te2 shows the semiconducting nature of the compound with a direct band gap of 1.0(2) eV, consistent with its black color. DFT calculations give a pseudo bandgap with a weak value of the DOS at the Fermi level.

Grafting of iron on amorphous silica surfaces from ab initio calculations.

Iron over silica catalytic systems have attracted considerable attention due to their activity and selectivity in different reactions, for instance, in the hydrodeoxygenation process. Here, the grafting mechanisms of iron under various forms (one atom, two atoms, or a cluster) on silica surfaces are studied using ab initio calculations. Various geometries with different locations of iron on the silica structure have been investigated, and it is found that a strong interaction between iron and the silanol groups takes place, mostly driven by the formation of Fe-O-Si bonds, and in few cases by nearby surface OH groups, creating Fe-OH-Si bonds. For the cluster, we show that the most favorable adsorption mode induces a strong effect on the silica surface accompanied with a large charge transfer, making it very stable and promising for a large panel of applications.

Modulated Linear Tellurium Chains in Ba3ScTe5: Synthesis, Crystal Structure, Optical and Resistivity Studies, and Electronic Structure.

A new ternary telluride, Ba3ScTe5, with a pseudo-one-dimensional structure, was synthesized at 1173 K by standard solid-state methods. A single-crystal X-ray diffraction study at 100(2) K shows the structure to be modulated. The structure of the subcell of Ba3ScTe5 crystallizes with two formula units in the hexagonal space group D6h3-P63/mcm with unit cell dimensions of a = b = 10.1190(5) Å and c = 6.8336(3) Å. The asymmetric unit of the subcell structure consists of four crystallographically independent sites: Ba1 (site symmetry: m2m), Sc1 (-3.m), Te1 (m2m), and Te2 (3.2). Its structure is made up of chains of ∞1[ScTe33-] that are separated by Ba2+ cations. The Sc atoms are bonded to six Te1 atoms that form a slightly distorted octahedral geometry. The structure of the subcell also contains linear infinite chains of Te2 with intermediate Te···Te interactions. The superstructure of Ba3ScTe5 is incommensurate and was solved in the hexagonal superspace group P-6(00γ)0 with a = 10.1188(3) Å and c = 6.8332(3) Å and a modulation vector of q = 0.3718(2)c*. The arrangement and coordination geometries of the atoms in the superstructure are very similar to those in the substructure. However, the main difference is that the infinite chains of Te atoms in the superstructure are distorted owing to the formation of long- and short-bonded pairs of Te atoms. The presence of these chains with intermediate Te···Te interactions makes assignment of the formal oxidation states arbitrary. The optical absorption study of a polycrystalline sample of Ba3ScTe5 that was synthesized by the stoichiometric reaction of elements at 1173 K reveals a direct band gap of 1.1(2) eV. The temperature-dependent resistivity study of polycrystalline Ba3ScTe5 shows semiconducting behavior corroborating the optical studies, while density functional theory calculations report a pseudo band gap of 1.3 eV.