Giant Carrier Mobility in a Room-Temperature Ferromagnetic VSi2N4 Monolayer.

Using density functional theory (DFT), we investigate that two possible phases of VSi2N4 (VSN) may be realized, one called the "H phase" corresponding to what is known from calculation and herein the other new "T phase" being stabilized by a biaxial tensile strain of 3%. Significantly, the H phase is predicted to display a giant carrier mobility of 1 × 106 cm2 V-1 s-1, which exceeds that for most 2D magnetic materials, with a Curie temperature (TC) exceeding room temperature and a band gap of 2.01 eV at the K point. Following the H-T phase transition, the direct band gap shifts to the Γ point and increases to 2.59 eV. The Monte Carlo (MC) simulations also indicate that TC of the T phase VSN can be effectively modulated by strain, reaching room temperature under a biaxial strain of -4%. These results show that VSN should be a promising functional material for future nanoelectronics.

Controlling piezoresistance in single molecules through the isomerisation of bullvalenes.

Nanoscale electro-mechanical systems (NEMS) displaying piezoresistance offer unique measurement opportunities at the sub-cellular level, in detectors and sensors, and in emerging generations of integrated electronic devices. Here, we show a single-molecule NEMS piezoresistor that operates utilising constitutional and conformational isomerisation of individual diaryl-bullvalene molecules and can be switched at 850 Hz. Observations are made using scanning tunnelling microscopy break junction (STMBJ) techniques to characterise piezoresistance, combined with blinking (current-time) experiments that follow single-molecule reactions in real time. A kinetic Monte Carlo methodology (KMC) is developed to simulate isomerisation on the experimental timescale, parameterised using density-functional theory (DFT) combined with non-equilibrium Green's function (NEGF) calculations. Results indicate that piezoresistance is controlled by both constitutional and conformational isomerisation, occurring at rates that are either fast (equilibrium) or slow (non-equilibrium) compared to the experimental timescale. Two different types of STMBJ traces are observed, one typical of traditional experiments that are interpreted in terms of intramolecular isomerisation occurring on stable tipped-shaped metal-contact junctions, and another attributed to arise from junction‒interface restructuring induced by bullvalene isomerisation.

Simulating optical linear absorption for mesoscale molecular aggregates: An adaptive hierarchy of pure states approach.

In this paper, we present dyadic adaptive HOPS (DadHOPS), a new method for calculating linear absorption spectra for large molecular aggregates. This method combines the adaptive HOPS (adHOPS) framework, which uses locality to improve computational scaling, with the dyadic HOPS method previously developed to calculate linear and nonlinear spectroscopic signals. To construct a local representation of dyadic HOPS, we introduce an initial state decomposition that reconstructs the linear absorption spectra from a sum over locally excited initial conditions. We demonstrate the sum over initial conditions can be efficiently Monte Carlo sampled and that the corresponding calculations achieve size-invariant [i.e., O(1)] scaling for sufficiently large aggregates while trivially incorporating static disorder in the Hamiltonian. We present calculations on the photosystem I core complex to explore the behavior of the initial state decomposition in complex molecular aggregates as well as proof-of-concept DadHOPS calculations on an artificial molecular aggregate inspired by perylene bis-imide to demonstrate the size-invariance of the method.

Density functionals with asymptotic-potential corrections are required for the simulation of spectroscopic properties of materials.

Five effects of correction of the asymptotic potential error in density functionals are identified that significantly improve calculated properties of molecular excited states involving charge-transfer character. Newly developed materials-science computational methods are used to demonstrate how these effects manifest in materials spectroscopy. Connection is made considering chlorophyll-a as a paradigm for molecular spectroscopy, 22 iconic materials as paradigms for 3D materials spectroscopy, and the VN - defect in hexagonal boron nitride as an example of the spectroscopy of defects in 2D materials pertaining to nanophotonics. Defects can equally be thought of as being "molecular" and "materials" in nature and hence bridge the relms of molecular and materials spectroscopies. It is concluded that the density functional HSE06, currently considered as the standard for accurate calculations of materials spectroscopy, should be replaced, in most instances, by the computationally similar but asymptotically corrected CAM-B3LYP functional, with some specific functionals for materials-use only providing further improvements.

Accurate prediction of the properties of materials using the CAM-B3LYP density functional.

Density functionals with asymptotic corrections to the long-range potential provide entry-level methods for calculations on molecules that can sustain charge transfer, but similar applications in materials science are rare. We describe an implementation of the CAM-B3LYP range-separated functional within the Vienna Ab-initio Simulation Package (VASP) framework, together with its analytical functional derivatives. Results obtained for eight representative materials: aluminum, diamond, graphene, silicon, NaCl, MgO, 2D h-BN, and 3D h-BN, indicate that CAM-B3LYP predictions embody mean-absolute deviations (MAD) compared to HSE06 that are reduced by a factor of six for lattice parameters, four for quasiparticle band gaps, three for the lowest optical excitation energies, and six for exciton binding energies. Further, CAM-B3LYP appears competitive compared to ab initio G0 W0 and Bethe-Salpeter equation approaches. The CAM-B3LYP implementation in VASP was verified by comparison of optimized geometries and reaction energies for isolated molecules taken from the ACCDB database, evaluated in large periodic unit cells, to analogous results obtained using Gaussian basis sets. Using standard GW pseudopotentials and energy cutoffs for the plane-wave calculations and the aug-cc-pV5Z basis set for the atomic-basis ones, the MAD in energy for 1738 chemical reactions was 0.34 kcal mol-1 , while for 480 unique bond lengths this was 0.0036 Å; these values reduced to 0.28 kcal mol-1 (largest error 0.94 kcal mol-1 ) and 0.0009 Å by increasing the plane-wave cutoff energy to 850 eV.

Identifying carbon as the source of visible single-photon emission from hexagonal boron nitride.

Single-photon emitters (SPEs) in hexagonal boron nitride (hBN) have garnered increasing attention over the last few years due to their superior optical properties. However, despite the vast range of experimental results and theoretical calculations, the defect structure responsible for the observed emission has remained elusive. Here, by controlling the incorporation of impurities into hBN via various bottom-up synthesis methods and directly through ion implantation, we provide direct evidence that the visible SPEs are carbon related. Room-temperature optically detected magnetic resonance is demonstrated on ensembles of these defects. We perform ion-implantation experiments and confirm that only carbon implantation creates SPEs in the visible spectral range. Computational analysis of the simplest 12 carbon-containing defect species suggest the negatively charged [Formula: see text] defect as a viable candidate and predict that out-of-plane deformations make the defect environmentally sensitive. Our results resolve a long-standing debate about the origin of single emitters at the visible range in hBN and will be key to the deterministic engineering of these defects for quantum photonic devices.

Silicon - single molecule - silicon circuits.

In 2020, silicon - molecule - silicon junctions were fabricated and shown to be on average one third as conductive as traditional junctions made using gold electrodes, but in some instances to be even more conductive, and significantly 3 times more extendable and 5 times more mechanically stable. Herein, calculations are performed of single-molecule junction structure and conductivity pertaining to blinking and scanning-tunnelling-microscopy (STM) break junction (STMBJ) experiments performed using chemisorbed 1,6-hexanedithiol linkers. Some strikingly different characteristics are found compared to analogous junctions formed using the metals which, to date, have dominated the field of molecular electronics. In the STMBJ experiment, following retraction of the STM tip after collision with the substrate, unterminated silicon surface dangling bonds are predicted to remain after reaction of the fresh tips with the dithiol solute. These dangling bonds occupy the silicon band gap and are predicted to facilitate extraordinary single-molecule conductivity. Enhanced junction extendibility is attributed to junction flexibility and the translation of adsorbed molecules between silicon dangling bonds. The calculations investigate a range of junction atomic-structural models using density-functional-theory (DFT) calculations of structure, often explored at 300 K using molecular dynamics (MD) simulations. These are aided by DFT calculations of barriers for passivation reactions of the dangling bonds. Thermally averaged conductivities are then evaluated using non-equilibrium Green's function (NEGF) methods. Countless applications through electronics, nanotechnology, photonics, and sensing are envisaged for this technology.

Asymmetry in the Q y Fluorescence and Absorption Spectra of Chlorophyll a Pertaining to Exciton Dynamics.

Significant asymmetry found between the high-resolution Q y emission and absorption spectra of chlorophyll-a is herein explained, providing basic information needed to understand photosynthetic exciton transport and photochemical reactions. The Q y spectral asymmetry in chlorophyll has previously been masked by interference in absorption from the nearby Q x transition, but this effect has recently been removed using extensive quantum spectral simulations or else by analytical inversion of absorption and magnetic circular dichroism data, allowing high-resolution absorption information to be accurately determined from fluorescence-excitation spectra. To compliment this, here, we measure and thoroughly analyze the high-resolution differential fluorescence line narrowing spectra of chlorophyll-a in trimethylamine and in 1-propanol. The results show that vibrational frequencies often change little between absorption and emission, yet large changes in line intensities are found, this effect also being strongly solvent dependent. Among other effects, the analysis in terms of four basic patterns of Duschinsky-rotation matrix elements, obtained using CAM-B3LYP calculations, predicts that a chlorophyll-a molecule excited into a specific vibrational level, may, without phase loss or energy relaxation, reemit the light over a spectral bandwidth exceeding 1,000 cm-1 (0.13 eV) to influence exciton-transport dynamics.

Covalent Linkages of Molecules and Proteins to Si-H Surfaces Formed by Disulfide Reduction.

Thiols and disulfide contacts have been, for decades, key for connecting organic molecules to surfaces and nanoclusters as they form self-assembled monolayers (SAMs) on metals such as gold (Au) under mild conditions. In contrast, they have not been similarly deployed on Si owing to the harsh conditions required for monolayer formation. Here, we show that SAMs can be simply formed by dipping Si-H surfaces into dilute solutions of organic molecules or proteins comprising disulfide bonds. We demonstrate that S-S bonds can be spontaneously reduced on Si-H, forming covalent Si-S bonds in the presence of traces of water, and that this grafting can be catalyzed by electrochemical potential. Cyclic disulfide can be spontaneously reduced to form complete monolayers in 1 h, and the reduction can be catalyzed electrochemically to form full surface coverages within 15 min. In contrast, the kinetics of SAM formation of the cyclic disulfide molecule on Au was found to be three-fold slower than that on Si. It is also demonstrated that dilute thiol solutions can form monolayers on Si-H following oxidation to disulfides under ambient conditions; the supply of too much oxygen, however, inhibits SAM formation. The electron transfer kinetics of the Si-S-enabled SAMs on Si-H is comparable to that on Au, suggesting that Si-S contacts are electrically transmissive. We further demonstrate the prospect of this spontaneous disulfide reduction by forming a monolayer of protein azurin on a Si-H surface within 1 h. The direct reduction of disulfides on Si electrodes presents new capabilities for a range of fields, including molecular electronics, for which highly conducting SAM-electrode contacts are necessary and for emerging fields such as biomolecular electronics as disulfide linkages could be exploited to wire proteins between Si electrodes, within the context of the current Si-based technologies.

Source code, input data, and sample output concerning the application of multistate density functional theory to the singdoublet and tripdoublet states of the ethylene cation.

This is the data and associated new software required to run multi-state density-functional theory (MSDDFT) calculations by the GAMESS programme. Also, data and software needed to drive GAMESS based on output from the Gaussian-16 package is included. Sample input and output files are included, as well as Perl scripts and Fortran source code. A separate execution of the scripts is required to create the input specifications for each state to be included in the MSDFT, then after GAMESS is run more software is included to calculate the final state energies. The associated basic theory and results are described in "Multistate density functional theory applied with 3 unpaired electrons in 3 orbitals: the singdoublet and tripdoublet states of the ethylene cation" [1].

Spontaneous S-Si bonding of alkanethiols to Si(111)-H: towards Si-molecule-Si circuits.

We report the synthesis of covalently linked self-assembled monolayers (SAMs) on silicon surfaces, using mild conditions, in a way that is compatible with silicon-electronics fabrication technologies. In molecular electronics, SAMs of functional molecules tethered to gold via sulfur linkages dominate, but these devices are not robust in design and not amenable to scalable manufacture. Whereas covalent bonding to silicon has long been recognized as an attractive alternative, only formation processes involving high temperature and/or pressure, strong chemicals, or irradiation are known. To make molecular devices on silicon under mild conditions with properties reminiscent of Au-S ones, we exploit the susceptibility of thiols to oxidation by dissolved O2, initiating free-radical polymerization mechanisms without causing oxidative damage to the surface. Without thiols present, dissolved O2 would normally oxidize the silicon and hence reaction conditions such as these have been strenuously avoided in the past. The surface coverage on Si(111)-H is measured to be very high, 75% of a full monolayer, with density-functional theory calculations used to profile spontaneous reaction mechanisms. The impact of the Si-S chemistry in single-molecule electronics is demonstrated using STM-junction approaches by forming Si-hexanedithiol-Si junctions. Si-S contacts result in single-molecule wires that are mechanically stable, with an average lifetime at room temperature of 2.7 s, which is five folds higher than that reported for conventional molecular junctions formed between gold electrodes. The enhanced "ON" lifetime of this single-molecule circuit enables previously inaccessible electrical measurements on single molecules.

Single-photon emitters in hexagonal boron nitride: a review of progress.

This report summarizes progress made in understanding properties such as zero-phonon-line energies, emission and absorption polarizations, electron-phonon couplings, strain tuning and hyperfine coupling of single photon emitters in hexagonal boron nitride. The primary aims of this research are to discover the chemical nature of the emitting centres and to facilitate deployment in device applications. Critical analyses of the experimental literature and data interpretation, as well as theoretical approaches used to predict properties, are made. In particular, computational and theoretical limitations and challenges are discussed, with a range of suggestions made to overcome these limitations, striving to achieve realistic predictions concerning the nature of emitting centers. A symbiotic relationship is required in which calculations focus on properties that can easily be measured, whilst experiments deliver results in a form facilitating mass-produced calculations.

Absorption-emission symmetry breaking and the different origins of vibrational structures of the 1Qy and 1Qx electronic transitions of pheophytin a.

The vibrational structure of the optical absorption and fluorescence spectra of the two lowest-energy singlet electronic states (Qy and Qx) of pheophytin a were carefully studied by combining low-resolution and high-resolution spectroscopy with quantum chemical analysis and spectral modeling. Large asymmetry was revealed between the vibrational structures of the Qy absorption and fluorescence spectra, integrally characterized by the total Huang-Rhys factor and reorganization energy in absorption of Svib A = 0.43 ± 0.06, λA = 395 cm-1 and in emission of Svib E = 0.35 ± 0.06, λE = 317 cm-1. Time-dependent density-functional theory using the CAM-B3LYP, ωB97XD, and MN15 functionals could predict and interpret this asymmetry, with the exception of one vibrational mode per model, which was badly misrepresented in predicted absorption spectra; for CAM-B3LYP and ωB97XD, this mode was a Kekulé-type mode depicting aromaticity. Other computational methods were also considered but performed very poorly. The Qx absorption spectrum is broad and could not be interpreted in terms of a single set of Huang-Rhys factors depicting Franck-Condon allowed absorption, with Herzberg-Teller contributions to the intensity being critical. For it, CAM-B3LYP calculations predict that Svib A (for modes >100 cm-1) = 0.87 and λA = 780 cm-1, with effective x and y polarized Herzberg-Teller reorganization energies of 460 cm-1 and 210 cm-1, respectively, delivering 15% y-polarized intensity. However, no method was found to quantitatively determine the observed y-polarized contribution, with contributions of up to 50% being feasible.

Competition between charge migration and charge transfer induced by nuclear motion following core ionization: Model systems and application to Li2.

Attosecond and femtosecond spectroscopies present opportunities for the control of chemical reaction dynamics and products, as well as for quantum information processing; we address the somewhat unique situation of core-ionization spectroscopy which, for dimeric chromophores, leads to strong valence charge localization and hence tightly paired potential-energy surfaces of very similar shape. Application is made to the quantum dynamics of core-ionized Li2 +. This system is chosen as Li2 is the simplest stable molecule facilitating both core ionization and valence ionization. First, the quantum dynamics of some model surfaces are considered, with the surprising result that subtle differences in shape between core-ionization paired surfaces can lead to dramatic differences in the interplay between electronic charge migration and charge transfer induced by nuclear motion. Then, equation-of-motion coupled-cluster calculations are applied to determine potential-energy surfaces for 8 core-excited state pairs, calculations believed to be the first of their type for other than the lowest-energy core-ionized molecular pair. While known results for the lowest-energy pair suggest that Li2 + is unsuitable for studying charge migration, higher-energy pairs are predicted to yield results showing competition between charge migration and charge transfer. Central is a focus on the application of Hush's 1975 theory for core-ionized X-ray photoelectron spectroscopy to understand the shapes of the potential-energy surfaces and hence predict key features of charge migration.

Bioferroelectric Properties of Glycine Crystals.

Biological ferroelectric materials have great potential in biosensing and disease diagnosis and treatment. Glycine crystals form the simplest bioferroelectric materials, and here we investigate the polarizations of its ß- and γ-phases. Using density functional theory, we predict that glycine crystals can develop polarizations  even larger than those of conventional inorganic ferroelectrics. Further, using systematic molecular dynamics simulations utilizing polarized crystal charges, we predict the Curie temperature of γ-glycine to be 630 K, with a required coercive field to switch its polarization states of 1 V·nm-1, consistent with experimental evidence. This work sheds light on the microscopic mechanism of electric dipole ordering in biomaterials, helping in the material design of novel bioferroelectrics.

van der Waals forces control ferroelectric-antiferroelectric ordering in CuInP2S6 and CuBiP2Se6 laminar materials.

We show how van der Waals (vdW) forces outcompete covalent and ionic forces to control ferroelectric ordering in CuInP2S6 nanoflakes as well as in CuInP2S6 and CuBiP2Se6 crystals. While the self-assembly of these 2D layered materials is clearly controlled by vdW effects, this result indicates that the internal layer structure is also similarly controlled. Using up to 14 first-principles computational methods, we predict that the bilayers of both materials should be antiferroelectric. However, antiferroelectric nanoflakes and bulk materials are shown to embody two fundamentally different types of inter-layer interactions, with vdW forces strongly favouring one and strongly disfavouring the other compared to ferroelectric ordering. Strong specific vdW interactions involving the Cu atoms control this effect. Thickness-dependent significant cancellation of these two large opposing vdW contributions results in a small net effect that interacts with weak ionic contributions to control ferroelectric ordering.

Cyanine platelet single crystals: growth, crystal structure and optical spectra.

Crystalline organic semiconducting materials are much in demand for multiple electronic and optoelectronic device applications. Here, solution grown ultrathin rhombic crystals of a trimethine carbocyanine anionic dye are used to establish relationships between structural and optical properties. The dye crystallized in the monoclinic space group P21/c featuring alternating layers of molecules in two different herringbone type patterns, with perchlorate counterions located mostly within one of the two layers. Micro transmittance spectroscopy revealed a broadened spectrum compared to those obtained in solution and in an amorphous thin film. Using polarized light, transmission spectroscopy revealed strong low-energy and weak high-energy bands polarized along the crystallographic b- and c-axis, respectively. Using the extended dipole approximation, significant exciton couplings are predicted between neighboring molecules in the crystal, of the order of the intrinsic monomer reorganization energies associated with nuclear relaxation after excitation, depicting a complex spectral scenario. The exciton coupling pattern explains the relative energies of the b- and c-polarized components but the observed intensities are opposite to expectations based on chromophore alignment within the crystal.

Faraday cage screening reveals intrinsic aspects of the van der Waals attraction.

General properties of the recently observed screening of the van der Waals (vdW) attraction between a silica substrate and silica tip by insertion of graphene are predicted using basic theory and first-principles calculations. Results are then focused on possible practical applications, as well as an understanding of the nature of vdW attraction, considering recent discoveries showing it competing against covalent and ionic bonding. The traditional view of the vdW attraction as arising from pairwise-additive London dispersion forces is considered using Grimme's "D3" method, comparing results to those from Tkatchenko's more general many-body dispersion (MBD) approach, all interpreted in terms of Dobson's general dispersion framework. Encompassing the experimental results, MBD screening of the vdW force between two silica bilayers is shown to scale up to medium separations as 1.25 de/d, where d is the bilayer separation and de is its equilibrium value, depicting antiscreening approaching and inside de Means of unifying this correlation effect with those included in modern density functionals are urgently required.