Simulating Polymerization by Boltzmann Inversion Force Field Approach and Dynamical Nonequilibrium Reactive Molecular Dynamics.

The radical polymerization process of acrylate compounds is, nowadays, numerically investigated using classical force fields and reactive molecular dynamics, with the aim to probe the gel-point transition as a function of the initial radical concentration. In the present paper, the gel-point transition of the 1,6-hexanediol dimethacrylate (HDDMA) is investigated by a coarser force field which grants a reduction in the computational costs, thereby allowing the simulation of larger system sizes and smaller radical concentrations. Hence, the polymerization is investigated using reactive classical molecular dynamics combined with a dynamical approach of the nonequilibrium molecular dynamics (D-NEMD). The network structures in the polymerization process are probed by cluster analysis tools, and the results are critically compared with the similar all-atom system, showing a good agreement.

Capturing Free-Radical Polymerization by Synergetic Ab Initio Calculations and Topological Reactive Molecular Dynamics.

Photocurable polymers are used ubiquitously in 3D printing, coatings, adhesives, and composite fillers. In the present work, the free radical polymerization of photocurable compounds is studied using reactive classical molecular dynamics combined with a dynamical approach of the nonequilibrium molecular dynamics (D-NEMD). Different concentrations of radicals and reaction velocities are considered. The mechanical properties of the polymer resulting from 1,6-hexanediol dimethacrylate systems are characterized in terms of viscosity, diffusion constant, and activation energy, whereas the topological ones through the number of cycles (polymer loops) and cyclomatic complexity. Effects like volume shrinkage and delaying of the gel point for increasing monomer concentration are also predicted, as well as the stress-strain curve and Young's modulus. Combining ab initio, reactive molecular dynamics, and the D-NEMD method might lead to a novel and powerful tool to describe photopolymerization processes and to original routes to optimize additive manufacturing methods relying on photosensitive macromolecular systems.

Controlling electroosmotic flows by polymer coatings: A joint experimental-theoretical investigation.

We analyze the electroosmotic flow (EOF) of an electrolytic solution in a polymer coated capillary electrophoresis tube. The polymeric density, charge, thickness, and the capillary tube charge vary as a function of pH and produce a non-trivial modulation of the EOF, including a flow reversal at acid pH conditions. By means of a theoretical argument and numerical simulations, we recover the experimental curve for the EOF, providing a firm approach for predictive analysis of electroosmosis under different polymeric coating conditions. A proposed application of the approach is to determine the near-wall charge of the coating to be used for further quantitative analysis of the electroosmotic flow and mobility.

Electro-osmotic flow in coated nanocapillaries: a theoretical investigation.

Motivated by recent experiments, we present a theoretical investigation of how the electro-osmotic flow occurring in a capillary is modified when its charged surfaces are coated with charged polymers. The theoretical treatment is based on a three-dimensional model consisting of a ternary fluid-mixture, representing the solvent and two species for the ions, confined between two parallel charged plates decorated with a fixed array of scatterers representing the polymer coating. The electro-osmotic flow, generated by a constant electric field applied in a direction parallel to the plates, is studied numerically by means of Lattice Boltzmann simulations. In order to gain further understanding we performed a simple theoretical analysis by extending the Stokes-Smoluchowski equation to take into account the porosity induced by the polymers in the region adjacent to the walls. We discuss the nature of the velocity profiles by focusing on the competing effects of the polymer charges and the frictional forces they exert. We show evidence of the flow reduction and of the flow inversion phenomenon when the polymer charge is opposite to the surface charge. By using the density of polymers and the surface charge as control variables, we propose a phase diagram that discriminates the direct and the reversed flow regimes and determines their dependence on the ionic concentration.

Lattice Boltzmann method for mixtures at variable Schmidt number.

When simulating multicomponent mixtures via the Lattice Boltzmann Method, it is desirable to control the mutual diffusivity between species while maintaining the viscosity of the solution fixed. This goal is herein achieved by a modification of the multicomponent Bhatnagar-Gross-Krook evolution equations by introducing two different timescales for mass and momentum diffusion. Diffusivity is thus controlled by an effective drag force acting between species. Numerical simulations confirm the accuracy of the method for neutral binary and charged ternary mixtures in bulk conditions. The simulation of a charged mixture in a charged slit channel show that the conductivity and electro-osmotic mobility exhibit a departure from the Helmholtz-Smoluchowski prediction at high diffusivity.

Quantum dynamical structure factor of liquid neon via a quasiclassical symmetrized method.

We apply the phase integration method for quasiclassical quantum time correlation functions [M. Monteferrante, S. Bonella, and G. Ciccotti, Mol. Phys. 109, 3015 (2011)] to compute the dynamic structure factor of liquid neon. So far the method had been tested only on model systems. By comparing our results for neon with experiments and previous calculations, we demonstrate that the scheme is accurate and efficient also for a realistic model of a condensed phase system showing quantum behavior.

The quantum free energy barrier for hydrogen vacancy diffusion in Na3AlH6.

The path integral single sweep method is used to assess quantum effects on the free energy barrier for hydrogen vacancy diffusion in a defective Na(3)AlH(6) crystal. This process has been investigated via experiments and simulations due to its potential relevance in the H release mechanism in sodium alanates, prototypical materials for solid state hydrogen storage. Previous computational studies, which used density functional methods for the electronic structure, were restricted to a classical treatment of the nuclear degrees of freedom. We show that, although they do not change the qualitative picture of the process, nuclear quantum effects reduce the free energy barrier height by about 18% with respect to the classical calculation improving agreement with available neutron scattering data.

Short range hydrogen diffusion in Na3AlH6.

Ab initio free energy and rate calculations are performed to investigate two activated mobility processes observed, respectively, in neutron scattering and anelastic spectroscopy experiments on sodium alanates. The system is modeled as a Na(3)AlH(6) crystal hosting one hydrogen vacancy. We identify the process observed via neutron scattering with a positively charged hydrogen vacancy diffusing from the AlH to one of the AlH groups. As for the anelastic spectroscopy experiments, our calculations negate the current hypothesis on the process, i.e. local rearrangement of the H vacancy around the pentacoordinated Al group.