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.

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.

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.

Free-Energy Landscape and Isomerization Rates of Au 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.

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 (Δ) on temperature in both the solid and the molecule, and a substantial effect of packing. While the Δ vanishes in the perfect crystal, it is of the order of ∼0.

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.

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 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.

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.

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.

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.

A Finite-Field Approach for GW Calculations beyond the Random Phase Approximation.

We describe a finite-field approach to compute density response functions, which allows for efficient G W and G WΓ calculations beyond the random phase approximation. The method is easily applicable to density functional calculations performed with hybrid functionals. We present results for the electronic properties of molecules and solids, and we discuss a general scheme to overcome slow convergence of quasiparticle energies obtained from G WΓ calculations, as a function of the basis set used to represent the dielectric matrix.

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 WO, 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.

First-Principles Simulations of Liquid Water Using a Dielectric-Dependent Hybrid Functional.

We carried out first-principles simulations of liquid water under ambient conditions using a dielectric-dependent hybrid functional, where the fraction of exact exchange is set equal to the inverse of the high-frequency dielectric constant of the liquid. We found excellent agreement with experiment for the oxygen-oxygen partial correlation function at the experimental equilibrium density and 311 ± 3 K. Other structural and dynamical properties, such as the diffusion coefficient, molecular dipole moments, and vibrational spectra, are also in good agreement with experiment.

Hierarchical Coupling of First-Principles Molecular Dynamics with Advanced Sampling Methods.

We present a seamless coupling of a suite of codes designed to perform advanced sampling simulations, with a first-principles molecular dynamics (MD) engine. As an illustrative example, we discuss results for the free energy and potential surfaces of the alanine dipeptide obtained using both local and hybrid density functionals (DFT), and we compare them with those of a widely used classical force field, Amber99sb. In our calculations, the efficiency of first-principles MD using hybrid functionals is augmented by hierarchical sampling, where hybrid free energy calculations are initiated using estimates obtained with local functionals.

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.

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.

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.

Reproducibility in density functional theory calculations of solids.

The widespread popularity of density functional theory has given rise to an extensive range of dedicated codes for predicting molecular and crystalline properties. However, each code implements the formalism in a different way, raising questions about the reproducibility of such predictions. We report the results of a community-wide effort that compared 15 solid-state codes, using 40 different potentials or basis set types, to assess the quality of the Perdew-Burke-Ernzerhof equations of state for 71 elemental crystals.

Performance and accuracy of recursive subspace bisection for hybrid DFT calculations in inhomogeneous systems.

The high cost of computing the Hartree-Fock exchange energy has resulted in a limited use of hybrid density functionals in solid-state and condensed phase calculations. Approximate methods based on the use of localized orbitals have been proposed as a way to reduce this computational cost. In particular, Boys orbitals (or maximally localized Wannier functions in solids) were recently used in plane wave, first-principles molecular dynamics simulations of water.

Density and Compressibility of Liquid Water and Ice from First-Principles Simulations with Hybrid Functionals.

We determined the equilibrium density and compressibility of water and ice from first-principles molecular dynamics simulations using gradient-corrected (PBE) and hybrid (PBE0) functionals. Both functionals predicted the density of ice to be larger than that of water, by 15 (PBE) and 35% (PBE0). The PBE0 functional yielded a lower density of both ice and water with respect to PBE, leading to better agreement with experiment for ice but not for liquid water.

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.