The Community Multiscale Air Quality (CMAQ) model is a comprehensive multipollutant air quality modeling system developed and maintained by the US Environmental Protection Agency’s (EPA) Office of Research and Development (ORD). Recently, version 5.1 of the CMAQ model (v5.1) was released to the public, incorporating a large number of science updates and extended capabilities over the previous release version of the model (v5.0.2). These updates include the following: improvements in the meteorological calculations in both CMAQ and the Weather Research and Forecast (WRF) model used to provide meteorological fields to CMAQ, updates to the gas and aerosol chemistry, revisions to the calculations of clouds and photolysis, and improvements to the dry and wet deposition in the model. Sensitivity simulations isolating several of the major updates to the modeling system show that changes to the meteorological calculations result in enhanced afternoon and early evening mixing in the model, periods when the model historically underestimates mixing. This enhanced mixing results in higher ozone (O3) mixing ratios on average due to reduced NO titration, and lower fine particulate matter (PM2.5) concentrations due to greater dilution of primary pollutants (e.g., elemental and organic carbon). Updates to the clouds and photolysis calculations greatly improve consistency between the WRF and CMAQ models and result in generally higher O3 mixing ratios, primarily due to reduced cloudiness and attenuation of photolysis in the model. Updates to the aerosol chemistry result in higher secondary organic aerosol (SOA) concentrations in the summer, thereby reducing summertime PM2.5 bias (PM2.5 is typically underestimated by CMAQ in the summer), while updates to the gas chemistry result in slightly higher O3 and PM2.5 on average in January and July. Overall, the seasonal variation in simulated PM2.5 generally improves in CMAQv5.1 (when considering all model updates), as simulated PM2.5 concentrations decrease in the winter (when PM2.5 is generally overestimated by CMAQ) and increase in the summer (when PM2.5 is generally underestimated by CMAQ). Ozone mixing ratios are higher on average with v5.1 vs. v5.0.2, resulting in higher O3 mean bias, as O3 tends to be overestimated by CMAQ throughout most of the year (especially at locations where the observed O3 is low); however, O3 correlation is largely improved with v5.1. Sensitivity simulations for several hypothetical emission reduction scenarios show that v5.1 tends to be slightly more responsive to reductions in NOx (NO + NO2), VOC and SOx (SO2 + SO4) emissions than v5.0.2, representing an improvement as previous studies have shown CMAQ to underestimate the observed reduction in O3 due to large, widespread reductions in observed emissions.

Simulations of the initial oxidation process of a SiC surface exposed to O2 and H2O molecules was studied with ReaxFF, an atomically detailed reactive molecular dynamics method that naturally models the breaking and forming of bonds. In this work, the ReaxFF forcefield was first expanded by training it with new quantum mechanics data of the binding energy, equation of state, and heat of formation of the SiC crystal, along with data from earlier studies that describes Si – Si, Si – O, and Si – H interactions. This expanded ReaxFF forcefield is capable of simultaneously describing both Si–C–O and Si–O–H bonding interactions. Using the forcefield, oxidation simulations were performed at various temperatures (in the range of 500 to 5000 K), and the trajectories were analyzed. The analyses showed that SiC gradually transforms into the oxides of silicon with simultaneous formation of a graphite-like layer. In presence of excess O2, the graphite-like layer was further oxidized to CO and CO2. We also analyzed the trajectories with two-atom and three-atom clusters to quantitatively track the oxidation progress. This analysis clearly showed Si–O and C–C bond formation at the expense of O–O and Si–C bond consumption indicating SiC oxidation with simultaneous formation for carbon-like structure. Consumption of SiC with O2 was found to be faster than that with H2O. We have also reported the oxidation simulation of SiC with a mixture of H2O and O2. Oxidation proceeded effectively as a two-part sequence, with O2 first oxidizing the SiC, followed then by H2O.

Several models were used to describe the partitioning of ammonia, water, and organic compounds between the gas and particle phases for conditions in the southeastern US during summer 2013. Existing equilibrium models and frameworks were found to be sufficient, although additional improvements in terms of estimating pure-species vapor pressures are needed. Thermodynamic model predictions were consistent, to first order, with a molar ratio of ammonium to sulfate of approximately 1.6 to 1.8 (ratio of ammonium to 2× sulfate, RN/2S ≈ 0.8 to 0.9) with approximately 70% of total ammonia and ammonium (NHx) in the particle. Southeastern Aerosol Research and Characterization Network (SEARCH) gas and aerosol and Southern Oxidant and Aerosol Study (SOAS) Monitor for AeRosols and Gases in Ambient air (MARGA) aerosol measurements were consistent with these conditions. CMAQv5.2 regional chemical transport model predictions did not reflect these conditions due to a factor of 3 overestimate of the nonvolatile cations. In addition, gas-phase ammonia was overestimated in the CMAQ model leading to an even lower fraction of total ammonia in the particle. Chemical Speciation Network (CSN) and aerosol mass spectrometer (AMS) measurements indicated less ammonium per sulfate than SEARCH and MARGA measurements and were inconsistent with thermodynamic model predictions. Organic compounds were predicted to be present to some extent in the same phase as inorganic constituents, modifying their activity and resulting in a decrease in [H+]air (H+ in μgm−3 air), increase in ammonia partitioning to the gas phase, and increase in pH compared to complete organic vs. inorganic liquid–liquid phase separation. In addition, accounting for nonideal mixing modified the pH such that a fully interactive inorganic–organic system had a pH roughly 0.7 units higher than predicted using traditional methods (pH = 1.5 vs. 0.7). Particle-phase interactions of organic and inorganic compounds were found to increase partitioning towards the particle phase (vs. gas phase) for highly oxygenated (O : C≥0.6) compounds including several isoprene-derived tracers as well as levoglu-cosan but decrease particle-phase partitioning for low O: C, monoterpene-derived species.