Shear Modes in a 2D Polar Metal.

Low-frequency shear and breathing modes are important Raman signatures of two-dimensional (2D) materials, providing information on the number of layers and insights into interlayer interactions. We elucidate the nature of low-frequency modes in a 2D polar metal-2D Ga covalently bonded to a SiC substrate, using a first-principles Green's function-based approach. The low-frequency Raman modes are dominated by surface resonance modes, consisting primarily of out-of-phase shear modes in Ga, coupled to SiC phonons. Breathing modes are strongly coupled to the substrate and do not give rise to peaks in the phonon spectra. The highest-frequency shear mode blue-shifts significantly with increasing thickness, reflecting both an increase in the number of Ga layers and an increase in the effective interlayer force constant. The surface resonance modes evolve into localized 2D Ga modes as the phonon momentum increases. The predicted low-frequency modes are consistent with Raman measurements on 2D polar Ga.

Organic-2D Material Heterostructures: A Promising Platform for Exciton Condensation and Multiplication.

We predict that high temperature Bose-Einstein condensation of charge transfer excitons can be achieved in organic-two-dimensional (2D) material heterostructures, at ∼50-100 K. Unlike 2D bilayers that can be angle-misaligned, organic-2D systems generally have momentum-direct low-energy excitons, thus favoring condensation. Our predictions are obtained for ZnPc-MoS2 using state-of-the-art first-principles calculations with the GW-BSE approach. The exciton energies we predict are in excellent agreement with recent experiments. The lowest energy charge transfer excitons in ZnPc-MoS2 are strongly bound with a spatial extent of ∼1-2 nm and long lifetimes (τ0 ∼ 1 ns), making them ideal for exciton condensation. We also predict the emergence of inter-ZnPc excitons that are stabilized by the interaction of the molecules with the 2D substrate. This novel way of stabilizing intermolecular excitons by indirect substrate mediation suggests design strategies for singlet fission and exciton multiplication, which are important to overcome the Shockley-Queisser efficiency limit in solar cells.

Chiral Phonons and Giant Magneto-Optical Effect in CrBr3 2D Magnet.

Phonons with chirality determine the optical helicity of inelastic light scattering processes due to their nonzero angular momentum. Here it is shown that 2D magnetic CrBr3 hosts chiral phonons at the Brillouin-zone center. These chiral phonons are linear combinations of the doubly-degenerate Eg phonons, and the phonon eigenmodes exhibit clockwise and counterclockwise rotational vibrations corresponding to angular momenta of l = ± 1. Such Eg chiral phonons completely switch the polarization of incident circularly polarized light. On the other hand, the non-degenerate non-chiral Ag phonons display a giant magneto-optical effect under an external out-of-plane magnetic field, rotating the plane of polarization of the scattered linearly polarized light. The corresponding degree of polarization of the scattered light changes from 91% to -68% as the magnetic field strength increases from 0 to 5 T. In contrast, the chiral Eg modes display no field dependence. The results lay a foundation for the study of phonon chirality and magneto-optical phenomena in 2D magnetic materials, as well as their related applications, such as the phonon Hall effect, topological photonics, and Raman lasing.

Adsorption of methane on single metal atoms supported on graphene: Role of electron back-donation in binding and activation.

We consider single metal atoms supported on graphene as possible candidate systems for on-board vehicular storage of methane or for methane activation. We use density functional theory to study the adsorption of one and two molecules of methane on such graphene-supported single atoms, where the metal atom M is a 3d-transition metal (Sc to Zn). Our results suggest that M = Sc, Ti, and V are the best candidates for gas storage applications, while Ni and Co seem particularly promising with respect to activation of the C-H bond in methane. We find a strong and linear correlation between the adsorption energy of methane and the degree of back-donation of electrons from occupied metal d-states to antibonding methane states. A similar correlation is found between the elongation of C-H bonds and electron back-donation. An important role is played by the graphene substrate in enhancing the binding of methane on metal atoms, compared to the negligible binding observed on isolated metal atoms.

Understanding the electrochemical double layer at the hematite/water interface: A first principles molecular dynamics study.

Using first principles molecular dynamics simulations, we probe the electrochemical double layer formed at the interface between the hematite surface and water. We consider two terminations of the (001) surface, viz., the fully hydroxylated (OH) and the stoichiometric (FeO3Fe) termination. We explicitly incorporate the counterions (Na+ and F-) in the solution, and model both specific and nonspecific adsorption of F- ions. We find that F- ions prefer to bind directly to the Fe ions (specific adsorption), with a substantial energy gain (0.75 eV/ion). We investigate the effect of the interface and the counterions on the dipole of individual water molecules. We find significant deviations of +0.2/-0.15 D for dipoles of the first solvation shell water molecules of F-/Na+ ions, respectively. Additionally, the hydration layers at the interface show an enhancement in the dipole moment resulting from stronger hydrogen bonding interactions between the water molecules and surface charged species. Furthermore, we analyze the electrostatic potential profile at the solid/liquid interface as a function of the kind of counterion present in the double layer and compute the capacitance of the compact (Helmholtz) layer. We find that our results (40.3 ± 3.5 µF/cm2 for the OH termination and 51 ± 5 µF/cm2 for the FeO3Fe termination) compare favorably with values reported by potentiometric titration based experimental studies (10-100 µF/cm2).

Characterization of peroxo reaction intermediates in the water oxidation process on hematite surfaces.

We use density functional theory-based calculations to study structural, electronic, and magnetic properties of two key reaction intermediates on a hematite, [Formula: see text]-Fe2O3, photoanode during the solar-driven water splitting reaction. Both intermediates contain an oxygen atom bonded to a surface iron atom. In one case, the adsorbed oxygen also forms a peroxo bond with a lattice oxygen from hematite; in the second case no such bond is formed. Both configurations are energetically equivalent and are related to the overpotential-determining step in the oxygen evolution reaction. The calculated reaction path for the breaking of the peroxo bond shows a barrier of about 0.86 eV for the transformation between the two intermediates. We explain this high barrier with the drastically different electronic and magnetic structure, which we also analyze using maximally localized Wannier functions. Photo-generated electron holes are shown to localize preferentially close to the reaction center at the surface in both configurations. In the case of the oxo species, this localization favors subsequent electron transfer steps during the oxygen evolution cycle. In the case of the peroxo configuration, this fact together with the high barrier for breaking the oxygen-oxygen bond indicates a possible loss mechanism due to hole trapping. Graphical Abstract Calculated spin density at a hematite surface with peroxo intermediate.

Quantum mechanical effects in zwitterionic amino acids: The case of proline, hydroxyproline, and alanine in water.

In this work, we use ab initio molecular dynamics simulations to elucidate the electronic properties of three hydrated zwitterionic amino acids, namely proline, hydroxyproline, and alanine, the former two forming an important constituent of collagen. In all three systems, we find a substantial amount of charge transfer between the amino acids and surrounding solvent, which, rather surprisingly, also involves the reorganization of electron density near the hydrophobic non-polar groups. Water around proline appears to be slightly more polarized, as reflected by the enhanced water dipole moment in its hydration shell. This observation is also complemented by an examination of the IR spectra of the three systems where there is a subtle red and blue shift in the O-H stretch and bend regions, respectively, for proline. We show that polarizability of these amino acids as revealed by a dipole moment analysis involves a significant enhancement from the solvent and that this also involves non-polar groups. Our results suggest that quantum mechanical effects are likely to be important in understanding the coupling between biomolecules and water in general and in hydrophobic interactions.

Magnetoelectric ϵ-Fe2O3: DFT study of a potential candidate for electrode material in photoelectrochemical cells.

The metastable iron oxide ϵ-Fe2O3 is rare but known for its magnetoelectric properties. While the more common alpha phase has been recognized for a long time as a suitable material for photoelectrochemical cells, its use is limited because of the electron-hole recombination problem when exposed to light. The indirect bandgap of the epsilon phase with its spontaneous polarization may offer a better potential for the application in photoelectrochemistry. Here, we report a detailed study of the electronic and structural features of the epsilon phase of iron oxide, its stability in thin films, and possible water dissociation reactions. Our studies are performed using density functional theory with a Hubbard-U correction. We observe that the stable ϵ-Fe2O3 surfaces favor the dissociation of water. The average difference in the energies of the states when water is adsorbed and when it is dissociated is roughly found to be -0.40 eV. Our results compare with the available experimental results where the epsilon phase is reported to be more efficient for the release of hydrogen from renewable oxygenates when exposed to sunlight.

Passivation of surface states of α-Fe2O3(0001) surface by deposition of Ga2O3 overlayers: A density functional theory study.

There is a big debate in the community regarding the role of surface states of hematite in the photoelectrochemical water splitting. Experimental studies on non-catalytic overlayers passivating the hematite surface states claim a favorable reduction in the overpotential for the water splitting reaction. As a first step towards understanding the effect of these overlayers, we have studied the system Ga2O3 overlayers on hematite (0001) surfaces using first principles computations in the PBE+U framework. Our computations suggest that stoichiometric terminations of Ga2O3 overlayers are energetically more favored than the bare surface, at ambient oxygen chemical potentials. Energetics suggest that the overlayers prefer to grow via a layer-plus-island (Stranski-Krastanov) growth mode with a critical layer thickness of 1-2 layers. Thus, a complete wetting of the hematite surface by an overlayer of gallium oxide is thermodynamically favored. We establish that the effect of deposition of the Ga2O3 overlayers on the bare hematite surface is to passivate the surface states for the stoichiometric termination. For the oxygen terminated surface which is the most stable termination under photoelectrochemical conditions, the effect of deposition of the Ga2O3 overlayer is to passivate the hole-trapping surface state.

Graphene oxide as an optimal candidate material for methane storage.

Methane, the primary constituent of natural gas, binds too weakly to nanostructured carbons to meet the targets set for on-board vehicular storage to be viable. We show, using density functional theory calculations, that replacing graphene by graphene oxide increases the adsorption energy of methane by 50%. This enhancement is sufficient to achieve the optimal binding strength. In order to gain insight into the sources of this increased binding, that could also be used to formulate design principles for novel storage materials, we consider a sequence of model systems that progressively take us from graphene to graphene oxide. A careful analysis of the various contributions to the weak binding between the methane molecule and the graphene oxide shows that the enhancement has important contributions from London dispersion interactions as well as electrostatic interactions such as Debye interactions, aided by geometric curvature induced primarily by the presence of epoxy groups.

Physical origins of weak H2 binding on carbon nanostructures: insight from ab initio studies of chemically functionalized graphene nanoribbons.

We have performed ab initio density functional theory calculations, incorporating London dispersion corrections, to study the absorption of molecular hydrogen on zigzag graphene nanoribbons whose edges have been functionalized by OH, NH2, COOH, NO2, or H2PO3. We find that hydrogen molecules always preferentially bind at or near the functionalized edge, and display induced dipole moments. Binding is generally enhanced by the presence of polar functional groups. The largest gains are observed for groups with oxygen lone pairs that can facilitate local charge reorganization, with the biggest single enhancement in adsorption energy found for "strong functionalization" by H2PO3 (115 meV/H2 versus 52 meV/H2 on bare graphene). We show that for binding on the "outer edge" near the functional group, the presence of the group can introduce appreciable contributions from Debye interactions and higher-order multipole electrostatic terms, in addition to the dominant London dispersion interactions. For those functional groups that contain the OH moiety, the adsorption energy is linearly proportional to the number of lone pairs on oxygen atoms. Mixed functionalization with two different functional groups on a graphene edge can also have a synergistic effect, particularly when electron-donating and electron-withdrawing groups are combined. For binding on the "inner edge" somewhat farther from the functional group, most of the binding again arises from London interactions; however, there is also significant charge redistribution in the π manifold, which directly reflects the electron donating or withdrawing capacity of the functional group. Our results offer insight into the specific origins of weak binding of gas molecules on graphene, and suggest that edge functionalization could perhaps be used in combination with other strategies to increase the uptake of hydrogen in graphene. They also have relevance for the storage of hydrogen in porous carbon materials, such as activated carbons.

Tuning spin transport properties and molecular magnetoresistance through contact geometry.

Molecular spintronics seeks to unite the advantages of using organic molecules as nanoelectronic components, with the benefits of using spin as an additional degree of freedom. For technological applications, an important quantity is the molecular magnetoresistance. In this work, we show that this parameter is very sensitive to the contact geometry. To demonstrate this, we perform ab initio calculations, combining the non-equilibrium Green's function method with density functional theory, on a dithienylethene molecule placed between spin-polarized nickel leads of varying geometries. We find that, in general, the magnetoresistance is significantly higher when the contact is made to sharp tips than to flat surfaces. Interestingly, this holds true for both resonant and tunneling conduction regimes, i.e., when the molecule is in its "closed" and "open" conformations, respectively. We find that changing the lead geometry can increase the magnetoresistance by up to a factor of ∼5. We also introduce a simple model that, despite requiring minimal computational time, can recapture our ab initio results for the behavior of magnetoresistance as a function of bias voltage. This model requires as its input only the density of states on the anchoring atoms, at zero bias voltage. We also find that the non-resonant conductance in the open conformation of the molecule is significantly impacted by the lead geometry. As a result, the ratio of the current in the closed and open conformations can also be tuned by varying the geometry of the leads, and increased by ∼400%.