Extremely anisotropic van der Waals thermal conductors.

The densification of integrated circuits requires thermal management strategies and high thermal conductivity materials1-3. Recent innovations include the development of materials with thermal conduction anisotropy, which can remove hotspots along the fast-axis direction and provide thermal insulation along the slow axis4,5. However, most artificially engineered thermal conductors have anisotropy ratios much smaller than those seen in naturally anisotropic materials. Here we report extremely anisotropic thermal conductors based on large-area van der Waals thin films with random interlayer rotations, which produce a room-temperature thermal anisotropy ratio close to 900 in MoS2, one of the highest ever reported. This is enabled by the interlayer rotations that impede the through-plane thermal transport, while the long-range intralayer crystallinity maintains high in-plane thermal conductivity. We measure ultralow thermal conductivities in the through-plane direction for MoS2 (57 ± 3 mW m-1 K-1) and WS2 (41 ± 3 mW m-1 K-1) films, and we quantitatively explain these values using molecular dynamics simulations that reveal one-dimensional glass-like thermal transport. Conversely, the in-plane thermal conductivity in these MoS2 films is close to the single-crystal value. Covering nanofabricated gold electrodes with our anisotropic films prevents overheating of the electrodes and blocks heat from reaching the device surface. Our work establishes interlayer rotation in crystalline layered materials as a new degree of freedom for engineering-directed heat transport in solid-state systems.

Machine Learning for Polaritonic Chemistry: Accessing Chemical Kinetics.

Altering chemical reactivity and material structure in confined optical environments is on the rise, and yet, a conclusive understanding of the microscopic mechanisms remains elusive. This originates mostly from the fact that accurately predicting vibrational and reactive dynamics for soluted ensembles of realistic molecules is no small endeavor, and adding (collective) strong light-matter interaction does not simplify matters. Here, we establish a framework based on a combination of machine learning (ML) models, trained using density-functional theory calculations and molecular dynamics to accelerate such simulations. We then apply this approach to evaluate strong coupling, changes in reaction rate constant, and their influence on enthalpy and entropy for the deprotection reaction of 1-phenyl-2-trimethylsilylacetylene, which has been studied previously both experimentally and using ab initio simulations. While we find qualitative agreement with critical experimental observations, especially with regard to the changes in kinetics, we also find differences in comparison with previous theoretical predictions. The features for which the ML-accelerated and ab initio simulations agree show the experimentally estimated kinetic behavior. Conflicting features indicate that a contribution of dynamic electronic polarization to the reaction process is more relevant than currently believed. Our work demonstrates the practical use of ML for polaritonic chemistry, discusses limitations of common approximations, and paves the way for a more holistic description of polaritonic chemistry.

GPAW: An open Python package for electronic structure calculations.

We review the GPAW open-source Python package for electronic structure calculations. GPAW is based on the projector-augmented wave method and can solve the self-consistent density functional theory (DFT) equations using three different wave-function representations, namely real-space grids, plane waves, and numerical atomic orbitals. The three representations are complementary and mutually independent and can be connected by transformations via the real-space grid. This multi-basis feature renders GPAW highly versatile and unique among similar codes. By virtue of its modular structure, the GPAW code constitutes an ideal platform for the implementation of new features and methodologies. Moreover, it is well integrated with the Atomic Simulation Environment (ASE), providing a flexible and dynamic user interface. In addition to ground-state DFT calculations, GPAW supports many-body GW band structures, optical excitations from the Bethe-Salpeter Equation, variational calculations of excited states in molecules and solids via direct optimization, and real-time propagation of the Kohn-Sham equations within time-dependent DFT. A range of more advanced methods to describe magnetic excitations and non-collinear magnetism in solids are also now available. In addition, GPAW can calculate non-linear optical tensors of solids, charged crystal point defects, and much more. Recently, support for graphics processing unit (GPU) acceleration has been achieved with minor modifications to the GPAW code thanks to the CuPy library. We end the review with an outlook, describing some future plans for GPAW.

Hot-Carrier Transfer across a Nanoparticle-Molecule Junction: The Importance of Orbital Hybridization and Level Alignment.

While direct hot-carrier transfer can increase photocatalytic activity, it is difficult to discern experimentally and competes with several other mechanisms. To shed light on these aspects, here, we model from first-principles hot-carrier generation across the interface between plasmonic nanoparticles and a CO molecule. The hot-electron transfer probability depends nonmonotonically on the nanoparticle-molecule distance and can be effective at long distances, even before a strong chemical bond can form; hot-hole transfer on the other hand is limited to shorter distances. These observations can be explained by the energetic alignment between molecular and nanoparticle states as well as the excitation frequency. The hybridization of the molecular orbitals is the key predictor for hot-carrier transfer in these systems, emphasizing the necessity of ground state hybridization for accurate predictions. Finally, we show a nontrivial dependence of the hot-carrier distribution on the excitation energy, which could be exploited when optimizing photocatalytic systems.

GPUMD: A package for constructing accurate machine-learned potentials and performing highly efficient atomistic simulations.

We present our latest advancements of machine-learned potentials (MLPs) based on the neuroevolution potential (NEP) framework introduced in Fan et al. [Phys. Rev. B 104, 104309 (2021)] and their implementation in the open-source package gpumd. We increase the accuracy of NEP models both by improving the radial functions in the atomic-environment descriptor using a linear combination of Chebyshev basis functions and by extending the angular descriptor with some four-body and five-body contributions as in the atomic cluster expansion approach. We also detail our efficient implementation of the NEP approach in graphics processing units as well as our workflow for the construction of NEP models and demonstrate their application in large-scale atomistic simulations. By comparing to state-of-the-art MLPs, we show that the NEP approach not only achieves above-average accuracy but also is far more computationally efficient. These results demonstrate that the gpumd package is a promising tool for solving challenging problems requiring highly accurate, large-scale atomistic simulations. To enable the construction of MLPs using a minimal training set, we propose an active-learning scheme based on the latent space of a pre-trained NEP model. Finally, we introduce three separate Python packages, viz., gpyumd, calorine, and pynep, that enable the integration of gpumd into Python workflows.

Investigating the Chemical Ordering in Quaternary Clathrate Ba8AlxGa16-xGe30.

Recently, there has been an increased interest in quaternary clathrate systems as promising thermoelectric materials. Because of their increased complexity, however, the chemical ordering in the host framework of quaternary clathrates has not yet been comprehensively analyzed. Here, we have synthesized a prototypical quaternary type-I clathrate Ba8AlxGa16-xGe30 by Czochralski and flux methods, and we employed a combination of X-ray and neutron diffraction along with atomic scale simulations to investigate chemical ordering in this material. We show that the site occupancy factors of trivalent elements at the 6c site differ, depending on the synthesis method, which can be attributed to the level of equilibration. The flux-grown samples are consistent with the simulated high-temperature disordered configuration, while the degree of ordering for the Czochralski sample lies between the ground state and the high-temperature state. Moreover, we demonstrate that the atomic displacement parameters of the Ba atoms in the larger tetrakaidecahedral cages are related to chemical ordering. Specifically, Ba atoms are either displaced toward the periphery or localized at the cage centers. Consequently, this study reveals key relationships between the chemical ordering in the quaternary clathrates Ba8AlxGa16-xGe30 and the structural properties, thereby offering new perspectives on designing these materials and optimizing their thermoelectric properties.

Dipolar coupling of nanoparticle-molecule assemblies: An efficient approach for studying strong coupling.

Strong light-matter interactions facilitate not only emerging applications in quantum and non-linear optics but also modifications of properties of materials. In particular, the latter possibility has spurred the development of advanced theoretical techniques that can accurately capture both quantum optical and quantum chemical degrees of freedom. These methods are, however, computationally very demanding, which limits their application range. Here, we demonstrate that the optical spectra of nanoparticle-molecule assemblies, including strong coupling effects, can be predicted with good accuracy using a subsystem approach, in which the response functions of different units are coupled only at the dipolar level. We demonstrate this approach by comparison with previous time-dependent density functional theory calculations for fully coupled systems of Al nanoparticles and benzene molecules. While the present study only considers few-particle systems, the approach can be readily extended to much larger systems and to include explicit optical-cavity modes.

Tunable Phases of Moiré Excitons in van der Waals Heterostructures.

Stacking monolayers of transition metal dichalcogenides into a heterostructure with a finite twist-angle gives rise to artificial moiré superlattices with a tunable periodicity. As a consequence, excitons experience a periodic potential, which can be exploited to tailor optoelectronic properties of these materials. Whereas recent experimental studies have confirmed twist-angle-dependent optical spectra, the microscopic origin of moiré exciton resonances has not been fully clarified yet. Here, we combine first-principles calculations with the excitonic density matrix formalism to study transitions between different moiré exciton phases and their impact on optical properties of the twisted MoSe2/WSe2 heterostructure. At angles smaller than 2°, we find flat, moiré-trapped states for inter- and intralayer excitons. This moiré exciton phase changes into completely delocalized states at 3°. We predict a linear and quadratic twist-angle dependence of excitonic resonances for the moiré-trapped and delocalized exciton phases, respectively.

To Every Rule There is an Exception: A Rational Extension of Loewenstein's Rule.

Loewenstein's rule, which states that Al-O-Al motifs are energetically unstable, is fundamental to the understanding and design of zeolites. Here, using a combination of electronic structure calculations and lattice models, we show under which circumstances this rule becomes invalid and how it can be rationally extended using the chabasite framework for demonstration.

Non-Bonded Radii of the Atoms Under Compression.

We present quantum mechanical estimates for non-bonded, van der Waals-like, radii of 93 atoms in a pressure range from 0 to 300 gigapascal. Trends in radii are largely maintained under pressure, but atoms also change place in their relative size ordering. Multiple isobaric contractions of radii are predicted and are explained by pressure-induced changes to the electronic ground state configurations of the atoms. The presented radii are predictive of drastically different chemistry under high pressure and permit an extension of chemical thinking to different thermodynamic regimes. For example, they can aid in assignment of bonded and non-bonded contacts, for distinguishing molecular entities, and for estimating available space inside compressed materials. All data has been made available in an interactive web application.

Optimizing photon upconversion by decoupling excimer formation and triplet triplet annihilation.

Perylene is a promising annihilator candidate for triplet-triplet annihilation photon upconversion, which has been successfully used in solar cells and in photocatalysis. Perylene can, however, form excimers, reducing the energy conversion efficiency and hindering further development of TTA-UC systems. Alkyl substitution of perylene can suppress excimer formation, but decelerate triplet energy transfer and triplet-triplet annihilation at the same time. Our results show that mono-substitution with small alkyl groups selectively blocks excimer formation without severly compromising the TTA-UC efficiency. The experimental results are complemented by DFT calculations, which demonstrate that excimer formation is suppressed by steric repulsion. The results demonstrate how the chemical structure can be modified to block unwanted intermolecular excited state relaxation pathways with minimal effect on the preferred ones.

Direct hot-carrier transfer in plasmonic catalysis.

Plasmonic metal nanoparticles can concentrate optical energy and enhance chemical reactions on their surfaces. Plasmons can interact with adsorbate orbitals and decay by directly exciting a carrier from the metal to the adsorbate in a process termed the direct-transfer process. Although this process could be useful for enhancing the efficiency of a chemical reaction, it remains poorly understood. Here, we report a preliminary investigation employing time-dependent density-functional theory (TDDFT) calculations to capture this process at a model metal-adsorbate interface formed by a silver nanoparticle (Ag147) and a carbon monoxide molecule (CO). Direct hot-electron transfer is observed to occur from the occupied states of Ag to the unoccupied molecular orbitals of CO. We determine the probability of this process and show that it depends on the adsorption site of CO. Our results are expected to aid the design of more efficient metal-molecule interfaces for plasmonic catalysis.

Dithiafulvene derivatized donor-acceptor norbornadienes with redshifted absorption.

Photoisomerization of norbornadiene (N) to its metastable isomer quadricyclane (Q) has attracted interest as a strategy for harvesting and storing solar energy. For this strategy to mature the absorption maximum of N has to be moved from the UV to the visible region. Here we show that functionalization of the system with dithiafulvene (DTF) electron donors causes remarkable redshifts of various N derivatives. Thus, some derivatives were found to absorb light with an absorption onset up to 556 nm. The incorporation of DTF units comes, however, with a drawback with regard to achieving reversible N-to-Q and Q-to-N isomerizations. For some derivatives, the photoisomerization was completely quenched. The compounds were subjected to a computational study to shed light on the underlying reason for this reluctance to undergo photoisomerization. The computational study revealed that in these systems, the first excited state (S1) is positioned close to or lower than the transition state for photoconversion, effectively blocking a possible conversion to Q, thus revealing a practical challenge for the future design of N-Q energy storage systems with an improved solar spectrum match.

Norbornadiene-Based Photoswitches with Exceptional Combination of Solar Spectrum Match and Long-Term Energy Storage.

Norbornadiene-quadricyclane (NBD-QC) photoswitches are candidates for applications in solar thermal energy storage. Functionally, they rely on an intramolecular [2+2] cycloaddition reaction, which couples the S0 landscape on the NBD side to the S1 landscape on the QC side of the reaction and vice-versa. This commonly results in an unfavourable correlation between the first absorption maximum and the barrier for thermal back-conversion. This work demonstrates that this correlation can be counteracted by using steric repulsion to hamper the rotational motion of the side groups along the back-conversion path. It is shown that this modification reduces the correlation between the effective back-conversion barrier and the first absorption maximum and also increases the back-conversion entropy. The resulting molecules exhibit exceptionally long half-lives for their metastable forms without significantly affecting other properties, most notably solar spectrum match and storage density.

Beyond Magic Numbers: Atomic Scale Equilibrium Nanoparticle Shapes for Any Size.

In the pursuit of complete control over morphology in nanoparticle synthesis, knowledge of the thermodynamic equilibrium shapes is a key ingredient. While approaches exist to determine the equilibrium shape in the large size limit (≳10-20 nm) as well as for very small particles (≲2 nm), the experimentally increasingly important intermediate size regime has largely remained elusive. Here, we present an algorithm, based on atomistic simulations in a constrained thermodynamic ensemble, that efficiently predicts equilibrium shapes for any number of atoms in the range from a few tens to many thousands of atoms. We apply the algorithm to Cu, Ag, Au, and Pd particles with diameters between approximately 1 and 7 nm and reveal an energy landscape that is more intricate than previously suggested. The thus obtained particle type distributions demonstrate that the transition from icosahedral particles to decahedral and further into truncated octahedral particles occurs only very gradually, which has implications for the interpretation of experimental data. The approach presented here is extensible to alloys and can in principle also be adapted to represent different chemical environments.

Loss channels in triplet-triplet annihilation photon upconversion: importance of annihilator singlet and triplet surface shapes.

Triplet-triplet annihilation photon upconversion (TTA-UC) can, through a number of energy transfer processes, efficiently combine two low frequency photons into one photon of higher frequency. TTA-UC systems consist of one absorbing species (the sensitizer) and one emitting species (the annihilator). Herein, we show that the structurally similar annihilators, 9,10-diphenylanthracene (DPA, 1), 9-(4-phenylethynyl)-10-phenylanthracene (2) and 9,10-bis(phenylethynyl)anthracene (BPEA, 3) have very different upconversion efficiencies, 15.2 ± 2.8%, 15.9 ± 1.3% and 1.6 ± 0.8%, respectively (of a maximum of 50%). We show that these results can be understood in terms of a loss channel, previously unaccounted for, originating from the difference between the BPEA singlet and triplet surface shapes. The difference between the two surfaces results in a fraction of the triplet state population having geometries not energetically capable of forming the first singlet excited state. This is supported by TD-DFT calculations of the annihilator excited state surfaces as a function of phenyl group rotation. We thereby highlight that the commonly used "spin-statistical factor" should be used with caution when explaining TTA-efficiencies. Furthermore, we show that the precious metal free zinc octaethylporphyrin (ZnOEP) can be used for efficient sensitization and that the upconversion quantum yield is maximized when sensitizer-annihilator spectral overlap is minimized (ZnOEP with 2).

Low Molecular Weight Norbornadiene Derivatives for Molecular Solar-Thermal Energy Storage.

Molecular solar-thermal energy storage systems are based on molecular switches that reversibly convert solar energy into chemical energy. Herein, we report the synthesis, characterization, and computational evaluation of a series of low molecular weight (193-260 g mol(-1) ) norbornadiene-quadricyclane systems. The molecules feature cyano acceptor and ethynyl-substituted aromatic donor groups, leading to a good match with solar irradiation, quantitative photo-thermal conversion between the norbornadiene and quadricyclane, as well as high energy storage densities (396-629 kJ kg(-1) ). The spectroscopic properties and energy storage capability have been further evaluated through density functional theory calculations, which indicate that the ethynyl moiety plays a critical role in obtaining the high oscillator strengths seen for these molecules.

Quantitative Predictions of the Thermal Conductivity in Transition Metal Dichalcogenides: Impact of Point Defects in MoS2 and WS2 Monolayers.

Transition metal dichalcogenides are investigated for various applications at the nanoscale because of their unique combination of properties and dimensionality. For many of the anticipated applications, heat conduction plays an important role. At the same time, these materials often contain relatively large amounts of point defects. Here, we provide a systematic analysis of the impact of intrinsic and selected extrinsic defects on the lattice thermal conductivity of MoS2 and WS2 monolayers. We combine Boltzmann transport theory and Green's function-based T-matrix approach for the calculation of scattering rates. The force constants for the defect configurations are obtained from density functional theory calculations via a regression approach, which allows us to sample a rather large number of defects at a moderate computational cost and to systematically enforce both the translational and rotational acoustic sum rules. The calculated lattice thermal conductivity is in quantitative agreement with the experimental data for heat transport and defect concentrations for both MoS2 and WS2. Crucially, this demonstrates that the strong deviation from a 1/T temperature dependence of the lattice thermal conductivity observed experimentally can be fully explained by the presence of point defects. We furthermore predict the scattering strengths of the intrinsic defects to decrease in the sequence VMo ≈ V2S= > V2S⊥ > VS > Sad in both materials, while the scattering rates for the extrinsic (adatom) defects decrease with increasing mass such that Liad > Naad > Kad. Compared with earlier work, we find that both intrinsic and extrinsic adatoms are relatively weak scatterers. We attribute this difference to the treatment of the translational and rotational acoustic sum rules, which, if not enforced, can lead to spurious contributions in the zero-frequency limit.