Supramolecular Organization of Heptapyrenotide Oligomers—An in Depth Investigation by Molecular Dynamics Simulations

We present a molecular modeling study based on ab initio and classical molecular dynamics calculations, for the investigation of the tridimensional structure and supramolecular assembly formation of heptapyrenotide oligomers in water solution. Our calculations show that free oligomers self-assemble in helical structures characterized by an inner core formed by π- stacked pyrene units, and external grooves formed by the linker moieties. The coiling of the linkers has high ordering, dominated by hydrogen-bond interactions among the phosphate and amide groups. Our models support a mechanism of longitudinal supramolecular oligomerization based on interstrand pyrene intercalation. Only a minimal number of pyrene units intercalate at one end, favoring formation of very extended longitudinal chains, as also detected by AFM experiment. Our results provide a structural explanation of the mechanism of chirality amplification in 1:1 mixtures of standard heptapyrenotides and modified oligomers with covalently linked deoxycytidine, based on selective molecular recognition and binding of the nucleotide to the groove of the left-wound helix.

Phases of many flavors QCD: Lattice results

This note is based on our recent results on QCD with varying number of flavors of fundamental fermions. Topics include unusual, strong dynamics in the preconformal, confining phase, the physics of the conformal window and the role of ab-initio lattice simulations in establishing our current knowledge of the phases of many flavor QCD.

Peptide nucleic acids: Mechanism of tight binding and application in artificial nuclease

Peptide nucleic acids (PNA) are mimics of nucleic acids with a peptidic backbone. Duplexes and triplexes formed between PNA and DNA or RNA possess remarkable thermal stability, they are resistant to nuclease cleavage and can better discriminate mismatches. Understanding the mechanism for the tight binding between PNA and oligonucleotides is important for the design and development of better PNA-based drugs.^ We have performed molecular dynamics (MD) simulations of 8-mer PNA/DNA duplex and two analogous duplexes with chiral modification of PNA strand (D- or L-Alanine modification). MD simulations were performed with explicit water and Na$\sp{+}$ counter ions. The 1.5-ns simulations were carried out with AMBER using periodic boundary and particle mesh Ewald summation. The point charges for PNA monomers were derived from fitting electrostatic potentials, obtained from ab initio calculation, to atomic centers using RESP. Derived charges reveal significantly altered charge distribution on the PNA bases and predict the Watson-Crick H-bonds involving PNA to be stronger. Results from NMR studies investigating H-bond interactions between DNA-DNA and DNA-PNA base pairs in non-polar environment are consistent with this prediction. MD simulations demonstrated that the PNA strand is more flexible than the DNA strand in the same duplex. That this flexibility might be important for the duplex stability is tested by introducing modification into the PNA backbones. Results from MD simulation revealed dramatically altered structures for the modified PNA-DNA duplexes. Consistent with previous NMR results, we also found no intrachain hydrogen bonds between O7$\sp\prime$ and N1$\sp\prime$ of the neighboring residues in our MD study. Our study reveals that in addition to the lack of charge repulsion, stronger Watson-Crick hydrogen bonds together with flexible backbone are important factors for the enhanced stability of the PNA-DNA duplex.^ In a related study, we have developed an application of Gly-Gly-His-(Gly)$\sb3$-PNA conjugate as an artificial nuclease. We were able to demonstrate cleavage of single stranded DNA at a single site upon Ni(II) binding to Gly-Gly-His tripeptide and activation of nuclease with monoperoxyphthalic acid. ^

The calculation of Afρ and mass loss rate for comets

Ab initio calculations of Afρ are presented using Mie scattering theory and a Direct Simulation Monte Carlo (DSMC) dust outflow model in support of the Rosetta mission and its target 67P/Churyumov-Gerasimenko (CG). These calculations are performed for particle sizes ranging from 0.010 μm to 1.0 cm. The present status of our knowledge of various differential particle size distributions is reviewed and a variety of particle size distributions is used to explore their effect on Afρ , and the dust mass production View the MathML sourcem˙. A new simple two parameter particle size distribution that curtails the effect of particles below 1 μm is developed. The contributions of all particle sizes are summed to get a resulting overall Afρ. The resultant Afρ could not easily be predicted a priori and turned out to be considerably more constraining regarding the mass loss rate than expected. It is found that a proper calculation of Afρ combined with a good Afρ measurement can constrain the dust/gas ratio in the coma of comets as well as other methods presently available. Phase curves of Afρ versus scattering angle are calculated and produce good agreement with observational data. The major conclusions of our calculations are: – The original definition of A in Afρ is problematical and Afρ should be: qsca(n,λ)×p(g)×f×ρqsca(n,λ)×p(g)×f×ρ. Nevertheless, we keep the present nomenclature of Afρ as a measured quantity for an ensemble of coma particles.– The ratio between Afρ and the dust mass loss rate View the MathML sourcem˙ is dominated by the particle size distribution. – For most particle size distributions presently in use, small particles in the range from 0.10 to 1.0 μm contribute a large fraction to Afρ. – Simplifying the calculation of Afρ by considering only large particles and approximating qsca does not represent a realistic model. Mie scattering theory or if necessary, more complex scattering calculations must be used. – For the commonly used particle size distribution, dn/da ∼ a−3.5 to a−4, there is a natural cut off in Afρ contribution for both small and large particles. – The scattering phase function must be taken into account for each particle size; otherwise the contribution of large particles can be over-estimated by a factor of 10. – Using an imaginary index of refraction of i = 0.10 does not produce sufficient backscattering to match observational data. – A mixture of dark particles with i ⩾ 0.10 and brighter silicate particles with i ⩽ 0.04 matches the observed phase curves quite well. – Using current observational constraints, we find the dust/gas mass-production ratio of CG at 1.3 AU is confined to a range of 0.03–0.5 with a reasonably likely value around 0.1.

Microwave Spectroscopy and Structures of Several Substituted Acetylenes and Fluorocarbons

The structures of two substituted acetylene compounds have been characterized from their microwave rotational spectra. In the first study, two structures of 6-methyl-3-heptyne have been determined. This compound can be thought of as an ethyl group separated from an isobutyl group by a C≡C spacer. Both structures have the ethyl and isobutyl groups eclipsed, consistent with the dominant interaction determining the orientation about the acetylene axis being the weak dispersion attraction between the end groups. One structure is with the isobutyl group in a symmetric conformation and the other with the isobutyl group asymmetric. In addition, the microwave spectrum of the butane analogue 3,5-octadiyne has been observed. This compound consists of two ethyl groups separated by two C≡C spacers. The study is still in progress, but it appears that the ethyl end groups are freely rotating. Therefore, it seems that the dispersion attractions between the end groups are too weak at this longer distance of about 7 Å. The structures of several fluorocarbons have also been studied by microwave spectroscopy. The structures of perfluoropentane and perfluorohexane have been shown to be helical, like the polymer polytetrafluoroethylene (Teflon©). The structure of perfluoropropane and two conformers of 1H-heptafluoropropane have been determined to be non-helical. It is apparent that the steric and dipole repulsions between fluorine atoms that have been attributed to the helical structure of longer fluorocarbon chains are not sufficient in a three carbon chain to cause a twist in the structures.

Estudio a escala laboratorio y planta piloto, de la adsorción de NO y SO2, sobre CR2O3/AL2O3 a altas temperaturas, provenientes de fuentes fijas

Resumen: Se planificaron las experiencias con el objeto de analizar el comportamiento del catalizador en la columna metálica de mayor diámetro. Se modificaron las masas usadas para verificar la eficiencia de retención respecto de la masa. Se realizaron ciclos de adsorción, desorción y readsorción sobre una misma muestra para determinar variaciones en la eficiencia del catalizador. En otra fase, en colaboración con el Dr. V. A. Ranea y el Prof. E. E. Mola (INIFTA, UNLP), se desarrolló el estudio teórico de la adsorción de moléculas de SO2, CH4, CO2, O2 y CO sobre Cr2O3(0001) mediante Teoría del Funcional Densidad (programa VASP, Vienna Ab-initio Simulation Package), y el estudio de la cinética de la reacción entre CH4, SO2 y el O2 junto con la presencia de especies sulfito y sulfato. Este estudio permitió hallar los sitios preferenciales de adsorción de Sº y la posible competencia con SO2 experimentalmente y por cálculos teóricos. Dentro del marco de la presente línea de investigación, la Ing. Sabrina Hernández Guiance continúa realizando experiencias en el marco del proyecto conjunto con el INIFTA, las cuales forman parte del desarrollo de su tesis doctoral. Experimentalmente, se observa que la eficiencia de adsorción del catalizador respecto al SO2 es cercana al 100%. Se observa un pico de termodesorción a 1120 K. Luego, se estudió la oxidación de CH4 con SO2. Se observa que hay producción de CO2 desde temperatura inicial, seguida de un aumento significativo en la formación de CO2 hasta 330-340 K. Luego, la producción de CO2 se mantiene aproximadamente constante. Mediante el empleo de la ecuación de Arrhenius y resultados experimentales, se obtuvo la energía de activación de la reacción global, de 7 Kcal/mol. Mediante estudios teóricos, se determinó que la energía de quimisorción del SO2 sobre el Cr2O3 es de -3.09 eV para la configuración más estable, una energía de adsorción de O2 en estado disociativo de -1.567 eV, una energía para CH4 sobre O2 adsorbido previamente de -0.335 eV, y - 0.812 eV para la configuración más estable de CO2 sobre el sustrato.

Quantum chemical investigation on iodine oxides and their role in the formation of atmospheric aerosols

En esta contribución se presenta un estudio teórico de diferentes reacciones químicas entre óxidos de iodo y agua que pueden contribuir a la formación de partículas en la atmósfera. Mediante el uso de cálculos quimicocuánticos ab initio con tratamiento de la correlación electrónica se han obtenido propiedades termodinámicas para caracterizar esas reacciones

Simulación atomística de metales líquidos y aleaciones para reactores de fusión nuclear = Atomistic simulation of liquid metals and alloys for nuclear fusion reactors

La presente tesis comprende un estudio de metales líquidos, Li, Pb y eutéctico Li17Pb en el ámbito de la tecnología de fusión nuclear. Uno de los problemas fundamentales en futuros reactores de fusión es la producción y extracción de tritio en la denominada envoltura regeneradora (blanket en inglés). Dicho blanket tendrá dos propósitos, la extracción del calor generado por las reacciones de fusión para su posterior conversión en energía eléctrica así como la producción de tritio para realimentar el proceso. Dicha producción se realizará mediante el “splitting” del Li con los neutrones provenientes de la fusión. Esta reacción produce tritio y helio por lo que la interacción del T y el He con el metal líquido, con los materiales estructurales así como con el He es un problema fundamental aun no bien entendido y de gran importancia para futuros diseños. Los capítulos 1 2 y 3 presentan una introducción a dichos problemas. El capítulo 1 introduce al lector en la tecnología de fusión nuclear. El segundo capítulo explica en mayor detalle el uso de metales líquidos en reactores de fusión, no solo en blankets sino también como primera pared, divertor etc, lo que se denomina en general “plasma facing materials”. Por último se ofrece una breve introducción a las técnicas de dinámica molecular clásica (CMD) y un breve resumen de los potenciales más usados. El estudio se ha llevado a cabo utilizando simulación atomística mediante potenciales semi-empíricos del tipo átomo embebido (EAM). La Tesis consta de 4 partes bien definidas. En primer lugar se verificó la idoneidad de los potenciales tipo EAM para simular las propiedades de los metales Li y Pb en fase líquida. Dicho estudio se detalla en el Capítulo 4 y en su extensión, el Apéndice 1, en el que se estudia los límites de validez de esta aproximación. Los resultados de dicho estudio han sido publicados y presentados en diversos congresos internacionales. Un resumen de la metodología seguida fue publicado como capítulo de libro en Technofusión 2011. Los resultados se presentaron en diversos congresos internacionales, entre ellos ICENES 2011, (Artículo en ICENES Proceedings) ICOPS-SOFE 2011, en una presentación oral etc. El trabajo ha sido aceptado recientemente en Journal of Nuclear Materiales (Fraile et al 2012). La segunda parte y más importante comprende el desarrollo de un potencial para el estudio de la mezcla de ambos metales. Éste es el trabajo más novedoso e importante dado que no existía en la literatura un potencial semejante. Se estudiaron dos aproximaciones distintas al problema, un potencial tipo EAM/cd y un potencial EAM/alloy. Ambos potenciales dan resultados satisfactorios para la simulación del eutéctico (y concentraciones de Li menores que el 17%). Sin embargo el sistema LiPb en todas las concentraciones es un sistema que se aparta enormemente de una solución ideal y dicho potencial no daba buenos resultados para mezclas PbLi con concentraciones de Li grandes. Este problema fue solventado mediante el desarrollo de un segundo potencial, esta vez tipo EAM/alloy (segunda parte del Capítulo 5). Dicho trabajo será enviado a Physical Review Letters o a Phys. Rev. B, y una extensión junto con un estudio detallado de las propiedades del eutéctico de acuerdo con nuestras simulaciones se enviará a continuación a Phys. Rev. B. En tercer lugar se estudió el problema de la difusividad del H en metales líquidos aprovechando distintos potenciales existentes en la literatura. El problema del H en metales líquidos es importante en metalurgia. En dicho capítulo se estudió la difusividad del H en Pd, Ni y Al con potenciales tipo EAM, y también con un potencial más sofisticado que tiene en cuenta la dependencia angular de las interacciones (ADP por sus siglas en inglés). De este modo disponemos de un estudio detallado del problema con diferentes modelos y diferentes metales. La conclusión apunta a que si se compara con los resultados experimentales (muy escasos) los resultados obtenidos mediante CMD dan valores bajos de la difusividad del H. Las razones de dicho desacuerdo entre simulación y experimentos se detallan en el Capítulo 6. Este trabajo ha sido presentado en una presentación oral en el reciente congreso internacional “Trends on Nanothecnology” TNT 2012 celebrado en Madrid. El trabajo será publicado en un futuro próximo. Por último, como se dijo anteriormente, el estudio del He, la formación de burbujas en metales líquidos, su difusión nucleación y cavitación es otro problema deseable de ser estudiado mediante técnicas atomísticas. Para ello es necesario el desarrollo de diversos potenciales, He-Li, He-Pb y un potencial ternario Pb-Li-He. Para ello se han realizado simulación ab initio de los sistemas Pb+He y Li+He. Dicho estudio pretende calcular las fuerzas entre los átomos del metal (Pb o Li) con intersticiales de He. De este modo aplicaremos el “force matching method” (FMM) para el desarrollo de dichos potenciales. En el Capítulo 7 se detallan los resultados obtenidos referidos a las posiciones más favorables de las impurezas de He dentro de redes cristalinas de Pb y Li así como el efecto de tener en cuenta el acoplo spin-orbita (SOC en inglés). El análisis de los resultados en términos de transferencia de carga y análisis de las densidades electrónicas, así como la creación de los potenciales mencionados está en progreso. En conjunto la tesis presenta un estudio de los diversos problemas relacionados con el uso de metales líquidos en reactores de fusión y representa un primer paso en la determinación de parámetros de gran importancia para el diseño de blankets y sistemas de primera pared. Con la simulación MD de dichos problemas mediante, importante, potenciales realistas, valores de difusión, solubilidad etc de especies ligeras, H (o sus isotopos) y He en metales líquidos podrá ser calculada complementando así la base de datos que presenta enormes incertidumbres.

Raman and infrared spectra of dimethyl ether 13C-isotopologue (CH3O13CH3) from a CCSD(T) potential energy surface

So far, no experimental data of the infrared and Raman spectra of 13C isotopologue of dimethyl ether are available. With the aim of providing some clues of its low-lying vibrational bands and with the hope of contributing in a next spectral analysis, a number of vibrational transition frequencies below 300 cm−1 of the infrared spectrum and around 400 cm−1 of the Raman spectrum have been predicted and their assignments were proposed. Calculations were carried out through an ab initio three dimensional potential energy surface based on a previously reported one for the most abundant dimethyl ether isotopologue (M. Villa et al., J. Phys. Chem. A 115 (2011) 13573). The potential function was vibrationally corrected and computed with a highly correlated CCSD(T) method involving the COC bending angle and the two large amplitude CH3 internal rotation degrees of freedom. Also, the Hamiltonian parameters could represent a support for the spectral characterization of this species. Although the computed vibrational term values are expected to be very accurate, an empirical adjustment of the Hamiltonian has been performed with the purpose of anticipating some workable corrections to any possible divergence of the vibrational frequencies. Also, the symmetry breaking derived from the isotopic substitution of 13C in the dimethyl ether was taken into account when the symmetrization procedure was applied.

CCSD(T) Study of CD3-O-CD3 and CH3-O-CD3 Far-Infrared Spectra

From a vibrationally corrected 3D potential energy surface determined with highly correlated ab initio calculations (CCSD(T)), the lowest vibrational energies of two dimethyl-ether isotopologues, 12CH3–16O–12CD3 (DME-d3) and 12CD3–16O–12CD3 (DME-d6), are computed variationally. The levels that can be populated at very low temperatures correspond to the COC-bending and the two methyl torsional modes. Molecular symmetry groups are used for the classification of levels and torsional splittings. DME-d6 belongs to the G36 group, as the most abundant isotopologue 12CH3–16O–12CH3 (DME-h6), while DME-d3 is a G18 species. Previous assignments of experimental Raman and far-infrared spectra are discussed from an effective Hamiltonian obtained after refining the ab initio parameters. Because a good agreement between calculated and experimental transition frequencies is reached, new assignments are proposed for various combination bands corresponding to the two deuterated isotopologues and for the 020 → 030 transition of DME-d6. Vibrationally corrected potential energy barriers, structural parameters, and anharmonic spectroscopic parameters are provided. For the 3N – 9 neglected vibrational modes, harmonic and anharmonic fundamental frequencies are obtained using second-order perturbation theory by means of CCSD and MP2 force fields. Fermi resonances between the COC-bending and the torsional modes modify DME-d3 intensities and the band positions of the torsional overtones.

Band gap control via tuning of inversion degree in CdIn2S4 spinel

Based on theoretical arguments, we propose a possible route for controlling the band-gap in the promising photovoltaic material CdIn2S4. Our ab initio calculations show that the experimental degree of inversion in this spinel (fraction of tetrahedral sites occupied by In) corresponds approximately to the equilibrium value given by the minimum of the theoretical inversion free energy at a typical synthesis temperature. Modification of this temperature, or of the cooling rate after synthesis, is then expected to change the inversion degree, which in turn sensitively tunes the electronic band-gap of the solid, as shown here by screened hybrid functional calculations.

Estructura electrónica de fotocatalizadores compuestos ZnO-TiO2

Recientes publicaciones han mostrado propiedades fotocatalíticas interesantes en sistemas basados en TiO2 y ZnO. En unos casos hay presentes fases de ambos óxidos binarios en íntimo contacto [1] y en otros se tienen óxidos mixtos (titanatos de Zn) de distintas estequiometrías [2]; estos últimos, además, se han podido dopar con nitrógeno para obtener actividad con luz visible [3]. Las características electrónicas relevantes de estos sistemas (posición relativa de los niveles de ambas fases en el primero, estructura de bandas para los titanatos con o sin N) se conocen muy poco. Aquí se realiza un estudio teórico cuántico de estos materiales, usando para mayor exactitud funcionales híbridos (pues es sabido que la DFT estándar predice mal los bandgaps). Además se tienen en cuenta desarrollos teóricos recientes que permiten determinar ab initio, para semiconductores de gap alto, el coeficiente más adecuado de mezcla de intercambio HF [4, 5], y formular reglas para obtener con más exactitud el alineamiento de bandas que se establece a través de una interfaz [5, 6].

A theoretical study on the formation of iodine oxide aggregates and monohydrates

Biotic and abiotic emissions of molecular iodine and iodocarbons from the sea or the ice surface and the intertidal zone to the coastal/polar marine boundary layer lead to the formation of iodine oxides, which subsequently nucleate forming iodine oxide particles (IOPs). Although the link between coastal iodine emissions and ultrafine aerosol bursts is well established, the details of the nucleation mechanism have not yet been elucidated. In this paper, results of a theoretical study of a range of potentially relevant aggregation reactions of different iodine oxides, as well as complexation with water molecules, are reported. Thermochemical properties of these reactions are obtained from high level ab initio correlated calculations including spin–orbit corrections. The results show that the nucleation path most likely proceeds through dimerisation of I2O4. It is also shown that water can hinder gas-to-particle conversion to some extent, although complexation with key iodine oxides does not remove enough of these to stop IOP formation. A consistent picture of this process emerges from the theoretical study presented here and the findings of a new laboratory study reported in the accompanying paper (Gomez Martin et al., 2013).

Theoretical investigation on iodine oxides formation in the production of atmospheric aerosols

In this contribution, results of a theoretical study on different reactions that odine oxides, in the presence of water, can undergo to form iodine oxides particles in the atmosphere. Thermodynamic and kinetic properties of these reactions have been obtained at high level ab initio correlated calculations.

Double ion implantation and pulsed laser melting processes for third generation solar cells

In the framework of the third generation of photovoltaic devices, the intermediate band solar cell is one of the possible candidates to reach higher efficiencies with a lower processing cost. In this work, we introduce a novel processing method based on a double ion implantation and, subsequently, a pulsed laser melting (PLM) process to obtain thicker layers of Ti supersaturated Si. We perform ab initio theoretical calculations of Si impurified with Ti showing that Ti in Si is a good candidate to theoretically form an intermediate band material in the Ti supersaturated Si. From time-of-flight secondary ion mass spectroscopy measurements, we confirm that we have obtained a Ti implanted and PLM thicker layer of 135 nm. Transmission electron microscopy reveals a single crystalline structure whilst the electrical characterization confirms the transport properties of an intermediate band material/Si substrate junction. High subbandgap absorption has been measured, obtaining an approximate value of 104 cm−1 in the photons energy range from 1.1 to 0.6 eV.