Abstract: In this contribution, we present an overview of the various techniques for combining atomistic molecular dynamics with Monte Carlo simulations, mainly in the context of condensed matter systems, as well as a brief summary of the main accelerated dynamics techniques. Special attention is given to the force bias Monte Carlo technique and its combination with molecular dynamics, in view of promising recent developments, including a definable timescale. Various examples of the application of combined molecular dynamics / Monte Carlo simulations are given, in order to demonstrate the enhanced simulation efficiency with respect to either pure molecular dynamics or Monte Carlo.

Abstract: The behavior of hydrocarbon species at step edges of diamond terraces is investigated by means of combined molecular dynamics−Metropolis Monte Carlo simulations. The results show that the formation of ballas-like diamond films (like UNCD) and well-faceted diamond films (like NCD) can be related to the gas phase concentrations of CxHy in a new manner: Species that have high concentrations above the growing UNCD films suppress the extension of step edges through defect formation. The species that are present above the growing NCD film, however, enhance the extension of diamond terraces, which is believed to result in well-faceted diamond films. Furthermore, it is shown that, during UNCD growth, CxHy species with x ≥ 2 play an important role, in contrast to the currently adopted CVD diamond growth mechanism. Finally, the probabilities for the extension of the diamond (100) terrace are much higher than those for the diamond (111) terrace, which is in full agreement with the experimental observation that diamond (100) facets are more favored than diamond (111) facets during CVD diamond growth.

Abstract: Computer simulations are performed for an argon inductively coupled plasma (ICP) with a capacitive radio-frequency bias power, to investigate the bias effect on the discharge mode transition and on the plasma characteristics at various ICP currents, bias voltages, and bias frequencies. When the bias frequency is fixed at 13.56 MHz and the ICP current is low, e.g., 6A, the spatiotemporal averaged plasma density increases monotonically with bias voltage, and the bias effect is already prominent at a bias voltage of 90 V. The maximum of the ionization rate moves toward the bottom electrode, which indicates clearly the discharge mode transition in inductive/capacitive discharges. At higher ICP currents, i.e., 11 and 13 A, the plasma density decreases first and then increases with bias voltage, due to the competing mechanisms between the ion acceleration power dissipation and the capacitive power deposition. At 11 A, the bias effect is still important, but it is noticeable only at higher bias voltages. At 13 A, the ionization rate is characterized by a maximum at the reactor center near the dielectric window at all selected bias voltages, which indicates that the ICP power, instead of the bias power, plays a dominant role under this condition, and no mode transition is observed. Indeed, the ratio of the bias power to the total power is lower than 0.4 over a wide range of bias voltages, i.e., 0300V. Besides the effect of ICP current, also the effect of various bias frequencies is investigated. It is found that the modulation of the bias power to the spatiotemporal distributions of the ionization rate at 2MHz is strikingly different from the behavior observed at higher bias frequencies. Furthermore, the minimum of the plasma density appears at different bias voltages, i.e., 120V at 2MHz and 90V at 27.12 MHz.

Abstract: The mechanism of interaction of cold nonequilibrium plasma jets with mammalian cells in physiologic liquid is reported. The major biological active species produced by an argon RF plasma jet responsible for cell viability reduction are analyzed by experimental results obtained through physical, biological, and chemical diagnostics. This is complemented with chemical kinetics modeling of the plasma source to assess the dominant reactive gas phase species. Different plasma chemistries are obtained by changing the feed gas composition of the cold argon based RF plasma jet from argon, humidified argon (0.27%), to argon/oxygen (1%) and argon/air (1%) at constant power. A minimal consensus physiologic liquid was used, providing isotonic and isohydric conditions and nutrients but is devoid of scavengers or serum constituents. While argon and humidified argon plasma led to the creation of hydrogen peroxide dominated action on the mammalian cells, argonoxygen and argonair plasma created a very different biological action and was characterized by trace amounts of hydrogen peroxide only. In particular, for the argonoxygen (1%), the authors observed a strong negative effect on mammalian cell proliferation and metabolism. This effect was distance dependent and showed a half life time of 30 min in a scavenger free physiologic buffer. Neither catalase and mannitol nor superoxide dismutase could rescue the cell proliferation rate. The strong distance dependency of the effect as well as the low water solubility rules out a major role for ozone and singlet oxygen but suggests a dominant role of atomic oxygen. Experimental results suggest that O reacts with chloride, yielding Cl2 − or ClO−. These chlorine species have a limited lifetime under physiologic conditions and therefore show a strong time dependent biological activity. The outcomes are compared with an argon MHz plasma jet (kinpen) to assess the differences between these (at least seemingly) similar plasma sources.

Abstract: In this paper, a new implementation of the Metropolis Monte Carlo (MMC) algorithm is presented. When combining the MMC model with a molecular dynamics (MD) code, crystal growth by plasma-enhanced chemical vapor deposition can be simulated. As the MD part simulates impacts of growth species onto the surface on a time scale of picoseconds, the MMC algorithm simulates the slower adatom surface processes. The implementation includes a criterion for the selection of atoms that are allowed to be displaced during the simulation, and a criterion of after how many MMC cycles the simulation is stopped. We performed combined MD-MMC simulations for hydrocarbon species that are important for the growth of ultrananocrystalline diamond (UNCD) films at partially hydrogenated diamond surfaces, since this implementation is part of a study of the growth mechanisms of (ultra)nanocrystalline diamond films. Exemplary for adatom arrangements during the growth of UNCD, the adatom surface behavior of C and C2H2 at diamond (111)1 × 1, C and C4H2 at diamond (111)1 × 1 and C3 at diamond (100)2 × 1 has been investigated. For all cases, the diamond crystal structure is pursued under the influence of MMC simulation. Additional longer time-scale MD simulations put forward very similar structures, verifying the MMC algorithm. Nevertheless, the MMC simulation time is typically one order of magnitude shorter than the MD simulation time.

Abstract: Dry reforming of methane has gained significant interest over the years. A novel reforming technique with great potential is plasma technology. One of its drawbacks, however, is energy consumption. Therefore, we performed an extensive computational study, supported by experiments, aiming to identify the influence of the operating parameters (gas mixture, power, residence time and frequency) of a dielectric barrier discharge plasma on the conversion and energy efficiency, and to investigate which of these parameters lead to the most promising results and whether these are eventually sufficient for industrial implementation. The best results, in terms of both energy efficiency and conversion, are obtained at a specific energy input (SEI) of 100 J cm−3, a 1090 CH4CO2 ratio, 10 Hz, a residence time of 1 ms, resulting in a total conversion of 84% and an energy efficiency of 8.5%. In general, increasing the CO2 content in the gas mixture leads to a higher conversion and energy efficiency. The SEI couples the effect of the power and residence time, and increasing the SEI always results in a higher conversion, but somewhat lower energy efficiencies. The effect of the frequency is more complicated: we observed that the product of frequency (f) and residence time (τ), being a measure for the total number of micro-discharge filaments which the gas molecules experience when passing through the reactor, was critical. For most cases, a higher number of filaments yields higher values for conversion and energy efficiency. To benchmark our model predictions, we also give an overview of measured conversions and energy efficiencies reported in the literature, to indicate the potential for improvement compared to the state-of-the art. Finally, we identify the limitations as well as the benefits and future possibilities of plasma technology.

Abstract: The tumor suppressor protein p53 can lose its function upon DNA-contact mutations (R273C and R273H) in the core DNA-binding domain. The activity can be restored by second-site suppressor or rescue mutations (R273CT284R, R273HT284R, and R273HS240R). In this paper, we elucidate the structural and functional consequence of p53 proteins upon DNA-contact mutations and rescue mutations and the underlying mechanisms at the atomic level by means of molecular dynamics simulations. Furthermore, we also apply the docking approach to investigate the binding phenomena between the p53 protein and DNA upon DNA-contact mutations and rescue mutations. This study clearly illustrates that, due to DNA-contact mutants, the p53 structure loses its stability and becomes more rigid than the native protein. This structural loss might affect the p53-DNA interaction and leads to inhibition of the cancer suppression. Rescue mutants (R273CT284R, R273HT284R and R273HS240R) can restore the functional activity of the p53 protein upon DNA-contact mutations and show a good interaction between the p53 protein and a DNA molecule, which may lead to reactivate the cancer suppression function. Understanding the effects of p53 cancer and rescue mutations at the molecular level will be helpful for designing drugs for p53 associated cancer diseases. These drugs should be designed so that they can help to inhibit the abnormal function of the p53 protein and to reactivate the p53 function (cell apoptosis) to treat human cancer.

Abstract: It is known that film structure may change dramatically with the extent of surface diffusion during the film growth process. In the present work, surface diffusion, induced thermally or activated by energetic impacts, is investigated theoretically under conditions appropriate for magnetron sputter-deposition of MgAlO thin films with varying stoichiometry. The distribution of surface diffusion energy barriers available to the system was determined for each stoichiometry, which allowed assessing in a qualitative way how much surface diffusion will take place on the time scale available between deposition events. The activation energy barriers increase with the Al concentration in the film, and therefore, the surface diffusion rates in the time frame of typical deposition rates drop, which can explain the decrease in crystallinity in the film structure and the transition to amorphous structure. The deposition process and the immediate surface diffusion enhanced by the energetic adatoms are simulated by means of a molecular dynamics model. The longer-time thermal surface diffusion and the energy landscape are studied by the temperature accelerated dynamics method, applied in an approximate way. The surface diffusion enhanced by the energetic impacts appears to be very important for the film structure in the low-temperature deposition regime.

Abstract: We studied the influence of dense, spherical packing materials, with different chemical compositions, on the dry reforming of methane (DRM) in a dielectric barrier discharge (DBD) reactor. Although not catalytically activated, a vast effect on the conversion and product selectivity could already be observed, an influence which is often neglected when catalytically activated plasma packing materials are being studied. The alpha-Al2O3 packing material of 2.0-2.24 mm size yields the highest total conversion (28%), as well as CO2 (23%) and CH4 (33%) conversion and a high product fraction towards CO (similar to 70%) and ethane (similar to 14%), together with an enhanced CO/H-2 ratio of 9 in a 4.5 mm gap DBD at 60 W and 23 kHz. gamma-Al2O3 is only slightly less active in total conversion (22%) but is even more selective in products formed than alpha-Al2O3 BaTiO3 produces substantially more oxygenated products than the other packing materials but is the least selective in product fractions and has a clear negative impact on CO2 conversion upon addition of CH4. Interestingly, when comparing to pure CO2 splitting and when evaluating differences in products formed, significantly different trends are obtained for the packing materials, indicating a complex impact of the presence of CH4 and the specific nature of the packing materials on the DRM process.