Methods for implementing the stream boundary condition in DSMC computations

In the direct simulation Monte‐Carlo (DSMC) method for simulating rarefied gas flows, the velocities of simulator particles that cross a simulation boundary and enter the simulation space are typically generated using the acceptance–rejection procedure that samples the velocities from a truncated theoretical velocity distribution that excludes low and high velocities. This paper analyses an alternative technique, where the velocities of entering particles are obtained by extending the simulation procedures to a region adjacent to the simulation space, and considering the movement of particles generated within that region during the simulation time step. The alternative method may be considered as a form of acceptance–rejection procedure, and permits the generation of all possible velocities, although the population of high velocities is depleted with respect to the theoretical distribution. Nevertheless, this is an improvement over the standard acceptance–rejection method. Previous implementations of the alternative method gave a number flux lower than the theoretical number required. Two methods for obtaining the correct number flux are presented. For upstream boundaries in high‐speed flows, the alternative method is more computationally efficient than the acceptance–rejection method. However, for downstream boundaries, the alternative method is extremely inefficient. The alternative method, with the correct theoretical number flux, should therefore be used in DSMC computations in favour of the acceptance–rejection method for upstream boundaries in high‐speed flows. Copyright © 2003 John Wiley & Sons, Ltd.     

Monte Carlo Hamiltonian: Linear Potentials

We further study the validity of the Monte Carlo Hamiltonian method. The advantage of the method,in comparison with the standard Monte Carlo Lagrangian approach, is its capability to study the excited states. Weconsider two quantum mechanical models: a symmetric one V(x) = |x|/2; and an asymmetric one V(x) = ∞, forx ＜ 0 and V(x) = x, for x ≥ 0. The results for the spectrum, wave functions and thermodynamical observables are inagreement with the analytical or Runge-Kutta calculations.     

The phase diagram of a single polymer chain: New insights from a new simulation method

We present simulation results for the phase behavior of a single chain for a flexible lattice polymer model using the Wang-Landau sampling idea. Applying this new algorithm to the problem of the homopolymer collapse allows us to investigate not only the high temperature coil–globule transition but also an ensuing crystallization at lower temperature. Performing a finite size scaling analysis on the two transitions, we show that they coincide for our model in the thermodynamic limit corresponding to a direct collapse of the random coil into the crystal without intermediate coil–globule transition. As a consequence, also the many chain phase diagram of this model can be predicted to consist only of gas and crystal phase in the limit of infinite chain length. This behavior is in agreement with findings on the phase behavior of hard-sphere systems with a relatively short-ranged attractive square well. © 2006 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 44: 2542–2555, 2006     

New Monte Carlo method for the self-avoiding walk

We introduce a new Monte Carlo algorithm for the self-avoiding walk (SAW), and show that it is particularly efficient in the critical region (long chains). We also introduce new and more efficient statistical techniques. We employ these methods to extract numerical estimates for the critical parameters of the SAW on the square lattice. We find=2.63820 ± 0.00004 ± 0.00030=1.352 ± 0.006 ± 0.025v=0.7590 ± 0.0062 ± 0.0042 where the first error bar represents systematic error due to corrections to scaling (subjective 95% confidence limits) and the second bar represents statistical error (classical 95% confidence limits). These results are based on SAWs of average length 166, using 340 hours CPU time on a CDC Cyber 170–730. We compare our results to previous work and indicate some directions for future research.     

Dynamics of the formation of two-dimensional ordered structures

The growth of ordered domains in lattice gas models, which occurs after the system is quenched from infinite temperature to a state below the critical temperatureT _c, is studied by Monte Carlo simulation. For a square lattice with repulsion between nearest and next-nearest neighbors, which in equilibrium exhibits fourfold degenerate (2×1) superstructures, the time-dependent energy E(t), domain size L(t), and structure functionS(q, t) are obtained, both for Glauber dynamics (no conservation law) and the case with conserved density (Kawasaki dynamics). At late times the energy excess and halfwidth of the structure factor decrease proportional tot ^–x, whileL(t) t ^x, where the exponent x=1/2 for Glauber dynamics and x1/3 for Kawasaki dynamics. In addition, the structure factor satisfies a scaling lawS(k,t)=t ^2xS(kt^x). The smaller exponent for the conserved density case is traced back to the excess density contained in the walls between ordered domains which must be redistributed during growth. Quenches toT>T _c, T=T_c (where we estimate dynamic critical exponents) andT=0 are also considered. In the latter case, the system becomes frozen in a glasslike domain pattern far from equilibrium when using Kawasaki dynamics. The generalization of our results to other lattices and structures also is briefly discussed.     

Kinetics of CO oxidation on the surface of heterophase catalysts. Modeling by the Monte Carlo method

The Monte Carlo method has been used to simulate CO oxidation on a lattice consisting of various alternating patches: M1, where s(CO)>s(O₂) and M2, where s(CO)2). The reaction is shown to proceed over all the surface at low temperature as CO_ads spillover from M1 to M2 and backwards.     

Investigation of a cylindrical chemosorptive denuder for sampling and phase separation of toluene diisocyanate aerosols

A cylindrical chemosorptive denuder in series with a glass fibre filter has been evaluated for sampling toluene diisocyanate (TDI) aerosols. The sampler is designed for measuring personal exposure to diisocyanates. Several denuder coatings and derivatising reagents were investigated. Dimethylpolysiloxane (SE-30) and 5% phenyl dimethylpolysiloxane (SE-54) with either dibutylamine (DBA) or dipentylamine (DPeA) as derivatising reagents yielded the lowest vapour breakthrough (the amount (%) of the vapour that passes through the denuder), close to values predicted by theory. Immobilisation of the SE-30 denuder coating by in-situ cross-linking yielded comparable results. With an SE-30/DBA-coated denuder operating within an airflow range of 100–500 mL min^–1, the phase separation was shown to be consistent with theoretical predictions derived by use of the Gormley–Kennedy equation. This provides a means of calculating the vapour breakthrough and correcting experimentally obtained values with regard to vapour–particulate phase distribution, suggesting that the denuder can provide accurate phase-distribution measurements. The SE-30/DBA denuder can be used over a concentration range spanning nearly six orders of magnitude. Its capacity is sufficient to perform 15-min exposure measurements of a TDI aerosol with air concentrations as high as 1,700 g m^–3, 40 times higher than the Swedish occupational exposure limit (OEL). At the other end of the range, the estimated limit of detection (LOD) was less than 2 ng m^–3 for both the vapour and the aerosol phases when LC–ESI–MS–MS was used for chemical analysis.Electronic Supplementary Material Supplementary material for this article is available at     

Numerical Study of φ4 Model by Potential Importance Sampling Method

We investigate the phenomena of spontaneous symmetry breaking for φ⁴ model on a square lattice in the parameter space by using the potential importance sampling method, which was proposed by Milchev, Heermann, and Binder [J. Stat. Phys. 44 (1986) 749]. The critical values of the parameters allow us to determine the phase diagram of the model. At the same time, some relevant quantities such as susceptibility and specific heat are also obtained.     

Surface Modification Strategy for Enhanced NO2 Capture in Metal–Organic Frameworks

The interaction strength of nitrogen dioxide (NO₂) with a set of 43 functionalized benzene molecules was investigated by performing density functional theory (DFT) calculations. The functional groups under study were strategically selected as potential modifications of the organic linker of existing metal–organic frameworks (MOFs) in order to enhance their uptake of NO₂ molecules. Among the functional groups considered, the highest interaction energy with NO₂ (5.4 kcal/mol) was found for phenyl hydrogen sulfate (-OSO₃H) at the RI-DSD-BLYP/def2-TZVPP level of theory—an interaction almost three times larger than the corresponding binding energy for non-functionalized benzene (2.0 kcal/mol). The groups with the strongest NO₂ interactions (-OSO₃H, -PO₃H₂, -OPO₃H₂) were selected for functionalizing the linker of IRMOF-8 and investigating the trend in their NO₂ uptake capacities with grand canonical Monte Carlo (GCMC) simulations at ambient temperature for a wide pressure range. The predicted isotherms show a profound enhancement of the NO₂ uptake with the introduction of the strongly-binding functional groups in the framework, rendering them promising modification candidates for improving the NO₂ uptake performance not only in MOFs but also in various other porous materials.     

Use Of Diffusion Controlled Preconcentration Systems With Chemically Immobilized Coatings For The Determination Of Atmospheric Nitric Acid

Abstract

A new coating procedure for different diffusion controlled preconcentration methods (tubular denuder, annular denuder, diffusion screen) is described for the determination of nitric acid in air. In this study, a silanization reaction is applied to obtain a chemically fixed coating with an end placed functional NH₂-group, which can be used as a sink to collect acidic compounds from air. This coating must be carried out only once and can be used for a long time. The denuders, coated with this procedures, were compared with conventionally coated NaF denuders.