Theoretical evaluation of flotation performance of carboxyl hydroxamic acids with different number of polar groups on the surfaces of diaspore (010) and kaolinite (001)

The adsorption behaviors of three carboxyl hydroxamic acids on diaspore (010) and kaolinite (001) have been studied by density functional theory (DFT) and molecular dynamics (MD) method. The results indicated that carboxyl hydroxamic acids could adsorb on diaspore surface by ionic bonds and hydrogen bonds, and adsorb on kaolinite surface by hydrogen bonds. The models of carboxyl hydroxamic acids adsorbed on diaspore and kaolinite surfaces are proposed.

Molecular dynamics simulation on miscibility of trans-1,4,5,8-tetranitro-1,4,5,8 -tetraazadecalin (TNAD) with some propellants

The solubility parameters of TNAD, HMX, RDX, DINA, DNP propellants were predicted by molecular dynamics (MD) simulation in order to evaluate the miscibility of TNAD and the other four propellants. The results show that the order of miscibility is TNAD/DINA > TNAD/DNP > TNAD/RDX > TNAD/HMX from the analysis of miscibility. The densities and binding energies of TNAD/propellants blends were further investigated. The results indicate that the better the miscibility between TNAD and the propellants, the smaller the variation of the density rate. The larger the intermolecular interaction, the better the miscibility between components. The analysis of radial distribution function shows that the main interaction ways between TNAD and other energetic components are short-range interactions.

Relationship between structure and entropy contributions in an anthraquinone mercapto derivative

The structural and thermodynamic properties of an anthraquinone derivative were studied by means of quantum-chemical calculations. Conformational analysis using ab initio and density functional theory methods revealed 14 low-energy conformers. In order to discuss similarities and differences in entropy of the conformers, the rotational and vibrational contributions to entropy were correlated with changes in conformer structure. The component of the moment of inertia perpendicular to the molecular plane gives significant input to ΔS _rot , whereas the largest contributions to the ΔS _vib have vibrations associated with the τ _S1C20 coordinate.

Heavy periodane

The potential energy surface of the hypothetical NaMgAlSiPSCl system (heavy periodane) is exhaustively analyzed via the gradient embedded genetic algorithm (GEGA) in combination with density functional theory (DFT) computations. The electronegativity differences among the elements in both the second and third rows of the periodic table indicate that low-energy heavy periodane structures are obtained when highly electronegative and electropositive elements are bound together, but the global minimum of the heavy periodane system is completely different to its second-row analog (LiBeBCNOF).

Natures of benzene-water and pyrrole-water interactions in the forms of σ and π types: theoretical studies from clusters to liquid mixture

A combined and sequential use of quantum mechanical (QM) calculations and classical molecular dynamics (MD) simulations was made to investigate the σ and π types of hydrogen bond (HB) in benzene-water and pyrrole-water as clusters and as their liquid mixture, respectively. This paper aims at analyzing similarities and differences of these HBs resulted from QM and MD on an equal footing. Based on the optimized geometry at ωb97xD/aug-cc-pVTZ level of theory, the nature and property of σ and π types of HBs are unveiled by means of atoms in molecules (AIM), natural bond orbital (NBO) and energy decomposition analysis (EDA). In light of the above findings, MD simulation with OPLS-AA and SPC model was applied to study the liquid mixture at different temperatures. The MD results further characterize the behavior and structural properties of σ and π types HBs, which are somewhat different but reasonable for the clusters by QM. Finally, we provide a reasonable explanation for the different solubility between benzene/water and pyrrole/water.

Simulations on the possibility of formation of complexes between fluorouracil drug and cucurbit[n]urils: ab initio van der Waals DFT study

The binding geometry of fluorouracil/cucurbit[n]urils (CB[n]s) complexes with n?=?5–8 is investigated using the first-principles van der Waals density functional (vdW-DF) method, involving full geometry optimization. Such host-guest complexes are typically calculated using conventional DFT method, which significantly underestimates non-local dispersion forces (or vdW contributions) and therefore affects interactions between respected entities. We address here the role of vdW forces for the fluorouracil and CB[n]s molecules which can form directional hydrogen bonds with each other. It was found that the inclusion of dispersion interactions significantly affects the host-guest binding properties and the hydrogen bonding between the molecules provides the main binding mechanism, while results in the same geometries for the considered complexes. The 0.84 eV binding energy, for the thermodynamically favorable state, reveals that the interaction of fluorouracil with CB[n]s is an exothermic interaction and typical for strong hydrogen bonding. Furthermore, we have investigated the binding nature of these host-guest systems in aqueous solution with ab initio MD simulations adopting vdW-DF method. These findings afford evidence for the applicability of the vdW-DF approach and provide a realistic benchmark for the investigation of the host-guest complexes.

Characterization of the structures and dynamics of phosphoric acid doped benzimidazole mixtures: a molecular dynamics study

Benzimidazole-based polymer membranes like poly(2,5-benzimidazole) (ABPBI) doped with phosphoric acid (PA) are electrolytes that exhibit high proton conductivity in fuel cells at elevated temperatures. The benzimidazole (BI) moiety is an important constituent of these membranes, so the present work was performed in order to achieve a molecular understanding of the BI–PA interactions in the presence of varying levels of the PA dopant, using classical molecular dynamics (MD) simulations. The various hydrogen-bonding interactions, as characterized based on structural properties and hydrogen-bond lifetime calculations, show that both BI and PA molecules exhibit dual proton-acceptor/donor functionality. An examination of diffusion coefficients showed that the diffusion of BI decreases with increasing PA uptake, whereas the diffusion of PA slightly increases. The hydrogen-bond lifetime calculations pointed to the existence of competitive hydrogen bonding between various sites in BI and PA.

Molecular simulation study of PAMAM dendrimer composite membranes

Pure polysulfone (PSF) and its composites with chitosan (CST), hyaluronic acid (HA), conventional poly(amidoamine), and hydroxyl poly(amidoamine) dendrimers as the membranes for separation of the gases, methane, carbon dioxide, hydrogen sulfide, nitrogen, and oxygen have been studied by molecular dynamics (MD) and grand canonical Monte Carlo (GCMC) simulations. The transport properties (solubility, diffusivity, and permeability) of pure and gas mixtures in the membranes were calculated and the results of the simulations were compared with the available experimental data. The simulated structural properties of the pure and composite PSF membranes including occupied volume, free volume, surface area, fractional free volume (FFV), and radius of gyration (R _g ) were evaluated and their effects on the separability of the gases by the membranes were analyzed and interpreted by the obtained results.

Nucleus-independent chemical shift criterion for aromaticity in π-extended tetraoxa[8]circulenes

Recently synthesized π-extended symmetrical tetraoxa[8]circulenes that exhibit electroluminescent properties were calculated at the density functional theory (DFT) level using the quantum theory of atoms in molecules (QTAIM) approach to electron density distribution analysis. Nucleus-independent chemical shift (NICS) indices were used to characterize the aromaticity of the studied molecules. The tetraoxa[8]circulene molecules were found to consist of two antiaromatic perimeters (according to the Hückel “4n” antiaromaticity rule) that include 8 and 24 π-electrons. Conversely, NICS calculations demonstrated the existence of a common π-extended system (distributed like a flat ribbon) in the studied tetraoxa[8]circulene molecules. Thus, these symmetrical tetraoxa[8]circulene molecules provide examples of diatropic systems characterized by the presence of induced diatropic ring currents.

Stereoselectivity of chalcone isomerase with chalcone derivatives: a computational study

Chalcone isomerase (CHI) catalyzes the intramolecular cyclization of chalcones into flavonoids. The activity of CHI is essential for the biosynthesis of flavonoids precursors of floral pigments and phenylpropanoid plant defense compounds. In the present study, we explored the detailed binding structures and binding free energies for two different active site conformations of CHI with s-cis/s-trans conformers of three chalcone compounds by performing molecular dynamics (MD) simulations and binding free energy calculations. The computational results indicate that s-cis/s-trans conformers of chalcone compounds are orientated in the similar binding position in the active site of CHI and stabilized by the different first hydrogen bond network and the same second hydrogen bond network. The first hydrogen bond network results in much lower binding affinity of s-trans conformer of chalcone compound with CHI than that of s-cis conformer. The conformational change of the active site residue T48 from indirectly interacting with the substrate via the second hydrogen bond network to directly forming the hydrogen bond with the substrates cannot affect the binding mode of both conformers of chalcone compounds, but remarkably improves the binding affinity. These results show that CHI has a strong stereoselectivity. The calculated binding free energies for three chalcone compounds with CHI are consistent with the experimental activity data. In addition, several valuable insights are suggested for future rational design and discovery of high-efficiency mutants of CHI.

Why tazobactam and sulbactam have different intermediates population with SHV-1 β-lactamase: a molecular dynamics study

The imine intermediates of tazobactam and sulbactam bound to SHV-1 β-lactamase were investigated by molecular dynamics (MD) simulation respectively. Hydrogen bond networks around active site were found different between tazobactam and sulbactam acyl-enzymes. In tazobactam imine intermediate, it was observed that the triazolyl ring formed stable hydrogen bonds with Asn170 and Thr167. The results suggest that conformation of imine determined the population of intermediates. In imine intermediate of tazobactam, the triazolyl ring is trapped in Thr_Asn pocket, and it restricts the rotation of C5-C6 bond so that tazobactam can only generate trans enamine intermediate. Further, conformational cluster analyses are performed to substantiate the results. These findings provide an explanation for the corresponding experimental results, and will be potentially useful in the development of new inhibitors.

Conformational entropy of a polymer chain grafted to rough surfaces

A polymer molecule (represented by a statistical chain) end-grafted to a topologically rough surface was studied by static MC simulations. A modified self-avoiding walk on a cubic lattice was used to model the polymer in an athermal solution. Different statistical models of surface roughness were applied. Conformational entropies of chains attached to uncorrelated Gaussian, Brownian, and fractional Brownian surfaces were calculated. Results were compared with the predictions of a simple analytical model of a macromolecule end-grafted to a fractal surface.

CO2 adsorption on polar surfaces of ZnO

Physical and chemical adsorption of CO₂ on ZnO surfaces were studied by means of two different implementations of periodic density functional theory. Adsorption energies were computed and compared to values in the literature. In particular, it was found that the calculated equilibrium structure and internuclear distances are in agreement with previous work. CO₂ adsorption was analyzed by inspection of the density of states and electron localization function. Valence bands, band gap and final states of adsorbed CO₂ were investigated and the effect of atomic displacements analyzed. The partial density of states (PDOS) of chemical adsorption of CO₂ on the ZnO(0001) surface show that the p orbitals of CO₂ were mixed with the ZnO valence band state appearing at the top of the valence band and in regions of low-energy conduction band.

Monte carlo study of the percolation in two-dimensional polymer systems

The structure of a two-dimensional film formed by adsorbed polymer chains was studied by means of Monte Carlo simulations. The polymer chains were represented by linear sequences of lattice beads and positions of these beads were restricted to vertices of a two-dimensional square lattice. Two different Monte Carlo methods were employed to determine the properties of the model system. The first was the random sequential adsorption (RSA) and the second one was based on Monte Carlo simulations with a Verdier-Stockmayer sampling algorithm. The methodology concerning the determination of the percolation thresholds for an infinite chain system was discussed. The influence of the chain length on both thresholds was presented and discussed. It was shown that the RSA method gave considerably lower thresholds for longer chains. This behavior can be explained by a different pool of chain conformations used in the calculations in both methods under consideration.

Quantitative analysis of molecular surfaces: areas, volumes, electrostatic potentials and average local ionization energies

We describe a procedure for performing quantitative analyses of fields f(r) on molecular surfaces, including statistical quantities and locating and evaluating their local extrema. Our approach avoids the need for explicit mathematical representation of the surface and can be implemented easily in existing graphical software, as it is based on the very popular representation of a surface as collection of polygons. We discuss applications involving the volumes, surface areas and molecular surface electrostatic potentials, and local ionization energies of a group of 11 molecules.

Effect of varying the 1–4 intramolecular scaling factor in atomistic simulations of long-chain N-alkanes with the OPLS-AA model

A comprehensive molecular dynamics simulation study of n-alkanes using the optimized potential for liquid simulation with all-atoms (OPLS-AA) force field at ambient condition has been performed. Our results indicate that while simulations with the OPLS-AA force field accurately predict the liquid state mass density for n-alkanes with carbon number equal or less than 10, for n-alkanes with carbon number equal or exceeding 12, the OPLS-AA force field with the standard scaling factor for the 1–4 intramolecular Van der Waals and electrostatic interaction gives rise to a quasi-crystalline structure. We found that accurate predictions of the liquid state properties are obtained by successively reducing the aforementioned scaling factor for each increase of the carbon number beyond n-dodecane. To better understand the effects of reducing the scaling factor, its influence on the torsion potential profile, and the corresponding gauche-trans conformer distribution, heat of vaporization, melting point, and self-diffusion coefficient for n-dodecane were investigated. This relatively simple procedure enables more accurate predictions of the thermo-physical properties of longer n-alkanes.

Entropy versus aromaticity in the conformational dynamics of aromatic rings

Comparison of the results of Car-Parrinello molecular dynamics simulations of isolated benzene, pyrimidine and 1,2,4-triazine molecules reveals that the unusually low population of planar geometry of the benzene ring is caused by entropy effects despite its high aromaticity. The decrease in symmetry of the molecule results in smaller changes in entropy and Gibbs free energy due to out-of-plane deformations of the ring, leading to an increase in the population of planar geometry of the ring. This leads to differences in the topology of potential energy and Gibbs free energy surfaces.

Cyclo-hexa-peptides at the water/cyclohexane interface: a molecular dynamics simulation

Molecular dynamic (MD) simulations have been performed to study the behaviors of ten kinds of cyclo-hexa-peptides (CHPs) composed of amino acids with the diverse hydrophilic/hydrophobic side chains at the water/cyclohexane interface. All the CHPs take the “horse-saddle” conformations at the interface and the hydrophilicity/hydrophobicity of the side chains influences the backbones’ structural deformations. The orientations and distributions of the CHPs at the interface and the differences of interaction energies (ΔΔE) between the CHPs and the two liquid phases have been determined. RDF analysis shows that the H-bonds were formed between the O_C atoms of the CHPs’ backbones and H_w atoms of water molecules. N atoms of the CHPs’ backbones formed the H-bonds or van der Waals interactions with the water solvent. It was found that there is a parallel relationship between ΔΔE and the lateral diffusion coefficients (D _xy ) of the CHPs at the interface. The movements of water molecules close to the interface are confined to some extent, indicating that the dynamics of the CHPs and interfacial water molecules are strongly coupled.

Assessing the quantum mechanical level of theory for prediction of linear and nonlinear optical properties of push-pull organic molecules

In this paper, we assessed the quantum mechanical level of theory for prediction of linear and nonlinear optical (NLO) properties of push-pull organic molecules. The electric dipole moment (μ), mean polarizability (〈α〉) and total static first hyperpolarizability (β_t) were calculated for a set of benzene, styrene, biphenyl and stilbene derivatives using HF, MP2 and DFT (31 different functionals) levels and over 71 distinct basis sets. In addition, we propose two new basis sets, NLO-V and aNLO-V, for NLO properties calculations. As the main outcomes it is shown that long-range corrected DFT functionals such as M062X, ωB97, cam-B3LYP, LC-BLYP and LC-ωPBE work satisfactorily for NLO properties when appropriate basis sets such as those proposed here (NLO-V or aNLO-V) are used. For most molecules with β ranging from 0 to 190 esu, the average absolute deviation was 13.2 esu for NLO-V basis sets, compared to 27.2 esu for the standard 6-31 G(2d) basis set. Therefore, we conclude that the new basis sets proposed here (NLO-V and aNLO-V), together with the cam-B3LYP functional, make an affordable calculation scheme to predict NLO properties of large organic molecules.