A new molecular dynamics (MD) approach for computer simulations of irradiation driven chemical transformations of complex molecular systems is suggested. The approach is based on the fact that irradiation induced quantum transformations can often be treated as random, fast and local processes involving small molecules or molecular fragments. We advocate that the quantum transformations, such as molecular bond breaks, creation and annihilation of dangling bonds, electronic charge redistributions, changes in molecular topologies, etc, could be incorporated locally into the molecular force fields that describe the classical MD of complex molecular systems under irradiation. The proposed irradiation driven molecular dynamics (IDMD) methodology is designed for the molecular level description of the irradiation driven chemistry. The IDMD approach is implemented into the MBN Explorer software package capable to operate with a large library of classical potentials, many-body force fields and their combinations. IDMD opens a broad range of possibilities for modelling of irradiation driven modifications and chemistry of complex molecular systems ranging from radiotherapy cancer treatments to the modern technologies such as focused electron beam deposition (FEBID). As an example, the new methodology is applied for studying the irradiation driven chemistry caused by FEBID of tungsten hexacarbonyl W(CO)6 precursor molecules on a hydroxylated SiO2 surface. It is demonstrated that knowing the interaction parameters for the fragments of the molecular system arising in the course of irradiation one can reproduce reasonably well experimental observations and make predictions about the morphology and molecular composition of nanostructures that emerge on the surface during the FEBID process.

No abstract

Focused electron beam induced deposition (FEBID) is a powerful technique for 3D-printing of complex nanodevices. However, for resolutions below 10 nm, it struggles to control size, morphology and composition of the structures, due to a lack of molecular-level understanding of the underlying irradiation-driven chemistry (IDC). Computational modeling is a tool to comprehend and further optimize FEBID-related technologies. Here we utilize a novel multiscale methodology which couples Monte Carlo simulations for radiation transport with irradiation-driven molecular dynamics for simulating IDC with atomistic resolution. Through an in depth analysis of $$\hbox {W(CO)}_6$$ W(CO) 6 deposition on $$\hbox {SiO}_2$$ SiO 2 and its subsequent irradiation with electrons, we provide a comprehensive description of the FEBID process and its intrinsic operation. Our analysis reveals that simulations deliver unprecedented results in modeling the FEBID process, demonstrating an excellent agreement with available experimental data of the simulated nanomaterial composition, microstructure and growth rate as a function of the primary beam parameters. The generality of the methodology provides a powerful tool to study versatile problems where IDC and multiscale phenomena play an essential role.

Abstract. The concept of molecular mechanics force field has been widely accepted nowadays for studying various processes in biomolecular systems. In this paper, we suggest a modification for the standard CHARMM force field that permits simulations of systems with dynamically changing molecular topologies. The implementation of the modified force field was carried out in the popular program MBN Explorer, and, to support the development, we provide several illustrative case studies where dynamical topology is necessary. In particular, it is shown that the modified molecular mechanics force field can be applied for studying processes where rupture of chemical bonds plays an essential role, e.g., in irradiation-or collision-induced damage, and also in transformation and fragmentation processes involving biomolecular systems.