Ab initio quantum chemical studies on neutral-radical reactions of ethynyl (C2H) and cyano (CN) with unsaturated hydrocarbons

An Ab Initio/RRKM study of the reaction mechanism and product branching ratios of neutral-radical ethynyl (C2H) and cyano (CN) radical species with unsaturated hydrocarbons is performed. The reactions studied apply to cold conditions such as planetary atmospheres including Titan, the Interstellar Medium (ISM), icy bodies and molecular clouds. The reactions of C2H and CN additions to gaseous unsaturated hydrocarbons are an active area of study. NASA's Cassini/Huygens mission found a high concentration of C2H and CN from photolysis of ethyne (C2H2) and hydrogen cyanide (HCN), respectively, in the organic haze layers of the atmosphere of Titan. The reactions involved in the atmospheric chemistry of Titan lead to a vast array of larger, more complex intermediates and products and may also serve as a chemical model of Earth's primordial atmospheric conditions. The C2H and CN additions are rapid and exothermic, and often occur barrierlessly to various carbon sites of unsaturated hydrocarbons. The reaction mechanism is proposed on the basis of the resulting potential energy surface (PES) that includes all the possible intermediates and transition states that can occur, and all the products that lie on the surface. The B3LYP/6-311g(d,p) level of theory is employed to determine optimized electronic structures, moments of inertia, vibrational frequencies, and zero-point energy. They are followed by single point higher-level CCSD(T)/cc-vtz calculations, including extrapolations to complete basis sets (CBS) of the reactants and products. A microcanonical RRKM study predicts single-collision (zero-pressure limit) rate constants of all reaction paths on the potential energy surface, which is then used to compute the branching ratios of the products that result. These theoretical calculations are conducted either jointly or in parallel to experimental work to elucidate the chemical composition of Titan's atmosphere, the ISM, and cold celestial bodies.<.

Ab Initio Quantum Chemical Studies on Neutral-Radical Reactions of Ethynyl (C2H) and Cyano (CN) with Unsaturated Hydrocarbons

An Ab Initio/RRKM study of the reaction mechanism and product branching ratios of neutral-radical ethynyl (C2H) and cyano (CN) radical species with unsaturated hydrocarbons is performed. The reactions studied apply to cold conditions such as planetary atmospheres including Titan, the Interstellar Medium (ISM), icy bodies and molecular clouds. The reactions of C2H and CN additions to gaseous unsaturated hydrocarbons are an active area of study. NASA’s Cassini/Huygens mission found a high concentration of C2H and CN from photolysis of ethyne (C2H2) and hydrogen cyanide (HCN), respectively, in the organic haze layers of the atmosphere of Titan. The reactions involved in the atmospheric chemistry of Titan lead to a vast array of larger, more complex intermediates and products and may also serve as a chemical model of Earth’s primordial atmospheric conditions. The C2H and CN additions are rapid and exothermic, and often occur barrierlessly to various carbon sites of unsaturated hydrocarbons. The reaction mechanism is proposed on the basis of the resulting potential energy surface (PES) that includes all the possible intermediates and transition states that can occur, and all the products that lie on the surface. The B3LYP/6-311g(d,p) level of theory is employed to determine optimized electronic structures, moments of inertia, vibrational frequencies, and zero-point energy. They are followed by single point higher-level CCSD(T)/cc-vtz calculations, including extrapolations to complete basis sets (CBS) of the reactants and products. A microcanonical RRKM study predicts single-collision (zero-pressure limit) rate constants of all reaction paths on the potential energy surface, which is then used to compute the branching ratios of the products that result. These theoretical calculations are conducted either jointly or in parallel to experimental work to elucidate the chemical composition of Titan’s atmosphere, the ISM, and cold celestial bodies.

The role of disulfide bridges in the 3-D structures of the antimicrobial peptides gomesin and protegrin-1: a molecular dynamics study

Some antimicrobial peptides have a broad spectrum of action against many different kinds of microorganisms. Gomesin and protegrin-1 are examples of such antimicrobial peptides, and they were studied by molecular dynamics in this research. Both have a beta-hairpin conformation stabilized by two disulfide bridges and are active against Gram-positive and Gram-negative bacteria, as well as fungi. In this study, the role of the disulfide bridge in the maintenance of the tertiary peptide structure of protegrin-1 and gomesin is analyzed by the structural characteristics of these peptides and two of their respective variants, gomy4 and proty4, in which the four cysteines are replaced by four tyrosine residues. The absence of disulfide bridges in gomy4 and proty4 is compensated by overall reinforcement of the original hydrogen bonds and extra attractive interactions between the aromatic rings of the tyrosine residues. The net effects on the variants with respect to the corresponding natural peptides are: i) maintenance of the original beta-hairpin conformation, with great structural similarities between the mutant and the corresponding natural peptide; ii) combination of positive F and. Ramachandran angles within the hairpin head region with a qualitative change to a combination of positive (F) and negative (.) angles, and iii) significant increase in structural flexibility. Experimental facts about the antimicrobial activity of the gomesin and protegrin-1 variants have also been established here, in the hope that the detailed data provided in the present study may be useful for understanding the mechanism of action of these peptides.

Study of the antimicrobial peptide indolicidin and a mutant in micelle medium by molecular dynamics simulation

The antimicrobial peptide indolicidin (IND) and the mutant CP10A in hydrated micelles were studied using molecular dynamics simulations in order to observe whether the molecular dynamics and experimental data could be sufficiently correlated and a detailed description of the interaction of the antimicrobial peptides with a model of the membrane provided by a hydrated micelle system could be obtained. In agreement with the experiments, the simulations showed that the peptides are located near the surface of the micelles. Peptide insertions agree with available experimental data, showing deeper insertion of the mutant compared with the peptide IND. Major insertion into the hydrophobic core of the micelle by all tryptophan and mutated residues of CP10A in relation to IND was observed. The charged residues of the terminus regions of both peptides present similar behavior, indicating that the major differences in the interactions with the micelles of the peptides IND and CP10A occur in the case of the hydrophobic residues.

Dynamic polarizability, Cauchy moments, and the optical absorption spectrum of liquid water: A sequential molecular dynamics/quantum mechanical approach

The dynamic polarizability and optical absorption spectrum of liquid water in the 6-15 eV energy range are investigated by a sequential molecular dynamics (MD)/quantum mechanical approach. The MD simulations are based on a polarizable model for liquid water. Calculation of electronic properties relies on time-dependent density functional and equation-of-motion coupled-cluster theories. Results for the dynamic polarizability, Cauchy moments, S(-2), S(-4), S(-6), and dielectric properties of liquid water are reported. The theoretical predictions for the optical absorption spectrum of liquid water are in good agreement with experimental information.

Electronic properties of liquid ammonia: A sequential molecular dynamics/quantum mechanics approach

The electronic properties of liquid ammonia are investigated by a sequential molecular dynamics/quantum mechanics approach. Quantum mechanics calculations for the liquid phase are based on a reparametrized hybrid exchange-correlation functional that reproduces the electronic properties of ammonia clusters [(NH(3))(n); n=1-5]. For these small clusters, electron binding energies based on Green's function or electron propagator theory, coupled cluster with single, double, and perturbative triple excitations, and density functional theory (DFT) are compared. Reparametrized DFT results for the dipole moment, electron binding energies, and electronic density of states of liquid ammonia are reported. The calculated average dipole moment of liquid ammonia (2.05 +/- 0.09 D) corresponds to an increase of 27% compared to the gas phase value and it is 0.23 D above a prediction based on a polarizable model of liquid ammonia [Deng , J. Chem. Phys. 100, 7590 (1994)]. Our estimate for the ionization potential of liquid ammonia is 9.74 +/- 0.73 eV, which is approximately 1.0 eV below the gas phase value for the isolated molecule. The theoretical vertical electron affinity of liquid ammonia is predicted as 0.16 +/- 0.22 eV, in good agreement with the experimental result for the location of the bottom of the conduction band (-V(0)=0.2 eV). Vertical ionization potentials and electron affinities correlate with the total dipole moment of ammonia aggregates. (c) 2008 American Institute of Physics.

Mn dimers on graphene nanoribbons: An ab initio study

We show, using ab initio density functional theory calculations, that Mn dimers adsorbed on graphene nanoribbons (Mn(2)/GNRs) present a magnetic bistability, as does the isolated Mn dimer. Our total energy results indicate that Mn dimers lying along the edge sites of zigzag GNRs represent the most likely configuration. We find that similar to the isolated Mn(2) molecule, the antiferromagnetic coupling represents the ground state for Mn(2)/GNR, and the spin density configuration of the GNR does not play an important role on the net magnetic moment of Mn(2), which makes GNRs an ideal substrate for adsorption of these molecules. The ground state and the excited state configuration of the Mn dimer, viz., low-spin (LS) and high-spin (HS), are maintained in the face of changes in the spin density configuration of the substrate. Here we find that the Mn(2)/GNR systems exhibit a LS <-> HS binary behavior, which can be considered as a useful property in the development of nanomemories based upon metallic clusters. (C) 2011 American Institute of Physics. [doi:10.1063/1.3553849]

Bilayer graphene dual-gate nanodevice: An ab initio simulation

We study the electronic transport properties of a dual-gated bilayer graphene nanodevice via first-principles calculations. We investigate the electric current as a function of gate length and temperature. Under the action of an external electrical field we show that even for gate lengths up 100 angstrom, a nonzero current is exhibited. The results can be explained by the presence of a tunneling regime due the remanescent states in the gap. We also discuss the conditions to reach the charge neutrality point in a system free of defects and extrinsic carrier doping.

Accurate prediction of the Si/SiO(2) interface band offset using the self-consistent ab initio DFT/LDA-1/2 method

We use the density functional theory/local-density approximation (DFT/LDA)-1/2 method [L. G. Ferreira , Phys. Rev. B 78, 125116 (2008)], which attempts to fix the electron self-energy deficiency of DFT/LDA by half-ionizing the whole Bloch band of the crystal, to calculate the band offsets of two Si/SiO(2) interface models. Our results are similar to those obtained with a ""state-of-the-art"" GW approach [R. Shaltaf , Phys. Rev. Lett. 100, 186401 (2008)], with the advantage of being as computationally inexpensive as the usual DFT/LDA. Our band gap and band offset predictions are in excellent agreement with experiments.

Edge effects in bilayer graphene nanoribbons: Ab initio total-energy density functional theory calculations

We show that the ground state of zigzag bilayer graphene nanoribbons is nonmagnetic. It also possesses a finite gap, which has a nonmonotonic dependence with the width as a consequence of the competition between bulk and strongly attractive edge interactions. All results were obtained using ab initio total-energy density functional theory calculations with the inclusion of parametrized van der Waals interactions.

Spin filtering and disorder-induced magnetoresistance in carbon nanotubes: Ab initio calculations

Nitrogen-doped carbon nanotubes can provide reactive sites on the porphyrin-like defects. It is well known that many porphyrins have transition-metal atoms, and we have explored transition-metal atoms bonded to those porphyrin-like defects inN-doped carbon nanotubes. The electronic structure and transport are analyzed by means of a combination of density functional theory and recursive Green's function methods. The results determined the heme B-like defect (an iron atom bonded to four nitrogens) is the most stable and has a higher polarization current for a single defect. With randomly positioned heme B defects in nanotubes a few hundred nanometers long, the polarization reaches near 100%, meaning they are effective spin filters. A disorder-induced magnetoresistance effect is also observed in those long nanotubes, and values as high as 20 000% are calculated with nonmagnectic eletrodes.

Barrier-free substitutional doping of graphene sheets with boron atoms: Ab initio calculations

Using ab initio methods, we propose a simple and effective way to substitutionally dope graphene sheets with boron. The method consists of selectively exposing each side of the graphene sheet to different elements. We first expose one side of the membrane to boron while the other side is exposed to nitrogen. Proceeding this way, the B atoms will be spontaneously incorporated into the graphene membrane without any activation barrier. In a second step, the system should be exposed to a H-rich environment, which will remove the CN radical from the layer and form HCN, leading to a perfect substitutional doping.

Polarization effects in molecular dynamics simulations of glass-formers Ca(NO(3))(2)center dot nH(2)O, n=4, 6, and 8

Thermodynamics, equilibrium structure, and dynamics of glass-forming liquids Ca(NO(3))(2)center dot nH(2)O, n=4, 6, and 8, have been investigated by molecular dynamics (MD) simulations. A polarizable model was considered for H(2)O and NO(3)- on the basis of previous fluctuating charge models for pure water and the molten salt 2Ca(NO(3))(2)center dot 3KNO(3). Similar thermodynamic properties have been obtained with nonpolarizable and polarizable models. The glass transition temperature, T(g), estimated from MD simulations was dependent on polarization, in particular the dependence of T(g) with electrolyte concentration. Significant polarization effects on equilibrium structure were observed in cation-cation, cation-anion, and water-water structures. Polarization increases the diffusion coefficient of H(2)O, but does not change significantly the diffusion coefficients of ions. Viscosity decreases upon inclusion of polarization, but the conductivity calculated with the polarizable model is smaller than the nonpolarizable model because polarization enhances anion-cation interactions.

Fragility and glassy dynamics of 2Ca(NO(3))(2)center dot 3KNO(3) under pressure: Molecular dynamics simulations

Molecular dynamics simulations of the glass-forming liquid 2Ca(NO(3))(2)center dot 3KNO(3) (CKN) were performed from high temperature liquid states down to low temperature glassy states at six different pressures from 10(-4) to 5.0 GPa. The temperature dependence of the structural relaxation time indicates that the fragility of liquid CKN changes with pressure. In line with recent proposal [Scopigno , Science 302, 849 (2003)], the change on liquid fragility is followed by a proportional change of the nonergodicity factor of the corresponding glass at low temperature. (c) 2008 American Institute of Physics.

Molecular dynamics, density functional, ADMET predictions, virtual screening, and molecular interaction field studies for identification and evaluation of novel potential CDK2 inhibitors in cancer therapy

In this work, we have used molecular dynamics, density functional theory, virtual screening, ADMET predictions, and molecular interaction field studies to design and propose eight novel potential inhibitors of CDK2. The eight molecules proposed showed interesting structural characteristics that are required for inhibiting the CDK2 activity and show potential as drug candidates for the treatment of cancer. The parameters related to the Rule of Five were calculated, and only one of the molecules violated more than one parameter. One of the proposals and one of the drug-like compounds selected by virtual screening indicated to be promising candidates for CDK2-based cancer therapy.