Interfacing CRYSTAL/AMBER to Optimize QM/MM Lennard⁻Jones Parameters for Water and to Study Solvation of TiO₂ Nanoparticles

A review on microfluidic-assisted nanoparticle synthesis, and their applications using multiscale simulation methods

Recent years have witnessed an increased interest in the development of nanoparticles (NPs) owing to their potential use in a wide variety of biomedical applications, including drug delivery, imaging agents, gene therapy, and vaccines, where recently, lipid nanoparticle mRNA-based vaccines were developed to prevent SARS-CoV-2 causing COVID-19. NPs typically fall into two broad categories: organic and inorganic. Organic NPs mainly include lipid-based and polymer-based nanoparticles, such as liposomes, solid lipid nanoparticles, polymersomes, dendrimers, and polymer micelles. Gold and silver NPs, iron oxide NPs, quantum dots, and carbon and silica-based nanomaterials make up the bulk of the inorganic NPs. These NPs are prepared using a variety of top-down and bottom-up approaches. Microfluidics provide an attractive synthesis alternative and is advantageous compared to the conventional bulk methods. The microfluidic mixing-based production methods offer better control in achieving the desired size, morphology, shape, size distribution, and surface properties of the synthesized NPs. The technology also exhibits excellent process repeatability, fast handling, less sample usage, and yields greater encapsulation efficiencies. In this article, we provide a comprehensive review of the microfluidic-based passive and active mixing techniques for NP synthesis, and their latest developments. Additionally, a summary of microfluidic devices used for NP production is presented. Nonetheless, despite significant advancements in the experimental procedures, complete details of a nanoparticle-based system cannot be deduced from the experiments alone, and thus, multiscale computer simulations are utilized to perform systematic investigations. The work also details the most common multiscale simulation methods and their advancements in unveiling critical mechanisms involved in nanoparticle synthesis and the interaction of nanoparticles with other entities, especially in biomedical and therapeutic systems. Finally, an analysis is provided on the challenges in microfluidics related to nanoparticle synthesis and applications, and the future perspectives, such as large-scale NP synthesis, and hybrid formulations and devices.

Transferable Potential Function for Flexible H2O Molecules Based on the Single-Center Multipole Expansion

A potential function is presented for describing a system of flexible H₂O molecules based on the single-center multipole expansion (SCME) of the electrostatic interaction. The model, referred to as SCME/f, includes the variation of the molecular quadrupole moment as well as the dipole moment with changes in bond length and angle so as to reproduce results of high-level electronic structure calculations. The multipole expansion also includes fixed octupole and hexadecapole moments, as well as anisotropic dipole-dipole, dipole-quadrupole, and quadrupole-quadrupole polarizability tensors. The model contains five adjustable parameters related to the repulsive interaction and damping functions in the electrostatic and dispersion interactions. Their values are adjusted to reproduce the lowest energy isomers of small clusters, (H₂O)_n with n = 2-6, as well as measured properties of the ice Ih crystal. Subsequent calculations of the energy difference between the various isomer configurations of the clusters show that SCME/f gives good agreement with results of electronic structure calculations and represents a significant improvement over the previously presented rigid SCME potential function. Analysis of the vibrational frequencies of the clusters and structural properties of ice Ih crystal show the importance of accurately describing the variation of the quadrupole moment with molecular structures.

Dopamine-Decorated TiO2 Nanoparticles in Water: A QM/MM vs an MM Description

Nanoparticle functionalization is a modern strategy in nanotechnology to build up devices for several applications. Modeling fully decorated metal oxide nanoparticles of realistic size (few nanometers) in an aqueous environment is a challenging task. In this work, we present a case study relevant for solar-light exploitation and for biomedical applications, i.e., a dopamine-functionalized TiO₂ nanoparticle (1700 atoms) in bulk water, for which we have performed an extensive comparative investigation with both MM and QM/MM approaches of the structural properties and of the conformational dynamics. We have used a combined multiscale protocol for a more efficient exploration of the complex conformational space. On the basis of the results of this study and of some QM and experimental data, we have defined strengths and limitations of the existing force field parameters. Our findings will be useful for an improved modeling and simulation of many other similar hybrid bioinorganic nanosystems in an aqueous environment that are pivotal in a broad range of nanotechnological applications.

Polarizable Embedding with a Transferable H2O Potential Function II: Application to (H2O)n Clusters and Liquid Water

The incorporation of polarization in multiscale quantum-mechanics/molecular-mechanics (QM/MM) simulations is important for a variety of applications, for example, charge-transfer reactions. A recently developed formalism based on a density functional theory description of the QM region and a potential energy function for H₂O molecules that includes quadrupole as well as dipole polarizability of the MM region is used to simulate liquid water and water clusters. Analysis of the energy, atomic forces, MM polarization, and structure is presented. A quantitative assessment of the QM/MM-MM/MM interaction energy differences of all possible QM/MM configurations of (H₂O)_n clusters shows that the interquartile range of the distributions of the QM/MM binding energies is never more than 20 meV/molecule higher or lower than the binding energies produced with either of the single-model results. Comparing these interaction energy differences with the QM/MM induction differences show that they are not systematically caused by the induced MM moments of our polarizable embedding scheme. Optimized hexamer geometries as well as the liquid water structure are shown to be improved in comparison with results obtained using point-charge based embedding models neglecting polarization.

Polarizable Embedding with a Transferable H2O Potential Function I: Formulation and Tests on Dimer

The incorporation of mutual polarization in multiscale simulations where different regions of the system are treated at different level of theory is important in studies of, for example, electronic excitations and charge transfer processes. We present here an energy functional for describing a quantum mechanics/molecular mechanics (QM/MM) scheme that includes reciprocal polarization between the two subsystems. The inclusion of polarization alleviates shortcomings inherent in electrostatic embedding QM/MM models based on point-charge force fields. A density functional theory (DFT) description of the QM subsystem is coupled to a single center multipole expansion (SCME) description of H₂O molecules in the MM subsystem that includes anisotropic dipole and quadrupole polarizability as well as static multipoles up to and including the hexadecapole. The energy functional and the coupling scheme is general and can be extended to arbitrary order in terms of both the static and induced moments. Tests of the energy surface for the H₂O dimer show that the QM/MM results lie in between the pure DFT and pure SCME values. The consistency of the many-body contributions to the energy and analytical forces is demonstrated for an H₂O pentamer.