Search for: All records

Understanding, predicting, and ultimately controlling exciton band structure and exciton dynamics are central to diverse chemical and materials problems. Here, we have developed a first-principles method to determine exciton dispersion and exciton–phonon interaction in semiconducting and insulating solids based on time-dependent density functional theory. The first-principles method is formulated in planewave bases and pseudopotentials and can be used to compute exciton band structures, exciton charge density, ionic forces, the non-adiabatic coupling matrix between excitonic states, and the exciton–phonon coupling matrix. Based on the spinor formulation, the method enables self-consistent noncollinear calculations to capture spin-orbital coupling. Hybrid exchange-correlation functionals are incorporated to deal with long-range electron–hole interactions in solids. A sub-Hilbert space approximation is introduced to reduce the computational cost without loss of accuracy. For validations, we have applied the method to compute the exciton band structure and exciton–phonon coupling strength in transition metal dichalcogenide monolayers; both agree very well with the previous GW-Bethe–Salpeter equation and experimental results. This development paves the way for accurate determinations of exciton dynamics in a wide range of solid-state materials.

 

During geomagnetic storms a large amount of energy is transferred into the ionosphere-thermosphere (IT) system, leading to local and global changes in e.g., the dynamics, composition, and neutral density. The more steady energy from the lower atmosphere into the IT system is in general much smaller than the energy input from the magnetosphere, especially during geomagnetic storms, and therefore details of the lower atmosphere forcing are often neglected in storm time simulations. In this study we compare the neutral density observed by Swarm-C during the moderate geomagnetic storm of 31 January to 3 February 2016 with the Thermosphere-Ionosphere-Electrodynamics-GCM (TIEGCM) finding that the model can capture the observed large scale neutral density variations better in the southern than northern hemisphere. The importance of more realistic lower atmospheric (LB) variations as specified by the Whole Atmosphere Community Climate Model eXtended (WACCM-X) with specified dynamics (SD) is demonstrated by improving especially the northern hemisphere neutral density by up to 15% compared to using climatological LB forcing. Further analysis highlights the importance of the background atmospheric condition in facilitating hemispheric different neutral density changes in response to the LB perturbations. In comparison, employing observationally based field-aligned current (FAC) versus using an empirical model to describe magnetosphere-ionosphere (MI) coupling leads to an 7–20% improved northern hemisphere neutral density. The results highlight the importance of the lower atmospheric variations and high latitude forcing in simulating the absolute large scale neutral density especially the hemispheric differences. However, focusing on the storm time variation with respect to the quiescent time, the lower atmospheric influence is reduced to 1–1.5% improvement with respect to the total observed neutral density. The results provide some guidance on the importance of more realistic upper boundary forcing and lower atmospheric variations when modeling large scale, absolute and relative neutral density variations. 

Electrochemical conversion of CO 2 into value-added chemicals continues to draw interest in renewable energy applications. Although many metal catalysts are active in the CO 2 reduction reaction (CO 2 RR), their reactivity and selectivity are nonetheless hindered by the competing hydrogen evolution reaction (HER). The competition of the HER and CO 2 RR stems from the energy scaling relationship between their reaction intermediates. Herein, we predict that bimetallic monolayer electrocatalysts (BMEs) – a monolayer of transition metals on top of extended metal substrates – could produce dual-functional active sites that circumvent the scaling relationship between the adsorption energies of HER and CO 2 RR intermediates. The antibonding interaction between the adsorbed H and the metal substrate is revealed to be responsible for circumventing the scaling relationship. Based on extensive density functional theory (DFT) calculations, we identify 11 BMEs which are highly active and selective toward the formation of formic acid with a much suppressed HER. The H–substrate antibonding interaction also leads to superior CO 2 RR performance on monolayer-coated penta-twinned nanowires. 

Abstract Recent experiments observed significant energy transfer in type-II van der Waals (vdW) heterostructures, such as WS 2 /MoSe 2 , which is surprising due to their staggered band alignment and weak spectral overlap. In this work, we carry out first-principles calculations to shed light on energy and charge transfer in WS 2 /MoSe 2 heterostructure. Incorporating excitonic effect in nonadiabatic electronic dynamics, our first-principles calculations uncover a two-step process in competing energy and charge transfer, unravel their relative efficiencies and explore the means to control their competition. While both Dexter and Förster mechanisms can be responsible for energy transfer, they are shown to operate at different conditions. The excitonic effect is revealed to drive ultrafast energy and charge transfer in type-II WS 2 /MoSe 2 heterostructure. Our work provides a comprehensive picture of exciton dynamics in vdW heterostructures and paves the way for rational design of novel vdW heterostructures for optoelectronic and photovoltaic applications. 

Excitons can be trapped by moiré potentials in van der Waals (vdW) heterostructures, forming ordered arrays of quantum dots. Excitons can also be trapped by defect potentials as single photon emitters. While the moiré and defect potentials in vdW heterostructures have been studied separately, their interplay remains largely unexplored. Here, we perform first-principles calculations to elucidate the interplay of the two potentials in determining the optoelectronic properties of twisted MoS 2 /WS 2 heterobilayers. The binding energy, charge density, localization, and hybridization of the moiré excitons can be modulated by the competition and cooperation of the two potentials. Their interplay can also be tuned by vertical electric fields, which can either de-trap the excitons or strongly localize them. One can further tailor the interplay of the two potentials via defect engineering to create one-dimensional exciton lattices with tunable orientations. Our work establishes defect engineering as a promising strategy to realize on-demand optoelectronic responses.