Size-Dependent Ab Initio Atomistic Thermodynamics from Cluster to Bulk: Application to Hydration of Titania Nanoparticles.

Ab initio atomistic thermodynamics (AIAT) has become an indispensable tool to estimate Gibbs free energy changes for solid surfaces interacting with gaseous species relative to pressure (p) and temperature (T). For such systems, AIAT assumes that solid vibrational contributions to Gibbs free energy differences cancel out. However, the validity of this assumption is unclear for nanoscale systems. Using hydrated titania nanoparticles (NPs) as an example, we estimate the vibrational contributions to the Gibbs free energy of hydration (ΔGhyd(T,p)) for arbitrary NP size and degree of hydration. Comparing ΔGhyd(T,p) phase diagrams for NPs when considering these contributions (AIATnano) relative to a standard AIAT approach reveals significant qualitative and quantitative differences, which only become negligible for large systems. By constructing a size-dependent ΔGhyd(T,p) phase diagram, we illustrate how our approach can provide deeper insights into how nanosytems interact with their environments, with many potential applications (e.g., catalytic nanoparticles, biological colloids, nanoparticulate pollutants).

Gas-Phase Production of Hydroxylated Silicon Oxide Cluster Cations: Structure, Infrared Spectroscopy, and Astronomical Relevance.

The interaction of free cationic silicon oxide clusters, Si x O y + (x = 2-5, y ≥ x), with dilute water vapor, was investigated in a flow tube reactor. Product mass distributions indicate cluster size-dependent dissociative water adsorption. To probe the structure and vibrational spectra of the resulting Si x O y H2 + (x = 2-4) clusters, we employed infrared multiple photon dissociation spectroscopy and density functional theory calculations. The planar rhombic cluster core of the disilicon oxides (x = 2) appears to be retained upon dissociative adsorption of one H2O unit, whereas a significant structural transformation of the tri- and tetra-silicon oxides (x = 3 and 4) is induced, resulting in an increased coordination of the Si atoms and more 3D cluster structures. In an astronomical context, we discuss the potential relevance of Si x O y H z + clusters as seeds for dust nucleation and catalysts for carbon-based chemistry in diffuse or translucent interstellar clouds, where all the necessary conditions for producing these species are found. In the produced clusters, the frequency of the isolated silanol Si-OH stretching vibrational mode is considerably blue-shifted compared to that in hydroxylated bulk silica and small inorganic compounds. This mode has a characteristic frequency range between 1200 cm-1 (8.3 µm) and 1090 cm-1 (9.2 µm) and is associated with the anomalously small Si-OH bond lengths in these ionised species. In infrared observations such high frequency Si-O stretching modes are usually associated with a pure bulk silica component of silicate cosmic dust. The presence of Si x O y H2 + clusters in low silica astrophysical environments could thus potentially be detected via their signature Si-O band using the James Webb space telescope.

Tuning electronic levels in photoactive hydroxylated titania nanosystems: combining the ligand dipole effect and quantum confinement.

Reducing the size of titania (TiO2) to the nanoscale promotes the photoactive anatase phase for use in a range of applications from industrial catalysis to environment remediation. The nanoscale dimensions of these systems affect the magnitude of the electronic energy gap by quantum confinement. Upon interaction with aqueous environments or water vapour, the surfaces of these systems will also be hydroxylated to some degree. In turn, this affects the electronic energy levels due to the cumulative electrostatic effect of the dipolar hydroxyl (-OH) ligands (i.e. the ligand dipole effect). Using accurate density functional calculations, we investigate the combined effects of quantum confinement and the hydration-induced ligand dipole effect on a set of realistic titania nanosystems over a wide range of hydroxylation. Our detailed investigation reveals that, contrary to previous models, the ligand dipole effect does not-linearly depend on the ligand coverage due to the formation of inter-ligand OH⋯OH hydrogen bonds. To account for the resulting effects, we propose a refined model, which describes the ligand dipole effect more accurately in our systems. We show that both hydroxylation (by the ligand dipole effect) and size (by quantum confinement) have significant but distinct impacts on the electronic energy levels in nanotitania. As an example, we discuss how variations in these effects can be used to tune the highest unoccupied energy level in nanotitania for enhancing the efficiency of the hydrogen evolution reaction. Overall, we show that any specific energy shift can be achieved by a range of different combinations of nanosystem size and degree of hydroxylation, thus providing options for energy-level tuning while also allowing consideration of practical constraints (e.g. synthetic limitations, operating conditions) for photochemical applications.

Minimum conditions for accurate modeling of urea production via co-electrolysis.

Co-electrolysis of carbon oxides and nitrogen oxides promise to simultaneously help restore the balance of the C and N cycles while producing valuable chemicals such as urea. However, co-electrolysis processes are still largely inefficient and numerous knowledge voids persist. Here, we provide a solid thermodynamic basis for modelling urea production via co-electrolysis. First, we determine the energetics of aqueous urea produced under electrochemical conditions based on experimental data, which enables an accurate assessment of equilibrium potentials and overpotentials. Next, we use density functional theory (DFT) calculations to model various co-electrolysis reactions producing urea. The calculated reaction free energies deviate significantly from experimental values for well-known GGA, meta-GGA and hybrid functionals. These deviations stem from errors in the DFT-calculated energies of molecular reactants and products. In particular, the error for urea is approximately -0.25 ± 0.10 eV. Finally, we show that all these errors introduce large inconsistencies in the calculated free-energy diagrams of urea production via co-electrolysis, such that gas-phase corrections are strongly advised.

Predicting observable infrared signatures of nanosilicates in the diffuse interstellar medium.

The destruction time scale of dust in the diffuse interstellar medium is estimated to be an order of magnitude shorter than its residence time. Nevertheless, dust is observed in the interstellar medium, leading to the conclusion that reformation and grain growth must take place. Direct observations of nanometre-sized silicate grains, the main constituent of interstellar dust, would provide a smoking gun for the occurrence of grain condensation in the diffuse interstellar medium. Here we employ quantum chemical calculations to obtain the mid-infrared (IR) optical properties of a library of Mg-end member silicate nanoparticles with olivine (Mg2SiO4) and pyroxene (MgSiO3) stoichiometries. We use this library as an input for a foreground-screen model to predict the spectral appearance of the absorption profile due to mixtures of bulk and nanoparticle silicates towards bright background sources. The mid-IR spectrum observed towards an O8V star or a carbon-rich Wolf-Rayet star starts to change when ∼3% of the silicate mass is in the form of nanosilicates. We predict that a 3-10% nanosilicate fraction can be detected with the James Webb Space Telescope (JWST) using the mid-IR instrument (MIRI). With our upcoming JWST observations using MIRI, we will be able to detect or place limits on the nanosilicate content in the diffuse interstellar medium, and thus potentially directly confirm interstellar dust formation.

A machine learning potential for simulating infrared spectra of nanosilicate clusters.

The use of machine learning (ML) in chemical physics has enabled the construction of interatomic potentials having the accuracy of ab initio methods and a computational cost comparable to that of classical force fields. Training an ML model requires an efficient method for the generation of training data. Here, we apply an accurate and efficient protocol to collect training data for constructing a neural network-based ML interatomic potential for nanosilicate clusters. Initial training data are taken from normal modes and farthest point sampling. Later on, the set of training data is extended via an active learning strategy in which new data are identified by the disagreement between an ensemble of ML models. The whole process is further accelerated by parallel sampling over structures. We use the ML model to run molecular dynamics simulations of nanosilicate clusters with various sizes, from which infrared spectra with anharmonicity included can be extracted. Such spectroscopic data are needed for understanding the properties of silicate dust grains in the interstellar medium and in circumstellar environments.

Buckletronics: how compression-induced buckling affects the mechanical and electronic properties of sp2-based two-dimensional materials.

We explore the mechanical and electronic response of sp2-based two-dimensional materials under in-plane compression employing first principles density functional theory-based calculations. Taking two carbon-based graphynes (α-graphyne and γ-graphyne) as example systems, we show that the structures of both two-dimensional materials are susceptible to out-of-plane buckling, which emerges for modest in-plane biaxial compression (1.5-2%). Out-of-plane buckling is found to be more energetically stable than in-plane scaling/distortion and significantly lowers the in-plane stiffness of both graphenes. The buckling also gives rise to in-plane auxetic behaviour in both two-dimensional materials. Under compression, the induced in-plane distortions and out-of-plane buckling also lead to modulations of the electronic band gap. Our work highlights the possibility of using in-plane compression to induce out-of-plane buckling in, otherwise planar, sp2-based two-dimensional materials (e.g. graphynes, graphdiynes). We suggest that controllable compression-induced buckling in planar two-dimensional materials (as opposed to two-dimensional materials, which are buckled due to sp3 hybridization) could provide a route to a new 'buckletronics' approach for tuning the mechanical and electronic properties of sp2-based systems. This article is part of a discussion meeting issue 'Supercomputing simulations of advanced materials'.

Emergent Spin Frustration in Neutral Mixed-Valence 2D Conjugated Polymers: A Potential Quantum Materials Platform.

Two-dimensional conjugated polymers (2DCPs)─organic 2D materials composed of arrays of carbon sp2 centers connected by π-conjugated linkers─are attracting increasing attention due to their potential applications in device technologies. This interest stems from the ability of 2DCPs to host a range of correlated electronic and magnetic states (e.g., Mott insulators). Substitution of all carbon sp2 centers in 2DCPs by nitrogen or boron results in diamagnetic insulating states. Partial substitution of C sp2 centers by B or N atoms has not yet been considered for extended 2DCPs but has been extensively studied in the analogous neutral mixed-valence molecular systems. Here, we employ accurate first-principles calculations to predict the electronic and magnetic properties of a new class of hexagonally connected neutral mixed-valence 2DCPs in which every other C sp2 nodal center is substituted by either a N or B atom. We show that these neutral mixed-valence 2DCPs significantly energetically favor a state with emergent superexchange-mediated antiferromagnetic (AFM) interactions between C-based spin-1/2 centers on a triangular sublattice. These AFM interactions are surprisingly strong and comparable to those in the parent compounds of cuprate superconductors. The rigid and covalently linked symmetric triangular AFM lattice in these materials thus provides a highly promising and robust basis for 2D spin frustration. As such, extended mixed-valence 2DCPs are a highly attractive platform for the future bottom-up realization of a new class of all-organic quantum materials, which could host exotic correlated electronic states (e.g., unusual magnetic ordering, quantum spin liquids).

Crystal properties without crystallinity? Influence of surface hydroxylation on the structure and properties of small TiO2 nanoparticles.

Titania (TiO2) nanoparticles (NPs) are widely employed in applications that take advantage of their photochemical properties (e.g. pollutant degradation, photocatalysis). Here, we study the interrelation between crystallinity, surface hydroxylation and electronic structure in titania NPs with 1.4-2.3 nm diameters using all electron density functional theory-based calculations. We show how the distribution of local coordination environments of the atoms in thermally annealed quasi-spherical non-crystalline NPs converge to those in correspondingly sized faceted crystalline anatase NPs upon increasing hydroxylation. When highly hydroxylated, annealed NPs also possess electronic energy gaps with very similar energies and band edge orbital characters to those of the crystalline anatase NPs. We refer to the crystallite-mimicking non-crystalline annealed NPs as "crystalikes". Small stable crystalike NPs could allow for photochemical applications of titania in the size range where crystalline anatase NPs tend to become thermodynamically unfavoured (<3-5 nm). Our work implies the anatase crystal structure may not be as essential as previously assumed for TiO2 NP applications and generally suggests that crystalikes could be possible in other nanomaterials.

Stable Organic Radical for Enhancing Metal-Monolayer-Semiconductor Junction Performance.

The preparation of monolayers based on an organic radical and its diamagnetic counterpart has been pursued on hydrogen-terminated silicon surfaces. The functional monolayers have been investigated as solid-state metal/monolayer/semiconductor (MmS) junctions showing a characteristic diode behavior which is tuned by the electronic characteristics of the organic molecule. The eutectic gallium-indium liquid metal is used as a top electrode to perform the transport measurements and the results clearly indicate that the SOMO-SUMO molecular orbitals impact the device performance. The junction incorporating the radical shows an almost two orders of magnitude higher rectification ratio (R(|J1V/J-1V|) = 104.04) in comparison with the nonradical one (R(|J1V/J-1V|) = 102.30). The high stability of the fabricated MmS allows the system to be interrogated under irradiation, evidencing that at the wavelength where the photon energy is close to the band gap of the radical there is a clear enhancement of the photoresponse. This is translated into an increase of the photosensitivity (Sph) value from 68.7 to 269.0 mA/W for the nonradical and radical based systems, respectively.

How graphenic are graphynes? Evidence for low-lying correlated gapped states in graphynes.

Graphynes can be structurally envisioned as 2D extensions to graphene, whereby linearly bonded carbon linkages increase the distance between trigonal carbon nodes. Many graphynes have been predicted to exhibit a Dirac-like semimetallic (SEM) graphenic electronic structure, which could potentially make them competitive with graphene for applications. Currently, most graphynes remain as attractive synthetic targets, and their properties are still unconfirmed. Here, we demonstrate that the electronic structure of hexagonal α-graphyne is analogous to that of biaxially strained graphene. By comparison with accurate quantum Monte Carlo results on strained graphene, we show that the relative energetic stability of electronic states in this correlated 2D system can be captured by density functional theory (DFT) calculations using carefully tailored hybrid functionals. Our tuned hybrid DFT approach confirms that α-graphyne has a low energy correlated Mott-like antiferromagnetic insulating (AFI) state, which competes with the SEM state. Our work shows that the AFI-SEM crossover in α-graphyne could be tunable by in-plane biaxial strain. Applying our approach to other graphynes shows that they should also exhibit correlated AFI states, which could be dominant even at zero strain. Calculations using an onsite Coulombic repulsive term (i.e., DFT + U) also confirm the predictions of our hybrid DFT calculations. Overall, our work strongly suggests that graphynes are not as graphenic (i.e., Dirac-like) as often previously predicted by DFT calculations using standard generalized gradient approximation functionals. However, due to the greater electronic versatility (e.g., tunable semiconducting bandgaps and accessible spin polarized states) implied by our study, graphynes could have novel device applications that are complementary to those of graphene.

An unconstrained approach to systematic structural and energetic screening of materials interfaces.

From grain boundaries and heterojunctions to manipulating 2D materials, solid-solid interfaces play a key role in many technological applications. Understanding and predicting properties of these complex systems present an ongoing and increasingly important challenge. Over the last few decades computer simulation of interfaces has become vastly more powerful and sophisticated. However, theoretical interface screening remains based on largely heuristic methods and is strongly biased to systems that are amenable to modelling within constrained periodic cell approaches. Here we present an unconstrained and generally applicable non-periodic screening approach for systematic exploration of material's interfaces based on extracting and aligning disks from periodic reference slabs. Our disk interface method directly and accurately describes how interface structure and energetic stability depends on arbitrary relative displacements and twist angles of two interacting surfaces. The resultant detailed and comprehensive energetic stability maps provide a global perspective for understanding and designing interfaces. We confirm the power and utility of our method with respect to the catalytically important TiO2 anatase (101)/(001) and TiO2 anatase (101)/rutile (110) interfaces.

Cluster Beam Study of (MgSiO3)+-Based Monomeric Silicate Species and Their Interaction with Oxygen: Implications for Interstellar Astrochemistry.

Silicates are ubiquitously found as small dust grains throughout the universe. These particles are frequently subject to high-energy processes and subsequent condensation in the interstellar medium (ISM), where they are broken up into many ultrasmall silicate fragments. These abundant molecular-sized silicates likely play an important role in astrochemistry. By approximately mimicking silicate dust grain processing occurring in the diffuse ISM by ablation/cooling of a Mg/Si source material in the presence of O2, we observed the creation of stable clusters based on discrete pyroxene monomers (MgSiO3 +), which traditionally have only been considered possible as constituents of bulk silicate materials. Our study suggests that such pyroxene monomer-based clusters could be highly abundant in the ISM from the processing of larger silicate dust grains. A detailed analysis, by infrared multiple-photon dissociation (IR-MPD) spectroscopy and density functional theory (DFT) calculations, reveals the structures and properties of these monomeric silicate species. We find that the clusters interact strongly with oxygen, with some stable cluster isomers having a silicate monomeric core bound to an ozone-like moiety. The general high tendency of these monomeric silicate species to strongly adsorb O2 molecules also suggests that they could be relevant to the observed and unexplained depletion of oxygen in the ISM. We further find clusters where a Mg atom is bound to the MgSiO3 monomer core. These species can be considered as the simplest initial step in monomer-initiated nucleation, indicating that small ionized pyroxenic clusters could also assist in the reformation of larger silicate dust grains in the ISM.

Efficiency of Interstellar Nanodust Heating: Accurate Bottom-up Calculations of Nanosilicate Specific Heat Capacities.

Ultrasmall nanosized silicate grains are likely to be highly abundant in the interstellar medium. From sporadically absorbing energy from ultraviolet photons, these nanosilicates are subjected to significant instantaneous temperature fluctuations. These stochastically heated nanograins subsequently emit in the infrared. Previous estimates of the extent of the heating and emission have relied on empirical fits to bulk silicate heat capacities. The heat capacity of a system depends on the range of available vibrational modes, which for nanosized solids is dramatically affected by the constraints of finite size. Although experimental vibrational spectra of nanosilicates is not yet available, we directly take these finite size effects into account by using accurate vibrational spectra of low-energy nanosilicate structures from quantum chemical density functional theory calculations. Our results indicate that the heat capacities of ultrasmall nanosilicates are smaller than previously estimated, which would lead to a higher temperature and more intense infrared emission during stochastic heating. Specifically, we find that stochastically heated grains ultrasmall nanosilicates could be up to 35-80 K hotter than previously predicted. Our results could help to improve the understanding of infrared emission from ultrasmall nanosilicates in the ISM, which could be observed by the James Webb Space Telescope.

Can calculated harmonic vibrational spectra rationalize the structure of TiC-based nanoparticles?

Nanoscale titanium carbide (TiC) is widely used in composites and energy applications. In order to design and optimize these systems and to gain a fundamental understanding of these nanomaterials, it is important to understand the atomistic structure of nano-TiC. Cluster beam experiments have provided detailed infrared vibrational spectra of numerous TixCy nanoparticles with well defined masses. However, these spectra have yet to be convincingly assigned to TixCy nanoparticle structures. Herein, using accurate density functional theory based calculations, we perform a systematic survey of likely candidate nanoparticle structures with masses corresponding to those in experiment. We calculate harmonic infrared vibrational spectra for a range of nanoparticles up to 100 atoms in size, with a focus on systems based on removing either four carbon atoms or a single titanium atom from rocksalt-structured stoichiometric TiC nanoparticles. Our calculations clearly show that Ti-deficient nanoparticles are unlikely candidates to explain the experimental spectra as such structures are highly susceptible to C-C bonding, whose characteristic frequencies are not observed in experiment. However, our calculated infrared spectra for C-deficient nanoparticles have some matching features with the experimental spectra but tend to have more complex infrared spectra with more peaks than those obtained from experiment. We suggest that the discrepancy between experiment and theory may be largely due to thermally induced anharmonicities and broadening in the latter nanoparticles, which are not be accounted for in harmonic vibrational calculations.

Understanding the nature and location of hydroxyl groups on hydrated titania nanoparticles.

TiO2 nanoparticles (NPs) are intensively studied and widely used due to their huge potential in numerous applications involving their interaction with ultraviolet light (e.g., photocatalysis and sunscreens). Typically, these NPs are in water-containing environments and thus tend to be hydrated. As such, there is a growing need to better understand the physicochemical properties of hydrated TiO2 NPs in order to improve their performance in photochemical applications (e.g., photocatalytic water splitting) and to minimise their environmental impact (e.g., potential biotoxicity). To help address the need for reliable and detailed data on how nano-titania interacts with water, we present a systematic experimental and theoretical study of surface hydroxyl (OH) groups on photoactive anatase TiO2 NPs. Employing well-defined experimentally synthesised NPs and detailed realistic NP models, we obtain the measured and computed infrared spectra of the surface hydroxyls, respectively. By comparing the experimental and theoretical spectra we are able to identify the type and location of different OH groups in these NP systems. Specifically, our study allows us to provide unprecedented and detailed information about the coverage-dependent distribution of hydroxyl groups on the surface of experimental titania NPs, the degree of their H-bonding interactions and their associated assigned vibrational modes. Our work promises to lead to new routes for developing new and safe nanotechnologies based on hydrated TiO2 NPs.

Controlling pairing of π-conjugated electrons in 2D covalent organic radical frameworks via in-plane strain.

Controlling the electronic states of molecules is a fundamental challenge for future sub-nanoscale device technologies. π-conjugated bi-radicals are very attractive systems in this respect as they possess two energetically close, but optically and magnetically distinct, electronic states: the open-shell antiferromagnetic/paramagnetic and the closed-shell quinoidal diamagnetic states. While it has been shown that it is possible to statically induce one electronic ground state or the other by chemical design, the external dynamical control of these states in a rapid and reproducible manner still awaits experimental realization. Here, via quantum chemical calculations, we demonstrate that in-plane uniaxial strain of 2D covalently linked arrays of radical units leads to smooth and reversible conformational changes at the molecular scale that, in turn, induce robust transitions between the two kinds of electronic distributions. Our results pave a general route towards the external control, and thus technological exploitation, of molecular-scale electronic states in organic 2D materials.

Twistable dipolar aryl rings as electric field actuated conformational molecular switches.

The ability to control the chemical conformation of a system via external stimuli is a promising route for developing molecular switches. For eventual deployment as viable sub-nanoscale components that are compatible with current electronic device technology, conformational switching should be controllable by a local electric field (i.e. E-field gateable) and accompanied by a rapid and significant change in conductivity. In organic chemical systems the degree of π-conjugation is linked to the degree of electronic delocalisation, and thus largely determines the conductivity. Here, by means of accurate first principles calculations, we study the prototypical biphenyl based molecular system in which the dihedral angle between the two rings determines the degree of conjugation. In order to make this an E-field gateable system we create a net dipole by asymmetrically functionalising one ring with: (i) electron withdrawing (F, Br and CN), (ii) electron donating (NH2), and (iii) mixed (NH2/NO2) substituents. In this way, the application of an E-field interacts with the dipolar system to influence the dihedral angle, thus controlling the conjugation. For all considered substituents we consider a range of E-fields, and in each case extract conformational energy profiles. Using this data we obtain the minimum E-field required to induce a barrierless switching event for each system. We further extract the estimated switching speeds, the conformational probabilities at finite temperatures, and the effect of applied E-field on electronic structure. Consideration of these data allow us to assess which factors are most important in the design of efficient gateable electrical molecular switches.

Neutral Organic Radical Formation by Chemisorption on Metal Surfaces.

Organic radical monolayers (r-MLs) bonded to metal surfaces are potential materials for the development of molecular (spin)electronics. Typically, stable radicals bearing surface anchoring groups are used to generate r-MLs. Following a recent theoretical proposal based on a model system, we report the first experimental realization of a metal surface-induced r-ML, where a rationally chosen closed-shell precursor 3,5-dichloro-4-[bis(2,4,6-trichlorophenyl)methylen]cyclohexa-2,5-dien-1-one (1) transforms into a stable neutral open-shell species (1•) via chemisorption on the Ag(111) surface. X-ray photoelectron spectroscopy reveals that the >C═O group of 1 reacts with the surface, forming a C-O-Ag linkage that induces an electronic rearrangement that transforms 1 to 1•. We further show that surface reactivity is an important factor in this process whereby Au(111) is inert towards 1, whereas the Cu(111) surface leads to dehalogenation reactions. The radical nature of the Ag(111)-bound monolayer was further confirmed by angle-resolved photoelectron spectroscopy and electronic structure calculations, which provide evidence of the emergence of the singly occupied molecular orbital (SOMO) of 1•.