Diffusion Tensors of Arbitrary-Shaped Nanoparticles in Fluid by Molecular Dynamics Simulation

Effective action for dissipative and nonholonomic systems

We show that the action of a dynamical system can be supplemented by an effective action for its environment to reproduce arbitrary coordinate dependent ohmic dissipation and gyroscopic forces. The action is a generalization of the harmonic bath model and describes a set of massless interacting scalar fields in an auxiliary space coupled to the original system at the boundary. A certain limit of the model implements nonholonomic constraints. In the case of dynamics with nonlinearly realized symmetries the effective action takes the form of a two-dimensional nonlinear σ model. It provides a basis for application of path integral methods to general dissipative and nonholonomic systems.

Orientation of graphene nanosheets in suspension under an electric field: theoretical model and molecular dynamic simulations

Orientation regulation of nanoparticles in a suspension by an electric field is a powerful tool to tune its mechanical, thermal, optical, electrical properties etc. However, how molecular modification can affect the orientation of two-dimensional nanoparticles is still unclear. In this paper, the influence of molecular modification on the orientation of graphene nanosheets (GNS) in water was investigated through theoretical analyses and molecular dynamics (MD) simulations. Firstly, a new orientation angle model was proposed, which considers hydration effects, dipole moments and resistance torque. Then, MD simulations were conducted to investigate the effects of position, direction, type, and number of functional groups on the orientation of GNS. The trend observed in MD simulations is consistent with the proposed theoretical model. The results reveal that, under the combined influence of the dipole moment and hydration effects, the modification with hydrophilic functional groups can reduce the orientation angle from 21.31° to 8.34°, while the modification with hydrophobic functional groups increases it to 26.43°. Among the hydrophilic functional groups, orientation of hydroxylated GNS is the best. With an increase in the number of hydroxyl groups, orientation angle is decreased from 12.61° to 8.34°. This work can provide valuable guidance for the design of high-performance suspensions and composites, such as thermal smart materials with adjustable thermal conductivity and intelligent devices with tailored capabilities.

Computational modeling of passive transport of functionalized nanoparticles

Functionalized nanoparticles (NPs) are complex objects present in a variety of systems ranging from synthetic grafted nanoparticles to viruses. The morphology and number of the decorating groups can vary widely between systems. Thus, the modeling of functionalized NPs typically considers simplified spherical objects as a first-order approximation. At the nanoscale label, complex hydrodynamic interactions are expected to emerge as the morphological features of the particles change, and they can be further amplified when the NPs are confined or near walls. Direct estimation of these variations can be inferred via diffusion coefficients of the NPs. However, the evaluation of the coefficients requires an improved representation of the NPs morphology to reproduce important features hidden by simplified spherical models. Here, we characterize the passive transport of free and confined functionalized nanoparticles using the Rigid Multi-Blob (RMB) method. The main advantage of RMB is its versatility to approximate the mobility of complex structures at the nanoscale with significant accuracy and reduced computational cost. In particular, we investigate the effect of functional groups' distribution, size, and morphology over nanoparticle translational and rotational diffusion. We identify that the presence of functional groups significantly affects the rotational diffusion of the nanoparticles; moreover, the morphology of the groups and number induce characteristic mobility reduction compared to non-functionalized nanoparticles. Confined NPs also evidenced important alterations in their diffusivity, with distinctive signatures in the off-diagonal contributions of the rotational diffusion. These results can be exploited in various applications, including biomedical, polymer nanocomposite fabrication, drug delivery, and imaging.

All-Atom Simulation Method for Zeeman Alignment and Dipolar Assembly of Magnetic Nanoparticles

Magnetic nanoparticles (MNPs) can organize into novel structures in solutions with excellent order and unique geometries. However, studies of the self-assembly of smaller MNPs are challenging due to a complicated interplay between external magnetic fields and van der Waals, electrostatic, dipolar, steric, and hydrodynamic interactions. Here, we present a novel all-atom molecular dynamics simulation method to enable detailed studies of the dynamics, self-assembly, structure, and properties of MNPs as a function of core sizes and shapes, ligand chemistry, solvent properties, and external field. We demonstrate the use and effectiveness of the model by simulating the self-assembly of oleic acid ligand-functionalized magnetite (Fe₃O₄) nanoparticles, with spherical and cubic shapes, into rings, lines, chains, and clusters under a uniform external magnetic field. We found that the long-range electrostatic interactions can favor the formation of a chain over a ring, the ligands promote MNP cluster growth, and the solvent can reduce the rotational diffusion of the MNPs. The algorithm has been parallelized to take advantage of multiple processors of a modern computer and can be used as a plugin for the popular simulation software LAMMPS to study the behavior of small MNPs and gain insights into the physics and chemistry of different magnetic assembly processes with atomistic details.

Anisotropic diffusion of membrane proteins at experimental timescales

Single-particle tracking (SPT) experiments of lipids and membrane proteins provide a wealth of information about the properties of biomembranes. Careful analysis of SPT trajectories can reveal deviations from ideal Brownian behavior. Among others, this includes confinement effects and anomalous diffusion, which are manifestations of both the nanoscale structure of the underlying membrane and the structure of the diffuser. With the rapid increase in temporal and spatial resolution of experimental methods, a new aspect of the motion of the particle, namely, anisotropic diffusion, might become relevant. This aspect that so far received only little attention is the anisotropy of the diffusive motion and may soon provide an additional proxy to the structure and topology of biomembranes. Unfortunately, the theoretical framework for detecting and interpreting anisotropy effects is currently scattered and incomplete. Here, we provide a computational method to evaluate the degree of anisotropy directly from molecular dynamics simulations and also point out a way to compare the obtained results with those available from SPT experiments. In order to probe the effects of anisotropic diffusion, we performed coarse-grained molecular dynamics simulations of peripheral and integral membrane proteins in flat and curved bilayers. In agreement with the theoretical basis, our computational results indicate that anisotropy can persist up to the rotational relaxation time [τ=(2D_r)^-1], after which isotropic diffusion is observed. Moreover, the underlying topology of the membrane bilayer can couple with the geometry of the particle, thus extending the spatiotemporal domain over which this type of motion can be detected.