Field-theoretic analysis of directed percolation: Three-loop approximation.

The directed bond percolation is a paradigmatic model in nonequilibrium statistical physics. It captures essential physical information on the nature of continuous phase transition between active and absorbing states. In this paper, we study this model by means of the field-theoretic formulation with a subsequent renormalization group analysis. We calculate all critical exponents needed for the quantitative description of the corresponding universality class to the third order in perturbation theory. Using dimensional regularization with minimal subtraction scheme, we carry out perturbative calculations in a formally small parameter ɛ, where ɛ=4-d is a deviation from the upper critical dimension d_{c}=4. We use a nontrivial combination of analytical and numerical tools in order to determine ultraviolet divergent parts of Feynman diagrams.

Surface Tension of Organophosphorus Compounds: Sarin and its Surrogates.

While the production and stockpiling of organophosphorus chemical warfare agents (CWAs), such as sarin, was banned three decades ago, CWAs have remained a threat. New approaches for decontamination and destruction of CWAs require detailed knowledge of their various physicochemical properties. In particular, surface tension is needed to describe the formation and evolution of hazardous aerosols when CWA liquids are dispersed in the air. Due to the extreme toxicity of sarin, most experimental studies are carried out using its surrogates─organophosphorus compounds which, while having similar structures, are much less toxic, e.g., dimethyl methylphosphonate (DMMP) and diisopropyl methylphosphonate (DIMP). However, not only for sarin, but also for its surrogates, literature data on the surface tension are scarce. Here we present experimental measurements and computational predictions of the surface tension of DMMP and DIMP. Classical molecular dynamics simulations using the Transferable Potentials for Phase Equilibria (TraPPE) force field produced an excellent agreement with the experimental results in the temperature range from 3 to 60 °C, validating the predictive capability of TraPPE. Consequently, we applied the TraPPE force field to sarin. Our modeled values for the sarin surface tension cover the range of temperatures from 0 to 85 °C, and the four experimental data points from the literature measured between 20 and 35 °C agree perfectly with our predictions. The temperature-dependent surface tension values for sarin and its surrogates obtained in our study can be used in models predicting the formation and evolution of aerosols made of these chemicals. Furthermore, our results justify the use of the TraPPE force field to derive the thermodynamic properties of other organophosphorus compounds with structures similar to the ones studied here.

Force Fields for Molecular Modeling of Sarin and its Simulants: DMMP and DIMP.

Even three decades after signing the Chemical Weapons Convention, organophosphorus chemical warfare agents (CWAs), such as sarin, remain a threat. The development of novel methods for the detection of CWAs, protection from CWAs, and CWA decontamination motivates research on their physicochemical properties. Due to the extreme toxicity of sarin, most of the experimental studies are carried out using less toxic simulant compounds. In addition to experimental studies of sarin simulants, both sarin and simulants can be studied using in silico experiments-molecular simulations. The results of classical molecular modeling of the compounds and their agreement with experimental data rely on the force field used to describe the system. In recent years, there have been several force fields proposed for sarin and its most common simulant dimethyl methylphosphonate (DMMP). However, other simulants frequently used in experiments received less attention from the molecular simulation perspective, for example, to date, there is no force field and no simulation data for diisopropyl methylphosphonate (DIMP). Here, we compare the literature force fields for sarin and DMMP, focusing specifically on the vapor-liquid equilibrium for the pure species. We carried out Monte Carlo and molecular dynamics simulations using the existing literature force fields from which we predicted the liquid densities and vapor pressures developing the entire binodal curves. We compared the predictions to the experimental data and showed that the TraPPE-UA force field outperformed the other force fields. Thus, we extended TraPPE-UA for DIMP, utilized this force field in molecular simulations, and predicted the liquid densities and vapor pressures for a range of temperatures (binodal curve), which agreed well with the published experimental data. From the binodal, we calculated the critical properties of DIMP and demonstrated that these parameters can be used in the Peng-Robinson equation of state for this compound.

Single Parameter for Predicting the Morphology of Atmospheric Black Carbon.

Black carbon (BC) from fuel combustion is an effective light absorber that contributes significantly to direct climate forcing. The forcing is altered when BC combines with other substances, which modify its mixing state and morphology, making the evaluation of its atmospheric lifetime and climate impact a challenge. To elucidate the associated mechanisms, we exposed BC aerosol to supersaturated vapors of different chemicals to form thin coatings and measured the coating mass required to induce the restructuring of BC aggregates. We found that studied chemicals fall into two distinct groups based on a single dimensionless parameter, χ, which depends on the diameter of BC monomer spheres and the coating material properties, including vapor supersaturation, molar volume, and surface tension. We show that when χ is small (low-volatility chemicals), the highly supersaturated vapor condenses uniformly over aggregates, including convex monomers and concave junctions in between monomers, but when χ is large (intermediate-volatility chemicals), junctions become preferred. The aggregates undergo prompt restructuring when condensation in the junctions dominates over condensation on monomer spheres. For a given monomer diameter, the coating distribution is mostly controlled by vapor supersaturation. The χ factor can be incorporated straightforwardly into atmospheric models to improve simulations of BC aging.