Variational calculation of vibrational linear and nonlinear optical properties

A variational approach for reliably calculating vibrational linear and nonlinear optical properties of molecules with large electrical and/or mechanical anharmonicity is introduced. This approach utilizes a self-consistent solution of the vibrational Schrödinger equation for the complete field-dependent potential-energy surface and, then, adds higher-level vibrational correlation corrections as desired. An initial application is made to static properties for three molecules of widely varying anharmonicity using the lowest-level vibrational correlation treatment (i.e., vibrational Møller-Plesset perturbation theory). Our results indicate when the conventional Bishop-Kirtman perturbation method can be expected to break down and when high-level vibrational correlation methods are likely to be required. Future improvements and extensions are discussed

Conformations of poly{G}-poly{C} π stacks with high hole mobility

Charge transfer properties of DNA depend strongly on the π stack conformation. In the present paper, we identify conformations of homogeneous poly-{G}-poly-{C} stacks that should exhibit high charge mobility. Two different computational approaches were applied. First, we calculated the electronic coupling squared, V2, between adjacent base pairs for all 1 ps snapshots extracted from 15 ns molecular dynamics trajectory of the duplex G15. The average value of the coupling squared 〈 V2 〉 is found to be 0.0065 eV2. Then we analyze the base-pair and step parameters of the configurations in which V2 is at least an order of magnitude larger than 〈 V2 〉. To obtain more consistent data, ∼65 000 configurations of the (G:C)2 stack were built using systematic screening of the step parameters shift, slide, and twist. We show that undertwisted structures (twist<20°) are of special interest, because the π stack conformations with strong electronic couplings are found for a wide range of slide and shift. Although effective hole transfer can also occur in configurations with twist=30° and 35°, large mutual displacements of neighboring base pairs are required for that. Overtwisted conformation (twist38°) seems to be of limited interest in the context of effective hole transfer. The results may be helpful in the search for DNA based elements for nanoelectronics

Estimates of electronic coupling for excess electron transfer in DNA

Electronic coupling Vda is one of the key parameters that determine the rate of charge transfer through DNA. While there have been several computational studies of Vda for hole transfer, estimates of electronic couplings for excess electron transfer (ET) in DNA remain unavailable. In the paper, an efficient strategy is established for calculating the ET matrix elements between base pairs in a π stack. Two approaches are considered. First, we employ the diabatic-state (DS) method in which donor and acceptor are represented with radical anions of the canonical base pairs adenine-thymine (AT) and guanine-cytosine (GC). In this approach, similar values of Vda are obtained with the standard 6-31 G* and extended 6-31++ G* basis sets. Second, the electronic couplings are derived from lowest unoccupied molecular orbitals (LUMOs) of neutral systems by using the generalized Mulliken-Hush or fragment charge methods. Because the radical-anion states of AT and GC are well reproduced by LUMOs of the neutral base pairs calculated without diffuse functions, the estimated values of Vda are in good agreement with the couplings obtained for radical-anion states using the DS method. However, when the calculation of a neutral stack is carried out with diffuse functions, LUMOs of the system exhibit the dipole-bound character and cannot be used for estimating electronic couplings. Our calculations suggest that the ET matrix elements Vda for models containing intrastrand thymine and cytosine bases are essentially larger than the couplings in complexes with interstrand pyrimidine bases. The matrix elements for excess electron transfer are found to be considerably smaller than the corresponding values for hole transfer and to be very responsive to structural changes in a DNA stack

Charge transfer in DNA: hole charge is confined to a single base pair due to solvation effects

We include solvation effects in tight-binding Hamiltonians for hole states in DNA. The corresponding linear-response parameters are derived from accurate estimates of solvation energy calculated for several hole charge distributions in DNA stacks. Two models are considered: (A) the correction to a diagonal Hamiltonian matrix element depends only on the charge localized on the corresponding site and (B) in addition to this term, the reaction field due to adjacent base pairs is accounted for. We show that both schemes give very similar results. The effects of the polar medium on the hole distribution in DNA are studied. We conclude that the effects of polar surroundings essentially suppress charge delocalization in DNA, and hole states in (GC)n sequences are localized on individual guanines

Electronic and vibrational contributions to first hyperpolarizability of donor-acceptor-substituted azobenzene

In this study we report on the electronic and vibrational (hyper)polarizabilities of donor–acceptorsubstituted azobenzene. It is observed that both electronic and vibrational contributions to the electric dipole first hyperpolarizability of investigated photoactive molecule substantially depend on the conformation. The contributions to the nuclear relaxation first hyperpolarizability are found to be quite important in the case of two considered isomers (cis and trans). Although the double-harmonic term is found to be the largest in terms of magnitude, it is shown that the total value of the nuclear relaxation contribution to vibrational first hyperpolarizability is a result of subtle interplay of higher-order contributions. As a part of the study, we also assess the performance of long-range-corrected density functional theory in determining vibrational contributions to electric dipole (hyper)polarizabilities. In most cases, the applied long-range-corrected exchange correlation potentials amend the drawbacks of their conventional counterparts

Mechanism of action of cyclic antimicrobial peptides

This PhD thesis is the result of the combination of experimental and computational techniques with the aim of understanding the mechanism of action of de novo cyclic decapeptides with high antimicrobial activity. By experimental techniques the influence of the replacement of the phenylalanine for tryptophan residue in their antimicrobial activity was tested and the stability in human serum was also analyzed, in order to evaluate their potential therapeutic application as antitumor agents. On the other hand, the interaction amongst the peptide BPC194 c(KKLKKFKKLQ), the best candidate from the whole library of cyclic peptides, and a model anionic membrane was simulated. The results showed a structure-function relationship derived from the stable conformation of the peptides involved in the membrane permeabilization. As a result, a rational design was performed being BPC490 the peptide with best antimicrobial activity compared with the best active peptide from the original library.

Significance of cross correlations in the stress relaxation of polymer melts

According to linear response theory, all relaxation functions in the linear regime can be obtained using time correlation functions calculated under equilibrium. In this paper, we demonstrate that the cross correlations make a significant contribution to the partial stress relaxation functions in polymer melts. We present two illustrations in the context of polymer rheology using (1) Brownian dynamics simulations of a single chain model for entangled polymers, the slip-spring model, and (2) molecular dynamics simulations of a multichain model. Using the single chain model, we analyze the contribution of the confining potential to the stress relaxation and the plateau modulus. Although the idea is illustrated with a particular model, it applies to any single chain model that uses a potential to confine the motion of the chains. This leads us to question some of the assumptions behind the tube theory, especially the meaning of the entanglement molecular weight obtained from the plateau modulus. To shed some light on this issue, we study the contribution of the nonbonded excluded-volume interactions to the stress relaxation using the multichain model. The proportionality of the bonded/nonbonded contributions to the total stress relaxation (after a density dependent "colloidal" relaxation time) provides some insight into the success of the tube theory in spite of using questionable assumptions. The proportionality indicates that the shape of the relaxation spectrum can indeed be reproduced using the tube theory and the problem is reduced to that of finding the correct prefactor. (c) 2007 American Institute of Physics

Hierarchical modeling of entangled polymers

We demonstrate that it is possible to link multi-chain molecular dynamics simulations with the tube model using a single chain slip-links model as a bridge. This hierarchical approach allows significant speed up of simulations, permitting us to span the time scales relevant for a comparison with the tube theory. Fitting the mean-square displacement of individual monomers in molecular dynamics simulations with the slip-spring model, we show that it is possible to predict the stress relaxation. Then, we analyze the stress relaxation from slip-spring simulations in the framework of the tube theory. In the absence of constraint release, we establish that the relaxation modulus can be decomposed as the sum of contributions from fast and longitudinal Rouse modes, and tube survival. Finally, we discuss some open questions regarding possible future directions that could be profitable in rendering the tube model quantitative, even for mildly entangled polymers

Deforming glassy polystyrene: Influence of pressure, thermal history, and deformation mode on yielding and hardening

The toughness of a polymer glass is determined by the interplay of yielding, strain softening, and strain hardening. Molecular-dynamics simulations of a typical polymer glass, atactic polystyrene, under the influence of active deformation have been carried out to enlighten these processes. It is observed that the dominant interaction for the yield peak is of interchain nature and for the strain hardening of intrachain nature. A connection is made with the microscopic cage-to-cage motion. It is found that the deformation does not lead to complete erasure of the thermal history but that differences persist at large length scales. Also we find that the strain-hardening modulus increases with increasing external pressure. This new observation cannot be explained by current theories such as the one based on the entanglement picture and the inclusion of this effect will lead to an improvement in constitutive modeling.

Microscopic Mechanisms of Strain Hardening in Glassy Polymers

The mechanisms underlying the increase in stress for large mechanical strains of a polymer glass, quantified by the strain-hardening modulus, are still poorly understood. In the present paper we aim to elucidate this matter and present new mechanisms. Molecular-dynamics simulations of two polymers with very different strain-hardening moduli (polycarbonate and polystyrene) have been carried out. Nonaffine displacements occur because of steric hindrances and connectivity constraints. We argue that it is not necessary to introduce the concept of entanglements to understand strain hardening, but that hardening is rather coupled with the increase in the rate of nonaffine particle displacements. This rate increases faster for polycarbonate, which has the higher strain-hardening modulus. Also more nonaffine chain stretching is present for polycarbonate. It is shown that the inner distances of such a nonaffinely deformed chain can be well described by the inner distances of the worm-like chain, but with an effective stiffness length (equal to the Kuhn length for an infinite worm-like chain) that increases during deformation. It originates from the finite extensibility of the chain. In this way the increase in nonaffine particle displacement can be understood as resulting from an increase in the effective stiffness length of the perturbed chain during deformation, so that at larger strains a higher rate of plastic events in terms of nonaffine displacement is necessary, causing in turn the observed strain hardening in polymer glasses.

Developments in near infrared spectroscopy

Vibrational spectroscopy at high excitation is an important research frontier for two reasons. Firstly, the near infrared is proving to be an important area for the analytical applications of spectroscopy, and we would therefore like to understand how the spectra we observe relate to the molecular structure of the absorbing species. Secondly, there is a fundamental interest in understanding molecular dynamics and energy flow within a polyatomic molecule at high excitation, because this is the boundary between spectroscopy and chemistry through which we try to understand the details of a chemical reaction. In this presentation I shall survey recent progress in this field.

The vibrational spectrum and ultimate modulus of polyethylene

We have performed the first completely ab initio lattice dynamics calculation of the full orthorhombic cell of polyethylene using periodic density functional theory in the local density approximation (LDA) and the generalized gradient approximation (GGA). Contrary to current perceptions, we show that LDA accurately describes the structure whereas GGA fails. We emphasize that there is no parametrization of the results. We then rigorously tested our calculation by computing the phonon dispersion curves across the entire Brillouin zone and comparing them to the vibrational spectra, in particular the inelastic neutron scattering (INS) spectra, of polyethylene (both polycrystalline and aligned) and perdeuteriopolyethylene. The F-point frequencies (where the infrared and Raman active modes occur) are in good agreement with the latest low temperature data. The near-perfect reproduction of the INS spectra, gives confidence in the results and allows Lis to deduce a number of physical properties including the elastic moduli, parallel and perpendicular to the chain. We find that the Young's modulus for an infinitely long, perfectly crystalline polyethylene is 360.2 GPa at 0 K. The highest experimental value is 324 GPa, indicating that current high modulus fibers are similar to 90% of their maximum possible strength.

Dinuclear copper and zinc complexes of a hexaazamacrocycle containing p-xylyl spacers and bridging anions: theoretical and spectroscopic studies

The hexaazamacrocycle 7,22-dimethyl-3,7,11,18,22,26-hexaazatricyclo[26.2.2.2(13,16)] tetratriaconta-1(30), 13,15,28,31,33- hexaene (Me-2[30] pbz(2)N(6)) was synthesized and characterised by single crystal X-ray diffraction. The macrocycle adopts a conformation with the two aromatic rings almost parallel at a distance of ca. 4.24 Angstrom, but displaced relative to each other by ca. 1.51 Angstrom. The protonation constants of this compound and the stability constants of its complexes with Cu2+ and Zn2+, were determined in water - methanol (9 : 1 v/v) at 25 degreesC with ionic strength 0.10 mol dm(-3) in KCl. The potentiometric and spectroscopic studies (NMR of zinc, cadmium and lead complexes, and EPR of the copper complexes) indicate the formation of only dinuclear complexes. The association constants of the dinuclear copper complex with anions ( thiocyanate, terephthalate and glyphosate) and neutral molecules (1,4-benzenedimethanol, p-xylylenediamine and terephthalic acid) were determined at 20 degreesC in methanol. The structural preferences of this ligand and of its dinuclear copper(II) complex with a variety of bridging ligands were evaluated theoretically by molecular mechanics calculations (MM) and molecular dynamics (MD) using quenching techniques.

Binding studies of a protonated dioxatetraazamacrocycle with carboxylate substrates

The tetraprotonated form of the dioxatetraazamacrocycle, 6,19-dioxa-3,9,16,22-tetraaza[22.2.2.2(11,14)]-triaconta-1(26),11,13,24, 27,29-hexaene, (H4L1)(4+), was used as the receptor for binding studies with carboxylate anionic substrates of different shapes, sizes, and charges [succinate (suc(2-)), cyclo- hexanetricarboxylate (cta(3-)), phthalate (ph(2-)), isophthalate (iph(2-)), terephthalate (tph(2-)), and benezenetricarboxylate (btc(3-))]. Association constants were determined by potentiometry in aqueous solution at 298.2 K and 0.10 M KCl and by H-1 NMR titration in D2O. The strongest association was found for the btc3- anion at 5-7 pH region. From both techniques it was possible to establish the binding preference trend of the receptor for the different substrates, and the H-1 NMR spectroscopy gave important suggestions about the type of interactions between partners and the location of the substrates in the supramolecular entities formed. The effective binding constants at pH 6 follow the order: btc(3-)>iph(2-)>cta(3-) =ph(2-)>tph(2-)>suc(2-). All the studies suggest that the anionic substrates bind to the receptor via N-H center dot center dot center dot O = C hydrogen bonds and electrostatic interactions, and the aromatic substrates can also establish pi-pi stacking interactions. The crystal structures of (H4L1)(4+) and its supramolecular assemblies with ph(2-) and tph(2-) were determined by X-ray diffraction. The last two structures showed that the association process in solid state occurs via multiple N-H center dot center dot center dot O = C hydrogen bonds with the anionic substrate located outside the macrocyclic cavity of the receptor. Molecular dynamics simulations carried out for the association of (H4L1)(4+) with tph(2-) and btC(3-) in water solution established at atomic level the existence of all interactions suggested by the experimental studies, which act cooperatively in the binding process. Furthermore, the binding free energies were estimated and the values are in agreement with the experimental ones, indicating that the binding of these two anionic substrates occurs into the receptor cavity. However, the tph(2-) has also propensity to leave the macrocyclic cavity and its molecular recognition can also happen at the top of the receptor. (C) 2008 Elsevier Ltd. All rights reserved.

Metal complexes of a dipyridine octaazamacrocycle: stability constants, structural and modelling studies

Two 28-membered octaazamacrocycles, [28]py(2)N(6) and Me-2[28]py(2)N(6), have been synthesized. The protonation constants of the N-methyl. derivative and the stability constants of its complexes with Ni2+, Cu2+, Zn2+, Cd2+, and Pb2+ were determined at 25degreesC in 0.10 mol dm(-3) KNO3. The high overall basicity of Me-2[28]py(2)N(6) is ascribed to the weaker repulsion between protonated contiguous charged ammonium sites separated by propyl chains. These studies together with NMR, UV-vis and EPR spectroscopies indicated the presence of mono- and di-nuclear species, The single crystal structure of the complex [Ni-2([28]py(2)N(6))(H2O)(4)]Cl-4.3H(2)O was determined, and showed each nickel centre in a distorted octahedral co-ordination environment. The nickel centres are held within the macrocycle at a large distance of 6.991(g) Angstrom from each other. The formation of mononuclear complexes was evaluated theoretically via molecular mechanics (MM) and molecular dynamics (MD) calculations and showed that these large macrocycles have sufficient flexibility to encapsulate metal ions with different stereo-electronic sizes. Structures for small and large metal ions are proposed.