Second-order Møller-Plesset perturbation theory without basis set superposition error : II : open-shell systems

The basis set superposition error-free second-order MØller-Plesset perturbation theory of intermolecular interactions was studied. The difficulties of the counterpoise (CP) correction in open-shell systems were also discussed. The calculations were performed by a program which was used for testing the new variants of the theory. It was shown that the CP correction for the diabatic surfaces should be preferred to the adiabatic ones

Energy partitioning for "fuzzy" atoms

The total energy of molecule in terms of 'fuzzy atoms' presented as sum of one- and two-atomic energy components is described. The divisions of three-dimensional physical space into atomic regions exhibit continuous transition from one to another. The energy components are on chemical energy scale according to proper definitions. The Becke's integration scheme and weight function determines realization of method which permits effective numerical integrations

Counterpoise-corrected geometries and harmonic frequencies of N-body clusters: application to (HF)n (n=3,4)

A study was conducted on the methods of basis set superposition error (BSSE)-free geometry optimization and frequency calculations in clusters larger than a dimer. In particular, three different counterpoise schemes were critically examined. It was shown that the counterpoise-corrected supermolecule energy can be easily obtained in all the cases by using the many-body partitioning of energy

Are the maximum hardness and minimum polarizability principles always obeyed in nontotally symmetric vibrations?

An overview is given on a study which showed that not only in chemical reactions but also in the favorable case of nontotally symmetric vibrations where the chemical and external potentials keep approximately constant, the generalized maximum hardness principle (GMHP) and generalized minimum polarizability principle (GMPP) may not be obeyed. A method that allows an accurate determination of the nontotally symmetric molecular distortions with more marked GMPP or anti-GMPP character through diagonalization of the polarizability Hessian matrix is introduced

Relations among several nuclear and electronic density functional reactivity indexes

A set of connections among several nuclear and electronic indexes of reactivity in the framework of the conceptual Density Functional Theory by using an expansion ofthe energy functional in terms of the total number of electrons and the normal coordinates within a canonical ensemble was derived. The relations obtained provided explicit links between important quantities related to the chemical reactivity of a system. This paper particularly demonstrates that the derivative of the electronic energy with respect to the external potential of a system in its equilibrium geometry was equal to the negative of the nuclear repulsion derivative with respect to the external potential

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

Basis set and electron correlation effects on initial convergence for vibrational nonlinear optical properties of conjugated organic molecules

The level of ab initio theory which is necessary to compute reliable values for the static and dynamic (hyper)polarizabilities of three medium size π-conjugated organic nonlinear optical (NLO) molecules is investigated. With the employment of field-induced coordinates in combination with a finite field procedure, the calculations were made possible. It is stated that to obtain reasonable values for the various individual contributions to the (hyper)polarizability, it is necessary to include electron correlation. Based on the results, the convergence of the usual perturbation treatment for vibrational anharmonicity was examined

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.