Quadratic Nonlinearity of One- and Two-Electron Oxidized Metalloporphyrins and Their Switching in Solution

We report the quadratic nonlinearity of one- and two-electron oxidation products of the first series of transition metal complexes of meso-tetraphenylporphyrin (TPP). Among many MTPP complexes, only CuTPP and ZnTPP show reversible oxidation/reduction cycles as seen from cyclic voltammetry experiments. While centrosymmetric neutral metalloporphyrins have zero first hyperpolarizability, β, as expected, the cation radicals and dications of CuTPP and ZnTPP have very high β values. The one- and two-electron oxidation of the MTPPs leads to symmetry-breaking of the metal−porphyrin core, resulting in a large β value that is perhaps aided in part by contributions from the two-photon resonance enhancement. The calculated static first hyperpolarizabilities, β0, which are evaluated in the framework of density functional theory by a coupled perturbed Hartree−Fock method, support the experimental trend. The switching of optical nonlinearity has been achieved between the neutral and the one-electron oxidation products but not between the one- and the two-electron oxidation products since dications that are electrochemically reversible are unstable due to the formation of stable isoporphyrins in the presence of nucleophiles such as halides.

Molecular theory of solvation and solvation dynamics of a classical ion in a dipolar liquid

A microscopic theory of equilibrium solvation and solvation dynamics of a classical, polar, solute molecule in dipolar solvent is presented. Density functional theory is used to explicitly calculate the polarization structure around a solvated ion. The calculated solvent polarization structure is different from the continuum model prediction in several respects. The value of the polarization at the surface of the ion is less than the continuum value. The solvent polarization also exhibits small oscillations in space near the ion. We show that, under certain approximations, our linear equilibrium theory reduces to the nonlocal electrostatic theory, with the dielectric function (c(k)) of the liquid now wave vector (k) dependent. It is further shown that the nonlocal electrostatic estimate of solvation energy, with a microscopic c(k), is close to the estimate of linearized equilibrium theories of polar liquids. The study of solvation dynamics is based on a generalized Smoluchowski equation with a mean-field force term to take into account the effects of intermolecular interactions. This study incorporates the local distortion of the solvent structure near the ion and also the effects of the translational modes of the solvent molecules.The latter contribution, if significant, can considerably accelerate the relaxation of solvent polarization and can even give rise to a long time decay that agrees with the continuum model prediction. The significance of these results is discussed.

Density-wave theory of dislocations in crystals

The density-wave theory of Ramakrishnan and Yussouff is extended to provide a scheme for describing dislocations and other topological defects in crystals. Quantitative calculations are presented for the order-parameter profiles, the atomic configuration, and the free energy of a screw dislocation with Burgers vector b=(a/2, a/2, a/2) in a bcc solid. These calculations are done using a simple parametrization of the direct correlation function and a gradient expansion. It is conventional to express the free energy of the dislocation in a crystal of size R as (λb2/4π)ln(αR/‖b‖), where λ is the shear elastic constant, and α is a measure of the core energy. Our results yield for Na the value α≃1.94a/(‖c1’’‖)1/2 (≃1.85) at the freezing temperature (371 K) and α≃2.48a/(‖c1’’‖)1/2 at 271 K, where c1’’ is the curvature of the first peak of the direct correlation function c(q). Detailed results for the density distribution in the dislocation, particularly the core region, are also presented. These show that the dislocation core has a columnar character. To our knowledge, this study represents the first calculation of dislocation structure, including the core, within the framework of an order-parameter theory and incorporating thermal effects.

Topological defects in crystals: A density-wave theory

A new approach for describing dislocations and other topological defects in crystals, based on the density wave theory of Ramakrishnan and Yussouff is presented. Quantitative calculations are discussed in brief for the order parameter profiles, the atomic configuration and the free energy of a screw dislocation with Burgers vector b = (a/2, a/2,a/2 ) in a bcc solid. Our results for the free energy of the dislocation in a crystal of sizeR, when expressed as (λb 2/4π) ln (αR/|b|) whereλ is the shear elastic constant, yield, for example, the valueα ⋍ 1·85 for sodium at its freezing temperature (371°K). The density distribution in the presence of the dislocation shows that the dislocation core has a columnar character. To our knowledge, this study represents the first calculation of dislocation structure, including the core, within the framework of an order parameter theory incorporating thermal effects.

Armchair or Zigzag? A tool for characterizing graphene edge

Electronic, magnetic, and structural properties of graphene flakes depend sensitively upon the type of edge atoms. We present a simple software tool for determining the type of edge atoms in a honeycomb lattice. The algorithm is based on nearest neighbor counting. Whether an edge atom is of armchair or zigzag type is decided by the unique pattern of its nearest neighbors. Particular attention is paid to the practical aspects of using the tool, as additional features such as extracting out the edges from the lattice could help in analyzing images from transmission microscopy or other experimental probes. Ultimately, the tool in combination with density-functional theory or tight-binding method can also be helpful in correlating the properties of graphene flakes with the different armchair-to-zigzag ratios. Program summary Program title: edgecount Catalogue identifier: AEIA_v1_0 Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AEIA_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data, etc.: 66685 No. of bytes in distributed program, including test data, etc.: 485 381 Distribution format: tar.gz Programming language: FORTRAN 90/95 Computer: Most UNIX-based platforms Operating system: Linux, Mac OS Classification: 16.1, 7.8 Nature of problem: Detection and classification of edge atoms in a finite patch of honeycomb lattice. Solution method: Build nearest neighbor (NN) list; assign types to edge atoms on the basis of their NN pattern. Running time: Typically similar to second(s) for all examples. (C) 2010 Elsevier B.V. All rights reserved.

Activation of Lattice Oxygen of TiO2 by Pd2+ Ion: Correlation of Low-Temperature CO and Hydrocarbon Oxidation with Structure of Ti1-xPdxO2-x (x=0.01-0.03)

Lattice oxygen of TiO2 is activated by the substitution of Pd ion in its lattice. Ti1-xPdxO2-x (x = 0.01-0.03) have been synthesized by solution combustion method crystallizing in anatase TiO2 structure. Pd is in +2 oxidation state and Ti is in +4 oxidation state in the catalyst. Pd is more ionic in TiO2 lattice compared to Pd in PdO. Oxygen storage capacity defined by ``amount of oxygen that is used reversibly to oxidize CO'' is as high as 5100 mu mol/g of Ti0.97Pd0.03O1.97. Oxygen is extracted by CO to CO2 in absence of feed oxygen even at room temperature which is more than 20 times compared to pure TiO2. Rate of CO oxidation is 2.75 mu mol g(-1) s(-1) at 60 degrees C over Ti0.97Pd0.03O1.97 and C2H2 gets oxidized to CO2 and H2O at room temperature. Catalyst is not poisoned on long time operation of the reactor. Such high catalytic activity is due to activated lattice oxygen created by the substitution of Pd ion as seen from first-principles density functional theory (DFT) calculations with 96 atom supercells of Ti32O64, Ti31Pd1O63, Ti30Pd2O62, and Ti29Pd3O61. The compounds crystallize in anatase TiO2 structure with Pd2+ ion in nearly square planar geometry and TiO6 octahedra are distorted by the creation of weakly bound oxygens. Structural analysis of Ti31Pd1O63 which is close to 3% Pd ion substituted TiO2 shows that oxygens associated with both Ti and Pd ions in the lattice show bond valence sum of 1.87, a low value characteristic of weak oxygen in the lattice compared to oxygens with valence 2 and above in the same lattice. Exact positions of activated oxygens have been identified in the lattice from DFT calculations.

Exploring the color of transition metal ions in irregular coordination geometries: new colored inorganic oxides based on the spiroffite structure, Zn2-xMxTe3O8 (M = Co, Ni, Cu)

We describe the synthesis, crystal structures, and optical absorption spectra of transition metal substituted spiroffite derivatives, Zn2-xMxTe3O8 (M-II = Co, Ni, Cu; 0 < x <= 1.0). The oxides are readily synthesized by solid state reaction of stoichiometric mixtures of the constituent binaries at 620 degrees C. Reitveld refinement of the crystal structures from powder X-ray diffraction (XRD) data shows that the Zn/MO6 octahedra are strongly distorted, as in the parent Zn2Te3O8 structure, consisting of five relatively short Zn/M-II-O bonds (1.898-2.236 angstrom) and one longer Zn/M-II-O bond (2.356-2.519 angstrom). We have interpreted the unique colors and the optical absorption/diffuse reflectance spectra of Zn2-xMxTe3O8 in the visible, in terms of the observed/irregular coordination geometry of the Zn/M-II-O chromophores. We could not however prepare the fully substituted M2Te3O8 (M-II = Co, Ni, Cu) by the direct solid state reaction method. Density Functional Theory (DFT) modeling of the electronic structure of both the parent and the transition metal substituted derivatives provides new insights into the bonding and the role of transition metals toward the origin of color in these materials. We believe that transition metal substituted spiroffites Zn2-xMxTe3O8 reported here suggest new directions for the development of colored inorganic materials/pigments featuring irregular/distorted oxygen coordination polyhedra around transition metal ions.

Performance Analysis of Strained Monolayer MoS2 MOSFET

We present a computational study on the impact of tensile/compressive uniaxial (epsilon(xx)) and biaxial (epsilon(xx) = epsilon(yy)) strain on monolayer MoS2, n-, and p-MOSFETs. The material properties like band structure, carrier effective mass, and the multiband Hamiltonian of the channel are evaluated using the density functional theory. Using these parameters, self-consistent Poisson-Schrodinger solution under the nonequilibrium Green's function formalism is carried out to simulate the MOS device characteristics. 1.75% uniaxial tensile strain is found to provide a minor (6%) ON current improvement for the n-MOSFET, whereas same amount of biaxial tensile strain is found to considerably improve the p-MOSFET ON currents by 2-3 times. Compressive strain, however, degrades both n-MOS and p-MOS devices performance. It is also observed that the improvement in p-MOSFET can be attained only when the channel material becomes indirect gap in nature. We further study the performance degradation in the quasi-ballistic long-channel regime using a projected current method.

Platinum Ion-Doped TiO2: High Catalytic Activity of Pt2+ with Oxide Ion Vacancy in Ti1-x4+Ptx2+O2-x Compared to Pt4+ without Oxide Ion Vacancy in Ti1-x4+Ptx4+O2

Ti0.97Pt0.032+O1.97 and Ti0.97Pt0.034+O2 have been synthesized by a solution combustion method using alanine and glycine as the fuels, respectively. Both crystallize in anatase TiO2 structure with 15 nm average crystallite size. X-ray photoelectron spectroscopy (XPS) confirmed Pt ions are in the 2+ state in Ti0.97Pt0.03O1.97 (alanine) and 4+ state in Ti0.97Pt0.03O2 (glycine). The rate of CO oxidation occurring over Ti0.97Pt0.032+O1.97 (0.76 mu mol.g(-1).s(-1)) is similar to 10, times more than that over Ti0.97Pt0.034+O2 at 60 degrees C (0.08 mu mol.g(-1).s(-1)). A large shift in 100% hydrocarbons conversion to lower temperature was observed for Pt2+ ion-substituted TiO2 relative 10 that for Pt4+ ion-substituted TiO2. After reoxidation of the reduced compound by H-2 as well as CO, Pt ions are stabilized in mixed valences, 2+ and 4+ states. The role of oxide ion vacancy has been demonstrated by CO oxidation and H-2 + O-2 recombination reactions in the presence and absence of O-2. We analyze the activated lattice oxygens upon substitution of Pt2+ and Pt4+ ions in TiO2, using first-principles density functional theory (DFT) calculations with supercells of Ti31Pt1O63, Ti30Pt2O62, and Ti29Pt3O61 for Pt2+ ion substitution and Ti31Pt1O64, Ti30Pt2O62, and Ti29Pt3O61 for Pt4+ ion substitution in TiO2. We find that the local structure of Pt2+ ion has a distorted square planar geometry and that of Pt4+ ion has an octahedral geometry similar to that of Ti4+ ion in pure TiO2. The change in coordination of Pt2+ ion gives rise to weakly bonded oxygens, and these oxygens are involved in high rates of catalytic reaction. Thus, the high catalytic activity results from synergistic roles of Pt2+ ion and oxide ion vacancy and weakly bonded lattice oxygen.

Quantum chemical study on influence of intermolecular hydrogen bonding on the geometry, the atomic charges and the vibrational dynamics of 2,6-dichlorobenzonitrile

FT-IR (4000-400 cm(-1)) and FT-Raman (4000-200 cm(-1)) spectral measurements on solid 2,6-dichlorobenzonitrile (2,6-DCBN) have been done. The molecular geometry, harmonic vibrational frequencies and bonding features in the ground state have been calculated by density functional theory at the B3LYP/6-311++G (d,p) level. A comparison between the calculated and the experimental results covering the molecular structure has been made. The assignments of the fundamental vibrational modes have been done on the basis of the potential energy distribution (PED). To investigate the influence of intermolecular hydrogen bonding on the geometry, the charge distribution and the vibrational spectrum of 2,6-DCBN; calculations have been done for the monomer as well as the tetramer. The intermolecular interaction energies corrected for basis set superposition error (BSSE) have been calculated using counterpoise method. Based on these results, the correlations between the vibrational modes and the structure of the tetramer have been discussed. Molecular electrostatic potential (MEP) contour map has been plotted in order to predict how different geometries could interact. The Natural Bond Orbital (NBO) analysis has been done for the chemical interpretation of hyperconjugative interactions and electron density transfer between occupied (bonding or lone pair) orbitals to unoccupied (antibonding or Rydberg) orbitals. UV spectrum was measured in methanol solution. The energies and oscillator strengths were calculated by Time Dependent Density Functional Theory (TD-DFT) and matched to the experimental findings. TD-DFT method has also been used for theoretically studying the hydrogen bonding dynamics by monitoring the spectral shifts of some characteristic vibrational modes involved in the formation of hydrogen bonds in the ground and the first excited state. The C-13 nuclear magnetic resonance (NMR) chemical shifts of the molecule were calculated by the Gauge independent atomic orbital (GIAO) method and compared with experimental results. Standard thermodynamic functions have been obtained and changes in thermodynamic properties on going from monomer to tetramer have been presented. (C) 2013 Elsevier B.V. All rights reserved.

Density Functional Theoretical Modeling, Electrostatic Surface Potential and Surface Enhanced Raman Spectroscopic Studies on Biosynthesized Silver Nanoparticles: Observation of 400 pM Sensitivity to Explosives

Interaction of adsorbate on charged surfaces, orientation of the analyte on the surface, and surface enhancement aspects have been studied. These aspects have been explored in details to explain the surface-enhanced Raman spectroscopic (SERS) spectra of 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane (HNIW or CL-20), a well-known explosive, and 2,4,6-trinitrotoluene (TNT) using one-pot synthesis of silver nanoparticles via biosynthetic route using natural precursor extracts of clove and pepper. The biosynthesized silver nanoparticles (bio Ag Nps) have been characterized using UV-vis spectroscopy, scanning electron microscopy and atomic force microscopy. SERS studies conducted using bio Ag Nps on different water insoluble analytes, such as CL-20 and TNT, lead to SERS signals at concentration levels of 400 pM. The experimental findings have been corroborated with density functional computational results, electrostatic surface potential calculations, Fukui functions and potential measurements.

Intermolecular interactions, charge-density distribution and the electrostatic properties of pyrazinamide anti-TB drug molecule: an experimental and theoretical charge-density study

An experimental charge-density analysis of pyrazinamide (a first line antitubercular drug) was performed using high-resolution X-ray diffraction data (sin theta/lambda)(max) = 1.1 angstrom(-1)] measured at 100 (2) K. The structure was solved by direct methods using SHELXS97 and refined by SHELXL97. The total electron density of the pyrazinamide molecule was modeled using the Hansen-Coppens multipole formalism implemented in the XD software. The topological properties of electron density determined from the experiment were compared with the theoretical results obtained from CRYSTAL09 at the B3LYP/6-31G** level of theory. The crystal structure was stabilized by N-H center dot center dot center dot N and N-H center dot center dot center dot O hydrogen bonds, in which the N3-H3B center dot center dot center dot N1 and N3-H3A center dot center dot center dot O1 interactions form two types of dimers in the crystal. Hirshfeld surface analysis was carried out to analyze the intermolecular interactions. The fingerprint plot reveals that the N center dot center dot center dot H and O center dot center dot center dot H hydrogen-bonding interactions contribute 26.1 and 18.4%, respectively, of the total Hirshfeld surface. The lattice energy of the molecule was calculated using density functional theory (B3LYP) methods with the 6-31G** basis set. The molecular electrostatic potential of the pyrazinamide molecule exhibits extended electronegative regions around O1, N1 and N2. The existence of a negative electrostatic potential (ESP) region just above the upper and lower surfaces of the pyrazine ring confirm the pi-electron cloud.

Improved quasiparticle wave functions and mean field for G(0)W(0) calculations: Initialization with the COHSEX operator

The GW approximation to the electron self-energy has become a standard method for ab initio calculation of excited-state properties of condensed-matter systems. In many calculations, the G W self-energy operator, E, is taken to be diagonal in the density functional theory (DFT) Kohn-Sham basis within the G0 W0 scheme. However, there are known situations in which this diagonal Go Wo approximation starting from DFT is inadequate. We present two schemes to resolve such problems. The first, which we called sc-COHSEX-PG W, involves construction of an improved mean field using the static limit of GW, known as COHSEX (Coulomb hole and screened exchange), which is significantly simpler to treat than GW W. In this scheme, frequency-dependent self energy E(N), is constructed and taken to be diagonal in the COHSEX orbitals after the system is solved self-consistently within this formalism. The second method is called off diagonal-COHSEX G W (od-COHSEX-PG W). In this method, one does not self-consistently change the mean-field starting point but diagonalizes the COHSEX Hamiltonian within the Kohn-Sham basis to obtain quasiparticle wave functions and uses the resulting orbitals to construct the G W E in the diagonal form. We apply both methods to a molecular system, silane, and to two bulk systems, Si and Ge under pressure. For silane, both methods give good quasiparticle wave functions and energies. Both methods give good band gaps for bulk silicon and maintain good agreement with experiment. Further, the sc-COHSEX-PGW method solves the qualitatively incorrect DFT mean-field starting point (having a band overlap) in bulk Ge under pressure.

Stacking Interactions in RNA and DNA: Roll-Slide Energy Hyperspace for Ten Unique Dinucleotide Steps

Understanding dinucleotide sequence directed structures of nuleic acids and their variability from experimental observation remained ineffective due to unavailability of statistically meaningful data. We have attempted to understand this from energy scan along twist, roll, and slide degrees of freedom which are mostly dependent on dinucleotide sequence using ab initio density functional theory. We have carried out stacking energy analysis in these dinucleotide parameter phase space for all ten unique dinucleotide steps in DNA and RNA using DFT-D by B97X-D/6-31G(2d,2p), which appears to satisfactorily explain conformational preferences for AU/AU step in our recent study. We show that values of roll, slide, and twist of most of the dinucleotide sequences in crystal structures fall in the low energy region. The minimum energy regions with large twist values are associated with the roll and slide values of B-DNA, whereas, smaller twist values correspond to higher stability to RNA and A-DNA like conformations. Incorporation of solvent effect by CPCM method could explain the preference shown by some sequences to occur in B-DNA or A-DNA conformations. Conformational preference of BII sub-state in B-DNA is preferentially displayed mainly by pyrimidine-purine steps and partly by purine-purine steps. The purine-pyrimidine steps show largest effect of 5-methyl group of thymine in stacking energy and the introduction of solvent reduces this effect significantly. These predicted structures and variabilities can explain the effect of sequence on DNA and RNA functionality. (c) 2014 Wiley Periodicals, Inc. Biopolymers 103: 134-147, 2015.

The Self-Association of Graphane Is Driven by London Dispersion and Enhanced Orbital Interactions

We investigated the nature of the cohesive energy between graphane sheets via multiple CH center dot center dot center dot HC interactions, using density functional theory (DFT) including dispersion correction (Grimmes D3 approach) computations of n]graphane sigma dimers (n = 6-73). For comparison, we also evaluated the binding between graphene sheets that display prototypical pi/pi interactions. The results were analyzed using the block-localized wave function (BLW) method, which is a variant of ab initio valence bond (VB) theory. BLW interprets the intermolecular interactions in terms of frozen interaction energy (Delta E-F) composed of electrostatic and Pauli repulsion interactions, polarization (Delta E-pol), charge-transfer interaction (Delta E-CT), and dispersion effects (Delta E-disp). The BLW analysis reveals that the cohesive energy between graphane sheets is dominated by two stabilizing effects, namely intermolecular London dispersion and two-way charge transfer energy due to the sigma CH -> sigma*(HC) interactions. The shift of the electron density around the nonpolar covalent C-H bonds involved in the intermolecular interaction decreases the C-H bond lengths uniformly by 0.001 angstrom. The Delta E-CT term, which accounts for similar to 15% of the total binding energy, results in the accumulation of electron density in the interface area between two layers. This accumulated electron density thus acts as an electronic glue for the graphane layers and constitutes an important driving force in the self-association and stability of graphane under ambient conditions. Similarly, the double faced adhesive tape style of charge transfer interactions was also observed among graphene sheets in which it accounts for similar to 18% of the total binding energy. The binding energy between graphane sheets is additive and can be expressed as a sum of CH center dot center dot center dot HC interactions, or as a function of the number of C-H bonds.