Electronic structure calculations with GPAW: a real-space implementation of the projector augmented-wave method

Halide-sodalites: thermal behavior at low temperatures and local deviations from the average structure

Sodalites of the general type |Na8X2|[T1T2O4]6 with X = Cl−, Br−, I− have been synthesized for Al–Si, Ga–Si, Al–Ge and Ga–Ge as T1–T2 frameworks. The structures were examined using in-house and synchrotron X-ray diffraction, Raman spectroscopy, force-field structure optimizations and DFT based ab-initio molecular dynamics (MD) computations. Calculated phonon density of states (PDOS) of the 12 compounds show only minor differences within a framework composition with a lowering of certain phonon energies with increasing anion size. Earlier published Debye and Einstein temperatures obtained with a Debye-Einstein-anharmonicity (DEA) model approach are confirmed using the determined low-temperature lattice parameters (18 K–293 K) and show no correlation with the respective PDOS. Small-box refinements against radial pair distribution functions (PDF) allowed the determination of anisotropic displacement ellipsoids (ADP) for Na+ and O2−, indicating a strong dependency of the ADP of Na+ on the chemical composition. Significantly lower thermal displacements from MD calculations suggested an influence of structural displacements. For compounds with an aspherical ADP for sodium, structural models could be refined in which the sodium is located on two 8e or one 24i site (both partially occupied), and also temperature-dependent (100 K–300 K) for the compounds with Ga–Ge framework. 3D-plots of the bond-valence sums of Na+ further validate the structural differences. These results imply that the local structure of halide-sodalites in many cases is not best described by the known average structure and may even not be cubic.

Electronic and Optical Properties of Rocksalt Mg1-xZnxO and Wurtzite Zn1-xMgxO with Varied Concentrations of Magnesium and Zinc

The structural, electronic and optical properties of rocksalt Mg1-xZnxO and wurtzite Zn1-xMgxO with the concentration of Zn and Mg varying from 0.125 to 0.875 were investigated using density functional theory (DFT), DFT+U, linear response theory and the Bethe-Salpeter equation. According to the experimental band gap for varied concentrations of magnesium and zinc, modeling the supercell was utilized for the varied concentrations of Mg/Zn/O compounds in order to not only avoid constructing the complicated interface systems that are observed in the experiments but also take into account the excitonic effects that usually require huge computational resources. From the calculated density of states, the Zn states are highly related to the edge of the conduction band minimum and responsible for the width of bandgap. In addition, the contribution of Zn-d states is below expectations as they are located away from the VBM. As for the optical response, an increase in Zn concentration would cause a red-shifted spectrum, on the whole. In contrast, the higher concentration of Mg also triggers the blue-shift of the optical spectrum. In addition, anisotropic properties could be found in the spectrum with consideration of the excitonic effects, whereas there is no apparent difference in optical response based on linear response theory. In addition, the optical features of this work reflect the characteristic peaks of the literature around the absorption onset.

Density-of-states similarity descriptor for unsupervised learning from materials data

We develop a materials descriptor based on the electronic density-of-states (DOS) and investigate the similarity of materials based on it. As an application example, we study the Computational 2D Materials Database (C2DB) that hosts thousands of two-dimensional materials with their properties calculated by density-functional theory. Combining our descriptor with a clustering algorithm, we identify groups of materials with similar electronic structure. We introduce additional descriptors to characterize these clusters in terms of crystal structures, atomic compositions, and electronic configurations of their members. This allows us to rationalize the found (dis)similarities and to perform an automated exploratory and confirmatory analysis of the C2DB data. From this analysis, we find that the majority of clusters consist of isoelectronic materials sharing crystal symmetry, but we also identify outliers, i.e., materials whose similarity cannot be explained in this way.

Indirect Band Gap Semiconductors for Thin-Film Photovoltaics: High-Throughput Calculation of Phonon-Assisted Absorption

Discovery of high-performance materials remains one of the most active areas in photovoltaics (PV) research. Indirect band gap materials form the largest part of the semiconductor chemical space, but predicting their suitability for PV applications from first-principles calculations remains challenging. Here, we propose a computationally efficient method to account for phonon-assisted absorption across the indirect band gap and use it to screen 127 experimentally known binary semiconductors for their potential as thin-film PV absorbers. Using screening descriptors for absorption, carrier transport, and nonradiative recombination, we identify 28 potential candidate materials. The list, which contains 20 indirect band gap semiconductors, comprises well-established (3), emerging (16), and previously unexplored (9) absorber materials. Most of the new compounds are anion-rich chalcogenides (TiS₃ and Ga₂Te₅) and phosphides (PdP₂, CdP₄, MgP₄, and BaP₃) containing homoelemental bonds and represent a new frontier in PV materials research. Our work highlights the previously underexplored potential of indirect band gap materials for optoelectronic thin-film technologies.

Force Dependent Barriers from Analytic Potentials within Elastic Environments

Bond rupture under the action of external forces is usually induced by temperature fluctuations, where the key quantity is the force dependent barrier that needs to be overcome. Using analytic potentials we find that these barriers are fully determined by the dissociation energy and the maximal force the potential can withstand. The barrier shows a simple dependence on these two quantities that allows for a re-interpretation of the Eyring-Zhurkov-Bell length Δ x ‡ and the expressions in theories going beyond that. It is shown that solely elastic environments do not change this barrier in contrast to the predictions of constraint geometry simulate external force (COGEF) strategies. The findings are confirmed by explicit calculations of bond rupture in a polydimethylsiloxane model.

Influence of Stress on Electronic and Optical Properties of Rocksalt and Wurtzite MgO-ZnO Nanocomposites with Varying Concentrations of Magnesium and Zinc

The structural, electronic and optical properties of stressed MgO-ZnO nanocomposite alloys with concentrations of Zn and Mg varying from 0.125 to 0.875 were studied using ab initio simulations. Two crystal structures are considered for the initial MgO-ZnO alloys: the rocksalt Mg1-xZnxO and wurtzite Zn1-xMgxO phases. For rocksalt Mg1-xZnxO, the optimized structures are stable at pressures below 10 GPa. The larger the Mg concentration and pressure, the wider the Eg of the rocksalt phase. In contrast, the optimal geometries of wurtzite Zn1-xMgxO reveal a diversity of possibilities, including rocksalt, wurtzite and mixed phases. These effects lead to the fact that the optical properties of wurtzite Zn1-xMgxO not only demonstrate the properties of the wurtzite phase but also indicate the optical features of the rocksalt phase. In addition, mixed phases of Zn1-xMgxO simultaneously provide the characteristics of both wurtzite and rocksalt phases with the same structures in different dielectric matrices.

Nanofluidic Trapping of Faceted Colloidal Nanocrystals for Parallel Single-Particle Catalysis

Catalyst activity can depend distinctly on nanoparticle size and shape. Therefore, understanding the structure sensitivity of catalytic reactions is of fundamental and technical importance. Experiments with single-particle resolution, where ensemble-averaging is eliminated, are required to study it. Here, we implement the selective trapping of individual spherical, cubic, and octahedral colloidal Au nanocrystals in 100 parallel nanofluidic channels to determine their activity for fluorescein reduction by sodium borohydride using fluorescence microscopy. As the main result, we identify distinct structure sensitivity of the rate-limiting borohydride oxidation step originating from different edge site abundance on the three particle types, as confirmed by first-principles calculations. This advertises nanofluidic reactors for the study of structure-function correlations in catalysis and identifies nanoparticle shape as a key factor in borohydride-mediated catalytic reactions.

Tuning the Optical Properties and Conductivity of Bundles in Networks of Single-Walled Carbon Nanotubes

The films of single-walled carbon nanotubes (SWCNTs) are a promising material for flexible transparent electrodes, which performance depends not only on the properties of individual nanotubes but also, foremost, on bundling of individual nanotubes. This work investigates the impact of densification on optical and electronic properties of SWCNT bundles and fabricated films. Our ab initio analysis shows that the optimally densified bundles, consisting of a mixture of quasi-metallic and semiconducting SWCNTs, demonstrate quasi-metallic behavior and can be considered as an effective conducting medium. Our density functional theory calculations indicate the band curving and bandgap narrowing with the reduction of the distance between nanotubes inside bundles. Simulation results are consistent with the observed conductivity improvement and shift of the absorption peaks in SWCNT films densified in isopropyl alcohol. Therefore, not only individual nanotubes but also the bundles should be considered as building blocks for high-performance transparent conductive SWCNT-based films.

Direct correlation between void formation and lithium dendrite growth in solid-state electrolytes with interlayers

Solid-state Li-ion batteries with lithium anodes offer higher energy densities and are safer than conventional liquid electrolyte-based Li-ion batteries. However, the growth of lithium dendrites across the solid-state electrolyte layer leads to the premature shorting of cells and limits their practical viability. Here, using solid-state Li half-cells with metallic interlayers between a garnet-based lithium-ion conductor and lithium, we show that interfacial void growth precedes dendrite nucleation and growth. Specifically, void growth was observed at a current density of around two-thirds of the critical current density for dendrite growth. Computational calculations reveal that interlayer materials with higher critical current densities for dendrite growth also have the largest thermodynamic and kinetic barriers for lithium vacancy accumulation at their interfaces with lithium. Our results suggest that interfacial modification with suitable metallic interlayers decreases the tendency for void growth and improves dendrite growth tolerance in solid-state electrolytes, even in the absence of high stack pressures.

Momentum-Dependent Oscillator Strength Crossover of Excitons and Plasmons in Two-Dimensional PtSe2

The 1T-phase layered PtX₂ chalcogenide has attracted widespread interest due to its thickness dependent metal-semiconductor transition driven by strong interlayer coupling. While the ground state properties of this paradigmatic material system have been widely explored, its fundamental excitation spectrum remains poorly understood. Here we combine first-principles calculations with momentum (q) resolved electron energy loss spectroscopy (q-EELS) to study the collective excitations in 1T-PtSe₂ from the monolayer limit to the bulk. At finite momentum transfer, all the spectra are dominated by two distinct interband plasmons that disperse to higher energy with increasing q. Interestingly, the absence of long-range screening in the two-dimensional (2D) limit inhibits the formation of long wavelength plasmons. Consequently, in the small-q limit, excitations in monolayer PtSe₂ are exclusively of excitonic nature, and the loss spectrum coincides with the optical spectrum. The qualitatively different momentum dependence of excitons and plasmons enables us to unambiguously disentangle their spectral fingerprints in the excited state spectrum of layered 1T-PtSe₂. This will help to discern the charge carrier plasmon and locally map the optical conductivity and trace the layer-dependent semiconductor to metal transition in 1T-PtSe₂ and other 2D materials.

Rational Designing of Bimetallic/Trimetallic Hydrogen Evolution Reaction Catalysts Using Supervised Machine Learning

Cost-efficient electrocatalysts to replace precious platinum group metals- (PGMs-) based catalysts for the hydrogen evolution reaction (HER) carry significant potential for sustainable energy solutions. Machine learning (ML) methods have provided new avenues for intelligent screening and predicting efficient heterogeneous catalysts in recent years. We coalesce density functional theory (DFT) and supervised ML methods to discover earth-abundant active heterogeneous NiCoCu-based HER catalysts. An intuitive generalized microstructure model was designed to study the adsorbate's surface coverage and generate input features for the ML process. The study utilizes optimized eXtreme Gradient Boost Regression (XGBR) models to screen NiCoCu alloy-based catalysts for HER. We show that the most active HER catalysts can be screened from an extensive set of catalysts with this approach. Therefore, our approach can provide an efficient way to discover novel heterogeneous catalysts for various electrochemical reactions.

WSe2 as Transparent Top Gate for Infrared Near-Field Microscopy

Independent control of carrier density and out-of-plane displacement field is essential for accessing novel phenomena in two-dimensional (2D) material heterostructures. While this is achieved with independent top and bottom metallic gate electrodes in transport experiments, it remains a challenge for near-field optical studies as the top electrode interferes with the optical path. Here, we characterize the requirements for a material to be used as the top-gate electrode and demonstrate experimentally that few-layer WSe₂ can be used as a transparent, ambipolar top-gate electrode in infrared near-field microscopy. We carry out nanoimaging of plasmons in a bilayer graphene heterostructure tuning the plasmon wavelength using a trilayer WSe₂ gate, achieving a density modulation amplitude exceeding 2 × 10¹² cm^-2. The observed ambipolar gate-voltage response allows us to extract the energy gap of WSe₂, yielding a value of 1.05 eV. Our results provide an additional tuning knob to cryogenic near-field experiments on emerging phenomena in 2D materials and moiré heterostructures.

Ultra-flat and long-lived plasmons in a strongly correlated oxide

Plasmons in strongly correlated systems are attracting considerable attention due to their unconventional behavior caused by electronic correlation effects. Recently, flat plasmons with nearly dispersionless frequency-wave vector relations have drawn significant interest because of their intriguing physical origin and promising applications. However, these flat plasmons exist primarily in low-dimensional materials with limited wave vector magnitudes (q < ~0.7 Å⁻¹). Here, we show that long-lived flat plasmons can propagate up to ~1.2 Å⁻¹ in α-Ti₂O₃, a strongly correlated three-dimensional Mott-insulator, with an ultra-small energy fluctuation (<40 meV). The strong correlation effect renormalizes the electronic bands near Fermi level with a small bandwidth, which is responsible for the flat plasmons in α-Ti₂O₃. Moreover, these flat plasmons are not affected by Landau damping over a wide range of wave vectors (q < ~1.2 Å⁻¹) due to symmetry constrains on the electron wavefunctions. Our work provides a strategy for exploring flat plasmons in strongly correlated systems, which in turn may give rise to novel plasmonic devices in which flat and long-lived plasmons are desirable.

Dispersionless plasmons could find important practical applications, but previous demonstrations have been limited to 2D materials and small momentum range. Here the authors report ultra-flat plasmons propagating over a wide range of momenta in a 3D strongly correlated oxide α-Ti₂O₃.

Static Electronic Density Response of Warm Dense Hydrogen: Ab Initio Path Integral Monte Carlo Simulations

The properties of hydrogen under extreme conditions are important for many applications, including inertial confinement fusion and astrophysical models. A key quantity is given by the electronic density response to an external perturbation, which is probed in x-ray Thomson scattering experiments-the state of the art diagnostics from which system parameters like the free electron density n_{e}, the electronic temperature T_{e}, and the charge state Z can be inferred. In this work, we present highly accurate path integral Monte Carlo results for the static electronic density response of hydrogen. We obtain the static exchange-correlation (XC) kernel K_{XC}, which is of central relevance for many applications, such as time-dependent density functional theory. This gives us a first unbiased look into the electronic density response of hydrogen in the warm-dense matter regime, thereby opening up a gamut of avenues for future research.

Polaritonic Chemistry from First Principles via Embedding Radiation Reaction

The coherent interaction of a large collection of molecules with a common photonic mode results in strong light-matter coupling, a feature that has proven highly beneficial for chemistry and has introduced the research topics polaritonic and QED chemistry. Here, we demonstrate an embedding approach to capture the collective nature while retaining the full ab initio representation of single molecules─an approach ideal for polaritonic chemistry. The accuracy of the embedding radiation-reaction ansatz is demonstrated for time-dependent density-functional theory. Then, by virtue of a simple proton-tunneling model, we illustrate that the influence of collective strong coupling on chemical reactions features a nontrivial dependence on the number of emitters and can alternate between strong catalyzing and an inhibiting effect. Bridging classical electrodynamics, quantum optical descriptions, and the ab initio description of realistic molecules, this work can serve as a guiding light for future developments and investigations in the quickly growing fields of QED chemistry and QED material design.

Accurate Computational Prediction of Core-Electron Binding Energies in Carbon-Based Materials: A Machine-Learning Model Combining Density-Functional Theory and GW

We present a quantitatively accurate machine-learning (ML) model for the computational prediction of core-electron binding energies, from which X-ray photoelectron spectroscopy (XPS) spectra can be readily obtained. Our model combines density functional theory (DFT) with GW and uses kernel ridge regression for the ML predictions. We apply the new approach to disordered materials and small molecules containing carbon, hydrogen, and oxygen and obtain qualitative and quantitative agreement with experiment, resolving spectral features within 0.1 eV of reference experimental spectra. The method only requires the user to provide a structural model for the material under study to obtain an XPS prediction within seconds. Our new tool is freely available online through the XPS Prediction Server.

Heteroanionic Control of Exemplary Second-Harmonic Generation and Phase Matchability in 1D LiAsS2-xSex

The isostructural heteroanionic compounds β-LiAsS_2-xSe_x (x = 0, 0.25, 1, 1.75, 2) show a positive correlation between selenium content and second-harmonic response and greatly outperform the industry standard AgGaSe₂. These materials crystallize in the noncentrosymmetric space group Cc as one-dimensional ¹/_∞ [AsQ₂]^- (Q = S, Se, S/Se) chains consisting of corner-sharing AsQ₃ trigonal pyramids with charge-balancing Li⁺ atoms interspersed between the chains. LiAsS_2-xSe_x melts congruently for 0 ≤ x ≤ 1.75, but when the Se content exceeds x = 1.75, crystallization is complicated by a phase transition. This behavior is attributed to the β- to α-phase transition present in LiAsSe₂, which is observed in the Se-rich compositions. The band gap decreases with increasing Se content, starting at 1.63 eV (LiAsS₂) and reaching 1.06 eV (β-LiAsSe₂). Second-harmonic generation measurements as a function of wavelength on powder samples of β-LiAsS_2-xSe_x show that these materials exhibit significantly higher nonlinearity than AgGaSe₂ (d₃₆ = 33 pm/V), reaching a maximum of 61.2 pm/V for LiAsS₂. In comparison, single-crystal measurements for LiAsSSe yielded a d_eff = 410 pm/V. LiAsSSe, LiAsS_0.25Se_1.75, and β-LiAsSe₂ show phase-matching behavior for incident wavelengths exceeding 3 μm. The laser-induced damage thresholds from two-photon absorption processes are on the same order of magnitude as AgGaSe₂, with S-rich materials slightly outperforming AgGaSe₂ and Se-rich materials slightly underperforming AgGaSe₂.

Learning Grain-Boundary Segregation: From First Principles to Polycrystals

The segregation of solute atoms at grain boundaries (GBs) can strongly impact the structural and functional properties of polycrystals. Yet, due to the limited availability of simulation tools to study polycrystals at the atomistic scale (i.e., interatomic potentials), there is a minimal understanding of the variation of solute segregation tendencies across the very complex space of GB microenvironments and the large range of alloys in which it can occur. Here, we develop an algorithmic framework that can directly learn the full spectrum of segregation energies for a metal solute atom in a metal polycrystal from ab initio methods, bypassing the need for alloy interatomic potentials. This framework offers a pathway to a comprehensive catalog of GB solute segregation with quantum accuracy, for the entire alloy space. As an initial demonstration in this pursuit, we build an extensive GB segregation database for aluminum-based alloys across the periodic table, including dozens of alloys for which there are substantially no prior data.

Beam-driven Dynamics of Aluminium Dopants in Graphene

Substituting heteroatoms into graphene can tune its properties for applications ranging from catalysis to spintronics. The further recent discovery that covalent impurities in graphene can be manipulated at atomic precision using a focused electron beam may open avenues towards sub-nanometer device architectures. However, the preparation of clean samples with a high density of dopants is still very challenging. Here, we report vacancy-mediated substitution of aluminium into laser-cleaned graphene, and without removal from our ultra-high vacuum apparatus, study their dynamics under 60 keV electron irradiation using aberration-corrected scanning transmission electron microscopy and spectroscopy. Three- and four-coordinated Al sites are identified, showing excellent agreement with ab initio predictions including binding energies and electron energy-loss spectrum simulations. We show that the direct exchange of carbon and aluminium atoms predicted earlier occurs under electron irradiation, although unexpectedly it is less probable than the same process for silicon. We also observe a previously unknown nitrogen-aluminium exchange that occurs at Al─N double-dopant sites at graphene divacancies created by our plasma treatment.

Machine-learned model Hamiltonian and strength of spin-orbit interaction in strained Mg2X (X = Si, Ge, Sn, Pb)

Machine-learned multi-orbital tight-binding (MMTB) Hamiltonian models have been developed to describe the electronic characteristics of intermetallic compounds Mg₂Si, Mg₂Ge, Mg₂Sn, and Mg₂Pb subject to strain. The MMTB models incorporate spin-orbital mediated interactions and they are calibrated to the electronic band structures calculated via density functional theory by a massively parallelized multi-dimensional Monte-Carlo search algorithm. The results show that a machine-learned five-band tight-binding (TB) model reproduces the key aspects of the valence band structures in the entire Brillouin zone. The five-band model reveals that compressive strain localizes the contribution of the 3sorbital of Mg to the conduction bands and the outer shellporbitals of X (X = Si, Ge, Sn, Pb) to the valence bands. In contrast, tensile strain has a reversed effect as it weakens the contribution of the 3sorbital of Mg and the outer shellporbitals of X to the conduction bands and valence bands, respectively. Theπbonding in the Mg₂X compounds is negligible compared to theσbonding components, which follow the hierarchy|σsp|>|σpp|>|σss|, and the largest variation against strain belongs toσ_pp. The five-band model allows for estimating the strength of spin-orbit coupling (SOC) in Mg₂X and obtaining its dependence on the atomic number of X and strain. Further, the band structure calculations demonstrate a significant band gap tuning and band splitting due to strain. A compressive strain of-10%can open a band gap at the Γ point in metallic Mg₂Pb, whereas a tensile strain of+10%closes the semiconducting band gap of Mg₂Si. A tensile strain of+5%removes the three-fold degeneracy of valence bands at the Γ point in semiconducting Mg₂Ge. The presented MMTB models can be extended for various materials and simulations (band structure, transport, classical molecular dynamics), and the obtained results can help in designing devices made of Mg₂X.

Machine Learning Assisted Exploration of High Entropy Alloy-Based Catalysts for Selective CO2 Reduction to Methanol

Catalytic conversion of CO₂ to carbon neutral fuels can be ecofriendly and allow for economic replacement of fossil fuels. Here, we have investigated high-throughput screening of high entropy alloy (Cu, Co, Ni, Zn, and Sn) based catalysts through machine learning (ML) for CO₂ hydrogenation to methanol. Stability and catalytic activity studies of these catalysts have been performed for all possible combinations, where different elemental, compositional, and surface microstructural features were used as input parameters. Adsorption energy values of CO₂ reduction intermediates on the CuCoNiZnMg- and CuCoNiZnSn-based catalysts have been used to train the ML models. Successful prediction of adsorption energies of the adsorbates using CuCoNiZnMg-based training data is achieved except for two intermediates. Hence, we show that activity and selectivity of these catalysts can be successfully predicted for CO₂ hydrogenation to methanol and have screened a series of high entropy-based catalysts (from 36750 considered catalysts) which could be promising for methanol synthesis.

Atomistic Global Optimization X: A Python package for optimization of atomistic structures

Modelling and understanding properties of materials from first principles require knowledge of the underlyingatomistic structure. This entails knowing the individual chemical identity and position of all atoms involved.Obtaining such information for macro-molecules, nano-particles, clusters, and for the surface, interface, andbulk phases of amorphous and solid materials represents a difficult high-dimensional global optimizationproblem. The rise of machine learning techniques in materials science has, however, led to many compellingdevelopments that may speed up structure searches. The complexity of such new methods has prompted aneed for an efficient way of assembling them into global optimization algorithms that can be experimentedwith. In this paper, we introduce the Atomistic Global Optimization X (AGOX) framework and code, asa customizable approach that enables efficient building and testing of global optimization algorithms. Amodular way of expressing global optimization algorithms is described and modern programming practicesare used to enable that modularity in the freely available AGOX python package. A number of examplesof global optimization approaches are implemented and analyzed. This ranges from random search andbasin-hopping to machine learning aided approaches with on-the-fly learnt surrogate energy landscapes. Themethods are show-cased on problems ranging from supported clusters over surface reconstructions to largecarbon clusters and metal-nitride clusters incorporated into graphene sheets.

Observation of Biradical Spin Coupling through Hydrogen Bonds

Investigation of intermolecular electron spin interaction is of fundamental importance in both science and technology. Here, radical pairs of all-trans retinoic acid molecules on Au(111) are created using an ultralow temperature scanning tunneling microscope. Antiferromagnetic coupling between two radicals is identified by magnetic-field-dependent spectroscopy. The measured exchange energies are from 0.1 to 1.0 meV. The biradical spin coupling is mediated through O─H⋯O hydrogen bonds, as elucidated from analysis combining density functional theory calculation and a modern version of valence bond theory.

DNA sequencing by Förster resonant energy transfer

We propose a new DNA sequencing concept based on nonradiative Förster resonant energy transfer (FRET) from a donor quantum dot (QD) to an acceptor molecule. The FRET mechanism combined with the nanopore-based DNA translocation is suggested as a novel concept for sequencing DNA molecules. A recently-developed hybrid quantum/classical method is employed, which uses time-dependent density functional theory and quasistatic finite difference time domain calculations. Due to the significant absorbance of DNA bases for photon energies higher than 4 eV, biocompatibility, and stability, we use Zinc-Oxide (ZnO) QD as a donor in the FRET mechanism. The most sensitivity for the proposed method to DNA is achieved for the Hoechst fluorescent-dye acceptor and 1 nm ZnO-QD. Results show that the insertion of each type of DNA nucleobases between the donor and acceptor changes the frequency of the emitted light from the acceptor molecule between 0.25 to 1.6 eV. The noise analysis shows that the method can determine any unknown DNA nucleobases if the signal-to-noise ratio is larger than 5 dB. The proposed concept and excellent results shed light on a new promising class of DNA sequencers.

Breaking with the Principles of Coreduction to Form Stoichiometric Intermetallic PdCu Nanoparticles

Intermetallic nanoparticles (NPs) have shown enhanced catalytic properties as compared to their disordered alloy counterparts. To advance their use in green energy, it is crucial to understand what controls the formation of intermetallic NPs over alloy structures. By carefully selecting the additives used in NP synthesis, it is here shown that monodisperse, intermetallic PdCu NPs can be synthesized in a controllable manner. Introducing the additives iron(III) chloride and ascorbic acid, both morphological and structural control can be achieved. Combined, these additives provide a synergetic effect resulting in precursor reduction and defect-free growth; ultimately leading to monodisperse, single-crystalline, intermetallic PdCu NPs. Using in situ X-ray total scattering, a hitherto unknown transformation pathway is reported that diverges from the commonly reported coreduction disorder-order transformation. A Cu-rich structure initially forms, which upon the incorporation of Pd(0) and atomic ordering forms intermetallic PdCu NPs. These findings underpin that formation of stoichiometric intermetallic NPs is not limited by standard reduction potential matching and coreduction mechanisms, but is instead driven by changes in the local chemistry. Ultimately, using the local chemistry as a handle to tune the NP structure might open new opportunities to expand the library of intermetallic NPs by exploiting synthesis by design.

Rational Design of Biaxial Tensile Strain for Boosting Electronic and Ionic Conductivities of Na2 MnSiO4 for Rechargeable Sodium-Ion Batteries

Using first-principles calculations, biaxial tensile (ϵ=2 and 4 %) and compressive (ϵ=-2 and -4 %) straining of Na2 MnSiO4 lattices resulted into radial distance cut offs of 1.65 and 2 Å, respectively, in the first and second nearest neighbors shell from the center. The Si-O and Mn-O bonds with prominent probability density peaks validated structural stability. Wide-band gap of 2.35 (ϵ=0 %) and 2.54 eV (ϵ=-4 %), and narrow bandgap of 2.24 eV (ϵ=+4 %) estimated with stronger coupling of p-d σ bond than that of the p-d π bond, mainly contributed from the oxygen p-state and manganese d-state. Na+ -ion diffusivity was found to be enhanced by three orders of magnitude as the applied biaxial strain changed from compressive to tensile. According to the findings, the rational design of biaxial strain would improve the ionic and electronic conductivity of Na2 MnSiO4 cathode materials for advanced rechargeable sodium-ion batteries.

Variational Density Functional Calculations of Excited States: Conical Intersection and Avoided Crossing in Ethylene Bond Twisting

Theoretical studies of photochemical processes require a description of the energy surfaces of excited electronic states, especially near degeneracies, where transitions between states are most likely. Systems relevant to photochemical applications are typically too large for high-level multireference methods, and while time-dependent density functional theory (TDDFT) is efficient, it can fail to provide the required accuracy. A variational, time-independent density functional approach is applied to the twisting of the double bond and pyramidal distortion in ethylene, the quintessential model for photochemical studies. By allowing for symmetry breaking, the calculated energy surfaces exhibit the correct topology around the twisted-pyramidalized conical intersection even when using a semilocal functional approximation, and by including explicit self-interaction correction, the torsional energy curves are in close agreement with published multireference results. The findings of the present work point to the possibility of using a single determinant time-independent density functional approach to simulate nonadiabatic dynamics, even for large systems where multireference methods are impractical and TDDFT is often not accurate enough.

Understanding the complementarities of surface-enhanced infrared and Raman spectroscopies in CO adsorption and electrochemical reduction

In situ/operando surface enhanced infrared and Raman spectroscopies are widely employed in electrocatalysis research to extract mechanistic information and establish structure-activity relations. However, these two spectroscopic techniques are more frequently employed in isolation than in combination, owing to the assumption that they provide largely overlapping information regarding reaction intermediates. Here we show that surface enhanced infrared and Raman spectroscopies tend to probe different subpopulations of adsorbates on weakly adsorbing surfaces while providing similar information on strongly binding surfaces by conducting both techniques on the same electrode surfaces, i.e., platinum, palladium, gold and oxide-derived copper, in tandem. Complementary density functional theory computations confirm that the infrared and Raman intensities do not necessarily track each other when carbon monoxide is adsorbed on different sites, given the lack of scaling between the derivatives of the dipole moment and the polarizability. Through a comparison of adsorbed carbon monoxide and water adsorption energies, we suggest that differences in the infrared vs. Raman responses amongst metal surfaces could stem from the competitive adsorption of water on weak binding metals. We further determined that only copper sites capable of adsorbing carbon monoxide in an atop configuration visible to the surface enhanced infrared spectroscopy are active in the electrochemical carbon monoxide reduction reaction.

Infrared and Raman spectroscopies are often assumed to provide similar insights into heterogeneous reaction mechanisms. This study shows that these techniques provide similar data when CO is strongly bound to a surface, yet distinct subpopulations of CO are probed when binding is weaker.

Density Functional Theory Perspective on the Nonlinear Response of Correlated Electrons across Temperature Regimes

We explore a new formalism to study the nonlinear electronic density response based on Kohn-Sham density functional theory (KS-DFT) at partially and strongly quantum degenerate regimes. It is demonstrated that the KS-DFT calculations are able to accurately reproduce the available path integral Monte Carlo simulation results at temperatures relevant for warm dense matter research. The existing analytical results for the quadratic and cubic response functions are rigorously tested. It is demonstrated that the analytical results for the quadratic response function closely agree with the KS-DFT data. Furthermore, the performed analysis reveals that currently available analytical formulas for the cubic response function are not able to describe simulation results, neither qualitatively nor quantitatively, at small wavenumbers q < 2qF, with qF being the Fermi wavenumber. The results show that KS-DFT can be used to describe warm dense matter that is strongly perturbed by an external field with remarkable accuracy. Furthermore, it is demonstrated that KS-DFT constitutes a valuable tool to guide the development of the nonlinear response theory of correlated quantum electrons from ambient to extreme conditions. This opens up new avenues to study nonlinear effects in a gamut of different contexts at conditions that cannot be accessed with previously used path integral Monte Carlo methods.

Low-Pressure Electrochemical Synthesis of Complex High-Pressure Superconducting Superhydrides

There is great current interest in multicomponent superhydrides due to their unique quantum properties under pressure. A remarkable example is the ternary superhydride Li_{2}MgH_{16} computationally identified to have an unprecedented high superconducting critical temperature T_{c} of ∼470  K at 250 GPa. However, the very high synthesis pressures required remains a significant hurdle for detailed study and potential applications. In this Letter, we evaluate the feasibility of synthesizing ternary Li-Mg superhydrides by the recently proposed pressure-potential (P^{2}) method that uniquely combines electrochemistry and applied pressure to control synthesis and stability. The results indicate that it is possible to synthesize Li-Mg superhydrides at modest pressures by applying suitable electrode potentials. Using pressure alone, no Li-Mg ternary hydrides are predicted to be thermodynamically stable, but in the presence of electrode potentials, both Li_{2}MgH_{16} and Li_{4}MgH_{24} can be stabilized at modest pressures. Three polymorphs are predicted as ground states of Li_{2}MgH_{16} below 300 GPa, with transitions at 33 and 160 GPa. The highest pressure phase is superconducting, while the two at lower pressures are not. Our findings point out the potentially important role of the P^{2} method in controlling phase stability of complex multicomponent superhydrides.

Prospects of Improving Molecular Solar Energy Storage of the Norbornadiene/Quadricyclane System through Bridgehead Modifications

We have investigated novel bicyclic diene molecular solar thermal energy storage systems that presently are the ones with the highest predicted energy density. Using a variety of different ab initio quantum chemical methods, we report storage energies, absorption spectra, and reaction barriers for the release of stored energy for a series of bicyclic dienes. The bicyclic dienes are all constructed by modifying the bridgehead of the well-known norbornadiene/quadricyclane (NBD/QC) system. In conclusion, we find it promising that it is possible to significantly amplify the storage energy of the NBD/QC system without seriously compromising other crucial properties by introducing simple modifications to the bridgehead.

NHC-Stabilized Au10 Nanoclusters and Their Conversion to Au25 Nanoclusters

Herein, we describe the synthesis of a toroidal Au10 cluster stabilized by N-heterocyclic carbene and halide ligands via reduction of the corresponding NHC-Au-X complexes (X = Cl, Br, I). The significant effect of the halide ligands on the formation, stability, and further conversions of these clusters is presented. While solutions of the chloride derivatives of Au10 show no change even upon heating, the bromide derivative readily undergoes conversion to form a biicosahedral Au25 cluster at room temperature. For the iodide derivative, the formation of a significant amount of Au25 was observed even upon the reduction of NHC-Au-I. The isolated bromide derivative of the Au25 cluster displays a relatively high (ca. 15%) photoluminescence quantum yield, attributed to the high rigidity of the cluster, which is enforced by multiple CH-π interactions within the molecular structure. Density functional theory computations are used to characterize the electronic structure and optical absorption of the Au10 cluster. 13C-Labeling is employed to assist with characterization of the products and to observe their conversions by NMR spectroscopy.

Shortcut to Self-Consistent Light-Matter Interaction and Realistic Spectra from First Principles

We introduce a simple approach to how an electromagnetic environment can be efficiently embedded into state-of-the-art electronic structure methods, taking the form of radiation-reaction forces. We demonstrate that this self-consistently provides access to radiative emission, natural linewidth, Lamb shifts, strong coupling, electromagnetically induced transparency, Purcell-enhanced and superradiant emission. As an example, we illustrate its seamless integration into time-dependent density-functional theory with virtually no additional cost, presenting a convenient shortcut to light-matter interactions.

Substituent Control of σ-Interference Effects in the Transmission of Saturated Molecules

The single-molecule conductance of saturated molecules can potentially be fully suppressed by destructive quantum interference in their σ-system. However, only few molecules with σ-interference have been identified, and the structure-property relationship remains to be elucidated. Here, we explore the role of substituents in modulating the electronic transmission of saturated molecules. In functionalized bicyclo[2.2.2]octanes, the transmission is suppressed by σ-interference when fluorine substituents are applied. For bicyclo[2.2.2]octasilane and -octagermanes, the transmission is suppressed when carbon-based substituents are used, and such molecules are likely to be highly insulating. For the carbon-based substituents, we find a strong correlation between the appropriate Hammett constants and the transmission. The substituent effect enables systematic optimization of the insulating properties of saturated molecular cores.

Dimerization of dehydrogenated polycyclic aromatic hydrocarbons on graphene

Dimerization of polycyclic aromatic hydrocarbons (PAHs) is an important, yet poorly understood, step in the on-surface synthesis of graphene (nanoribbon), soot formation, and growth of carbonaceous dust grains in the interstellar medium (ISM). The on-surface synthesis of graphene and the growth of carbonaceous dust grains in the ISM require the chemical dimerization in which chemical bonds are formed between PAH monomers. An accurate and cheap method of exploring structure rearrangements is needed to reveal the mechanism of chemical dimerization on surfaces. This work has investigated the chemical dimerization of two dehydrogenated PAHs (coronene and pentacene) on graphene via an evolutionary algorithm augmented by machine learning surrogate potentials and a set of customized structure operators. Different dimer structures on surfaces have been successfully located by our structure search methods. Their binding energies are within the experimental errors of temperature programmed desorption measurements. The mechanism of coronene dimer formation on graphene is further studied and discussed.

Suspension electrolyte with modified Li+ solvation environment for lithium metal batteries

Designing a stable solid-electrolyte interphase on a Li anode is imperative to developing reliable Li metal batteries. Herein, we report a suspension electrolyte design that modifies the Li⁺ solvation environment in liquid electrolytes and creates inorganic-rich solid-electrolyte interphases on Li. Li₂O nanoparticles suspended in liquid electrolytes were investigated as a proof of concept. Through theoretical and empirical analyses of Li₂O suspension electrolytes, the roles played by Li₂O in the liquid electrolyte and solid-electrolyte interphases of the Li anode are elucidated. Also, the suspension electrolyte design is applied in conventional and state-of-the-art high-performance electrolytes to demonstrate its applicability. Based on electrochemical analyses, improved Coulombic efficiency (up to ~99.7%), reduced Li nucleation overpotential, stabilized Li interphases and prolonged cycle life of anode-free cells (~70 cycles at 80% of initial capacity) were achieved with the suspension electrolytes. We expect this design principle and our findings to be expanded into developing electrolytes and solid-electrolyte interphases for Li metal batteries.

High-Quality Graphene Using Boudouard Reaction

Following the game-changing high-pressure CO (HiPco) process that established the first facile route toward large-scale production of single-walled carbon nanotubes, CO synthesis of cm-sized graphene crystals of ultra-high purity grown during tens of minutes is proposed. The Boudouard reaction serves for the first time to produce individual monolayer structures on the surface of a metal catalyst, thereby providing a chemical vapor deposition technique free from molecular and atomic hydrogen as well as vacuum conditions. This approach facilitates inhibition of the graphene nucleation from the CO/CO2 mixture and maintains a high growth rate of graphene seeds reaching large-scale monocrystals. Unique features of the Boudouard reaction coupled with CO-driven catalyst engineering ensure not only suppression of the second layer growth but also provide a simple and reliable technique for surface cleaning. Aside from being a novel carbon source, carbon monoxide ensures peculiar modification of catalyst and in general opens avenues for breakthrough graphene-catalyst composite production.

Addressing Dynamics at Catalytic Heterogeneous Interfaces with DFT-MD: Anomalous Temperature Distributions from Commonly Used Thermostats

Density functional theory-based molecular dynamics (DFT-MD) has been widely used for studying the chemistry of heterogeneous interfacial systems under operational conditions. We report frequently overlooked errors in thermostated or constant-temperature DFT-MD simulations applied to study (electro)catalytic chemistry. Our results demonstrate that commonly used thermostats such as Nosé-Hoover, Berendsen, and simple velocity-rescaling methods fail to provide a reliable temperature description for systems considered. Instead, nonconstant temperatures and large temperature gradients within the different parts of the system are observed. The errors are not a "feature" of any particular code but are present in several ab initio molecular dynamics implementations. This uneven temperature distribution, due to inadequate thermostatting, is well-known in the classical MD community, where it is ascribed to the failure in kinetic energy equipartition among different degrees of freedom in heterogeneous systems (Harvey et al. J. Comput. Chem. 1998, 726-740) and termed the flying ice cube effect. We provide tantamount evidence that interfacial systems are susceptible to substantial flying ice cube effects and demonstrate that the traditional Nosé-Hoover and Berendsen thermostats should be applied with care when simulating, for example, catalytic properties or structures of solvated interfaces and supported clusters. We conclude that the flying ice cube effect in these systems can be conveniently avoided using Langevin dynamics.

Influence of Au, Ag, and Cu Adatoms on Optical Properties of TiO2 (110) Surface: Predictions from RT-TDDFT Calculations

In this paper, real-time time-dependent density-functional theory (RT-TDDFT) calculations are performed to analyze the optical property and charge transitions of a single noble metal atom deposited on rutile TiO2 (110) surface. The model structures are built reflecting the equilibrium positions of deposited adatoms atop the TiO2 surface. The absorption spectra are calculated for all model structures under study. To provide deeper insight into photo-absorption processes, the transition contribution maps are computed for the states of deposited adatoms involved in transitions. Assuming the photon energy is enough to overcome the band gap of TiO2 (∼3 eV), the photogenerated electrons of TiO2 seem to be partly accumulated around deposited Au atoms. In contrast, this is rarely observed for deposited Ag and Cu atoms. Based on our calculations, we have identified the transition state mechanism that is important for the design strategy of future photocatalytic materials.

Hexacene on Cu(110) and Ag(110): Influence of the Substrate on Molecular Orientation and Interfacial Charge Transfer

Hexacene, composed of six linearly fused benzene rings, is an organic semiconductor material with superior electronic properties. The fundamental understanding of the electronic and chemical properties is prerequisite to any possible application in devices. We investigate the orientation and interface properties of highly ordered hexacene monolayers on Ag(110) and Cu(110) with X-ray photoemission spectroscopy (XPS), photoemission orbital tomography (POT), X-ray absorption spectroscopy (XAS), low-energy electron diffraction (LEED), scanning tunneling microscopy (STM), and density functional theory (DFT). We find pronounced differences in the structural arrangement of the molecules and the electronic properties at the metal/organic interfaces for the two substrates. While on Cu(110) the molecules adsorb with their long molecular axis parallel to the high symmetry substrate direction, on Ag(110), hexacene adsorbs in an azimuthally slightly rotated geometry with respect to the metal rows of the substrate. In both cases, molecular planes are oriented parallel to the substrate. A pronounced charge transfer from both substrates to different molecular states affects the effective charge of different C atoms of the molecule. Through analysis of experimental and theoretical data, we found out that on Ag(110) the LUMO of the molecule is occupied through charge transfer from the metal, whereas on Cu(110) even the LUMO+1 receives a charge. Interface dipoles are determined to a large extent by the push-back effect, which are also found to differ significantly between 6A/Ag(110) and 6A/Cu(110).