Triarylmethyl-based 2D covalent networks: virtual screening of chemical functionalisation for optimising strain-induced property control.

Two-dimensional covalent networks based on triarylmethyl (TAM) radical monomers have been proposed as versatile materials whose unpaired electrons may be externally localised/delocalised through the application of external uniaxial strain. This phenomenon arises through the strain-induced variance of the dihedral twist angles of the aryl rings within the network, and allows the control of important physico-chemical properties (e.g. magnetic interactions, electronic band gap). In order to experimentally realise such materials, one must find a compromise between the kinetic stability of the TAM monomers (through sterically protecting the radical centre with the appropriate aryl ring functionalisation) and the structural flexibility of the resulting material (provided by low intra-ring steric hindrance). In this work, through an efficient search procedure based on force field-based screening, employing ∼1750 calculations, followed by selected accurate electronic structure calculations, we provide support for the experimental viability of TAM-based 2D networks with highly controllable properties.

Exchange Coupling Inversion in a High-Spin Organic Triradical Molecule.

The magnetic properties of a nanoscale system are inextricably linked to its local environment. In adatoms on surfaces and inorganic layered structures, the exchange interactions result from the relative lattice positions, layer thicknesses, and other environmental parameters. Here, we report on a sample-dependent sign inversion of the magnetic exchange coupling between the three unpaired spins of an organic triradical molecule embedded in a three-terminal device. This ferro-to-antiferromagnetic transition is due to structural distortions and results in a high-to-low spin ground-state change in a molecule traditionally considered to be a robust high-spin quartet. Moreover, the flexibility of the molecule yields an in situ electric tunability of the exchange coupling via the gate electrode. These findings open a route to the controlled reversal of the magnetic states in organic molecule-based nanodevices by mechanical means, electrical gating, or chemical tailoring.

Challenges in modelling the reaction chemistry of interstellar dust.

Studies aiming to understand the physicochemical properties of interstellar dust and the chemical reactions that occur on and in it have traditionally been the preserve of astronomical observation and experimental attempts to mimic astronomically relevant conditions in the laboratory. Increasingly, computational modelling in its various guises is establishing a complementary third pillar of support to this endeavour by providing detailed insights into the complexities of interstellar dust chemistry. Inherently, the basis of computational modelling is to be found in the details (e.g. atomic structure/composition, reaction barriers) that are difficult to probe accurately from observation and experiment. This bottom-up atom-based theoretical approach, often itself based on deeper quantum mechanical principles, although extremely powerful, also has limitations when systems become too large or complex. In this Perspective, after first providing a general background to the current state of observational-based knowledge, we introduce a number of computational modelling methods with reference to recent state-of-the-art studies, in order to highlight the capabilities of such approaches in this field. Specifically, we first outline the use of computational chemistry methods for dust nucleation, structure, and individual reactions on bare and icy dust surfaces. Later, we review kinetic modelling of networks of reactions relevant to dust chemistry and how to take into account quantum tunnelling effects in the low temperature reactions in the interstellar medium. Finally, we point to the future challenges that need to be overcome for computational modelling to provide even more detailed and encompassing perspectives on the nature and reaction chemistry of interstellar dust.

One-dimensional embedded cluster approach to modeling CdS nanowires.

We present an embedded cluster model to treat one-dimensional nanostructures, using a hybrid quantum mechanical/molecular mechanical (QM/MM) approach. A segment of the nanowire (circa 50 atoms) is treated at a QM level of theory, using density functional theory (DFT) with a hybrid exchange-correlation functional. This segment is then embedded in a further length of wire, treated at an MM level of theory. The interaction between the QM and MM regions is provided by an embedding potential located at the interface. Point charges are placed beyond the ends of the wire segment in order to reproduce the Madelung potential of the infinite system. We test our model on the ideal system of a CdS linear chain, benchmarking our results against calculations performed on a periodic system using a plane-wave DFT approach, with electron exchange and correlation treated at the same level of approximation in both methods. We perform our tests on pure CdS and, importantly, the system containing a single In or Cu impurity. We find excellent agreement in the determined electronic structure using the two approaches, validating our embedded cluster model. As the hybrid QM/MM model avoids spurious interactions between charged defects, it will be of benefit to the analysis of the role of defects in nanowire materials, which is currently a major challenge using a plane-wave DFT approach. Other advantages of the hybrid QM/MM approach over plane-wave DFT include the ability to calculate ionization energies with an absolute reference and access to high levels of theory for the QM region which are not incorporated in most plane-wave codes. Our results concur with available experimental data.

Long range coupling between defect centres in inorganic nanostructures: valence alternation pairs in nanoscale silica.

Valence alternation pair (VAP) states are formed by a closed-shell combination of two space- and charge-separated topological defect centres. These pairs of defects, although historically invoked to explain the electronic properties of bulk inorganic glassy materials (e.g., amorphous silicon dioxide) via the concept of negative-U defects, have more recently been found in a number of theoretical studies of silica surfaces and nanoscale silica clusters. Using density functional theory we systematically probe the structure and internal stability of VAPs in a number of silica nanoclusters with respect to the separation of the two constituent defect centres. We find that VAP states in nanosilica are strongly stabilised by the attractive electrostatic interaction between their separated oppositely charged component defects such that VAPs can persist up to an internal separation of a least 1.5 nanometres. Beyond this distance VAPs become unstable with respect to an open-shell combination of topological defects, virtually indistinguishable from two isolated open-shell defect centres. Finally, we theoretically analyse the possibility of experimental observation of VAP states through their infra-red vibrational spectra.

The effect of local environment on photoluminescence: a time-dependent density functional theory study of silanone groups on the surface of silica nanostructures.

The optical absorption spectrum and lowest photoluminescence (PL) signal for silanone terminated silica nanostructures are studied using time-dependent density functional theory calculations on a range of realistic low energy silica nanocluster models. We show that the broad experimental absorption spectrum for silanone centers [V. A. Radtsig and I. M. Senchenya Russ. Chem. Bull. 45, 1849 (1996)] is most likely the result of a synergetic combination of inhomogeneous broadening, thermal broadening and the small energy differences between different excitations. We further demonstrate that upon relaxation of the excited state the excited electron and hole localize on only one silanone center, and that there is a clear and distinct link between the local environment of a silanone center and its absorption and PL spectra. Finally, we provide strong evidence that the silanone center does not have a double bond between the constituent silicon and oxygen atoms but rather can be probably more aptly described as the =Si(+)-O(-) charge-transfer species.

Optical excitations of defects in realistic nanoscale silica clusters: comparing the performance of density functional theory using hybrid functionals with correlated wavefunction methods.

Optical excitations of low energy silica (SiO(2))(4) clusters obtained by global optimization, as opposed to constructed by hand, are studied using a range of theoretical methods. By focusing on the lowest energy silica clusters we hope to capture at least some of the characteristic ways by which the dry surfaces of silica nanosystems preferentially terminate. Employing the six lowest energy (SiO(2))(4) cluster isomers, we show that they exhibit a surprisingly wide range of geometries, defects, and associated optical excitations. Some of the clusters show excitations localized on isolated defects, which are known from previous studies using hydrogen-terminated versions of the defect in question. Other clusters, however, exhibit novel charge-transfer excitations in which an electron transfers between two spatially separated defects. In these cases, because of the inherent proximity of the constituent defects due to the small cluster dimensions, the excitation spectrum is found to be very different from that of the same defects in isolation. Excitation spectra of all clusters were calculated using time-dependent density functional theory (TD-DFT) and delta-SCF DFT (DeltaDFT) methods employing two different hybrid density functionals (B3LYP and BB1K) differing essentially in the amount of incorporated Hartree-Fock-like exchange (HFLE). In all cases the results were compared with CASPT2 calculated values which are taken as a benchmark standard. In line with previous work, the spatially localized excitations are found to be well described by TD-DFT/B3LYP but which gives excitation energies that are significantly underestimated in the case of the charge-transfer excitations. The TD-DFT/BB1K combination in contrast is found to give generally good excitation energies for the lowest excited states of both localized and charge-transfer excitations. Finally, our calculations suggest that the increased quality of the predicted excitation spectra by adding larger amounts of HFLE is mainly due to an increased localization of the excited state associated with the elimination of spurious self-interaction inherent to (semi-)local DFT functionals.

Self-diffusion of molecular hydrogen in clathrasils compared: dodecasil 3C versus sodalite.

The self-diffusion coefficient of molecular hydrogen through the all-silica microporous dodecasil 3C structure is calculated by means of molecular-dynamics (MD) calculations, allowing for full framework flexibility, in order to assess the material's feasibility as a hydrogen storage medium. The hydrogen uptake rate into dodecasil 3C is compared to that previously calculated for sodalite and it is found that the latter performs significantly better. The reason for this variation in performance is found to lie in intrinsic topological differences between each framework type. This is explicitly demonstrated by means of a simplified version of transition state theory helping to succinctly rationalize the MD data.

Novel structures and energy spectra of hydroxylated (SiO2)8-based clusters: searching for the magic (SiO2)8O2H3- cluster.

The prominent (SiO(2))(8)O(2)H(3) (-) mass peak resulting from the laser ablation of hydroxylated silica, attributed to magic cluster formation, is investigated employing global optimization with a dedicated interatomic potential and density functional calculations. The low-energy spectra of cluster isomers are calculated for the closed shell clusters: (SiO(2))(8)OH(-) and (SiO(2))(8)O(2)H(3) (-) giving the likely global minima in each case. Based upon our calculated cluster structures and energetics, and further on the known experimental details, it is proposed that the abundant formation of (SiO(2))(8)O(2)H(3) (-) clusters is largely dependent on the high stability of the (SiO(2))(8)OH(-) ground state cluster. Both the (SiO(2))(8)O(2)H(3) (-) and (SiO(2))(8)OH(-) ground state clusters are found to exhibit cagelike structures with the latter containing a particularly unusual tetrahedrally four-coordinated oxygen center not observed before in either bulk silica or silica clusters. The bare ground state (SiO(2))(8)O(2-) cluster ion core is also found to have four tetrahedrally symmetric Si==O terminations making it a possible candidate, when combined with suitable cations, for extended cluster-based structures/materials.

Effect of cation distribution on self-diffusion of molecular hydrogen in Na3Al3Si3O12 sodalite: a molecular dynamics study.

The diffusion of hydrogen in sodium aluminum sodalite (NaAlSi-SOD) is modeled using classical molecular dynamics, allowing for full flexibility of the host framework, in the temperature range 800-1200 K. From these simulations, the self-diffusion coefficient is determined as a function of temperature and the hydrogen uptake at low equilibrium hydrogen concentration is estimated at 573 K. The influence of the cation distribution over the framework on the hydrogen self-diffusion is investigated by comparing results employing a low energy fully ordered cation distribution with those obtained using a less ordered distribution. The cation distribution is found to have a surprisingly large influence on the diffusion, which appears to be due to the difference in framework flexibility for different cation distributions, the occurrence of correlated hopping in case of the ordered distribution, and the different nature of the diffusion processes in both systems. Compared to our previously reported calculations on all silica sodalite (all-Si-SOD), the hydrogen diffusion coefficient of sodium aluminum sodalite is higher in the case of the ordered distribution and lower in case of the disordered distribution. The hydrogen uptake rates of all-Si-SOD and NaSiAl-SOD are comparable at high temperatures (approximately 1000 K) and lower for all-Si-SOD at lower temperatures (approximately 400 K).

Molecular-dynamics analysis of the diffusion of molecular hydrogen in all-silica sodalite.

In order to investigate the technical feasibility of crystalline porous silicates as hydrogen storage materials, the self-diffusion of molecular hydrogen in all-silica sodalite is modeled using large-scale classical molecular-dynamics simulations employing full lattice flexibility. In the temperature range of 700-1200 K, the diffusion coefficient is found to range from 1.610(-10) to 1.810(-9) m(2)/s. The energy barrier for hydrogen diffusion is determined from the simulations allowing the application of transition state theory, which, together with the finding that the pre-exponential factor in the Arrhenius-type equation for the hopping rate is temperature-independent, enables extrapolation of our results to lower temperatures. Estimates based on mass penetration theory calculations indicate a promising hydrogen uptake rate at 573 K.

Fully coordinated silica nanoclusters: (SiO2)N molecular rings.

A new form of finite silica with edge-sharing SiO2 units connected in a ring is proposed. High-level density-functional calculations for (SiO2)(N), N=4-14, show the rings to be energetically more stable than the corresponding (SiO2)(N) linear chains for N>11. The rings display frequency modes in remarkable agreement with infrared bands measured on dehydrated silica surfaces indicating their potential as models of strained extended silica systems. Silica rings, if synthesized, may also be useful precursors for new bulk-silica polymorphs with tubular or porous morphologies.