Probing the effect of nitro groups in nitramine based energetic molecules: a DFT and charge density study

The structure, electron density distribution, energetic and electrostatic properties of simple nitramine based energetic TMA, DMNA, MDA and TNA molecules were determined using density functional theory (B3LYP) with the 6-311G^** and aug-cc-pVDZ basis sets coupled with Bader's theory of atoms in molecules. In the NO₂ group substituted molecules, the N–N bond distance increases with the increase of NO₂ groups, whereas in C–N bonds, this effect is relatively less, and the distances are almost equal. The topological analysis of electron density reveals that the electron density ρ_bcp(r) of C–N and N–N bonds are significantly decreasing with the increase of NO₂ groups in the nitramine molecules. The Laplacian of electron density ▽²ρ_bcp(r) of N–NO₂ bonds [DMNA: ? 16.7 eÅ^? 5, MDA: ? 12.8 eÅ^? 5 and TNA: ? 7.9 eÅ^? 5] of the molecules are relatively less negative, and the values also decrease with the increase of NO₂ groups; this implies that the charge concentration decreases with the increase of NO₂ groups, which leads to weakening the N–N bonds of the molecules. The isosurface of molecular electrostatic potential displays high electronegative regions around the NO₂ groups. The oxygen balance OB₁₀₀ of the molecules increases as the number of NO₂ group increases in the molecules, in which, the TNA molecule having maximum OB₁₀₀ value [+7.89]. The band gap, heat of detonation, bond dissociation energy and charge imbalance are predominantly depends on the number of NO₂ group present in the molecule. The charge imbalance parameter (ν) has been calculated for all molecules, which reveals that TNA is a highly sensitive molecule, the corresponding ν value is 0.047.     

Comparison of the structures and electronic properties of sodalites containing alkali metals and alkali-earth metals and their hydrates

The effects of different alkali and alkali-earth metal ions on the electronic structures and properties of sodalite M_n[AlSiO₄]₆ (M-SOD) and their hydrates M_n[AlSiO₄]₆?8H₂O (M=Li, Na, K, n = 6; M=Ca, n = 3) were studied using density functional theory method. Theoretical calculations predicted that the Al–O–Si bond angle and cation-framework oxide distance in sodalites with alkali metal cations are correlated with cell volumes. The reduced bandwidths in M-SOD (M=Li, Na and K) show that the inter-atomic orbital overlap in sodalites is weaker than those in the hydrate phases. Frontier molecular orbital analysis indicated that oxygen atoms in the frameworks and most metal ions of SOD and their corresponding hydrates exhibit high reactivity. The interactions existing in sodalites and hydrates were qualitative described. The calculated combination energies of metal ions with framework of sodalites are in the order of K⁺< Na⁺< Li⁺< Ca²⁺. This finding confirms the experimental observation for ion exchange.     

Density functional theory and ab initio studies of vibrational spectra of 2-bis (2-chloroethyl) aminoperhydro-1,3,2-oxazaphosphorinane-2-oxide

The Fourier transform Raman (FTR) and Fourier transform infrared (FTIR) spectra of 2-bis (2-chloroethyl) aminoperhydro-1,3,2-oxazaphosphorinane-2-oxide were recorded in the regions 4000–100 cm^? 1 and 4000–400 cm¹, respectively, in the solid phase. Molecular electronic energy, geometrical structure, harmonic vibrational spectra, infrared intensities and Raman scattering activities, highest occupied molecular orbital, lowest unoccupied molecular orbital energy, energy gaps and thermodynamical properties such as zero-point vibrational energies, rotational constants, entropies and dipole moment were computed at the Hartree–Fock/6-31G(d,p) and three parameter hybrid functional Lee–Yang–Parr/6-31G(d,p) levels of theory. The vibrational studies were interpreted in terms of potential energy distribution. The results were compared with experimental values with the help of scaling procedures. The observed wave number in FTIR and FTR spectra was analysed and assigned to different normal modes of the molecule. Most of the modes have wave numbers in the expected range and are in good agreement with computed values.     

Molecular structure and spectroscopic (FT-IR,FT-Raman, 13C, 1H NMR and UV) studies of 3,4-dihydroxy-l-phenylalanine using density functional theory

Fourier transform infrared (FT-IR) and Fourier transform Raman (FT-Raman) spectra of 3,4-dihydroxy-l-phenylalanine (3,4-DPA) in solid phase were recorded and analysed in this research. Along with this, the IR spectra in CHCl₃ and the use of acetone as solvents of 3,4-DPA were also recorded. The equilibrium geometry, bonding features and harmonic vibrational frequencies were investigated with the help of density functional theory (DFT) method. The ¹H and ¹³C nuclear magnetic resonance chemical shifts of the molecule were calculated using the gauge including atomic orbital method and compared with experimental results. Stability of the molecule arising from hyperconjugative interactions and charge delocalisation was analysed using natural bond orbital analysis. The results show that charge in electron density (E _D) in the σ* and π* antibonding orbitals and second-order delocalisation energies E(2) confirms the occurrence of intramolecular charge transfer within the molecule. UV–vis spectrum of the compound was recorded and the electronic properties, such as HOMO and LUMO energies, were analysed using the time-dependent (TD)-DFT approach. Finally, the calculation results were applied to simulate infrared and Raman spectra of the title compound, which showed good agreement with the observed spectra.     

Ab initio and molecular dynamics studies of solid β-HMX: effects of hydrostatic pressure and high temperature

Using first-principles density functional theory and classical molecular dynamics (MD), the structural, electronic and mechanical properties of the energetic material β-HMX have been studied. The crystal structure optimised by the local density approximation calculations compares reasonably with the experimental data. Electronic band structure and density of states indicate that β-HMX is an insulator with a band gap of 3.059 eV. The pressure effect on the crystal structure and physical properties has been investigated in the range of 0–40 GPa. The crystal structure and electronic properties change slightly as the pressure increases from 0 to 2.5 GPa; when the pressure is above 2.5 GPa, further increment of the pressure results in significant changes in crystal structure. There is a larger compression along the b-axis than along the a- and c-axes. Isothermal–isobaric MD simulations on β-HMX were performed in the temperature range of 5–400 K. Phase transition at 360 K, corresponding to a volume interrupt, was found. The computed thermal expansion coefficients show anisotropic behaviour with a slightly larger expansion along the b- and c-axes than along the a-axis. In the temperature range of 5–360 K, β-HMX possesses good plasticity and its stiffness decreases with increasing the temperature.     

Large variation in electrostatic contours upon addition of steric parameters and the effect of charge calculation schemes in CoMFA on mutagenicity of MX analogues

Partial atomic charge is a useful concept to describe physico-chemical properties of a molecule. For this, various schemes have been devised to get reasonable values. Mutagen X is an ideal set to test the effect of partial atomic charge variation. Therefore, we collected data from previous reports and studied various charge schemes. Our systematic study covers 26 charge calculation schemes along with a broad range of levels of theory. Charge calculation schemes include charges from charge equalisation, electrostatic potential fitting, molecular orbital and atomic polar tensor. Calculation levels span from empirical, semi-empirical, Hartree–Fock, density functional and Møller–Plesset 2. We also used two validation statistics for internal prediction. To observe the electrostatic effect accurately during comparative molecular field analysis (CoMFA) modelling, we first studied isolated electrostatic parameters to avoid interaction effect with steric parameters. The results clearly show that adding steric parameters does change statistical conclusions as well as CoMFA maps. Although there was a weak trend that quantum mechanical (QM)-derived charges gave better statistical values, it is not apparent statistically (alpha = 0.05). Particularly, Mülliken population analysis (MPA) did not produce better results. Therefore, when we excluded MPA schemes from QM calculation, the QM-derived charges were found to be significant, i.e. sophisticated charge schemes other than MPA with QM methods were found to be superior to simple empirical charge schemes. In addition, we demonstrated that in order to test charge schemes properly, excluding steric parameter is more important. This work exemplifies Occam's theorem of parsimony. A simpler model is a better model.     

Intramolecular proton transfer in adenine imino tautomers

The electronic structural impact on intramolecular proton transfer in the cis- and trans-imino N7 and N9 tautomers of adenine (A) has been studied quantum mechanically, using density functional theory (B3LYP/TZVP, SAOP/TZ2P, LB94/TZ2P) and Green function (OVGF/TZVP) models. It is found that proton transfer does not significantly change isotropic properties but has profound impact on electron distributions of the species through anisotropic properties. The relative energies with respect to the canonical A tautomer (amino-9H), ΔE, for imino 7Hcis, imino 7Htrans, imino 9Hcis and imino 9Htrans are calculated as 16.15, 16.43, 18.46 and 13.80 kcal mol^? 1 (B3LYP/TZVP model) and some minor changes in perimeters of the purine ring is also observed. The Hirshfeld atomic charges indicate that whether a proton attached to N₍₇₎ or N₍₉₎ causes a significant local charge redistribution. However, these charges are insensitive to cis–trans proton transfer. Condensed Fukui function reveals N₍₁₀₎ and C₍₈₎ as the most electrophilic reactive site among N- and C-atom sites, respectively. We also found that proton transfer significantly alters in-plane σ orbitals, rather than out of plane π orbitals including the frontier orbital 6a″. Moreover, orbital based responses to various proton transfers are presented: the orbital 29a′ (HOMO-1) is a signature orbital differentiating all the four tautomers. Orbital 27a′ is a site (N₍₇₎ and N₍₉₎) sensitive orbital, whereas orbital 22a′ is only sensitive to proton orientation on the imino group = N–H.     

Electronic structure simulations of 2,6-dimethyl-2,5-heptadien-4-one by FTIR,FT-Raman,NMR, UV–vis,NBO and density functional theory

Density functional theory (DFT) (B3LYP and B3PW91) calculations have been carried out for 2,6-dimethyl-2,5-heptadien-4-one (DMHD4O) using 6–311++ G** basis set. Complete vibrational assignment and analysis of the fundamental modes of the compound were carried out from the FTIR and FT-Raman spectral data. The theoretical electronic absorption has been calculated by using time-dependent DFT (TD-DFT) methods and compared with the experimental spectra. The theoretically computed Frontier energy gaps and TD-DFT calculations are in good agreement with the experimental UV–vis spectral absorption. The chemical hardness measured from the Frontier molecular orbital energies of DMHD4O is 0.0693 eV. Electronic stability of the compound arising from hyperconjugative interactions and charge delocalisation were also investigated based on the natural bond orbital (NBO) analysis. Effective stabilisation energy E ⁽²⁾ associated with the interactions of the π and the lone pair of electrons was determined by the NBO analysis. ¹³C and ¹H NMR chemical shifts of the compound have been calculated by means of Gauge-Invariant Atomic Orbital using B3LYP/6–311++ G** method. The partial ionic character of the carbonyl group due to resonance render a partially positive charge to the carbonyl carbon, and thus C4 chemical shift lie in the very downfield 191.6 ppm. Comparison between the experimental and the theoretical results indicates that B3LYP method is able to provide satisfactory results for predicting vibrational, electronic and NMR properties.     

Molecular dynamics study of methane in water: diffusion and structure

In this study, theoretical analysis on the geometries and electronic properties of various conjugated oligomers based on thiophene (Th) or bicyclic non-classical Th units is reported. The dihedral angle, bond length, bond-length alternation, bond critical point (BCP) properties, nucleus-independent chemical shift (NICS) and Wiberg bond index (WBIs) are analysed and correlated with conduction properties. The changes of bond length, BCP properties, NICS and WBIs all show that the conjugational degree is increased systematically with main chain extension. As a result, the highest occupied molecular orbital–lowest unoccupied molecular orbital energy separation (E _g) is decreased upon chain elongation. The E _g of oligomers based on bicyclic non-classical Th unit is much lower than that of Th-based oligomers due to the narrower E _g of bicyclic non-classical Ths compared to Th, which indicates that the narrow E _g of the bicyclic non-classical Ths can be carried over to their polymers by using them as building blocks for the polymers. The band structures and density of states analysis show that the four polymers all have small band gap ( < 0.9 eV), wide highest occupied bandwidth and relatively small effective mass of hole, which indicate that those proposed polymers may be potential conductors.     

Spectroscopic investigation on glutamic acid by Coulomb-attenuating and double hybrid density functional theory methods

This study deals with the identification of glutamic acid by means of quantum chemical approach. FT-IR, FT-Raman and UV–vis spectra were recorded in the region 4000–400, 4000–50 cm^? 1 and 200–600 nm, respectively. CAM-B3LYP/6-31G(d,p) and B2PLYP/6-31G(d,p) calculations were performed to obtain the optimised molecular structures, vibrational frequencies and corresponding vibrational assignment, thermodynamic properties and natural bonding orbital (NBO) analysis. The results show that the obtained optimised geometric parameters (bond lengths, bond angles and bond dihedrals) and vibrational frequencies were found to be in good agreement with the experimental results. The calculations of the electronic spectra were compared with the experimental ones. Furthermore, highest occupied molecular orbital and lowest unoccupied molecular orbital analyses and UV–vis spectral analysis were also performed to determine the energy band gaps and transition states. NBO analysis, calculated using density functional theory methods (CAM-B3LYP/6-31G(d,p) and B2PLYP/6-31G(d,p)), was induced to find inter-molecular atoms. ¹³C and ¹H NMR isotropic chemical shifts were calculated and the assignments made were compared with the ChemDraw Ultra values.     

Electronic structure and band gap engineering of ZnO-based semiconductor alloy films

We have conducted first-principles total-energy density functional calculations to study the atomic structures, band structures and electronic structures of Zn_1 ? xM_xO (M = Be, Mg, Cd, Ag, Cu) semiconductor alloys. The Heyd–Scuseria–Ernzerhof hybrid functional has been performed to yield lattice constants and band gaps of Zn_1 ? xM_xO semiconductor in much better agreement with experimental data than with the standard local exchange correlation functional. We found that the strong coupling between O 2p and Cu 3d or Ag 4d bands plays a key role in narrowing of band gaps and leading to the half-metallic behaviour interpreted with the unique spatial distribution pattern between the highest occupied molecular orbital and the lowest unoccupied molecular orbital.     

Molecular dynamics and quantum chemical studies of solvent effects on cyclo glycylglycine and glycylalanine dipeptides

Hydrogen bond (H-bond) interactions between the two cyclo dipeptides, cyclo(glycyl-glycine) (CGG) and cyclo(glycyl-alanine) (CGA), and water have been studied using molecular dynamics (MD) and quantum chemical methods. The MD studies have been carried out on CGG and CGA in water using fixed charge force field (AMBER ff03) for over 10 ns with a MD time step of 2 fs. The results of this study show that the solvation pattern influences the conformations of the cyclo dipeptides. Following molecular simulations, post Hartree–Fock and density functional theory methods have been used to explore the molecular properties of the cyclo dipeptides in gaseous and aqueous phase environments. The self-consistent reaction field theory has been used to optimise the cyclopeptides in diethyl ether (? = 4.3) and water (? = 78.5), and the solvent effects have been analysed. A cluster of eight water molecules leads to the formation of first solvation shell of CGG and CGA and the strong H-bonding mainly contributes to the interaction energies. The H-bond interactions have been analysed by the calculation of electron density ρ(r) and its Laplacian ▽²ρ(r) at bond critical points using atoms in molecules theory. The natural bond orbital analysis was carried out to reveal the nature of H-bond interactions. In the solvated complexes, the keto carbons registered the maximum NMR chemical shifts.     

Shielding effect of a PEG molecule of a mono-PEGylated peptide varies with PEG chain length

Abstract

‘Shielding’ effect of a conjugated PEG molecule could cause a change in the electrostatic interaction characteristics of a PEGylate. We investigated how PEG chain length (or molecular weight) alters the electrostatic interaction potential of exenatide variants using their mono-PEGylates in a branched and linear form as model PEGylates. First, we performed the experiments to demonstrate the elution time changes of the mono-PEGylates conjugated with various MW PEGs (5, 10, 20, and 40 kD) using cation exchange chromatography (HiTrap^® SP) at various pHs (2.5, 3.0, 3.5, and 4.0). Then, we calculated the net surface charge of each mono-PEGylate to propose the PEG molecule’s shielding range in terms of the number of amino acids adjacent to the conjugation residue, assuming that a PEG molecule in solution sweeps out a spherical space and an exenatide molecule have a secondary structure. The net charge calculation result was well-correlated with the experimental elution time data, where 5, 10, 20, and 40 kD PEG hindered the electrostatic potential of 5, 8, 12, and 17 amino acid residues in maximum, respectively, on each side of the conjugation point.     

Ab initio investigations on planar (MgO)n clusters (n = 1–5) and their hydrogen adsorption behaviour

Ab initio calculations (B3LYP and PBE-D3) of the structures, stabilities, vibrational, electronic and hydrogen adsorption behaviour of (MgO)_n clusters are performed using 6-311+ + G(d,p) basis set. The planar (MgO)_n clusters are found to be global minima for n ≤ 3 and local minima for n = 4 and 5. In addition, we have also analysed global minimum structures of (MgO)₄ and (MgO)₅. The binding energies suggest that their stabilities increase successively. Vibrational frequencies and IR intensities further support the enhanced stability with an increase in the size of (MgO)_n clusters. The highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) surfaces are used to explain and discuss the electronic properties. Finally, we have demonstrated hydrogen storage capacity of (MgO)_n clusters, considering hydrogen adsorption on planar as well as global minimum (MgO)₄ and (MgO)₅ clusters. We have noticed that four and five H₂ molecules can be easily adsorbed by (MgO)₄ and (MgO)₅ clusters having adsorption energy of 0.13–0.14 eV with mass ratio of 4.76%. Thus, the present study is expected to motivate further the applications of small clusters for efficient hydrogen energy storage.     

Modeling the interaction between anti-cancer drug penicillamine and pristine and functionalized carbon nanotubes for medical applications: density functional theory investigation and a molecular dynamics simulation

Abstract

The present study focuses on the prediction and investigation of binding properties of penicillamine with pure (5,5) single-walled carbon nanotube (SWCNT) and functionalized SWCNT (f-SWCNT) through the B3LYP and M06-2X functionals using the 6-31G** basis set. The electronic and structural properties, adsorption energy and frontier molecular orbitals of various configurations are examined. Our theoretical results indicated that the interaction of the nanotubes with penicillamine molecule is weak so that the drug adsorption process is typically physisorption. Also, results of theoretical calculations show that the adsorption of the drug molecule on f-SWCNT is stronger with shorter intermolecular distances in comparison to pure SWCNT. The natural bond orbital (NBO) analysis of studied systems demonstrates that the charge is transferred from penicillamine molecule to the nanotubes. Furthermore, molecular dynamics (MD) simulation is employed to evaluate the dynamic and diffusion behavior of drug molecule on SWCNT and f-SWCNT. Energy results show that drug molecule spontaneously moves toward the carriers, and the van der Waals energy contributions in drug adsorption are more than electrostatic interaction. The obtained results from MD simulation confirm that the functionalization of SWCNT leads to increase in the solubility of the carrier in aqueous solution.

Communicated by Ramaswamy H. Sarma     

Complexes of DNA bases and Watson–Crick base pairs interaction with neutral silver Agn (n = 8, 10, 12) clusters: a DFT and TDDFT study

We study the binding of the neutral Ag_n (n = 8, 10, 12) to the DNA base-adenine (A), guanine (G) and Watson–Crick –adenine-thymine, guanine-cytosine pairs. Geometries of complexes were optimized at the DFT level using the hybrid B3LYP functional. LANL2DZ effective core potential was used for silver and 6–31 + G^** was used for all other atoms. NBO charges were analyzed using the Natural population analysis. The absorption properties of Ag_n–A,G/WC complexes were also studied using time-dependent density functional theory. The absorption spectra for these complexes show wavelength in the visible region. It was revealed that silver clusters interact more strongly with WC pairs than with isolated DNA complexes. Furthermore, it was found that the electronic charge transferred from silver to isolated DNA clusters are less than the electronic charge transferred from silver to the Ag_n–WC complexes. The vertical ionization potential, vertical electron affinity, hardness, and electrophilicity index of Ag_n–DNA/WC complexes have also been discussed.     

A DFT study on the second-order nonlinear optical properties of the plenary mixed-metal polyoxometalate

The redox and second-order nonlinear optical (NLO) properties of organic–inorganic hybrid polyoxoanions [LNbOEMe₃]^3 ? , [LNbOEPh₃]^3 ?  (L = α-{PW₁₁O₃₉}^7 ? , E = Si, Ge, Sn, Pb) are investigated by density functional theory method. The element substitution effects on the molecular nonlinear response have been analysed. The results show that the computed β₀ values depend on both E (E = Si < Ge < Sn < Pb) and the organic groups connected with E. For [LNbOEPh₃]^3 ?  (L = α-{PW₁₁O₃₉}^7 ? , E = Si, Ge, Sn, Pb), the analysis of major contributions to β₀ value suggests that the charge transfer from organic group to Keggin polyanion plays the key role in NLO response; the polyanion acts as a donor, whereas the organic group acts as an acceptor in [LNbOEPh₃]^3 ? . Our results show that this kind of organic–inorganic hybrid compound possesses larger molecular second-order polarisability and might be a potential NLO material.     

Spectroscopic (FTIR,FT-Raman,NMR and UV) and molecular structure investigations of 1,5-diphenylpenta-2,4-dien-1-one: a combined experimental and theoretical approach

The title molecule 1,5-diphenylpenta-2,4-dien-1-one (cinnamylideneacetophenone, CA) has been synthesised and characterised by FTIR, FT-Raman, NMR and UV–vis spectral analyses. The possible stable conformers of the CA molecule were searched by potential energy surface scan at B3LYP level of theory. The molecular geometry from X-ray determination of the CA molecule in the ground state has been compared using the density functional theory (DFT) with 6-31G(d,p) basis set. The harmonic vibrational modes, the corresponding wavenumbers and IR and Raman intensities of most stable conformer were calculated by the DFT method. The assignments of the fundamentals were proposed on the basis of total energy distribution calculations. The calculated ¹³C and ¹H NMR chemical shifts using gauge including atomic orbitals approach are in good agreement with the observed chemical shifts. The molecular stability and bond strength have been investigated by applying natural bond orbital analysis. Using the time-dependent DFT method, the electronic absorption spectrum of the title compound has been predicted and the electronic transitions within the molecule have been interpreted. The molecular electrostatic potential map was used for predicting possible hydrogen and oxygen bonding sites in the CA molecule.