Influence of a Counterion on the Ion Atmosphere of an Anion: A Molecular Dynamics Study of LiX and CsX (X = F-, Cl-, I-) in Methanol

We report molecular dynamics (MD) simulations to explore the influence of a counterion on the structure and dynamics of cationic and anionic solvation shells for various ions in methanol at 298 K. We show that the variation in ionic size of either the cation or the anion in an ion pair influences the solvation structure of the other ion as well as the diffusivity in an electrolyte solution of methanol. The extent of ionic association between the cation and its counteranion of different ionic sizes has been investigated by analyzing the radial distribution functions (RDFs) and the orientation of methanol molecules in the first solvation shell (FSS) of ions. It is shown that the methanol in the FSS of the anion as well the cation exhibit quite different radial and orientational structures as compared to methanol which lie in the FSS of either the anion or the cation but not both. We find that the coordination number (CN) of F-, Cr-, and I- ions decreases with increasing size of the anion which is contrary to the trend reported for the anions in H2O. The mean residence time (MRT) of methanol molecules in the FSS of ions has been calculated using the stable states picture (SSP) approach. It is seen that the ion-counterion interaction has a considerable influence on the MRT of methanol molecules in the FSS of ions. We also discuss the stability order of the ion-counterion using the potentials of mean force (PMFs) for ion pairs with ions of different sizes. The PMF plots reveal that the Li+-F- pair (small-small) is highly stable and the Li+-I- pair is least stable (small-large) in electrolyte solutions.

Electroluminescence in Molecular Junctions: A Diagrammatic Approach

We compute electroluminescent signal in a current carrying single molecule junction using a superoperator formalism. Liouville space loop diagrams are used to identify all density matrix pathways that emit photons via the electroluminescence process. A frequency resolved spectrum is expressed in terms of the various Pock space states of the isolated molecule that participate in the creation and subsequent recombination of exciton. Application is made to a multilevel Coulomb blockade model system and to a gold-benzene-1,4-dithiol-gold molecular junction.

Ab Initio Molecular Dynamics Simulations of Amino Acids in Aqueous Solutions: Estimating pK(a) Values from Metadynamics Sampling

Changes in the protonation and deprotonation of amino acid residues in proteins play a key role in many biological processes and pathways. Here, we report calculations of the free-energy profile for the protonation deprotonation reaction of the 20 canonical alpha amino acids in aqueous solutions using ab initio Car-Parrinello molecular dynamics simulations coupled with metad-ynamics sampling. We show here that the calculated change in free energy of the dissociation reaction provides estimates of the multiple pK(a) values of the amino acids that are in good agreement with experiment. We use the bond-length-dependent number of the protons coordinated to the hydroxyl oxygen of the carboxylic and the amine groups as the collective variables to explore the free-energy profiles of the Bronsted acid-base chemistry of amino acids in aqueous solutions. We ensure that the amino acid undergoing dissociation is solvated by at least three hydrations shells with all water molecules included in the simulations. The method works equally well for amino acids with neutral, acidic and basic side chains and provides estimates of the multiple pK(a) values with a mean relative error, with respect to experimental results, of 0.2 pK(a) units.

A Pd-8 Tetrafacial Molecular Barrel as Carrier for Water Insoluble Fluorophore

A new carbazole-based tetraimidazole ligand 1,3,6,8-tetra(1H-imidazol-1-yl)-9-methyl-9H-carbazole (L) has been synthesized. The unsymmetrical nature of L as well as the rotational freedom of imidazole donor moieties around C-N bond make it a special building unit, which upon treatment with cis-(tmeda)Pd(NO3)(2) produced an unprecedented single linkage-isomeric Pd-8 tetrafacial molecular nanobarrel (PSMBR-1) tmeda N,N,N',N'-tetramethylethane-1,2-diamine]. Unlike closed architectures, open barrel architecture of water-soluble PSMBR-1 makes it an ideal host for some water insoluble polyaromatic hydrocarbons in aqueous medium; one such inclusion complex coroneneCPSMBR-1 was characterized by X-ray diffraction study. Moreover, the potential application of PSMER-1 as carrier in aqueous medium for the transportation of water insoluble fluorophore (perylene) for live cell imaging is explored.

Molecular Dynamics Study of the Structure, Flexibility, and Hydrophilicity of PETIM Dendrimers: A Comparison with PAMAM Dendrimers

A new class of dendrimers, the poly(propyl ether imine) (PETIM) dendrimer, has been shown to be a novel hyperbranched polymer having potential applications as a drug delivery vehicle. Structure and dynamics of the amine terminated PETIM dendrimer and their changes with respect to the dendrimer generation are poorly understood. Since most drugs are hydrophobic in nature, the extent of hydrophobicity of the dendrimer core is related to its drug encapsulation and retention efficacy. In this study, we carry out fully atomistic molecular dynamics (MD) simulations to characterize the structure of PETIM (G2-G6) dendrimers in salt solution as a function of dendrimer generation at different protonation levels. Structural properties such as radius of gyration (R-g), radial density distribution, aspect ratio, and asphericity are calculated. In order to assess the hydrophilicity of the dendrimer, we compute the number of bound water molecules in the interior of dendrirner as well as the number of dendrimer-water hydrogen bonds. We conclude that PETIM dendrimers have relatively greater hydrophobicity and flexibility when compared with their extensively investigated PAMAM counterparts. Hence PETIM dendrimers are expected to have stronger interactions with lipid membranes as well as improved drug encapsulation and retention properties when compared with PAMAM dendrimers. We compute the root-mean-square fluctuation of dendrimers as well as their entropy to quantify the flexibility of the dendrimer. Finally we note that structural and solvation properties computed using force field parameters derived based on the CHARMM general purpose force field were in good quantitative agreement with those obtained using the generalized Amber force field (GAFF).

Molecular structure of the discotic liquid crystalline phase of hexa-peri-hexabenzocoronene/oligothiophene hybrid and their charge transport properties

Using atomistic molecular dynamics simulation, we study the discotic columnar liquid crystalline (LC) phases formed by a new organic compound having hexa-peri-Hexabenzocoronene (HBC) core with six pendant oligothiophene units recently synthesized by Nan Hu et al. Adv. Mater. 26, 2066 (2014)]. This HBC core based LC phase was shown to have electric field responsive behavior and has important applications in organic electronics. Our simulation results confirm the hexagonal arrangement of columnar LC phase with a lattice spacing consistent with that obtained from small angle X-ray diffraction data. We have also calculated various positional and orientational correlation functions to characterize the ordering of the molecules in the columnar arrangement. The molecules in a column are arranged with an average twist of 25 degrees having an average inter-molecular separation of similar to 5 angstrom. Interestingly, we find an overall tilt angle of 43 degrees between the columnar axis and HBC core. We also simulate the charge transport through this columnar phase and report the numerical value of charge carrier mobility for this liquid crystal phase. The charge carrier mobility is strongly influenced by the twist angle and average spacing of the molecules in the column. (C) 2015 AIP Publishing LLC.

Towards an exact factorization of the molecular wave function

An exact single-product factorisation of the molecular wave function for the timedependent Schrodinger equation is investigated by using an ansatz involving a phasefactor. By using the Frenkel variational method, we obtain the Schrodinger equations for the electronic and nuclear wave functions. The concept of a potential energy surface (PES) is retained by introducing a modified Hamiltonian as suggested earlier by Cederbaum. The parameter in the phase factor is chosen such that the equations of motion retain the physically appealing Born- Oppenheimer-like form, and is therefore unique.

Crystal chemistry and photomechanical behavior of 3,4-dimethoxycinnamic acid: correlation between maximum yield in the solid-state topochemical reaction and cooperative molecular motion

A new monoclinic polymorph, form II (P2(1)/c, Z = 4), has been isolated for 3,4-dimethoxycinnamic acid (DMCA). Its solid-state 2 + 2 photoreaction to the corresponding alpha-truxillic acid is different from that of the first polymorph, the triclinic form I (P (1) over bar, Z = 4) that was reported in 1984. The crystal structures of the two forms are rather different. The two polymorphs also exhibit different photomechanical properties. Form I exhibits photosalient behavior but this effect is absent in form II. These properties can be explained on the basis of the crystal packing in the two forms. The nanoindentation technique is used to shed further insights into these structure-property relationships. A faster photoreaction in form I and a higher yield in form II are rationalized on the basis of the mechanical properties of the individual crystal forms. It is suggested that both Schmidt-type and Kaupp-type topochemistry are applicable for the solid-state trans-cinnamic acid photodimerization reaction. Form I of DMCA is more plastic and seems to react under Kaupp-type conditions with maximum molecular movements. Form II is more brittle, and its interlocked structure seems to favor Schmidt-type topochemistry with minimum molecular movement.

Rational molecular dynamics scheme for predicting optimum concentration loading of nano-additive in phase change materials

The present study deals with the diffusion and phase transition behaviour of paraffin reinforced with carbon nano-additives namely graphene oxide (GO) and surface functionalized single walled carbon nanotubes (SWCNT). Bulk disordered systems of paraffin hydrocarbons impregnated with carbon nano-additives have been generated in realistic equilibrium conformations for potential application as latent heat storage systems. Ab initio molecular dynamics(MD) in conjugation with COMPASS forcefield has been implemented using periodic boundary conditions. The proposed scheme allows determination of optimum nano-additive loading for improving thermo-physical properties through analysis of mass, thermal and transport properties; and assists in determination of composite behaviour and related performance from microscopic point of view. It was observed that nanocomposites containing 7.8% surface functionalised SWCNT and 55% GO loading corresponds to best latent heat storage system. The propounded methodology could serve as a by-pass route for economically taxing and iterative experimental procedures required to attain the optimum composition for best performance. The results also hint at the large unexplored potential of ab-initio classical MD techniques for predicting performance of new nanocomposites for potential phase change material applications. (C) 2015 Author(s).