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.

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.

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.

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.

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.

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.

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.

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.

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.

Gold surfaces and nanoparticles are protected by Au(0)-thiyl species and are destroyed when Au(I)-thiolates form.

The synthetic chemistry and spectroscopy of sulfur-protected gold surfaces and nanoparticles is analyzed, indicating that the electronic structure of the interface is Au(0)-thiyl, with Au(I)-thiolates identified as high-energy excited surface states. Density-functional theory indicates that it is the noble character of gold and nanoparticle surfaces that destabilizes Au(I)-thiolates. Bonding results from large van der Waals forces, influenced by covalent bonding induced through s-d hybridization and charge polarization effects that perturbatively mix in some Au(I)-thiolate character. A simple method for quantifying these contributions is presented, revealing that a driving force for nanoparticle growth is nobleization, minimizing Au(I)-thiolate involvement. Predictions that Brust-Schiffrin reactions involve thiolate anion intermediates are verified spectroscopically, establishing a key feature needed to understand nanoparticle growth. Mixing of preprepared Au(I) and thiolate reactants always produces Au(I)-thiolate thin films or compounds rather than monolayers. Smooth links to O, Se, Te, C, and N linker chemistry are established.

How Equilibrium Gets Mimicked During Kinetic and Thermodynamic Control in Porphyrin and Phthalocyanine Self-Assembled Monolayers.

The recent review of Hipps and Mazur is extended to consider the origins and significance of their conclusion that "surface structures vary with relative component concentration in a way that may mimic equilibria but is not". How this situation can arise during self-assembly is discussed, as well as a range of qualitative and quantitative observations and high-level free-energy calculations that document the effect for meso-tetraalkylporphyrins self-assembled monolayer (SAM) polymorphs. This leads to a discussion of modern challenges facing quantification of the effects caused by kinetic control, as well as to the usefulness of equilibrium mimicking in the design and synthesis of SAMs.

Molecular quantum cellular automata cell design trade-offs: latching vs. power dissipation.

The use of molecules to enact quantum cellular automata (QCA) cells has been proposed as a new way for performing electronic logic operations at sub-nm dimensions. A key question that arises concerns whether chemical or physical processes are to be exploited. The use of chemical reactions allows the state of a switch element to be latched in molecular form, making the output of a cell independent of its inputs, but costs energy to do the reaction. Alternatively, if purely electronic polarization is manipulated then no internal latching occurs, but no power is dissipated provided the fields from the inputs change slowly compared to the molecular response times. How these scenarios pan out is discussed by considering calculated properties of the 1,4-diallylbutane cation, a species often used as a paradigm for molecular electronic switching. Utilized are results from different calculation approaches that depict the ion either as a charge-localized mixed-valence compound functioning as a bistable switch, or else as an extremely polarizable molecule with a delocalized electronic structure. Practical schemes for using molecular cells in QCA and other devices emerge.

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.

A priori calculations of the free energy of formation from solution of polymorphic self-assembled monolayers.

Modern quantum chemical electronic structure methods typically applied to localized chemical bonding are developed to predict atomic structures and free energies for meso-tetraalkylporphyrin self-assembled monolayer (SAM) polymorph formation from organic solution on highly ordered pyrolytic graphite surfaces. Large polymorph-dependent dispersion-induced substrate-molecule interactions (e.g., -100 kcal mol(-1) to -150 kcal mol(-1) for tetratrisdecylporphyrin) are found to drive SAM formation, opposed nearly completely by large polymorph-dependent dispersion-induced solvent interactions (70-110 kcal mol(-1)) and entropy effects (25-40 kcal mol(-1) at 298 K) favoring dissolution. Dielectric continuum models of the solvent are used, facilitating consideration of many possible SAM polymorphs, along with quantum mechanical/molecular mechanical and dispersion-corrected density functional theory calculations. These predict and interpret newly measured and existing high-resolution scanning tunnelling microscopy images of SAM structure, rationalizing polymorph formation conditions. A wide range of molecular condensed matter properties at room temperature now appear suitable for prediction and analysis using electronic structure calculations.

Challenges facing an understanding of the nature of low-energy excited states in photosynthesis.

While the majority of the photochemical states and pathways related to the biological capture of solar energy are now well understood and provide paradigms for artificial device design, additional low-energy states have been discovered in many systems with obscure origins and significance. However, as low-energy states are naively expected to be critical to function, these observations pose important challenges. A review of known properties of low energy states covering eight photochemical systems, and options for their interpretation, are presented. A concerted experimental and theoretical research strategy is suggested and outlined, this being aimed at providing a fully comprehensive understanding.

Intermixed adatom and surface-bound adsorbates in regular self-assembled monolayers of racemic 2-butanethiol on Au(111).

In situ scanning tunneling microscopy combined with density functional theory molecular dynamics simulations reveal a complex structure for the self-assembled monolayer (SAM) of racemic 2-butanethiol on Au(111) in aqueous solution. Six adsorbate molecules occupy a (10×√3)R30° cell organized as two RSAuSR adatom-bound motifs plus two RS species bound directly to face-centered-cubic and hexagonally close-packed sites. This is the first time that these competing head-group arrangements have been observed in the same ordered SAM. Such unusual packing is favored as it facilitates SAMs with anomalously high coverage (30%), much larger than that for enantiomerically resolved 2-butanethiol or secondary-branched butanethiol (25%) and near that for linear-chain 1-butanethiol (33%).

Electron-vibration entanglement in the Born-Oppenheimer description of chemical reactions and spectroscopy.

Entanglement is sometimes regarded as the quintessential measure of the quantum nature of a system and its significance for the understanding of coupled electronic and vibrational motions in molecules has been conjectured. Previously, we considered the entanglement developed in a spatially localized diabatic basis representation of the electronic states, considering design rules for qubits in a low-temperature chemical quantum computer. We extend this to consider the entanglement developed during high-energy processes. We also consider the entanglement developed using adiabatic electronic basis, providing a novel way for interpreting effects of the breakdown of the Born-Oppenheimer (BO) approximation. We consider: (i) BO entanglement in the ground-state wavefunction relevant to equilibrium thermodynamics, (ii) BO entanglement associated with low-energy wavefunctions relevant to infrared and tunneling spectroscopies, (iii) BO entanglement in high-energy eigenfunctions relevant to chemical reaction processes, and (iv) BO entanglement developed during reactive wavepacket dynamics. A two-state single-mode diabatic model descriptive of a wide range of chemical phenomena is used for this purpose. The entanglement developed by BO breakdown correlates simply with the diameter of the cusp introduced by the BO approximation, and a hierarchy appears between the various BO-breakdown correction terms, with the first-derivative correction being more important than the second-derivative correction which is more important than the diagonal correction. This simplicity is in contrast to the complexity of BO-breakdown effects on thermodynamic, spectroscopic, and kinetic properties. Further, processes poorly treated at the BO level that appear adequately treated using the Born-Huang adiabatic approximation are found to have properties that can only be described using a non-adiabatic description. For the entanglement developed between diabatic electronic states and the nuclear motion, qualitatively differently behavior is found compared to traditional properties of the density matrix and hence entanglement provides new information about system properties. For chemical reactions, this type of entanglement simply builds up as the transition-state region is crossed. It is robust to small changes in parameter values and is therefore more attractive for making quantum qubits than is the related fragile ground-state entanglement, provided that coherent motion at the transition state can be sustained.

A unified diabatic description for electron transfer reactions, isomerization reactions, proton transfer reactions, and aromaticity.

While diabatic approaches are ubiquitous for the understanding of electron-transfer reactions and have been mooted as being of general relevance, alternate applications have not been able to unify the same wide range of observed spectroscopic and kinetic properties. The cause of this is identified as the fundamentally different orbital configurations involved: charge-transfer phenomena involve typically either 1 or 3 electrons in two orbitals whereas most reactions are typically closed shell. As a result, two vibrationally coupled electronic states depict charge-transfer scenarios whereas three coupled states arise for closed-shell reactions of non-degenerate molecules and seven states for the reactions implicated in the aromaticity of benzene. Previous diabatic treatments of closed-shell processes have considered only two arbitrarily chosen states as being critical, mapping these states to those for electron transfer. We show that such effective two-state diabatic models are feasible but involve renormalized electronic coupling and vibrational coupling parameters, with this renormalization being property dependent. With this caveat, diabatic models are shown to provide excellent descriptions of the spectroscopy and kinetics of the ammonia inversion reaction, proton transfer in N2H7(+), and aromaticity in benzene. This allows for the development of a single simple theory that can semi-quantitatively describe all of these chemical phenomena, as well as of course electron-transfer reactions. It forms a basis for understanding many technologically relevant aspects of chemical reactions, condensed-matter physics, chemical quantum entanglement, nanotechnology, and natural or artificial solar energy capture and conversion.

Non-adiabatic effects in thermochemistry, spectroscopy and kinetics: the general importance of all three Born-Oppenheimer breakdown corrections.

Using a simple model Hamiltonian, the three correction terms for Born-Oppenheimer (BO) breakdown, the adiabatic diagonal correction (DC), the first-derivative momentum non-adiabatic correction (FD), and the second-derivative kinetic-energy non-adiabatic correction (SD), are shown to all contribute to thermodynamic and spectroscopic properties as well as to thermal non-diabatic chemical reaction rates. While DC often accounts for >80% of thermodynamic and spectroscopic property changes, the commonly used practice of including only the FD correction in kinetics calculations is rarely found to be adequate. For electron-transfer reactions not in the inverted region, the common physical picture that diabatic processes occur because of surface hopping at the transition state is proven inadequate as the DC acts first to block access, increasing the transition state energy by (ℏω)(2)λ/16J(2) (where λ is the reorganization energy, J the electronic coupling and ω the vibration frequency). However, the rate constant in the weakly-coupled Golden-Rule limit is identified as being only inversely proportional to this change rather than exponentially damped, owing to the effects of tunneling and surface hopping. Such weakly-coupled long-range electron-transfer processes should therefore not be described as "non-adiabatic" processes as they are easily described by Born-Huang ground-state adiabatic surfaces made by adding the DC to the BO surfaces; instead, they should be called just "non-Born-Oppenheimer" processes. The model system studied consists of two diabatic harmonic potential-energy surfaces coupled linearly through a single vibration, the "two-site Holstein model". Analytical expressions are derived for the BO breakdown terms, and the model is solved over a large parameter space focusing on both the lowest-energy spectroscopic transitions and the quantum dynamics of coherent-state wavepackets. BO breakdown is investigated pertinent to: ammonia inversion, aromaticity in benzene, the Creutz-Taube ion, the bacterial photosynthetic reaction centre, BNB, the molecular conductor Alq3, and inverted-region charge recombination in a ferrocene-porphyrin-fullerene triad photosynthetic model compound. Throughout, the fundamental nature of BO breakdown is linked to the properties of the cusp catastrophe: the cusp diameter is shown to determine the magnitudes of all couplings, numerical basis-set and trajectory-integration requirements, and to determine the transmission coefficient κ used to understand deviations from transition-state theory.

Bond angle variations in XH3 [X = N, P, As, Sb, Bi]: the critical role of Rydberg orbitals exposed using a diabatic state model.

Ammonia adopts sp(3) hybridization (HNH bond angle 108°) whereas the other members of the XH3 series PH3, AsH3, SbH3, and BiH3 instead prefer octahedral bond angles of 90-93°. We use a recently developed general diabatic description for closed-shell chemical reactions, expanded to include Rydberg states, to understand the geometry, spectroscopy and inversion reaction profile of these molecules, fitting its parameters to results from Equation of Motion Coupled-Cluster Singles and Doubles (EOM-CCSD) calculations using large basis sets. Bands observed in the one-photon absorption spectrum of NH3 at 18.3 eV, 30 eV, and 33 eV are reassigned from Rydberg (formally forbidden) double excitations to valence single-excitation resonances. Critical to the analysis is the inclusion of all three electronic states in which two electrons are placed in the lone-pair orbital n and/or the symmetric valence σ* antibonding orbital. An illustrative effective two-state diabatic model is also developed containing just three parameters: the resonance energy driving the high-symmetry planar structure, the reorganization energy opposing it, and HXH bond angle in the absence of resonance. The diabatic orbitals are identified as sp hybrids on X; for the radical cations XH3(+) for which only 2 electronic states and one conical intersection are involved, the principle of orbital following dictates that the bond angle in the absence of resonance is acos(-1/5) = 101.5°. The multiple states and associated multiple conical intersection seams controlling the ground-state structure of XH3 renormalize this to acos[3 sin(2)(2(1/2)atan(1/2))/2 - 1/2] = 86.7°. Depending on the ratio of the resonance energy to the reorganization energy, equilibrium angles can vary from these limiting values up to 120°, and the anomalously large bond angle in NH3 arises because the resonance energy is unexpectedly large. This occurs as the ordering of the lowest Rydberg orbital and the σ* orbital swap, allowing Rydbergization to compresses σ* to significantly increase the resonance energy. Failure of both the traditional and revised versions of the valence-shell electron-pair repulsion (VSEPR) theory to explain the ground-state structures in simple terms is attributed to exclusion of this key physical interaction.