Interface motion in a two-dimensional Ising model with a field

We determine by Monte Carlo simulations the width of an interface between the stable phase and the metastable phase in a two-dimensional Ising model with a magnetic field, in the case of nonconversed order parameter (Glauber dynamics). At zero temperature, the width increases ast with–1/3, as predicted by earlier theories. As temperature increases, the value of the effective exponent that we measure decreases toward the value 1/4, which is the value in the absence of magnetic field.     

Metastable systems driven by colored noise: The stationary state

The influence of the Bardeen-Herring back-jump correlations on the Fermi-Dirac statistics of the one-dimensional nonhomogeneous fermionic lattice gas is studied by the Monte Carlo simulation technique and semianalytically. The resulting distribution is obtained, exhibiting increased population of the lower levels in comparison to the Fermi-Dirac statistics.     

Monte Carlo simulation of phase separation and clustering in the ABV model

As a model for a binary alloy undergoing an unmixing phase transition, we consider a square lattice where each site can be either taken by an A atom, a B atom, or a vacancy (V), and there exists a repulsive interaction between AB nearest neighbor pairs. Starting from a random initial configuration, unmixing proceeds via random jumps of A atoms or B atoms to nearest neighbor vacant sites. In the absence of any interaction, these jumps occur at jump rates _A and _B, respectively. For a small concentration of vacancies (c _v=0.04) the dynamics of the structure factorS(k,t) and its first two momentsk ₁(t),k ₂ ² (t) is studied during the early stages of phase separation, for several choices of concentrationc _B of B atoms. Forc _B=0.18 also the time evolution of the cluster size distribution is studied. Apart from very early times, the mean cluster sizel(t) as well as the moments of the structure function depend on timet and the ratio of the jump rates (= _B/ _A) only via a scaled timet/(). Qualitatively, the behavior is very similar to the direct exchange model containing no vacancies. Consequences for phase separation of real alloys are briefly discussed.     

Monte carlo calculations of the conformal charge

The conformal charge is an important quantity which characterizes the nature of the two-dimensional phase transition. We report a first attempt to use a new numerical method to calculate the conformal charge. In this paper, we apply our method to the 2-dimensional, ⁴, continuous-spin Ising model. By varying the parameters in the Hamiltonian, one can change continuously from the known Gaussian limit to the Ising limit. It is well known that the critical points for these two systems are not in the same universality class. We study this behavior for the Gaussian model, the single-well ⁴ model, the border model, and the double-well ⁴ model for a large lattice. Our results, while giving a good general picture, are not so far sufficient to differentiate whether the non-Gaussian cases studied belong to the Ising model universality class or not. Further studies of other lattice sizes should serve to improve greatly our conclusions.     

Monte Carlo simulation of spatial distribution of X rays in multifilm targets. III: optimum design as applied to W/Al‐film targets

The brightness of the X‐ray source in a W/Al‐film target used for X‐ray projection microscopy was studied by an approach using the Monte Carlo simulation. Since continuous X rays generated in a thin film have a specific angular distribution of emission, the brightness of the continuous X‐ray source cannot simply be estimated on the assumption that the angular distribution is homogeneous. The newly developed approach using the Monte Carlo simulation enables the evaluation of the effective source size, angular distribution, and brightness of a continuous X‐ray source with sufficient accuracy that it leads to the optimum design of a high‐brightness X‐ray source for uses such as X‐ray projection microscopy. The Monte Carlo calculations were performed for W(Δz)/Al (200 µm)‐film targets with different thicknesses of W film, Δz, under bombardment of 60 kV electrons. The results have suggested an optimum design consisting of a W (2 µm)/Al (200 µm)‐film target as most promising for providing an X‐ray source of higher brightness than the W (5 µm)/Al (200 µm)‐film target, which has already been in practice for X‐ray projection microscopy. Copyright © 2006 John Wiley & Sons, Ltd.     

Monte Carlo simulation of a lattice model for ternary polymer mixtures

Monte Carlo studies of symmetrical polymer mixturesAB, modelled by selfavoiding walks withN _A =N _B =N steps on a simple cubic lattice, are presented for arbitrary concentrations of vacancies _v in the range from _v =0.2 to _v =0.8 and chain lengthsN64. We obtained the phase diagrams and the equation of state for three choices of the ratio / _AB ( being the energy between monomers of the same kind, _AB being the energy between different monomers). Flory-Huggins theory provides only a qualitative understanding of these results. If the equation of state is fitted with an effective Flory-Huggins parameter _eff , the latter turns out to be strongly dependent on both concentration and temperature.Contributed paper delivered at the Tagung der Deutschen Physikalischen Gesellschaft, Fachausschuß Polymerphysik, Berlin, March 30–April 3, 1987.     

张量积型Said-Ball曲面的预处理渐近迭代逼近法

为加快张量积型Said-Ball曲面渐近迭代逼近法的收敛速度,探讨了张量积型Said-Ball曲面渐近迭代逼近法的预处理技术。首先利用对角补偿约化技术构造了预处理子,然后结合矩阵Kronecker积性质,采取预处理渐近迭代逼近法求解张量积型Said-Ball曲面。为进一步降低计算量并提高算法的稳定性,利用广义极小残差法求解预处理方程,得到预处理渐近迭代逼近法的非精确求解方法。分析了预处理渐近迭代逼近法及非精确求解方法的收敛性。最后用数值实例说明预处理子能大大减小迭代矩阵的谱半径,令预处理技术及其非精确求解方法的计算效率明显提高。此外,由于对角补偿预处理子能改善配置矩阵的谱分布,因此也可用于对广义极小残差法的预处理,以改善其收敛性。     

MLIMC: Machine Learning-Based Implicit-Solvent Monte Carlo

Monte Carlo (MC) methods are important computational tools for molecular structure optimizations and predictions. When solvent effects are explicitly considered, MC methods become very expensive due to the large degree of freedom associated with the water molecules and mobile ions. Alternatively implicit-solvent MC can largely reduce the computational cost by applying a mean field approximation to solvent effects and meanwhile maintains the atomic detail of the target molecule. The two most popular implicit-solvent models are the Poisson-Boltzmann (PB) model and the Generalized Born (GB) model in a way such that the GB model is an approximation to the PB model but is much faster in simulation time. In this work, we develop a machine learning-based implicit-solvent Monte Carlo (MLIMC) method by combining the advantages of both implicit solvent models in accuracy and efficiency. Specifically, the MLIMC method uses a fast and accurate PB-based machine learning (PBML) scheme to compute the electrostatic solvation free energy at each step. We validate our MLIMC method by using a benzene-water system and a protein-water system. We show that the proposed MLIMC method has great advantages in speed and accuracy for molecular structure optimization and prediction.     

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.     

Spectral Analysis for HSS Preconditioners

In this paper,we are interested in HSS preconditioners for saddle point lin- ear systems with a nonzero(2,2)-th block.We study an approximation of the spectra of HSS preconditioned matrices and use these results to illustrate and explain the spectra obtained from numerical examples,where the previous spectral analysis of HSS precon- ditioned matrices does not cover.