Quantum simulations of thermally activated delayed fluorescence in an all-organic emitter.

We investigate the prototypical NAI-DMAC thermally activated delayed fluorescence (TADF) emitter in the gas phase- and high-packing fraction limits at finite temperature, by combining first principles molecular dynamics with a quantum thermostat to account for nuclear quantum effects (NQE). We find a weak dependence of the singlet-triplet energy gap (ΔEST) on temperature in both the solid and the molecule, and a substantial effect of packing. While the ΔEST vanishes in the perfect crystal, it is of the order of ∼0.3 eV in the molecule, with fluctuations ranging from 0.1 to 0.4 eV at 300 K. The transition probability between the HOMOs and LUMOs has a stronger dependence on temperature than the singlet-triplet gap, with a desirable effect for thermally activated fluorescence; such temperature effect is weaker in the condensed phase than in the molecule. Our results on ΔEST and oscillator strengths, together with our estimates of direct and reverse intersystem crossing rates, show that optimization of packing and geometrical conformation is critical to increase the efficiency of TADF compounds. Our findings highlight the importance of considering thermal fluctuations and NQE to obtain robust predictions of the electronic properties of NAI-DMAC.

Free-Energy Landscape and Isomerization Rates of Au4 Clusters at Finite Temperatures.

In metallic nanoparticles, the geometry of atomic positions controls the particle's electronic band structure, polarizability, and catalytic properties. Analyzing the structural properties is a complex problem; the structure of an assembled cluster changes from moment to moment due to thermal fluctuations. Conventional structural analyses based on spectroscopy or diffraction cannot determine the instantaneous structure exactly and can merely provide an averaged structure. Molecular simulations offer an opportunity to examine the assembly and evolution of metallic clusters, as the preferred assemblies and conformations can easily be visualized and explored. Here, we utilize the adaptive biasing force algorithm applied to first-principles molecular dynamics to demonstrate the exploration of a relatively simple system, which permits a comprehensive study of the small metal cluster Au4 in both neutral and charged configurations. Our simulation work offers a quantitative understanding of these clusters' dynamic structure, which is significant for single-site catalytic reactions on metal clusters and provides a starting point for a detailed quantitative understanding of more complex pure metal and alloy clusters' dynamic properties.

Capturing the Complete Reaction Profile of a C-H Bond Activation.

The activation of C-H bonds requires the generation of extremely reactive species, which hinders the study of this reaction and its key intermediates. To overcome this challenge, we synthesized an iron(III) chloride-pyridinediimine complex that generates a chlorine radical proximate to reactive C-H bonds upon irradiation with light. Transient spectroscopy confirms the formation of a Cl·|arene complex, which then activates C-H bonds on the PDI ligand to yield HCl and a carbon-centered radical as determined by photocrystallography. First-principles molecular dynamics-density functional theory calculations reveal the trajectory for the formation of a Cl·|arene intermediate. Together, these experimental and computational results show the complete reaction profile for the preferential activation of a C-H bond in the solid state.

The role of defects and excess surface charges at finite temperature for optimizing oxide photoabsorbers.

Computational screening of materials for solar to fuel conversion technologies has mostly focused on bulk properties, thus neglecting the structure and chemistry of surfaces and interfaces with water. We report a finite temperature study of WO3, a promising anode for photoelectrochemical cells, carried out using first-principles molecular dynamics simulations coupled with many-body perturbation theory. We identified three major factors determining the chemical reactivity of the material interfaced with water: the presence of surface defects, the dynamics of excess charge at the surface, and finite temperature fluctuations of the surface electronic orbitals. These general descriptors are essential for the understanding and prediction of optimal oxide photoabsorbers for water oxidation.

Finite-Field Approach to Solving the Bethe-Salpeter Equation.

We present a method to compute optical spectra and exciton binding energies of molecules and solids based on the solution of the Bethe-Salpeter equation and the calculation of the screened Coulomb interaction in a finite field. The method does not require either the explicit evaluation of dielectric matrices or of virtual electronic states, and can be easily applied without resorting to the random phase approximation. In addition, it utilizes localized orbitals obtained from Bloch states using bisection techniques, thus greatly reducing the complexity of the calculation and enabling the efficient use of hybrid functionals to obtain single particle wave functions. We report exciton binding energies of several molecules and absorption spectra of condensed systems of unprecedented size, including water and ice samples with hundreds of atoms.

Ensemble first-principles molecular dynamics simulations of water using the SCAN meta-GGA density functional.

We present an ensemble of 16 independent first-principles molecular dynamics simulations of water performed using the Strongly Constrained and Appropriately Normed (SCAN) meta-generalized gradient approximation exchange-correlation functional. These simulations were used to compute the structural and electronic properties of liquid water, as well as polarizabilities, Raman and infrared spectra. Overall, we find that the SCAN functional used at a simulation temperature of 330 K provides an accurate description of the structural and electronic properties of water while incurring a moderate computational cost. The availability of an ensemble of independent simulations provides a quantitative estimate of the uncertainty in computed structural and electronic properties. Results are also compared with a similar dataset generated using the Perdew, Burke, and Ernzerhof exchange-correlation functional at a temperature of 400 K. All simulation data and trajectories are available at http://quantum-simulation.org.

Roadmap on semiconductor-cell biointerfaces.

This roadmap outlines the role semiconductor-based materials play in understanding the complex biophysical dynamics at multiple length scales, as well as the design and implementation of next-generation electronic, optoelectronic, and mechanical devices for biointerfaces. The roadmap emphasizes the advantages of semiconductor building blocks in interfacing, monitoring, and manipulating the activity of biological components, and discusses the possibility of using active semiconductor-cell interfaces for discovering new signaling processes in the biological world.

Equilibration and analysis of first-principles molecular dynamics simulations of water.

First-principles molecular dynamics (FPMD) simulations based on density functional theory are becoming increasingly popular for the description of liquids. In view of the high computational cost of these simulations, the choice of an appropriate equilibration protocol is critical. We assess two methods of estimation of equilibration times using a large dataset of first-principles molecular dynamics simulations of water. The Gelman-Rubin potential scale reduction factor [A. Gelman and D. B. Rubin, Stat. Sci. 7, 457 (1992)] and the marginal standard error rule heuristic proposed by White [Simulation 69, 323 (1997)] are evaluated on a set of 32 independent 64-molecule simulations of 58 ps each, amounting to a combined cumulative time of 1.85 ns. The availability of multiple independent simulations also allows for an estimation of the variance of averaged quantities, both within MD runs and between runs. We analyze atomic trajectories, focusing on correlations of the Kohn-Sham energy, pair correlation functions, number of hydrogen bonds, and diffusion coefficient. The observed variability across samples provides a measure of the uncertainty associated with these quantities, thus facilitating meaningful comparisons of different approximations used in the simulations. We find that the computed diffusion coefficient and average number of hydrogen bonds are affected by a significant uncertainty in spite of the large size of the dataset used. A comparison with classical simulations using the TIP4P/2005 model confirms that the variability of the diffusivity is also observed after long equilibration times. Complete atomic trajectories and simulation output files are available online for further analysis.

Blind test of density-functional-based methods on intermolecular interaction energies.

In the past decade, a number of approaches have been developed to fix the failure of (semi)local density-functional theory (DFT) in describing intermolecular interactions. The performance of several such approaches with respect to highly accurate benchmarks is compared here on a set of separation-dependent interaction energies for ten dimers. Since the benchmarks were unknown before the DFT-based results were collected, this comparison constitutes a blind test of these methods.

Solar nanocomposites with complementary charge extraction pathways for electrons and holes: Si embedded in ZnS.

We propose that embedding silicon nanoparticles (NP) into amorphous, nonstoichiometric ZnS leads to promising nanocomposites for solar energy conversion. Using ab initio molecular dynamics simulations we show that, upon high temperature amorphization of the host chalcogenide, sulfur atoms are drawn to the NP surface. We find that the sulfur content may be engineered to form a type II heterojunction, with complementary charge transport channels for electrons and holes, and that sulfur capping is beneficial to lower the nanoparticle gap, with respect to that of NPs embedded in oxide matrices. Our analysis is conducted using density functional theory with local and hybrid functionals and many body perturbation theory at the GW level.

Raman Spectra of Electrified Si-Water Interfaces: First-Principles Simulations.

We investigate the Raman spectra of liquid water in contact with a semiconductor surface using first-principles molecular dynamics simulations. We focus on a hydrogenated silicon-water interface and compute the Raman spectra from time correlation functions of the polarizability. We establish a relationship between Raman spectral signatures and structural properties of the liquid at the interface, and we identify the vibrational impacts of an applied electric field. We show that negative bias leads to a reduction of the number of hydrogen bonds (HBs) formed between the surface and the topmost water layer and an enhancement of the HB interactions between water molecules. Instead, positive bias leads to an enhancement of both the HB interactions between water and the surface and between water molecules, creating a semi-ordered interfacial layer. Our work provides molecular-level insights into electrified semiconductor/water interfaces and the identification of specific structural features through Raman spectroscopy.

Communication: electronic structure of the solvated chloride anion from first principles molecular dynamics.

We present first principles molecular dynamics simulations of the chloride anion in liquid water performed using gradient-corrected and hybrid density functionals. We show that it is necessary to use hybrid functionals both for the generation of molecular dynamics trajectories and for the calculation of electronic states in order to obtain a qualitatively correct description of the electronic properties of the solution. In particular, it is only with hybrid functionals that the highest occupied molecular orbital of the anion is found above the valence band maximum of water, consistent with photoelectron detachment measurements. Similar results were obtained using many body perturbation theory within the G0W0 approximation.

All-Electron Plane-Wave Electronic Structure Calculations.

We demonstrate the use of the plane wave basis for all-electron electronic structure calculations. The approach relies on the definition of an analytic, norm-conserving, regularized Coulomb potential, and a scalable implementation of the plane wave method capable of handling large energy cutoffs (up to 80 kRy in the examples shown). The method is applied to the computation of electronic properties of isolated atoms as well as the diamond and silicon crystals, MgO, solid argon, and a configuration of 64 water molecules extracted from a first-principles molecular dynamics simulation. The computed energies, band gaps, ionic forces, and stress tensors provide reference results for the validation of pseudopotentials and/or localized basis sets. A calculation of the all-electron band structure of diamond and silicon using the SCAN meta-GGA density functional allows for a validation of calculations based on pseudopotentials derived using the PBE exchange-correlation functional. In the case of (H2O)64, the computed ionic forces provide a reference from which the errors incurred in pseudopotential calculations and in localized Gaussian basis sets calculations can be estimated.

Engineering the formation of spin-defects from first principles.

The full realization of spin qubits for quantum technologies relies on the ability to control and design the formation processes of spin defects in semiconductors and insulators. We present a computational protocol to investigate the synthesis of point-defects at the atomistic level, and we apply it to the study of a promising spin-qubit in silicon carbide, the divacancy (VV). Our strategy combines electronic structure calculations based on density functional theory and enhanced sampling techniques coupled with first principles molecular dynamics. We predict the optimal annealing temperatures for the formation of VVs at high temperature and show how to engineer the Fermi level of the material to optimize the defect's yield for several polytypes of silicon carbide. Our results are in excellent agreement with available experimental data and provide novel atomistic insights into point defect formation and annihilation processes as a function of temperature.

A simplified implementation of van der Waals density functionals for first-principles molecular dynamics applications.

We present a simplified implementation of the non-local van der Waals correlation functional introduced by Dion et al. [Phys. Rev. Lett. 92, 246401 (2004)] and reformulated by Román-Pérez et al. [Phys. Rev. Lett. 103, 096102 (2009)]. The proposed numerical approach removes the logarithmic singularity of the kernel function. Complete expressions of the self-consistent correlation potential and of the stress tensor are given. Combined with various choices of exchange functionals, five versions of van der Waals density functionals are implemented. Applications to the computation of the interaction energy of the benzene-water complex and to the computation of the equilibrium cell parameters of the benzene crystal are presented. As an example of crystal structure calculation involving a mixture of hydrogen bonding and dispersion interactions, we compute the equilibrium structure of two polymorphs of aspirin (2-acetoxybenzoic acid, C(9)H(8)O(4)) in the P2(1)/c monoclinic structure.

Tailored nanoheterojunctions for optimized light emission.

We present coupled classical and quantum simulations of 1 to 2 nm Si nanocrystals (NCs) embedded in amorphous SiO(2) and we show that by tuning the density of the oxide matrix one may change the relative alignment of Si NC and SiO(2) electronic states at the interface. We find that interfacial strain plays a key role in determining the variation of the nanaoparticle gap as a function of size, as well as of conduction band offsets with the oxide. In particular, our results show that it is the variation of the valence band offset with size that is responsible for the gap change. Our findings suggest that the elastic properties of the embedding matrix may be tuned to tailor the energy levels of small Si NCs so as to optimize their performance in optoelectronic devices and solar cells.

Melting of ice under pressure.

The melting of ice under pressure is investigated with a series of first-principles molecular dynamics simulations. In particular, a two-phase approach is used to determine the melting temperature of the ice-VII phase in the range of 10-50 GPa. Our computed melting temperatures are consistent with existing diamond anvil cell experiments. We find that for pressures between 10 and 40 GPa, ice melts as a molecular solid. For pressures above approximately 45 Gpa, there is a sharp increase in the slope of the melting curve because of the presence of molecular dissociation and proton diffusion in the solid before melting. The onset of significant proton diffusion in ice-VII as a function of increasing temperature is found to be gradual and bears many similarities to that of a type-II superionic solid.

Neural Network Sampling of the Free Energy Landscape for Nitrogen Dissociation on Ruthenium.

In heterogeneous catalysis, free energy profiles of reactions govern the mechanisms, rates, and equilibria. Energetics are conventionally computed using the harmonic approximation (HA), which requires determination of critical states a priori. Here, we use neural networks to efficiently sample and directly calculate the free energy surface (FES) of a prototypical heterogeneous catalysis reaction-the dissociation of molecular nitrogen on ruthenium-at density-functional-theory-level accuracy. We find that the vibrational entropy of surface atoms, often neglected in HA for transition metal catalysts, contributes significantly to the reaction barrier. The minimum free energy path for dissociation reveals an "on-top" adsorbed molecular state prior to the transition state. While a previously reported flat-lying molecular metastable state can be identified in the potential energy surface, it is absent in the FES at relevant reaction temperatures. These findings demonstrate the importance of identifying critical points self-consistently on the FES for reactions that involve considerable entropic effects.

Dissociation of salts in water under pressure.

The investigation of salts in water at extreme conditions is crucial to understanding the properties of aqueous fluids in the Earth. We report first principles (FP) and classical molecular dynamics simulations of NaCl in the dilute limit, at temperatures and pressures relevant to the Earth's upper mantle. Similar to ambient conditions, we observe two metastable states of the salt: the contact (CIP) and the solvent-shared ion-pair (SIP), which are entropically and enthalpically favored, respectively. We find that the free energy barrier between the CIP and SIP minima increases at extreme conditions, and that the stability of the CIP is enhanced in FP simulations, consistent with the decrease of the dielectric constant of water. The minimum free energy path between the CIP and SIP becomes smoother at high pressure, and the relative stability of the two configurations is affected by water self-dissociation, which can only be described properly by FP simulations.