Asymptotic behavior and interpretation of virtual states: The effects of confinement and of basis sets

A real-space high order finite difference method is used to analyze the effect of spherical domain size on the Hartree-Fock (and density functional theory) virtual eigenstates. We show the domain size dependence of both positive and negative virtual eigenvalues of the Hartree-Fock equations for small molecules. We demonstrate that positive states behave like a particle in spherical well and show how they approach zero. For the negative eigenstates, we show that large domains are needed to get the correct eigenvalues. We compare our results to those of Gaussian basis sets and draw some conclusions for real-space, basis-sets, and plane-waves calculations. (C) 2016 AIP Publishing LLC.

Molecular dynamics and simulations study on the vibrational and electronic solvatochromism of benzophenone

Solvent plays a key role in diverse physico-chemical and biological processes. Therefore, understanding solute-solvent interactions at the molecular level of detail is of utmost importance. A comprehensive solvatochromic analysis of benzophenone (Bzp) was carried out in various solvents using Raman and electronic spectroscopy, in conjunction with Density Functional Theory (DFT) calculations of supramolecular solute-solvent clusters generated using classical Molecular Dynamics Simulations (c-MDSs). The >C=O stretching frequency undergoes a bathochromic shift with solvent polarity. Interestingly, in protic solvents this peak appears as a doublet: c-MDS and ad hoc explicit solvent ab initio calculations suggest that the lower and higher frequency peaks are associated with the hydrogen bonded and dangling carbonyl group of Bzp, respectively. Additionally, the dangling carbonyl in methanol (MeOH) solvent is 4 cm(-1) blue-shifted relative to acetonitrile solvent, despite their similar dipolarity/polarizability. This suggests that the cybotactic region of the dangling carbonyl group in MeOH is very different from its bulk solvent structure. Therefore, we propose that this blue-shift of the dangling carbonyl originates in the hydrophobic solvation shell around it resulting from extended hydrogen bonding network of the protic solvents. Furthermore, the 1(1)n pi* (band I) and 1(1)pi pi* (band II) electronic transitions show a hypsochromic and bathochromic shift, respectively. In particular, these shifts in protic solvents are due to differences in their excited state-hydrogen bonding mechanisms. Additionally, a linear relationship is obtained for band I and the >C=O stretching frequency (cm(-1)), which suggests that the different excitation wavelengths in band I correspond to different solvation states. Therefore, we hypothesize that the variation in excitation wavelengths in band I could arise from different solvation states leading to varying solvation dynamics. This will have implications for ultrafast processes associated with electron-transfer, charge transfer, and also the photophysical aspects of excited states. (C) 2016 AIP Publishing LLC.

Primary and rechargeable zinc-air batteries using ceramic and highly stable TiCN as an oxygen reduction reaction electrocatalyst

Primary and secondary zinc-air batteries based on ceramic, stable, one dimensional titanium carbonitride (TiCN) nanostructures are reported. The optimized titanium carbonitride composition by density functional theory reveals their good activity towards the oxygen reduction reaction (ORR). Electrochemical measurements show their superior performance for the ORR in alkaline media coupled with favourable kinetics. The nanostructured TiCN lends itself amenable to be used as an air cathode material in primary and rechargeable zinc-air batteries. The battery performance and cyclability are found to be good. Further, we have demonstrated a gel-based electrolyte for rechargeable zinc-air batteries based on a TiCN cathode under ambient, atmospheric conditions without any oxygen supply from a cylinder. The present cell can work at current densities of 10-20 mA cm(2) (app. 10 000 mA g(-1) of TiCN) for several hours (63 h in the case of 10 mA cm(-2)) with a charge retention of 98%. The low cost, noble metal-free, mechanically stable and corrosion resistant TiCN is a very good alternative to Pt for metal-air battery chemistry.

Na2M2(SO4)(3) (M = Fe, Mn, Co and Ni): towards high-voltage sodium battery applications

Sodium-ion-based batteries have evolved as excellent alternatives to their lithium-ion-based counterparts due to the abundance, uniform geographical distribution and low price of Na resources. In the pursuit of sodium chemistry, recently the alluaudite framework Na2M2(SO4)(3) has been unveiled as a high-voltage sodium insertion system. In this context, the framework of density functional theory has been applied to systematically investigate the crystal structure evolution, density of states and charge transfer with sodium ions insertion, and the corresponding average redox potential, for Na2M2(SO4)(3) (M = Fe, Mn, Co and Ni). It is shown that full removal of sodium atoms from the Fe-based device is not a favorable process due to the 8% volume shrinkage. The imaginary frequencies obtained in the phonon dispersion also reflect this instability and the possible phase transition. This high volume change has not been observed in the cases of the Co- and Ni-based compounds. This is because the redox reaction assumes a different mechanism for each of the compounds investigated. For the polyanion with Fe, the removal of sodium ions induces a charge reorganization at the Fe centers. For the Mn case, the redox process induces a charge reorganization of the Mn centers with a small participation of the oxygen atoms. The Co and Ni compounds present a distinct trend with the redox reaction occurring with a strong participation of the oxygen sublattice, resulting in a very small volume change upon desodiation. Moreover, the average deintercalation potential for each of the compounds has been computed. The implications of our findings have been discussed both from the scientific perspective and in terms of technological aspects.

Triplet excited electronic state switching induced by hydrogen bonding: A transient absorption spectroscopy and time-dependent DFT study

The solvent plays a decisive role in the photochemistry and photophysics of aromatic ketones. Xanthone (XT) is one such aromatic ketone and its triplet-triplet (T-T) absorption spectra show intriguing solvatochromic behavior. Also, the reactivity of XT towards H-atom abstraction shows an unprecedented decrease in protic solvents relative to aprotic solvents. Therefore, a comprehensive solvatochromic analysis of the triplet-triplet absorption spectra of XT was carried out in conjunction with time dependent density functional theory using the ad hoc explicit solvent model approach. A detailed solvatochromic analysis of the T-T absorption bands of XT suggests that the hydrogen bonding interactions are different in the corresponding triplet excited states. Furthermore, the contributions of non-specific and hydrogen bonding interactions towards differential solvation of the triplet states in protic solvents were found to be of equal magnitude. The frontier molecular orbital and electron density difference analysis of the T-1 and T-2 states of XT indicates that the charge redistribution in these states leads to intermolecular hydrogen bond strengthening and weakening, respectively, relative to the S-0 state. This is further supported by the vertical excitation energy calculations of the XT-methanol supra-molecular complex. The intermolecular hydrogen bonding potential energy curves obtained for this complex in the S-0, T-1, and T-2 states support the model. In summary, we propose that the different hydrogen bonding mechanisms exhibited by the two lowest triplet excited states of XT result in a decreasing role of the n pi* triplet state, and are thus responsible for its reduced reactivity towards H-atom abstraction in protic solvents. (C) 2016 AIP Publishing LLC.

Simultaneous tunability of the electronic and phononic gaps in SnS2 under normal compressive strain

Controlled variation of the electronic properties of. two-dimensional (2D) materials by applying strain has emerged as a promising way to design materials for customized applications. Using density functional theory (DFT) calculations, we show that while the electronic structure and indirect band gap of SnS2 do not change significantly with the number of layers, they can be reversibly tuned by applying biaxial tensile (BT), biaxial compressive (BC), and normal compressive (NC) strains. Mono to multilayered SnS2 exhibit a reversible semiconductor to metal (S-M) transition with applied strain. For bilayer (2L) SnS2, the S-Mtransition occurs at the strain values of 17%,-26%, and -24% under BT, BC, and NC strains, respectively. Due to weaker interlayer coupling, the critical strain value required to achieve the S-Mtransition in SnS2 under NC strain is much higher than for MoS2. From a stability viewpoint, SnS2 becomes unstable at very low strain values on applying BC (-6.5%) and BT strains (4.9%), while it is stable even up to the transition point (-24%) in the case of NC strain. In addition to the reversible tuning of the electronic properties of SnS2, we also show tunability in the phononic band gap of SnS2, which increases with applied NC strain. This gap increases three times faster than for MoS2. This simultaneous tunability of SnS2 at the electronic and phononic levels with strain, makes it a potential candidate in field effect transistors (FETs) and sensors as well as frequency filter applications.

Visualisation of charge-transfer excitations in donor-acceptor molecules using the particle-hole map: a case study

Charge-transfer (CT) excitations are essential for photovoltaic phenomena in organic solar cells. Owing to the complexity of molecular geometries and orbital coupling, a detailed analysis and spatial visualisation of CT processes can be challenging. In this paper, a new detail-oriented visualisation scheme, the particle-hole map (PHM), is applied and explained for the purpose of spatial analysis of excitations in organic molecules. The PHM can be obtained from the output of a time-dependent density-functional theory calculation with negligible additional computational cost, and provides a useful physical picture for understanding the origins and destinations of electrons and holes during an excitation process. As an example, we consider intramolecular CT excitations in Diketopyrrolopyrrole-based molecules, and relate our findings to experimental results.

Theoretical insights on the electro-thermal transport properties of monolayer MoS2 with line defects

Two dimensional (2D) materials demonstrate several novel electrical, mechanical, and thermal properties which are quite distinctive to those of their bulk form. Among many others, one important potential application of the 2D material is its use in the field of energy harvesting. Owing to that, here we present a detailed study on electrical as well as thermal transport of monolayer MoS2, in quasi ballistic regime. Besides the perfect monolayer in its pristine form, we also consider various line defects which have been experimentally observed in mechanically exfoliated MoS2 samples. For calculating various parameters related to the electrical transmission, we employ the non-equilibrium Green's function-density functional theory combination. However, to obtain the phonon transmission, we take help of the parametrized Stillinger-Weber potential which can accurately delineate the inter-atomic interactions for the monolayer MoS2. Due to the presence of line defects, we observed significant reductions in both the charge carrier and the phonon transmissions through a monolayer MoS2 flake. Moreover, we also report a comparative analysis showing the temperature dependency of the thermoelectric figure of merit values, as obtained for the perfect as well as the other defective 2D samples. (C) 2016 AIP Publishing LLC.

Elasticity, stability, and ideal strength of beta-SiC in plane-wave-based ab initio calculations

On the basis of the pseudopotential plane-wave method and the local-density-functional theory, this paper studies energetics, stress-strain relation, stability, and ideal strength of beta-SiC under various loading modes, where uniform uniaxial extension and tension and biaxial proportional extension are considered along directions [001] and [111]. The lattice constant, elastic constants, and moduli of equilibrium state are calculated and the results agree well with the experimental data. As the four SI-C bonds along directions [111], [(1) over bar 11], [11(1) over bar] and [111] are not the same under the loading along [111], internal relaxation and the corresponding internal displacements must be considered. We find that, at the beginning of loading, the effect of internal displacement through the shuffle and glide plane diminishes the difference among the four Si-C bonds lengths, but will increase the difference at the subsequent loading, which will result in a crack nucleated on the {111} shuffle plane and a subsequently cleavage fracture. Thus the corresponding theoretical strength is 50.8 GPa, which agrees well with the recent experiment value, 53.4 GPa. However, with the loading along [001], internal relaxation is not important for tetragonal symmetry. Elastic constants during the uniaxial tension along [001] are calculated. Based on the stability analysis with stiffness coefficients, we find that the spinodal and Born instabilities are triggered almost at the same strain, which agrees with the previous molecular-dynamics simulation. During biaxial proportional extension, stress and strength vary proportionally with the biaxial loading ratio at the same longitudinal strain.

Development of a classical force field for the oxidized Si surface: application to hydrophilic wafer bonding.

We have developed a classical two- and three-body interaction potential to simulate the hydroxylated, natively oxidized Si surface in contact with water solutions, based on the combination and extension of the Stillinger-Weber potential and of a potential originally developed to simulate SiO(2) polymorphs. The potential parameters are chosen to reproduce the structure, charge distribution, tensile surface stress, and interactions with single water molecules of a natively oxidized Si surface model previously obtained by means of accurate density functional theory simulations. We have applied the potential to the case of hydrophilic silicon wafer bonding at room temperature, revealing maximum room temperature work of adhesion values for natively oxidized and amorphous silica surfaces of 97 and 90 mJm(2), respectively, at a water adsorption coverage of approximately 1 ML. The difference arises from the stronger interaction of the natively oxidized surface with liquid water, resulting in a higher heat of immersion (203 vs 166 mJm(2)), and may be explained in terms of the more pronounced water structuring close to the surface in alternating layers of larger and smaller densities with respect to the liquid bulk. The computed force-displacement bonding curves may be a useful input for cohesive zone models where both the topographic details of the surfaces and the dependence of the attractive force on the initial surface separation and wetting can be taken into account.

Ab initio investigation of the elasticity and stability of aluminium

On the basis of the pseudopotential plane-wave (PP-PW) method in combination with the local density functional theory (LDFT), complete stress-strain curves for the uniaxial loading and uniaxial deformation along the [001] and [111] directions, and the biaxial proportional extension along [010] and [001] for aluminium are obtained. During the uniaxial loading, certain general behaviours of the energy versus the stretch and the load versus the stretch are confirmed; in each case, there exist three special unstressed structures: f.c.c., b.c.c., and f.c.t. for [001]; f.c.c., s.c., and b.c.c. for [111]. Using stability criteria, we find that all of these states are unstable, and always occur together with shear instability, except the natural f.c.c. structure. A Pain transformation from the stable f.c.c. structure to the stable b.c.c. configuration cannot be obtained by uniaxial compression along any equivalent [001] and [111] direction. The tensile strengths are similar for the two directions. For the higher energy barrier of the [111] direction, the compressive strength is greater than that for the [001] direction. With increase in the ratio of the biaxial proportional extension, the stress and tensile strength increase; however, the critical strain does not change significantly. Our results add to the existing ab initio database for use in fitting and testing interatomic potentials.

Graphite intercalation compounds under pressure

Motivated by recent experimental work, we use first-principles density functional theory methods to conduct an extensive search for low enthalpy structures of C$_6$Ca under pressure. As well as a range of buckled structures, which are energetically competitive over an intermediate range of pressures, we show that the high pressure system ($\gtrsim 18$ GPa) is unstable towards the formation of a novel class of layered structures, with the most stable compound involving carbon sheets containing five- and eight-membered rings. As well as discussing the energetics of the different classes of low enthalpy structures, we comment on the electronic structure of the high pressure compound and its implications for superconductivity.

Growth Mechanism And Joint Structure Of Zno Tetrapods

Regular zinc oxide (ZnO) tetrapods with a flat plane have been obtained on Si(1 0 0) substrate via the chemical vapour deposition approach. The x-ray diffraction result suggests that these tetrapods are all single crystals with a wurtzite structure that grow along the (0 0 0 1) direction and corresponding electron backscatter diffraction analysis reveals the crystal orientation of growth and exposed surface. Furthermore, we find some ZnO tetrapods with some legs off and the angles between every two legs are measured with the aid of scanning electron microscopy and image analysis, which benefit to reveal the structure of ZnO tetrapods joint. The structure model and growth mechanism of ZnO tetrapods are proposed. Besides, the stable model of the interface was obtained through the density-functional theory calculation and the energy needed to break the twin plane junction was calculated as 5.651 J m(-2).

Effect Of Load Triaxiality On Polymorphic Transitions In Zinc Oxide

Molecular dynamics (MD) simulations and first-principles calculations are carried out to analyze the stability of both newly discovered and previously known phases of ZnO under loading of various triaxialities. The analysis focuses on a graphite-like phase (FIX) and a body-centered-tetragonal phase (BCT-4) that were observed recently in [0 1 (1) over bar 0]- and [0 0 0 1]-oriented nanowires respectively under uniaxial tensile loading as well as the natural state of wurtzite (WZ) and the rocksalt (RS) phase which exists under hydrostatic pressure loading. Equilibrium critical stresses for the transformations are obtained. The WZ -> HX transformation is found to be energetically favorable above a critical tensile stress of 10 GPa in [0 1 (1) over tilde 0] nanowires. The BCT-4 phase can be stabilized at tensile stresses above 7 GPa in [0 0 0 1] nanowires. The RS phase is stable at hydrostatic pressures above 8.2 GPa. The identification and characterization of these phase transformations reveal a more extensive polymorphism of ZnO than previously known. A crystalline structure-load triaxiality map is developed to summarize the new understanding. (c) 2007 Elsevier Ltd. All rights reserved.

Unconventional spin texture of a topologically nontrivial semimetal Sb(110)

9 p.