Prediction of ionic conductivity from adiabatic heating in non-equilibrium molecular dynamics on various test systems.

CONTEXT: The evaluation of ionic conductivity through atomistic modeling typically involves calculating diffusion coefficients, which often necessitates simulations spanning several hundreds of nanoseconds. This study introduces a less computationally demanding approach based on non-equilibrium molecular dynamics applicable to a wide range of systems. METHOD: Ionic conductivity is determined by evaluating the Joule heating effect recorded during non-equilibrium molecular dynamics (NEMD) simulations. These simulations which involve applying a uniform electric field using classical force fields in LAMMPS are conducted within the MedeA software environment. The conductivity value for a specific temperature can thus be obtained from a single simulation together with an estimation of the associated uncertainty. Guidelines for selecting NEMD parameters such as electric field intensity and initial temperature are proposed to satisfy linear irreversible transport. RESULTS: The protocol presented in this study is applied to four different types of systems, namely, (i) molten NaCl, (ii) NaCl and LiCl aqueous solutions, (iii) solution of ionic liquid with two solvents, and (iv) NaX zeolites in the anhydrous and hydrated states. The main advantages of the proposed protocol are simplicity of implementation (eliminating the need to store individual ion trajectories), reliability (low electric field, linear response, no perturbation of the equations of motion by a thermostat), and a wide range of applications. The estimated contribution of field-induced drift motion of ions to kinetic energy appears very low, justifying the use of standard kinetic energy in the method. For each system, the reported influence of temperature, ion concentration, solvent nature, or hydration is correctly predicted.

Amorphous cis-1,4-polybutadiene P-V-T properties from atomistic simulations.

CONTEXT: As a result of the diversity of microstructures encountered in cis-1,4-polybutadiene and the variety of measurement methods used, experimental values of variation of glass transition temperature (Tg) with pressure are relatively dispersed. However, atomistic simulations enable access to valuable information for very well-controlled chemistry and structures with a well-defined and systematic acquisition protocol. By varying the temperature and pressure, the specific volume of the melt was computed, yielding results that deviated by only 2% from experimental data. A linear relationship between Tg and pressure was observed, with Tg predicted to be 162 K at zero pressure and a rate of change of Tg with respect to pressure (dTg/dP) of 0.24 K/MPa. METHOD: The atomistic dilatometry experiments were conducted on a model of amorphous cis-1,4 polybutadiene with an approximate molecular weight of 5400 g/mol using the LAMMPS code and the all-atom forcefield pcff + . The dilatometry process involved cooling and heating at a rate of 9 × 1012 K/min. The specific volume was calculated by averaging over seven independent configurations for each temperature. The Tait equation was employed to fit the specific volume evolution within the temperature range of 10 to 700 K under different pressures of 0, 60, and 100 MPa.

Computing thermomechanical properties of dry homopolymers used as raw materials for formulation of biomedical hydrogels.

Different static properties have been calculated with COMPASS force field for polyacrylamide, poly(2-hydroxyethylacrylate) (HEA), poly(2-hydroxyethylmethacrylate) (HEMA), poly(glycidylmethacrylate) (GMA), polyethylene glycol (PEG), and poly(2,2,2-trifluoroethylmethacrylate) (TFEM). For each polymers, the calculated values were averaged on five equilibrated configurations of amorphous cell composed of one atactic chain containing 100 repeat units. The ranking obtained from the densities calculated at 300 K is TFEM > HEA ≈ xpolycrylamide > HEMA ≈ GMA > PEG. Concerning the glass transition temperature we have obtained polyacrylamide > HEMA ≈ GMA ≈ HEA > PEG, and polyacrylamide > HEMA ≈ HEA > GMA ≈ PEG > TFEM for the bulk modulus. The calculated results, when available, have been compared with experimental data coming from literature.

Density functional theory calculations on azobenzene derivatives: a comparative study of functional group effect.

Density functional theory (DFT) calculations have been used to investigate the structural properties, dipole moments, polarizabilities, Gibbs energies, hardness, electronegativity, HOMO/LUMO energies, and chemical potentials of trans and cis configurations of eight para-substituted azobenzene derivatives. All properties have been obtained using the B3LYP functional and 6-31++G(d,p) basis set. The planar structures have been obtained for all optimized trans configurations. The energy difference between trans and cis configurations for considered derivatives was found to be between 64.2-73.1 kJ/mole. It has been obtained that the p-aminodiazo-benzene (ADAB) has the difference in the dipole moments between trans and cis forms higher than for trans and cis azobenzene.

Decomposition mechanisms of trinitroalkyl compounds: a theoretical study from aliphatic to aromatic nitro compounds.

The chemical mechanisms involved in the decomposition of trinitroethyl compounds were studied for both aliphatic and aromatic derivatives using density functional theory calculations. At first, in the case of 1,1,1-trinitrobutane, used as a reference molecule, two primary channels were highlighted among the five investigated ones: the breaking of the C-N bond and the HONO elimination. Then, the influence of various structural parameters was studied for these two reactions by changing the length of the carbon chain, adding substituents or double bonds along the carbon chain. If some slight changes in activation energies were observed for most of these features, no modification of the competition between the two investigated reactions was highlighted and the breaking of the C-N bond remained the favoured mechanism. At last, the reactions involving the trinitroalkyl fragments were highlighted to be more competitive than reactions involving nitro groups linked to aromatic cycles in two aromatic systems (4-(1,1,1-trinitrobutyl)-nitrobenzene and 2-(1,1,1-trinitrobutyl)-nitrobenzene). This showed that aromatic nitro compounds with trinitroalkyl derivatives decompose from their alkyl part and may be considered more likely as aliphatic than as aromatic regarding the initiation of their decomposition process.

A DFT study on the initial stage of thermal degradation of Poly(methyl methacrylate)/carbon nanotube system.

DFT calculations, with VWN exchange correlation functional and double numeric basis set, were used to evaluate the energies required for the scission reactions taking place in the initial stage of the thermal degradation of Poly(methyl methacrylate) (PMMA) in the presence of a carbon nanotube (CNT). Side group and main chain scissions were investigated. The results averaged from five configurations of pure PMMA (DP=5) were used as references and compared to the results obtained for the five same configurations of PMMA grafted on three carbon nanotubes of similar diameter (1.49 nm). The bond dissociation energies (BDE) of main chain scission evaluated for grafted PMMA was 4 % less endothermic than for pure PMMA. These results seemed independent of the tested chirality (11,11); (12,10) and (16,5) of the carbon nanotubes. Comparisons with the BDE of the weakest bonds due to intrinsic defaults (head to head and unsaturated end chain) were performed.

Calculations of the charge distribution in dodecyltrimethylammonium: a quantum chemical investigation.

In this work we report the atomic partial charges evaluated on dodecyltrimethylammonium ion. The values obtained from 17 quantum methods [CHELP, CHELPG, MK, NPA at (HF, LDA, PBE, B3LYP)//6-31G++(d,p) level and APT at B3LYP//6-31G++(d,p)] on the molecule optimised at B3LYP/6-31G++(d,p) level were compared to two semiempirical methods (Gasteiger and Qeq) and the commercial force field PCFF. All methods based on quantum calculation gave a positive charge delocalised on at least the first four alkyl groups of the tail. However, those deriving partial charges from the electrostatic potential gave an unrealistic set of alternative positive and negative alkyl group charges along the tail. In comparison, the NPA and APT methods lead to a steady decrease in the partial charges from the third alkyl group, and agreed closely with the representation of the electrostatic potential mapped onto the 0.002 au isodensity surface. The choice of the exchange correlation treatment does not drastically influence the atomic partial charges.

Density functional theory characterisation of 4-hydroxyazobenzene.

We report the structural properties, infrared (IR) and Raman spectra, dipole moment, polarisability, hardness and chemical potential of the trans and cis configurations of 4-hydroxyazobenzene calculated using the B3LYP functionals. All calculations were performed with the following basis sets: 6-31G, 6-31++G, 6-31G(d,p), 6-31++G(d,p), 6-31G(2d,2p), 6-31++G(2d,2p) and 6-311++G(2d,2p). We observed that 6-31++G(d,p) gives similar results to 6-311++G(2d,2p). Consequently, SVWN and PW91 methods were also used in association with 6-31++G(d,p) to test the influence of the different models of exchange and correlation functionals. A planar structure was obtained for all the optimised trans configuration structures. In both isomers, the presence of the hydroxyl group leads to an asymmetry in certain structural parameters. From these results, two IR or Raman active frequencies can be used to easily distinguish trans and cis configurations. The trans configuration was found to be more stable than the cis configuration by 67 +/- 2 kJ mol(-1) at 0 K. The difference of the dipole moment between trans and cis for 4-hydroxyazobenzene was found to be lower than for trans and cis azobenzene.