No abstract

No abstract

To make practical the molecular dynamics simulation of large scale reactive chemical systems (1000s of atoms), we developed ReaxFF, a force field for reactive systems. ReaxFF uses a general relationship between bond distance and bond order on one hand and between bond order and bond energy on the other hand that leads to proper dissociation of bonds to separated atoms. Other valence terms present in the force field (angle and torsion) are defined in terms of the same bond orders so that all these terms go to zero smoothly as bonds break. In addition, ReaxFF has Coulomb and Morse (van der Waals) potentials to describe nonbond interactions between all atoms (no exclusions). These nonbond interactions are shielded at short range so that the Coulomb and van der Waals interactions become constant as R ij f 0. We report here the ReaxFF for hydrocarbons. The parameters were derived from quantum chemical calculations on bond dissociation and reactions of small molecules plus heat of formation and geometry data for a number of stable hydrocarbon compounds. We find that the ReaxFF provides a good description of these data. Generally, the results are of an accuracy similar or better than PM3, while ReaxFF is about 100 times faster. In turn, the PM3 is about 100 times faster than the QC calculations. Thus, with ReaxFF we hope to be able to study complex reactions in hydrocarbons.

We report here an approach for predicting charge distributions in molecules for use in molecular dynamics simulations. The input data are experimental atomic ionization potentials, electron affinities, and atomic radii. An atomic chemical potential is constructed by using these quantities plus shielded electrostatic interactions between all charges. Requiring equal chemical potentials leads to equilibrium charges that depend upon geometry. This charge equilibration (QEq) approach leads to charges in excellent agreement with experimental dipole moments and with the atomic charges obtained from the electrostatic potentials of accurate ab initio calculations. QEq can be used to predict charges for any polymer, ceramic, semiconductor, or biological system, allowing extension of molecular dynamics studies to broad classes of new systems. The charges depend upon environment and change during molecular dynamics calculations. We indicate how this approach can also be used to predict infrared intensities, dielectric constants, and other charge-related properties.

In memoriam Piero PinoStarburst dendrimers are three-dimensional, highly ordered oligomeric and polymeric compounds formed by reiterative reaction sequences starting from smaller molecules-"initiator cores" such as ammonia or pentaerythritol. Protecting group strategies are crucial in these syntheses, which proceed via discrete "Aufbau" stages referred to as generations. Critical molecular design parameters (CMDPs) such as size, shape, and surface chemistry may be controlled by the reactions and synthetic building blocks used. Starburst dendrimers can mimic certain properties of micelles and liposomes and even those of biomolecules and the still more complicated, but highly organized, building blocks of biological systems. Numerous applications of these compounds are conceivable, particularly in mimicking the functions of large biomolecules as drug carriers and immunogens. This new branch of "supramolecular chemistry" should spark new developments in both organic and macromolecular chemistry.

A summary of the technical advances that are incorporated in the fourth major release of the Q-Chem quantum chemistry program is provided, covering approximately the last seven years. These include developments in density functional theory methods and algorithms, nuclear magnetic resonance (NMR) property evaluation, coupled cluster and perturbation theories, methods for electronically excited and openshell species, tools for treating extended environments, algorithms for walking on potential surfaces, analysis tools, energy and electron transfer modelling, parallel computing capabilities, and graphical user interfaces. In addition, a selection of example case studies that illustrate these capabilities is given. These include extensive benchmarks of the comparative accuracy of modern density functionals for bonded and non-bonded interactions, tests of attenuated second order Møller-Plesset (MP2) methods for intermolecular interactions, a variety of parallel performance benchmarks, and tests of the accuracy of implicit solvation models. Some specific chemical examples include calculations on the strongly correlated Cr 2 dimer, exploring zeolitecatalysed ethane dehydrogenation, energy decomposition analysis of a charged ter-molecular complex arising from glycerol photoionisation, and natural transition orbitals for a Frenkel exciton state in a nine-unit model of a self-assembling nanotube.Keywords quantum chemistry, software, electronic structure theory, density functional theory, electron correlation, computational modelling, Q-Chem Disciplines Chemistry CommentsThis article is from Molecular Physics: An International Journal at the Interface Between Chemistry and Physics 113 (2015): 184, doi:10.1080/00268976.2014. RightsWorks produced by employees of the U.S. Government as part of their official duties are not copyrighted within the U.S. The content of this document is not copyrighted. Authors 185A summary of the technical advances that are incorporated in the fourth major release of the Q-CHEM quantum chemistry program is provided, covering approximately the last seven years. These include developments in density functional theory methods and algorithms, nuclear magnetic resonance (NMR) property evaluation, coupled cluster and perturbation theories, methods for electronically excited and open-shell species, tools for treating extended environments, algorithms for walking on potential surfaces, analysis tools, energy and electron transfer modelling, parallel computing capabilities, and graphical user interfaces. In addition, a selection of example case studies that illustrate these capabilities is given. These include extensive benchmarks of the comparative accuracy of modern density functionals for bonded and non-bonded interactions, tests of attenuated second order Møller-Plesset (MP2) methods for intermolecular interactions, a variety of parallel performance benchmarks, and tests of the accuracy of implicit solvation models. Some specific chemical examples include calculations on the strongly corre...

To investigate the initial chemical events associated with high-temperature gas-phase oxidation of hydrocarbons, we have expanded the ReaxFF reactive force field training set to include additional transition states and chemical reactivity of systems relevant to these reactions and optimized the force field parameters against a quantum mechanics (QM)-based training set. To validate the ReaxFF potential obtained after parameter optimization, we performed a range of NVT-MD simulations on various hydrocarbon/O 2 systems. From simulations on methane/O 2 , o-xylene/O 2 , propene/O 2 , and benzene/O 2 mixtures, we found that ReaxFF obtains the correct reactivity trend (propene > o-xylene > methane > benzene), following the trend in the C-H bond strength in these hydrocarbons. We also tracked in detail the reactions during a complete oxidation of isolated methane, propene, and o-xylene to a CO/CO 2 /H 2 O mixture and found that the pathways predicted by ReaxFF are in agreement with chemical intuition and our QM results. We observed that the predominant initiation reaction for oxidation of methane, propene, and o-xylene under fuel lean conditions involved hydrogen abstraction of the methyl hydrogen by molecular oxygen forming hydroperoxyl and hydrocarbon radical species. While under fuel rich conditions with a mixture of these hydrocarbons, we observed different chemistry compared with the oxidation of isolated hydrocarbons including a change in the type of initiation reactions, which involved both decomposition of the hydrocarbon or attack by other radicals in the system. Since ReaxFF is capable of simulating complicated reaction pathways without any preconditioning, we believe that atomistic modeling with ReaxFF provides a useful method for determining the initial events of oxidation of hydrocarbons under extreme conditions and can enhance existing combustion models.

An activity lift for platinum Platinum is an excellent but expensive catalyst for the oxygen reduction reaction (ORR), which is critical for fuel cells. Alloying platinum with other metals can create shells of platinum on cores of less expensive metals, which increases its surface exposure, and compressive strain in the layer can also boost its activity (see the Perspective by Stephens et al. ). Bu et al. produced nanoplates—platinum-lead cores covered with platinum shells—that were in tensile strain. These nanoplates had high and stable ORR activity, which theory suggests arises from the strain optimizing the platinum-oxygen bond strength. Li et al. optimized both the amount of surface-exposed platinum and the specific activity. They made nanowires with a nickel oxide core and a platinum shell, annealed them to the metal alloy, and then leached out the nickel to form a rough surface. The mass activity was about double the best reported values from previous studies. Science , this issue p. 1410 , p. 1414 ; see also p. 1378