Dynamic arrest in polymer melts: competition between packing and intramolecular barriers

We present molecular dynamics simulations of a simple model for polymer melts with intramolecular barriers. We investigate structural relaxation as a function of the barrier strength. Dynamic correlators can be consistently analyzed within the framework of the mode coupling theory of the glass transition. Control parameters are tuned in order to induce a competition between general packing effects and polymer-specific intramolecular barriers as mechanisms for dynamic arrest. This competition yields unusually large values of the so-called mode coupling theory exponent parameter and rationalizes qualitatively different observations for simple bead-spring and realistic polymers. The systematic study of the effect of intramolecular barriers presented here also establishes a fundamental difference between the nature of the glass transition in polymers and in simple glass formers.

Cluster glasses of ultrasoft particles

We present molecular dynamics (MD) simulations results for dense fluids of ultrasoft, fully penetrable particles. These are a binary mixture and a polydisperse system of particles interacting via the generalized exponential model, which is known to yield cluster crystal phases for the corresponding monodisperse systems. Because of the dispersity in the particle size, the systems investigated in this work do not crystallize and form disordered cluster phases. The clusteringtransition appears as a smooth crossover to a regime in which particles are mostly located in clusters, isolated particles being infrequent. The analysis of the internal cluster structure reveals microsegregation of the big and small particles, with a strong homo-coordination in the binary mixture. Upon further lowering the temperature below the clusteringtransition, the motion of the clusters" centers-of-mass slows down dramatically, giving way to a cluster glass transition. In the cluster glass, the diffusivities remain finite and display an activated temperature dependence, indicating that relaxation in the cluster glass occurs via particle hopping in a nearly arrested matrix of clusters. Finally we discuss the influence of the microscopic dynamics on the transport properties by comparing the MD results with Monte Carlo simulations.

Control and femtosecond time-resolved imaging of torsion in a chiral molecule

We study how the combination of long and short laser pulses can be used to induce torsion in an axially chiral biphenyl derivative (3,5-difluoro-3 ,5 -dibromo-4 -cyanobiphenyl). A long, with respect to the molecular rotational periods, elliptically polarized laser pulse produces 3D alignment of the molecules, and a linearly polarized short pulse initiates torsion about the stereogenic axis. The torsional motion is monitored in real-time by measuring the dihedral angle using femtosecond time-resolved Coulomb explosion imaging. Within the first 4 picoseconds (ps), torsion occurs with a period of 1.25 ps and an amplitude of 3◦ in excellent agreement with theoretical calculations. At larger times, the quantum states of the molecules describing the torsional motion dephase and an almost isotropic distribution of the dihedral angle is measured.We demonstrate an original application of covariance analysis of two-dimensional ion images to reveal strong correlations between specific ejected ionic fragments from Coulomb explosion. This technique strengthens our interpretation of the experimental data

Structure and dynamics of the membrane attaching nitric oxide transporter nitrophorin 7

Nitrophorins represent a unique class of heme proteins that are able to perform the delicate transportation and release of the free-radical gaseous messenger nitric oxide (NO) in a pH-triggered manner. Besides its ability to bind to phospholipid membranes, the N-terminus contains an additional Leu-Pro-Gly stretch, which is a unique sequence trait, and the heme cavity is significantly altered with respect to other nitrophorins. These distinctive features encouraged us to solve the X-ray crystallographic structures of NP7 at low and high pH and bound with different heme ligands (nitric oxide, histamine, imidazole). The overall fold of the lipocalin motif is well preserved in the different X-ray structures and resembles the fold of other nitrophorins. However, a chain-like arrangement in the crystal lattice due to a number of head-to-tail electrostatic stabilizing interactions is found in NP7. Furthermore, the X-ray structures also reveal ligand-dependent changes in the orientation of the heme, as well as in specific interactions between the A-B and G-H loops, which are considered to be relevant for the biological function of nitrophorins. Fast and ultrafast laser triggered ligand rebinding experiments demonstrate the pH-dependent ligand migration within the cavities and the exit route. Finally, the topological distribution of pockets located around the heme as well as from inner cavities present at the rear of the protein provides a distinctive feature in NP7, so that while a loop gated exit mechanism to the solvent has been proposed for most nitrophorins, a more complex mechanism that involves several interconnected gas hosting cavities is proposed for NP7.

Speed of travelling fronts: Two-dimensional and ballistic dispersal probability distributions

The speed of traveling fronts for a two-dimensional model of a delayed reactiondispersal process is derived analytically and from simulations of molecular dynamics. We show that the one-dimensional (1D) and two-dimensional (2D) versions of a given kernel do not yield always the same speed. It is also shown that the speeds of time-delayed fronts may be higher than those predicted by the corresponding non-delayed models. This result is shown for systems with peaked dispersal kernels which lead to ballistic transport

Diffusive transport properties in monovalent and divalent metal-ion halide melts: A computer simulation study

Self- and cross-velocity correlation functions and related transport coefficients of molten salts are studied by molecular-dynamics simulation. Six representative systems are considered, i.e., NaCl and KCl alkali halides, CuCl and CuBr noble-metal halides, and SrCl2 and ZnCl2 divalent metal-ion halides. Computer simulation results are compared with experimental self-diffusion coefficients and electrical conductivities. Special attention is paid to dynamic cross correlations and their dependence on the Coulomb interactions as well as on the size and mass differences between anions and cations.

Collective dynamic properties in simple crystals: Influence of the structural disorder

Collective dynamic properties in Lennard-Jones crystals are investigated by molecular dynamics simulation. The study is focused on properties such as the dynamic structure factors, the longitudinal and transverse currents and the density of states. The influence on these properties of the structural disorder is analyzed by comparing the results for one-component crystals with those for liquids and supercooled liquids at analogous conditions. The effects of species-disorder on the collective properties of binary crystals are also discussed.

Fases de la materia, átomos y deslocalización cuántica.

Este artículo pertenece a un monográfico sobre química

Fronts from complex two-dimensional dispersal kernels: theory and application to Reid's paradox

Bimodal dispersal probability distributions with characteristic distances differing by several orders of magnitude have been derived and favorably compared to observations by Nathan [Nature (London) 418, 409 (2002)]. For such bimodal kernels, we show that two-dimensional molecular dynamics computer simulations are unable to yield accurate front speeds. Analytically, the usual continuous-space random walks (CSRWs) are applied to two dimensions. We also introduce discrete-space random walks and use them to check the CSRW results (because of the inefficiency of the numerical simulations). The physical results reported are shown to predict front speeds high enough to possibly explain Reid's paradox of rapid tree migration. We also show that, for a time-ordered evolution equation, fronts are always slower in two dimensions than in one dimension and that this difference is important both for unimodal and for bimodal kernels

Additional compact formulas for vibrational dynamic dipole polarizabilities and hyperpolarizabilities

Compact expressions, complete through second order in electrical and/or mechanical anharmonicity, are given for the dynamic dipole vibrational polarizability and dynamic first and second vibrational hyperpolarizabilities. Certain contributions not previously formulated are now included

Critical thoughts on computing atom condensed Fukui functions

Different procedures to obtain atom condensed Fukui functions are described. It is shown how the resulting values may differ depending on the exact approach to atom condensed Fukui functions. The condensed Fukui function can be computed using either the fragment of molecular response approach or the response of molecular fragment approach. The two approaches are nonequivalent; only the latter approach corresponds in general with a population difference expression. The Mulliken approach does not depend on the approach taken but has some computational drawbacks. The different resulting expressions are tested for a wide set of molecules. In practice one must make seemingly arbitrary choices about how to compute condensed Fukui functions, which suggests questioning the role of these indicators in conceptual density-functional theory

Critical analysis and extension of the Hirshfeld atoms in molecules

The computational approach to the Hirshfeld [Theor. Chim. Acta 44, 129 (1977)] atom in a molecule is critically investigated, and several difficulties are highlighted. It is shown that these difficulties are mitigated by an alternative, iterative version, of the Hirshfeld partitioning procedure. The iterative scheme ensures that the Hirshfeld definition represents a mathematically proper information entropy, allows the Hirshfeld approach to be used for charged molecules, eliminates arbitrariness in the choice of the promolecule, and increases the magnitudes of the charges. The resulting "Hirshfeld-I charges" correlate well with electrostatic potential derived atomic charges

Selected configuration interaction with truncation energy error and application to the Ne atom

Selected configuration interaction (SCI) for atomic and molecular electronic structure calculations is reformulated in a general framework encompassing all CI methods. The linked cluster expansion is used as an intermediate device to approximate CI coefficients BK of disconnected configurations (those that can be expressed as products of combinations of singly and doubly excited ones) in terms of CI coefficients of lower-excited configurations where each K is a linear combination of configuration-state-functions (CSFs) over all degenerate elements of K. Disconnected configurations up to sextuply excited ones are selected by Brown's energy formula, ΔEK=(E-HKK)BK2/(1-BK2), with BK determined from coefficients of singly and doubly excited configurations. The truncation energy error from disconnected configurations, Δdis, is approximated by the sum of ΔEKS of all discarded Ks. The remaining (connected) configurations are selected by thresholds based on natural orbital concepts. Given a model CI space M, a usual upper bound ES is computed by CI in a selected space S, and EM=E S+ΔEdis+δE, where δE is a residual error which can be calculated by well-defined sensitivity analyses. An SCI calculation on Ne ground state featuring 1077 orbitals is presented. Convergence to within near spectroscopic accuracy (0.5 cm-1) is achieved in a model space M of 1.4× 109 CSFs (1.1 × 1012 determinants) containing up to quadruply excited CSFs. Accurate energy contributions of quintuples and sextuples in a model space of 6.5 × 1012 CSFs are obtained. The impact of SCI on various orbital methods is discussed. Since ΔEdis can readily be calculated for very large basis sets without the need of a CI calculation, it can be used to estimate the orbital basis incompleteness error. A method for precise and efficient evaluation of ES is taken up in a companion paper

Anharmonicity contributions to the vibrational second hyperpolarizability of conjugated oligomers

Restricted Hartree-Fock 6-31G calculations of electrical and mechanical anharmonicity contributions to the longitudinal vibrational second hyperpolarizability have been carried out for eight homologous series of conjugated oligomers - polyacetylene, polyyne, polydiacetylene, polybutatriene, polycumulene, polysilane, polymethineimine, and polypyrrole. To draw conclusions about the limiting infinite polymer behavior, chains containing up to 12 heavy atoms along the conjugated backbone were considered. In general, the vibrational hyperpolarizabilities are substantial in comparison with their static electronic counterparts for the dc-Kerr and degenerate four-wave mixing processes (as well as for static fields) but not for electric field-induced second harmonic generation or third harmonic generation. Anharmonicity terms due to nuclear relaxation are important for the dc-Kerr effect (and for the static hyperpolarizability) in the σ-conjugated polymer, polysilane, as well as the nonplanar π systems polymethineimine and polypyrrole. Restricting polypyrrole to be planar, as it is in the crystal phase, causes these anharmonic terms to become negligible. When the same restriction is applied to polymethineimine the effect is reduced but remains quantitatively significant due to the first-order contribution. We conclude that anharmonicity associated with nuclear relaxation can be ignored, for semiquantitative purposes, in planar π-conjugated polymers. The role of zero-point vibrational averaging remains to be evaluated

Linear response functions for a vibrational configuration interaction state

Linear response functions are implemented for a vibrational configuration interaction state allowing accurate analytical calculations of pure vibrational contributions to dynamical polarizabilities. Sample calculations are presented for the pure vibrational contributions to the polarizabilities of water and formaldehyde. We discuss the convergence of the results with respect to various details of the vibrational wave function description as well as the potential and property surfaces. We also analyze the frequency dependence of the linear response function and the effect of accounting phenomenologically for the finite lifetime of the excited vibrational states. Finally, we compare the analytical response approach to a sum-over-states approach