Molecular dynamics simulations of heat conduction in multi-walled carbon nanotubes

Heat conduction in multi-walled carbon nanotubes (MWNTs) was studied using non-equilibrium molecular dynamics simulations. This research focuses on the effects of the multi-wall structure of the MWNTs on the heat conduction. The results show that the thermal conductivity of a MWNT is almost the same as that of the corresponding single-walled carbon nanotubes (SWNTs) rather than much smaller as has been suggested. Thus, the multi-wall structure does not significantly affect the thermal conduction in the MWNTs. Analysis of the temperature profiles and the phonon density of states confirms that there is almost no heat transport between the MWNT layers and that each layer conducts heat nearly independently along parallel channels. This is physically reasonable since the weak inter-wall interactions and large interfacial thermal resistances make the MWNT layers behave like parallel thermal circuits.     

Molecular dynamics simulations of carbon nanotube oscillators deformed by encapsulated copper nanowires

Pure carbon nanotube (CNT) oscillators are compared to the corresponding CNT oscillators encapsulating copper nanowires (Cu@CNTs) by molecular dynamics simulations. The classical oscillation theory provides a fairly good estimate of the mass dependence of the operating frequency when the CNT surface is not deformed by the Cu nanowire. The structural deformations of the CNT induced by the encapsulated copper nanowire have a greater effect on the oscillation frequency than the mass of the copper nanowire. The excess forces of the Cu@CNT oscillator are slightly higher than those of the CNT oscillator and the excess van der Waals forces induced by the inter-wall interactions are 17 times higher than the excess forces induced by the Cu nanowire–CNT interactions.     

Molecular dynamics simulation of doxorubicin adsorption on a bundle of functionalized CNT

In this study, molecular dynamics simulation is used to investigate the adsorption of an anticancer drug, doxorubicin, on bundles of functionalized single-walled carbon nanotubes (SWNTs) in an aqueous solution. Carboxylic group has been selected as the functional group. Molecular dynamics (MD) simulations are performed for both separated systems containing a SWNT bundle and a functionalized carbon nanotube bundle, and results are compared with existing experimental data. MD results show that doxorubicin can be adsorbed on CNTs using different methods such as entrapment within CNT bundle, attachment to the side wall of the CNT, and adsorption on the CNT inner cavity. For functionalized CNT, the adsorption of drugs on the functional groups is essential for predicting the enhancement of drug loading on the functionalized nanotubes. Furthermore, the adsorption behavior of doxorubicin on CNTs is fitted with Langmuir and Freundlich isotherm models. The results show that Langmuir model can predict the adsorption behavior of doxorubicin on CNTs more accurately than Freundlich model does. As predicted by this isotherm model, the adsorption process of doxorubicin on CNTs is relatively difficult, but it can be improved by increasing the functional groups on the CNTs surface.     

Interfacial adhesion properties of graphene sheet on nanoscale corrugated surface: a molecular dynamics simulation study

Adhesive contacts between graphene sheet (GS) and corrugated substrates made of an ordered array of atomic pillars with variable geometries were investigated by molecular dynamics simulations. Depending on the height and interval distance of the pillars, GS can conformably coat the surface, partially adhere, or remain flat on top of the pillars. The relationship between the geometries of the pillar and the final adhesion configurations of GS was partially established. A critical adsorption energy was determined to achieve stable adsorption configuration of GS on corrugated substrates made of ordered pillar arrays. Besides the geometries of pillars, the effects of initial coating angle of GS were also considered as an important factor that affects the final adsorption configuration. We observed two interesting morphologies of GS, ‘I shape’ and ‘L shape’, which were determined by the initial coating angles.     

Molecular dynamics simulations of the aggregation behaviour of overlapped graphene sheets in linear aliphatic hydrocarbons

Molecular dynamics simulations were used to investigate the aggregation of two partially overlapped graphene sheets in hexane, dodecane and eicosane. When partially overlapped graphene sheets are adjacent to one another, they will expel the adsorbed layers of the solvent molecules on the graphene surface, and the amount of overlap will increase. When the overlapped regions of the graphene sheets are separated by solvent molecules, they cannot expel the adsorption layers between them, and so the sheets remain separated. The driving force for aggregation is the van der Waals interaction between the two graphene sheets, while the van der Waals interaction between the graphene sheets and the solvent molecules inhibits graphene aggregation. The diffusion rate of the hydrocarbon molecules with shorter chain lengths is higher. Thus, they diffuse faster during graphene aggregation, which leads to a higher rate of graphene overlapping in the shorter hydrocarbons. This work provides useful insights into graphene aggregation in linear hydrocarbon solvents of varying lengths at the nanoscale.     

Molecular dynamics simulation of lactate dehydrogenase adsorption onto pristine and carboxylic-functionalized graphene

ABSTRACT

Lactate dehydrogenase (LDH) is a tetrameric enzyme which is composed of two subunits known as LDHA and LDHB, which are encoded by the LDHA and LDHB genes respectively. LDH catalyses the last step in anaerobic glycolysis through the reversible conversion of pyruvate to lactate via coupled oxidation of NADH cofactor. The LDHA plays an important regulatory role in anaerobic glycolysis, by catalysing the final step of the process. Therefore, it is likely that increases in the expression level of LDHA in cancer cells could facilitate the efficiency of anaerobic glycolysis. Measuring the level of serum LDHA is a key step in the diagnosis of many cancer types. In this study, the adsorption, stability, and dynamics of LDHA on the surface of pristine graphene (PG) and carboxylated graphene (COOH-Graphene) were investigated using its molecular dynamics simulation. Variations in root mean square deviation, root mean square fluctuation, solvent accessible surface area and adsorption energy of the LDHA during the simulation were calculated to analyse the effect of PG and COOH-Graphene on the overall conformation of LDHA. Results showed that the adsorption of LDHA on COOH-Graphene is mostly mediated by electrostatic interactions, whereas on the PG, both Van der Waals and π-π interactions are prominent.     

Molecular dynamics simulation of dispersion improvement of graphene sheets in nanofluids by steric hindrance resulting from functional groups

Nanofluids are candidate materials for thermal management of heat transfer equipment. Practical applications of thermally enhanced nanofluids contribute to the reduction of weight of systems, leading to improved energy efficiency. Microsize particles sink into the systems because of gravity, therefore rendering the addition meaningless in terms of improving thermal properties. However, nanoparticles can be buoyant, leading to Brownian motion in the fluid, when they do not aggregate with each other. The most important factor in nanofluids is long-term stability of the dispersion in the fluid. Numerous studies have reported the dispersion stability; functional groups attached to nanoparticles play a role in causing steric hindrance and have an affinity for the surrounding fluid, resulting in preserving the dispersion. We investigate the structural effects on dispersion by molecular dynamics simulations of nanofluid containing graphene sheets with functional groups of varying lengths at the surface. The results demonstrate that short functional groups were too short to cause significant steric hindrance, while relatively longer functional groups tended to stack onto the graphene sheets, leading to trapping due to strong van der Waals interactions. Additionally, we discuss the minimum number of functional groups necessary for maintaining dispersion through calculations of the area of a single functional group.     

Molecular dynamics simulations of thermal expansion properties of single layer graphene sheets

The area coefficients of thermal expansion (CTEs) of perfect single layer graphene sheet (SLGS) and SLGS with vacancy defects of different distributions were calculated in this work through molecular dynamics (MD) simulations. The effects of some parameters such as temperature, SLGS size, sample area size, vacancy fraction and vacancy distribution on CTE were investigated extensively. Numerical results clearly revealed that for both perfect and defective SLGSs, the area CTEs are negative and nonlinear with the temperature variation within a wide temperature range. Moreover, the area CTEs tend to be more insensitive to the temperature when temperature is higher than 600 K. The area CTE of a perfect SLGS converges only when the SLGS size and the ratio of the sample size to the SLGS size is above a critical value. When the SLGS size or the sample size is small, the area CTE shows distinct size-dependence. In addition, a set of empirical formulations is proposed for evaluating the area CTEs of perfect SLGSs within a wide temperature range. For the SLGS with vacancy defects, the area CTE decreases with the increase of vacancy fraction within the temperature range considered. Furthermore, compared with a decentralised distribution of vacancy defects, a concentrated distribution leads to a smaller value of area CTE of SLGS, especially for the case of high vacancy fraction.     

Molecular dynamics simulation of replacement of methane hydrate with carbon dioxide

Molecular dynamics simulation was performed to analyse the phenomena of replacement of methane hydrate with carbon dioxide (CO₂) at 270 K and 5.0 MPa for 5300 ps. The methane hydrate phase was constructed with 16 unit cells of hydrate. Every cage in the hydrate was occupied by one methane molecule. The methane hydrate phase was sandwiched between two CO₂ phases. During the simulation the hydrate partially melted and liquid water phase appeared, and CO₂ dissolved in the liquid water phase. The replacements were observed three times at the hydrate–liquid water interface during the simulation. In the first case, the replacement occurred at a S-cage without changing the structure. In the second case, an M-cage of methane hydrate partially collapsed, and methane and CO₂ molecules exchanged. After the exchange, the cage occupied by CO₂ remained in the M-cage structure. In the third case, a S-cage of methane hydrate partially collapsed, and methane and CO₂ molecules exchanged. After the exchange, the cage occupied by CO₂ changed to an M-cage-like structure.     

Molecular insight into adsorption affinities of Carmustine drug on boron and nitrogen doped functionalized single-walled carbon nanotubes using density functional theory including dispersion correction calculations and molecular dynamics simulation

Abstract

We report a quantum mechanics calculation and molecular dynamics simulation study of Carmustine drug (BNU) adsorption on the surface of nitrogen (N) and boron (B) doped-functionalized single-walled carbon nanotubes. The stability of the optimized complexes is determined on the basis of relative adsorption energy (ΔE_ads). The ΔE_ads results claim that drug molecule tends to adsorb on the nitrogen and boron doped functionalized tubes with the energy values in the range of ?61.177 to ?95.806?kJ/mol. Based on the obtained results, it is observed that N-doping compared with B-doping has improved more effectively drug absorption on the surface of functionalized nanotube. The results of Atoms in Molecule calculations indicate that drug adsorbs molecularly via hydrogen bonds interactions on the surface doped-functionalized carbon nanotubes. Moreover, molecular dynamics simulation is performed to investigate the dynamics behavior of the drug molecules on the nitrogen-doped functionalized carbon nanotube (f-NNT) and functionalized carbon nanotube (f-CNT). The higher average calculated electrostatic and van der Waals energies as well as higher number of intermolecular hydrogen bonds in BNU-f-NNT in comparison with BNU-f-CNT model suggest the more effectual interaction between drug molecules and nitrogen-doped functionalized carbon nanotube.

Communicated by Ramaswamy H. Sarma     

The influence of nicotine on pioglitazone encapsulation into carbon nanotube: the investigation of molecular dynamic and density functional theory

In this work, molecular dynamics simulations of the insertion of pioglitazone into the nanotube with chirality (10, 10) at 400 K and 1 bar in the presence and absence of nicotine molecules and in different drug concentrations have been studied. The main aim is consideration of the effect of nicotine in the drug encapsulation process. The results indicate that encapsulation of pioglitazone could be attributed to the water flow via van der Waals and hydrophilic interactions. Because of the existence of the partial π–π interactions between aromatic rings of pioglitazone and the conjugated aromatic rings of nanotube, pioglitazone molecule can enter inside the nanotube. Some physical properties such as hydrogen bonding, number of contacts, also, the diffusion coefficient of the pioglitazone and water molecules, and variation of the center of mass have been calculated during the simulation. Furthermore, computing the electronic structure has also been done on model systems for quantitative determination of the adsorption energy (E_ads). The B3LYP/6-31G* level calculations on four different configurations of pioglitazone/carbon nanotube (CNT) and nicotine/CNT show that the interaction of drug with the inside of the nanotube is stronger than the other forms.     

Interfacial layering and orientation ordering of ionic liquid around single-walled carbon nanotube: a molecular dynamics study

Single-walled carbon nanotubes (SWNTs) tend to aggregate to heavily tangled bundles due to the strong van der Waals attraction. Ionic liquids (ILs) are a kind of newly proposed solvents in which SWNT can be physically well dispersed. In this article, the cylindrical interface has been investigated by molecular dynamics simulation between IL of 1-butyl-3-methylimidazolium tetrafluoroborate ([Bmim][BF4]) and an infinite long armchair (6,6) SWNT. The highly ordered structure of the cations and anions is elucidated by the simulation results. Two evident dense layers are found for both the cations and anions along the surface normal direction of the SWNT. In addition, we have observed two different orientation patterns of the cations in the first layer. In sublayer 1A, which is the nearest to the surface, the imidazolium rings of the cations prefer to be parallel to the surface, with a slight tilt angle less than 15°. In sublayer 1B, they tend to be perpendicular to the surface, with their butyl chains appearing in sublayer 1A. The [BF₄]^?  anions are found to cling to the nanotube surface with three fluoride atoms, also indicating a highly ordered orientation. The simulation results in this work provide a clue to understand the stabilisation and dispersion of SWNT bundles in ILs.     

Predicting doxorubicin drug delivery by single-walled carbon nanotube through cell membrane in the absence and presence of nicotine molecules: a molecular dynamics simulation study

Abstract

In order to study the interaction of the anticancer agent Doxorubicin with the single-walled carbon nanotubes with different diameters as drug delivery systems, the molecular dynamics (MD) simulations have been used. Also, for design and development of intracellular Doxorubicin drug delivery systems, a series of steered MD simulations are applied to explore the possibility of encapsulated Doxorubicin–carbon nanotube penetration through a lipid bilayer in presence and absence of Nicotine molecules at different pulling rates. Our simulation results showed that in spite of the adsorption of drug molecules on the outer sidewall of the nanotubes, the spontaneous localization of one Doxorubicin molecule into the cavity of the nanovectors with larger diameters is observed. It is found that the presence of Nicotine molecules in extracellular medium increases the required force for pulling nanotube-encapsulated drug as well as the required time for penetration process, especially at higher velocity. Also, the entering process of the Nicotine molecules into the carbon nanotube causes that the encapsulated drug molecule is fully released in the hydrophobic phase of the lipid bilayer.

Communicated by Ramaswamy H. Sarma     

Assessment of the chitosan-functionalized graphene oxide as a carrier for loading thioguanine,an antitumor drug and effect of urea on adsorption process: Combination of DFT computational and molecular dynamics simulation studies

In this study, the interaction thioguanine (TG) anticancer drug with the functionalized graphene oxide (GO) nanosheet surface is theoretically studied in both gas phase and separately in physiological media using the density functional theory (DFT) calculations. DFT calculations indicated the adsorption and solvation energies are negative for f-GONS/TG complexes which propose the adsorption process of TG molecule onto the f-GONS surface is possible from the energetic viewpoint. QTAIM calculations confirm the nature of partially covalent-partially electrostatic between drug and nanosheet. These results are sorely relevant that an approach for loading of TG molecule is the chemical modification of GO using covalent functionalization which can serve as a nanocarrier to load drug molecules. Moreover, to understand the effect of urea on the nature of the interaction between TG and f-GONS, molecular dynamics (MD) simulation was employed. The results indicated that in the presence of urea the adsorption process gets affected and leads to instability of system, while the affinity of the TG for adsorption onto GO surface is increased in pure water.

Communicated by Ramaswamy H. Sarma     

Effect of charge redistribution on the binding of DNA nucleotide to carbon nanotube in molecular dynamics

All-atom molecular dynamics (MD) simulations are performed to study the binding of DNA nucleotides with two carbon nanotubes (CNTs) with similar diameters but different chiralities. Two schemes for assigning partial atomic charges (PACs) are adopted: (I) using PACs obtained from isolated DNA nucleotide and CNT optimised in vacuum, and (II) using PACs obtained from optimising nucleotide-CNT hybrid in solution. The former approach is what most MD simulations have used in the study of DNA-CNT hybrids, while in the latter approach, a redistribution of the PACs has occurred upon the hybridisation. Our results show that the charge redistribution has a profound effect on the dynamics of binding. In particular, PACs obtained from (II) lead to more stable binding structures in the MD simulations. The findings suggest that care should be taken in simulating DNA-CNT interactions using the classical force field approach.     

Molecular dynamics study of carbon nanotube oscillator on gold surface

We investigated the substrate effect of carbon nanotube (CNT) oscillators using classical molecular dynamics simulations. Double-walled CNT oscillators on {100} gold surface were considered. The nanotube–gold interactions induced the compressive deformations of the outer nanotube and affected the transitional velocity and the energy dissipation of the nanotube oscillator. When the inner nanotube was extruded from the outer nanotube, the central regions of the outer nanotube were compressed by the nanotube–gold interactions and then, these compressive forces pushed out the inner nanotube and finally, the transitional velocity of the inner nanotube was slightly increased at the edges regions. Since the energy dissipation of the nanotube oscillator on gold surface was higher than that in vapor, the decrease of the transitional velocity for the nanotube oscillator on gold surface was greater than that for the nanotube oscillator in vapor.     

Molecular dynamics simulations of ribonuclease T1

Molecular dynamics simulations in vacuum and with a water sphere around the active site were performed on the 2GMP-RNase T1 complex. The presence of water led to the maintenance of the 2-GMP-RNase T1 interactions as compared to the X-ray structure, including the hydrogen bonds implicated in the enzyme-inhibitor recognition process. The sidechain of His92 in the molecular dynamics water simulation, however, hydrogen bonds directly to the phosphate of 2GMP in contrast to the X-ray structure but in support of the role of that residue in the enzyme's catalytic mechanism. Fluctuations of activesite residues are not strongly influenced by water, possibly owing to the exclusion of water by the bound 2GMP, which did show an increase in mobility. Analysis of the 2GMP-RNase T1 interactions versus time reveal an equilibrium fluctuation in the presence of water, leading to a less favorable 2GMP-RNase T1 interaction energy, suggesting a possible relationship between picosecond fluctuations and inhibitor dissociation occurring in the millisecond time domain.Abbreviations RNase T1 Ribonuclease T1 (EC.3.1.27.3) - 2GMP Guanosine-2-monophosphate - SBS Stochastic Bondary Simulation - VS Vacuum Simulation - MD Molecular Dynamics     

Molecular dynamics modelling of tethered organics in confined spaces

A computational method for constructing and evaluating the dynamic behaviour of functionalised hexagonal mesoporous silica (HMS) MCM-41 models is reported. HMS with three pore diameters (1.7, 2.2 and 2.9 nm) were prepared, and, from these, two series of derivative structures were constructed – one with 1,3-diphenylpropyl (DPP) tethers and the other with smaller dimethylsilyl (DMS) tethers attached to the mesopores' internal surfaces. Comparison with experimental data shows that simulation results correctly predict the maximum tether density that can be achieved for each tether and each pore diameter. For the smaller pore models, the extent of DPP functionalisation that can be achieved is limited by the available pore volume. However, for the larger pore model, the extent of functionalisation is limited by access to potentially reactive sites on the pore surface. The dynamic behaviour of the models was investigated over a range of temperatures (240–648 K). At lower temperatures ( < 400 K), the mobility of DPP tethers in the 2.9 nm model is actually less than that observed in either the 2.2 nm model or the 1.7 nm model due to the extensive non-bonded interactions that are able to develop between tethers and the silica surface at this diameter. At higher temperatures, the free ends of these tethers break away from the surface, extend further into the pore space and the DPP mobility in the 2.9 nm model is higher than in the smaller pore systems.     

Molecular dynamics simulations of Alzheimer's beta-amyloid protofilaments

Filamentous amyloid aggregates are central to the pathology of Alzheimer's disease. We use all-atom molecular dynamics (MD) simulations with explicit solvent and multiple force fields to probe the structural stability and the conformational dynamics of several models of Alzheimer's beta-amyloid fibril structures, for both wild-type and mutated amino acid sequences. The structural models are based on recent solid state NMR data. In these models, the peptides form in-register parallel beta-sheets along the fibril axis, with dimers of two U-shaped peptides located in layers normal to the fibril axis. Four different topologies are explored for stacking the beta-strand regions against each other to form a hydrophobic core. Our MD results suggest that all four NMR-based models are structurally stable, and we find good agreement with dihedral angles estimated from solid-state NMR experiments. Asp23 and Lys28 form buried salt-bridges, resulting in an alternating arrangement of the negatively and positively charged residues along the fibril axis that is reminiscent of a one-dimensional ionic crystal. Interior water molecules are solvating the buried salt-bridges. Based on data from NMR measurements and MD simulations of short amyloid fibrils, we constructed structural models of long fibrils. Calculated X-ray fiber diffraction patterns show the characteristics of packed beta-sheets seen in experiments, and suggest new experiments that could discriminate between various fibril topologies.     

Molecular dynamics simulation of intrinsically disordered proteins

Intrinsically disordered proteins are biomolecules that do not have a definite 3D structure; therefore, their dynamical simulation cannot start from a known list of atomistic positions, such as a Protein Data Bank file. We describe a method to start a computer simulation of these proteins. The first step of the procedure is the creation of a multi-rod configuration of the molecule, derived from its primary sequence. This structure is dynamically evolved in vacuo until its gyration radius reaches the experimental average value; at this point solvent molecules, in explicit or implicit implementation, are added to the protein and a regular molecular dynamics simulation follows. We have applied this procedure to the simulation of tau, one of the largest totally disordered proteins.