Bal K (2018) New ways to bridge the gap between microscopic simulations and macroscopic chemistry. Antwerpen
Keywords: Doctoral thesis; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
|
“Establishing uniform acceptance in force biased Monte Carlo simulations”. Neyts EC, Thijsse BJ, Mees MJ, Bal KM, Pourtois G, Journal of chemical theory and computation 8, 1865 (2012). http://doi.org/10.1021/ct2008268
Abstract: Uniform acceptance force biased Monte Carlo (UFMC) simulations have previously been shown to be a powerful tool to simulate atomic scale processes, enabling one to follow the dynamical path during the simulation. In this contribution, we present a simple proof to demonstrate that this uniform acceptance still complies with the condition of detailed balance, on the condition that the characteristic parameter lambda = 1/2 and that the maximum allowed step size is chosen to be sufficiently small. Furthermore, the relation to Metropolis Monte Carlo (MMC) is also established, and it is shown that UFMC reduces to MMC by choosing the characteristic parameter lambda = 0 [Rao, M. et al. Mol. Phys. 1979, 37, 1773]. Finally, a simple example compares the UFMC and MMC methods.
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 5.245
Times cited: 20
DOI: 10.1021/ct2008268
|
“Temperature influence on the reactivity of plasma species on a nickel catalyst surface : an atomic scale study”. Somers W, Bogaerts A, van Duin ACT, Huygh S, Bal KM, Neyts EC, Catalysis today 211, 131 (2013). http://doi.org/10.1016/j.cattod.2013.02.010
Abstract: In recent years, the potential use of hydrogen as a clean energy source has gained considerable attention. Especially H2 formation by Ni-catalyzed reforming of methane at elevated temperatures is an attractive process. However, a more fundamental knowledge at the atomic level is needed for a full comprehension of the reactions at the catalyst surface. In this contribution, we therefore investigate the H2 formation after CHx impacts on a Ni(1 1 1) surface in the temperature range 4001600 K, by means of reactive molecular dynamics (MD) simulations using the ReaxFF potential. While some H2 formation is already observed at the lower temperatures, substantial H2 formation is only obtained at elevated temperatures of 1400 K and above. At 1600 K, the H2 molecules are even the most frequently formed species. In direct correlation with the increasing dehydrogenation at elevated temperatures, an increased surface-to-subsurface C-diffusivity is observed as well. This study highlights the major importance of the temperature on the H2 formation.
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 4.636
Times cited: 27
DOI: 10.1016/j.cattod.2013.02.010
|
“On the time scale associated with Monte Carlo simulations”. Bal KM, Neyts EC, The journal of chemical physics 141, 204104 (2014). http://doi.org/10.1063/1.4902136
Abstract: Uniform-acceptance force-bias Monte Carlo (fbMC) methods have been shown to be a powerful technique to access longer timescales in atomistic simulations allowing, for example, phase transitions and growth. Recently, a new fbMC method, the time-stamped force-bias Monte Carlo (tfMC) method, was derived with inclusion of an estimated effective timescale; this timescale, however, does not seem able to explain some of the successes the method. In this contribution, we therefore explicitly quantify the effective timescale tfMC is able to access for a variety of systems, namely a simple single-particle, one-dimensional model system, the Lennard-Jones liquid, an adatom on the Cu(100) surface, a silicon crystal with point defects and a highly defected graphene sheet, in order to gain new insights into the mechanisms by which tfMC operates. It is found that considerable boosts, up to three orders of magnitude compared to molecular dynamics, can be achieved for solid state systems by lowering of the apparent activation barrier of occurring processes, while not requiring any system-specific input or modifications of the method. We furthermore address the pitfalls of using the method as a replacement or complement of molecular dynamics simulations, its ability to explicitly describe correct dynamics and reaction mechanisms, and the association of timescales to MC simulations in general.
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 2.965
Times cited: 26
DOI: 10.1063/1.4902136
|
“Merging Metadynamics into Hyperdynamics: Accelerated Molecular Simulations Reaching Time Scales from Microseconds to Seconds”. Bal KM, Neyts EC, Journal of chemical theory and computation 11, 4545 (2015). http://doi.org/10.1021/acs.jctc.5b00597
Abstract: The hyperdynamics method is a powerful tool to simulate slow processes at the atomic level. However, the construction of an optimal hyperdynamics potential is a task that is far from trivial. Here, we propose a generally applicable implementation of the hyperdynamics algorithm, borrowing two concepts from metadynamics. First, the use of a collective variable (CV) to represent the accelerated dynamics gives the method a very large flexibility and simplicity. Second, a metadynamics procedure can be used to construct a suitable history-dependent bias potential on-the-fly, effectively turning the algorithm into a self-learning accelerated molecular dynamics method. This collective variable-driven hyperdynamics (CVHD) method has a modular design: both the local system properties on which the bias is based, as well as the characteristics of the biasing method itself, can be chosen to match the needs of the considered system. As a result, system-specific details are abstracted from the biasing algorithm itself, making it extremely versatile and transparent. The method is tested on three model systems: diffusion on the Cu(001) surface and nickel-catalyzed methane decomposition, as examples of reactive processes with a bond-length-based CV, and the folding of a long polymer-like chain, using a set of dihedral angles as a CV. Boost factors up to 109, corresponding to a time scale of seconds, could be obtained while still accurately reproducing correct dynamics.
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 5.245
Times cited: 41
DOI: 10.1021/acs.jctc.5b00597
|
“Direct observation of realistic-temperature fuel combustion mechanisms in atomistic simulations”. Bal KM, Neyts EC, Chemical science 7, 5280 (2016). http://doi.org/10.1039/C6SC00498A
Abstract: Atomistic simulations can in principle provide an unbiased description of all mechanisms, intermediates, and products of complex chemical processes. However, due to the severe time scale limitation of conventional simulation techniques, unrealistically high simulation temperatures are usually applied, which are a poor approximation of most practically relevant low-temperature applications. In this work, we demonstrate the direct observation at the atomic scale of the pyrolysis and oxidation of n-dodecane at temperatures as low as 700 K through the use of a novel simulation technique, collective variable-driven hyperdynamics (CVHD). A simulated timescale of up to 39 seconds is reached. Product compositions and dominant mechanisms are found to be strongly temperature-dependent, and are consistent with experiments and kinetic models. These simulations provide a first atomic-level look at the full dynamics of the complicated fuel combustion process at industrially relevant temperatures and time scales, unattainable by conventional molecular dynamics simulations.
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 8.668
Times cited: 22
DOI: 10.1039/C6SC00498A
|
“Structures and spectroscopic properties of sulfur-nitrogen-pnictogen chains : R2P-N=S=N-PR2 and R2P-N=S=N-AsR2”. Bal KM, Cautereels J, Blockhuys F, Journal of molecular structure 1132, 102 (2017). http://doi.org/10.1016/J.MOLSTRUC.2016.08.008
Abstract: The conformational and configurational preferences of Me2PNSNPMe2 (3) and Me2PNSNAsMe2 (4) have been identified using quantum chemical calculations at the DFT/B3LYP/6-311+G* level of theory. An approach in which energetic, structural (geometries and bond orders), electronic (analysis of the electron density) and spectroscopic properties are combined leads to the conclusion that these sulfur-nitrogen-pnictogen chains share many of the properties of their chalcogen-nitrogen analogues but that the through-space intramolecular interactions favouring the Z,Z configuration are even weaker than in these latter compounds. The results of this analysis also lead to an unambiguous assignment of the variable-temperature 31P and 15N NMR spectra of these compounds and their structures both in solution and in the solid state.
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 1.753
DOI: 10.1016/J.MOLSTRUC.2016.08.008
|
“Effect of electric fields on plasma catalytic hydrocarbon oxidation from atomistic simulations”. Neyts EC, Bal KM, Plasma processes and polymers 14, e1600158 (2017). http://doi.org/10.1002/PPAP.201600158
Abstract: The catalytic oxidative dehydrogenation of hydrocarbons is an industrially important process, in which selectivity is a key issue. We here investigate the conversion of methanol to formaldehyde on a vanadia surface employing long timescale simulations, reaching a time scale of seconds. In particular, we compare the thermal process to the case where an additional external electric field is applied, as would be the case in a direct plasma-catalysis setup. We find that the electric field influences the retention time of the molecules at the catalyst surface. These simulations provide an atomic scale insight in the thermal catalytic oxidative dehydrogenation process, and in how an external electric field may affect this process.
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 2.846
Times cited: 2
DOI: 10.1002/PPAP.201600158
|
“Atomistic simulations of graphite etching at realistic time scales”. Aussems DUB, Bal K M, Morgan TW, van de Sanden MCM, Neyts EC, Chemical science 8, 7160 (2017). http://doi.org/10.1039/C7SC02763J
Abstract: Hydrogen–graphite interactions are relevant to a wide variety of applications, ranging from astrophysics to fusion devices and nano-electronics. In order to shed light on these interactions, atomistic simulation using Molecular Dynamics (MD) has been shown to be an invaluable tool. It suffers, however, from severe timescale
limitations. In this work we apply the recently developed Collective Variable-Driven Hyperdynamics (CVHD) method to hydrogen etching of graphite for varying inter-impact times up to a realistic value of 1 ms, which corresponds to a flux of 1020 m2 s1. The results show that the erosion yield, hydrogen surface coverage and species distribution are significantly affected by the time between impacts. This can be explained by the higher probability of C–C bond breaking due to the prolonged exposure to thermal stress and the subsequent transition from ion- to thermal-induced etching. This latter regime of thermal-induced etching – chemical erosion – is here accessed for the first time using atomistic simulations. In conclusion, this study demonstrates that accounting for long time-scales significantly affects ion bombardment simulations and should not be neglected in a wide range of conditions, in contrast to what is typically assumed.
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 8.668
Times cited: 3
DOI: 10.1039/C7SC02763J
|
“Effect of plasma-induced surface charging on catalytic processes: application to CO2activation”. Bal KM, Huygh S, Bogaerts A, Neyts EC, Plasma sources science and technology 27, 024001 (2018). http://doi.org/10.1088/1361-6595/aaa868
Abstract: Understanding the nature and effect of the multitude of plasma–surface interactions in plasma catalysis is a crucial requirement for further process development and improvement. A particularly intriguing and rather unique property of a plasma-catalytic setup is the ability of the plasma to modify the electronic structure, and hence chemical properties, of the catalyst through charging, i.e. the absorption of excess electrons. In this work, we develop a quantum chemical model based on density functional theory to study excess negative surface charges in a heterogeneous catalyst exposed to a plasma. This method is specifically applied to investigate plasma-catalytic CO2 activation on supported M/Al2O3 (M=Ti, Ni, Cu) single atom catalysts. We find that (1) the presence of a negative surface charge dramatically improves the reductive power of the catalyst, strongly promoting the splitting of CO2 to CO and oxygen, and (2) the relative activity of the investigated transition metals is also changed upon charging, suggesting that controlled surface charging is a powerful additional parameter to tune catalyst activity and selectivity. These results strongly point to plasma-induced surface charging of the catalyst as an important factor contributing to the plasma-catalyst synergistic effects frequently reported for plasma catalysis.
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 3.302
Times cited: 19
DOI: 10.1088/1361-6595/aaa868
|
“Modelling molecular adsorption on charged or polarized surfaces: a critical flaw in common approaches”. Bal KM, Neyts EC, Physical chemistry, chemical physics 20, 8456 (2018). http://doi.org/10.1039/C7CP08209F
Abstract: A number of recent computational material design studies based on density functional theory (DFT) calculations have put forward a new class of materials with electrically switchable chemical characteristics that can be exploited in the development of tunable gas storage and electrocatalytic applications. We find systematic flaws in almost every computational study of gas adsorption on polarized or charged surfaces, stemming from an improper and unreproducible treatment of periodicity, leading to very large errors of up to 3 eV in some cases. Two simple corrective procedures that lead to consistent results are proposed, constituting a crucial course correction to the research in the field.
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 4.123
Times cited: 8
DOI: 10.1039/C7CP08209F
|
“High Coke Resistance of a TiO2Anatase (001) Catalyst Surface during Dry Reforming of Methane”. Huygh S, Bogaerts A, Bal KM, Neyts EC, Journal Of Physical Chemistry C 122, 9389 (2018). http://doi.org/10.1021/acs.jpcc.7b10963
Abstract: The resistance of a TiO2 anatase (001) surface to coke formation was studied in the context of dry reforming of methane using density functional theory (DFT) calculations. As carbon atoms act as precursors for coke formation, the resistance to coke formation can be measured by the carbon coverage of the surface. This is related to the stability of different CHx (x = 0−3) species and their rate of hydrogenation and dehydrogenation on the TiO2 surface. Therefore, we studied the reaction mechanisms and their corresponding rates as a function of the temperature for the dehydrogenation of the species on the surface. We found that the stabilities of C and CH are significantly lower than those of CH3 and CH2. The hydrogenation rates of the different species are significantly higher than the dehydrogenation rates in a temperature range of 300−1000 K. Furthermore, we found that dehydrogenation of CH3, CH2, and CH will only occur at appreciable rates starting from 600, 900, and 900 K, respectively. On the basis of these results, it is clear that the anatase (001) surface has a high coke resistance, and it is thus not likely that the surface will become poisoned by coke during dry reforming of methane. As the rate limiting step in dry reforming is the dissociative adsorption of CH4, we studied an alternative approach to thermal catalysis. We found that the temperature threshold for dry reforming is at least 700 K. This threshold temperature may be lowered by the use of plasma-catalysis, where the appreciable rates of adsorption of plasma-generated CHx radicals result in bypassing the rate limiting step of the reaction.
Keywords: A1 Journal Article; Plasma, laser ablation and surface modeling Antwerp (PLASMANT) ;
Impact Factor: 4.536
Times cited: 1
DOI: 10.1021/acs.jpcc.7b10963
|
“Mechanisms of elementary hydrogen ion-surface interactions during multilayer graphene etching at high surface temperature as a function of flux”. Aussems DUB, Bal KM, Morgan TW, van de Sanden MCM, Neyts EC, Carbon 137, 527 (2018). http://doi.org/10.1016/j.carbon.2018.05.051
Abstract: In order to optimize the plasma-synthesis and modification process of carbon nanomaterials for applications such as nanoelectronics and energy storage, a deeper understanding of fundamental hydrogengraphite/graphene interactions is required. Atomistic simulations by Molecular Dynamics have proven to be indispensable to illuminate these phenomena. However, severe time-scale limitations restrict them to very fast processes such as reflection, while slow thermal processes such as surface diffusion and molecular desorption are commonly inaccessible. In this work, we could however reach these thermal processes for the first time at time-scales and surface temperatures (1000 K) similar to high-flux plasma exposure experiments during the simulation of multilayer graphene etching by 5 eV H ions. This was achieved by applying the Collective Variable-Driven Hyperdynamics biasing technique, which extended the inter-impact time over a range of six orders of magnitude, down to a more realistic ion-flux of 1023m2s1. The results show that this not only causes a strong shift from predominant ion-to thermally induced interactions, but also significantly affects the hydrogen uptake and surface evolution. This study thus elucidates H ion-graphite/graphene interaction mechanisms and stresses the importance of including long time-scales in atomistic simulations at high surface temperatures to understand the dynamics of the ion-surface system.
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 6.337
Times cited: 4
DOI: 10.1016/j.carbon.2018.05.051
|
“CO2 activation on TiO2-supported Cu5 and Ni5 nanoclusters : effect of plasma-induced surface charging”. Jafarzadeh A, Bal KM, Bogaerts A, Neyts EC, The journal of physical chemistry: C : nanomaterials and interfaces 123, 6516 (2019). http://doi.org/10.1021/ACS.JPCC.8B11816
Abstract: Surface charging is an often overlooked factor in many plasma-surface interactions and in particular in plasma catalysis. In this study, we investigate the effect of excess electrons induced by a plasma on the adsorption properties of CO2 on titania-supported Cu-5 and Ni-5 clusters using spin-polarized and dispersion-corrected density functional theory calculations. The effect of excess electrons on the adsorption of Ni and Cu pentamers as well as on CO2 adsorption on a pristine anatase TiO2(101) slab is studied. Our results indicate that adding plasma-induced excess electrons to the system leads to further stabilization of the bent CO2 structure. Also, dissociation of CO2 on charged clusters is energetically more favorable than on neutral clusters. We hypothesize that surface charge is a plausible cause for the synergistic effects sometimes observed in plasma catalysis.
Keywords: A1 Journal article; Engineering sciences. Technology; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 4.536
Times cited: 4
DOI: 10.1021/ACS.JPCC.8B11816
|
“Overcoming Old Scaling Relations and Establishing New Correlations in Catalytic Surface Chemistry: Combined Effect of Charging and Doping”. Bal KM, Neyts EC, The journal of physical chemistry: C : nanomaterials and interfaces 123, 6141 (2019). http://doi.org/10.1021/acs.jpcc.9b01216
Abstract: Optimization of catalytic materials for a given application is greatly constrained by linear scaling relations. Recently, however, it has been demonstrated that it is possible to reversibly modulate the chemisorption of molecules on nanomaterials by charging (i.e., injection or removal of electrons) and hence reversibly and selectively modify catalytic activity beyond structure−activity correlations. The fundamental physical relation between the properties of the material, the charging process, and the chemisorption energy, however, remains unclear, and a systematic exploration and optimization of charge-switchable sorbent materials is not yet possible. Using hybrid DFT calculations of CO2 chemisorption on hexagonal boron nitride nanosheets with several types of defects and dopants, we here reveal the existence of fundamental correlations between the electron affinity of a material and charge-induced chemisorption, show how defect engineering can be used to modulate the strength and efficiency of the adsorption process, and demonstrate that excess electrons stabilize many topological defects. We then show how these insights could be exploited in the development of new electrocatalytic materials and the synthesis of doped nanomaterials. Moreover, we demonstrate that calculated chemical properties of charged materials are highly sensitive to the employed computational methodology because of the self-interaction error, which underlines the theoretical challenge posed by such systems.
Keywords: A1 Journal article; Engineering sciences. Technology; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 4.536
Times cited: 5
DOI: 10.1021/acs.jpcc.9b01216
|
“Ensemble-Based Molecular Simulation of Chemical Reactions under Vibrational Nonequilibrium”. Bal KM, Bogaerts A, Neyts EC, Journal Of Physical Chemistry Letters 11, 401 (2020). http://doi.org/10.1021/acs.jpclett.9b03356
Abstract: We present an approach to incorporate the effect of vibrational nonequilibrium in molecular dynamics (MD) simulations. A perturbed canonical ensemble, in which selected modes are excited to higher temperature while all others remain equilibrated at low temperature, is simulated by applying a specifically tailored bias potential. Our method can be readily applied to any (classical or quantum mechanical) MD setup at virtually no additional computational cost and allows the study of reactions of vibrationally excited molecules in nonequilibrium environments such as plasmas. In combination with enhanced sampling methods, the vibrational efficacy and mode selectivity of vibrationally stimulated reactions can then be quantified in terms of chemically relevant observables, such as reaction rates and apparent free energy barriers. We first validate our method for the prototypical hydrogen exchange reaction and then show how it can capture the effect of vibrational excitation on a symmetric SN2 reaction and radical addition on CO2.
Keywords: A1 Journal article; Engineering sciences. Technology; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 5.7
DOI: 10.1021/acs.jpclett.9b03356
|
“Accelerated molecular dynamics simulation of large systems with parallel collective variable-driven hyperdynamics”. Fukuhara S, Bal KM, Neyts EC, Shibuta Y, Computational Materials Science 177, 109581 (2020). http://doi.org/10.1016/j.commatsci.2020.109581
Abstract: The limitation in time and length scale is a major issue of molecular dynamics (MD) simulation. Although several methods have been developed to extend the MD time scale, their performance usually deteriorates with increasing system size. Therefore, an acceleration method which is applicable to large systems is required to bridge the gap between the MD simulations and target phenomena. In this study, an accelerated MD method for large system is developed based on the collective variable-driven hyperdynamics (CVHD) method [K.M. Bal and E.C. Neyts, 2015]. The key idea is to run CVHD in parallel with rate control and accelerate multiple possible events simultaneously. Using this novel method, carbon diffusion in bcc-iron bicrystal with grain boundary is examined as an application for practical materials. Carbon atoms reaching at the grain boundary are trapped whereas carbon atoms in the bulk region diffuse randomly, and both dynamic regimes can be simultaneously accelerated with the parallel CVHD technique.
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 3.3
DOI: 10.1016/j.commatsci.2020.109581
|
“Activation of CO2on Copper Surfaces: The Synergy between Electric Field, Surface Morphology, and Excess Electrons”. Jafarzadeh A, Bal KM, Bogaerts A, Neyts EC, Journal Of Physical Chemistry C 124, 6747 (2020). http://doi.org/10.1021/acs.jpcc.0c00778
Abstract: In this work, we use density functional theory calculations to study the combined effect of external electric fields, surface morphology, and surface charge on CO2 activation over Cu(111), Cu(211), Cu(110), and Cu(001) surfaces. We observe that the binding energy of the CO2 molecule on Cu surfaces increases significantly upon increasing the applied electric field strength. In addition, rougher surfaces respond more effectively to the presence of the external electric field toward facilitating the formation of a carbonate-like CO2 structure and the transformation of the most stable adsorption mode from physisorption to chemisorption. The presence of surface charges further strengthens the electric field effect and consequently causes an improved bending of the CO2 molecule and C−O bond length elongation. On the other hand, a net charge in the absence of an externally applied electric field shows only a marginal effect on CO2 binding. The chemisorbed CO2 is more stable and further activated when the effects of an external electric field, rough surface, and surface charge are combined. These results can help to elucidate the underlying factors that control CO2 activation in heterogeneous and plasma catalysis, as well as in electrochemical processes.
Keywords: A1 Journal article; Engineering sciences. Technology; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 3.7
DOI: 10.1021/acs.jpcc.0c00778
|
“Free energy barriers from biased molecular dynamics simulations”. Bal KM, Fukuhara S, Shibuta Y, Neyts EC, Journal Of Chemical Physics 153, 114118 (2020). http://doi.org/10.1063/5.0020240
Abstract: Atomistic simulation methods for the quantification of free energies are in wide use. These methods operate by sampling the probability density of a system along a small set of suitable collective variables (CVs), which is, in turn, expressed in the form of a free energy surface (FES). This definition of the FES can capture the relative stability of metastable states but not that of the transition state because the barrier height is not invariant to the choice of CVs. Free energy barriers therefore cannot be consistently computed from the FES. Here, we present a simple approach to calculate the gauge correction necessary to eliminate this inconsistency. Using our procedure, the standard FES as well as its gauge-corrected counterpart can be obtained by reweighing the same simulated trajectory at little additional cost. We apply the method to a number of systems—a particle solvated in a Lennard-Jones fluid, a Diels–Alder reaction, and crystallization of liquid sodium—to demonstrate its ability to produce consistent free energy barriers that correctly capture the kinetics of chemical or physical transformations, and discuss the additional demands it puts on the chosen CVs. Because the FES can be converged at relatively short (sub-ns) time scales, a free energy-based description of reaction kinetics is a particularly attractive option to study chemical processes at more expensive quantum mechanical levels of theory.
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 4.4
DOI: 10.1063/5.0020240
|
“Entropic and enthalpic factors determining the thermodynamics and kinetics of carbon segregation from transition metal nanoparticles”. Fukuhara S, Bal KM, Neyts EC, Shibuta Y, Carbon 171, 806 (2021). http://doi.org/10.1016/j.carbon.2020.09.059
Abstract: The free energy surface (FES) for carbon segregation from nickel nanoparticles is obtained from advanced molecular dynamics simulations. A suitable reaction coordinate is developed that can distinguish dissolved carbon atoms from segregated dimers, chains and junctions on the nanoparticle surface. Because of the typically long segregation time scale (up to ms), metadynamics simulations along the developed reaction coordinate are used to construct FES over a wide range of temperatures and carbon concentrations. The FES revealed the relative stability of different stages in the segregation process, and free energy barriers and rates of the individual steps could then be calculated and decomposed into enthalpic and entropic contributions. As the carbon concentration in the nickel nanoparticle increases, segregated carbon becomes more stable in terms of both enthalpy and entropy. The activation free energy of the reaction also decreases with the increase of carbon concentration, which can be mainly attributed to entropic effects. These insights and the methodology developed to obtain them improve our understanding of carbon segregation process across materials science in general, and the nucleation and growth of carbon nanotube in particular.
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 6.337
DOI: 10.1016/j.carbon.2020.09.059
|
“Reaction mechanisms of C(3PJ) and C+(2PJ) with benzene in the interstellar medium from quantum mechanical molecular dynamics simulations”. Izadi ME, Bal KM, Maghari A, Neyts EC, Physical Chemistry Chemical Physics 23, 4205 (2021). http://doi.org/10.1039/D0CP04542J
Abstract: While spectroscopic data on small hydrocarbons in interstellar media in combination with crossed molecular beam (CMB) experiments have provided a wealth of information on astrochemically relevant species, much of the underlying mechanistic pathways of their formation remain elusive. Therefore, in this work, the chemical reaction mechanisms of C(<sup>3</sup>P<sub>J</sub>) + C<sub>6</sub>H<sub>6</sub>and C<sup>+</sup>(<sup>2</sup>P) + C<sub>6</sub>H<sub>6</sub>systems using the quantum mechanical molecular dynamics (QMMD) technique at the PBE0-D3(BJ) level of theory is investigated, mimicking a CMB experiment. Both the dynamics of the reactions as well as the electronic structure for the purpose of the reaction network are evaluated. The method is validated for the first reaction by comparison to the available experimental data. The reaction scheme for the C(<sup>3</sup>P<sub>J</sub>) + C<sub>6</sub>H<sub>6</sub>system covers the literature data,<italic>e.g.</italic>the major products are the 1,2-didehydrocycloheptatrienyl radical (C<sub>7</sub>H<sub>5</sub>) and benzocyclopropenyl radical (C<sub>6</sub>H<sub>5</sub>–CH), and it reveals the existence of less common pathways for the first time. The chemistry of the C<sup>+</sup>(<sup>2</sup>P<sub>J</sub>) + C<sub>6</sub>H<sub>6</sub>system is found to be much richer, and we have found that this is because of more exothermic reactions in this system in comparison to those in the C(<sup>3</sup>P<sub>J</sub>) + C<sub>6</sub>H<sub>6</sub>system. Moreover, using the QMMD simulation, a number of reaction paths have been revealed that produce three distinct classes of reaction products with different ring sizes. All in all, at all the collision energies and orientations, the major product is the heptagon molecular ion for the ionic system. It is also revealed that the collision orientation has a dominant effect on the reaction products in both systems, while the collision energy mostly affects the charged system. These simulations both prove the applicability of this approach to simulate crossed molecular beams, and provide fundamental information on reactions relevant for the interstellar medium.
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 4.123
DOI: 10.1039/D0CP04542J
|
“Reactive plasma cleaning and restoration of transition metal dichalcogenide monolayers”. Marinov D, de Marneffe J-F, Smets Q, Arutchelvan G, Bal KM, Voronina E, Rakhimova T, Mankelevich Y, El Kazzi S, Nalin Mehta A, Wyndaele P-J, Heyne MH, Zhang J, With PC, Banerjee S, Neyts EC, Asselberghs I, Lin D, De Gendt S, npj 2D Materials and Applications 5, 17 (2021). http://doi.org/10.1038/s41699-020-00197-7
Abstract: The cleaning of two-dimensional (2D) materials is an essential step in the fabrication of future devices, leveraging their unique physical, optical, and chemical properties. Part of these emerging 2D materials are transition metal dichalcogenides (TMDs). So far there is limited understanding of the cleaning of “monolayer” TMD materials. In this study, we report on the use of downstream H<sub>2</sub>plasma to clean the surface of monolayer WS<sub>2</sub>grown by MOCVD. We demonstrate that high-temperature processing is essential, allowing to maximize the removal rate of polymers and to mitigate damage caused to the WS<sub>2</sub>in the form of sulfur vacancies. We show that low temperature in situ carbonyl sulfide (OCS) soak is an efficient way to resulfurize the material, besides high-temperature H<sub>2</sub>S annealing. The cleaning processes and mechanisms elucidated in this work are tested on back-gated field-effect transistors, confirming that transport properties of WS<sub>2</sub>devices can be maintained by the combination of H<sub>2</sub>plasma cleaning and OCS restoration. The low-damage plasma cleaning based on H<sub>2</sub>and OCS is very reproducible, fast (completed in a few minutes) and uses a 300 mm industrial plasma etch system qualified for standard semiconductor pilot production. This process is, therefore, expected to enable the industrial scale-up of 2D-based devices, co-integrated with silicon technology.
Keywords: A1 Journal article; Engineering sciences. Technology; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
DOI: 10.1038/s41699-020-00197-7
|
Bal KM, Neyts EC (2021) Quantifying the impact of vibrational nonequilibrium in plasma catalysis: insights from a molecular dynamics model of dissociative chemisorption. 394004
Abstract: The rate, selectivity and efficiency of plasma-based conversion processes is strongly affected by nonequilibrium phenomena. High concentrations of vibrationally excited molecules are such a plasma-induced effect. It is frequently assumed that vibrationally excited molecules are important in plasma catalysis because their presence lowers the apparent activation energy of dissociative chemisorption reactions and thus increases the conversion rate. A detailed atomic-level understanding of vibrationally stimulated catalytic reactions in the context of plasma catalysis is however lacking. Here, we couple a recently developed statistical model of a plasma-induced vibrational nonequilibrium to molecular dynamics simulations, enhanced sampling methods, and machine learning techniques. We quantify the impact of a vibrational nonequilibrium on the dissociative chemisorption barrier of H2 and CH4 on nickel catalysts over a wide range of vibrational temperatures. We investigate the effect of surface structure and compare the role of different vibrational modes of methane in the dissociation process. For low vibrational temperatures, very high vibrational efficacies are found, and energy in bend vibrations appears to dominate the dissociation of methane. The relative impact of vibrational nonequilibrium is much higher on terrace sites than on surface steps. We then show how our simulations can help to interpret recent experimental results, and suggest new paths to a better understanding of plasma catalysis.
Keywords: A1 Journal Article;plasma catalysis; vibrational nonequilibrium; dissociative chemisorption; free energy barriers; Plasma, laser ablation and surface modeling Antwerp (PLASMANT) ;
Impact Factor: 2.588
DOI: 10.1088/1361-6463/ac113a
|
“Nucleation rates from small scale atomistic simulations and transition state theory”. Bal KM, Journal Of Chemical Physics 155, 144111 (2021). http://doi.org/10.1063/5.0063398
Abstract: The evaluation of nucleation rates from molecular dynamics trajectories is hampered by the slow nucleation time scale and impact of finite size effects. Here, we show that accurate nucleation rates can be obtained in a very general fashion relying only on the free energy barrier, transition state theory, and a simple dynamical correction for diffusive recrossing. In this setup, the time scale problem is overcome by using enhanced sampling methods, in casu metadynamics, whereas the impact of finite size effects can be naturally circumvented by reconstructing the free energy surface from an appropriate ensemble. Approximations from classical nucleation theory are avoided. We demonstrate the accuracy of the approach by calculating macroscopic rates of droplet nucleation from argon vapor, spanning 16 orders of magnitude and in excellent agreement with literature results, all from simulations of very small (512 atom) systems.
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 2.965
DOI: 10.1063/5.0063398
|
“Reweighted Jarzynski sampling : acceleration of rare events and free energy calculation with a bias potential learned from nonequilibrium work”. Bal KM, Journal Of Chemical Theory And Computation 17, 6766 (2021). http://doi.org/10.1021/ACS.JCTC.1C00574
Abstract: We introduce a simple enhanced sampling approach for the calculation of free energy differences and barriers along a one-dimensional reaction coordinate. First, a small number of short nonequilibrium simulations are carried out along the reaction coordinate, and the Jarzynski equality is used to learn an approximate free energy surface from the nonequilibrium work distribution. This free energy estimate is represented in a compact form as an artificial neural network and used as an external bias potential to accelerate rare events in a subsequent molecular dynamics simulation. The final free energy estimate is then obtained by reweighting the equilibrium probability distribution of the reaction coordinate sampled under the influence of the external bias. We apply our reweighted Jarzynski sampling recipe to four processes of varying scales and complexities.spanning chemical reaction in the gas phase, pair association in solution, and droplet nucleation in supersaturated vapor. In all cases, we find reweighted Jarzynski sampling to be a very efficient strategy, resulting in rapid convergence of the free energy to high precision.
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 5.245
DOI: 10.1021/ACS.JCTC.1C00574
|
“Extending and validating bubble nucleation rate predictions in a Lennard-Jones fluid with enhanced sampling methods and transition state theory”. Bal KM, Neyts EC, Journal Of Chemical Physics 157, 184113 (2022). http://doi.org/10.1063/5.0120136
Abstract: We calculate bubble nucleation rates in a Lennard-Jones fluid through explicit molecular dynamics simulations. Our approach-based on a recent free energy method (dubbed reweighted Jarzynski sampling), transition state theory, and a simple recrossing correction-allows us to probe a fairly wide range of rates in several superheated and cavitation regimes in a consistent manner. Rate predictions from this approach bridge disparate independent literature studies on the same model system. As such, we find that rate predictions based on classical nucleation theory, direct brute force molecular dynamics simulations, and seeding are consistent with our approach and one another. Published rates derived from forward flux sampling simulations are, however, found to be outliers. This study serves two purposes: First, we validate the reliability of common modeling techniques and extrapolation approaches on a paradigmatic problem in materials science and chemical physics. Second, we further test our highly generic recipe for rate calculations, and establish its applicability to nucleation processes.
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 4.4
DOI: 10.1063/5.0120136
|