Towards the identification of the gold binding region within trypsin stabilized nanoclusters using microwave synthesis routes.

Elucidating the location of stabilized nanoclusters within their protein hosts is an existing challenge towards the optimized development of functional protein-nanoclusters. While nanoclusters of various metal compositions can be readily synthesized within a wide array of protein hosts and exhibit tailorable properties, the inability to identify the cluster stabilization region prevents controllable property manipulation of both metallic and protein components. Additionally, the ability to synthesize protein-nanoclusters in a consistent and high-throughput fashion is also highly desirable. In this effort, trypsin stabilized gold nanoclusters are synthesized through standard and microwave-enabled methodologies to determine the impact of processing parameters on the materials physical and functional properties. Density functional theory simulations are employed to localize high probability regions within the trypsin enzyme for Au25 cluster stabilization, which reveal that cluster location is likely within close proximity of the trypsin active region. Trypsin activity measurements support our findings from DFT, as trypsin enzymatic activity is eliminated following cluster growth and stabilization. Moreover, studies on the reactivity of Au NCs and synchrotron characterization measurements further reveal that clusters made by microwave-based techniques exhibit slight structural differences to those made via standard methodologies, indicating that microwave-based syntheses largely maintain the native structural attributes despite the faster synthetic conditions. Overall, this work illustrates the importance of understanding the connections between synthetic conditions, atomic-scale structure, and materials properties that can be potentially used to further tune the properties of metal cluster-protein materials for future applications.

Dissipative particle dynamics with reactions: Application to RDX decomposition.

We present a general, flexible framework for a constant-energy variant of the dissipative particle dynamics method that allows chemical reactions (DPD-RX). In our DPD-RX approach, reaction progress variables are assigned to each particle that monitor the time evolution of an extent-of-reaction associated with the prescribed reaction mechanisms and kinetics assumed to occur within the particle, where chemistry can be modeled using complex or reduced reaction mechanisms. We demonstrate our DPD-RX method by considering thermally initiated unimolecular decomposition of the energetic material, cyclotrimethylene trinitramine (RDX), into a molecular gas mixture. Studies are performed to demonstrate the effect of a spatially averaged particle internal temperature and a local reaction volume term in the chemical kinetics expressions, where both provide implicit mechanisms for capturing condensed phase reactivity. We also present an analysis of the expansion of the product gas mixture during decomposition. Finally, a discussion of other potential applications and extensions of the DPD-RX method is given.

Freestanding organogels by molecular velcro of unsaturated amphiphiles.

A simple amphiphile, N-cardanyltaurine amide (NCT) with different degrees of cis-unsaturation in its tail resulted in the formation of strong organogels. Interestingly, this is in contrast to the commonly accepted notion that introducing unsaturation in alkyl chains enhances fluidity in lipid assemblies. The physico-chemical and first-principles DFT calculations confirmed the pegging of 'kinked' unsaturated side chains, where the hydrophobic interlocking as in Velcro fasteners leads to a network of cylindrical micelles, resulting in self-standing organogels. Textural profile analysis and spectroscopic details substantiated the dynamic assembly to resemble a 3D network of gelators rather than being a cross-linked or polymerized matrix of monomers.

A coarse-grain force field for RDX: Density dependent and energy conserving.

We describe the development of a density-dependent transferable coarse-grain model of crystalline hexahydro-1,3,5-trinitro-s-triazine (RDX) that can be used with the energy conserving dissipative particle dynamics method. The model is an extension of a recently reported one-site model of RDX that was developed by using a force-matching method. The density-dependent forces in that original model are provided through an interpolation scheme that poorly conserves energy. The development of the new model presented in this work first involved a multi-objective procedure to improve the structural and thermodynamic properties of the previous model, followed by the inclusion of the density dependency via a conservative form of the force field that conserves energy. The new model accurately predicts the density, structure, pressure-volume isotherm, bulk modulus, and elastic constants of the RDX crystal at ambient pressure and exhibits transferability to a liquid phase at melt conditions.

Free-energy calculations using classical molecular simulation: application to the determination of the melting point and chemical potential of a flexible RDX model.

We present an extension of various free-energy methodologies to determine the chemical potential of the solid and liquid phases of a fully-flexible molecule using classical simulation. The methods are applied to the Smith-Bharadwaj atomistic potential representation of cyclotrimethylene trinitramine (RDX), a well-studied energetic material, to accurately determine the solid and liquid phase Gibbs free energies, and the melting point (Tm). We outline an efficient technique to find the absolute chemical potential and melting point of a fully-flexible molecule using one set of simulations to compute the solid absolute chemical potential and one set of simulations to compute the solid-liquid free energy difference. With this combination, only a handful of simulations are needed, whereby the absolute quantities of the chemical potentials are obtained, for use in other property calculations, such as the characterization of crystal polymorphs or the determination of the entropy. Using the LAMMPS molecular simulator, the Frenkel and Ladd and pseudo-supercritical path techniques are adapted to generate 3rd order fits of the solid and liquid chemical potentials. Results yield the thermodynamic melting point Tm = 488.75 K at 1.0 atm. We also validate these calculations and compare this melting point to one obtained from a typical superheated simulation technique.

Lattice strain due to an atomic vacancy.

Volumetric strain can be divided into two parts: strain due to bond distance change and strain due to vacancy sources and sinks. In this paper, efforts are focused on studying the atomic lattice strain due to a vacancy in an FCC metal lattice with molecular dynamics simulation (MDS). The result has been compared with that from a continuum mechanics method. It is shown that using a continuum mechanics approach yields constitutive results similar to the ones obtained based purely on molecular dynamics considerations.