Computational strategies for a four-component Dirac-Kohn-Sham program: implementation and first applications

An implementation of the generalized gradient approximation within the four-component formulation of relativistic density-functional theory using G-spinor basis sets is presented. This approach is based on the direct evaluation of the relativistic density and its gradient from the G-spinor amplitudes and gradients without explicit reference to the total density matrix. This proves to be a particularly efficient scheme, with an intrinsic computational cost that scales linearly with the number of G-spinor basis functions. In order to validate this new implementation, incorporated in the parallel version of the program BERTHA, a detailed study of the diatomic system CsAu is also reported. The spectroscopic constants D(e),r(e),omega(e), and x(e)omega(e) and the dipole moment mu have been calculated and compared with the best available theoretical and experimental data. The sensitivity of our results to the details of the numerical schemes used to evaluate the matrix elements is analyzed in detail. Also presented is a comparative study of molecular properties in the alkali auride series which have been obtained using several standard non-relativistic density functionals.     

Performance of nonrelativistic and quasi-relativistic hybrid DFT for the prediction of electric and magnetic hyperfine parameters in 57Fe Mössbauer spectra

57Fe electric and magnetic hyperfine parameters were calculated for a series of 10 iron model complexes, covering a wide range of oxidation and spin states. Employing the B3LYP hybrid method, results from nonrelativistic density functional theory (DFT) and quasi-relativistic DFT within the zero-order regular approximation (ZORA) were compared. Electron densities at the iron nuclei were calculated and correlated with experimental isomer shifts. It was shown that the fit parameters do not depend on a specific training set of iron complexes and are, therefore, more universal than might be expected. The nonrelativistic and quasi-relativistic electron densities gave fit parameters of similar quality; the ZORA densities are only shifted by a factor of 1.32, upward in the direction of the four-component Dirac-Fock value. From a correlation of calculated electric field gradients and experimental quadrupole splittings, the value of the 57Fe nuclear quadrupole moment was redetermined to a value of 0.16 barn, in good agreement with other studies. The ZORA approach gave no additional improvement of the calculated quadrupole splittings in comparison to the nonrelativistic approach. The comparison of the calculated and measured 57Fe isotropic hyperfine coupling constants (hfcc's) revealed that both the ZORA approach and the inclusion of spin-orbit contributions lead to better agreement between theory and experiment in comparison to the nonrelativistic results. For all iron complexes with small spin-orbit contributions (high-spin ferric and ferryl systems), a distinct underestimation of the isotropic hfcc's was found. Scaling factors of 1.81 (nonrelativistic DFT) and 1.69 (ZORA) are suggested. The calculated 57Fe isotropic hfcc's of the remaining model systems (low-spin ferric and high-spin ferrous systems) contain 10-50% second-order contributions and were found to be in reasonable agreement with the experimental results. This is assumed to be the consequence of error cancellation because g-tensor calculations for these systems are of poor quality with the existing DFT approaches. Excellent agreement between theory and experiment was found for the 57Fe anisotropic hfcc's. Finally, all of the obtained fit parameters were used for an application study of the [Fe(H2O)6]3+ ion. The calculated spectroscopic data are in good agreement with the Mossbauer and electron paramagnetic resonance results discussed in detail in a forthcoming paper.     

High temperature structures and orientational disorder in compressed solid nitrogen

An experimental study of the phase transitions at high temperature in compressed solid nitrogen has been performed using Raman spectroscopy. Knowledge of the equilibrium phase diagram in the region of the ordered epsilon phase and the two disordered delta and deltaloc phases, at pressures between 10 and 20 GPa, has been extended up to 500 K. The Raman scattering line shape and line width of the active vibrons has been measured accurately, along isobaric scans, across the phase transitions. Analysis of the width and of its different behavior with increasing temperature in the three phases led to more precise conclusions about the nature of the disorder in the different phases. Observation of an evident shoulder in the nu2 band of the deltaloc phase suggests the possibility that sites of two different symmetries may be occupied by the disk molecules in this structure.     

How well does molecular simulation reproduce environment-specific conformations of the intrinsically disordered peptides PLP,TP2 and ONEG?

Understanding the conformational ensembles of intrinsically disordered proteins and peptides (IDPs) in their various biological environments is essential for understanding their mechanisms and functional roles in the proteome, leading to a greater knowledge of, and potential treatments for, a broad range of diseases. To determine whether molecular simulation is able to generate accurate conformational ensembles of IDPs, we explore the structural landscape of the PLP peptide (an intrinsically disordered region of the proteolipid membrane protein) in aqueous and membrane-mimicking solvents, using replica exchange with solute scaling (REST2), and examine the ability of four force fields (ff14SB, ff14IDPSFF, CHARMM36 and CHARMM36m) to reproduce literature circular dichroism (CD) data. Results from variable temperature (VT) ¹H and Rotating frame Overhauser Effect SpectroscopY (ROESY) nuclear magnetic resonance (NMR) experiments are also presented and are consistent with the structural observations obtained from the simulations and CD. We also apply the optimum simulation protocol to TP2 and ONEG (a cell-penetrating peptide (CPP) and a negative control peptide, respectively) to gain insight into the structural differences that may account for the observed difference in their membrane-penetrating abilities. Of the tested force fields, we find that CHARMM36 and CHARMM36m are best suited to the study of IDPs, and accurately predict a disordered to helical conformational transition of the PLP peptide accompanying the change from aqueous to membrane-mimicking solvents. We also identify an α-helical structure of TP2 in the membrane-mimicking solvents and provide a discussion of the mechanistic implications of this observation with reference to the previous literature on the peptide. From these results, we recommend the use of CHARMM36m with the REST2 protocol for the study of environment-specific IDP conformations. We believe that the simulation protocol will allow the study of a broad range of IDPs that undergo conformational transitions in different biological environments.

A protocol for simulating intrinsically disordered peptides in aqueous and hydrophobic solvents is proposed. Results from four force fields are compared with experiment. CHARMM36m performs the best for the simulated IDPs in all environments.     

Thermal and viscosity characterization of two asphalt emulsions

Journal of Thermal Analysis and Calorimetry - Asphalt emulsion is a versatile product that can be used for different waterproofing applications in civil construction and highway constructions....     

Accelerating molecular modeling applications with graphics processors

Molecular mechanics simulations offer a computational approach to study the behavior of biomolecules at atomic detail, but such simulations are limited in size and timescale by the available computing resources. State-of-the-art graphics processing units (GPUs) can perform over 500 billion arithmetic operations per second, a tremendous computational resource that can now be utilized for general purpose computing as a result of recent advances in GPU hardware and software architecture. In this article, an overview of recent advances in programmable GPUs is presented, with an emphasis on their application to molecular mechanics simulations and the programming techniques required to obtain optimal performance in these cases. We demonstrate the use of GPUs for the calculation of long-range electrostatics and nonbonded forces for molecular dynamics simulations, where GPU-based calculations are typically 10-100 times faster than heavily optimized CPU-based implementations. The application of GPU acceleration to biomolecular simulation is also demonstrated through the use of GPU-accelerated Coulomb-based ion placement and calculation of time-averaged potentials from molecular dynamics trajectories. A novel approximation to Coulomb potential calculation, the multilevel summation method, is introduced and compared with direct Coulomb summation. In light of the performance obtained for this set of calculations, future applications of graphics processors to molecular dynamics simulations are discussed.     

Nanoscale conformational ordering in polyanilines investigated by SAXS and AFM

Understanding the adsorption mechanisms in nanostructured polymer films has become crucial for their use in technological applications, since film properties vary considerably with the experimental conditions utilized for film fabrication. In this paper, we employ small-angle X-ray scattering (SAXS) to investigate solutions of polyanilines and correlate the chain conformations with morphological features of the nanostructured films obtained with atomic force microscopy (AFM). It is shown that aggregates formed already in solution affect the film morphology; in particular, at early stages of adsorption film morphology appears entirely governed by the chain conformation in solution and adsorption of aggregates. We also use SAXS data for modeling poly(o-ethoxyaniline) (POEA) particle shape through an ab initio procedure based on simulated annealing using the dummy atom model (DAM), which is then compared to the morphological features of POEA films fabricated with distinct pHs and doping acids. Interestingly, when the derivative POEA is doped with p-toluene sulfonic acid (TSA), the resulting films exhibit a fibrillar morphology-seen with atomic force microscopy and transmission electron microscopy-that is consistent with the cylindrical shape inferred from the SAXS data. This is in contrast with the globular morphology observed for POEA films doped with other acids.     

Study of ionizing radiation on the properties of polyamide 6 with fiberglass reinforcement

The use of polymers reinforced with fiberglass is becoming more and more common in the switches for household industries. These compounds perform a good tension resistance to the impact and the humidity absorption being used at the present time and also are in the automobile industry in parts underneath the hood, especially in the radiator frames. The aim of this work is to study the effect of ionizing radiation on the properties of polyamide 6 with fiberglass reinforcement and undergone to different irradiation doses. Samples were prepared and irradiated on JOB 188 accelerator with an electron beam energy of 1.5 MeV in air with different doses and a dose rate of 27.99 kGy/h. Afterward, the properties of the non-irradiated and irradiated polyamide 6 with fiberglass reinforcement were evaluated.